论文部分内容阅读
直翅目属于节肢动物门、六足总纲、昆虫纲,是起源比较早的类群之一,已知有2万多种,广泛分布于世界各地,以热带地区种类较多。直翅目由蝗亚目(Caelifera)和螽亚目(Ensifera)两个亚目组成。目前在GenBank中虽然已经积累了47个直翅目物种的全线粒体基因组,但是主要集中于蝗亚目的蝗总科物种,螽亚目物种相对较少。迄今为止,利用比较基因组手段来探讨直翅目线粒体基因组的结构及进化特征的研究甚少;采用不同的数据集划分策略来探讨线粒体DNA(mtDNA)基因在直翅目系统研究中的作用的研究也很少,且大多研究者集中于直翅目内高级阶元的系统关系即两个亚目是否为单系群,但是利用mtDNA基因数据集对直翅目内中级阶元和低级阶元系统发育的研究甚少。直翅目中癞蝗科和槌角蝗科均属于蝗亚目物种。癞蝗科中疙蝗属同短鼻蝗属在形态上的差别微乎甚微,且均分布于我国的荒漠草原地区;而槌角蝗科大足蝗属内西伯利亚蝗同李氏大足蝗的形态上也仅仅依据前后肘脉是否合并来判定,且二者在地域分布上有交叉。这种形态相似且地域相近的物种在线粒体基因组水平上是否也很保守?本研究测定了蝗总科中3种癞蝗和1种槌角蝗的线粒体基因组全序列,并对其进行了注释和分析。在此基础上,联合GenBank中已测直翅目线粒体基因组全序列,采用生物信息学同比较基因组学相结合的方法分析了直翅目内不同阶元物种mtDNA基因的演化模式。在此基础上,基于线粒体蛋白基因密码子3位点及不同种类蛋白和rRNA不同区域构建的线粒体基因数据集对直翅目内不同分类阶元的系统发育关系进行了研究。主要研究结果如下:1.西伯利亚蝗(Gomphocerus sibiricus)、红胫波腿蝗(Asiotmethis zacharjini)、贺兰短鼻蝗(Filchnerella helanshanensis)和红缘疙蝗(Pseudotmethis rubimarginis)的mtDNA总长度分别为15590bp、15660bp、15657bp和15661bp,均编码昆虫线粒体基因组中典型的37个基因,即13个蛋白编码基因,2个rRNA基因和22个tRNA基因。4种蝗虫的mtDNA基因排序及转录方向同蝗亚目中已测物种的一致。2.从整体mtDNA、22个tRNA基因和2个rRNA基因序列计算大足蝗属内西伯利亚蝗同李氏大足蝗之间的P-距离最小,但从13个蛋白基因计算,西伯利亚蝗和西藏大足蝗间的P-距离最小。无论从核苷酸和氨基酸水平变异率还是从非同义/同义置换率的比率来分析,线粒体蛋白基因cox1、cox2、cox3和nad4L在槌角蝗科内最为保守,进化最慢,受到的净化选择功能约束最少,而atp6、atp8和nad6进化速率最快,受到自然选择压的影响较少。在4个槌角蝗trnThr的TΨC臂中第3对碱基后存在2个碱基A的插入,此特征是否为槌角蝗科的共源性状还有待证明。rrnS在西伯利亚蝗和李氏大足蝗中几乎完全一致。然而在同属的西藏大足蝗物种中rrnS结构域Ⅰ和Ⅱ却存在较大的变异区域。3.无论从线粒体全基因组整体序列还是从13个蛋白基因、22个tRNA基因或2个rRNA基因计算,4个癞蝗物种中,贺兰短鼻蝗同红缘疙蝗的P-距离均为最近。线粒体蛋白基因cox1、cox2、cytb和atp6在癞蝗科内最为保守,而nad5和nad3的进化速率相对比较快。4.在槌角蝗科和癞蝗科两个类群中,nad6和atp8都显示了非常高的AT%,在癞蝗中nad4L的A+T含量也非常高,而3个cox基因的A+T含量都为最低。除ATP8为中度的A-偏斜外,分布于J链的每个线粒体蛋白基因几乎不存在AT-偏斜性。然而由N链编码的4个蛋白基因都有显著的T-偏斜值。在4个槌角蝗和4个癞蝗物种的A+T丰富区中均找到了茎环结构,但是该二级结构在红拟棒角蝗中同其他3个槌角蝗物种中差别较大;宽纹蠢蝗同其他3个癞蝗的茎环结构差异也比较大。而且在癞蝗物种中发现了位于茎环J链序列内部的T簇。5.昆虫纲所选6个目物种的A+T含量集中在75%~80%之间;AT-skew集中于0~0.05之间,GC-skew基本为负值(C>G),且集中于-0.3~0之间。其中,直翅目物种的A+T含量和AT-skew分布均比较分散,但二者似乎存在正相关性。直翅目碱基组成异质性分析表明,A+T含量相似、AT百分比接近平均值或亲缘关系较近的物种间ID值较小;线粒体蛋白不同密码子位点的差异指数由低到高顺序为:第2位点<第1位点<第3位点。6.线粒体13个蛋白编码基因中,cox1和cox3在直翅目大多数物种中的A+T含量普遍偏低,而atp8, nad4L和nad6的A+T含量基本偏高。在碱基偏斜方面,N-链蛋白基因的AT偏斜均为显著的T偏斜。J-链蛋白基因在螽亚目物种中基本为T偏斜,而在蝗亚目多数物种中为A偏斜。对GC偏斜而言,N链蛋白基因均为正值,而J链蛋白基因则均为负值。基于JC和K2P两种模型对直翅目内两个亚目13个线粒体蛋白基因的平均遗传距离进行分析后发现,在两个亚目中cox2变异程度都为最小,在蝗亚目中nad5变异最大,而在螽亚目中nad1变异最大。13个蛋白基因无论从核苷酸序列还是氨基酸序列的差异性进行分析,cytb基因变异程度都为最小,其次为3个cox基因,而nad6基因变异度最大。通常,氨基酸序列比核酸序列应该更加保守。然而,在13个线粒体蛋白基因中,仅在cox1基因中发现氨基酸序列相似度高于核酸序列,多数蛋白基因的氨基酸序列差异性高于或等于核酸序列,在nad6和atp8中尤为明显。7.22种线粒体tRNA基因在蝗总科各科的保守位点比例表明tRNA基因的核酸保守性有链间偏向性。16种位于J-链的tRNAs在蝗总科中核酸一致度均超过50%,另外,核酸保守度很高的trnLeuUUR, trnAsn和trnLys基因均位于J-链上。tRNA家族在癞蝗科内最为保守,其次为槌角蝗科。rrnS的3个结构域中,结构域Ⅰ在去除5’端后最为保守,其次为结构域Ⅲ,结构域Ⅱ进化速率最快。在rrnL的6个结构域中,结构域Ⅳ和Ⅴ最为保守,而结构域Ⅰ和Ⅱ进化最快。8. PCG1、PCG2、COX、COX+cytb和rRNA (C)这类保守基因或区段主要适用于分析直翅目内高级阶元的系统发育关系,如亚目及总科阶元。进化中等的PCG23及NADH基因则适用于从亚科至亚目任意一个阶元的系统关系研究,且建树方法也至关重要。rRNA (V)这种可变位点较多的区段比较适用于对亚科级阶元的系统关系研究。但是像ATP这种进化速率超快且总长度比较短的基因似乎不太适用于系统关系研究。