Correlated rates of synonymous site evolution across plant genomes

Synonymous substitution rates have been shown to vary among evolutionary lineages of both nuclear and organellar genes across a broad range of taxonomic groups. In animals, rate heterogeneity does not appear to be correlated across nuclear and mitochondrial genes. In this paper, we contrast substitution rates in two plant groups and show that grasses evolve more rapidly than palms at synonymous sites in a mitochondrial, a nuclear, and a plastid gene. Furthermore, we show that the relative rates of synonymous substitution between grasses and palms are similar at the three loci. The correlation in synonymous substitution rates across genes is particularly striking because the three genes evolve at very different absolute rates. In contrast, relative rates of nonsynonymous substitution are not conserved among the three genes.      

Influence of secondary structure of polyglycine on some conformational properties

Values of four conformational properties, namely unperturbed dimension [r2]0, dipole moment [mu 2], mean squared optical anisotropy [gamma 2], and molar Kerr constant [mK], have been calculated for polyglycine chains allowing several combinations of the secondary structure with the aim of studying the dependence of these magnitudes on the secondary structure of the chain. Two different approaches to the secondary structure have been used. In the first, chains with all their units in a given conformation (random coil, alpha-helix or beta-sheet) are interrupted at several positions by one unit in a different conformation. In the second, chains with varying composition of two conformations alpha-helix/beta-sheet and beta-sheet/random coil were allowed and the results obtained compared with previous work for alpha-helix/random coil chains.     

Patterns of nucleotide substitution among simultaneously duplicated gene pairs in Arabidopsis thaliana

We characterized rates and patterns of synonymous and nonsynonymous substitution in 242 duplicated gene pairs on chromosomes 2 and 4 of Arabidopsis thaliana. Based on their collinear order along the two chromosomes, the gene pairs were likely duplicated contemporaneously, and therefore comparison of genetic distances among gene pairs provides insights into the distribution of nucleotide substitution rates among plant nuclear genes. Rates of synonymous substitution varied 13.8-fold among the duplicated gene pairs, but 90% of gene pairs differed by less than 2.6-fold. Average nonsynonymous rates were approximately fivefold lower than average synonymous rates; this rate difference is lower than that of previously studied nonplant lineages. The coefficient of variation of rates among genes was 0.65 for nonsynonymous rates and 0.44 for synonymous rates, indicating that synonymous and nonsynonymous rates vary among genes to roughly the same extent. The causes underlying rate variation were explored. Our analyses tentatively suggest an effect of physical location on synonymous substitution rates but no similar effect on nonsynonymous rates. Nonsynonymous substitution rates were negatively correlated with GC content at synonymous third codon positions, and synonymous substitution rates were negatively correlated with codon bias, as observed in other systems. Finally, the 242 gene pairs permitted investigation of the processes underlying divergence between paralogs. We found no evidence of positive selection, little evidence that paralogs evolve at different rates, and no evidence of differential codon usage or third position GC content.     

The predicted secondary structure of enolase.

The results of several secondary-structure prediction programs were combined to produce an estimate of the regions of alpha-helix, beta-sheet and reverse turn for both chicken skeletal-muscle and yeast enolase sequences. The predicted secondary-structure content of the chicken enzyme is 27% alpha-helix and less than 10% beta-sheet, whereas in the yeast enolase a similar helix content but virtually no sheet are predicted. These results are in fair agreement with published experimental estimates of the amount of secondary structure in the yeast enzyme. The enzyme appears to be formed from three domains.     

Estimating absolute rates of synonymous and nonsynonymous nucleotide substitution in order to characterize natural selection and date species divergences

The rate of molecular evolution can vary among lineages. Sources of this variation have differential effects on synonymous and nonsynonymous substitution rates. Changes in effective population size or patterns of natural selection will mainly alter nonsynonymous substitution rates. Changes in generation length or mutation rates are likely to have an impact on both synonymous and nonsynonymous substitution rates. By comparing changes in synonymous and nonsynonymous rates, the relative contributions of the driving forces of evolution can be better characterized. Here, we introduce a procedure for estimating the chronological rates of synonymous and nonsynonymous substitutions on the branches of an evolutionary tree. Because the widely used ratio of nonsynonymous and synonymous rates is not designed to detect simultaneous increases or simultaneous decreases in synonymous and nonsynonymous rates, the estimation of these rates rather than their ratio can improve characterization of the evolutionary process. With our Bayesian approach, we analyze cytochrome oxidase subunit I evolution in primates and infer that nonsynonymous rates have a greater tendency to change over time than do synonymous rates. Our analysis of these data also suggests that rates have been positively correlated.     

Secondary conformations and temperature effect on structural transformation of amyloid beta (1-28), (1-40) and (1-42) peptides

Secondary structure of three amyloid b-peptides [A beta(1-28), A beta(1-40) and A beta(1-42)] in the solid state was respectively determined by Fourier transform infrared (FT-IR) microspectroscopy. Their thermal-dependent structural transformation were also investigated by FT-IR microspectroscopy equipped with a thermal analyzer. The present result demonstrates that the solid-state A beta(1-28), A beta(1-40) and A beta(1-42) peptides showed a significant IR spectral difference in the amide I and II bands. The secondary conformation of A beta(1-28) peptide was the combination of major beta-sheet and minor alpha-helix with little random coil structures, but A beta(1-40) peptide showed the co-existence of major beta-sheet and minor random coil with little alpha-helix structures. A beta(1-42) peptide mainly consisted of the predominant b-sheet structure. Although the intact A beta(1-28), A beta(1-40) or A beta(1-42) peptide exhibits a different secondary structure, a similar beta-conformation may form after thermal treatment. A thermal-dependent transition was found for solid A beta(1-28) and A beta(1-40) peptides near 40 degrees C and 45 degrees C, respectively. There was no transition temperature for solid A beta(1-42) peptide, however, due to only a very little level of alpha-helix and random coil structure containing in the solid A beta(1-42) peptide. The thermal denaturation plays an important role in the structural transformation from alpha-helix/random coil to beta-sheet.     

Substitution rates in hepatitis delta virus

Substitution rates were estimated for the coding and noncoding regions of the hepatitis delta virus (HDV). The estimated rates of synonymous substitution in HDV were lower than the rates of substitution at nonsynonymous sites and in the noncoding region. HDV has lower synonymous substitution rates than the hepatitis C virus, though both are RNA viruses. The relatively low rate of synonymous substitution in HDV may be due to a strong preference of G and C nucleotides at third codon positions. Variation in substitution rate among HDV lineages may be correlated with the clinical development of the HDV-induced hepatitis. The phylogenetic tree inferred for 24 HDV strains reveals similarities between lineages isolated from the same geographic region. Correspondence to: W.-H. Li     

Uncorrected nucleotide bias in mtDNA can mimic the effects of positive Darwinian selection

The relative rates of nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions have been widely used to infer the action of natural selection from comparative sequence data. It is known, however, that mutational and repair biases can affect rates of evolution at both synonymous and nonsynonymous sites. More importantly, it is also known that synonymous sites are particularly prone to the effects of nucleotide bias. This means that nucleotide biases may affect the calculated ratio of substitution rates at synonymous and nonsynonymous sites. Using a large data set of animal mitochondrial sequences, we demonstrate that this is, in fact, the case. Highly biased nucleotide sequences are characterized by significantly elevated dN/dS ratios, but only when the nucleotide frequencies are not taken into account. When the analysis is repeated taking the nucleotide frequencies at each codon position into account, such elevated ratios disappear. These results suggest that the recently reported differences in dN/dS ratios between vertebrate and invertebrate mitochondrial sequences could be explained by variations in mitochondrial nucleotide frequencies rather than the effects of positive Darwinian selection.     

A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome

A model of DNA sequence evolution applicable to coding regions is presented. This represents the first evolutionary model that accounts for dependencies among nucleotides within a codon. The model uses the codon, as opposed to the nucleotide, as the unit of evolution, and is parameterized in terms of synonymous and nonsynonymous nucleotide substitution rates. One of the model's advantages over those used in methods for estimating synonymous and nonsynonymous substitution rates is that it completely corrects for multiple hits at a codon, rather than taking a parsimony approach and considering only pathways of minimum change between homologous codons. Likelihood-ratio versions of the relative-rate test are constructed and applied to data from the complete chloroplast DNA sequences of Oryza sativa, Nicotiana tabacum, and Marchantia polymorpha. Results of these tests confirm previous findings that substitution rates in the chloroplast genome are subject to both lineage-specific and locus-specific effects. Additionally, the new tests suggest tha the rate heterogeneity is due primarily to differences in nonsynonymous substitution rates. Simulations help confirm previous suggestions that silent sites are saturated, leaving no evidence of heterogeneity in synonymous substitution rates.      

Secondary structure of Escherichia coli glucosamine-6-phosphate deaminase from amino acid sequence and circular dichroism spectroscopy

The secondary structure of the purified glucosamine-6-phosphate deaminase from Escherichia coli K12 was investigated by both circular dichroism (CD) spectroscopy and empirical prediction methods. The enzyme was obtained by allosteric-site affinity chromatography from an overproducing strain bearing a pUC18 plasmid carrying the structural gene for the enzyme. From CD analysis, 34% of alpha-helix, 9% of parallel beta-sheet, 11% of antiparallel beta-sheet, 15% turns and 35% of non-repetitive structures, were estimated. A joint prediction scheme, combining six prediction methods with defined rules using several physicochemical indices, gave the following values: alpha-helix, 37%; beta-sheet, 22%; turns, 18% and coil, 23%. The structure predicted showed also a considerable degree of alternacy of alpha and beta structures; 64% of helices are amphipathic and 90% of beta-sheets are hydrophobic. Overall, the data suggest that deaminase has as dominant motif, an alpha/beta structure.     

Rates of nucleotide substitution and mammalian nuclear gene evolution. Approximate and maximum-likelihood methods lead to different conclusions

Rates and patterns of synonymous and nonsynonymous substitutions have important implications for the origin and maintenance of mammalian isochores and the effectiveness of selection at synonymous sites. Previous studies of mammalian nuclear genes largely employed approximate methods to estimate rates of nonsynonymous and synonymous substitutions. Because these methods did not account for major features of DNA sequence evolution such as transition/transversion rate bias and unequal codon usage, they might not have produced reliable results. To evaluate the impact of the estimation method, we analyzed a sample of 82 nuclear genes from the mammalian orders Artiodactyla, Primates, and Rodentia using both approximate and maximum-likelihood methods. Maximum-likelihood analysis indicated that synonymous substitution rates were positively correlated with GC content at the third codon positions, but independent of nonsynonymous substitution rates. Approximate methods, however, indicated that synonymous substitution rates were independent of GC content at the third codon positions, but were positively correlated with nonsynonymous rates. Failure to properly account for transition/transversion rate bias and unequal codon usage appears to have caused substantial biases in approximate estimates of substitution rates.     

pH- and thermal-dependent conformational transition of PGAIPG, a repeated hexapeptide sequence from tropoelastin

The secondary structure of PGAIPG (Pro-Gly-Ala-IIe-Pro-Gly), a repeated hexapeptide of tropoelastin, in buffer solution of different pH was determined by using attenuated total reflection-Fourier transform infrared (ATR-FTIR) spectroscopy. The thermal-dependent structural change of PGAIPG in aqueous solution or in solid state was also examined by thermal FTIR microspectroscopy. The conformation of PGAIPG in aqueous solution exhibited a pH-dependent structural characterization. A predominant peak at 1614 cm(-1) (aggregated beta-sheet) with a shoulder near 1560 cm(-1) (beta-sheet) appeared in pH 5.5-8.5 buffer solutions. A new broad shoulder at 1651 cm(-1) (random coil and/or alpha-helix) with 1614 cm(-1) was observed in the pH 4.5 buffer solution. However, the broad shoulder at 1651 cm(-1) was converted to a maximum peak at 1679 cm(-1) (beta-turn/antiparallel beta-sheet) when the pH shifted from 4.5 to 3.5, but the original pronounced peak at 1614 cm(-1) became a shoulder. Once the pH was lowered to 2.5, the IR spectrum of PGAIPG was dominated by major absorption at 1679 cm(-1) with a minor peak at 1552 cm(-1) (alpha-helix/random coil). The result indicates that the pH was a predominant factor to transform PGAIPG structure from aggregated beta-sheet (pH 8.5) to beta-turn/intermolecular antiparallel beta-sheet (pH 2.5). Moreover, a partial conformation of PGAIPG with minor alpha-helix/random coil structures was also explored in the lower pH buffer solution. There was no thermal-dependent structural change for solid-state PGAIPG. The thermal-induced formation of aggregated beta-sheet for PGAIPG in aqueous solution was found from 28 to 30 degrees C, however, which might be correlated with the formation of an opaque gel that turned from clear solution. The formation of aggregated beta-sheet structure for PGAIPG beyond 30 degrees C might be due to the intermolecular hydrogen bonded interaction between the hydrophobic PGAIPG fragments induced by coacervation.     

Genetic code redundancy and the evolutionary stability of protein secondary structure

The genetic code has an inherent bias towards some amino acids because of the variable number of synonymous codons per amino acid. The extent to which these biases are expressed in protein secondary structure is described through the analysis of the overall amino acid compositions of the alpha-helix, beta-sheet, beta-turn and random coil segments elucidated by X-ray crystallography. Given the concept of neutral mutation in proteins, the allocation of synonyms in the genetic code appears to protect secondary structures from amino acid changes and discourages the appearance of chemically complex residues. The level of protection is similar for each structural form, despite their clear preferences for certain amino acids. The organization of the code is therefore relevant to the preservation of conformation seen in the evolution of many protein families.     

Evolution of the cetacean mitochondrial D-loop region.

We sequenced the mitochondrial DNA D-loop regions from two cetacean species and compared these with the published D-loop sequences of several other mammalian species, including one other cetacean. Nucleotide substitution rates, DNA sequence simplicity, possible open reading frames (ORFs), and potential RNA secondary structure were investigated. The substitution rate is an order of magnitude lower than would be expected on the basis of reports on human sequence variation in this region but are consistent with interspecific primate and rodent D-loop sequence variation and with estimates of substitution rates from whole mitochondrial genomes. Deletions/insertions are less common in the cetacean D-loop than in other vertebrate species. Areas of high sequence simplicity (clusters of short repetitive motifs) across the region correspond to areas of high sequence divergence. Three regions predicted to form secondary structures are homologous to such putative structures in other species; however, the presumptive structures most conserved in cetaceans are different from those reported for other taxa. While all three species have possible long ORFs, only a short sequence of seven amino acids is shared with other mammalian species, and those changes that had occurred within it are all nonsynonymous. We conclude that DNA slippage, in addition to point mutation, contributes to the evolution of the D-loop and that regions of conserved secondary structure in cetaceans and an ORF are unlikely to contribute significantly to the conservation of the central region.     

Evolution of nucleic acids coding for ribonucleases: the mRNA sequence of mouse pancreatic ribonuclease

The cDNA of mouse pancreatic mRNA has been cloned. After the library was screened with a rat ribonuclease cDNA probe, the positive clones were isolated and sequenced. There were no differences from the previously determined protein sequence. The mRNA codes for a preribonuclease of 149 amino acid residues including a signal peptide of 25 amino acids. The 3' noncoding region has a length of 260 bp, and the total mRNA length is approximately 940 bp. Comparison with the rat pancreatic ribonuclease sequence showed a high rate of nucleotide substitution. Within the coding region, nonsynonymous and synonymous substitution rates are 4.3 X 10(-9) and 15 X 10(-9) nucleotide substitutions/site/year, respectively. The latter value is one of the highest rates observed in the molecular evolution of mammalian nuclear genes. In the signal sequences the synonymous substitution rate is much lower and about the same as the nonsynonymous rate. Signal sequences of other mouse and rat proteins also exhibit little difference between synonymous and nonsynonymous rates. The sequences of rat and mouse pancreatic ribonuclease messengers were compared with those of bovine pancreatic, seminal, and brain ribonuclease. While the 3' noncoding regions of rat and mouse are very similar, as are those of the three bovine messengers, there is no significant similarity between both rodent and the three bovine messengers for the greater part of these regions. There is a duplication of approximately 50 nucleotides in the 3' noncoding region of the bovine messengers, with a region rich in A and C in between. The presence of this structural feature may be correlated with recent gene duplications that have occurred in the bovine genome.     

Analysis of substitution events in HIV-1 vif gene of the Korean clade

Nucleotide and amino acid substitution pattern in vif gene of the Korean clade of HIV-1 isolated from Koreans were analyzed using consensus sequences. At nucleotide level, transition/transversion substitution ratio was 1.88, and nonsynonymous/synonymous substitution ratio was 2.67, suggesting a divergent evolution in the Korean clade. At amino acid level, there were 17 substitutions and G-->E substitution at position 37 may be responsible for change in predicted secondary structure.     

Synonymous and Nonsynonymous Substitutions in Mammalian Genes: Intragenic Correlations

Previous investigations indicated that synonymous and nonsynonymous substitution rates are correlated in mammalian genes. In the present work, this correlation has been studied at the intragenic level using a dataset of 48 orthologous genes from species belonging to at least four different mammalian orders. The results obtained show that the intragenic variability in synonymous rates is correlated with that of nonsynonymous rates. Moreover, the variation in GC level (and especially of C level) of silent positions along each gene is correlated with the variation in synonymous rate. These results reinforce the previous conclusions that synonymous and nonsynonymous rates as well as GC levels of silent positions are to some extent under common selective constraints. Received: 10 July 1997 / Accepted: 13 August 1997     

Not so different after all: a comparison of methods for detecting amino acid sites under selection

We consider three approaches for estimating the rates of nonsynonymous and synonymous changes at each site in a sequence alignment in order to identify sites under positive or negative selection: (1) a suite of fast likelihood-based "counting methods" that employ either a single most likely ancestral reconstruction, weighting across all possible ancestral reconstructions, or sampling from ancestral reconstructions; (2) a random effects likelihood (REL) approach, which models variation in nonsynonymous and synonymous rates across sites according to a predefined distribution, with the selection pressure at an individual site inferred using an empirical Bayes approach; and (3) a fixed effects likelihood (FEL) method that directly estimates nonsynonymous and synonymous substitution rates at each site. All three methods incorporate flexible models of nucleotide substitution bias and variation in both nonsynonymous and synonymous substitution rates across sites, facilitating the comparison between the methods. We demonstrate that the results obtained using these approaches show broad agreement in levels of Type I and Type II error and in estimates of substitution rates. Counting methods are well suited for large alignments, for which there is high power to detect positive and negative selection, but appear to underestimate the substitution rate. A REL approach, which is more computationally intensive than counting methods, has higher power than counting methods to detect selection in data sets of intermediate size but may suffer from higher rates of false positives for small data sets. A FEL approach appears to capture the pattern of rate variation better than counting methods or random effects models, does not suffer from as many false positives as random effects models for data sets comprising few sequences, and can be efficiently parallelized. Our results suggest that previously reported differences between results obtained by counting methods and random effects models arise due to a combination of the conservative nature of counting-based methods, the failure of current random effects models to allow for variation in synonymous substitution rates, and the naive application of random effects models to extremely sparse data sets. We demonstrate our methods on sequence data from the human immunodeficiency virus type 1 env and pol genes and simulated alignments.     

The Bacillus thuringiensis delta-endotoxin. Evidence for a two domain structure of the minimal toxic fragment

The conformational characteristics of the minimal toxic fragment of the delta-endotoxin from Bacillus thuringiensis berliner 1715 were examined by fluorescence and circular dichroism spectroscopy. This insecticidal protein, specifically toxic to lepidopteran species, was found to consist of two structural domains. Experimental evidence for this conclusion was provided by biphasic guanidine hydrochloride unfolding curves at different pH values and electrophoretic patterns of protease digests. Two stable fragments of comparable molecular weight were obtained using four different broad specificity proteolytic enzymes. A secondary structure model was constructed using seven B. thuringiensis toxin sequences. These toxins were selected on the basis of their limited sequence homology and represent all known insecticidal specificities. Despite this divergence, a consensus secondary structure pattern was obtained, confirming the structural homology among the toxins. The N-terminal halves of all toxins are predicted to be relatively rich in alpha-helix structure and the C-terminal parts to contain alternating beta-strand and coil structures. The latter seems characteristic for a beta-sheet conformation. Comparing this model to the unfolding data obtained by circular dichroism, whose far UV signal gives a measure of the alpha-helix content, allowed us to delineate the structural domains into the primary structure.     

Secondary structures of rat lipolytic enzymes: circular dichroism studies and relation to hydrophobic moments

To explore the secondary structures of lingual and pancreatic lipases, circular dichroism measurements were performed. Maximum average ellipticities were used to calculate the percentage of alpha-helices, beta-sheets, and random coils. Lingual lipase had an ellipticity of -20235 +/- 140 deg cm2/dmol (mean +/- SE) at 220 nm suggesting 60% alpha-helix, 20% beta-sheet and 20% random coil structure, but the mean ellipticity for pancreatic lipase was -14093 +/- 82 deg cm2/dmol (mean +/- SE) at 210 nm suggesting a 34.8% alpha-helical, 25% beta-sheet and 40% random coil secondary structure. An alpha-helical stretch of residues with a large hydrophobic moment ("globular" alpha-helix by hydrophobic moment plot) from amino acids 382 through 389 at the COOH-terminal end of lingual lipase was noted. This sequence, absent in pancreatic lipase, may account for the avid binding of lingual lipase to fat emulsion particles.