High Rate Of Nucleotide Substitution Research Articles

Rates of evolutionary change in viruses: patterns and determinants

Understanding the factors that determine the rate at which genomes generate and fix mutations provides important insights into key evolutionary mechanisms. We review our current knowledge of the rates of mutation and substitution, as well as their determinants, in RNA viruses, DNA viruses and retroviruses. We show that the high rate of nucleotide substitution in RNA viruses is matched by some DNA viruses, suggesting that evolutionary rates in viruses are explained by diverse aspects of viral biology, such as genomic architecture and replication speed, and not simply by polymerase fidelity.

Molecular phylogeny of Banza (Orthoptera: Tettigoniidae), the endemic katydids of the Hawaiian Archipelago

The extant endemic katydids (Orthoptera: Tettigoniidae) of the Hawaiian Archipelago include one to three species per high island and a single species on Nihoa, all currently placed in the genus Banza. These acoustic insects provide an excellent opportunity for investigating the evolution of reproductive isolation and speciation, but such studies require an understanding of phylogenetic relationships within the group. We use maximum parsimony, likelihood-based Bayesian inference, and maximum likelihood to infer phylogenetic relationships among these taxa, based on ∼2 kb of mitochondrial cytochrome oxidase I and cytochrome b. Our results strongly support two distinct high island clades: one clade (“Clade I”) composed of species from Kauai, Oahu, Molokai, and Lanai and another clade (“Clade II”) composed of species from Maui and Hawaii ( Banza unica, from Oahu, may be basal to both these clades, but its placement is not well resolved). Within these clades, some inferred relationships are strongly supported, such as the sister status of B. kauaiensis (Kauai) and B. parvula (Oahu) within Clade I, but other relationships remain more ambiguous, such as the relative position of B. brunnea (Maui) within Clade II. Although a detailed reconstruction of the historical biogeography of the Hawaiian katydids is difficult, we use our genetic data combined with the known geological history of the Hawaiian Islands to set limits on plausible historical scenarios for diversification of this group. Beyond these historical biogeographic inferences, our results indicate possible cryptic speciation on both Oahu and Hawaii, as well as what may be unusually high average rates of nucleotide substitution. The present work sets the stage for future genetic and experimental investigations of this group.

Genetic variability of measles virus in acute and persistent infections

RNA viruses have high nucleotide substitution rates, and therefore the potential to mutate rapidly. In the case of vaccine preventable RNA viruses, this may potentially lead to emergence of vaccine escape mutants. The WHO has targeted measles virus (MV) for elimination in many regions, and its genetic variability is monitored to estimate appearance of such mutants. Phylogenetic analysis of partial N or H genes of 230 MV strains circulating in the UK over a 10-year period was performed. Substitution rates in three outbreaks were determined to be 3.9 × 10 −3 to 6.7 × 10 −3 per nucleotide per annum. This is an order of magnitude higher than previously reported for circulating MV. Analysis of virus detected sporadically in the UK between 1992 and 2000 lead to a slightly higher substitution rate of 7.8 × 10 −3 per site per year. Additionally, genetic variability of persistent MV, isolated from subacute sclerosing panencephalitis (SSPE) patients, was investigated and appeared more stable than circulating viruses. Profiles of nucleotide changes in acute and persistent virus were compared. In acute virus, 33% of all mutation events occurred from A-to-G, which contrasts the predominant U-to-C mutations found in persistent infections. Mutations do not seem to be driven by positive selection and no association with known biological functions could be found. We conclude that substitution rates in circulating virus may be higher than in persistent, hypermutated virus and that the high substitution rate of MV may allow evolution of escape. Diversity of circulating strains should be closely monitored in the future.

Relationships among the three major lineages of the Acari (Arthropoda:Arachnida) inferred from small subunit rRNA: paraphyly of the Parasitiformes with respect to the Opilioacariformes and relative rates of nucleotide substitution

We inferred phylogeny among the three major lineages of the Acari (mites) from the small subunit rRNA gene. Our phylogeny indicates that the Opilioacariformes is the sister-group to the Ixodida+Holothyrida, not the Ixodida+Mesostigmata+Holothyrida, as previously thought. Support for this relationship increased when sites with the highest rates of nucleotide substitution, and thus the greatest potential for saturation with nucleotide substitutions, were removed. Indeed, the increase in support (and resolution) was despite a 70% reduction in the number of parsimony-informative sites from 408 to 115. This shows that rather than ‘noisy’ sites having no impact on resolution of deep branches, ‘noisy’ sites have the potential to obscure phylogenetic relationships. The arrangement, Ixodida+Holothyrida+Opilioacariformes, however, may be an artefact of long-branch attraction since relative-rate tests showed that the Mesostigmata have significantly faster rates of nucleotide substitution than other parasitiform mites. Thus, the fast rates of nucleotide substitution of the Mesostigmata might have caused the Mesostigmata to be attracted to the outgroup in our trees. We tested the hypothesis that the high rate of nucleotide substitution in some mites was related to their short generation times. The Acari species that have high nucleotide substitution rates usually have short generation times; these mites also tend to be more active and thus have higher metabolic rates than other mites. Therefore, more than one factor may affect the rate of nucleotide substitution in these mites.

Slow algae, fast fungi: exceptionally high nucleotide substitution rate differences between lichenized fungi Omphalina and their symbiotic green algae Coccomyxa

Omphalina basidiolichens are obligate mutualistic associations of a fungus of the genus Omphalina (the exhabitant) and a unicellular green alga of the genus Coccomyxa (the inhabitant). It has been suggested that symbiotic inhabitants have a lower rate of genetic change compared to exhabitants because the latter are more exposed to abiotic environmental variation and competition from other organisms. In order to test this hypothesis we compared substitution rates in the nuclear ribosomal internal transcribed spacer region (ITS1, 5.8S, ITS2) among fungal species with rates among their respective algal symbionts. To ensure valid comparisons, only taxon pairs (12) with a common evolutionary history were used. On average, substitution rates in the ITS1 portion of Omphalina pairs were 27.5 times higher than rates in the corresponding pairs of Coccomyxa since divergence from their respective ancestor at the base of the Omphalina/Coccomyxa lineage. Substitution rates in the 5.8S and the ITS2 portions were 2.4 and 18.0 times higher, respectively. The highest rate difference (43.0) was found in the ITS1 region. These are, to our knowledge, the highest differences of substitution rates reported for symbiotic organisms. We conclude that the Omphalina model system conforms to the proposed hypothesis of lower substitution rates in the inhabitant, but that the mode of transmission of the inhabitant (vertical versus horizontal) could be a prevailing factor in the regulation of unequal rates of nucleotide substitution between co-evolving symbionts. Our phylogenetic study of Coccomyxa revealed three main lineages within this genus, corresponding to free-living Coccomyxa, individuals isolated from basidiolichens Omphalina and Coccomyxa isolated from ascolichens belonging to the Peltigerales.

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects.

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hemipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one another. We found that the correlation was positive and statistically significant (R2 = 0.73, P = 0.01; Rs = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Isolation of novel olfactory receptor genes in marmosets ( Callithrix): insights into pseudogene formation and evidence for functional degeneracy in non-human primates

Nineteen olfactory receptor (OR) genes were isolated from three OR subfamilies in two species of marmoset ( Callithrix). Olfactory receptor 912-93 has high sequence similarity among marmosets and between marmosets and humans, suggesting strong conservation of function. All of the remaining seventeen OR genes identified from subfamilies 3A and 1E were pseudogenes. Following pseudogene formation, marmoset OR genes in both 1E and 3A subfamilies underwent duplications, indel events and a high rate of nucleotide substitution. These results provide a contrast to previous studies, and show that in spite of the keen olfactory sense of marmosets, they harbour many OR pseudogenes. A high rate of in vitro recombination using Pfu polymerase but not Taq polymerase was confirmed. The rapid molecular evolution of OR pseudogenes suggests that they do not provide a useful source of sequence variation for conversion to intact OR genes over evolutionary timescales. The overall pattern of OR evolution in marmosets is comparable to the ‘birth-and-death’ model of gene family evolution. An unbiased view on the evolutionary timing of the reduction of the functional olfactory repertoire in humans must await more data.

Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria).

Mitochondrial genes have been used extensively in population genetic and phylogeographical analyses, in part due to a high rate of nucleotide substitution in animal mitochondrial DNA (mtDNA). Nucleotide sequences of anthozoan mitochondrial genes, however, are virtually invariant among conspecifics, even at third codon positions of protein-coding sequences. Hence, mtDNA markers are of limited use for population-level studies in these organisms. Mitochondrial gene sequence divergence among anthozoan species is also low relative to that exhibited in other animals, although higher level relationships can be resolved with these markers. Substitution rates in anthozoan nuclear genes are much higher than in mitochondrial genes, whereas nuclear genes in other metazoans usually evolve more slowly than, or similar to, mitochondrial genes. Although several mechanisms accounting for a slow rate of sequence evolution have been proposed, there is not yet a definitive explanation for this observation. Slow evolution and unique characteristics may be common in primitive metazoans, suggesting that patterns of mtDNA evolution in these organisms differ from that in other animal systems.

Origin of sphinx, a young chimeric RNA gene in Drosophila melanogaster.

Non-protein-coding RNA genes play an important role in various biological processes. How new RNA genes originated and whether this process is controlled by similar evolutionary mechanisms for the origin of protein-coding genes remains unclear. A young chimeric RNA gene that we term sphinx (spx) provides the first insight into the early stage of evolution of RNA genes. spx originated as an insertion of a retroposed sequence of the ATP synthase chain F gene at the cytological region 60DB since the divergence of Drosophila melanogaster from its sibling species 2-3 million years ago. This retrosequence, which is located at 102F on the fourth chromosome, recruited a nearby exon and intron, thereby evolving a chimeric gene structure. This molecular process suggests that the mechanism of exon shuffling, which can generate protein-coding genes, also plays a role in the origin of RNA genes. The subsequent evolutionary process of spx has been associated with a high nucleotide substitution rate, possibly driven by a continuous positive Darwinian selection for a novel function, as is shown in its sex- and development-specific alternative splicing. To test whether spx has adapted to different environments, we investigated its population genetic structure in the unique "Evolution Canyon" in Israel, revealing a similar haplotype structure in spx, and thus similar evolutionary forces operating on spx between environments.

Mosaic structure and molecular evolution of the leukotoxin operon (lktCABD) in Mannheimia (Pasteurella) haemolytica, Mannheimia glucosida, and Pasteurella trehalosi.

The mosaic structure and molecular evolution of the leukotoxin operon (lktCABD) was investigated by nucleotide sequence comparison of the lktC, lktB, and lktD genes in 23 Mannheimia (Pasteurella) haemolytica, 6 Mannheimia glucosida, and 4 Pasteurella trehalosi strains. Sequence variation in the lktA gene has been described previously (R. L. Davies et al., J. Bacteriol. 183:1394-1404, 2001). The leukotoxin operon of M. haemolytica has a complex mosaic structure and has been derived by extensive inter- and intraspecies horizontal DNA transfer and intragenic recombination events. However, the pattern of recombination varies throughout the operon and among the different evolutionary lineages of M. haemolytica. The lktA and lktB genes have the most complex mosaic structures with segments derived from up to four different sources, including M. glucosida and P. trehalosi. In contrast, the lktD gene is highly conserved in M. haemolytica. The lktC, lktA, and lktB genes of strains representing the major ovine lineages contain recombinant segments derived from bovine or bovine-like serotype A2 strains. These findings support the previous conclusion that host switching of bovine A2 strains from cattle to sheep has played a major role in the evolution of the leukotoxin operon in ovine strains of M. haemolytica. Homologous segments of donor and recipient alleles are identical, or nearly identical, indicating that the recombinational exchanges occurred relatively recent in evolutionary terms. The 5' and 3' ends of the operon are highly conserved in M. haemolytica, which suggests that multiple horizontal exchanges of the complete operon have occurred by a common mechanism such as transduction. Although the lktA and lktB genes both have complex mosaic structures and high nucleotide substitution rates, the amino acid diversity of LktB is significantly lower than that of LktA due to a higher degree of evolutionary constraint against amino acid replacement. The recombinational exchanges within the leukotoxin operon have had greatest effect on LktA and probably provide an adaptive advantage against the host antibody response by generating novel antigenic variation at surface-exposed sites.

Comparative Analysis of a Mitochondrial DNA Control Region Fragment Amplified from Three Adriatic Flatfish Species and Molecular Phylogenesis of Pleuronectiformes.

: The 5'-end of the mitochondrial control region of three Pleuronectiformes from the Adriatic Sea, Platichthys flesus italicus (Adriatic flounder), Solea vulgaris (common sole), and Solea kleini (Klein's sole), was sequenced and compared with that of six other flatfish species from the families Pleuronectidae and Bothidae. The sequence structures of all flatfishes appear very similar and consist of alternate short segments with low, medium, and high rates of nucleotide substitution. Four conserved 19-bp repeats occur at the beginning of the European and Adriatic flounder sequences. The common occurrence of tandem arrays in fish control regions could be related to a stable secondary structure. Molecular phylogenetic relationships among Pleuronectiformes agree well with previous morphologic data at all taxonomic levels. Molecular analyses could therefore contribute to resolving phylogenetic and taxonomic debates within the Pleuronectiformes.

A molecular phylogeny of apiaceae subfamily Apioideae: evidence from nuclear ribosomal DNA internal transcribed spacer sequences

Phylogenetic relationships among 40 New World and Old World members of Apiaceae subfamily Apioideae, representing seven of the eight tribes and eight of the ten subtribes commonly recognized in the subfamily, were inferred from nucleotide sequence variation in the internal transcribed spacer (ITS) regions of 18‐26S nuclear ribosomal DNA. Although the sequences are alignable, with only 11% of sites excluded from the analyses because of alignment ambiguity, divergence values in pairwise comparisons of unambiguous positions among all taxa were high and ranged from 0.5 to 33.2% of nucleotides in ITS 1 and from 0 to 33.2% of nucleotides in ITS 2. Average sequence divergence across both spacer regions was 18.4% of nucleotides. Phylogenies derived from ITS sequences estimated using neighbor‐joining analysis of substitution rates, and maximum likelihood and parsimony methods give trees of essentially similar topology and indicate that: (1) there is little support for any existing system of classification of the subfamily that is based largely on morphological and anatomical features of the mericarp; (2) there is a major phylogenetic division within the subfamily, with one clade comprising the genusSmyrniumand those taxa belonging to Drude's tribes Dauceae, Scandiceae, and Laserpitieae and the other clade comprising all other examined taxa; and (3) the generaArracacia, Coaxana, Coulterophytum, Enantiophylla, Myrrhidendron, Prionosciadium, andRhodosciadium, all endemic to Mexico and Central America, comprise a clade but their relationships to other New World taxa are equivocal. A phylogeny derived from parsimony analysis of chloroplast DNArpoC1 intron sequences is consistent with, but considerably less resolved than, relationships derived from these ITS regions. This study affirms that ITS sequences are useful for phylogenetic inference among closely related members of Apioideae but, owing to high rates of nucleotide substitution, are less useful in resolving relationships among the more ancestral nodes of the phylogeny.

Initial mutational steps toward loss of opsin gene function in cavefish.

The river-dwelling fishes of Astyanax fasciatus and their derivatives from Pachon and Micos caves in Mexico were collected. By comparing red (r007) and green (g103) opsin gene sequences from the three populations, we found a high frequency of C-->T changes in r007 and g103 of the eyed Micos fishes presumably due to spontaneous cytosine methylation or enzymatic deamination of the CG dinucleotide, and a deletion of 12 consecutive nucleotides in g101 of the blind Pachon fishes. Thus, high rates of nucleotide substitution and C-->T transitions, usually found for pseudogenes, may occur even before they become nonfunctional.

Nucleotide sequence analysis of recent epidemic strains of enterovirus 70.

Nucleotide sequences of the genome of enterovirus 70 (EV70) isolated in Osaka in 1993 were determined, and compared with those of the past epidemic strains. Nucleotide substitution rates in 332 bp, 174 bp, and 178 bp of the genes encoding viral protein (VP)1, VP2, or VP3 were 9.0, 7.5, and 5.6% between Kumoi-2/93 and J670/71, respectively. Likewise, the putative amino acid substitution rates were 1.8, 0, and 0%. It seems that the epidemic strains of EV70 in Japan have been evolving at a constant high nucleotide substitution rate but almost all the substitutions were synonymous.

Extensive DNA sequence conservation throughout the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex includes the four species M. tuberculosis, M. bovis, M. africanum, and M. microti. We sequenced 13 M. tuberculosis complex strains in the 16S-to-23S rDNA internal transcribed spacer (ITS). The ITS has a high rate of nucleotide substitution. Previous reports found three nucleotide substitutions in the ITS between two M. tuberculosis complex strains. In contrast, we found the same ITS sequence in all 13 M. tuberculosis complex strains (including all four species and M. bovis BCG). This finding confirms the conservation of 16S rDNA sequence and the high DNA-DNA relatedness found in previous studies. By the usual criteria, the four species of the M. tuberculosis complex would be considered a single species. In a phylogenetic analysis based on the ITS sequence, the four species of the M. tuberculosis complex were distinct from nontuberculous mycobacteria. The ITS contains at least seven potential sites for oligonucleotide probes with specificity for the M. tuberculosis complex.

Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus).

Intraspecific sequence variation in the D-loop region of mtDNA in white sturgeon (Acipenser transmontanus), a relict North American fish species, was examined in 27 individuals from populations of the Columbia and Fraser rivers. Thirty-three varied nucleotide positions were present in a 462-nucleotide D-loop sequence, amplified using the polymerase chain reaction. Bootstrapped neighbor-joining and maximum-parsimony trees of sequences from 19 haplotypes suggest that the two populations have recently diverged. This is consistent with the hypothesis that the Columbia River, a Pleistocene refugium habitat, was the source of founders for the Fraser River after the last glacial recession. On the basis of a divergence time of 10-12 thousand years ago, the estimated substitution rate of the white sturgeon D-loop region is 1.1-1.3 x 10(-7) nucleotides/site/year, which is comparable to rates for hypervariable sequences in the human D-loop region. Furthermore, the ratio of mean percent nucleotide differences in the D-loop (2.27%) to that in whole mtDNA (0.54%, as estimated from restriction-enzyme data) is 4.3, which is similar to the fourfold-to-fivefold-higher substitution rate estimated for the human D-loop. The high nucleotide substitution rate of the hypervariable region indicates that the vertebrate D-loop has potential as a genetic marker in molecular population studies.

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.

Patterns of nucleotide substitutions and implications for the immunological diversity of human immunodeficiency virus

Human immunodeficiency virus (HIV) exhibits immunological hypervariability, which has been an obstacle to successful production of effective anti-HIV vaccines. In this study, we estimated patterns of nucleotide and amino acid substitutions in the env gene of HIVs, with the aim of finding characteristics of the mechanism which generates the immunological diversity of the env protein of HIVs. We found that nucleotide changes between A and G are predominant compared to those between other nucleotides. Since this feature is consistent with the pattern of nucleotide substitutions of other retroviral genes but is quite different from those of most eukaryotic genes, a high rate of nucleotide substitution between A and G appears to be specific for retroviruses including HIVs. We discuss the biological relationship between this biased substitution and the mechanism generating hypervariability of epitopes on the env protein of HIVs.

Molecular Mechanisms of Mitochondrial-Genome Evolution in Higher Plants

Mitochondria are organelles that carry out oxidative phosphorylation in all eukaryotic cells. An independent genetic system is contained within mitochondria that differs in origin and function from its nuclear and chloroplast counterparts. The mitochondrial DNA (mtDNA) codes for some of the enzymes that carry out oxidative phosphorylation and for some mitochondrial components needed to synthesize these proteins. Many of the proteins in the mitochondria are coded by the nucleus. These features of mitochondrial genomes are common to all eukaryotes and have been maintained over long periods of evolutionary time. The mitochondria of higher plants differ from those of animals and fungi in that they share the cytoplasm with the chloroplasts; consequently, they have an evolutionary history that reflects different selective pressures, genetic interactions, and mechanisms of molecular evolution. Plants contain the largest mitochondrial genomes of all the eukaryotes. The mtDNA of higher plants is more variable in structure and organization than the DNA of chloroplasts. One major cause of this variation is DNA rearrangement within and between mtDNA molecules. Rearrangement creates gene duplications, novel peptides, and mitochondrial chromosomes that differ greatly in size. Mitochondria also contain populations of small diverse linear and circular supernumerary DNA's. Compared with the high rate of nucleotide substitution found in mammalian mtDNA, the rate for plant mtDNA is low. Comparison of nucleotide sequences from plant mtDNA with those of eubacteria provides support for the hypothesis that mitochondria originated as endosymbionts.

Higher rates of nucleotide substitution in Drosophila than in mammals.

To examine whether the rate of nucleotide substitution is affected by generation time of the organism, I attempted to estimate an accurate rate of synonymous (silent) substitution in Drosophila lineages, using alcohol dehydrogenase (Adh) and heat shock protein 82 (hsp82) genes. The results obtained suggest that the rate of synonymous substitution in Drosophila lineages is roughly 10-8 per site per year. This rate is approximately two times higher than that of rodents and ten times greater than higher primates. The higher rate in Drosophila may be explained by the shorter generation times of the Drosophila species, though the possibility that the mutation mechanism in Drosophila may differ from that in mammals cannot be excluded.