Evolution of hominoid mitochondrial DNA with special reference to the silent substitution rate over the genome

Summary Focusing on the synonymous substitution rate, we carried out detailed sequence analyses of hominoid mitochondrial (mt) DNAs of ca. 5-kb length. Owing to the outnumbered transitions and strong biases in the base compositions, synonymous substitutions in mtDNA reach rapidly a rather low saturation level. The extent of the compositional biases differs from gene to gene. Such changes in base compositions, even if small, can bring about considerable variation in observed synonymous differences and may result in the region-dependent estimate of the synonymous substitution rate. We demonstrate that such a region dependency is due to a failure to take proper account of heterogeneous compositional biases from gene to gene but that the actual synonymous substitution rate is rather uniform. The synonymous substitution rate thus estimated is 2.37 ± 0.11 × 10^–8 per site per year and comparable to the overall rate for the noncoding region. On the other hand, the rate of nonsynonymous substitutions differs considerably from gene to gene, as expected under the neutral theory of molecular evolution. The lowest rate is 0.8 × 10^–9 per site per year forCOI and the highest rate is 4.5 × 10^–9 forATPase 8, the degree of functional constraints (measured by the ratio of the nonsynonymous to the synonymous substitution rate) being 0.03 and 0.19, respectively. Transfer RNA (tRNA) genes also show variability in the base contents and thus in the nucleotide differences. The average rate for 11 tRNAs contained in the 5-kb region is 3.9 × 10^–9 per site per year. The nucleotide substitutions in the genome suggest that the transition rate is about 17 times faster than the transversion rate.     

The Spontaneous Mutation Rate in the Fission Yeast Schizosaccharomyces pombe

The rate at which new mutations arise in the genome is a key factor in the evolution and adaptation of species. Here we describe the rate and spectrum of spontaneous mutations for the fission yeast Schizosaccharomyces pombe, a key model organism with many similarities to higher eukaryotes. We undertook an ∼1700-generation mutation accumulation (MA) experiment with a haploid S. pombe, generating 422 single-base substitutions and 119 insertion-deletion mutations (indels) across the 96 replicates. This equates to a base-substitution mutation rate of 2.00 × 10⁻¹⁰ mutations per site per generation, similar to that reported for the distantly related budding yeast Saccharomyces cerevisiae. However, these two yeast species differ dramatically in their spectrum of base substitutions, the types of indels (S. pombe is more prone to insertions), and the pattern of selection required to counteract a strong AT-biased mutation rate. Overall, our results indicate that GC-biased gene conversion does not play a major role in shaping the nucleotide composition of the S. pombe genome and suggest that the mechanisms of DNA maintenance may have diverged significantly between fission and budding yeasts. Unexpectedly, CpG sites appear to be excessively liable to mutation in both species despite the likely absence of DNA methylation.     

Rates of Transition and Transversion in Coding Sequences since the Human-Rodent Divergence

Protein-coding sequences of 337 human genes were compared with those of homologous genes from rodent (mouse or rat). A composite alignment containing 477,189 nucleotide positions was constructed, and 21,570 amino acid replacements were inferred. The rates of transitional and transversional silent substitutions in fourfold degenerate sites are estimated as 1.71 × 10^-9 and 1.22 × 10^-9 site^-1 year^-1, respectively. Rates of substitutions in replacement sites, subject to selective constraints mediated by the genetic code, are lower, but also reflect a transitional bias. The amino acid exchange rejected least often during evolution is Asp/Glu, which is fixed at 30% the rate of transversions in silent sites. The most mutable amino acids in this survey are threonine and serine; serine coded by AGY is more mutable than serine coded by TCN. A scoring matrix for evaluating amino acid similarity was derived from this study.     

Rates of synonymous substitution in plant nuclear genes

Summary The rate of synonymous nucleotide substitution in nuclear genes of higher plants has been estimated. The rate varies among genes by a factor of up to two, in a manner that is not immediately explicable in terms of base composition or codon usage bias. The average rate, in both monocots and dicots, is about four times higher than that in chloroplast genes. This leads to an estimated absolute silent substitution rate of 6 × 10^–9 substitutions per site per year that falls within the range of average rates (2–8 × 10^–9) seen in different mammalian nuclear genomes.     

Microevolution of Helicobacter pylori during Prolonged Infection of Single Hosts and within Families

Our understanding of basic evolutionary processes in bacteria is still very limited. For example, multiple recent dating estimates are based on a universal inter-species molecular clock rate, but that rate was calibrated using estimates of geological dates that are no longer accepted. We therefore estimated the short-term rates of mutation and recombination in Helicobacter pylori by sequencing an average of 39,300 bp in 78 gene fragments from 97 isolates. These isolates included 34 pairs of sequential samples, which were sampled at intervals of 0.25 to 10.2 years. They also included single isolates from 29 individuals (average age: 45 years) from 10 families. The accumulation of sequence diversity increased with time of separation in a clock-like manner in the sequential isolates. We used Approximate Bayesian Computation to estimate the rates of mutation, recombination, mean length of recombination tracts, and average diversity in those tracts. The estimates indicate that the short-term mutation rate is 1.4×10⁻⁶ (serial isolates) to 4.5×10⁻⁶ (family isolates) per nucleotide per year and that three times as many substitutions are introduced by recombination as by mutation. The long-term mutation rate over millennia is 5–17-fold lower, partly due to the removal of non-synonymous mutations due to purifying selection. Comparisons with the recent literature show that short-term mutation rates vary dramatically in different bacterial species and can span a range of several orders of magnitude.     

Molecular analysis of hprt mutations induced by chromium picolinate in CHO AA8 cells

Chromium picolinate (CrPic) is a popular dietary supplement, marketed to the public for weight loss, bodybuilding, and control of blood sugar. Recommendations for long-term use at high dosages have led to questions regarding its safety. Previous studies have reported that CrPic can cause chromosomal aberrations and mutations. The purpose of the current work was to compare the mutagenicity of CrPic as a suspension in acetone versus a solution in DMSO, and to characterize the hprt mutations induced by CrPic in CHO AA8 cells. Treatments of 2% acetone or 2% DMSO alone produced no significant increase in 6-thioguanine (6-TG)-resistant mutants after 48 h exposures. Mutants resistant to 6-TG were generated by exposing cells for 48 h to 80 μg/cm² CrPic in acetone or to 1.0 mM CrPic in DMSO. CrPic in acetone produced an average induced mutation frequency (MF) of 56 per 10⁶ surviving cells relative to acetone solvent. CrPic in acetone was 3.5-fold more mutagenic than CrPic in DMSO, which produced an MF of 16.2. Characterization of 61 total mutations in 48 mutants generated from exposure to CrPic in acetone showed that base substitutions comprised 33% of the mutations, with transversions being predominant; deletions made up 62% of the mutations, with one-exon deletions predominating; and 1–4 bp insertions made up 5% of the characterized mutations. CrPic induced a statistically greater number of deletions and a statistically smaller number of base substitutions than have been measured in spontaneously generated mutants. These data confirm previous studies showing that CrPic is mutagenic, and support the contention that further study is needed to verify the safety of CrPic for human consumption.     

Estimation of the Spontaneous Mutation Rate per Nucleotide Site in a Drosophila melanogaster Full-Sib Family

We employed deep genome sequencing of two parents and 12 of their offspring to estimate the mutation rate per site per generation in a full-sib family of Drosophila melanogaster recently sampled from a natural population. Sites that were homozygous for the same allele in the parents and heterozygous in one or more offspring were categorized as candidate mutations and subjected to detailed analysis. In 1.23 × 10⁹ callable sites from 12 individuals, we confirmed six single nucleotide mutations. We estimated the false negative rate in the experiment by generating synthetic mutations using the empirical distributions of numbers of nonreference bases at heterozygous sites in the offspring. The proportion of synthetic mutations at callable sites that we failed to detect was <1%, implying that the false negative rate was extremely low. Our estimate of the point mutation rate is 2.8 × 10⁻⁹ (95% confidence interval = 1.0 × 10⁻⁹ − 6.1 × 10⁻⁹) per site per generation, which is at the low end of the range of previous estimates, and suggests an effective population size for the species of ∼1.4 × 10⁶. At one site, point mutations were present in two individuals, indicating that there had been a premeiotic mutation cluster, although surprisingly one individual had a G→A transition and the other a G→T transversion, possibly associated with error-prone mismatch repair. We also detected three short deletion mutations and no insertions, giving a deletion mutation rate of 1.2 × 10⁻⁹ (95% confidence interval = 0.7 × 10⁻⁹ − 11 × 10⁻⁹).     

Variation in synonymous substitution rates among mammalian genes and the correlation between synonymous and nonsynonymous divergences

Using mammalian gene sequences, the variances in the numbers of synonymous and nonsynonymous substitutions among genes were estimated together with the correlation coefficient between the two. The expected correlation coefficient can be obtained under the neutral theory using these estimated values of the variances. The expected coefficient is found to often be one-half to two-thirds of the observed value. Possible causes for the disagreement were discussed, such as correlated selective constraints on the two types of substitutions and excess doublet mutations. The variance of mutation rate and that of selective constraint were also estimated. The results show that the coefficient of variation of the former is 0.2–0.3, whereas that of the latter is 0.7–0.9. Correspondence to: T. Ohta     

Chloroplast DNA base substitutions: an experimental assessment

An experimental assessment was carried out to determine directly the frequency and types of spontaneous base substitutions that occur in chloroplast DNA. A target site within the chloroplast 16S rRNA gene of the green alga Chlamydomonas reinhardtii was chosen for the assay. Mutations at this site were known to confer spectinomycin resistance and simultaneously result in the loss of an AatII cleavage site. In the experiments reported here, base substitutions at any individual base occurred at a frequency in the range of 0.9–11 per 10⁹ viable cells plated. Four new mutations that confer resistance to spectinomycin were identified at the target site in the Chlamydomonas chloroplast 16S rRNA gene. When the relative rates of transition and transversion mutations were quantified, a bias toward transversions was observed. The prominence of A/T C/G transversions in the observed mutation spectrum suggests that oxidative damage may be the major cause of base substitution mutations within the chloroplast.     

Avoiding Dangerous Missense: Thermophiles Display Especially Low Mutation Rates

Rates of spontaneous mutation have been estimated under optimal growth conditions for a variety of DNA-based microbes, including viruses, bacteria, and eukaryotes. When expressed as genomic mutation rates, most of the values were in the vicinity of 0.003–0.004 with a range of less than two-fold. Because the genome sizes varied by roughly 10⁴-fold, the mutation rates per average base pair varied inversely by a similar factor. Even though the commonality of the observed genomic rates remains unexplained, it implies that mutation rates in unstressed microbes reach values that can be finely tuned by evolution. An insight originating in the 1920s and maturing in the 1960s proposed that the genomic mutation rate would reflect a balance between the deleterious effect of the average mutation and the cost of further reducing the mutation rate. If this view is correct, then increasing the deleterious impact of the average mutation should be countered by reducing the genomic mutation rate. It is a common observation that many neutral or nearly neutral mutations become strongly deleterious at higher temperatures, in which case they are called temperature-sensitive mutations. Recently, the kinds and rates of spontaneous mutations were described for two microbial thermophiles, a bacterium and an archaeon. Using an updated method to extrapolate from mutation-reporter genes to whole genomes reveals that the rate of base substitutions is substantially lower in these two thermophiles than in mesophiles. This result provides the first experimental support for the concept of an evolved balance between the total genomic impact of mutations and the cost of further reducing the basal mutation rate.     

Germline Mutations of Tetranucleotide DNA Repeats in Families with Normal Children and Reproductive Pathology

We have previously reported a high rate of tetranucleotide DNA repeat mutations, including mutations of both germline and somatic origin, in spontaneous human abortions. To analyze in more detail mutational microsatellite (MS) variability in meiosis and its possible association with disturbed embryonic development, we have conducted a comparative study of mutation rates of a panel of 15 autosomal tetranucleotide MSs in 55 families with healthy children and in 103 families that have had spontaneous abortions with normal karyotypes. In the families with miscarriage, the gametic MS mutation rate was higher than in the families with normal reproductive function (4.36 × 10⁻³ versus 2.32 × 10⁻³ per locus per gamete per generation), but this difference was statistically nonsignificant (P = 0.25). No association of MS mutations with familiar miscarriage was found. Mutations at the MS loci studied were recorded almost 3 times as often in spermatogenesis as in oogenesis, which is likely to result from a greater number of DNA replication cycles in male germline cell precursors than in female ones. Mutations increasing and reducing the MS sequence length appeared at virtually the same rate. Changes in MS DNA sequence length per one repeated element, i.e., single-step mutations (93% of cases) exceeded all other events of allele length change. The highest number of mutations (81.2%) was found in longer alleles. This distribution of mutations by size, direction, and parental origin corresponds to the multistep mutation model of their emergence via mechanism of DNA strand slippage during replication.__________Translated from Genetika, Vol. 41, No. 7, 2005, pp. 943–953.Original Russian Text Copyright © 2005 by Nikitina, Lebedev, Sukhanova, Nazarenko.     

The molecular evolution of pancreatic ribonuclease

Summary The primary structures of pancreatic ribonucleases from 26 species (18 artiodactyls, horse, whale, 5 rodents and turtle) are known. Several species contain identical ribonucleases (cow/bison; sheep/goat), other species show polymorphism (arabian camel) or the presence of two structural gene loci (guinea pig pancreas contains two ribonucleases that differ at 31 positions). 26 different sequences (including the ribonuclease from bovine seminal plasma which is paralogous to the pancreatic ribonucleases) were used to construct a most parsimonious tree. A second tree that most closely approximates current biological opinion requires 402 whereas the most parsimonious tree requires 389 nucleotide substitutions. The artiodactyl part of the most parsimonious tree conforms quite well with the biological one of this order, except for the position of the giraffe which is placed with the pronghorn. Other parts of the most parsimonious tree agree less with the biological tree, probably as a result of the occurrence of many parallel and back substitutions. Bovine seminal ribonuclease was found to be the result of a gene duplication which occurred before the divergence of the true ruminants, but after the divergence of this group from the cameloids.The evolutionary rate of ribonuclease was found to be 390, 3.0 and 11 nucleotide substitutions per 10⁹ yrs per ribonuclease gene, codon and covarion respectively. However, there is much variation in evolutionary rate in different taxa. Values ranging from about 100 (in the bovidae) to about 700 (in the rodents) nucleotide substitutions per 10⁹ yrs per gene were found.A method for counting parallel and back mutations is presented. The 389 nucleotide substitutions in the most parsimonious tree occur at 88 codon positions; 154 of them are the result of parallel and back mutations. Parallel evolution to a similar structure, including the presence of 2 sites with carbohydrate, was demonstrated in an extensive region at the surface of pig and guinea pig ribonuclease B. The presence of carbohydrate probably is important in a number of species. A correlation between the presence of heavily glycosidated ribonucleases and coecal digestion was observed. Hypothetical sequences of ancestral ungulate ribonucleases contain many recognition sites for carbohydrate attachment; this suggests that herbivores with coecal digestion might have preceded the true ruminants in mammalian evolution.     

Quantifying Selection against Synonymous Mutations in HIV-1 env Evolution

Intrapatient evolution of human immunodeficiency virus type 1 (HIV-1) is driven by the adaptive immune system resulting in rapid change of HIV-1 proteins. When cytotoxic CD8⁺ T cells or neutralizing antibodies target a new epitope, the virus often escapes via nonsynonymous mutations that impair recognition. Synonymous mutations do not affect this interplay and are often assumed to be neutral. We test this assumption by tracking synonymous mutations in longitudinal intrapatient data from the C2-V5 part of the env gene. We find that most synonymous variants are lost even though they often reach high frequencies in the viral population, suggesting a cost to the virus. Using published data from SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) assays, we find that synonymous mutations that disrupt base pairs in RNA stems flanking the variable loops of gp120 are more likely to be lost than other synonymous changes: these RNA hairpins might be important for HIV-1. Computational modeling indicates that, to be consistent with the data, a large fraction of synonymous mutations in this genomic region need to be deleterious with a cost on the order of 0.002 per day. This weak selection against synonymous substitutions does not result in a strong pattern of conservation in cross-sectional data but slows down the rate of evolution considerably. Our findings are consistent with the notion that large-scale patterns of RNA structure are functionally relevant, whereas the precise base pairing pattern is not.     

Reduction of synonymous substitutions in the core protein gene of hepatitis C virus

Molecular evolutionary analyses were carried out to elucidate the phylogenetic relationships, the evolutionary rate, and the divergence times of hepatitis C viruses. Using the nucleotide sequences of the viruses isolated from various locations in the world, we constructed phylogenetic trees. The trees showed that strains isolated from a single location were not necessarily clustered as a group. This suggests that the viruses may be transferred with blood on a worldwide scale. We estimated the evolutionary rates at synonymous and nonsynonymous sites for all genes in the viral genome. We then found that the rate (1.35 × 10^–3 per site per year) at synonymous sites for the C gene was much smaller than those for the other genes (e.g., 6.29 × 10^–3 per site per year for the E gene). This indicates that a special type of functional constraint on synonymous substitutions may exist in the C gene. Because we found an open reading frame (ORF) with the C gene region, the possibility exists that synonymous substitutions for the C gene are constrained by the overlapping ORF whose reading frame is different from that of the C gene. Applying the evolutionary rates to the trees, we also suggest that major groups of hepatitis C viruses diverged from their common ancestor several hundred years ago. Correspondence to: T. Gojobori     

Genetic Variability among Complete Human Respiratory Syncytial Virus Subgroup A Genomes: Bridging Molecular Evolutionary Dynamics and Epidemiology

Human respiratory syncytial virus (RSV) is an important cause of severe lower respiratory tract infections in infants and the elderly. In the vast majority of cases, however, RSV infections run mild and symptoms resemble those of a common cold. The immunological, clinical, and epidemiological profile of severe RSV infections suggests a disease caused by a virus with typical seasonal transmission behavior, lacking clear-cut virulence factors, but instead causing disease by modifying the host’s immune response in a way that stimulates pathogenesis. Yet, the interplay between RSV-evoked immune responses and epidemic behavior, and how this affects the genomic evolutionary dynamics of the virus, remains poorly understood. Here, we present a comprehensive collection of 33 novel RSV subgroup A genomes from strains sampled over the last decade, and provide the first measurement of RSV-A genomic diversity through time in a phylodynamic framework. In addition, we map amino acid substitutions per protein to determine mutational hotspots in specific domains. Using Bayesian genealogical inference, we estimated the genomic evolutionary rate to be 6.47×10⁻⁴ (credible interval: 5.56×10⁻⁴, 7.38×10⁻⁴) substitutions/site/year, considerably slower than previous estimates based on G gene sequences only. The G gene is however marked by elevated substitution rates compared to other RSV genes, which can be attributed to relaxed selective constraints. In line with this, site-specific selection analyses identify the G gene as the major target of diversifying selection. Importantly, statistical analysis demonstrates that the immune driven positive selection does not leave a measurable imprint on the genome phylogeny, implying that RSV lineage replacement mainly follows nonselective epidemiological processes. The roughly 50 years of RSV-A genomic evolution are characterized by a constant population size through time and general co-circulation of lineages over many epidemic seasons – a conclusion that might be taken into account when developing future therapeutic and preventive strategies.     

Biochemical analysis of related,independently arising histocompatibility mutants: bm17 and KB-98 enlarge the “bg series” of H-2Kb mutants

The “bg” series of MHC mutations is the most prevalent type of mutations of K^b in C57BL/6 mice screened by reciprocal tail skin grafting. The basis for identification of this series of mutations is the incompatibility of grafts between the parental B6 and the mutant. This series takes the longest to reciprocally reject the skin grafts. The series can be subdivided into “bg 1” and “bg 2” groups based on K^b-restricted recognition of virus-infected mutant target cells. The biochemical basis for these mutations are amino acid substitutions at residues 116 and 121 of the K^b transplantation antigen. These substitutions do not alter monoclonal antibody binding sites. The structural basis of MAb binding and the genetic basis of the mutation are discussed. This study was supported in part by USPHS Grants AI-07289, AI-10702, NCI P30-CA-13330, American Cancer Society Grant IM-236, and American Cancer Society Fellowship PF-2126. Stanley G. Nathenson is a member of the Irvington House Institute for Medical Research.     

Phylogenetic Analysis Reveals Rapid Evolutionary Dynamics in the Plant RNA Virus Genus Tobamovirus

Early studies on the evolutionary dynamics of plant RNA viruses suggested that they may evolve more slowly than their animal counterparts, sometimes dramatically so. However, these estimates were often based on an assumption of virus–host codivergence over time-scales of many millions of years that is difficult to verify. An important example are viruses of the genus Tobamovirus, where the assumption of host–virus codivergence over 100 million years has led to rate estimates in the range of ~1 × 10⁻⁸ nucleotide substitutions per site, per year. Such a low evolutionary rate is in apparent contradiction with the ability of some tobamoviruses to quickly overcome inbred genetic resistance. To resolve how rapidly molecular evolution proceeds in the tobomaviruses, we estimated rates of nucleotide substitution, times to common ancestry, and the extent of congruence between virus and host phylogenies. Using Bayesian coalescent methods applied to time-stamped sequences, we estimated mean evolutionary rates at the nucleotide and amino acid levels of between 1 × 10⁻⁵ and 1.3 × 10⁻³ substitutions per site, per year, and hence similar to those seen in a broad range of animal and plant RNA viruses. Under these rates, a conservative estimate for the time of origin of the sampled tobamoviruses is within the last 100,000 years, and hence a far more recently than proposed assuming codivergence. This is supported by our cophylogeny analysis which revealed significantly discordant evolutionary histories between the tobamoviruses and the plant families they infect.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.