目前大多数生态位研究建立在物种水平，将物种作为单个无差异的实体来响应气候变化。一些研究认为物种水平的生态位是保守的，也有研究表明物种生态位发生了分化。物种是一系列在分化时间、进化独立性和遗传特殊性上不同的进化谱系的集合。对于具有系统地理结构的广泛分布的物种，应用物种水平上综合谱系的生态位模型（ENMs）难于准确描述适应不同气候条件下的不同遗传支系间的生态位进化动态。此外，种内生态位动态下不同支系的形态进化和表型可塑性亟待观测研究。 

 青藏高原受印度洋板块和亚洲板块碰撞的影响持续向北抬升造就了青藏高原周边山区错综复杂的高山峡谷及剧烈的气候波动，这些变化为物种内支系分化及生态位进化提供了多样化的生态机遇。西藏齿突蟾广泛分布于青藏高原东南边缘，跨喜马拉雅、青藏高原、横断山及岷山等毗邻山系，海拔从2400m到5270m都有分布，如此广泛分布的物种为研究种内支系间生态位进化及形态进化动态提供了很好的研究材料。

　　图1.各支系发生位点分布图. 

 成都生物所两栖爬行动物多样性与保护项目组江建平研究员和谢锋研究员带领的团队联合首都师范大学等多团队基于多年科考采样的生物学数据，结合世界气候数据（WorldClim database），利用生态位模型（ecological niche models, ENMs）、排序法（environment principal component analyses, PCA-env）及系统发生对比方法（phylogenetic comparative methods, PCMs）和系统发生最小二乘法（phylogenetic generalized least squares, PGLS）探究了西藏齿突蟾支系间生态位分化水平、生态位分化的具体形式、生态位分化的贡献最大的气候因素以及形态进化的响应。

　　图2. (a)谱系发生关系及基于最小凸边形法获得的地理范围; (b)各支系适宜分布区; (c)沿各气候主成分(PCs)及海拔(Elevation)的祖先关系重建. 

 结果表明：基于线粒体基因组在遗传明显分化为6个支系，其中有3个（支系 E. A、E. B和E. C）发生了显著的生态位分化（基于ENMs和PCA-env方法：Schoener's D小于零假设预期Schoener's D值，双向生态位相似性测试P < .05，则零假设被拒绝，即生态位发生分化）。在该案例中，生态位的分化主要是由生态位扩张导致的而不是生态位填充（生态位扩张：转移生态位vs. 本地生态位 E > U；生态位填充：转移生态位vs. 本地生态位U > E，详见原文Table 3）。环境主成分PC3（主要由BIO2和BIO3组成）、PC4和PC5（主要由BIO6和BIO7组成）维度环境因子遵循布朗运动沿系统发育关系促进了生态位的进化，这一结论与促进头体长进化的环境因子几乎一致。而且，由生态位扩张引起的生态位分化的确伴随着形态上的进化（例如：头体长（SVL）具有强系统发育信号）与表型可塑性的发生（如：前臂及手长（LAHL），腿全长（HLL）和足长（FL））。本研究为广泛分布种的种内各支系间生态位进化及形态特征变化提供了研究模型和理论框架。 

　　图3. (a)环境变量的3D-PCA; (b)各支系沿气候变量BIO3的生态位分布; (c)参与PCA-env分析的各因子相关性及主成分贡献; (d-e)生态位等效性测试结果; (f-g)生态位相似性双向测试结果; (h-i)生态位动态分析结果. 红色代表扩张生态位，绿色代表稳定生态位，蓝色代表扩张生态位与稳定生态位重叠区域. 

 本研究报道以题为“Climatic-niche evolution with key morphological innovations across clades within Scutiger boulengeri (Anura: Megophryidae)”，发表在国际著名生态与进化类期刊Ecology and Evolution [2021, DOI: 10.1002/ece3.7838]。中国科学院成都生物研究所在读博士生蔺秀琴为该论文第一作者，谢锋研究员为通讯作者。该研究得到第二次青藏高原综合科学考察研究计划（STEP, 2019QZKK05010503）、生态环境部生物多样性调查评估项目（2019HJ2096001006）和四川省科技厅基础条件平台建设（2019JDPT0020）及中国生物多样性观测网络（Sino BON）等项目资助。

 原文链接