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.     

Comparative analysis of codon usage patterns in chloroplast genomes of the Asteraceae family

Codon usage bias (CUB) is an important evolutionary feature in a genome and has been widely documented from prokaryotes to eukaryotes. However, the significance of CUB in the Asteraceae family has not been well understood, with no Asteraceae species having been analyzed for this characteristic. Here, we use bioinformatics approaches to comparatively analyze the general patterns and influencing factors of CUB in five Asteraceae chloroplast (cp) genomes. The results indicated that the five genomes had similar codon usage patterns, showing a strong bias towards a high representation of NNA and NNT codons. Neutrality analysis showed that these cp genomes had a narrow GC distribution and no significant correlation was observed between GC₁₂ and GC₃. Parity Rule 2 (PR2) plot analysis revealed that purines were used more frequently than pyrimidines. Effective number of codons (ENc)-plot analysis showed that most genes followed the parabolic line of trajectory, but several genes with low ENc values lying below the expected curve were also observed. Furthermore, correspondence analysis of relative synonymous codon usage (RSCU) yielded a first axis that explained only a partial amount of variation of codon usage. These findings suggested that both natural selection and mutational bias contributed to codon bias, while selection was the major force to shape the codon usage in these Asteraceae cp genomes. Our study, which is the first to investigate codon usage patterns in Asteraceae plastomes, will provide helpful information about codon distribution and variation in these species, and also shed light on the genetic and evolutionary mechanisms of codon biology within this family.     

Factors affecting mito-nuclear codon usage interactions in the OXPHOS system of Drosophila melanogaster

Codon usage bias varies considerably among genomes and even within the genes of the same genome.In eukaryotic organisms,energy production in the form of oxidative phosphorylation(OXPHOS)is the only process under control of both nuclear and mitochondrial genomes.Although factors affecting codon usage in a single genome have been studied,this has not occurred when both interactional genomes are involved.Consequently, we investigated whether or not other factors influence codon usage of coevolved genes.We used Drosophila melanogaster as a model organism.Our χ2 test on the number of codons of nuclear and mitochondrial genes involved in the OXPHOS system was significantly different (χ2=7945.16,P<0.01).A plot of effective number of codons against GC3s content of nuclear genes showed that few genes lie on the expected curve,indicating that codon usage was random.Correspondence analysis indicated a significant correlation between axis 1 and codon adaptation index(R=0.947,P<0.01)in every nuclear gene sequence.Thus,codon usage bias of nuclear genes appeared to be affected by translational selection.Correlation between axis 1 coordinates and GC content(R=0.814.P<0.01)indicated that the codon usage of nuclear genes was also affected by GC composition.Analysis of mitochondrial genes did not reveal a significant correlation between axis 1 and any parameter.Statistical analyses indicated that codon usages of both nDNA and mtDNA were subjected to context-dependent mutations.     

Patterns of Genomic Integration of Nuclear Chloroplast DNA Fragments in Plant Species

The transfer of organelle DNA fragments to the nuclear genome is frequently observed in eukaryotes. These transfers are thought to play an important role in gene and genome evolution of eukaryotes. In plants, such transfers occur from plastid to nuclear [nuclear plastid DNAs (NUPTs)] and mitochondrial to nuclear (nuclear mitochondrial DNAs) genomes. The amount and genomic organization of organelle DNA fragments have been studied in model plant species, such as Arabidopsis thaliana and rice. At present, publicly available genomic data can be used to conduct such studies in non-model plants. In this study, we analysed the amount and genomic organization of NUPTs in 17 plant species for which genome sequences are available. The amount and distribution of NUPTs varied among the species. We also estimated the distribution of NUPTs according to the time of integration (relative age) by conducting sequence similarity analysis between NUPTs and the plastid genome. The age distributions suggested that the present genomic constitutions of NUPTs could be explained by the combination of the rapidly eliminated deleterious parts and few but constantly existing less deleterious parts.     

普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻(Oryza rufipogon Griff.)线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻(O. sativa L.)的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示:普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数(Nc)分布于45.32～61.00之间,其密码子偏性较弱; Nc值仅与GC_3呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数(CBI)和最优密码子使用频率(Fop)的相关性均达到显著水平;而GC_3和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度(RSCU)、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻(O. sativa L. subsp. japonica Kato)和籼稻(O. sativa L. subsp.indica Kato)内具有同质性。     

Evolution of Plant Mitochondrial Intron-Encoded Maturases: Frequent Lineage-Specific Loss and Recurrent Intracellular Transfer to the Nucleus

Among land plants, mitochondrial and plastid group II introns occasionally encode proteins called maturases that are important for splicing. Angiosperm nuclear genomes also encode maturases that are targeted to the organelles, but it is not known whether nucleus-encoded maturases exist in other land plant lineages. To examine the evolutionary diversity and history of this essential gene family, we searched for maturase homologs in recently sequenced nuclear and mitochondrial genomes from diverse land plants. We found that maturase content in mitochondrial genomes is highly lineage specific, such that orthologous maturases are rarely shared among major land plant groups. The presence of numerous mitochondrial pseudogenes in the mitochondrial genomes of several species implies that the sporadic maturase distribution is due to frequent inactivation and eventual loss over time. We also identified multiple maturase paralogs in the nuclear genomes of the lycophyte Selaginella moellendorffii, the moss Physcomitrella patens, and the representative angiosperm Vitis vinifera. Phylogenetic analyses of organelle- and nucleus-encoded maturases revealed that the nuclear maturase genes in angiosperms, lycophytes, and mosses arose by multiple shared and independent transfers of mitochondrial paralogs to the nuclear genome during land plant evolution. These findings indicate that plant mitochondrial maturases have experienced a surprisingly dynamic history due to a complex interaction of multiple evolutionary forces that affect the rates of maturase gain, retention, and loss.     

普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻（Oryza rufipogon Griff.）线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻（O.sativa L.）的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示：普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数（Nc）分布于45.32~61.00之间,其密码子偏性较弱;Nc值仅与GC₃呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数（CBI）和最优密码子使用频率（Fop）的相关性均达到显著水平;而GC₃和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度（RSCU）、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻（O.sativa L.subsp.japonica Kato）和籼稻（O.sativa L.subsp.indica Kato）内具有同质性。     

Inheritance of organelle DNA sequences in a citrus-poncirus intergeneric cross

Many land plants deviate from the maternal pattern of organelle inheritance. In this study, heterologous mitochondrial and chloroplast probes were used to investigate the inheritance of organelle genomes in the progeny of an intergeneric cross. The seed parent was LB 1-18 (a hybrid of Citrus reticulata Blanco cv. Clementine x C. paradisi Macf. cv. Duncan) and the pollen parent was the cross-compatible species Poncirus trifoliata (L.) Raf. All 26 progeny examined exhibited maternal inheritance of plastid petA and petD loci. However, 17 of the 26 progeny exhibited an apparent biparental inheritance of mitochondrial atpA, cob, coxII, and coxIII restriction fragment length polymorphisms (RFLPs) and maternal inheritance of mitochondrial rrn26 and coxI RFLPs. The remaining nine progeny inherited only maternal mitochondrial DNA (mtDNA) configurations. Investigations of plant mitochondrial genome inheritance are complicated by the multipartite structure of this genome, nuclear gene control over mitochondrial genome organization, and transfer of mitochondrial sequences to the nucleus. In this study, paternal mtDNA configurations were not detected in purified mtDNA of progeny plants, but were present in progeny DNA preparations enriched for nuclear genome sequences. MtDNA sequences in the nuclear genome therefore produced an inheritance pattern that mimics biparental inheritance of mtDNA.     

Comparison of codon usage and tRNAs in mitochondrial genomes of Candida species

To gain insight into the nature of the mitochondrial genomes (mtDNA) of different Candida species, the synonymous codon usage bias of mitochondrial protein coding genes and the tRNAs in C. albicans, C. parapsilosis, C. stellata, C. glabrata and the closely related yeast Saccharomyces cerevisiae were analyzed. Common features of the mtDNA in Candida species are a strong A+T pressure on protein coding genes, and insufficient mitochondrial tRNA species are encoded to perform protein synthesis. The wobble site of the anticodon is always U for the NNR (NNA and NNG) codon families, which are dominated by A-ending codons, and always G for the NNY (NNC and NNU) codon families, which is dominated by U-ending codons, and always U for the NNN (NNA, NNU, NNC and NNG) codon families, which are dominated by A-ending codons and U-ending codons. Patterns of synonymous codon usage of Candida species can be classified into three groups: (1) optimal codon-anticodon usage, Glu, Lys, Leu (translated by anti-codon UAA), Gln, Arg (translated by anti-codon UCU) and Trp are containing NNR codons. NNA, whose corresponding tRNA is encoded in the mtDNA, is used preferentially. (2) Non-optimal codon-anticodon usage, Cys, Asp, Phe, His, Asn, Ser (translated by anti-codon GCU) and Tyr are containing NNY codons. The NNU codon, whose corresponding tRNA is not encoded in the mtDNA, is used preferentially. (3) Combined codon-anticodon usage, Ala, Gly, Leu (translated by anti-codon UAG), Pro, Ser (translated by anti-codon UGA), Thr and Val are containing NNN codons. NNA (tRNA encoded in the mtDNA) and NNU (tRNA not encoded in the mtDNA) are used preferentially. In conclusion, we propose that in Candida species, codons containing A or U at third position are used preferentially, regardless of whether corresponding tRNAs are encoded in the mtDNA. These results might be useful in understanding the common features of the mtDNA in Candida species and patterns of synonymous codon usage.     

Is GC bias in the nuclear genome of the carnivorous plant Utricularia?driven by ROS-based mutation and biased gene conversion?

At less than 90 Mbp, the tiny nuclear genome of the carnivorous bladderwort plant Utricularia is an attractive model system for studying molecular evolutionary processes leading to genome miniaturization. Recently, we reported that expression of genes encoding DNA repair and reactive oxygen species (ROS) detoxification enzymes is highest in Utricularia traps, and we argued that ROS mutagenic action correlates with the high nucleotide substitution rates observed in the Utricularia plastid, mitochondrial, and nuclear genomes. Here, we extend our analysis of 100 nuclear genes from Utricularia and related asterid eudicots to examine nucleotide substitution biases and their potential correlation with ROS-induced DNA lesions. We discovered an unusual bias toward GC nucleotides, most prominently in transition substitutions at the third position of codons, which are presumably silent with respect to adaptation. Given the general tendency of biased gene conversion to drive GC bias, and of ROS to induce double strand breaks requiring recombinational repair, we propose that some of the unusual features of the bladderwort and its genome may be more reflective of these nonadaptive processes than of natural selection.     

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.     

Nuclear–cytoplasmic balance: whole genome duplications induce elevated organellar genome copy number

The plant genome is partitioned across three distinct subcellular compartments: the nucleus, mitochondria, and plastids. Successful coordination of gene expression among these organellar genomes and the nuclear genome is critical for plant function and fitness. Whole genome duplication (WGD) events in the nucleus have played a major role in the diversification of land plants and are expected to perturb the relative copy number (stoichiometry) of nuclear, mitochondrial, and plastid genomes. Thus, elucidating the mechanisms whereby plant cells respond to the cytonuclear stoichiometric imbalance that follows WGDs represents an important yet underexplored question in understanding the evolutionary consequences of genome doubling. We used droplet digital PCR to investigate the relationship between nuclear and organellar genome copy numbers in allopolyploids and their diploid progenitors in both wheat and Arabidopsis. Polyploids exhibit elevated organellar genome copy numbers per cell, largely preserving the cytonuclear stoichiometry observed in diploids despite the change in nuclear genome copy number. To investigate the timescale over which cytonuclear stoichiometry may respond to WGD, we also estimated the organellar genome copy number in Arabidopsis synthetic autopolyploids and in a haploid-induced diploid line. We observed corresponding changes in organellar genome copy number in these laboratory-generated lines, indicating that at least some of the cellular response to cytonuclear stoichiometric imbalance is immediate following WGD. We conclude that increases in organellar genome copy numbers represent a common response to polyploidization, suggesting that maintenance of cytonuclear stoichiometry is an important component in establishing polyploid lineages.     

Codon Usage in Plastid Genes Is Correlated with Context, Position Within the Gene, and Amino Acid Content

Highly expressed plastid genes display codon adaptation, which is defined as a bias toward a set of codons which are complementary to abundant tRNAs. This type of adaptation is similar to what is observed in highly expressed Escherichia coli genes and is probably the result of selection to increase translation efficiency. In the current work, the codon adaptation of plastid genes is studied with regard to three specific features that have been observed in E. coli and which may influence translation efficiency. These features are (1) a relatively low codon adaptation at the 5′ end of highly expressed genes, (2) an influence of neighboring codons on codon usage at a particular site (codon context), and (3) a correlation between the level of codon adaptation of a gene and its amino acid content. All three features are found in plastid genes. First, highly expressed plastid genes have a noticeable decrease in codon adaptation over the first 10–20 codons. Second, for the twofold degenerate NNY codon groups, highly expressed genes have an overall bias toward the NNC codon, but this is not observed when the 3′ neighboring base is a G. At these sites highly expressed genes are biased toward NNT instead of NNC. Third, plastid genes that have higher codon adaptations also tend to have an increased usage of amino acids with a high G + C content at the first two codon positions and GNN codons in particular. The correlation between codon adaptation and amino acid content exists separately for both cytosolic and membrane proteins and is not related to any obvious functional property. It is suggested that at certain sites selection discriminates between nonsynonymous codons based on translational, not functional, differences, with the result that the amino acid sequence of highly expressed proteins is partially influenced by selection for increased translation efficiency. Received: 21 July 1999 / Accepted: 5 November 1999     

MutS HOMOLOG1 is a nucleoid protein that alters mitochondrial and plastid properties and plant response to high light

Mitochondrial-plastid interdependence within the plant cell is presumed to be essential, but measurable demonstration of this intimate interaction is difficult. At the level of cellular metabolism, several biosynthetic pathways involve both mitochondrial- and plastid-localized steps. However, at an environmental response level, it is not clear how the two organelles intersect in programmed cellular responses. Here, we provide evidence, using genetic perturbation of the MutS Homolog1 (MSH1) nuclear gene in five plant species, that MSH1 functions within the mitochondrion and plastid to influence organellar genome behavior and plant growth patterns. The mitochondrial form of the protein participates in DNA recombination surveillance, with disruption of the gene resulting in enhanced mitochondrial genome recombination at numerous repeated sequences. The plastid-localized form of the protein interacts with the plastid genome and influences genome stability and plastid development, with its disruption leading to variegation of the plant. These developmental changes include altered patterns of nuclear gene expression. Consistency of plastid and mitochondrial response across both monocot and dicot species indicate that the dual-functioning nature of MSH1 is well conserved. Variegated tissues show changes in redox status together with enhanced plant survival and reproduction under photooxidative light conditions, evidence that the plastid changes triggered in this study comprise an adaptive response to naturally occurring light stress.     

Evolution in Hawaiian cave-adapted isopods (Oniscidea: Philosciidae): vicariant speciation or adaptive shifts?

Diatom plastid genes are examined with respect to codon adaptation and rates of silent substitution (Ks). It is shown that diatom genes follow the same pattern of codon usage as other plastid genes studied previously. Highly expressed diatom genes display codon adaptation, or a bias toward specific major codons, and these major codons are the same as those in red algae, green algae, and land plants. It is also found that there is a strong correlation between Ks and variation in codon adaptation across diatom genes, providing the first evidence for such a relationship in the algae. It is argued that this finding supports the notion that the correlation arises from selective constraints, not from variation in mutation rate among genes. Finally, the diatom genes are examined with respect to variation in Ks among different synonymous groups. Diatom genes with strong codon adaptation do not show the same variation in synonymous substitution rate among codon groups as the flowering plant psbA gene which, previous studies have shown, has strong codon adaptation but unusually high rates of silent change in certain synonymous groups. The lack of a similar finding in diatoms supports the suggestion that the feature is unique to the flowering plant psbA due to recent relaxations in selective pressure in that lineage.     

Codon usage in Entamoeba histolytica.

The codon usage of 10 E. histolytica genes comprising 4455 codons was analysed. The codon usage revealed an extremely biased use of synonymous codons with a preference for NNU (44%) and NNA (41.4%) codons. Codons CGG (arg), AGG (arg) and CCG (pro) were absent in the E. histolytica genes examined. The codon usage of E. histolytica resembled that of Plasmodium falciparum.     

Evolutionary analysis of aspartate aminotransferases

Aspartate aminotransferase isoenzymes are located in both the cytosol and organelles of eukaryotes, but all are encoded in the nuclear genome. In the work described here, a phylogenetic analysis was made of aspartate aminotransferases from plants, animals, yeast, and a number of bacteria. This analysis suggested that five distinct branches are present in the aspartate aminotransferase tree. Mitochondrial forms of the enzyme form one distinct group, bacterial aspartate aminotransferase formed another, and the plant and vertebrate cytosolic isoenzymes each formed a distinct group. Plant cytosolic isozymes formed a further group of which the plastid sequences were a member. The yeast mitochondrial and cytosolic aspartate aminotransferases formed groups separate from other members of the family. Correspondence to: C.J. Marshall     

The mitochondrial genome of Arabidopsis is composed of both native and immigrant information

Plants contain large mitochondrial genomes, which are several times as complex as those in animals, fungi or algae. However, genome size is not correlated with information content. The mitochondrial genome (mtDNA) of Arabidopsis specifies only 58 genes in 367 kb, whereas the 184 kb mtDNA in the liverwort Marchantia polymorpha codes for 66 genes, and the 58 kb genome in the green alga Prototheca wickerhamii encodes 63 genes. In Arabidopsis’ mtDNA, genes for subunits of complex II, for several ribosomal proteins and for 16 tRNAs are missing, some of which have been transferred recently to the nuclear genome. Numerous integrated fragments originate from alien genomes, including 16 sequence stretches of plastid origin, 41 fragments of nuclear (retro)transposons and two fragments of fungal viruses. These immigrant sequences suggest that the large size of plant mitochondrial genomes is caused by secondary expansion as a result of integration and propagation, and is thus a derived trait established during the evolution of land plants.