论文部分内容阅读
昆虫线粒体基因组具有基因组短、基因排列相对稳定、进化速率快、普遍为母系遗传的特征，目前已成为研究物种进化和系统发育的重要分子标记。随着系统发育分析方法的日益成熟多样和测序技术的不断升级，以线粒体基因组为数据源的系统发育研究已成为焦点领域。基于线粒体基因组数据重建系统发育关系探讨昆虫系统发育的研究，主要集中在物种丰富度较高的目级和科级阶元水平，而针对于亚科及更低级阶元的研究甚少。通过比较发现，基于线粒体基因组的较高级阶元系统发育研究广泛存在无法避免的问题，如：长枝吸引、碱基替换饱和等，可导致分枝结构不理想。而亚科级及以下的低级阶元在分析过程中极少出现类似问题。因此，线粒体基因组对于解决低级阶元系统发育关系表现出更好的适应性。 毛蚜亚科是半翅目蚜科的一个重要类群，包括2族，即毛蚜族和伪毛蚜族和12属，营同寄主全周期生活，而且两个族分别与木本植物和草本植物相关联，具有比较严格的寄主专化性。前期的研究主要集中在该亚科的物种鉴定与形态特征描述，寄主植物与生物学习性观察和地理分布等方面。有关该类群基于分子数据的相关研究很少，涉及的物种数量也非常有限；同时蚜虫线粒体基因组的获得相对比较困难，导致该亚科有关线粒体基因组的测序工作相对延迟，目前也没有关于毛蚜线粒体基因组的相关报道。规模化获取线粒体基因组信息，基于线粒体基因组数据重建毛蚜亚科内部系统发育关系，对了解毛蚜亚科的演化历史，探讨毛蚜与寄主植物之间的关系，推动蚜虫类线粒体基因组的研究水平具有重要意义。 本研究选取毛蚜亚科2族8属22种，对毛蚜线粒体基因组进行测序、注释和比较分析，阐述毛蚜亚科线粒体基因组的基本特征；利用贝叶斯法(BI)和最大似然法(ML)构建系统发育树，重建毛蚜亚科的系统发育关系。主要研究结果如下： （1）获取21种蚜虫约87％-95%的基因组全长，序列长度13459-15384 bp，1个蚜虫（Trichaitophorus foliatus）获得了约65%的基因组，序列长度为10355 bp。对测序结果进行注释，表明基因排列方式与祖先排列方式一致，未发现基因重排现象。由于假基因的干扰和A+T含量高，大部分物种未能获得tRNAArg- nad5区域的序列。所有蛋白质编码基因的起始密码子都为ATN，终止密码子为TAATAGT。间隔区大小范围为10-71 bp，最大的重叠区位于atp6和atp8之间，大小为20 bp。 （2）对毛蚜的碱基组成进行分析可知，A+T含量大于G+C含量，22个毛蚜物种A+T含量为81.55%-84.96%。在蛋白编码基因中，AT偏斜和GC偏斜都为负值，说明A少于T，G少于C。从密码子使用情况来看，Phe、Leu和Ile是编码最多的氨基酸， TTT（Phe）、TTA（Leu）和ATT（Ile）是使用次数最多的密码子。 （3）以粉毛蚜Pterocomma pilosum（蚜科：粉毛蚜亚科）为参考，计算毛蚜物种的同义替换率Ks、非同义替换率Ka以及非同义替换率Ka与同义替换率Ks的比值Ka/Ks,所有的毛蚜Ka/Ks值小于1，说明受纯化选择。除三角枫多态毛蚜Periphyllus acerihabitans外，伪毛蚜族的进化速率要大于毛蚜族的进化速率。 （4）利用BI和ML方法构建了毛蚜亚科系统发育树，两种方法构建的系统发育树拓扑结构一致。结果表明，在族级水平上，毛蚜族和伪毛蚜族的单系性得到了支持。在伪毛蚜族内，小毛蚜属和伪毛蚜属的物种聚在一起，因此使其成为并系群。赖毛蚜属的单系性成立并且位于伪毛蚜族的基部位置；在毛蚜族中，多态毛蚜属的Periphyllus acerihabitans和Periphyllus acericola与毛蚜族的其他物种形成姐妹群，在毛蚜族的其他物种中，三毛蚜属和桠毛蚜属的单系性得到了支持并与多态毛蚜属的物种聚为一支，形成（（Trichaitophorus+Yamatochaitophorus）+Periphyllus）的姐妹群关系。朗伯毛蚜属的Lambersaphis pruinosa聚到毛蚜属内，使毛蚜属成为并系群。