两栖动物皮肤的功能包括离子及水交换(Bentley & Yorio，1979；Marrero & Hillyard，1985；Sullivan et al., 2000)、感觉(Kitson & Roberts，1983；Mearow & Diamond，1988)、机械与化学防御(Fox，1986；Evans & Brodie，1994；Woodhams et al., 2007；Jeckel et al., 2015)和呼吸(Duellman & Trueb，1994；Nussbaum & Wilkinson，1995；Bickford et al., 2008；Robischon，2017)。对皮肤各功能对应结构特征的量化常借助石蜡切片技术，涉及如分裂相表皮细胞比例(耿欣莲，1959)、皮肤各层厚度(Greven et al., 1995；Barrionuevo，2017；Yang et al., 2019)、腺体数量和大小(de Almeida et al., 2007；Gonçalves & de Brito-Gitirana，2008；Yang et al., 2019)等。特征值的准确获得受样本设计影响。目前在石蜡切片技术领域，即使针对相同特征，已发表工作所采用的设计(包括切片数、测量数等)也少见统一(Norris et al., 1989；刘满樱等，2007；李洋等，2009；Prates et al., 2012；Brunetti et al., 2016；Barrionuevo，2017；Ponssa et al., 2017；Wanninger et al., 2018)。例如皮肤各层厚度的测量，曹燕等(2011)对齿突蟾属Scutiger 4个物种背腹面等距测量5次，而Chammas等(2015)对疣蟾Rhinella granulosa背腹面测量30次，Ponssa等(2017)对细趾蟾属Leptodactylus物种背部不同部位测量5次，VanBuren等(2019)对白唇树蛙Litoria infrafrenata背腹及后肢腹面随机测量10次。

不同结构特征的量化，缘其自身特点，对样本设计的要求理当不同。实际工作中，同一批次石蜡切片样品往往为不同特征提供数据。通过直接测量即可获得数据的指标，如厚度可能需要较小的样本量。而另外一些无法直接测量的特征则可能需要更大工作量的取样，进而影响制备切片样品的策略和工作量。因此，对后者样本设计的探讨有利于实践工作。

皮肤气体交换主要发生于表皮下毛细血管，其密度为衡量呼吸能力的一项重要参数(Krogh，1904；Feder & Burggren，1985；Malvin, 1988, 1993)。通过向两栖动物活体心脏注射墨水，可染色毛细血管并测量其单位面积内的长度(Czopek，1965；de Saint-Aubain，1982)。但基于动物伦理的考量，该方法近年来已罕有使用，现常用石蜡切片计数评估皮肤毛细血管丰富度。基于毛细血管计数和切面长度，单条切片仅可提供一个丰富度值。为获得有代表性的结果，需要相当的样本量。在长时间完全没入水中的情况下，一些蛙类有战斗和鸣叫等依赖于高效率皮肤呼吸的行为，但其相应皮肤毛细血管特征尚待研究。这些蛙类中生活史和行为学研究开展较多的物种之一为峨眉髭蟾Leptobrachium boringii(Zheng et al., 2011；Hudson & Fu，2013)。本研究以峨眉髭蟾10个部位表皮下毛细血管丰富度为例，通过对样本进行重抽样，探讨取样与结果稳定性的关系，以期为两栖类皮肤切片的样品制备策略提供参考。

峨眉髭蟾隶属角蟾科Megophryidae，主要分布于四川盆地周边山地，每年2—4月于溪流中繁殖，其余时间陆栖。2只峨眉髭蟾标本于2019年5、6月采自四川省峨眉山，均为非繁殖季节雄性成体，编号为IOZCAS09432(吻肛距=80.6 mm，大个体)、IOZCAS09433(吻肛距=72.2 mm，小个体)。该物种非繁殖季节时活动隐秘，陆栖个体多年来少有获取，对其陆栖个体皮肤毛细血管丰富度的报道有利于国内外同行据此进行对比分析。气体交换是蛙类皮肤的重要功能，体型较大的个体因其体表面积/体质量比较小，毛细血管丰富度较高。2只个体体型大小不同，若符合该预期则支持实测结果的可靠性。麻醉采用低于5%乙醇以减小刺激，固定液为8%甲醛，固定后流水置换甲醛24 h再转70%乙醇长期保存。切片前，每只蛙取头背部、躯干背部、前肢背侧、后肢背侧、体侧、下颌、胸部、下腹、前肢腹侧、后肢腹侧10个身体右侧部位皮肤，均为5 mm×6 mm长方形。对皮肤组织进行常规石蜡包埋切片，切片厚度6 μm。每块皮肤尽量获得约15张载玻片的连续切片，每张载玻片上为3列、每列15条切片。苏木精-伊红染色(Preece，1965)，树胶封片。

使用徕卡DM2500显微镜加装DFC495相机观察、拍照。每块皮肤测量5张载玻片、每隔2张测量1张；同一载玻片测量6条组织、每隔7条测量1条；同一条组织获取1张显微图像测量、大小为1 255 μm×942 μm。图像内，计数表皮下毛细血管数量，利用ImageJ version 1.52a(Schneider et al., 2012)测量表皮长度，得出单位长度表皮下毛细血管数为1个丰富度样本值。由此得到跨度约3.5 mm的30个样本值，结果以x表示。

编写R脚本对30个样本值无放回抽样，模拟1载玻片×1图像、1载玻片×2图像、……、1载玻片×6图像、2载玻片×1图像、……、5载玻片×5图像共29种样本设计，每种设计1 000次重复。不同设计的样本大小共16种情况，包括1、2、3、4、5、6、8、9、10、12、15、16、18、20、24、25。不同载玻片被抽取的皮肤切片/图像位置相同，以契合常见的等距取样策略(曹燕等，2011；Brunetti et al., 2012；Prates et al., 2012)。1 000个模拟结果的均值几等于真实结果。以真实结果的5%为评判阈值，即样本均值×0.05，若其大于1 000个模拟结果的标准差，则大多数模拟结果和真实结果的相似性高于95%、该模拟设计可靠。以此相似性标准，统计各部位获得接近真实结果所需样本量，记录满足该标准的3个最小样本量的中位数。为比较相同样本量、不同工作量的设计，对“1载玻片×5图像、5载玻片×1图像”和“2载玻片×5图像、5载玻片×2图像”这2对设计的模拟结果标准差进行配对t检验。对毛细血管丰富度和可靠设计数量进行Pearson相关检验，以探索前者对所需样本量的影响。各检验在SPSS 11.0中进行，数据处理前用Kolmogorov-Smirnov检验确认用于配对t检验的变量的差值和用于Pearson检验的变量不显著偏离正态分布(P值均大于0.05)。

2只个体表皮下毛细血管丰富度分布模式基本一致。与其理论上较低的体表面积/体质量比值相符，大个体各部位毛细血管丰富度多高于小个体(图 1)。两者在体侧和后肢腹侧都有较高密度的毛细血管，为10.5~15.3条/mm；而头背部是毛细血管密度均较低的部位，分别为4.6条/mm和6.1条/mm。

图 1 2只雄性峨眉髭蟾身体右侧10个部位表皮下毛细血管丰富度 Fig. 1 Densities of capillaries beneath the epidermis in 10 parts on the right side of 2 male Leptobrachium boringii

图选项

和预期相符，模拟结果随着样本量的增加趋近实测值(图 2：A)。各部位模拟结果有差异(表 1)：体侧、下颌、下腹、后肢腹侧获得接近真实结果所需最小样本量约为8，密度较低的头背部则需20或25，密度最低的小个体胸部无法获得接近真实值的模拟结果。大个体所需样本量更小、可靠样本设计共95个，小个体为46个；除前肢背侧相同外，大个体各部位可靠样本设计均多于小个体。Pearson相关检验表明，2只个体各部位毛细血管丰富度和可靠设计数量均呈显著正相关(图 2：B，C)，大、小个体相关系数分别为0.855和0.874(P值均小于0.01，n=10)。配对t检验显示，当样本量固定为5或10时，不同工作量的抽样设计对结果影响不明显，模拟结果标准差之间的差异无统计学意义(P值均大于0.1，n=10)。

图 2 模拟结果示例和毛细血管丰富度对所需样本量的影响 Fig. 2 An example of the simulation results and the relationship between capillary densities and sample sizes required A.小个体模拟结果，横轴显示模拟参数，如2*5代表观察2张载玻片、每载玻片5个图像，数值为样本量、与黑色柱形高度一致，箱形为模拟结果(x ± s和范围)，竖线代表模拟结果范围；B.大个体毛细血管丰富度与可靠设计数的关系；C.小个体毛细血管丰富度与可靠设计数的关系 A. Simulation results for the density of capillaries beneath the epidermis in the flank of the small male individual, horizontal coordinates are sampling designs, e.g., 2*5 means sampling 2 slides and 5 images per slide, numbers refer to sample sizes which are proportional to the heights of black bars; box plots show x ± s and range; B. the correlation between capillary densities and numbers of successful sampling designs observed in the large male individual; C. the correlation between capillary densities and numbers of successful sampling designs observed in the small male individual

图选项

本研究系统中，对一些部位如头背部、前肢背侧表皮下毛细血管丰富度的量化需要较多的切片。对比大、小个体，各部位表皮下毛细血管相对丰富度分布基本一致，有较小体表面积/体质量比的大个体毛细血管丰富度更高。两者体侧和后肢腹侧较高的毛细血管密度与在蛙类中已观察到的模式相同(Czopek，1965；Toledo & Jared，1993)。无尾两栖类后肢腹侧是皮肤水、气交换的关键部位，毛细血管密度一般较高(Roth，1973；Bentley & Yorio，1979)。这些一致性提示本工作的实测结果可靠。基于实测数据的模拟显示，小个体雄性多数部位需要接近20条组织或更高的样本量，大个体多数部位需要接近15条或更高的样本量。对比相同样本量、不同工作量的取样设计，基于较大工作量即更多制片的取样未得到显著理想的结果。但2对设计的样本量均较小，分别为5和10，不足以支持较少制片的策略，有待于在更合适的研究系统中探讨。一套皮肤石蜡切片可用于不同特征的量化，实际工作中也常如此，制片量不能仅针对需少量样本即可准确量化的特征。常规石蜡切片实验周期较长，但更多的展片、染色、封片并不会大幅增加工作量。相反，如果后期发现样本量不够、追加取样则会成倍耗费时间。另外，以毛细血管密度为例，一些特征在两栖动物皮肤不同部位的分布很不均一(Moalli et al., 1980)，追加取样可能难以对应部位、无法补救。因此，建议在两栖类皮肤石蜡切片工作中对全部有效切片进行制片，以满足不同特征对样本量的要求。

对比不同个体相同部位以及同一个体不同部位，表皮下毛细血管丰富度与所需样本量呈负相关。这显示来自高血管丰富度部位的测量值方差更小，这些部位的血管分布可能更均匀。这样的规律是否在其他两栖动物中存在、在什么样的毛细血管丰富度范围内存在，仍需考查。如有其普遍性，则对实际工作有指导意义。如前所述，同一物种的较小个体一般有更高的皮肤面积/体质量比，故可能有较低的毛细血管密度。两栖动物皮肤毛细血管丰富度分布有一定规律可循，部分物种的该特征已有量化(Czopek，1965；de Saint-Aubain，1982；Malvin，1988；Toledo & Jared，1993)。某些生活史特征如陆栖、水栖与该特征关联(耿欣莲，1959；唐以杰等，1999；Barej et al., 2010)。这些已有信息可否用于抽样策略的制定，有待在其他两栖动物类群中进行与本工作类似的相关性分析。需要指出的是，本工作以模拟结果的标准差是否小于实测结果值的5%为判断依据，该阈值在低丰富度部位显然更小，但这样的设定更符合以相似性评判结果稳定性的实践标准。

致谢: 衷心感谢中国科学院成都生物研究所曾一唯师姐、周星梅老师、曾晓茂老师的帮助与建议。