Testing for Selection on Synonymous Sites in Plant Mitochondrial DNA: The Role of Codon Bias and RNA Editing

Since plant mitochondrial genomes exhibit some of the slowest known synonymous substitution rates, it is generally believed that they experience exceptionally low mutation rates. However, the use of synonymous substitution rates to infer mutation rates depends on the implicit assumption that synonymous sites are evolving neutrally (or nearly so). To assess the validity of this assumption in plant mitochondrial genomes, we examined coding sequence for footprints of selection acting at synonymous sites. We found that synonymous sites exhibit an AT rich and pyrimidine skewed nucleotide composition compared to both non-synonymous sites and non-coding regions. We also found some evidence for selection associated with both biased codon usage and conservation of regulatory sequences involved in mRNA processing, although some of these findings are subject to alternative non-adaptive interpretations. Regardless, the inferred strength of selection appears too weak to account for the variation in substitution rates between the mitochondrial genomes of plants and other multicellular eukaryotes. Therefore, these results are consistent with the interpretation that plant mitochondrial genomes experience a substantially lower mutation rate rather than increased functional constraints acting on synonymous sites. Nevertheless, there are important nucleotide composition patterns (particularly the differences between synonymous sites and non-coding DNA) that remain largely unexplained.     

人群变异的分子基础:从单核苷酸多态性到单氨基酸多态性

单核苷酸多态性可以划分为位于基因编码区的SNP和非编码区的SNP两大种类;而在基因编码区的SNP还可以进一步划分为两个亚类:不改变氨基酸序列的同义SNP和改变氨基酸序列的非同义SNP.显然,非同义SNP将导致氨基酸序列的改变,即形成单氨基酸多态性.基于蛋白质组学方法,对亚洲人群血浆样本中的SAP进行了系统研究,发现某一特定SAP在纯合人群和杂合人群中可能与生理或病理性状有着不同的关联.更为重要的是,近期有研究发现,在生物体中广泛存在着RNA序列与DNA序列不一致的现象.导致这种差异的主要原因是在转录水平上存在着规模化的RNA编辑(被称为RNA编辑组,RNA editome).该发现表明,个体拥有的SAP中可能有一部分与基因组SNP无关,而是源于RNA编辑组.进一步推论,可能在翻译水平上存在着不依赖DNA和RNA序列的全新的SAP.     

An Integrated Sequence-Structure Database incorporating matching mRNA sequence, amino acid sequence and protein three-dimensional structure data.

We have constructed a non-homologous database, termed the Integrated Sequence-Structure Database (ISSD) which comprises the coding sequences of genes, amino acid sequences of the corresponding proteins, their secondary structure and straight phi,psi angles assignments, and polypeptide backbone coordinates. Each protein entry in the database holds the alignment of nucleotide sequence, amino acid sequence and the PDB three-dimensional structure data. The nucleotide and amino acid sequences for each entry are selected on the basis of exact matches of the source organism and cell environment. The current version 1.0 of ISSD is available on the WWW at http://www.protein.bio.msu.su/issd/ and includes 107 non-homologous mammalian proteins, of which 80 are human proteins. The database has been used by us for the analysis of synonymous codon usage patterns in mRNA sequences showing their correlation with the three-dimensional structure features in the encoded proteins. Possible ISSD applications include optimisation of protein expression, improvement of the protein structure prediction accuracy, and analysis of evolutionary aspects of the nucleotide sequence-protein structure relationship.     

Protein structure and the sequential structure of mRNA: α-Helix and β-sheet signals at the nucleotide level

A direct comparison of experimentally determined protein structures and their corresponding protein coding mRNA sequences has been performed. We examine whether real world data support the hypothesis that clusters of rare codons correlate with the location of structural units in the resulting protein. The degeneracy of the genetic code allows for a biased selection of codons which may control the translational rate of the ribosome, and may thus in vivo have a catalyzing effect on the folding of the polypeptide chain. A complete search for GenBank nucleotide sequences coding for structural entries in the Brookhaven Protein Data Bank produced 719 protein chains with matching mRNA sequence, amino acid sequence, and secondary structure assignment. By neural network analysis, we found strong signals in mRNA sequence regions surrounding helices and sheets. These signals do not originate from the clustering of rare codons, but from the similarity of codons coding for very abundant amino acid residues at the N- and C-termini of helices and sheets. No correlation between the positioning of rare codons and the location of structural units was found. The mRNA signals were also compared with conserved nucleotide features of 16S-like ribosomal RNA sequences and related to mechanisms for maintaining the correct reading frame by the ribosome. © 1996 Wiley-Liss, Inc.     

mRNA的序列、结构以及翻译速率与蛋白质结构的关系

mRNA所包含的核苷酸序列通过三联体密码子决定了蛋白质的氨基酸序列。但是, 由于对氨基酸同义密码使用频率上的差异, 密码子与反密码子相互作用效率上的不同, 以及密码子上下文关系和mRNA 不同区域二级结构上的差异, 造成了核糖体对mRNA 不同区域翻译速度上的差异, 加之共翻译折叠的作用, 使得mRNA 的序列和结构影响着蛋白质空间结构的形成。     

Genome-wide comparative analysis of simple sequence coding repeats among 25 insect species

We present a detailed genome-scale comparative analysis of simple sequence repeats within protein coding regions among 25 insect genomes. The repetitive sequences in the coding regions primarily represented single codon repeats and codon pair repeats. The CAG triplet is highly repetitive in the coding regions of insect genomes. It is frequently paired with the synonymous codon CAA to code for polyglutamine repeats. The codon pairs that are least repetitive code for polyalanine repeats. The frequency of hexanucleotide and dinucleotide motifs of codon pair repeats is significantly (p<0.001) different in the Drosophila species compared to the non-Drosophila species. However, the frequency of synonymous and non-synonymous codon pair repeats varies in a correlated manner (r(2)=0.79) among all the species. Results further show that perfect and imperfect repeats have significant association with the trinucleotide and hexanucleotide coding repeats in most of these insects. However, only select species show significant association between the numbers of perfect/imperfect hexamers and repeat coding for single amino acid/amino acid pair runs. Our data further suggests that genes containing simple sequence coding repeats may be under negative selection as they tend to be poorly conserved across species. The sequences of coding repeats of orthologous genes vary according to the known phylogeny among the species. In conclusion, the study shows that simple sequence coding repeats are important features of genome diversity among insects.     

Evolutionary pattern of protein architecture in mammal and fruit fly genomes

Mutations, which can alter amino acid constitution, contribute greatly to protein evolution. However, little is reported of their pattern during protein structural evolution. We investigated the distribution of non-synonymous single nucleotide polymorphisms (nsSNPs) and insertions/deletions (indels) along mammal and fruit fly proteins. We found the nsSNPs (and d(N)) and indels increased in protein boundary regions, and this pattern is inversely correlated with the distribution of protein domain density. Additionally, synonymous substitutions (and d(S)) are reduced in 5' and 3' regions, indicating more variable protein boundaries, compared with central interior. All evidence suggests that the inner part of coding sequences (CDSs) is comparatively conserved, whereas the 5' and 3' regions, with higher evolution rates, are more variable. We assumed that due to greater frequencies of nsSNPs and indels in adaptive regions of CDSs it could be easier to ultimately alter, gain, or lose amino acids, thus becoming the front line of protein evolution.     

Cloning of cDNA Sequences Encoding the Calcium-Binding Protein, Calmodulin, from Barley (Hordeum vulgare L.)

Full- and partial-length cDNAs encoding calmodulin mRNA have been cloned and sequenced from barley (Hordeum vulgare L.). Barley leaf mRNA, size-fractionated in sucrose density gradients, was used to synthesize double-stranded cDNA. The cDNA was cloned in λgt10 and screened with a synthetic, 14-nucleotide oligonucleotide probe, which was designed using the predicted coding sequences of the carboxy termini of spinach and wheat calmodulin proteins. The primary structure of barley calmodulin, predicted from DNA sequencing experiments, consists of 148 amino acids and differs from that of wheat calmodulin in only three positions. In two of the three positions, the amino acid changes are conservative, while the third change consists of an apparent deletion/insertion. The overall nucleotide sequence similarity between the amino acid coding regions of barley and vertebrate calmodulin mRNAs is approximately 77%. However, a region encoding 11 amino acids of the second Ca²⁺-binding domain is very highly conserved at the nucleotide level compared with the rest of the coding sequences (94% sequence identity between barley and chicken calmodulin mRNAs). Genomic Southern blots reveal that barley calmodulin is encoded by a single copy gene. This gene is expressed as a single size class of mRNA in all tissues of 7-day-old barley seedlings. In addition, these analyses indicate that a barley calmodulin cDNA coding region subclone is suitable as a probe for isolating calmodulin genes from other plants.     

Characterization of coding synonymous and non-synonymous variants in ADAMTS13 using ex vivo and in silico approaches

Synonymous variations, which are defined as codon substitutions that do not change the encoded amino acid, were previously thought to have no effect on the properties of the synthesized protein(s). However, mounting evidence shows that these "silent" variations can have a significant impact on protein expression and function and should no longer be considered "silent". Here, the effects of six synonymous and six non-synonymous variations, previously found in the gene of ADAMTS13, the von Willebrand Factor (VWF) cleaving hemostatic protease, have been investigated using a variety of approaches. The ADAMTS13 mRNA and protein expression levels, as well as the conformation and activity of the variants have been compared to that of wild-type ADAMTS13. Interestingly, not only the non-synonymous variants but also the synonymous variants have been found to change the protein expression levels, conformation and function. Bioinformatic analysis of ADAMTS13 mRNA structure, amino acid conservation and codon usage allowed us to establish correlations between mRNA stability, RSCU, and intracellular protein expression. This study demonstrates that variants and more specifically, synonymous variants can have a substantial and definite effect on ADAMTS13 function and that bioinformatic analysis may allow development of predictive tools to identify variants that will have significant effects on the encoded protein.     

Codon conservation in the influenza A virus genome defines RNA packaging signals

Genome segmentation facilitates reassortment and rapid evolution of influenza A virus. However, segmentation complicates particle assembly as virions must contain all eight vRNA species to be infectious. Specific packaging signals exist that extend into the coding regions of most if not all segments, but these RNA motifs are poorly defined. We measured codon variability in a large dataset of sequences to identify areas of low nucleotide sequence variation independent of amino acid conservation in each segment. Most clusters of codons showing very little synonymous variation were located at segment termini, consistent with previous experimental data mapping packaging signals. Certain internal regions of conservation, most notably in the PA gene, may however signify previously unidentified functions in the virus genome. To experimentally test the bioinformatics analysis, we introduced synonymous mutations into conserved codons within known packaging signals and measured incorporation of the mutant segment into virus particles. Surprisingly, in most cases, single nucleotide changes dramatically reduced segment packaging. Thus our analysis identifies cis-acting sequences in the influenza virus genome at the nucleotide level. Furthermore, we propose that strain-specific differences exist in certain packaging signals, most notably the haemagglutinin gene; this finding has major implications for the evolution of pandemic viruses.     

Comparison of dissimilarity patterns of E coli, yeast and mammalian tRNAs.

A number of experimental approaches have been developed for identification of recognition (identity) sites in tRNAs. Along with them a theoretical methodology has been proposed by McClain et al that is based on concomitant analysis of all tRNA sequences from a given species. This approach allows an evaluation of nucleotide combinations present in isoacceptor tRNAs specific for the given amino acid, and not present in equivalent positions in cloverleaf structure in other tRNAs of the same organism. These elements predicted from computer analysis of the databank could be tested experimentally for their participation in forming recognition sites. The correlation between theoretical predictions and experimental data appeared promising. The aim of the present work consisted of introducing further improvements into McClain's procedure by: i), introducing into analysis a variable region in tRNAs which had not been previously considered; to accomplish this, 'normalization' of variable nucleotides was suggested, based on primary and tertiary structures of tRNAs; ii), developing a new procedure for comparison of patterns for synonymous and non-synonymous tRNAs from different organisms; iii), analysis of 3- and 4-positional contacts between tRNAs and enzymes in addition to a formerly used 2-positional model. A systematic application of McClain's procedure to mammalian, yeast and E coli tRNAs led to the following results: i), imitancy patterns for non-synonymous tRNAs of any amino acid specificity and from any organisms analysed so far overlap by no more than 30%, providing a structural basis for discrimination with high fidelity between cognate and non-cognate tRNAs; ii), the predicted identity sites are non-randomly distributed within tRNA molecules; the dominant role is ascribed to only two regions--anticodon and amino acid stem which are located far apart from one another at extremes of all tRNA molecules; iii), the imitancy patterns for synonymous tRNAs in lower (yeast) and higher (mammalian) eukaryotes are similar but not identical; iv), distribution of predicted identity sites in the cloverleaf structure in prokaryotes and eukaryotes is essentially different: in eubacterial tRNAs the major role in recognition plays anticodon and/or amino acid acceptor stem, whereas in eukaryotic (both unicellular and multicellular) tRNAs the remaining part of the molecules is also involved in recognition; v), the imitancy patterns of synonymous tRNAs from prokaryotes and eukaryotes are dissimilar, this observation leads to the prediction that the tRNA identity sites for the same amino acid in prokaryotes and eukaryotes may differ.     

Selective pressures at a codon-level predict deleterious mutations in human disease genes

Deleterious mutations affecting biological function of proteins are constantly being rejected by purifying selection from the gene pool. The non-synonymous/synonymous substitution rate ratio (omega) is a measure of selective pressure on amino acid replacement mutations for protein-coding genes. Different methods have been developed in order to predict non-synonymous changes affecting gene function. However, none has considered the estimation of selective constraints acting on protein residues. Here, we have used codon-based maximum likelihood models in order to estimate the selective pressures on the individual amino acid residues of a well-known model protein: p53. We demonstrate that the number of residues under strong purifying selection in p53 is much higher than those that are strictly conserved during the evolution of the species. In agreement with theoretical expectations, residues that have been noted to be of structural relevance, or in direct association with DNA, were among those showing the highest signals of purifying selection. Conversely, those changing according to a neutral, or nearly neutral mode of evolution, were observed to be irrelevant for protein function. Finally, using more than 40 human disease genes, we demonstrate that residues evolving under strong selective pressures (omega<0.1) are significantly associated (p<0.01) with human disease. We hypothesize that non-synonymous change on amino acids showing omega<0.1 will most likely affect protein function. The application of this evolutionary prediction at a genomic scale will provide an a priori hypothesis of the phenotypic effect of non-synonymous coding single nucleotide polymorphisms (SNPs) in the human genome.     

Thioredoxin-catalyzed Refolding of Recombinant Protein: Refolding of Human Pro-urokinase

We have isolated four independent genomic DNA fragments encoding a rice storage protein, glutelin, by using its cDNA as a probe. Restriction mapping and sequencing analyses showed that all four clones encode type II mRNA of glutelin and have very similar nucleotide sequences. However, several differences were found among the four clones with respect to the nucleotide sequences of their coding and 5′-upstream regions, some of which confirmed corresponding changes in the amino acid sequences of the encoded proteins. These results indicate that the genes encoding type II mRNA species of glutelin constitute a multi-gene family of closely related members.     

Predicted functional RNAs within coding regions constrain evolutionary rates of yeast proteins

Functional RNAs (fRNAs) are being recognized as an important regulatory component in biological processes. Interestingly, recent computational studies suggest that the number and biological significance of functional RNAs within coding regions (coding fRNAs) may have been underestimated. We hypothesized that such coding fRNAs will impose additional constraint on sequence evolution because the DNA primary sequence has to simultaneously code for functional RNA secondary structures on the messenger RNA in addition to the amino acid codons for the protein sequence. To test this prediction, we first utilized computational methods to predict conserved fRNA secondary structures within multiple species alignments of Saccharomyces sensu strico genomes. We predict that as much as 5% of the genes in the yeast genome contain at least one functional RNA secondary structure within their protein-coding region. We then analyzed the impact of coding fRNAs on the evolutionary rate of protein-coding genes because a decrease in evolutionary rate implies constraint due to biological functionality. We found that our predicted coding fRNAs have a significant influence on evolutionary rates (especially at synonymous sites), independent of other functional measures. Thus, coding fRNA may play a role on sequence evolution. Given that coding regions of humans and flies contain many more predicted coding fRNAs than yeast, the impact of coding fRNAs on sequence evolution may be substantial in genomes of higher eukaryotes.