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.     

Nucleotide Content Gradients in Maternally and Paternally Inherited Mitochondrial Genomes of the Mussel <Emphasis Type="Italic">Mytilus</Emphasis>

Several studies have shown that in vertebrate mtDNAs the nucleotide content at fourfold degenerate sites is well correlated with the site’s time of exposure to the single-strand state, as predicted from the asymmetrical model of mtDNA replication. Here we examine whether the same explanation may hold for the regional variation in nucleotide content in the maternal and paternal mtDNAs of the mussel Mytilus galloprovincialis. The origin of replication of the heavy strand (O_H) of these genomes has been previously established. A systematic search of the two genomes for sequences that are likely to act as the origin of replication of the light strand (O_L) suggested that the most probable site lies within the ND3 gene. By adopting this O_L position we calculated times of exposure for 0_FD (nondegenerate), 2_FD (twofold degenerate), and 4_FD (fourfold degenerate) sites of the protein-coding part of the genome and for the rRNA, tRNA and noncoding parts. The presence of thymine and absence of guanine at 4_FD sites was highly correlated with the presumed time of exposure. Such an effect was not found for the 2_FD sites, the rRNA, the tRNA, or the noncoding parts. There was a trend for a small increase in cytosine at 0_FD sites with exposure time, which is explicable as the result of biased usage of 4_FD codons. The same analysis was applied to a recently sequenced mitochondrial genome of Mytilus trossulus and produced similar results. These results are consistent with the asymmetrical model of replication and suggest that guanine oxidation due to single-strand exposure is the main cause of regional variation of nucleotide content in Mytilus mitochondrial genomes. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Dr. David Pollock]     

Analysis of mitochondrial protein‐coding genes of Antheraea assamensis: Muga silkworm of Assam

To understand the synonymous codon usage pattern in mitochondrial genome of Antheraea assamensis, we analyzed the 13 mitochondrial protein‐coding genes of this species using a bioinformatic approach as no work was reported yet. The nucleotide composition analysis suggested that the percentages of A, T, G,and C were 33.73, 46.39, 9.7 and 10.17, respectively and the overall GC content was 19.86, that is, lower than 50% and the genes were AT rich. The mean effective number of codons of mitochondrial protein‐coding genes was 36.30 and it indicated low codon usage bias (CUB). Relative synonymous codon usage analysis suggested overrepresented and underrepresented codons in each gene and the pattern of codon usage was different among genes. Neutrality plot analysis revealed a narrow range of distribution for GC content at the third codon position and some points were diagonally distributed, suggesting both mutation pressure and natural selection influenced the CUB.     

Across bacterial phyla, distantly-related genomes with similar genomic GC content have similar patterns of amino acid usage

The GC content of bacterial genomes ranges from 16% to 75% and wide ranges of genomic GC content are observed within many bacterial phyla, including both gram negative and gram positive phyla. Thus, divergent genomic GC content has evolved repeatedly in widely separated bacterial taxa. Since genomic GC content influences codon usage, we examined codon usage patterns and predicted protein amino acid content as a function of genomic GC content within eight different phyla or classes of bacteria. We found that similar patterns of codon usage and protein amino acid content have evolved independently in all eight groups of bacteria. For example, in each group, use of amino acids encoded by GC-rich codons increased by approximately 1% for each 10% increase in genomic GC content, while the use of amino acids encoded by AT-rich codons decreased by a similar amount. This consistency within every phylum and class studied led us to conclude that GC content appears to be the primary determinant of the codon and amino acid usage patterns observed in bacterial genomes. These results also indicate that selection for translational efficiency of highly expressed genes is constrained by the genomic parameters associated with the GC content of the host genome.     

The mitochondrial genome of <Emphasis Type="Italic">Protohermes concolorus</Emphasis> Yang et Yang 1988 (Insecta: Megaloptera: Corydalidae)

The first complete mitochondrial genome of dobsonfly Protohermes concolorus Yang et Yang, 1988 (Megaloptera: Corydalidae) was sequenced in this study. The genome was a circular molecule of 15,851 bp containing the typical 37 genes that arranged in the same order as that of the putative ancestor of hexapods. Sequences overlaps were observed between several neighbor genes, which made the genome relatively compact. The tRNA-Ser (GCT) could not be folded into typical secondary structure because its DHU arm was replaced with a simple loop. Six of the 13 protein genes were terminated with a single T adjacent to a downstream tRNA gene in the same strand. The variation of GC content caused the different nucleotide substitution patterns of the protein genes. The genome was AT-biased with a total A + T content of 75.83% which was also demonstrated by the codon usage. The control region was the most AT-rich region with a sub-region of even higher A + T content. Protein genes of two strands presented opposite CG-skew trends which was also reflected by the codon usage. For most of the amino acids, the protein coding sequences did not prefer to use the cognate codons of corresponding tRNAs and the codon usage of the protein genes was not random. The variation of nucleotide substitution patterns of protein genes was significantly correlated with the GC content. The phylogenetic analyses based on all the 13 protein genes showed that Megaloptera was the sister group of other holometabolous insects except Coleoptera.     

Using the taxon-specific genes for the taxonomic classification of bacterial genomes

Background

The correct taxonomic assignment of bacterial genomes is a primary and challenging task. With the availability of whole genome sequences, the gene content based approaches appear promising in inferring the bacterial taxonomy. The complete genome sequencing of a bacterial genome often reveals a substantial number of unique genes present only in that genome which can be used for its taxonomic classification.

Results

In this study, we have proposed a comprehensive method which uses the taxon-specific genes for the correct taxonomic assignment of existing and new bacterial genomes. The taxon-specific genes identified at each taxonomic rank have been successfully used for the taxonomic classification of 2,342 genomes present in the NCBI genomes, 36 newly sequenced genomes, and 17 genomes for which the complete taxonomy is not yet known. This approach has been implemented for the development of a tool ‘Microtaxi’ which can be used for the taxonomic assignment of complete bacterial genomes.

Conclusion

The taxon-specific gene based approach provides an alternate valuable methodology to carry out the taxonomic classification of newly sequenced or existing bacterial genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1542-0) contains supplementary material, which is available to authorized users.     

双翅目昆虫线粒体基因组结构特点及其测序通用引物的设计和应用

研究双翅目昆虫线粒体基因组的结构特点, 并设计其测序的通用引物, 为今后双翅目昆虫线粒体基因组的研究提供参考和依据。利用比较基因组学和生物信息学方法, 分析了已经完全测序的26个双翅目昆虫线粒体基因组的结构特点、 碱基组成和保守区, 并据此设计了双翅目昆虫基因组测序的通用引物。结果表明： 双翅目昆虫线粒体基因组长14 503~19 517 bp, 其结构保守, 含有37个编码基因, 包括13个蛋白质编码基因, 22个tRNA编码基因和2个rRNA编码基因, 此外还包含一段长度差异很大的非编码区(AT富含区)。基因组内基因排列次序稳定, 除个别基因外, 其余都与黑腹果蝇Drosophila melanogaster基因排列次序一致。基因组的碱基组成不均衡, AT含量在72.59%~85.15%之间, 碱基使用存在偏向性, 偏好使用AC碱基。全基因组的核苷酸和氨基酸序列保守, 共鉴定了11个保守区。在保守区内共设计了26对双翅目线粒体基因组测序通用引物, 扩增的目标片段都在1 200 bp以内。将该套通用引物用于葱蝇Delia antiqua线粒体全基因组测序, 结果证明其高效、 合用。     

Nonneutral GC3 and Retroelement Codon Mimicry in Phytophthora

Phytophthora is a genus entirely comprised of destructive plant pathogens. It belongs to the Stramenopila, a unique branch of eukaryotes, phylogenetically distinct from plants, animals, or fungi. Phytophthora genes show a strong preference for usage of codons ending with G or C (high GC3). The presence of high GC3 in genes can be utilized to differentiate coding regions from noncoding regions in the genome. We found that both selective pressure and mutation bias drive codon bias in Phytophthora. Indicative for selection pressure is the higher GC3 value of highly expressed genes in different Phytophthora species. Lineage specific GC increase of noncoding regions is reminiscent of whole-genome mutation bias, whereas the elevated Phytophthora GC3 is primarily a result of translation efficiency-driven selection. Heterogeneous retrotransposons exist in Phytophthora genomes and many of them vary in their GC content. Interestingly, the most widespread groups of retroelements in Phytophthora show high GC3 and a codon bias that is similar to host genes. Apparently, selection pressure has been exerted on the retroelement’s codon usage, and such mimicry of host codon bias might be beneficial for the propagation of retrotransposons. Reviewing Editor: Dr. Yves van de Peer     

The complete mitochondrial genome of the yellow coaster, <Emphasis Type="Italic">Acraea</Emphasis> <Emphasis Type="Italic">issoria</Emphasis> (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini): sequence,gene organization and a unique tRNA translocation event

In this paper, the complete mitochondrial genome of Acraea issoria (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini) is reported; a circular molecule of 15,245 bp in size. For A. issoria, genes are arranged in the same order and orientation as the complete sequenced mitochondrial genomes of the other lepidopteran species, except for the presence of an extra copy of tRNA^Ile(AUR)b in the control region. All protein-coding genes of A. issoria mitogenome start with a typical ATN codon and terminate in the common stop codon TAA, except that COI gene uses TTG as its initial codon and terminates in a single T residue. All tRNA genes possess the typical clover leaf secondary structure except for tRNA^Ser(AGN), which has a simple loop with the absence of the DHU stem. The sequence, organization and other features including nucleotide composition and codon usage of this mitochondrial genome were also reported and compared with those of other sequenced lepidopterans mitochondrial genomes. There are some short microsatellite-like repeat regions (e.g., (TA)₉, polyA and polyT) scattered in the control region, however, the conspicuous macro-repeats units commonly found in other insect species are absent.     

Unparalleled GC content in the plastid DNA of Selaginella

One of the more conspicuous features of plastid DNA (ptDNA) is its low guanine and cytosine (GC) content. As of February 2009, all completely-sequenced plastid genomes have a GC content below 43% except for the ptDNA of the lycophyte Selaginella uncinata, which is 55% GC. The forces driving the S. uncinata ptDNA towards G and C are undetermined, and it is unknown if other Selaginella species have GC-biased plastid genomes. This study presents the complete ptDNA sequence of Selaginella moellendorffii and compares it with the previously reported S. uncinata plastid genome. Partial ptDNA sequences from 103 different Selaginella species are also described as well as a significant proportion of the S. moellendorffii mitochondrial genome. Moreover, S. moellendorffii express sequence tags are data-mined to estimate levels of plastid and mitochondrial RNA editing. Overall, these data are used to show that: (1) there is a genus-wide GC bias in Selaginella ptDNA, which is most pronounced in South American articulate species; (2) within the Lycopsida class (and among plants in general), GC-biased ptDNA is restricted to the Selaginella genus; (3) the cause of this GC bias is arguably a combination of reduced AT-mutation pressure relative to other plastid genomes and a large number of C-to-U RNA editing sites; and (4) the mitochondrial DNA (mtDNA) of S. moellendorffii is also GC biased (even more so than the ptDNA) and is arguably the most GC-rich organelle genome observed to date—the high GC content of the mtDNA also appears to be influenced by RNA editing. Ultimately, these findings provide convincing support for the earlier proposed theory that the GC content of land-plant organelle DNA is positively correlated and directly connected to levels of organelle RNA editing.     

Presence of two mitochondrial genomes in the mytilid Perumytilus purpuratus: Phylogenetic evidence for doubly uniparental inheritance

This study presents evidence, using sequences of ribosomal 16S and COI mtDNA, for the presence of two mitochondrial genomes in Perumytilus purpuratus. This may be considered evidence of doubly uniparental mtDNA inheritance. The presence of the two types of mitochondrial genomes differentiates females from males. The F genome was found in the somatic and gonadal tissues of females and in the somatic tissues of males; the M genome was found in the gonads and mantle of males only. For the mitochondrial 16S region, ten haplotypes were found for the F genome (nucleotide diversity 0.004), and 7 haplotypes for the M genome (nucleotide diversity 0.001), with a distance Dxy of 0.125 and divergence Kxy of 60.33%. For the COI gene 17 haplotypes were found for the F genome (nucleotide diversity 0.009), and 10 haplotypes for the M genome (nucleotide diversity 0.010), with a genetic distance Dxy of 0.184 and divergence Kxy of 99.97%. Our results report the presence of two well-differentiated, sex-specific types of mitochondrial genome (one present in the male gonad, the other in the female gonad), implying the presence of DUI in P. purpuratus. These results indicate that care must be taken in phylogenetic comparisons using mtDNA sequences of P. purpuratus without considering the sex of the individuals.     

Complete mitochondrial genome of the <Emphasis Type="Italic">Pleuronichthys lighti</Emphasis> (Pleuronectiformes,Pleuronectidae) with phylogenetic consideration

Mitochondrial genome has been used to shed light on many fields of both basic and applied research, including the study of molecular evolution. The complete mitochondrial genome sequence of 17368 bp nucleotides from the Pleuronichthys lighti was determined. It was a circular double-stranded DNA molecule with identical set of 22 transfer RNA genes, 2 ribosomal RNA genes, 13 protein-coding genes as well as a non-coding control region. Stand asymmetry in the nucleotide composition was reflected in the codon usage of genes oriented in opposite directions. In the control region, we identified the extended termination associated sequence domain, the central conserved sequence block domain and the conserved sequence block domain, and two complete repeat region. They were “TTACAATA” and “TGTTGTAA”, respectively. All known 12 mitochondrial genomes of Pleuronectinae fishes were downloaded and analyzed; there were 5570 variable sites in the consensus sequences of 15241 base pairs, calculation of total sites were 35.5%. The highest sequence divergence was 50% (ATP8) and the Kimura-2-parameter genetic distance was 0.235 (ND6), whereas the COIII had the lowest sequence divergence (28.8%) and genetic distance (0.128); the protein coding genes were mainly acted by purifying selection which was detected by selection tests. Analysis of confidence and the information content for per nucleotide revealed ND5, ATP6, COI and ND4 genes were suitable molecular markers for phylogenetic study of Pleuronectinae fishes. Phylogenetic analysis using Bayesian computational algorithms based on COI genes provided support for the taxonomic status of P. lighti, which was consistent with the traditional taxonomy.     

桃儿七叶绿体比较基因组学分析

为探究桃儿七(Sinopodophyllum hexandrum)不同叶绿体基因组特征,本研究以桃儿七5个叶绿体基因组为研究对象,借助生物信息学工具进行基因组图谱构建、重复序列分析、密码子偏好分析、反向重复序列区(inverted repeat, IR)/单拷贝区(single-copy, SC)边界分析、基因组序列比较分析及系统发育分析。结果表明:桃儿七5个叶绿体基因组全长为157 203–157 940 bp,为典型的叶绿体四分体结构,共注释出133–137个基因,说明桃儿七叶绿体基因组具有多样性。桃儿七不同叶绿体基因组简单重复序列(simple sequence repeat, SSR)位点不同,单核苷酸A/T占主要优势,散在重复序列包括正向重复、回文重复和反向重复3类。密码子偏好分析显示有效密码子(effective number of codon,ENc)值为51.14–51.17,密码子偏好性弱,GC与GC3s所占比例小于50%,密码子偏向使用A和U碱基并且以A和U碱基结尾。桃儿七5个叶绿体基因组IR/SC边界和基因组序列均比较保守。系统发育分析结果表明桃儿七和北美桃儿七亲...     

Rapid divergence of codon usage patterns within the rice genome

Background

Synonymous codon usage varies widely between genomes, and also between genes within genomes. Although there is now a large body of data on variations in codon usage, it is still not clear if the observed patterns reflect the effects of positive Darwinian selection acting at the level of translational efficiency or whether these patterns are due simply to the effects of mutational bias. In this study, we have included both intra-genomic and inter-genomic comparisons of codon usage. This allows us to distinguish more efficiently between the effects of nucleotide bias and translational selection.

Results

We show that there is an extreme degree of heterogeneity in codon usage patterns within the rice genome, and that this heterogeneity is highly correlated with differences in nucleotide content (particularly GC content) between the genes. In contrast to the situation observed within the rice genome, Arabidopsis genes show relatively little variation in both codon usage and nucleotide content. By exploiting a combination of intra-genomic and inter-genomic comparisons, we provide evidence that the differences in codon usage among the rice genes reflect a relatively rapid evolutionary increase in the GC content of some rice genes. We also noted that the degree of codon bias was negatively correlated with gene length.

Conclusion

Our results show that mutational bias can cause a dramatic evolutionary divergence in codon usage patterns within a period of approximately two hundred million years.The heterogeneity of codon usage patterns within the rice genome can be explained by a balance between genome-wide mutational biases and negative selection against these biased mutations. The strength of the negative selection is proportional to the length of the coding sequences. Our results indicate that the large variations in synonymous codon usage are not related to selection acting on the translational efficiency of synonymous codons.

     

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.     

How Mitochondria Redefine the Code

Annotated, complete DNA sequences are available for 213 mitochondrial genomes from 132 species. These provide an extensive sample of evolutionary adjustment of codon usage and meaning spanning the history of this organelle. Because most known coding changes are mitochondrial, such data bear on the general mechanism of codon reassignment. Coding changes have been attributed variously to loss of codons due to changes in directional mutation affecting the genome GC content (Osawa and Jukes 1988), to pressure to reduce the number of mitochondrial tRNAs to minimize the genome size (Anderson and Kurland 1991), and to the existence of transitional coding mechanisms in which translation is ambiguous (Schultz and Yarus 1994a). We find that a succession of such steps explains existing reassignments well. In particular, (1) Genomic variation in the prevalence of a codon's third-position nucleotide predicts relative mitochondrial codon usage well, though GC content does not. This is because A and T, and G and C, are uncorrelated in mitochondrial genomes. (2) Codons predicted to reach zero usage (disappear) do so more often than expected by chance, and codons that do disappear are disproportionately likely to be reassigned. However, codons predicted to disappear are not significantly more likely to be reassigned. Therefore, low codon frequencies can be related to codon reassignment, but appear to be neither necessary nor sufficient for reassignment. (3) Changes in the genetic code are not more likely to accompany smaller numbers of tRNA genes and are not more frequent in smaller genomes. Thus, mitochondrial codons are not reassigned during demonstrable selection for decreased genome size. Instead, the data suggest that both codon disappearance and codon reassignment depend on at least one other event. This mitochondrial event (leading to reassignment) occurs more frequently when a codon has disappeared, and produces only a small subset of possible reassignments. We suggest that coding ambiguity, the extension of a tRNA's decoding capacity beyond its original set of codons, is the second event. Ambiguity can act alone but often acts in concert with codon disappearance, which promotes codon reassignment. Received: 26 October 2000 / Accepted: 19 January 2001     

Mitochondrial cytochrome C oxidase subunit I of Manduca sexta and a comparison with other invertebrate genes

A cDNA encoding mitochondrial cytochrome c oxidase subunit I (mt COI) from Manduca sexta (Lepidoptera: Sphingidae) was cloned and sequenced. AT (adenine-thymine) content is high and codon usage is biased and likely reflects the role of mt COI in electron transport. The encoded protein is 514 amino acids long, contains seven invariant His residues observed in COIs in all organisms and would be predicted to be composed of 12 transmembrane regions.     

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.     

Comparison of the Yeast Proteome to Other Fungal Genomes to Find Core Fungal Genes

The purpose of this research was to search for evolutionarily conserved fungal sequences to test the hypothesis that fungi have a set of core genes that are not found in other organisms, as these genes may indicate what makes fungi different from other organisms. By comparing 6355 predicted or known yeast (Saccharomyces cerevisiae) genes to the genomes of 13 other fungi using Standalone TBLASTN at an e-value <1E-5, a list of 3340 yeast genes was obtained with homologs present in at least 12 of 14 fungal genomes. By comparing these common fungal genes to complete genomes of animals (Fugu rubripes, Caenorhabditis elegans), plants (Arabidopsis thaliana, Oryza sativa), and bacteria (Agrobacterium tumefaciens, Xylella fastidiosa), a list of common fungal genes with homologs in these plants, animals, and bacteria was produced (938 genes), as well as a list of exclusively fungal genes without homologs in these other genomes (60 genes). To ensure that the 60 genes were exclusively fungal, these were compared using TBLASTN to the major sequence databases at GenBank: NR (nonredundant), EST (expressed sequence tags), GSS (genome survey sequences), and HTGS (unfinished high-throughput genome sequences). This resulted in 17 yeast genes with homologs in other fungal genomes, but without known homologs in other organisms. These 17 core, fungal genes were not found to differ from other yeast genes in GC content or codon usage patterns. More intensive study is required of these 17 genes and other common fungal genes to discover unique features of fungi compared to other organisms.Reviewing Editor: Prof. David Gottman