静止代谢率(resting metabolic rate，RMR)是指动物维持生理稳态的最低能量消耗(Biro & Stamps，2008；孙儒泳等，2019)。在内温动物和外温动物中，RMR分别指代的是基础代谢率(basal metabolic rate，BMR)和标准代谢率(standard metabolic rate，SMR)，即成年动物处于非繁殖、静止和消化吸收后的状态时，在热中区(内温动物)或特定温度(外温动物)测得的最低代谢速率。无论外温动物或内温动物，RMR在个体、种群或物种水平上都存在显著差异，且与生长、存活等适合度成分关联密切(Biro & Stamps，2010；Auer et al., 2015；Pettersen et al., 2016；Arnold et al., 2021)。因而探究这些差异的形态、行为、生活史、系统发育方面的原因一直是动物生理生态学、进化生理学的研究重点。

RMR应该与心脏、肝脏、肾脏、小肠等高耗能器官的质量显著正相关(Daan et al., 1990；Hammond & Diamond，1997；Koteja，2000；Ksiazek et al., 2009)。代谢能力更强的有机体可以负荷更高的持续能量支出(如每日能量支出、持续代谢率或者最大代谢率)，但它们同时须承受更高的维持代价(即RMR)(Careau et al., 2008)。按照形态功能适应的原则，更高的持续能量支出需要更强大的高耗能器官(主要用于食物同化过程)与之匹配(又称“摄食增量假说”)。这个假说预测，在相关选择作用下，更大的高耗能器官需要更高的RMR以维持其基础运行。该预测在内温动物中得到较普遍的支持。尽管小鼠Mus musculus的小肠、肝脏、肾脏和心脏总重只占体质量的17%，但却贡献了个体间RMR 52%的差异(Konarzewski & Diamond，1995；Ksiazek et al., 2004)。此外，人工选择BMR更高或食物摄入量更多的小鼠品系会导致小肠、肝脏、肾脏和心脏的比重增加(Selman et al., 2001；Ksiazek et al., 2004)。扣除体质量的效应，布氏田鼠Lasiopodomys brandtii的BMR和心脏、肝脏、肾脏以及盲肠干重显著正相关(Song & Wang，2006)。在22种鸟类中，心脏和肾脏的组合相对质量也解释了其BMR 50%的变化(Daan et al., 1990)。然而，该假说很少在外温脊椎动物研究中进行验证，已有的少数研究结论也不一致(Biro & Stamps，2010)。

与器官形态特征不同，生理特征具有较强的可塑性。一般认为RMR的可重复性为0.3~0.4。因此，单次测定往往会与其真实值有较大偏离(Careau & Garland，2015)。进化生理学家极力主张采用重复测量的研究范式，结合随机回归模型分解生理特征的个体间和个体内方差及其相关性(Roche et al., 2016；Careau，2017)。然而，在过往RMR与器官大小的研究中，还未见该实验设计和统计方法的应用。

中华蟾蜍Bufo gargarizans是东亚地区演化最为成功的两栖类之一，分布在20°~50°N，海拔分布从太平洋岛屿延伸至青藏高原(Zhan & Fu，2011；Lu et al., 2021)。在外温脊椎动物中具有较强的代表性。本研究拟从野外采集中华蟾蜍成年个体，重复测定RMR 3次，利用马尔可夫-蒙泰卡罗双变量混合线性模型(MCMCglmm)(Hadfield，2010)，探究RMR和高耗能器官相对质量的个体间相关性。预测扣除体质量影响后，中华蟾蜍的RMR和小肠、肝脏、肾脏和心脏正相关，支持“摄食增量假说”。

1 研究方法 1.1 实验动物

成年中华蟾蜍为在繁殖季节采自川西山区多个海拔样点(雄性27只，雌性19只)。动物捕获后运送至中国科学院成都生物研究所动物房饲养。动物房内温度控制在20 ℃±1 ℃，光周期设定为12 h∶ 12 h(08∶ 00—20∶ 00开灯)。动物单笼饲养在长35.5 cm×宽25 cm×高15 cm的塑料盒内，盒内放置1块湿润的海绵提供水分，以及1块U形瓦片用于躲藏。饲喂蟾蜍的主要食物是沾有钙粉(Exo Terra)的黄粉虫Tenebrio molitor，每2 d喂食1次。此外，每周给蟾蜍饲喂蟋蟀1次(2~3只)。所有蟾蜍均在实验室驯化3个月以上。所有动物操作均按照中国科学院成都生物研究所动物福利委员会的要求执行。

1.2 呼吸代谢测定

RMR采用开放式呼吸代谢仪(TSE Systems，德国)在21 ℃环境下测定(Tan et al., 2021)。代谢测定前后均测量动物体质量(精确到0.1 g)。实验测定前，动物被转移到900 mL的亚力克材质的圆柱形呼吸室中适应2 h。呼吸室内均放置10 cm×10 cm湿纸巾以防止动物在测定过程中脱水。代谢测定时，首先利用气体流速控制器(G246)将压缩空气的气体流速控制在200 mL·min^-1。然后，气流开关(G244)将穿过呼吸室的气流与高氯酸锰的干燥柱连接(去除其中的水汽)，再将干燥过的气流导入O₂和CO₂分析装置中(S104氧气分析仪和QS151氧气分析仪)。因为中华蟾蜍为夜行性动物，选择在20∶ 30—06∶ 30测定RMR。代谢仪每批测定5只动物，每次耗时5 h。所有个体测定重复3次，每次测定个体的顺序保持一致。实验数据使用TSE Systems配备的专用软件进行采集。RMR以测定期内5 min内连续、最低耗氧速率的平均值计算。同时计算呼吸熵(CO₂生成和O₂消耗的比例)用于推断其有氧代谢的底物利用状态(Withers，2001)。

1.3 器官形态测定

代谢测定结束后，采用MS222(0.2%M/V)深度麻醉中华蟾蜍(镊子夹挤后腿肌肉无条件反射)。随即分离各内脏器官(心脏、肝脏、脾脏、肺脏、肾脏、性腺)和胴体，并将各器官间的结缔组织剥离干净。消化道用清水清理内容物后，分割成大肠、小肠和胃。各器官用滤纸吸干水分后，用铝箔纸包裹放置于干燥箱中，在65 ℃下干燥直至质量不再变化。

1.4 统计分析

数据采用SPSS 21.0和R 4.0.3进行分析。胴体干重和RMR均采用对数转换使其符合正态分布(Shapiro-Wilk检验)。采用一般线性模型分析性别和总干重(待分析的目标器官没有被包括在总干重中)对各器官干重的影响，其中固定效应包括性别、总干重、性别和总干重的交互作用。单变量混合线性模型用于分析重复测定的SMR和呼吸熵，固定效应设置为性别、总干重、性别和总干重的交互作用，个体标签设定为随机效应。随机效应的显著性基于最大似然法进行检验。

采用MCMCglmm双变量混合线性模型探究RMR和各器官质量的相关性。反应变量和解释变量在分析前都进行x=0、SD=1的数据转换以去除截距差异对随机效应的影响。对于每一个协变量分析，性别和代谢测定次数设置成固定效应，总干重和去除目标器官的总干重分别设定为RMR和目标器官的协变量。个体标签设定为随机效应。模型的残差和随机效应的先验设定为：V=diag(2)，nu=0.002；V=diag(2)×0.002，nu=1.002，alpha. mu=rep[0.2]，alpha. V=diag[25²，2，2]。为计算后验分布，模型共运行500 000次迭代，每100次迭代重复抽样一次，舍弃前50 000次重复。由于器官质量缺乏重复测定，因而只有随机方差而无残差。然而，方差成分必须为正值，因此将该变量的残差指定为一个极小的正值(0.000 1)。

2 实验结果 2.1 中华蟾蜍各器官相对干重的性别差异

中华蟾蜍雌雄个体存在二型性，雌性各器官的绝对干重均大于雄性(图 1：A)。除去体质量影响后，雌性的脾脏、肾脏干重显著大于雄性(P＜0.05)，雄性的心脏、肺脏干重显著大于雌性(P＜0.05)，而肝脏、胃、小肠、大肠、消化道和胴体的性别间差异不显著(P＞0.05；图 1：B；表 1)。

2.2 性别和体质量对RMR的影响及其可重复性

总干重显著影响RMR(P＜0.001)，但对性别和重复测定的时序没有影响(P＞0.05；表 2；图 2)。呼吸熵也不受性别和重复测定的时序影响(P＞0.05；表 2)。RMR和呼吸熵的重复性分别为0.321和0.201。

2.3 器官干重与RMR的相关性

MCMCglmm分析结果显示，在个体水平上，总干重校正的RMR和器官大小都无显著相关性(P＞0.05；表 3；图 3)。

3 讨论

中华蟾蜍的高耗能器官质量与RMR不存在相关性：1)在该物种中，RMR和持续能量支出间可能不存在相关选择或其选择压力较小。无论外温动物还是内温动物，种间比较都支持持续能量支出和RMR显著正相关(Taigen，1983；Bozinovic，1992；Merritt et al., 2013)。然而，种内分析二者相关性往往出现不一致的结果。例如，在个体水平上，根田鼠Microtus oeconomus和高原鼠兔Ochotona curzoniae夏季的BMR与最大代谢率并不相关，而冬季根田鼠BMR与最大代谢率显著相关(r=0.612)(王玉山等，2002)；花园石龙子Lampropholis delicate的RMR和最大代谢率不相关(Merritt et al., 2013)；在福勒蟾蜍Bufo woodhoursei fowleri中二者虽然相关，但相关系数仅为0.37(Walton，1988)。因而，何种生态、进化及生活史因素可以影响RMR和持续能量支出的关联性值得进一步研究。2)在该物种中，持续能量支出的个体差异不是由高耗能器官的质量大小变异来贡献的。在美洲豹蛙Rana pipiens中，仅见RMR和肾脏干重残差显著正相关，而与心脏、肝脏、小肠和腓肠肌干重不相关(Steyermark et al., 2005)。然而，该相关性受单一数据奇值点影响。这说明，即便是用肾脏来解释RMR的个体间变异，也存在较大的不确定性。在安第斯蟾蜍Bufo spinulosus中发现，具有更高RMR的个体拥有更大的小肠、肝脏、肾脏和心脏(Naya et al., 2009)，但该研究将冬眠组、短期禁食组和对照组样本合并进行了关联分析。而冬眠、短期禁食处理可能消耗高耗能器官中的脂肪等低代谢身体成分，潜在地增加了数据的极化程度，进而导致器官质量与RMR相关。因此，利用高耗能器官的大小来阐释持续能量支出的个体差异在两栖类中可能并不具有普遍性。将来的研究有必要聚焦于高耗能器官的生理生化特征(如柠檬酸循环或电子呼吸传递链过程中的限速酶)，来阐释高耗能器官与RMR的关联性(Bozinovic，1992；Trzcionka et al., 2008；Koch et al., 2021)。

本研究结果也证实，在开展整体生理学特征相关性研究时，进行重复测量设计的必要性。中华蟾蜍的RMR重复性仅为0.321，许多个体在3次测定过程中的RMR出现大幅度波动。如果仅依据单次测定结果进行相关分析，必然会偏离其真实值。实际上，利用单(多)变量随机回归模型分析重复测量的数据集可以解析何种因子限制特征稳定性、可塑性、进化权衡等一系列有趣的进化和生态问题，且日益受到学界关注和重视(Dingemanse & Dochtermann，2013；Sih et al., 2015；Videlier et al., 2019)。因此建议国内动物生理生态学者们在采用非入侵性生理学指标测定时充分利用这一设计理念。

本研究采用的动物来自不同海拔位点，而不同海拔的持续能量支出、RMR和体质量的尺度关系可能存在差异，从而影响器官大小与RMR的关系。然而，前期研究发现不同海拔中华蟾蜍的RMR没有差异，而且其与体质量的尺度关系也不随海拔变化(Tan et al., 2021)。而且，这些动物经过较长时间的室内驯化，很大程度上抵消了采用前即时环境的影响。因此，理论上来说，多个海拔位点取样不会对本研究结果产生显著影响。

本研究发现，中华蟾蜍的高耗能器官大小与其RMR不相关，提示该物种的持续能量支出与其高耗能器官质量不相关，或者持续能量支出与RMR不相关，因而不支持“摄食增量假说”。该研究为扩展外温动物持续能量支出理论奠定了实验基础。