The Role of Context-Dependent Mutations in Generating Compositional and Codon Usage Bias in Grass Chloroplast DNA

Abstract The influence of local base composition on mutations in chloroplast DNA (cpDNA) is studied in detail and the resulting, empirically derived, mutation dynamics are used to analyze both base composition and codon usage bias. A 4 × 4 substitution matrix is generated for each of the 16 possible flanking base combinations (contexts) using 17,253 noncoding sites, 1309 of which are variable, from an alignment of three complete grass chloroplast genome sequences. It is shown that substitution bias at these sites is correlated with flanking base composition and that the A+T content of these flanking sites as well as the number of flanking pyrimidines on the same strand appears to have general influences on substitution properties. The context-dependent equilibrium base frequencies predicted from these matrices are then applied to two analyses. The first examines whether or not context dependency of mutations is sufficient to generate average compositional differences between noncoding cpDNA and silent sites of coding sequences. It is found that these two classes of sites exist, on average, in very different contexts and that the observed mutation dynamics are expected to generate significant differences in overall composition bias that are similar to the differences observed in cpDNA. Context dependency, however, cannot account for all of the observed differences: although silent sites in coding regions appear to be at the equilibrium predicted, noncoding cpDNA has a significantly lower A+T content than expected from its own substitution dynamics, possibly due to the influence of indels. The second study examines the codon usage of low-expression chloroplast genes. When context is accounted for, codon usage is very similar to what is predicted by the substitution dynamics of noncoding cpDNA. However, certain codon groups show significant deviation when followed by a purine in a manner suggesting some form of weak selection other than translation efficiency. Overall, the findings indicate that a full understanding of mutational dynamics is critical to understanding the role selection plays in generating composition bias and sequence structure.     

鲨烯合酶基因的密码子偏性分析

为了分析鲨烯合酶(squalene synthase, SS)基因密码子的使用方式及其影响因素,利用codon W和SPSS 16.0软件对47条来自不同物种的SS基因进行多元统计分析、对应性分析.SS基因密码子1~3位碱基的GC含量(GC1, GC2和GC3)依次为51.33%、34.65%和54.37%,3个位点的GC含量均呈极显著相关关系(p<0.01),对应性分析的结果表明,第1轴显示30.71%的差异,有效密码子数和GC3、GC1和GC2的均值与GC3之间的相关性均达极显著水平(p<0.01).筛选出的26个最优密码子的第3位碱基均为G或C.以MEGA 5.0构建的基于SS蛋白质序列的进化树比基于RSCU的聚类更符合传统的系统发育观点.SS基因密码子偏好以G/C结尾,使用模式受选择和突变影响,突变对密码子偏好影响较大.     

Codon Usage in Pathogenic Entamoeba histolytica

Summary An analysis of 4680 codons expressed by pathogenic Entamoeba histolytica showed the A+U content of coding sequences to be 67%. The preference for A+U resulted in an unusual codon usage with an A+U content of 84% in the third codon position. The data show a remarkable similarity to those obtained for Plasmodium falciparum.     

Codon usage and base composition inRickettsia prowazekii

Codon usage and base composition in sequences from the A + T-rich genome ofRickettsia prowazekii, a member of the alpha Proteobacteria, have been investigated. Synonymous codon usage patterns are roughly similar among genes, even though the data set includes genes expected to be expressed at very different levels, indicating that translational selection has been ineffective in this species. However, multivariate statistical analysis differentiates genes according to their G + C contents at the first two codon positions. To study this variation, we have compared the amino acid composition patterns of 21R. prowazekii proteins with that of a homologous set of proteins fromEscherichia coli. The analysis shows that individual genes have been affected by biased mutation rates to very different extents: genes encoding proteins highly conserved among other species being the least affected. Overall, protein coding and intergenic spacer regions have G + C content values of 32.5% and 21.4%, respectively. Extrapolation from these values suggests thatR. prowazekii has around 800 genes and that 60–70% of the genome may be coding. Correspondence to: S.G.E. Andersson     

Codon bias of the gene for chloroplast glycerol-3-phosphate acyltransferase in Camellia sinensis (L.) O. Kuntze

The analysis on codon usage bias of GPAT gene of Camellia sinensis (L.) O. Kuntze may provide a basis for understanding the evolution relationship of C. sinensis and for selecting appropriate host expression systems to improve the expression of target genes. In the present study, the coding sequence of CsGPAT was analyzed with CodonW, CHIPS and CUSP programs, and compared with the genome of C. sinensis and GPAT genes of other 9 plant species. Our results showed that the cluster tree based on CDs could reveal the evolutional relations among the 10 plant species, whereas the cluster tree based on relative synonymous codon usage (RSCU) could not. There were 31 codons showing distinct usage differences between CsGPAT and genome of Escherichia coli, 21 between CsGPAT and yeast, but 13 between CsGPAT and Arabidopsis thaliana. But there were slightly fewer differences in codon usage between CsGPAT and A. thaliana. Therefore, the A. thaliana expression system may be more suitable for the expression of CsGPAT. These results may improve our understanding of the codon usage bias and functional studies of CsGPAT.     

基因及其表达调控中遗传密码选择的偏爱性

遗传密码在基因及其表达调控中具有明显的选择性.在低等生物及细胞器基因组中,同义密码优先选择A、T;在高等生物的核基因组中,同义密码首先考虑C、G;编码基因的邻近序列对基因转录调控影响很大.环境因素与遗传密码的选择有关.     

Selection Conflicts,Gene Expression,and Codon Usage Trends in Yeast

Synonymous codon usage in yeast appears to be influenced by natural selection on gene expression, as well as regional variation in compositional bias. Because of the large number of potential targets of selection (i.e., most of the codons in the genome) and presumed small selection coefficients, codon usage is an excellent model for studying factors that limit the effectiveness of selection. We use factor analysis to identify major trends in codon usage for 5836 genes in Saccharomyces cerevisiae. The primary factor is strongly correlated with gene expression, consistent with the model that a subset of codons allows for more efficient translation. The secondary factor is very strongly correlated with third codon position GC content and probably reflects regional variation in compositional bias. We find that preferred codon usage decreases in the face of three potential limitations on the effectiveness of selection: reduced recombination rate, increased gene length, and reduced intergenic spacing. All three patterns are consistent with the Hill–Robertson effect (reduced effectiveness of selection among linked targets). A reduction in gene expression in closely spaced genes may also reflect selection conflicts due to antagonistic pleiotropy.     

Analysis of Synonymous Codon Usage Bias in 09H1N1

A novel subtype of influenza A virus 09H1N1 has rapidly spread across the world. Evolutionary analyses of this virus have revealed that 09H1N1 is a triple reassortant of segments from swine, avian and human influenza viruses. In this study, we investigated factors shaping the codon usage bias of 09H1N1 and carried out cluster analysis of 60 strains of influenza A virus from different subtypes based on their codon usage bias. We discovered that more preferentially used codons of 09H1N1 are A-ended or U-ended...     

Evolutionary Basis of Codon Usage and Nucleotide Composition Bias in Vertebrate DNA Viruses

Understanding the extent and causes of biases in codon usage and nucleotide composition is essential to the study of viral evolution, particularly the interplay between viruses and host cells or immune responses. To understand the common features and differences among viruses we analyzed the genomic characteristics of a representative collection of all sequenced vertebrate-infecting DNA viruses. This revealed that patterns of codon usage bias are strongly correlated with overall genomic GC content, suggesting that genome-wide mutational pressure, rather than natural selection for specific coding triplets, is the main determinant of codon usage. Further, we observed a striking difference in CpG content between DNA viruses with large and small genomes. While the majority of large genome viruses show the expected frequency of CpG, most small genome viruses had CpG contents far below expected values. The exceptions to this generalization, the large gammaherpesviruses and iridoviruses and the small dependoviruses, have sufficiently different life-cycle characteristics that they may help reveal some of the factors shaping the evolution of CpG usage in viruses. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]     

Reduced Selection for Codon Usage Bias in <Emphasis Type="Italic">Drosophila miranda</Emphasis>

Biased codon usage in many species results from a balance among mutation, weak selection, and genetic drift. Here I show that selection to maintain biased codon usage is reduced in Drosophila miranda relative to its ancestor. Analyses of mutation patterns in noncoding DNA suggest that the extent of this reduction cannot be explained by changes in mutation bias or by biased gene conversion. Low levels of variability in D. miranda relative to its sibling species, D. pseudoobscura, suggest that it has a much smaller effective population size. Reduced codon usage bias in D. miranda may thus result from the reduced efficacy of selection against newly arising mutations to unpreferred codons. [Reviewing Editor: Dr. Richard Kliman]     

Codon Usage Decreases the Error Minimization Within the Genetic Code

The genetic code is not random but instead is organized in such a way that single nucleotide substitutions are more likely to result in changes between similar amino acids. This fidelity, or error minimization, has been proposed to be an adaptation within the genetic code. Many models have been proposed to measure this adaptation within the genetic code. However, we find that none of these consider codon usage differences between species. Furthermore, use of different indices of amino acid physicochemical characteristics leads to different estimations of this adaptation within the code. In this study, we try to establish a more accurate model to address this problem. In our model, a weighting scheme is established for mistranslation biases of the three different codon positions, transition/transversion biases, and codon usage. Different indices of amino acids physicochemical characteristics are also considered. In contrast to pervious work, our results show that the natural genetic code is not fully optimized for error minimization. The genetic code, therefore, is not the most optimized one for error minimization, but one that balances between flexibility and fidelity for different species.     

Codon usage patterns among genes for lepidopteran hemolymph proteins

Patterns in codon usage were examined for the coding regions of the 23 known lepidopteran hemolymph proteins. Coding triplets are GC rich at the third position and a significant linear relationship between GC content of silent and nonsilent (replacement) sites was demonstrated. Intron GC content was significantly lower than in coding regions and no relationship between intron GC content and the same at silent and nonsilent sites was found. Though hemolymph proteins are all produced by the same tissue—fat body—significantly less bias was observed when all moth sequences were pooled than when sequences of the two major species were analyzed separately, as predicted by the genome hypothesis. In cases where no statistically significant bias was observed, polar or acidic basic amino acids were almost exclusively involved. Calculation of codon adaptation indices (CAI) was of limited value in quantifying the degree of codon bias and probably reflects the complexity of multicellular-organism life cycles and the changing patterns of gene expression over different developmental stages. Correspondence to: D.R. Frohlich     

Synonymous Codon Usage, Accuracy of Translation, and Gene Length in Caenorhabditis elegans

In many unicellular organisms, invertebrates, and plants, synonymous codon usage biases result from a coadaptation between codon usage and tRNAs abundance to optimize the efficiency of protein synthesis. However, it remains unclear whether natural selection acts at the level of the speed or the accuracy of mRNAs translation. Here we show that codon usage can improve the fidelity of protein synthesis in multicellular species. As predicted by the model of selection for translational accuracy, we find that the frequency of codons optimal for translation is significantly higher at codons encoding for conserved amino acids than at codons encoding for nonconserved amino acids in 548 genes compared between Caenorhabditis elegans and Homo sapiens. Although this model predicts that codon bias correlates positively with gene length, a negative correlation between codon bias and gene length has been observed in eukaryotes. This suggests that selection for fidelity of protein synthesis is not the main factor responsible for codon biases. The relationship between codon bias and gene length remains unexplained. Exploring the differences in gene expression process in eukaryotes and prokaryotes should provide new insights to understand this key question of codon usage. Received: 18 June 2000 / Accepted: 10 November 2000     

A New Test for Selection Applied to Codon Usage in Drosophila simulans and D. mauritiana

In many organisms, synonymous codon usage is biased by a history of natural selection. However, codon bias, itself, does not indicate that selection is ongoing; it may be a vestige of past selection. Simple statistical tests have been devised to infer ongoing selection on codon usage by comparing the derived state frequency spectra at polymorphic sites segregating either derived preferred codons or derived unpreferred codons; if selection is effective, the frequency of derived states should be higher in the former. We propose a new test that uses the inferred degree of preference, essentially calculating the correlation of derived state frequency and the difference in preference between the derived and the ancestral states; the correlation should be positive if selection is effective. When implementing the test, derived and ancestral states can be assigned by parsimony or on the basis of relative probability. In either case, statistical significance is estimated by a simple permutation test. We explored the statistical power of the test by sampling polymorphism data from 14 loci in 16 strains of D. simulans, finding that the test retains 80% power even when quite a few of the data are discarded. The power of the test likely reflects better use of multiple features of the data, combining population frequencies of polymorphic variants and quantitative estimates of codon preferences. We also applied this novel test to 14 newly sequenced loci in five strains of D. mauritiana, showing for the first time ongoing selection on codon usage in this species.     

Phylogenetic Diversity of New Zealand Gelidiales as Revealed by rbcL Sequence Data

Diversity and phylogenetic relationships of New Zealand representatives of the red algal order Gelidiales have been examined using rbcL sequence data. Extensive field collections have been made from throughout the New Zealand region. Six genera have been reported previously from New Zealand (Capreolia, Gelidium, Pterocladia, Pterocladiella, Pterocladiastrum, Ptilophora). This research has revealed species with very restricted local distributions, as well as the discovery of several undescribed, cryptic taxa. The common and widespread Gelidium caulacantheum is confirmed to be more closely related to Capreolia than to other species of Gelidium. The generic concept of Capreolia, based on life history characters, will need to be modified to accommodate additional species possessing “Gelidium” life histories. A species endemic to New Zealand, Gelidium ceramoides, has been found to differ significantly from all other members of the Gelidiales and requires reclassification in another genus and order. Examination of field collections and herbarium specimens in addition to molecular sequence data have led us to conclude that specimens previously placed in the genera Ptilophora and Pterocladiastrum belong within Pterocladia lucida.     

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     

植物KCO基因密码子使用特性

以植物钾离子外排通道（K’channeloutward．rectifier,KCO）基因为研究对象,运用CodonW软件分析了75个植物KCO基因密码子的使用模式,探讨密码子的使用模式和影响密码子使用的各种可能因素。结果表明：碱基组成差异（r=0．961,P〈0．01）和自然选择（r=0．568,P〈0．01）是影响密码子使用的主要因素,并且高表达的基因强烈偏爱使用以G或C结尾的密码子。确定了UUC、CUC等26个均以G／C结尾的密码子为植物KcD基因的高表达优越密码子。