地面甲虫金属色变异的基因组景观：Nebria ingens复合体的遗传基础与生态适应

《Molecular Ecology》：The Genomic Landscape of Metallic Colour Variation in Ground Beetles

　　

甲虫的金属色变异是一种壮观的特征，启发了多种人类文化。然而，关于这一性状的遗传基础及其生态重要性所知甚少。在本研究中，研究人员描述了美国加利福尼亚州内华达山脉高山步甲Nebria ingens复合体中金属色变异的地理分布、光学机制、遗传基础以及生态与进化重要性。研究人员发现，鞘翅颜色在两个异域种（从黑色的N. ingens到绿色的N. riversi）之间连续变化，杂交种群呈现中间色；并证明该金属色由鞘翅上表皮中的多层反射器（multilayer reflector）产生。通过在高密度基因型变异中应用自然种群关联作图（wild-GWAS），研究人员鉴定了五个与金属色变异共变、且在角质层形成和色素沉着途径中具有已知功能的候选基因。这一发现，连同颜色变异与水分可利用性之间的显著相关性，表明金属色变异在N. ingens复合体中作为对环境变异的局部适应而进化。

金属色是甲虫中引人注目的表型特征，但其遗传基础与生态意义长期未被充分理解。以往研究多聚焦于金属色在种间识别、拟态和警戒中的作用，而较少关注其与局部环境适应之间的关系。在部分甲虫类群中，颜色变异已被证明与水分可利用性相关，提示金属色可能参与耐旱适应。然而，由于缺乏系统性的遗传学数据，控制金属色变异的基因和机制仍然不清。本研究以内华达山脉高山步甲Nebria ingens复合体为对象，该复合体包括两个地理隔离的异域种——鞘翅呈黑色的N. ingens与呈现绿色的N. riversi，并在分布区中部形成具中间颜色的杂交种群。微生境沿纬度梯度从北部的积雪环境过渡至南部的较干燥环境，为解析颜色变异的遗传基础及其生态功能提供了理想系统。研究人员利用全基因组重测序、高光谱成像、透射电子显微镜以及自然种群全基因组关联分析，系统研究了该复合体金属色变异的地理分布、光学机制、遗传基础及环境相关性，旨在鉴定控制金属色变异的候选基因，并检验其是否作为对水分可利用性变化的局部适应而进化。

研究人员共采集27个种群共404头N. ingens复合体成虫，使用374头个体的全基因组重测序数据进行变异鉴定，经质量过滤后获得355头个体的526,717个高质量单核苷酸多态性（SNP）和插入缺失（indel）变异用于关联分析。颜色表型通过高光谱相机（Specim IQ）测量368头个体的鞘翅反射光谱，并转换为CIELAB色度空间中的色度（chroma）值。采用透射电子显微镜观察绿色、深绿色和黑色鞘翅的上表皮微观结构，测量多层反射器的层数与厚度，并利用多层反射模型模拟反射光谱以验证光学机制。关联分析使用隐因子混合模型（LFMM）和偏冗余分析（pRDA），并通过设定潜在因子或主成分控制群体结构；同时使用线性混合效应模型检验雪水当量（SWE）与色度的关系，并使用线性回归检验色度与水分丧失率的关系。

**3.1 连续变异及金属色控制机制**
通过反射光谱测量，绿色鞘翅在500–600 nm波长处呈现反射峰，而黑色鞘翅振幅较低并失去特征峰。CIELAB色度转换显示，颜色变异从绿色到黑色连续分布，表明该性状为多基因控制的数量性状。透射电镜显示，绿色鞘翅的上表皮内含有1–3个电子致密层（厚50–100 nm）与电子透明层交替排列，构成多层反射器；深绿色鞘翅仅有一层电子致密层，黑色鞘翅中未观察到清晰的多层结构。基于层厚度计算的模拟光谱与观测光谱及实际颜色一致，且交替层数越多，绿色反射峰振幅越强，证明金属色由上表皮多层反射器产生，其变异属于由层数变化导致的饱和度变异。

**3.2 颜色相关基因及其作用**
在526,717个变异中，892个在LFMM分析中显著，573个在pRDA分析中为异常值，两者交集共195个变异，对应127个基因，其中65个具有功能注释。这些基因中，5个已知参与昆虫角质层形成和色素沉着途径：yellow-e和pdm3属于酪氨酸介导的鞣化途径，aristaless参与眼色素途径，chico编码胰岛素受体底物蛋白，属于胰岛素/类胰岛素生长因子信号通路，PK1-like编码焦激肽受体样蛋白，属于焦激肽信号通路。其中chico、PK1-like以及一个小G蛋白信号调节因子3（GPSM3）同源基因此前已在选择性扫描分析中被识别，提示其受到正选择。

**3.3 颜色变异、水分丧失率与SWE关联**
利用线性混合效应模型检验发现，鞘翅色度与采样点的雪水当量（SWE）显著正相关（t=5.58，p=2.51×10^?6），即来自SWE较高（积雪较多）区域的甲虫通常具有更鲜艳的金属色。水分丧失率与色度呈正相关趋势，但未达到传统显著阈值（p=0.0516，R²=0.068），表明色度只能解释水分丧失率6.8%的变异。这些结果支持金属色变异与水分可利用性存在关联，但水分保持能力可能还受其他因素影响。

讨论部分指出，N. ingens复合体的金属色由上表皮多层反射器产生，层数是决定反射峰振幅和颜色饱和度的关键。在候选基因中，yellow-e和pdm3最可能通过影响黑色素沉积参与多层反射器的形成，因为黑色素的分布决定了电子致密层与电子透明层的对比度。chico和PK1-like可能通过信号通路间接影响色素沉着，且具有选择信号，值得进一步研究；aristaless因参与不同色素途径而被认为可能性最低。生态方面，研究人员未发现该复合体存在基于颜色的择偶偏好，且作为高营养级夜行捕食者，其捕食压力较低，因此金属色更可能在湿润生境中起到伪装或水分调节的作用。黑色个体在较干燥生境中更常见，且倾向于具有更好的水分保持能力，但色度对水分丧失率的解释力有限。研究局限性包括颜色量化的测量误差以及大量候选基因缺乏功能注释，后续需通过表达调控或基因编辑实验验证候选基因功能。

综合而言，本研究表明N. ingens复合体的金属色变异由多层反射器的层数差异产生，并受多个已知色素途径基因的调控；颜色变异与雪水当量显著相关，提示金属色作为局部适应性状响应环境水分梯度而进化。这一发现为理解甲虫金属色的遗传基础和生态功能提供了重要证据，也凸显了野生全基因组关联分析在非模式物种表型遗传解析中的应用价值。