Patterns of DNA-Sequence Divergence Between <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. pseudoobscura</Emphasis>

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d _N) and codon usage bias (F _op), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.     

Variation in synonymous codon use and DNA polymorphism within the Drosophila genome

A strong negative correlation between the rate of amino-acid substitution and codon usage bias in Drosophila has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. To further explore this possibility we have investigated polymorphism and divergence at three kinds of sites: synonymous, nonsynonymous and intronic in relation to codon bias in D. melanogaster and D. simulans. We confirmed that protein evolution is one of the main explicative parameters for interlocus codon bias variation (r(2) approximately 40%). However, intron or synonymous diversities, which could have been expected to be good indicators of local interference [here defined as the additional increase of drift due to selection on tightly linked sites, also called 'genetic draft' by Gillespie (2000)] did not covary significantly with codon bias or with protein evolution. Concurrently, levels of polymorphism were reduced in regions of low recombination rates whereas codon bias was not. Finally, while nonsynonymous diversities were very well correlated between species, neither synonymous nor intron diversities observed in D. melanogaster were correlated with those observed in D. simulans. All together, our results suggest that the selective constraint on the protein is a stable component of gene evolution while local interference is not. The pattern of variation in genetic draft along the genome therefore seems to be instable through evolutionary times and should therefore be considered as a minor determinant of codon bias variance. We argue that selective constraints for optimal codon usage are likely to be correlated with selective constraints on the protein, both between codons within a gene, as previously suggested, and also between genes within a genome.     

Protein Evolutionary Rates Correlate with Expression Independently of Synonymous Substitutions in Helicobacter pylori

In free-living microorganisms, such as Escherichia coli and Saccharomyces cerevisiae, both synonymous and nonsynonymous substitution frequencies correlate with expression levels. Here, we have tested the hypothesis that the correlation between amino acid substitution rates and expression is a by-product of selection for codon bias and translational efficiency in highly expressed genes. To this end, we have examined the correlation between protein evolutionary rates and expression in the human gastric pathogen Helicobacter pylori, where the absence of selection on synonymous sites enables the two types of substitutions to be uncoupled. The results revealed a statistically significant negative correlation between expression levels and nonsynonymous substitutions in both H. pylori and E. coli. We also found that neighboring genes located on the same, but not on opposite strands, evolve at significantly more similar rates than random gene pairs, as expected by co-expression of genes located in the same operon. However, the two species differ in that synonymous substitutions show a strand-specific pattern in E. coli, whereas the weak similarity in synonymous substitutions for neighbors in H. pylori is independent of gene orientation. These results suggest a direct influence of expression levels on nonsynonymous substitution frequencies independent of codon bias and selective constraints on synonymous sites. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]     

The Evolution of Codon Preferences in Drosophila: A Maximum-Likelihood Approach to Parameter Estimation and Hypothesis Testing

Synonymous codon usage in related species may differ as a result of variation in mutation biases, differences in the overall strength and efficiency of selection, and shifts in codon preference—the selective hierarchy of codons within and between amino acids. We have developed a maximum-likelihood method to employ explicit population genetic models to analyze the evolution of parameters determining codon usage. The method is applied to twofold degenerate amino acids in 50 orthologous genes from D. melanogaster and D. virilis. We find that D. virilis has significantly reduced selection on codon usage for all amino acids, but the data are incompatible with a simple model in which there is a single difference in the long-term N _e, or overall strength of selection, between the two species, indicating shifts in codon preference. The strength of selection acting on codon usage in D. melanogaster is estimated to be |N _e s|≈ 0.4 for most CT-ending twofold degenerate amino acids, but 1.7 times greater for cysteine and 1.4 times greater for AG-ending codons. In D. virilis, the strength of selection acting on codon usage for most amino acids is only half that acting in D. melanogaster but is considerably greater than half for cysteine, perhaps indicating the dual selection pressures of translational efficiency and accuracy. Selection coefficients in orthologues are highly correlated (ρ= 0.46), but a number of genes deviate significantly from this relationship. Received: 20 December 1998 / Accepted: 17 February 1999     

Difference in larval type explains patterns of nonsynonymous substitutions in two ancient paralogs of the histone H3 gene in sea stars

SUMMARY Paralogous genes frequently show differences in patterns and rates of substitution that are typically attributed to different selection regimes, mutation rates, or local recombination rates. Here, two anciently diverged paralogous copies of the histone H3 gene in sea stars, the tandem‐repetitive early‐stage gene and a newly isolated gene with lower copy number that was termed the “putative late‐stage histone H3 gene” were analyzed in 69 species with varying mode of larval development. The two genes showed differences in relative copy number, overall substitution rates, nucleotide composition, and codon usage, but similar patterns of relative nonsynonymous substitution rates, when analyzed by the d_N/d_S ratio. Sea stars with a nonpelagic and nonfeeding larval type (i.e., brooding lineages) were observed to have d_N/d_S ratios that were larger than for nonbrooders but equal between the two paralogs. This finding suggested that demographic differences between brooding and nonbrooding lineages were responsible for the elevated d_N/d_S ratios observed for brooders and refuted a suggestion from a previous analysis of the early‐stage gene that the excess nonsynonymous substitutions were due to either (1) gene expression differences at the larval stage between brooders and nonbrooders or (2) the highly repetitive structure of the early‐stage histone H3 gene.     

A Survey of the Molecular Evolutionary Dynamics of Twenty-Five Multigene Families from Four Grass Taxa

We surveyed the molecular evolutionary characteristics of 25 plant gene families, with the goal of better understanding general processes in plant gene family evolution. The survey was based on 247 GenBank sequences representing four grass species (maize, rice, wheat, and barley). For each gene family, orthology and paralogy relationships were uncertain. Recognizing this uncertainty, we characterized the molecular evolution of each gene family in four ways. First, we calculated the ratio of nonsynonymous to synonymous substitutions (d _N/d _S) both on branches of gene phylogenies and across codons. Our results indicated that the d _N/d _S ratio was statistically heterogeneous across branches in 17 of 25 (68%) gene families. The vast majority of d _N/d _S estimates were <<1.0, suggestive of selective constraint on amino acid replacements, and no estimates were >1.0, either across phylogenetic lineages or across codons. Second, we tested separately for nonsynonymous and synonymous molecular clocks. Sixty-eight percent of gene families rejected a nonsynonymous molecular clock, and 52% of gene families rejected a synonymous molecular clock. Thus, most gene families in this study deviated from clock-like evolution at either synonymous or nonsynonymous sites. Third, we calculated the effective number of codons and the proportion of G+C synonymous sites for each sequence in each gene family. One or both quantities vary significantly within 18 of 25 gene families. Finally, we tested for gene conversion, and only six gene families provided evidence of gene conversion events. Altogether, evolution for these 25 gene families is marked by selective constraint that varies among gene family members, a lack of molecular clock at both synonymous and nonsynonymous sites, and substantial variation in codon usage. Received: 25 May 2000 / Accepted: 16 October 2000     

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]     

Likelihood Analysis of the Chalcone Synthase Genes Suggests the Role of Positive Selection in Morning Glories (<Emphasis Type="Italic">Ipomoea</Emphasis>)

Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoides, which are important for the pigmentation of flowers and act as attractants to pollinators. Genes encoding CHS constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. In morning glories (Ipomoea), five functional CHS genes (A–E) have been described. Phylogenetic analysis of the Ipomoea CHS gene family revealed that CHS A, B, and C experienced accelerated rates of amino acid substitution relative to CHS D and E. To examine whether the CHS genes of the morning glories underwent adaptive evolution, maximum-likelihood models of codon substitution were used to analyze the functional sequences in the Ipomoea CHS gene family. These models used the nonsynonymous/synonymous rate ratio ( = d_N/d_S) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio test suggested significant variation in selection pressure among amino acid sites, with a small proportion of them detected to be under positive selection along the branches ancestral to CHS A, B, and C. Positive Darwinian selection appears to have promoted the divergence of subfamily ABC and subfamily DE and is at least partially responsible for a rate increase following gene duplication.     

Patterns of synonymous codon usage in Drosophila melanogaster genes with sex-biased expression

The nonrandom use of synonymous codons (codon bias) is a well-established phenomenon in Drosophila. Recent reports suggest that levels of codon bias differ among genes that are differentially expressed between the sexes, with male-expressed genes showing less codon bias than female-expressed genes. To examine the relationship between sex-biased gene expression and level of codon bias on a genomic scale, we surveyed synonymous codon usage in 7276 D. melanogaster genes that were classified as male-, female-, or non-sex-biased in their expression in microarray experiments. We found that male-biased genes have significantly less codon bias than both female- and non-sex-biased genes. This pattern holds for both germline and somatically expressed genes. Furthermore, we find a significantly negative correlation between level of codon bias and degree of sex-biased expression for male-biased genes. In contrast, female-biased genes do not differ from non-sex-biased genes in their level of codon bias and show a significantly positive correlation between codon bias and degree of sex-biased expression. These observations cannot be explained by differences in chromosomal distribution, mutational processes, recombinational environment, gene length, or absolute expression level among genes of the different expression classes. We propose that the observed codon bias differences result from differences in selection at synonymous and/or linked nonsynonymous sites between genes with male- and female-biased expression.     

A Genomic Comparison of Faster-Sex, Faster-X, and Faster-Male Evolution Between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d_N) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d_N values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d_N. Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation. [Reviewing Editor: Dr. Willie J. Swanson]     

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.     

蕨类植物叶绿体基因ycf94的分子进化研究

ycf94基因是近年来在叶绿体基因组中新发现的一个基因,在蕨类植物中表现高度保守。该研究共选取94种蕨类植物,在系统发育背景下,对ycf94基因的结构特征、密码子偏好性、进化速率和适应性进化进行分析。结果表明, ycf94基因的密码子偏好性较弱,偏好使用以A/U结尾的密码子,且不同物种间的偏好性存在一定差异。密码子偏好性的形成主要受到突变压的影响,同时也存在其他因素的作用;基于凤尾蕨科和其他蕨类中ycf94基因的结构特征存在区别,对两者的分子替换速率进行了比较,表明颠换率、非同义替换率和ω值间存在显著差异;仅检测出1个正选择位点74A,强烈的负选择作用表明ycf94基因的结构和功能基本趋于稳定。这为蕨类系统发育分析提供了新依据,并提供了解析ycf94基因功能的线索。     

Synonymous Codon Usage,GC<Subscript>3</Subscript>, and Evolutionary Patterns Across Plastomes of Three Pooid Model Species: Emerging Grass Genome Models for Monocots

We have analyzed factors affecting the codon usage pattern of the chloroplasts genomes of representative species of pooid grass family. Correspondence analysis of relative synonymous codon usages (RSCU) showed that genes on secondary axis were correlated with their GC_3S values (all r > 0.3, p < 0.05), indicating mutational bias as an important selective force that shaped the variation in the codon usage among chloroplast genes. The Nc-plot showed that although a majority of the points with low-Nc values were lying below the expected curve, a few genes lied on the expected curve. Nc plot clearly showed that mutational bias plays a major role in codon biology across the monocot plastomes. The hydrophobicity and aromaticity of encoded proteins of each species were found to be other factors of codon usage variation. In the view of above light, besides natural selection, several other factors also likely to be involved in determining the selective constraints on codon bias in plastomes of pooid grass genomes. In addition, five codons (B. distachyon), seven codons (H. vulgare), and four codons (T. aestivum) were identified as optimal codons of the three grass chloroplasts. To identify genes evolving under positive selection, rates of nonsynonymous substitutions (Ka) and synonymous substitutions (Ks) were computed for all groups of orthologous gene pairs.     

Association Between Levels of Coding Sequence Divergence and Gene Misregulation in <Emphasis Type="Italic">Drosophila</Emphasis> Male Hybrids

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Carlo G. Artieri and Wilfried Haerty contributed equally to this publication.     

Gene expression and protein length influence codon usage and rates of sequence evolution in Populus tremula

Codon bias is generally thought to be determined by a balance between mutation, genetic drift, and natural selection on translational efficiency. However, natural selection on codon usage is considered to be a weak evolutionary force and selection on codon usage is expected to be strongest in species with large effective population sizes. In this paper, I study associations between codon usage, gene expression, and molecular evolution at synonymous and nonsynonymous sites in the long-lived, woody perennial plant Populus tremula (Salicaceae). Using expression data for 558 genes derived from expressed sequence tags (EST) libraries from 19 different tissues and developmental stages, I study how gene expression levels within single tissues as well as across tissues affect codon usage and rates sequence evolution at synonymous and nonsynonymous sites. I show that gene expression have direct effects on both codon usage and the level of selective constraint of proteins in P. tremula, although in different ways. Codon usage genes is primarily determined by how highly expressed a genes is, whereas rates of sequence evolution are primarily determined by how widely expressed genes are. In addition to the effects of gene expression, protein length appear to be an important factor influencing virtually all aspects of molecular evolution in P. tremula.     

Selection on Codon Usage for Error Minimization at the Protein Level

Given the structure of the genetic code, synonymous codons differ in their capacity to minimize the effects of errors due to mutation or mistranslation. I suggest that this may lead, in protein-coding genes, to a preference for codons that minimize the impact of errors at the protein level. I develop a theoretical measure of error minimization for each codon, based on amino acid similarity. This measure is used to calculate the degree of error minimization for 82 genes of Drosophila melanogaster and 432 rodent genes and to study its relationship with CG content, the degree of codon usage bias, and the rate of nucleotide substitution. I show that (i) Drosophila and rodent genes tend to prefer codons that minimize errors; (ii) this cannot be merely the effect of mutation bias; (iii) the degree of error minimization is correlated with the degree of codon usage bias; (iv) the amino acids that contribute more to codon usage bias are the ones for which synonymous codons differ more in the capacity to minimize errors; and (v) the degree of error minimization is correlated with the rate of nonsynonymous substitution. These results suggest that natural selection for error minimization at the protein level plays a role in the evolution of coding sequences in Drosophila and rodents.Reviewing Editor: Dr. Massimo Di Giulio     

Positive and Negative Selection in the β-Esterase Gene Cluster of the Drosophila melanogaster Subgroup

We examine the pattern of molecular evolution of the β-esterase gene cluster, including the Est-6 and ψEst-6 genes, in eight species of the Drosophila melanogaster subgroup. Using maximum likelihood estimates of nonsynonymous/synonymous rate ratios, we show that the majority of Est-6 sites evolves under strong (48% of sites) or moderate (50% of sites) negative selection and a minority of sites (1.5%) is under significant positive selection. Est-6 sites likely to be under positive selection are associated with increased intraspecific variability. One positively selected site is responsible for the EST-6 F/S allozyme polymorphism; the same site is responsible for the EST-6 functional divergence between species of the melanogaster subgroup. For ψEst-6 83.7% sites evolve under negative selection, 16% sites evolve neutrally, and 0.3% sites are under positive selection. The positively selected sites of ψEst-6 are located at the beginning and at the end of the gene, where there is reduced divergence between D. melanogaster and D. simulans; these regions of ψEst-6 could be involved in regulation or some other function. Branch-site-specific analysis shows that the evolution of the melanogaster subgroup underwent episodic positive selection. Collating the present data with previous results for the β-esterase genes, we propose that positive and negative selection are involved in a complex relationship that may be typical of the divergence of duplicate genes as one or both duplicates evolve a new function. [Reviewing Editor: Dr. Martin Kreitman]     

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.     

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