Comparative studies on codon usage pattern of chloroplasts and their host nuclear genes in four plant species

A detailed comparison was made of codon usage of chloroplast genes with their host (nuclear) genes in the four angiosperm speciesOryza sativa, Zea mays, Triticum aestivum andArabidopsis thaliana. The average GC content of the entire genes, and at the three codon positions individually, was higher in nuclear than in chloroplast genes, suggesting different genomic organization and mutation pressures in nuclear and chloroplast genes. The results of Nc-plots and neutrality plots suggested that nucleotide compositional constraint had a large contribution to codon usage bias of nuclear genes inO. sativa, Z. mays, andT. aestivum, whereas natural selection was likely to be playing a large role in codon usage bias in chloroplast genomes. Correspondence analysis and chi-test showed that regardless of the genomic environment (species) of the host, the codon usage pattern of chloroplast genes differed from nuclear genes of their host species by their AU-richness. All the chloroplast genomes have predominantly A- and/or U-ending codons, whereas nuclear genomes have G-, C- or U-ending codons as their optimal codons. These findings suggest that the chloroplast genome might display particular characteristics of codon usage that are different from its host nuclear genome. However, one feature common to both chloroplast and nuclear genomes in this study was that pyrimidines were found more frequently than purines at the synonymous codon position of optimal codons.     

Chloroplast DNA codon use: Evidence for selection at the psb A locus based on tRNA availability

Codon use in the three sequenced chloroplast genomes (Marchantia, Oryza, and Nicotiana) is examined. The chloroplast has a bias in that codons NNA and NNT are favored over synonymous NNC and NNG codons. This appears to be a consequence of an overall high A + T content of the genome. This pattern of codon use is not followed by the psb A gene of all three genomes and other psb A sequences examined. In this gene, the codon use favors NNC over NNT for twofold degenerate amino acids. In each case the only tRNA coded by the genome is complementary to the NNC codon. This codon use is similar to the codon use by chloroplast genes examined from Chlamydomonas reinhardtii. Since psb A is the major translation product of the chloroplast, this suggests that selection is acting on the codon use of this gene to adapt codons to tRNA availability, as previously suggested for unicellular organisms.     

Selection on the codon bias of chloroplast and cyanelle genes in different plant and algal lineages

In the plant chloroplast genome the codon usage of the highly expressed psbA gene is unique and is adapted to the tRNA population, probably due to selection for translation efficiency. In this study the role of selection on codon usage in each of the fully sequenced chloroplast genomes, in addition to Chlamydomonas reinhardtii, is investigated by measuring adaptation to this pattern of codon usage. A method is developed which tests selection on each gene individually by constructing sequences with the same amino acid composition as the gene and randomly assigning codons based on the nucleotide composition of noncoding regions of that genome. The codon bias of the actual gene is then compared to a distribution of random sequences. The data indicate that within the algae selection is strong in Cyanophora paradoxa, affecting a majority of genes, of intermediate intensity in Odontella sinensis, and weaker in Porphyra purpurea and Euglena gracilis. In the plants, selection is found to be quite weak in Pinus thunbergii and the angiosperms but there is evidence that an intermediate level of selection exists in the liverwort Marchantia polymorpha. The role of selection is then further investigated in two comparative studies. It is shown that average relative codon bias is correlated with expression level and that, despite saturation levels of substitution, there is a strong correlation among the algae genomes in the degree of codon bias of homologous genes. All of these data indicate that selection for translation efficiency plays a significant role in determining the codon bias of chloroplast genes but that it acts with different intensities in different lineages. In general it is stronger in the algae than the higher plants, but within the algae Euglena is found to have several unusual features which are noted. The factors that might be responsible for this variation in intensity among the various genomes are discussed. Received: 6 June 1997 / Accepted: 24 July 1997     

滇黄精叶绿体全基因组序列及其密码子使用偏性分析

为探究滇黄精(Polygonatum kingianum)叶绿体全基因组特征和密码子使用偏性,利用第二代测序技术对滇黄精嫩叶进行测序,再经组装与注释后得到其叶绿体基因组全序列,通过MISA、EMBOSS和CodonW等软件对滇黄精叶绿体全基因组的SSR位点、系统发育及密码子偏好性进行分析。结果表明,滇黄精完整叶绿体基因组长度为155 852 bp,基因组平均GC含量为37.7%,其大、小单拷贝区(LSC)长度分别为84 633和185 25 bp,反向重复区长度为26 347 bp,注释了132个基因,包括86个蛋白编码基因、38个tRNA基因和8个核糖rRNA基因。叶绿体基因组中共有69个SSR位点,绝大多数属于单碱基重复的A/T类型。系统发育分析表明滇黄精与格脉黄精(P. tessellatum)亲缘关系近,可能与分布地域有关。密码子偏好性分析表明,滇黄精叶绿体基因组密码子使用模式受到自然选择影响大于突变因素,最终确定9个最优密码子。因此, 滇黄精叶绿体基因组遗传结构和系统发育位置及其密码子偏倚的分析,为叶绿体基因工程研究提供理论依据。     

Complete mitochondrial genome of the boreal digging frog Kaloula borealis (Anura, Microhylidae)

We sequenced the complete mitochondrial genome from the boreal digging frog Kaloula borealis. The genome sequence was 17,173 bp in size, and the gene order and contents were identical to those of previously reported amphibian mitochondrial genomes. Of 13 protein-coding genes (PCGs), 5 genes (CO2, ATPase 6, CO3, ND3, and ND4) had incomplete stop codons. Also ND1 gene used GTG as a start codon, while CO1 and ND5 genes used AGG as a stop codon. The base composition of K. borealis mitogenome showed a strong anti-G bias (6.11%) on the 3rd position of PCGs.     

Selection on the codon bias ofChlamydomonas reinhardtii chloroplast genes and the plantpsbA gene

Plant chloroplast genes have a codon use that reflects the genome compositional bias of a high A+T content with the single exception of the highly translatedpsbA gene which codes for the photosystem II D1 protein. The codon usage of plantpsbA corresponds more closely to the limited tRNA population of the chloroplast and is very similar to the codon use observed in the chloroplast genes of the green algaChlamydomonas reinhardtii. This pattern of codon use may be an adaptation for increased translation efficiency. A correspondence between codon use of plantpsbA andChlamydomonas chloroplast genes and the tRNAs coded by the chloroplast genome, however, is not observed in all synonymous codon groups. It is shown here that the degree of correspondence between codon use and tRNA population in different synonymous groups is correlated with the second codon position composition. Synonymous groups with an A or T at the second codon position have a high representation of codons for which a complementary tRNA is coded by the chloroplast genome. Those with a G or C at the second position have an increased representation of codons that bind a chloroplast tRNA by wobble. It is proposed that the difference between synonymous groups in terms of codon adaptation to the tRNA population in plantpsbA andChlamydomonas chloroplast genes may be the result of differences in second position composition.     

Extensive Rearrangements in the Chloroplast Genome of <Emphasis Type="Italic">Trachelium caeruleum</Emphasis> Are Associated with Repeats and tRNA Genes

Chloroplast genome organization, gene order, and content are highly conserved among land plants. We sequenced the chloroplast genome of Trachelium caeruleum L. (Campanulaceae), a member of an angiosperm family known for highly rearranged genomes. The total genome size is 162,321 bp, with an inverted repeat (IR) of 27,273 bp, large single-copy (LSC) region of 100,114 bp, and small single-copy (SSC) region of 7,661 bp. The genome encodes 112 different genes, with 17 duplicated in the IR, a tRNA gene (trnI-cau) duplicated once in the LSC region, and a protein-coding gene (psbJ) with two duplicate copies, for a total of 132 putatively intact genes. ndhK may be a pseudogene with internal stop codons, and clpP, ycf1, and ycf2 are so highly diverged that they also may be pseudogenes. ycf15, rpl23, infA, and accD are truncated and likely nonfunctional. The most conspicuous feature of the Trachelium genome is the presence of 18 internally unrearranged blocks of genes inverted or relocated within the genome relative to the ancestral gene order of angiosperm chloroplast genomes. Recombination between repeats or tRNA genes has been suggested as a mechanism of chloroplast genome rearrangements. The Trachelium chloroplast genome shares with Pelargonium and Jasminum both a higher number of repeats and larger repeated sequences in comparison to eight other angiosperm chloroplast genomes, and these are concentrated near rearrangement endpoints. Genes for tRNAs occur at many but not all inversion endpoints, so some combination of repeats and tRNA genes may have mediated these rearrangements.     

药用植物华重楼(黑药花科)叶绿体全基因组研究

为探究华重楼(Paris polyphylla var. chinensis)的叶绿体基因组特征,利用叶绿体系统发育基因组学方法,对华重楼与其它百合目植物的叶绿体全基因组进行了比较。结果表明,华重楼的叶绿体全基因组长158307 bp,由4个区组成,包括2个反向重复区(IRA和IRB,27473 bp)、1个小单拷贝区(SSC,18175 bp)和1个大单拷贝区(LSC,85187 bp)。其叶绿体基因组有115个基因,包括81个编码蛋白质基因、30个转运RNA基因和4 个核糖体RNA基因。11种百合目植物的叶绿体全基因组的基因组成和基因顺序相似。华重楼的cemA基因是假基因,其起始密码子后有多聚核苷酸poly(A)及CA双核苷酸重复序列,编码序列中出现多个终止密码子, 且与北重楼(Paris verticillata)的cemA编码序列中的终止密码子位置不同。因此,华重楼叶绿体基因组比较保守;cemA结构及假基因化现象可能具有重要的进化与系统发育信息,其编码序列中的终止密码子可以区分华重楼和北重楼。     

Differential fate of plastid and mitochondrial genomes in Petunia somatic hybrids

Summary The chloroplast (cp) and mitochondrial (mt) DNAs of Petunia somatic hybrid plants, which were derived from the fusion of wild-type P. parodii protoplasts with albino P. inflata protoplasts, were analyzed by endonuclease restriction and Southern blot hybridization. Using ³²P-labelled probes that distinguished the two parental cpDNAs at a BamH1 site and at a HpaII site, only the P. parodii chloroplast genome was detected in the 10 somatic hybrid plants analyzed. To examine whether cytoplasmic mixing had resulted in rearrangement of the mitochondrial genome in the somatic hybrids, restriction patterns of purified somatic hybrid and parental mtDNAs were analyzed. Approximately 87% of those restriction fragments which distinguish the two parental genomes are P. inflata-specific. Restriction patterns of the somatic hybrid mtDNAs differ both from the parental patterns and from each other, suggesting that an interaction occurred between the parental mitochondrial genomes in the somatic fusion products which resulted in generation of the novel mtDNA patterns. Southern blot hybridization substantiates this conclusion. In addition, somatic hybrid lines derived from the same fusion product were observed to differ in mtDNA restriction pattern, reflecting a differential sorting-out of mitochondrial genomes at the time the plants were regenerated.     

Mitochondrial Genome of Prasinophyte Alga Pyramimonas parkeae

Prasinophytes are a paraphyletic assemblage of nine heterogeneous lineages in the Chlorophyta clade of Archaeplastida. Until now, seven complete mitochondrial genomes have been sequenced from four prasinophyte lineages. Here, we report the mitochondrial genome of Pyramimonas parkeae, the first representative of the prasinophyte clade I. The circular‐mapping molecule is 43,294 bp long, AT rich (68.8%), very compact and it comprises two 6,671 bp long inverted repeat regions. The gene content is slightly smaller than the gene‐richest prasinophyte mitochondrial genomes. The single identified intron is located in the cytochrome c oxidase subunit 1 gene (cox1). Interestingly, two exons of cox1 are encoded on the same strand of DNA in the reverse order and the mature mRNA is formed by trans‐splicing. The phylogenetic analysis using the data set of 6,037 positions assembled from 34 mtDNA‐encoded proteins of 48 green algae and plants is not in compliance with the branching order of prasinophyte clades revealed on the basis of 18S rRNA genes and cpDNA‐encoded proteins. However, the phylogenetic analyses based on all three genomic elements support the sister position of prasinophyte clades Pyramimonadales and Mamiellales.     

巨枝叶绿体基因组密码子偏好性分析

该文针对巨桉叶绿体基因组序列,选取其中长于300 nt且以AUG为起始密码子的43个非重复基因作为研究对象,采用CodonW1.4.2软件分析巨桉叶绿体基因组的密码子使用偏好性。结果表明:第3位密码子的平均GC含量为27.97%; ENC的变化范围为39.49～61.00,平均为47.04; RSCU1的密码子有31个,其中29个以A/U结尾;中性分析显示,GC12与GC3无显著相关;回归分析未达到显著性水平; ENCplot分析发现,大部分基因落在曲线上或附近;对应分析表明第1轴的贡献率为17.68%,第2轴的贡献率为11.49%,第3轴、第4轴的贡献率分别为8.00%和5.76%,前4轴累计贡献率达42.93%,第1轴与GC、ENC、CAI达到极显著相关。上述分析结果表明,巨桉叶绿体基因组的密码子偏好较弱,密码子第3位偏好以A或U结尾,选择和突变在巨桉叶绿体基因组密码子偏好中起相对均衡的作用,最终确定UUG、CUU、GUU、UCC、UCA、ACA、UAU、UAA、CAU、AAU、AGA和GGA 12个高频高表达密码子为最优密码子。这为转化叶绿体基因密码子优化,提高表达效率和改良巨桉目标性状奠定了坚实基础。     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

沙枣叶绿体全基因组序列及其使用密码子偏性分析

该研究以2株野生沙枣(Elaeagnus angustifolia Linn.)嫩枝经温室水培后的嫩叶为材料,采用CTAB法分别提取总DNA,并利用第二代测序技术进行总DNA从头测序,组装后得到2株沙枣叶绿体基因组全序列,并详细分析了其蛋白质编码基因密码子使用的偏好性及其原因,为沙枣叶绿体基因工程和分子系统进化等研究奠定基础。结果显示:(1)组装得到沙枣叶绿体基因组序列全长150 546 bp,由长度为81 113 bp的长单拷贝(LSC)区域和25 494 bp的短单拷贝(SSC)区域,以及1对分隔开它们的长18 445 bp的反向重复序列(IRS)组成;注释共得到132个基因,包括86个蛋白编码基因、38个tRNA基因和8个rRNA基因。(2)沙枣叶绿体基因组蛋白编码基因密码子的第三位碱基GC含量(GC_3)为28.47%,明显低于整个叶绿体基因组GC含量(37%),也低于第一位(GC_1)和第二位(GC_2)碱基的GC含量,说明密码子对AT碱基结尾有偏好性;其中, UCU、CCU、UGU、GCU、CUU、GAU、UCA和UAA为最优密码子。(3)同义密码子相对使用频率(RSCU)分析发现,影响密码子使用模式的因素并不单一,密码子的偏好性受到突变、选择及其他因素的共同影响,并且自然选择表达引起的序列差异比突变对密码子偏好性的影响要显著;中性绘图分析、有效密码子数(ENC-plot)分析和奇偶偏好性(PR2-plot)分析表明,沙枣叶绿体基因组使用密码子的偏性受选择的影响更大。(4)通过最大似然法、最大简约法和贝叶斯方法对胡颓子科6个物种和1个枣的叶绿体基因序列构建系统发育树,与它们使用密码子偏性聚类的结果一致,表明叶绿体基因组使用密码子偏性与物种的亲缘关系相关。     

Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer

The mitochondrial genome of grape (Vitis vinifera), the largestorganelle genome sequenced so far, is presented. The genomeis 773,279 nt long and has the highest coding capacity amongknown angiosperm mitochondrial DNAs (mtDNAs). The proportionof promiscuous DNA of plastid origin in the genome is also thelargest ever reported for an angiosperm mtDNA, both in absoluteand relative terms. In all, 42.4% of chloroplast genome of Vitishas been incorporated into its mitochondrial genome. In orderto test if horizontal gene transfer (HGT) has also contributedto the gene content of the grape mtDNA, we built phylogenetictrees with the coding sequences of mitochondrial genes of grapeand their homologs from plant mitochondrial genomes. Many incongruentgene tree topologies were obtained. However, the extent of incongruencebetween these gene trees is not significantly greater than thatobserved among optimal trees for chloroplast genes, the commonancestry of which has never been in doubt. In both cases, weattribute this incongruence to artifacts of tree reconstruction,insufficient numbers of characters, and gene paralogy. Thisfinding leads us to question the recent phylogenetic interpretationof Bergthorsson et al. (2003, 2004) and Richardson and Palmer(2007) that rampant HGT into the mtDNA of Amborella best explainsphylogenetic incongruence between mitochondrial gene trees forangiosperms. The only evidence for HGT into the Vitis mtDNAfound involves fragments of two coding sequences stemming fromtwo closteroviruses that cause the leaf roll disease of thisplant. We also report that analysis of sequences shared by bothchloroplast and mitochondrial genomes provides evidence fora previously unknown gene transfer route from the mitochondrionto the chloroplast.     

Complete nucleotide sequence of the chloroplast genome from a leptosporangiate fern, Adiantum capillus-veneris L.

We determined the complete nucleotide sequence of the chloroplast genome of the leptosporangiate fern, Adiantum capillus-veneris L. (Pteridaceae). The circular genome is 150,568 bp, with a large single-copy region (LSC) of 82,282 bp, a small-single copy region (SSC) of 21,392 bp and inverted repeats (IR) of 23,447 bp each. We compared the sequence to other published chloroplast genomes to infer the location of putative genes. When the IR is considered only once, we assigned 118 genes, of which 85 encode proteins, 29 encode tRNAs and 4 encode rRNAs. Four protein-coding genes, all four rRNA genes and six tRNA genes occur in the IR. Most (57) putative protein-coding genes appear to start with an ATG codon, but we also detected five other possible start codons, some of which suggest tRNA editing. We also found 26 apparent stop codons in 18 putative genes, also suggestive of RNA editing. We found all but one of the tRNA genes necessary to encode the complete repertoire required for translation. The missing trnK gene appears to have been disrupted by a large inversion, relative to other published chloroplast genomes. We detected several structural rearrangements that may provide useful information for phylogenetic studies.     

Analysis of synonymous codon usage patterns in different plant mitochondrial genomes

Codon usage in mitochondrial genome of the six different plants was analyzed to find general patterns of codon usage in plant mitochondrial genomes. The neutrality analysis indicated that the codon usage patterns of mitochondrial genes were more conserved in GC content and no correlation between GC12 and GC3. T and A ending codons were detected as the preferred codons in plant mitochondrial genomes. The Parity Rule 2 plot analysis showed that T was used more frequently than A. The EN_C-plot showed that although a majority of the points with low EN_C values were lying below the expected curve, a few genes lied on the expected curve. Correspondence analysis of relative synonymous codon usage yielded a first axis that explained only a partial amount of variation of codon usage. These findings suggest that natural selection is likely to be playing a large role in codon usage bias in plant mitochondrial genomes, but not only natural selection but also other several factors are likely to be involved in determining the selective constraints on codon bias in plant mitochondrial genomes. Meantime, 1 codon (P. patens), 6 codons (Z. mays), 9 codons (T. aestivum), 15 codons (A. thaliana), 15 codons (M. polymorpha) and 15 codons (N. tabacum) were defined as the preferred codons of the six plant mitochondrial genomes.     

蔓赤车叶绿体基因组结构及系统进化分析

蔓赤车(Pellionia scabra)属于荨麻科赤车属,是一种具有高营养价值的优质野菜。本研究以蔓赤车为材料,基于高通量技术方法完成叶绿体基因组测序、组装注释、结构解析及构建系统发育树,以此深入研究蔓赤车叶绿体基因组特征。结果表明,蔓赤车叶绿体基因组大小为153220 bp,GC含量为36.4%,属于典型的四分体结构。共注释到130个基因,其中85个蛋白编码基因,37个转运RNA基因,8个核糖体RNA基因;其中,有15个基因包含1个内含子,2个基因包含2个内含子,rps12存在反式剪接情况。蔓赤车叶绿体基因组可分为光合作用(43个)、自我复制(64个)、其他编码蛋白(7个)以及未知功能(4个)4大类基因。蔓赤车叶绿体基因组共检测出51073个密码子,其中编码亮氨酸(Leu)的密码子占比最大,密码子偏向使用A和U两种碱基。检测到72个简单重复序列位点,其中有58个单核苷酸、12个二核苷酸、1个三核苷酸和1个四核苷酸4种不同类型的简单重复序列(simple sequence repeats,SSRs)。蔓赤车IRb/SSC边界存在ycf1基因扩张现象。系统进化树显示,蔓赤车(Pellionia scabra OL800583)与庐山楼梯草(Elatostema stewardii MZ292972)、盘托楼梯草(Elatostema dissectum MK227819)、光叶楼梯草(Elatostema laevissimum var.laevissimum MN189961)亲缘关系最密切。基于蔓赤车叶绿体基因组的研究,旨在为蔓赤车物种鉴定、遗传进化及基因组学研究奠定理论基础。     

The Complete Mitochondrial Genome of the Entomopathogenic Nematode Steinernema carpocapsae: Insights into Nematode Mitochondrial DNA Evolution and Phylogeny

We determined the complete sequence of the mitochondrial DNA of the entomopathogenic nematode Steinernema carpocapsae and analyzed its structure and composition as well as the secondary structures predicted for its tRNAs and rRNAs. Almost the complete genome has been amplified in one fragment with long PCR and sequenced using a shotgun strategy. The 13,925-bp genome contains genes for 2 rRNAs, 22 tRNAs, and 12 proteins and lacks an ORF encoding ATPase subunit 8. Four initiation codons were inferred, TTT, TTA, ATA, and ATT, most of the genes ended with TAA or TAG, and only two had a T as an incomplete stop codon. All predicted tRNAs showed the nonconventional secondary structure typical of Secernentea. Although we were able to fold the sequences of trnN, trnD, and trnC into more conventional cloverleaf structures after adding adjacent nucleotides, northern blot experiments showed that the nonstandard tRNAs are actually expressed. Phylogenetic and comparative analyses showed that the mitochondrial genome of S. carpocapsae is more closely related to the genomes of A. suum and C. elegans than to that of Strongyloides stercoralis. This finding does not support the phylogeny based on nuclear small subunit ribosomal DNA sequences previously published. This discrepancy may result from differential reproductive strategies and/or differential selective pressure acting on nuclear and mitochondrial genes. The distinctive characteristics observed among mitochondrial genomes of Secernentea may have arisen to counteract the deleterious effects of Muller’s ratchet, which is probably enhanced by the reproductive strategies and selective pressures referred to above. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Rafael Zardoya]     

两种梧桐叶绿体基因组密码子使用偏性分析

为了分析美丽梧桐、云南梧桐叶绿体基因组密码子的使用偏性,该研究通过筛选美丽梧桐、云南梧桐叶绿体基因组中各52条蛋白编码序列,并利用CodonW、CUSP和SPSS软件对其密码子使用模式及偏性进行了分析。结果表明:(1)美丽梧桐、云南梧桐的GC含量分别为38.12%、38.05%,表明叶绿体基因组内富含A/T碱基。(2)有效密码子数(ENC)范围为36.91～56.46、36.55～58.04,表明多数密码子偏性较弱。(3)相对同义密码子(RSCU)分析显示,RSCU1的密码子各有29个,其中28个以A、U结尾。(4)中性绘图显示,GC_3与GC_(12)的相关性不显著,回归曲线斜率分别为0.195和0.304,说明密码子偏好性主要受到自然选择的影响。(5) ENC-plot分析中大部分基因分布于曲线的周围和下方,ENC比值多分布于-0.04～0.10之间,表明突变会影响密码子偏性的形成。此外,17、18个密码子分别被鉴定为美丽梧桐、云南梧桐的最优密码子。以上结果说明美丽梧桐、云南梧桐叶绿体基因组的密码子使用偏性可能受选择和突变共同作用,且使用模式较为相似,但具有一定的差异,可能与适应环境的进化机制有关。     

桃叶珊瑚属植物的叶绿体基因组结构特征及系统发育分析

为确定桃叶珊瑚属(Aucuba)植物叶绿体基因组的结构及其序列变异,揭示其属下种间亲缘关系,该研究对桃叶珊瑚(A. chinensis)、花叶青木(A. japonica var. variegata)等6种桃叶珊瑚属植物和丝缨花属植物黄杨叶丝缨花(Garrya buxifolia)进行二代测序,利用生物信息学软件对其叶绿体基因组序列进行组装和注释,并进行基本特征分析、序列比较以及系统发育分析。结果表明:(1)桃叶珊瑚属植物叶绿体基因组具典型的环状四分体结构,6条序列全长157 891～158 325 bp,均编码114个基因,包括80个蛋白质编码基因、30个tRNA基因和4个rRNA基因。(2)6种植物叶绿体基因组高频密码子数均为29个,偏好以A/U结尾,确定了这6条序列的最优密码子共100个,包含12个共有的最优密码子。(3)6条叶绿体基因组序列共检测到270条散在重复序列,133条串联重复序列以及412个SSR位点。(4)比较基因组学分析结果表明,该属植物叶绿体基因组序列高度保守。(5)从叶绿体基因组中筛选出10个高变片段。(6)系统发育分析结果显示支持桃叶珊瑚属为一个支持率较高的单系,与丝缨花属关系较近。该研究中的5种桃叶珊瑚属植物以及1种丝缨花属植物的叶绿体基因组均为首次测序组装,揭示了桃叶珊瑚属及其属下种间的系统发育关系,为桃叶珊瑚属植物的分类鉴定和系统发育提供了参考资料。