Evolutionary analyses of DNA sequences subject to constraints of secondary structure.

Evolutionary models appropriate for analyzing nucleotide sequences that are subject to constraints on secondary structure are developed. The models consider the evolution of pairs of nucleotides, and they incorporate the effects of base-pairing constraints on nucleotide substitution rates by introducing a new parameter to extensions of standard models of sequence evolution. To illustrate some potential uses of the models, a likelihood-ratio test is constructed for the null hypothesis that two (prespecified) regions of DNA evolve independently of each other. The sampling properties of the test are explored via simulation. The test is then incorporated into a heuristic method for identifying the location of unknown stems. The test and related procedures are applied to data from ribonuclease P RNA sequences of bacteria.

Testing for equality of evolutionary rates.

A likelihood ratio test is presented for comparing rates of evolutionary change in the paths of descent leading to two species. The test is compared to previous relative rate tests based on variances of estimated numbers of base substitutions. The likelihood approach allows for different transversion and transition rates, and when these rates are actually different, the likelihood ratio test can be much more powerful than the variance-based tests. For single-parameter mutation models, however, the two tests have similar power. The tests are applied to a set of chloroplast sequences from several species of grasses, and additional indications of significantly different rates leading to barley were found with the likelihood ratio test.

Comparing patterns of nucleotide substitution rates among chloroplast loci using the relative ratio test.

Even when several genetic loci are used in molecular evolutionary studies, each locus is typically analyzed independently of the others. This type of approach makes it difficult to study mechanisms and processes that affect multiple genes. In this work we develop a statistical approach for the joint analysis of two or more loci. The tests we propose examine whether or not nucleotide substitution rates across evolutionary lineages have the same relative proportions at two loci. Theses procedures are applied to 33 genes from the chloroplast genomes of rice, tobacco, pine, and liverwort. With the exception of five clearly distinct loci, we find that synonymous substitution rates tend to change proportionally across genes. We interpret these results to be consistent with a "lineage effect" acting on the entire chloroplast genome. In contrast, nonsynonymous rates do not change proportionally across genes, suggesting that locus-specific evolutionary effects dominate patterns of nonsynonymous substitution.

Constraints on intron evolution in the gene encoding the myosin alkali light chain in Drosophila.

Interspecific comparisons of intron sequences reveal conserved blocks of invariant nucleotides and several other departures from the strictly neutral model of molecular evolution. To distinguish the past action of evolutionary forces in introns known to have regulatory information, we examined nucleotide sequence variation at 991 sites in a random sample of 16 Drosophila melanogaster alleles of the gene encoding the myosin alkali light chain (Mlc1). The Mlc1 gene of D. melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. Analyses of these data reveal that introns 4 and 5, which flank the alternatively spliced exon 5, have reduced levels of both intraspecific polymorphism and interspecific divergence relative to intron 3. No polymorphism was observed in any of the exons examined in D. melanogaster. A genealogical analysis clearly demonstrates the occurrence of intragenic recombination in the ancestral history of Mlc1. Recombination events are estimated to be 13 times more likely than mutation events over the span of the sequenced region. Although there is little evidence for pairwise linkage disequilibrium in the Mlc1 region, higher order disequilibrium does seem to be present in the 5' half of the portion of the gene that was examined. Predictions of the folding free energy of the pre-mRNA reveal that sampled alleles have a significantly higher (less stable) free energy than do randomly permuted sequences. These results are consistent with the hypothesis that introns surrounding an alternatively spliced exon are subjected to additional constraints, perhaps due to specific aspects of secondary structure required for appropriate splicing of the pre-mRNA molecule.

Estimating synonymous and nonsynonymous substitution rates.

Partitioning the total substitution rate into synnonymous and nonsynonymous components is a key aspect of many analyses in molecular evolution. Numerous methods exist for estimating these rates. However, until recently none of the estimation procedures were based on a sound statistical footing. In this paper, the evolutionary model of Muse and Gaut (1994) is used as the basis for two sets of parameters quantifying silent and replacement substitution rates. The parameters are shown to be equal when the four nucleotides are equally frequent and unequal otherwise. Maximum-likelihood estimation of these parameters is described, and the performance of these estimates is compared to that of existing estimation procedures. It is shown that the estimates of Nei and Gojobori (1986) are not unbiased for either set of parameters, although they provide very good estimates for one set as long as sequence divergence is not too high. However, some disturbing properties are found for the Nei and Gojobori estimates. In particular, it is shown that the expected value of the Nei and Gojobori estimate of silent substitution rate is a function of both the silent and replacement substitution rates. The maximum-likelihood estimates have no such problems.

Examining rates and patterns of nucleotide substitution in plants.

Driven by rapid improvements in affordable computing power and by the even faster accumulation of genomic data, the statistical analysis of molecular sequence data has become an active area of interdisciplinary research. Maximum likelihood methods have become mainstream because of their desirable properties and, more importantly, their potential for providing statistically sound solutions in complex data analysis settings. In this chapter, a review of recent literature focusing on rates and patterns of nucleotide substitution rates in the nuclear, chloroplast, and mitochondrial genomes of plants demonstrates the power and flexibility of these new methods. The emerging picture of the nucleotide substitution process in plants is a complex one. Evolutionary rates are seen to be quite variable, both among genes and among plant lineages. However, there are hints, particularly in the chloroplast, that individual factors can have important effects on many genes simultaneously.

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.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

Relative rates of nucleotide substitution in the chloroplast genome.

Coding sequences from maize, rice, tobacco, and liverwort chloroplasts are aligned and subjected to relative rate tests. Results of rate tests suggest that coding sequences from maize and rice are evolving with homogeneous rates of nucleotide substitution while coding sequences from the grass lineages (i.e., maize and rice) are evolving at a faster rate than coding sequences from the tobacco chloroplast. Rate tests also suggest that particular loci evolve at significantly faster rates in grass chloroplast genomes than the tobacco chloroplast genome. These loci encode proteins important to RNA polymerase, the H(+)-ATPase complex, and the ribosomal proteins. Much of the variation at these loci can be attributed to differences in nonsynonymous substitution rates. Taken together, these studies suggest that the chloroplast DNA molecular clock varies both between evolutionary lineages and between protein coding loci.

Comparisons of the nucleotide substitution process among repetitive segments of the alpha- and beta-spectrin genes.

The actin-cross-linking protein spectrin is a prominent component of the membrane cytoskeleton. Spectrin is a tetramer of two antiparallel alphabeta-dimers which share a unique and ancient gene structure. The alpha-spectrin and beta-spectrin genes are composed primarily of tandemly repeated 106-amino-acid segments, each of which forms a triple alpha-helical coiled coil. Both the genes and the repeats themselves are homologous. The two genes are thought to be the result of a gene duplication event, and each gene is the product of duplications of the 106-amino-acid repeats. In this work we compare the process of molecular evolution across the repeated segments of the alpha- and beta-spectrin genes. We find that the alpha-spectrin segments have, for the most part, evolved in a homogeneous fashion, while considerable heterogeneity is found among beta-spectrin segments. Several segments with unique known functions are found to have evolved differently than the others. On the basis of heterogeneity of the evolutionary process, we suggest that at least one repeat has a unique function that has yet to be documented. We also present new statistical methods for comparing the evolutionary process between different regions of DNA sequences.

Length variation and secondary structure of introns in the Mlc1 gene in six species of Drosophila.

A nearly universal feature of intron sequences is that even closely related species exhibit a large number of insertion/deletion differences. The goal of the analysis described here is to test whether the observed pattern of insertion/deletion events in the genealogy of the myosin alkali light chain (Mlc1) gene is consistent with neutrality, and if not, to determine the underlying forces of evolutionary change. Mlc1 pre-mRNA is alternatively spliced, and one constraint is that signals necessary for tissue-specificity of directed splicing must be conserved. If the total length of an intron is functionally constrained, then the distribution of indels on branches of the gene genealogy should reflect a departure from randomness. Here we perform a phylogenetic analysis, inferring ancestral states wherever possible on a phylogeny of 29 alleles of Mlc1 from six species of Drosophila. Observed patterns of indels on the genealogy were compared to those from simulated data, with the result that we cannot reject the null hypothesis of neutrality. A clear departure from a neutral prediction was seen in the excess folding free energy predicted for the introns flanking the alternatively spliced exon. Relative rate tests also suggest a retardation in the rate of Mlc1 sequence evolution in the simulans clade.

Relative rates of nucleotide substitution at the rbcL locus of monocotyledonous plants.

We subjected 35 rbcL nucleotide sequences from monocotyledonous taxa to maximum likelihood relative rate tests and estimated relative differences in rates of nucleotide substitution between groups of sequences without relying on knowledge of divergence times between taxa. Rate tests revealed that there is a hierarchy of substitution rate at the rbcL locus within the monocots. Among the taxa analyzed the grasses have the most rapid substitution rate; they are followed in rate by the Orchidales, the Liliales, the Bromeliales, and the Arecales. The overall substitution rate for the rbcL locus of grasses is over 5 times the substitution rate in the rbcL of the palms. The substitution rate at the third codon positions in the rbcL of the grasses is over 8 times the third position rate in the palms. The pattern of rate variation is consistent with the generation-time-effect hypothesis. Heterogenous rates of substitution have important implications for phylogenetic reconstruction.

Comparisons of the molecular evolutionary process at rbcL and ndhF in the grass family (Poaceae).

We examine rate heterogeneity among evolutionary lineages of the grass family at two plasmid loci, ndhF and rbcL, and we introduce a method to determine whether patterns of rate heterogeneity are correlated between loci. We show both that rates of synonymous evolution are heterogeneous among grass lineages and that are heterogeneity is correlated between loci at synonymous sites. At nonsynonymous sites, the pattern of rate heterogeneity is not correlated between loci, primarily due to an aberrant pattern of rate heterogeneity at nonsynonymous sites of rbcL. We compare patterns of synonymous rate heterogeneity to predictors based on the generation time effect and the speciation rate hypotheses. Although there is some evidence for generation time effects, neither generation time effects nor speciation rates appear to be sufficient to explain patterns of rate heterogeneity in the grass plastid sequences.

Molecular evolution of type 1 serine/threonine protein phosphatases.

Type 1 serine/threonine protein phosphatases (PP1s) play key roles in many cellular processes. To understand the evolutionary relationships among PP1s from various kingdoms and to provide a valid basis to evaluate the structure-function relationships of these phosphatases, 44 PP1 sequences were aligned, revealing a high sequence similarity among PP1 homologs. About one-third of the total amino acids are conserved in all the sequences studied. Most of these conserved amino acids are located within a 270-amino-acid core region. They include most sites critical to the activity and regulation of PP1s based on three-dimensional structural studies of mammalian PP1s. Positional variation analysis using a sliding window approach revealed two variable blocks in the 270-amino-acid core region. The major variable block corresponds to a subdomain composed of three alpha-helices (alphaG, alphaH, and alphaI) and three beta-sheets (beta7, beta8, and beta9). Phylogenetic analyses suggested that plant and animal PP1s form distinct monophyletic groups. The plant PP1 family contains several subgroups that may be older than the monocot-dicot divergence. In the animal PP1 family, different vertebrate isoforms appear to form distinct subgroups. Relative substitution rate studies indicated that plant PP1s are more diverse than animal PP1s, with an average substitution rate 1.5 times as large as that of animal PP1s. The possible involvement of PP1s in the establishment of multicellularity is discussed.

Intragenic duplication and divergence in the spectrin superfamily of proteins.

Many structural, signaling, and adhesion molecules contain tandemly repeated amino acid motifs. The alpha-actinin/spectrin/dystrophin superfamily of F-actin-crosslinking proteins contains an array of triple alpha-helical motifs (spectrin repeats). We present here the complete sequence of the novel beta-spectrin isoform beta(Heavy)-spectrin (beta H). The sequence of beta H supports the origin of alpha- and beta-spectrins from a common ancestor, and we present a novel model for the origin of the spectrins from a homodimeric actin-crosslinking precursor. The pattern of similarity between the spectrin repeat units indicates that they have evolved by a series of nested, nonuniform duplications. Furthermore, the spectrins and dystrophins clearly have common ancestry, yet the repeat unit is of a different length in each family. Together, these observations suggest a dynamic period of increase in repeat number accompanied by homogenization within each array by concerted evolution. However, today, there is greater similarity of homologous repeats between species than there is across repeats within species, suggesting that concerted evolution ceased some time before the arthropod/vertebrate split. We propose a two-phase model for the evolution of the spectrin repeat arrays in which an initial phase of concerted evolution is subsequently retarded as each new protein becomes constrained to a specific length and the repeats diverge at the DNA level. This evolutionary model has general applicability to the origins of the many other proteins that have tandemly repeated motifs.