利用cvTree方法构建及分析tRNA序列的亲缘关系进化网络

本文利用cvTree和复杂网络理论相结合的方法,构建了3 420条tRNA序列亲缘关系进化网络.计算了相关网络的平均度、平均聚类系数和平均最短路径随进化距离Dis变化的关系.同时还分析和讨论了tRNA序列之间的亲缘关系及其进化机制,观察到一个与其它方法不同的结果:随着亲缘距离Dis的增加,网络的度分布先是从power-law分布转变为Gussian分布,然后又从Gussian分布转变为反power-law分布.本研究结果说明由cvTree方法和复杂网络理论相结合的方法研究tRNA基因序列的进化的结果也许更加符合其进化历程.     

球孢白僵菌线粒体序列分析及系统进化树构建

线粒体是真核细胞的重要的细胞器,在细胞能量代谢和细胞衰老方面扮演了重要角色,同时也是研究物种进化关系的重要资源。本文应用球孢白僵菌的线粒体序列分析了白僵菌与粪壳菌纲其他11属的13种真菌间的进化关系,为阐明球孢白僵菌的进化地位提供了新的论据。线粒体编码的rnl,rns,25个tRNAs和14个蛋白质基因均由同一条链编码,大部分的tRNA分布为三个tRNA簇,这与粪壳菌纲的其他物种相似。与NCBI上公布的序列比对发现球孢白僵菌线粒体序列与蝇蚧霉线粒体序列最相似,相似率达73%,用14个蛋白编码基因通过邻接法(Neighbor-joint)和最大简约法(Maximum parsimony)构建的系统进化树也得到相同结论,其支持率均为100%。     

COP9信号传导体与调节盖子复合体的进化

COP9信号传导体和26S蛋白酶体的调节盖子复合体皆为含有8个亚基的蛋白复合体,在真核生物体中普遍存在,它们的相应亚基在大小和氨基酸序列上具有一一对应关系.从NCBI站点的所有数据库中获得了裂殖酵母、酿酒酵母、线虫、果蝇、哺乳动物和拟南芥等多种生物的复合体的亚基序列共8组.COP9信号传导体与调节盖子复合体相应亚基之间的氨基酸序列一致性大于12%,它们均具有一些保守的区域,而且保守位点分布均匀,表明它们来自于同一祖先.在基于氨基酸序列构建的系统发育树中,各组序列分别形成两个分支:一个分支由COP9信号传导体亚基和相似蛋白组成,另一分支由相应的调节盖子复合体亚基和相似蛋白构成.各个分支中单细胞生物的序列位于动、植物序列的根部,表明COP9信号传导体与调节盖子复合体的基因重复发生在真核单细胞生物和多细胞生物分化以前,并且二者的亚基基因沿各自的方向独立进化.几乎所有编码两个蛋白复合体的基因在基因组中均为单拷贝,第Ⅴ、Ⅵ组的亚基保守程度最高,暗示着它们在复合体中起着关键的作用.对COP9信号传导体和调节盖子复合体的相应亚基基因两两之间进行dN/dS的相关性分析,分别鉴定出21和15对亚基编码序列间具有显著的Pearson相关关系,推测其相应亚基间可能通过承担相互关联的重要的生物学功能而协同进化.     

COP9信号传导体与调节盖子复合体的进化

COP9信号传导体和26S蛋白酶体的调节盖子复合体皆为含有8个亚基的蛋白复合体,在真核生物体中普遍存在,它们的相应亚基在大小和氨基酸序列上具有一一对应关系。从NCBI站点的所有数据库中获得了裂殖酵母、酿酒酵母、线虫、果蝇、哺乳动物和拟南芥等多种生物的复合体的亚基序列共8组。COP9信号传导体与调节盖子复合体相应亚基之间的氨基酸序列一致性大于12％,它们均具有一些保守的区域,而且保守位点分布均匀,表明它们来自于同一祖先。在基于氨基酸序列构建的系统发育树中,各组序列分别形成两个分支：一个分支由COP9信号传导体亚基和相似蛋白组成,另一分支由相应的调节盖子复合体亚基和相似蛋白构成。各个分支中单细胞生物的序列位于动、植物序列的根部,表明COP9信号传导体与调节盖子复合体的基因重复发生在真核单细胞生物和多细胞生物分化以前,并且二者的亚基基因沿各自的方向独立进化。几乎所有编码两个蛋白复合体的基因在基因组中均为单拷贝,第Ⅴ、Ⅵ组的亚基保守程度最高,暗示着它们在复合体中起着关键的作用。对COP9信号传导体和调节盖子复合体的相应亚基基因两两之间进行dN／dS的相关性分析,分别鉴定出21和15对亚基编码序列间具有显著的Pearson相关关系,推测其相应亚基间可能通过承担相互关联的重要的生物学功能而协同进化。     

颈棱蛇蛇毒金属蛋白酶基因扩增及系统进化分析

[目的]扩增颈棱蛇蛇毒金属蛋白酶基因,并对其进行系统进化分析。[方法]通过PCR方法扩增颈棱蛇蛇毒金属蛋白酶基因,使用Lasergene Meg Align 8.1.3软件对多个蛇种的蛇毒金属蛋白酶序列进行比对,使用MEGA 6.06对颈棱蛇和其它蛇种的蛇毒金属蛋白酶氨基酸序列进行进化分析。[结果]颈棱蛇蛇毒金属蛋白酶编码区长1 854 bp,编码617个氨基酸残基。该酶的核酸序列与绿栖林蛇、银环蛇、舟山眼镜蛇、圆斑蝰蛇和尖吻蝮蛇的序列相似度分别为82.1%、79.2%、79.1%、75.2%和75.7%,氨基酸序列相似度分别为68.9%、63.3%、63.1%、60.3%和60.1%。进化分析显示,颈棱蛇和同属于游蛇科水游蛇亚科的波加丹游蛇亲缘关系最近,其次为游蛇科美洲钝头蛇亚科的绿栖林蛇,与眼镜蛇科和蝰科这两个科亲缘关系最远,与非洲蛇科的小鳞穴蝰亲缘关系介于游蛇科和眼镜蛇科之间。[结论]蛇毒金属蛋白酶可以作为蛇种进化分析的一种重要参考基因。     

小麂、赤麂、黑麂mtDNA序列变异性及反映的进化关系

利用已测定的鹿科麂亚科动物小麂、赤麂、黑麂的线粒体全基因组序列,统计它们各自连接在一起的13个蛋白编码基因、22个tRNA基因、2个rRNA基因和1个控制区序列的碱基长度和组成,计算rRNA基因遗传距离,估算分歧时间,比较蛋白编码基因的碱基水平和氨基酸水平上的差异,基于连接在一起的13个氨基酸序列,以羊为外群,通过邻位相连法和最大简约性法构建进化树,探讨小麂、赤麂、黑麂的进化关系。结果表明,小麂是较原始的物种,赤麂和黑麂较为近缘,是从类似小麂的祖先演化而来。     

第一类肽链释放因子研究进展

第一类释放因子是新生肽链释放所必需的因子,它能正确地识别终止信号,水解肽酰-tRNA酯键,释放出新合成的多肽链.它的高级结构与tRNA结构相似,从而解释了它们功能上的相似性,其氨基酸序列的高度保守区GGQ motif和三肽反密码子区在新生肽链的释放中分别执行重要的功能.     

中华大蟾蜍γ-肌动蛋白基因Bbgγ-Actin的克隆与生物信息学分析

旨在获得中华大蟾蜍肌动蛋白序列全长。以中华大蟾蜍(Bufo bufo gargarizans)耳后腺为材料,提取总RNA,通过转录组高通量测序快速得到中华大蟾蜍肌动蛋白基因(Bbg Actin)cDNA全长2 055 bp,通过RT-PCR技术对其ORF序列进行验证,并对Bbg Actin基因进行分析。序列分析表明,Bbg Actin基因开放阅读框长1 131 bp,编码376个氨基酸。通过BLAST P程序分析,所得序列与Gen Bank数据库中肌动蛋白基因序列相似度均在89%以上,氨基酸序列的相似性达90%以上。系统进化分析表明,Bbg Actin与γ-Actin聚为一类,所得中华大蟾蜍肌动蛋白属于γ-Actin。首次获得了中华大蟾蜍γ-Actin cDNA全长序列。     

BRCA1蛋白的序列及空间结构的分析

查询人的BRCA1蛋白的氨基酸序列,利用生物信息学的方法进行相似性搜索,获得一系列BRCA1蛋白的氨基酸序列。选择了其中的11条序列,对BRCA1蛋白进行了多重序列分析和进化分析,对BRCA1蛋白的BRCT结构域进行三维同源模型的构建与比较分析。分析结果表明:BRCA1中某些特定部位的氨基酸序列高度保守;确定氨基酸的保守位点并联合进化分析可对基因错义突变的致病性做初步地猜测;相近物种来源的BRCA1具有较近的亲缘关系,而且具有极其相似的三维空间结构。这些为研究BRCA1蛋白的结构与功能关系提供指导意义。     

美洲大蠊基因组重复序列分析

重复序列是真核生物基因组的重要组成部分。一些重复序列,如自主型的逆转录转座子LINE,在昆虫的系统进化和遗传多样性研究方面得到了广泛的应用。de novo从头预测和基于同源比对预测相结合的方法被用来搜索美洲大蠊Periplaneta americana基因组,共鉴定出大约占全基因组62%的重复序列。研究发现,散在重复序列中,DNA转座子占美洲大蠊基因组的16.18%;逆转座元件中LINE最多,占基因组的13.64%,SINE和LTR分别占基因组的3.52%和1.32%。LINEs中的Bov Bs亚家族在所有转座子亚家族中比例最高(约6.73%)。美洲大蠊与德国小蠊Blattella germanica相比,除LTR外,其他类型的转座子占基因组的比例均高于德国小蠊。通过分析逆转录转座子反转录酶完整度、氨基酸序列相似度及遗传距离,从美洲大蠊基因组中鉴定出一类BovBs:RTE-1_PAm。BovBs的反转录酶氨基酸序列的系统树表明,美洲大蠊与内华达古白蚁Zootermopsis nevadensis的进化关系比与其同属蜚蠊科Blattidae的德国小蠊的关系更近。昆虫中BovBs的进化关系与传统核基因进化关系的不同,表明转座子的进化相对宿主基因的进化具有一定的独立性。     

Studying the evolutionary relationships and phylogenetic trees of 21 groups of tRNA sequences based on complex networks

To find out the evolutionary relationships among different tRNA sequences of 21 amino acids, 22 networks are constructed. One is constructed from whole tRNAs, and the other 21 networks are constructed from the tRNAs which carry the same amino acids. A new method is proposed such that the alignment scores of any two amino acids groups are determined by the average degree and the average clustering coefficient of their networks. The anticodon feature of isolated tRNA and the phylogenetic trees of 21 group networks are discussed. We find that some isolated tRNA sequences in 21 networks still connect with other tRNAs outside their group, which reflects the fact that those tRNAs might evolve by intercrossing among these 21 groups. We also find that most anticodons among the same cluster are only one base different in the same sites when S ≥ 70, and they stay in the same rank in the ladder of evolutionary relationships. Those observations seem to agree on that some tRNAs might mutate from the same ancestor sequences based on point mutation mechanisms.     

Comparing two evolutionary mechanisms of modern tRNAs

All modern tRNA gene families have a high similarity in their primary structure, and share the same cloverleaf secondary structure and an inverted L tertiary structure, which provide the clues for the study of their origin and evolution. There are two important mechanisms of the tRNA sequences evolution. One is point mutation, another is complementary duplication method. Both of them are supported by some specific examples. To find out the superior one of the two mechanisms or find out the most suitable mechanism for modern tRNAs evolution, we constructed two types of networks, parallel and antiparallel networks, based on the two mechanisms respectively, and then compared the degree distribution and clustering coefficient of networks constructed by the tRNA sequences of the single anticodon group, single isoaccepting group, and the whole tRNAs group of the two types of networks. The result of the comparison seems consistent with the idea that modern tRNA sequences evolved primarily by the mechanism of complementary method, and point mutation is an important and indispensable auxiliary mechanism during the evolutionary event.     

氨酰-tRNA合成酶专一性识别的进化

氨酰-tRNA合成酶(AARS)是一类在蛋白质合成过程中起着重要作用的酶,它通过与tRNA及其相应氨基酸的专一性识别作用,使得基因序列能够被精确地翻译成蛋白质序列.然而,氨酰-tRNA合成酶的这种识别作用既有专一性,也具有“兼容性”.氨酰-tRNA合成酶的这种双重性质不仅与其结构的进化有关,而且还与其所处的各类生物的不同进化阶段有关.AARS似乎经历了一个由“模糊专一性”（多重专一性）到“精确专一性”（单一专一性）的演变历程.     

Sequence evidence in the archaeal genomes that tRNAs emerged through the combination of ancestral genes as 5' and 3' tRNA halves

The discovery of separate 5' and 3' halves of transfer RNA (tRNA) molecules-so-called split tRNA-in the archaeal parasite Nanoarchaeum equitans made us wonder whether ancestral tRNA was encoded on 1 or 2 genes. We performed a comprehensive phylogenetic analysis of tRNAs in 45 archaeal species to explore the relationship between the three types of tRNAs (nonintronic, intronic and split). We classified 1953 mature tRNA sequences into 22 clusters. All split tRNAs have shown phylogenetic relationships with other tRNAs possessing the same anticodon. We also mimicked split tRNA by artificially separating the tRNA sequences of 7 primitive archaeal species at the anticodon and analyzed the sequence similarity and diversity of the 5' and 3' tRNA halves. Network analysis revealed specific characteristics of and topological differences between the 5' and 3' tRNA halves: the 5' half sequences were categorized into 6 distinct groups with a sequence similarity of >80%, while the 3' half sequences were categorized into 9 groups with a higher sequence similarity of >88%, suggesting different evolutionary backgrounds of the 2 halves. Furthermore, the combinations of 5' and 3' halves corresponded with the variation of amino acids in the codon table. We found not only universally conserved combinations of 5'-3' tRNA halves in tRNA(iMet), tRNA(Thr), tRNA(Ile), tRNA(Gly), tRNA(Gln), tRNA(Glu), tRNA(Asp), tRNA(Lys), tRNA(Arg) and tRNA(Leu) but also phylum-specific combinations in tRNA(Pro), tRNA(Ala), and tRNA(Trp). Our results support the idea that tRNA emerged through the combination of separate genes and explain the sequence diversity that arose during archaeal tRNA evolution.     

The origin of the protein synthesis mechanism

The origin and development of the protein synthesis mechanism is considered in four successive steps. The genetic code is supposed to be controlled by the relative amount (availability) of various amino acids and nucleotides on the one hand, and utility on each amino acid in the polypeptide. on the other hand. Thus, more simple (inutile) and abundant amino acids tended to correspond to codons which were rich in the less frequent base species, G and C. Features of primitive tRNA in the discrimination of amino acid are discussed. Primitive tRNA is proposed to have a discriminator site for amino acid and, separated from it, an anticodon site for interaction with nucleotides. A hypothetical course of subdivision of various nucleic acid species is proposed. In the scheme, mRNA and ribosomal RNA (rRNA) were derived from more primitive insoluble RNA. DNA appeared in the late, not first, step of the development. Several other aspects of evolutionary development of the whole protein synthesis mechanism, e.g., role of the discriminator site on primitive tRNA, modification and subdivision of code catalogue into a more precise specification of amino acids, and possible primordial interactions between tRNA and tRNA-binding sites on insoluble rRNA, are discussed.     

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ‐stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNA^Tyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences “U‐U‐x2‐C” and “U‐x‐G‐x2‐T,” thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.     

Evolution of tRNA recognition systems and tRNA gene sequences

The aminoacylation of tRNAs by the aminoacyl-tRNA synthetases recapitulates the genetic code by dictating the association between amino acids and tRNA anticodons. The sequences of tRNAs were analyzed to investigate the nature of primordial recognition systems and to make inferences about the evolution of tRNA gene sequences and the evolution of the genetic code. Evidence is presented that primordial synthetases recognized acceptor stem nucleotides prior to the establishment of the three major phylogenetic lineages. However, acceptor stem sequences probably did not achieve a level of sequence diversity sufficient to faithfully specify the anticodon assignments of all 20 amino acids. This putative bottleneck in the evolution of the genetic code may have been alleviated by the advent of anticodon recognition. A phylogenetic analysis of tRNA gene sequences from the deep Archaea revealed groups that are united by sequence motifs which are located within a region of the tRNA that is involved in determining its tertiary structure. An association between the third anticodon nucleotide (N36) and these sequence motifs suggests that a tRNA-like structure existed close to the time that amino acid-anticodon assignments were being established. The sequence analysis also revealed that tRNA genes may evolve by anticodon mutations that recruit tRNAs from one isoaccepting group to another. Thus tRNA gene evolution may not always be monophyletic with respect to each isoaccepting group.Based on a presentation made at a workshop— Aminoacyl-tRNA Synthetases and the Evolution of the Genetic Code—held at Berkeley, CA, July 17–20, 1994 Correspondence to: M.E. Saks     

Aminoacyl-tRNA synthetase classes and groups in prokaryotes

Knowledge on the evolution of aminoacyl-tRNA synthetases is crucial to studies on the origins of life. The relationships between the different aminoacyl-tRNA synthetase specificities in prokaryotic organisms are studied in this work. We reconstructed the ancestor sequences and the phylogenetic relationships utilizing the Maximum Likelihood method. The results suggest that in class I the evolution of the N-terminal segment was strongly influenced by the amino acid hydropathy in both domains of prokaryotes. The results for the C-terminal segments of class I were different in the two domains, indicating that its evolution was strongly influenced by the specific types of tRNA modification in each domain. The class II groups in Archaea were more heterogeneous with respect to the hydropathy of amino acids, indicating the interference of other influences. In bacteria, the configuration was also complex but the overall consensual division in two groups was maintained, group IIa forming a single branch with the five hydroapathetic amino acid specificities and group IIb containing the specificities for the moderately hydrophobic together with the hydrophilic amino acids. It is indicated that the aminoacyl-tRNA synthetase in both domains were subjected to different selective forces in diverse parts of the proteins, resulting in complex phylogenetic patterns.     

Revisiting the operational RNA code for amino acids: Ensemble attributes and their implications

It has been suggested that tRNA acceptor stems specify an operational RNA code for amino acids. In the last 20 years several attributes of the putative code have been elucidated for a small number of model organisms. To gain insight about the ensemble attributes of the code, we analyzed 4925 tRNA sequences from 102 bacterial and 21 archaeal species. Here, we used a classification and regression tree (CART) methodology, and we found that the degrees of degeneracy or specificity of the RNA codes in both Archaea and Bacteria differ from those of the genetic code. We found instances of taxon-specific alternative codes, i.e., identical acceptor stem determinants encrypting different amino acids in different species, as well as instances of ambiguity, i.e., identical acceptor stem determinants encrypting two or more amino acids in the same species. When partitioning the data by class of synthetase, the degree of code ambiguity was significantly reduced. In cryptographic terms, a plausible interpretation of this result is that the class distinction in synthetases is an essential part of the decryption rules for resolving the subset of RNA code ambiguities enciphered by identical acceptor stem determinants of tRNAs acylated by enzymes belonging to the two classes. In evolutionary terms, our findings lend support to the notion that in the pre-DNA world, interactions between tRNA acceptor stems and synthetases formed the basis for the distinction between the two classes; hence, ambiguities in the ancient RNA code were pivotal for the fixation of these enzymes in the genomes of ancestral prokaryotes.     

An informational approach to recognizing the similarity of tRNA genes and their global classification]

Based on informational approach to the comparative analysis of DNA sequences and the conception of fuzzy relations, new classification schemes were suggested for two aggregates of tRNA genes from the EMBL-8 data bank connecting amino acids TRP and ASP. These classification schemes are presented as a dendrogram and reflect possible evolutionary interrelation between tRNA genes in both aggregates. Evolutionary interrelations between tRNA genes were analyzed from two standpoints: 1) taking into consideration the interrelation between DNA sequences by all four nucleotides; 2) taking into account possible evolutionary pressure to the divergence of DNA sequences towards transitions. The similarity and difference of dendrograms for each gene aggregate were considered in both instances. It was shown that the analysis of all possible directions of the evolutionary pressure to the divergence of DNA and RNA sequences provides a new approach to the construction of a global tRNA phylogeny.