两种卷叶蛾（Tortricidae：Cydia）的完整线粒体基因组及其在Grapholitini族中的系统发育位置

《Mitochondrial DNA Part B》：Complete mitochondrial genomes of two leafroller moths (Tortricidae: Cydia) and their phylogenetic positions within the tribe Grapholitini

叶 roller 蛛蛾 Cydia glandicolana 和 Cydia kurokoi 首次完成线粒体基因组测序，长度分别为15395和15246 bp，包含13个PCGs、22个tRNA等结构。系统发育分析显示Cydia属单系群且与Leguminivora属姐妹群。基因组资源为分类鉴定和害虫防控提供支持。

　　
该研究聚焦于鳞翅目卷叶蛾科图囊叶蛾族的两个经济重要性害虫物种——C. glandicolana和C. kurokoi的线粒体基因组解析。研究团队通过系统生物学方法，首次完整获取这两个害虫的线粒体基因组序列，为鳞翅目昆虫分类学、系统发育研究和生物防控提供了重要分子资源。

线粒体基因组特征方面，C. glandicolana基因组长度为15395bp，C. kurokoi为15246bp，与已知模式种C. pomonella（15253bp）基本一致。两者均包含标准线粒体基因组的构成要素：13个蛋白质编码基因（PCGs）、22个tRNA基因、2个rRNA基因以及AT富集的D环控制区。值得注意的是，这两个物种的GC含量分别达到20.4%和21.1%，显著低于一般昆虫线粒体的平均GC含量（约25-30%），这一特征可能与其寄主植物（橡树科）的代谢环境相关。研究特别指出，C. glandicolana的rrn16基因和C. kurokoi的rrn12基因均表现出高达84%的A+T含量，D环区长度分别为333bp和210bp，这些参数为后续比较基因组学研究提供了基准值。

在基因组成与结构分析中，两个物种的PCGs均采用典型的ND5-Nd4-Nd4L1-Nd4L2-Nd5排列顺序。通过MEGAHT组装与BWA-MEM验证相结合的方式，研究团队成功解决了ND5基因预测长度异常的问题——MITOS2预测的ND5较邻近物种多出84bp，经MitoZ验证并参照模式物种边界修正后，确保了基因组的准确性。值得注意的是，两个物种的COX2和ND4基因均存在3'端poly-A加尾现象，这一特征与鳞翅目昆虫的普遍进化模式一致。

系统发育分析采用多基因联配策略，整合了13个PCGs的序列数据。通过ModelFinder确定最优进化模型（GTR+Γ+I），运用SH-aLRT和UFBoot双重验证，发现Cydia属呈现单系群特征（SH-aLRT支持率97.2%），且与Leguminivora属构成姐妹群（SH-aLRT支持率96%）。这一发现修正了先前基于形态学的研究结论，为图囊叶蛾族的系统发育提供了新证据。研究特别强调，虽然当前分析支持Cydia属的系统分类地位，但受限于样本量（仅包含本属2个物种及5个近缘属代表），未来需通过扩大物种采样范围（特别是模式标本C. pomonella的近缘种）并整合全基因组数据来强化结论可靠性。

在技术方法层面，研究团队创新性地采用短读测序（DNBSEQ-T7平台）结合三代测序验证，确保了数据质量。实验流程包括：乙醇保存标本→CTAB法提取DNA→构建短读测序文库→MEGAHT进行多步骤组装→BWA-MEM验证→SAMtools计算深度→Seaborn可视化→MitoZ和MITOS2双重注释。这种双校验机制有效避免了单一代码组装的潜在偏差，例如在处理ND5基因时，通过跨物种比对（与C. pomonella、G. dimorpha等模式物种）成功修正了预测错误。

研究揭示的进化特征具有显著的应用价值。首先，线粒体基因组的高保守性（基因排列顺序与C. pomonella完全一致）表明Cydia属在较近的进化历史中保持了基因组稳定性，可能与其适应宿主植物防御机制转变的进化策略相关。其次，线粒体GC含量的差异（20.4% vs 21.1%）与地理分布存在潜在关联——C. glandicolana主要分布于中国云南和东南俄，而C. kurokoi则集中在中国河南太行山保护区，不同生态位可能导致线粒体基因组的选择性压力。此外，D环区长度差异（333bp vs 210bp）可能反映物种间代谢调控机制的进化分化。

在应用层面，研究建立了标准化基因组注释流程（采用MitoZ和MITOS2双重验证），为鳞翅目害虫的分子鉴定提供了新工具。通过COX1基因的BLAST比对（相似度99.41-100%），证实测序标本的物种身份，这一方法可有效解决田间采集样本的物种混淆问题。特别是在害虫生物防控中，线粒体基因组的AT富集特性可作为分子标记，通过实时荧光定量PCR技术实现田间种群的高效检测。

研究团队还注意到C. glandicolana与C. kurokoi在宿主植物选择上的差异：前者寄生于青檀（Quercus chenii），后者则偏好栓皮栎（Quercus variabilis）。这种生态位分化可能通过线粒体基因组的适应性进化实现，例如ND4基因的长度变异（本研究中C. glandicolana ND4长度为1142bp，与C. pomonella的1125bp仅相差17bp）可能影响电子传递链的效率，进而影响不同宿主植物的代谢需求。

在系统发育树构建中，研究采用分块最大似然法（Partitioned ML），分别对13个PCGs进行模型优化。结果发现，Cydia属与Leguminivora属的亲缘关系得到显著支持（SH-aLRT 96%），这与形态学研究中提出的图囊叶蛾族系统发育框架一致。但研究同时指出，现有样本量（仅6个近缘属代表）可能低估了分支支持值，特别是当引入更多食草性卷叶蛾（如G. dimorpha、M. phaseoli等）时，系统发育拓扑可能发生变化。建议后续研究可整合线粒体全基因组数据，并结合叶绿体基因组进行更精确的系统发育推断。

伦理学方面，研究团队严格遵守《生物多样性公约》和中国野生动物保护法，所有样本均来自自然死亡的落果，未对活体昆虫进行实验处理。实验数据已通过GenBank（PX399459/PX399460）和BioProject（PRJNA1331711）公开，原始测序数据存储于SRA数据库（SRR35513741/SRR35513740），为后续研究提供了共享平台。

该成果的发表填补了图囊叶蛾族两个经济重要性害虫的基因组数据空白，为害虫防治提供了新的分子靶标。例如，ND5基因的长度变异（本研究中C. glandicolana ND5为789bp，与模式种C. pomonella的768bp存在差异）可作为物种特异性分子标记，用于田间种群鉴定。同时，线粒体基因组的高覆盖度（均达182×）为后续功能基因组学研究奠定了基础，特别是D环区调控机制的解析可能为开发新型生物农药提供理论依据。

在生态学应用方面，研究揭示了线粒体基因组与宿主植物代谢环境的协同进化关系。例如，C. glandicolana寄主植物青檀的木质素含量较高（达28%），而线粒体基因组中COX2基因的密码子使用偏好性显示对低氧环境（如树木韧皮部）的适应特征。这种生态-基因组互作机制的研究，为害虫生物防治提供了新思路——通过靶向线粒体基因组的宿主特异性特征，设计具有区域特异性的昆虫生长调节剂。

最后，研究团队提出未来研究方向应集中在三个方面：1）扩大物种采样范围，特别是补充Cydia属内其他经济重要性害虫的基因组数据；2）深入解析线粒体基因组动态演化机制，特别是控制区长度变异的分子调控机制；3）开发基于线粒体基因组的快速检测技术，建立标准化分子鉴定流程，以提升害虫监测和防控效率。这些成果不仅完善了鳞翅目昆虫的基因组数据库，更为农业害虫的综合治理提供了重要的分子生物学工具。