蚧总科物种基因组进化与水平基因转移模式：揭示其环境适应性与害虫治理新视角

《Ecology and Evolution》：Genomic Evolution and Patterns of Horizontal Gene Transfer in Coccomorpha Species

本研究通过比较基因组学分析，揭示了蚧总科物种在进化过程中经历两次全基因组复制（WGD）事件，并通过水平基因转移（HGT）获得大量源自细菌和古菌的基因，这些基因主要参与营养代谢过程。研究系统鉴定了包括ABC、COE、P450和UGT在内的关键解毒基因家族，为理解其环境适应机制和开发靶向治理策略提供了重要资源。

　　

1 引言

蚧总科昆虫是一类高度特化的植物寄生虫，全球已描述物种超过8000种，广泛分布于热带、亚热带和温带地区。作为重要的农林害虫，它们通过刺吸式口器取食植物汁液，直接损害植物组织并影响寄主活力。其排泄的蜜露还会促进烟霉病菌生长，进一步削弱寄主植物。部分种类还可传播植物病原体，加剧危害程度。蚧总科在长期进化过程中形成了独特的适应机制，包括固着生活方式、蜡质保护层分泌以及高度的寄主植物特异性。在分子层面，水平基因转移（HGT）和解毒相关基因家族的进化在其适应演化中扮演着关键角色。

2 材料与方法

研究选取了五个染色体级别的蚧总科物种基因组数据，包括<1>Acanthococcus lagerstroemiae、<1>Coronaproctus castanopsis、<1>Planococcus citri、<1>Balanococcus diminutus和<1>Phenacoccus solenopsis。采用<1>Bactericera cockerelli和<1>Drosophila melanogaster作为外群。使用OrthoFinder鉴定单拷贝直系同源基因，通过MAFFT进行多序列比对，trimAl修剪低质量区域，RAxML构建最大似然系统发育树。分化时间估算使用r8s软件，以362百万年作为校准点。基因家族扩张与收缩分析采用CAFE5软件。HGT获得基因的鉴定使用HGTfinder v1.0计算Alien Index（AI）和out_pct值，并通过系统发育树验证。解毒相关基因通过HMMER搜索Pfam数据库中的特定结构域进行鉴定。

3 结果

3.1 蚧总科系统发育与进化见解

系统发育分析显示，五个蚧总科物种的分化时间介于333.18至84.22百万年前之间。其中<1>C. castanopsis最先从蚧总科谱系中分化出来（约333.18百万年前），其次是<1>A. lagerstroemiae（约226.93百万年前），而<1>B. diminutus和<1>P. citri的分化时间最晚（84.22百万年前）。同源基因的Ka/Ks分析表明，五个蚧总科物种经历了两次全基因组复制（WGD）事件。最早的WGD事件发生在88-95百万年前的白垩纪晚期，这一时期正值温室气候，大气CO₂水平升高，全球温度比现在高4°C–6°C，同时被子植物经历了快速的进化辐射。

3.2 蚧总科基因家族扩张与收缩

基因家族分析显示，五个蚧总科物种的基因家族发生了不同程度的扩张和收缩。其中<1>P. solenopsis显示出最高数量的基因家族扩张（955个），而<1>C. castanopsis则表现出最高数量的基因家族收缩（1698个）。五个物种共有的基因家族有5054个，这些可能负责基本的生物学功能。扩张的基因家族主要富集在氧化酸代谢过程、有机酸代谢过程和羧酸代谢过程等抗氧化相关过程中。

3.3 蚧总科中HGT获得基因的鉴定

研究共鉴定出260个HGT获得基因，这些基因主要来源于细菌和古菌。细菌来源中，假单胞菌门（包括<1>Wolbachia和<1>Rickettsia属）贡献了最高比例（53%），其次是拟杆菌门（13%）和芽孢杆菌门（8%）。古菌来源的基因占2%（4个基因），其中广古菌门贡献了75%。不同物种间HGT获得基因的数量和来源存在显著差异，<1>P. citri从不同类群中获得了最多的HGT基因，而<1>A. lagerstroemiae的HGT获得基因数量最少。

3.4 蚧总科HGT获得基因的功能分析

GO富集分析显示，这些HGT获得基因主要参与赖氨酸代谢过程、糖胺聚糖生物合成过程和水溶性维生素生物合成过程等生物学过程。KEGG通路富集突出了丙酮酸代谢、氨基酸代谢（如丙氨酸、天冬氨酸和谷氨酸代谢）和脂肪酸延伸等代谢通路。这些结果表明HGT获得基因可能与细胞能量生产和氨基酸合成密切相关。

3.5 蚧总科解毒相关基因

研究鉴定了解毒相关基因家族，发现ABC、COE、P450和UGT是蚧总科中最大的四个解毒相关基因家族，每个家族在五个物种中均含有超过200个成员。相比之下，EPHX基因家族规模较小，总共只有4个成员。不同物种间解毒基因的数量和组成存在显著差异，<1>P. citri拥有最多的解毒相关基因（554个），而<1>C. castanopsis则最少（217个）。COE基因家族在五个物种中的基因数量变异最为显著，其中<1>C. castanopsis最少（39个），<1>P. citri最多（208个）。

4 讨论

研究表明蚧总科物种在进化过程中经历了复杂的基因组重排事件，包括两次WGD事件和大量的HGT事件。这些进化事件为其环境适应提供了遗传基础。扩张的基因家族主要富集在抗氧化相关过程中，表明蚧总科可能进化出了扩展的抗氧化机制来应对农药和其他外源化合物引起的氧化应激。HGT获得基因主要参与营养代谢过程，这与蚧总科专性取食植物韧皮部的习性密切相关。韧皮部富含糖类但缺乏氨基酸和维生素，而HGT获得基因可能帮助蚧总科补偿其共生基因组中营养和结构合成途径的缺失。解毒相关基因家族的变异反映了不同物种在生态位、食性偏好或历史农药暴露方面的差异。特别值得注意的是，<1>C. castanopsis作为一个新出现的寡食性物种，其解毒基因数量显著较少，这可能与其狭窄的寄主范围和对农药选择压力的有限暴露有关。相比之下，泛食性物种如<1>P. citri拥有更多的解毒基因，可能在多样化的寄主植物或农药压力下具有生存优势。这些发现为了解蚧总科的解毒机制和农药抗性进化提供了重要见解，并为开发害虫控制策略提供了潜在的分子靶标。