Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

Comparing Evolutionary Patterns and Variability in the Mitochondrial Control Region and Cytochrome <Emphasis Type="Italic">b</Emphasis> in Three Species of Baleen Whales

The rapidly evolving mitochondrial control region remains an important source of information on phylogeography and demographic history for cetaceans and other vertebrates, despite great uncertainty in the rate of nucleotide substitution across both nucleotide positions and lineages. Patterns of variation in linked markers with slower rates of evolution can potentially be used to calibrate the rate of nucleotide substitution in the control region and to better understand the interplay of evolutionary and demographic forces across the mitochondrial genome above and below the species level. We have examined patterns of diversity within and between three baleen whale species (gray, humpback, and Antarctic minke whales) in order to determine how patterns of molecular evolution differ between cytochrome b and the control region. Our results show that cytochrome b is less variable than expected given the diversity in the control region for gray and humpback whales, even after functional differences are taken into account, but more variable than expected for minke whales. Differences in the frequency distributions of polymorphic sites and in best-fit models of nucleotide substitution indicate that these patterns may be the result of hypervariability in the control region in gray and humpback whales but, in minke whales, may result from a large, stable or expanding population size coupled with saturation at the control region. Using paired cytochrome b and control region data across individuals, we show that the average rate of nucleotide substitution in the control region may be on average 2.6 times higher than phylogenetically derived estimates in cetaceans. These results highlight the complexity of making inferences from control region data alone and suggest that applying simple rules of DNA sequence analyses across species may be difficult.     

Nucleotide distribution and the recognition of coding regions in DNA sequences: an information theory approach

A method for the recognition of coding-regions along DNA sequences is described. The method is based on the observation, made in several cases, that nucleotide distribution at the third position of the codon is more biased (less random) than that in the other two positions. It is suggested that since nucleotide distribution at the third position is only weakly influenced by the amino acid distribution in the coded protein, there must be some constraints at the DNA level which bias the nucleotide distribution at the third position. The distinction between DNA-level constraints and protein-level constraints is discussed in the frame of Information Theory, and the analysis of the Mitochondrial gene coding for subunit-1 of the yeast cytochrome oxidase is presented.     

Accelerated molecular evolution in Microtus (Rodentia) as assessed via complete mitochondrial genome sequences

Microtus is one of the most taxonomically diverse mammalian genera, including over 60 extant species. These rodents have evolved rapidly, as the genus originated less than 2 million years ago. If these numbers are taken at face value, then an average of 30 microtine speciation events have occurred every million years. One explanation for the rapid rate of cladogenesis in Microtus could be the karyotypic differentiation exhibited across the genus: diploid numbers range from 17 to 64. Despite the striking chromosomal variability within Microtus, phenotypic variation is unremarkable. To determine whether nucleotide substitution rates are also elevated in voles, we sequenced the entire mitochondrial DNA (mtDNA) genome of the Eurasian sibling vole (Microtus rossiaemeridionalis). We compared this genome to another previously sequenced vole mtDNA genome (Microtus kikuchii) and performed pairwise sequence comparisons with the mtDNA genomes of ten additional mammalian genera. We found that microtine mtDNA genomes are evolving more rapidly than any other mammalian lineage we sampled, as gauged by the rate of nucleotide substitution across the entire mtDNA genome as well as at each individual protein-coding gene. Additionally, we compared substitution rates within the cytochrome b gene to seven other rodent genera and found that Microtus mtDNA is evolving fastest. The root cause of accelerated evolution in Microtus remains uncertain, but merits further investigation. The mitochondrial genome sequence from this article has been deposited with the GenBank database under accession number DQ015676.     

Estimating the total number of nucleotide substitutions since the common ancestor of a pair of homologous genes: Comparison of several methods and three beta hemoglobin messenger RNA's

Summary The mRNA sequences of beta hemoglobin for human, mouse and rabbit were examined. Observations included the following: (1) there is a significant bias against the use of codons only one nucleotide different from terminating codons; (2) less than 4% of the codons end in adenine; (3), guanine is the most common third position nucleotide but it never follows a second position cytosine; (4) nearest neighbor (doublet) nucleotides are non-random with the greatest contributor to non-randomness being the third position suggesting that codon choice for a given amino acid rather than a choice among amino acids is the more important contributor; (5) the CG dinucleotide is even rarer in positions other than the first and second of the codon than it is in those two, suggesting that the need for arginine has in fact elevated the CG frequency in those positions; (6) 77 per cent of the nucleotides are unsubstituted among these three taxa, which could be a sampling effect, but there is strong evidence that about one-third of them are in fact unsubstitutable because of selective constrainsts; (7) the two longest stretches of unsubstituted nucleotides (32 and 35 consecutive nucleotides) surround the points of the two non-coding insertion sequences; (8) over half the substitutions occur in the third nucleotide position of the codons; (9) silent (non-amino acid changing) substitutions occur at about four times the rate of non-silent substitutions on the basis of their relative opportunity to occur; (10) silent substitutions occur slightly but significantly more often in codons that also have non-silent substitutions than independence of the two events would predict; (11) substitutions occur in adjacent nucleotides significantly more often than chance would predict; (12) among four-fold degenerate codons, third position transitions (principally cytosine-uracil interchanges) outnumber transversions by two to one although the reverse ratio would be expected.The analysis of these messengers provided an opportunity to evaluate the random evolutionary hit (REH) theory. I observed that: (1) the REH theory is premised upon five assumptions, all false; (2) the theory leads to contradictory estimates of the number of varions; (3) the REH values are underestimates; (4) the REH values frequently violate the triangle inequality; (5) the REH values, contrary to claim, are not concordant either with accepted point mutations (PAMs) or augmented distances; (6) the REH values are more likely than values uncorrected for multiple substitutions to give incorrect phylogenies; and (7) the REH values have statistical problems probably associated with a large variance in its fundamental parameter, r_e. From this I conclude that REH theory is not suitable for its intended purpose of estimating from protein sequences of nucleotide substitutions since the common ancestor of two gene products.     

The analysis of intron data and their use in the detection of short signals

Summary In order to examine whether certain short DNA sequences (putative splice signals) occurred in a certain region of an intron more often than would be expected by chance, intron data were examined to see what structure they took. There were significant departures from equal nucleotide frequency, and successive nucleotides clearly did not occur independently in the rat and mouse introns examined. The nonindependence was mainly due to a CG shortage and a less marked TA shortage. However the pairwise frequencies explained almost all the variability in triplet frequencies in the data and so the data could be approximately modeled by using nucleotide frequencies conditional on what the previous nucleotide was. Some coding DNA was also examined and the pairs in second and third positions, and third and first positions in a codon, showed similar departures from independence to those of the intron data. Using the probability model derived for intron data, expected frequencies of putative signals were derived and compared with the observed frequencies.Some of the work for this paper was done while the author was at the Department of Applied Statistics, University of Reading, England     

Evolution of proteins in mammalian cytoplasmic and mitochondrial ribosomes

Summary The proteins of cytoplasmic and mitochondrial ribosomes from the cow and the rat were analyzed by co-electrophoresis in two dimensional polyacrylamide gels to determine their relative evolutionary rates. In a pairwise comparison of individual ribosomal proteins (r-proteins) from the cow and the rat, over 85% of the cytoplasmic r-proteins have conserved electrophoretic properties in this system, while only 15% of the proteins of mitochondrial ribosomes from these animals fell into this category. These values predict that mammalian mitochondrial r-proteins are evolving about 13 times more rapidly than cytoplasmic r-proteins. Based on actual evolutionary rates for representative cytoplasmic r-proteins, this mitochondrial r-protein evolutionary rate corresponds to an amino acid substitution rate of 40×10^–10 per site per year, placing mitochondrial r-proteins in the category of rapidly evolving proteins. The mitochondrial r-proteins are apparently evolving at a rate comparable to that of the mitochondrial rRNA, suggesting that functional constraints act more or less equally on both kinds of molecules in the ribosome. It is significant that mammalian mitochondrial r-proteins are evolving more rapidly than cytoplasmic r-proteins in the same cell, since both sets of r-proteins are encoded by nuclear genes. Such a difference in evolutionary rates implies that the functional constraints operating on ribosomes are somewhat relaxed for mitochondrial ribosomes.Presented at the FEBS Symposium on Genome Organization and Evolution, held in Crete, Greece, September 1–5, 1986     

Variability of Evolutionary Rates of DNA

A statistical analysis of DNA sequences from four nuclear loci and five mitochondrial loci from different orders of mammals is described. A major aim of the study is to describe the variation in the rate of molecular evolution of proteins and DNA. A measure of rate variability is the statistic R, the ratio of the variance in the number of substitutions to the mean number. For proteins, R is found to be in the range 0.16 less than R less than 35.55, thus extending in both directions the values seen in previous studies. An analysis of codons shows that there is a highly significant excess of double substitutions in the first and second positions, but not in the second and third or first and third positions. The analysis of the dynamics of nucleotide evolution showed that the ergodic Markov chain models that are the basis of most published formulas for correcting for multiple substitutions are incompatible with the data. A bootstrap procedure was used to show that the evolution of the individual nucleotides, even the third positions, show the same variation in rates as seen in the proteins. It is argued that protein and silent DNA evolution are uncoupled, with the evolution at both levels showing patterns that are better explained by the action of natural selection than by neutrality. This conclusion is based primarily on a comparison of the nuclear and mitochondrial results.     

Divergent protein coding regions in otherwise closely related androgen-regulated mRNAs

Rat seminal vesicles serve as a model system for studying androgen action. We have sequenced and compared full length cDNAs for two major proteins (S and F) synthesised and secreted under hormonal control. Overall, mRNAS and mRNAF share 57% nucleotide sequence homology suggesting that their genes arose by duplication of a common ancestor. However, the mRNAs display a highly unusual regional distribution of sequence homology, with the untranslated regions (UTRs) being substantially more homologous than the protein-coding regions (PCRs). Detailed analysis of nucleotide substitutions at synonymous and replacement sites shows that the PCRs have evolved very rapidly. Evolutionary conservation of the UTRs is no higher than that of UTRs generally and thus provides no evidence of a specific regulatory role for the UTRs in androgen action. The primary sequences of proteins S and F have diverged so rapidly that they are the best examples of neutrally evolving proteins for which comparative nucleotide sequence data are available. However, despite their rapid divergence, the predicted higher order structures for both proteins consist largely of non-regular conformation. This is discussed in terms of their roles as structural components of the rodent copulatory plug.     

The restriction of codon ambiguity on the basis of known variants

Summary The genetic code may be used to formulate the nucleotide sequence of a messenger RNA from the known amino acid sequence of a protein. Unfortunately, the degeneracy of the code means that there will be ambiguity in the nucleotide assignments in a third or more of the positions. A simple procedure is given that utilizes the information of known genetic variants to reduce that ambiguity. Problems associated with silent polymorphism are treated. The human alpha and beta hemoglobins are used to exemplify the technique. A total of 68 nucleotides in the two sequences are thereby made less ambiguous. One reduction leads to a nucleotide inconsistent with the result of the recently published beta hemoglobin sequence.     

Molecular Evolution of the Plant R Regulatory Gene Family

Anthocyanin pigmentation patterns in different plant species are controlled in part by members of the myc-like R regulatory gene family. We have examined the molecular evolution of this gene family in seven plant species. Three regions of the R protein show sequence conservation between monocot and dicot R genes. These regions encode the basic helix-loop-helix domain, as well as conserved N-terminal and C-terminal domains; mean replacement rates for these conserved regions are 1.02 X 10(-9) nonsynonymous nucleotide substitutions per site per year. More than one-half of the protein, however, is diverging rapidly, with nonsynonymous substitution rates of 4.08 X 10(-9) substitutions per site per year. Detailed analysis of R homologs within the grasses (Poaceae) confirm that these variable regions are indeed evolving faster than the flanking conserved domains. Both nucleotide substitutions and small insertion/deletions contribute to the diversification of the variable regions within these regulatory genes. These results demonstrate that large tracts of sequence in these regulatory loci are evolving at a fairly rapid rate.     

Exploring patterns and extent of bias in estimating divergence time from mitochondrial DNA sequence data in a particular lineage: a case study of salamanders (order Caudata)

In the practice of molecular dating, substitution saturation will bias the results if not properly modeled. Date estimates based on commonly used mitochondrial DNA sequences likely suffer from this problem because of their high substitution rate. Nevertheless, the patterns and extent of such expected bias remain unknown for many major evolutionary lineages, which often differ in ages, available calibrations, and substitution rates of their mitochondrial genome. In this case study of salamanders, we used estimates based on multiple nuclear exons to assess the effects of saturation on dating divergences using mitochondrial genome sequences on a timescale of ~200-300 My. The results indicated that, due to saturation for older divergences and in the absence of younger effective calibration points, dates derived from the mitochondrial data were considerably overestimated and systematically biased toward the calibration point for the ingroup root. The overestimate might be as great as 3-10 times (about 20 My) older than actual divergence dates for recent splitting events and 40 My older for events that are more ancient. For deep divergences, dates estimated were strongly compressed together. Furthermore, excluding the third codon positions of protein-coding genes or only using the RNA genes or second codon positions did not considerably improve the performance. In the order Caudata, slowly evolving markers such as nuclear exons are preferred for dating a phylogeny covering a relatively wide time span. Dates estimated from these markers can be used as secondary calibrations for dating recent events based on rapidly evolving markers for which mitochondrial DNA sequences are attractive candidates due to their short coalescent time. In other groups, similar evaluation should be performed to facilitate the choice of markers for molecular dating and making inferences from the results.     

Allele-mining of rice blast resistance genes at AC134922 locus

The AC134922 locus is one of the most rapidly evolving nucleotide binding site-leucine-rich repeat (NBS-LRR) gene family in rice genome. Six rice blast resistance (R) genes have been cloned from this locus and other two resistance candidate genes, Pi34 and Pi47, are also mapped to this complex locus. Therefore, it seems that more functional R genes could be identified from this locus. In this study, we cloned 22 genes from 12 cultivars based on allele-mining strategy at this locus and identified 6 rice blast R genes with 4 of them recognizing more than one isolates. Our result suggests that gene stacking might be the evolutionary strategy for complex gene locus to interact with rapidly evolving pathogens, which might provide a potential way for the cloning of durable resistance genes. Moreover, the mosaic structure and ambiguous ortholog/paralog relationships of these homologous genes, caused by frequent recombination and gene conversion, indicate that multiple alleles of this complex locus may serve as a reservoir for the evolutionary novelty of these R genes.     

An Improved Codon Modeling Approach for Accurate Estimation of the Mutation Bias

Phylogenetic codon models are routinely used to characterize selective regimes in coding sequences. Their parametric design, however, is still a matter of debate, in particular concerning the question of how to account for differing nucleotide frequencies and substitution rates. This problem relates to the fact that nucleotide composition in protein-coding sequences is the result of the interactions between mutation and selection. In particular, because of the structure of the genetic code, the nucleotide composition differs between the three coding positions, with the third position showing a more extreme composition. Yet, phylogenetic codon models do not correctly capture this phenomenon and instead predict that the nucleotide composition should be the same for all three positions. Alternatively, some models allow for different nucleotide rates at the three positions, an approach conflating the effects of mutation and selection on nucleotide composition. In practice, it results in inaccurate estimation of the strength of selection. Conceptually, the problem comes from the fact that phylogenetic codon models do not correctly capture the fixation bias acting against the mutational pressure at the mutation–selection equilibrium. To address this problem and to more accurately identify mutation rates and selection strength, we present an improved codon modeling approach where the fixation rate is not seen as a scalar, but as a tensor. This approach gives an accurate representation of how mutation and selection oppose each other at equilibrium and yields a reliable estimate of the mutational process, while disentangling the mean fixation probabilities prevailing in different mutational directions.     

Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.     

Pathways of purine nucleotide synthesis in the human trophoblast early and late in gestation

Purine nucleotide and nucleic acid synthesis were studied in cultures of human first and third trimester trophoblastic cells. De novo synthesis was measured as incorporation of 14C-formate into purine nucleotides. Reutilization of purine bases was evaluated by the incorporation of 14C-adenine and -hypoxanthine. Utilization of 14C-adenine was also studied. The incorporation of formate was significantly (P less than 0.01) less active in the third trimester cells. Adenine incorporation was an order of magnitude higher than that of formate in both first and third trimester cells, and significantly (P less than 0.001) higher in the first than third trimester cells. No change in the reutilization of hypoxanthine was observed as a function of gestational age, and the rate was not increased by high extracellular inorganic phosphate. Both phosphorylation and deamination of adenosine increased as a function of concentration up to at least 60 microM, and the high concentration was more efficiently utilized in the first trimester cells. The major pathways of purine nucleotide synthesis function in the human trophoblast throughout gestation, but the contribution of reutilization seems larger than that of de novo synthesis. First trimester trophoblast appears more active in nucleotide and nucleic acid synthesis. Hypoxanthine utilization appears not be enhanced by increased extracellular phosphate. Hypoxanthine may be the major precursor utilized in trophoblastic purine nucleotide synthesis.     

Comparison of responses by bacteriophages and bacteria to pressures on the base composition of open reading frames

Differences in the base composition of genomes can occur because of GC pressure, purine-loading pressure (AG pressure) and RNY pressure, for which there are possible functional explanations, and because of the more abstract pressures exerted by individual bases. The graphical approach of Muto and Osawa was used to analyse how bacteriophages and bacteria balance potentially conflicting pressures on their genomes. Phages generally respond to AG pressure by increasing A while keeping T constant, and by decreasing C while keeping G constant. In contrast, bacteria generally increase both A and T, the former more so, and decrease both G and C, the latter more so. These differences largely occur at third codon positions, which are more responsive than first and second codon positions to AG pressure and GC pressure. Phages respond to AG pressure more in the third codon position than bacteria, whereas bacteria respond more in the first codon position than phages. Conversely, bacteria respond to GC pressure more in the third codon position than phages, whereas phages respond more in the first codon position than bacteria. As GC pressure increases, A is traded for C and AG pressure decreases; first and second codon positions, having more A than T, are most responsive to this negative effect of increased GC pressure; third positions either do not respond (phages) or respond weakly (bacteria). In a set of 48 phage-host pairs, degrees of purine loading were less correlated between phage and host than were GC percentages. These results suggest that pressures on conventional and genome phenotypes operate differentially in phages and bacteria, generating both general differences in base composition and specific differences characteristic of particular phage-host pairs. The reciprocal relationship between GC pressure and AG pressure implies that effects attributed to GC pressure may actually be due to AG pressure, and vice versa.