Evidence that purifying selection acts on promoter sequences

We tested whether functionally important sites in bacterial, yeast, and animal promoters are more conserved than their neighbors. We found that substitutions are predominantly seen in less important sites and that those that occurred tended to have less impact on gene expression than possible alternatives. These results suggest that purifying selection operates on promoter sequences.     

Relationship between G + C in silent sites of codons and amino acid composition of human proteins

Summary We have investigated the relationship between the G + C content of silent (synonymous) sites in codons and the amino acid composition of encoded proteins for approximately 1,600 human genes. There are positive correlations between silent site G + C and the proportions of codons for Arg, Pro, Ala, Trp, His, Gln, and Leu and negative ones for Tyr, Phe, Asn, Ile, Lys, Asp, Thr, and Glu. The median proteins coded by groups of genes that differ in silent-site G + C content also differ in amino acid composition, as do some proteins coded by homologous genes. The pattern of compositional change can be largely explained by directional mutation pressure, the genetic code, and differences in the frequencies of accepted amino acid substitutions; the shifts in protein composition are likely to be selectively neutral.Offprint requests to: D.W. Collins     

Codon usage changes and sequence dissimilarity between human and rat

Summary This paper reports on the relationship between the number of silent differences and the codon usage changes in the lineages leading to human and rat. Examination of 102 pairs of homologous genes gives rise to four main conclusions: (1) We have previously demonstrated the existence of a codon usage change (called the minor shift) between human and rat; this was confirmed here with a larger sample. For genes with extreme C+G frequencies, the C+G level in the third codon position is less extreme in rat than in human. (2) Protein similarity and percentage of positive differences are the two main factors that discriminate homologous genes when characterized by differences between rat and human. By definition, positive differences result from silent changes between A or T and C or G with a direction implying a C+G content variation in the same direction as the overall gene variation. (3) For genes showing both codon usage change and low protein similarity, a majority of amino acid replacements contributes to C+G level variation in positions I and II in the same direction as the variation in position III. This is thus a new example of protein evolution due to constraints acting at the DNA level. (4) In heavy isochores (high C+G content) no direct correlation exists between codon usage change (measured by the dissymmetry of differences) and silent dissimilarity. In light isochores the opposite situation is observed: modification of codon usage is associated with a high synonymous dissimilarity. This result shows that, in some cases, modification of constraints acting at the DNA level could accelerate divergence between genomes.     

Modulation of Base-Specific Mutation and Recombination Rates EnablesFunctional Adaptation Within the Context of the Genetic Code

The persistence of life requires populations to adapt at a rate commensurate with the dynamics of their environment. Successful populations that inhabit highly variable environments have evolved mechanisms to increase the likelihood of successful adaptation. We introduce a 64 × 64 matrix to quantify base-specific mutation potential, analyzing four different replicative systems, error-prone PCR, mouse antibodies, a nematode, and Drosophila. Mutational tendencies are correlated with the structural evolution of proteins. In systems under strong selective pressure, mutational biases are shown to favor the adaptive search of space, either by base mutation or by recombination. Such adaptability is discussed within the context of the genetic code at the levels of replication and codon usage.Supplementary material to this paper is available in electronic form.Reviewing Editor: Dr. Edward Trifonov     

Evidence that accumulation of mutants in a biofilm reflects natural selection rather than stress-induced adaptive mutation

The accumulation of mutant genotypes within a biofilm evokes the controversy over whether the biofilm environment induces adaptive mutation or whether the accumulation can be explained by natural selection. A comparison of the virulence of two strains of the dental pathogen Streptococcus mutans showed that rats infected with one of the strains accumulated a high proportion (average, 22%) of organisms that had undergone a deletion between two contiguous and highly homologous genes. To determine if the accumulation of deletion mutants was due to selection or to an increased mutation rate, accumulations of deletion mutants within in vitro planktonic and biofilm cultures and within rats inoculated with various proportions of deletion organisms were quantified. We report here that natural selection was the primary force behind the accumulation of the deletion mutants.     

Equal G and C contents in histone genes indicate selection pressures on mRNA secondary structure

Summary Protein-specific versus taxon-specific patterns of nucleotide frequencies were studied in histone genes. The third positions of codons have a (well-known) taxon-specific G+C level and a histone type-specific G/C ratio. This ratio counterbalances the G/C ratio in the first and second positions so that the overall G and C levels in the coding region become approximately equal. The compensation of the G/C ratio indicates a selection pressure at the mRNA level rather than a selection pressure or mutation bias at the DNA level or a selection pressure on codon usage. The structure of histone mRNAs is compatible with the hypothesis that the G/C compensation is due to selection pressures on mRNA secondary structure. Nevertheless, no specific motifs seem to have been selected, and the free energy of the secondary structures is only slightly lower than that expected on the basis of nucleotide frequencies.Offprint requests to: M. A. Huynen     

Codon usage in mitochondrial genomes: distinguishing context-dependent mutation from translational selection

We analyze the frequencies of synonymous codons in animal mitochondrial genomes, focusing particularly on mammals and fish. The frequencies of bases at 4-fold degenerate sites are found to be strongly influenced by context-dependent mutation, which causes correlations between pairs of neighboring bases. There is a pattern of excess of certain dinucleotides and deficit of others that is consistent across large numbers of species, despite the wide variation of single-nucleotide frequencies among species. In many bacteria, translational selection is an important influence on codon usage. In order to test whether translational selection also plays a role in mitochondria, we need to control for context-dependent mutation. Selection for translational accuracy can be detected by comparison of codon usage in conserved and variable sites in the same genes. We give a test of this type that works in the presence of context-dependent mutation. There is very little evidence for translational accuracy selection in the mitochondrial genes considered here. Selection for translational efficiency might lead to preference for codons that match the limited repertoire of anticodons on the mitochondrial tRNAs. This is difficult to detect because the effect would usually be in the same direction in comparable to codon families and so would not cause an observable difference in codon usage between families. Several lines of evidence suggest that this type of selection is weak in most cases. However, we found several cases where unusual bases occur at the wobble position of the tRNA, and in these cases, some evidence for selection on codon usage was found. We discuss the way that these unusual cases are associated with codon reassignments in the mitochondrial genetic code.     

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

ScopeSynonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.PurposeThis synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.Physiological and medical relevanceUnderstanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.     

Towards a resolution on the inherent methodological weakness of the "effective number of codons used by a gene"

Recently Anders Fuglsang provided a modified way for calculating N(c) when biased discrepancy is present in a gene [Biochem. Biophys. Res. Commun. 317 (2004) 957]. Instead of taking the average codon homozygosity for each synonymous family type (as proposed by Wright) [Gene 87 (1990) 23] Fuglsang considered codon homozygosity of each amino acid individually. Marsashi and Najafabadi [Biochem. Biophys. Res. Commun. 324 (2004) 1] in their recent article demonstrated that the readjustment for overestimation at the level of individual amino acids results in loss of considerable amount of information. Immediately after the publication of Marsashi and Najafabadi, Fuglsang proposed that codon homozygosities can be calculated based on the classical population genetics [Biochem. Biophys. Res. Commun. 327 (2005) 1]. Though Fuglsang's approach is a novel one, it fails when any of the amino acids are absent in a gene. However, the inherent cause of overestimation at the level of individual amino acids is still obscured in the literature. Here in this communication we have presented a general condition where effective number of codons is overestimated using Wright's formula and also we propose a new way to calculate N(c), which is independent of amino acid composition.     

基于相对密码子频率的芜菁花叶病毒基因组序列的进化分析

现有92株芜菁花叶病毒(TuMV)的全基因组序列已在GenBank报道,据分析报道其中58株不含重组序列。利用系统聚类法对92株TuMV的全基因组序列和58株TuMV全基因组序列的相对密码子频率RSCU值进行聚类分析。同时利用系统发育分析方法分析了这92株和58株TuMV全基因组序列。结果发现,92株芜菁花叶病毒株的密码子偏性聚类树与其系统进化树的一致度很低;而不含重组序列的58株芜菁花叶病毒株的密码子偏性聚类树与其系统进化树的一致度却非常高,且与寄生宿主类型基本对应。这表明在不存在重组的情况下,TuMV密码子频率的偏性可能是宿主内的一种选择压力,影响TuMV基因组的点突变进化方向,促使TuMV适应宿主内环境。     

Common silent mutations in all types of hereditary complement C1q deficiencies

Hereditary complete deficiency of complement component C1q is a rare genetic disorder that is associated with severe recurrent infections and a high prevalence of lupus-erythematosus-like symptoms. In the past, several single nucleotide polymorphisms have been identified in all three genes coding for the C1q A, B, and C chains. These point mutations which either lead to termination codons, frameshift, or amino acid exchanges were thought to be responsible for these defects as no other nonsense or missense mutations were found. As a result of the aberrations, either a nonfunctional C1q antigen is present or no C1q protein is detectable in the patients' sera. Screening 46 individuals from seven families with different forms of C1q deficiencies identified a homologous silent mutation at position Gly70 (GGG>GGA) of the C1q A gene of all 11 C1q-deficient patients. A high number of family members that were heterozygous for the coding mutations carried the silent mutation in the homozygous (18%) or heterozygous (36%) state. In addition to the Gly70 mutation in the A gene, another homozygous silent mutation (C gene at position Pro14, CCT>CCC) was detected in all C1q-deficient patients.     

Effects of codon modification on human BMP2 gene expression in tobacco plants

Bone morphogenetic protein 2 (BMP2) has great potential in therapeutic applications. We are working on generating transgenic plants as a bioreactor to produce BMP2. We have studied the effects of codon optimization on the expression of human BMP2 (hBMP2) in tobacco plants. Three modified hBMP2 genes were transformed into tobacco under the control of either cauliflower mosaic virus 35S (CaMV35S) promoter or double-CaMV35S promoter plus alfalfa mosaic virus (AMV) enhancer. The fused β-glucuronidase (GUS) reporter gene was used to facilitate the assay of protein expression. The results indicated that codon optimization could increase the protein expression level obviously under CaMV35S promoter. However, under relatively stronger initiation condition (double-CaMV35S promoter plus AMV enhancer), only the gene with the lowest degree of codon optimization could increase the protein expression level. Our findings suggest that the action of codon optimization may be influenced by the factors of promoter strength and A+T content in tobacco plants.     

Evidence for purifying selection acting on silent sites in BRCA1

In mammals, it is usually assumed that selection cannot be strong enough to act on nucleotide mutations that do not cause a change at the protein level (i.e. 'silent' or 'synonymous' mutations). Here we report the results of a molecular evolutionary analysis of BRCA1. We find a repeatable pronounced peak in the ratio of nonsynonymous to synonymous substitutions between codons 200-300. Unusually, this peak is caused by a plummet in the silent-site rate of evolution. The most parsimonious interpretation of these data is that purifying selection is acting on silent sites.     

Nonrandomness of point mutation as reflected in nucleotide substitutions in pseudogenes and its evolutionary implications

Summary We have obtained a revised estimate of the pattern of point mutation by considering more pseudogene sequences. Compared with our previous estimate, it agrees better with expectations based on the double-strand structure of DNA. The revised pattern, like the previous one, indicates that mutation occurs nonrandomly among the four nucleotides. In particular, the proportion of transitional mutations (59%) is almost twice as high as the value (33%) expected under random mutation. The same high proportion of transitions is observed in synonymous substitutions in genes. The proportion of transitional changes observed among electrophoretic variants of human hemoglobin is about the same as that predicted by the revised pattern of mutation. We also show that nonrandom mutation increases, by about 15%, the proportion of synonymous mutations due to single-nucleotide changes in the codon table, and increases, from 10% to 50%, the rate of synonymous mutation in the seven genes studied. However, nonrandom mutation reduces (by about 10%) the proportion of polar changes among nonsynonymous mutations in a gene. As far as single-nucleotide changes (in the codon table) are concerned, nonrandom mutation only slightly favors relatively conservative amino acid interchanges, and has virtually no effect on the proportions of radical changes and nonsense mutations.     

Reducing mutation load through sexual selection on males

Mutation load is a key parameter in evolutionary theories, but relatively little empirical information exists on the mutation load of populations, or the elimination of this load through selection. We manipulated the opportunity for sexual selection within a mutation accumulation divergence experiment to determine how sexual selection on males affected the accumulation of mutations contributing to sexual and nonsexual fitness. Sexual selection prevented the accumulation of mutations affecting male mating success, the target trait, as well as reducing mutation load on productivity, a nonsexual fitness component. Mutational correlations between mating success and productivity (estimated in the absence of sexual selection) were positive. Sexual selection significantly reduced these fitness component correlations. Male mating success significantly diverged between sexual selection treatments, consistent with the fixation of genetic differences. However, the rank of the treatments was not consistent across assays, indicating that the mutational effects on mating success were conditional on biotic and abiotic context. Our experiment suggests that greater insight into the genetic targets of natural and sexual selection can be gained by focusing on mutational rather than standing genetic variation, and on the behavior of trait variances rather than means.     

Effect of codon adaptation on codon-level and gene-level translation efficiency in vivo

Background

There is a significant difference between synonymous codon usage in many organisms, and it is known that codons used more frequently generally showed efficient decoding rate. At the gene level, however, there are conflicting reports on the existence of a correlation between codon adaptation and translation efficiency, even in the same organism.

Results

To resolve this issue, we cultured Escherichia coli under conditions designed to maintain constant levels of mRNA and protein and subjected the cells to ribosome profiling (RP) and mRNA-seq analyses. We showed that the RP results correlated more closely with protein levels generated under similar culture conditions than with the mRNA abundance from the mRNA-seq. Our result indicated that RP/mRNA ratio could be used as a measure of translation efficiency at gene level. On the other hand, the RP data showed that codon-specific ribosome density at the decoding site negatively correlated with codon usage, consistent with the hypothesis that preferred codons display lower ribosome densities due to their faster decoding rate. However, highly codon-adapted genes showed higher ribosome densities at the gene level, indicating that the efficiency of translation initiation, rather than higher elongation efficiency of preferred codons, exerted a greater effect on ribosome density and thus translation efficiency.

Conclusions

These findings indicate that evolutionary pressure on highly expressed genes influenced both codon bias and translation initiation efficiency and therefore explains contradictory findings that codon usage bias correlates with translation efficiency of native genes, but not with the artificially created gene pool, which was not subjected to evolution pressure.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1115) contains supplementary material, which is available to authorized users.     

Role in tumorigenesis of silent mutations in the TP53 gene

Over 10,000 mutations in the TP53 suppressor gene have been recorded in the International Agency for Research on Cancer (IARC) tumor data base. About 4% of these mutations are silent. It is a question whether these mutations play a role in tumor development. In order to approach this question, we asked whether the reported silent mutations are randomly distributed throughout the TP53 gene. The p53 data base was searched exon by exon. From the frequency of codons with no silent mutations, the average number of silent mutations per codon for each exon was calculated using the Poisson distribution. The results indicate the distribution to be non-random. About one-third of all silent mutations occur in “hot-spots” and after subtraction of these hot-spots, the remaining silent mutations are randomly distributed. In addition, the percentage of silent mutations among the total in the silent mutation hot-spots is close to that expected for random mutation. We conclude that most of the silent mutations recorded in tumors play no role in tumor development and that the percentage of silent mutation is an indication of the amount of random mutation during tumorigenesis. Silent mutations occur to a significantly different extent in different tumor types. Tumors of the esophagus and colon have a low frequency of silent mutations, tumors of the prostate have a high frequency.     

流感病毒基因的密码子偏好性及聚类分析

流行性感冒病毒是一种造成人类及动物患流行性感冒的RNA病毒,它造成急性上呼吸道感染,并由空气迅速传播,在世界各地常有周期性的大流行。根据该病毒的基因组CDS序列,探讨了基因组序列密码子的使用模式和特性,并进行了病毒间的聚类分析。结果表明:流感病毒的G+C含量均低于A+U含量,偏向使用以A、U结尾的密码子的程度比使用以G、C结尾的较高,CUG、UCA、AGU、AGC、AGA、AGG、GUG、CCA、ACA、GGA、GCA、AUU、UGA、CAU、CAA、AAU、AAA、GAA等18个密码子为流感病毒共有的偏好性密码子,且以A结尾的居多,尤其偏爱AGA、GGA。聚类结果表明首先亚洲流感病毒H2N2和香港流感病毒H2N2聚为一类,亚洲流感病毒H1N1和俄罗斯流感病毒H1N1聚为一类,1997年和2003年~2004年发生的人禽流感聚为一类,说明它们的密码子使用的偏好性相似;而2009年爆发的甲型H1N1流感和任何一个流感的距离都比较远,说明甲型H1N1流感病毒是一种新型的病毒,不同于以往任何一种流感病毒。     

A combination dimensionality reduction approach to codon position patterns of eubacteria based on their complete genomes

Graphical techniques have become powerful tools for the visualization and analysis of complicated biological systems. However, we cannot give such a graphical representation in a 2D/3D space when the dimensions of the represented data are more than three dimensions. The proposed method, a combination dimensionality reduction approach (CDR), consists of two parts: (i) principal component analysis (PCA) with a newly defined parameter ρ and (ii) locally linear embedding (LLE) with a proposed graphical selection for its optional parameter k. The CDR approach with ρ and k not only avoids loss of principal information, but also sufficiently well preserves the global high-dimensional structures in low-dimensional space such as 2D or 3D. The applications of the CDR on characteristic analysis at different codon positions in genome show that the method is a useful tool by which biologists could find useful biological knowledge.