Whole Genome Sequencing of Mutation Accumulation Lines Reveals a Low Mutation Rate in the Social Amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10⁻⁹, with a Poisson confidence interval of 4.1×10⁻⁹ − 9.5×10⁻⁹, per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10⁻¹¹, with a Poisson confidence interval ranging from 7.4×10⁻¹³ to 1.6×10⁻¹⁰, is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.     

Kinetics of Leptin Binding to the Q223R Leptin Receptor

Studies in human populations and mouse models of disease have linked the common leptin receptor Q223R mutation to obesity, multiple forms of cancer, adverse drug reactions, and susceptibility to enteric and respiratory infections. Contradictory results cast doubt on the phenotypic consequences of this variant. We set out to determine whether the Q223R substitution affects leptin binding kinetics using surface plasmon resonance (SPR), a technique that allows sensitive real-time monitoring of protein-protein interactions. We measured the binding and dissociation rate constants for leptin to the extracellular domain of WT and Q223R murine leptin receptors expressed as Fc-fusion proteins and found that the mutant receptor does not significantly differ in kinetics of leptin binding from the WT leptin receptor. (WT: k_a 1.76×10⁶±0.193×10⁶ M⁻¹ s⁻¹, k_d 1.21×10⁻⁴±0.707×10⁻⁴ s⁻¹, K_D 6.47×10⁻¹¹±3.30×10⁻¹¹ M; Q223R: k_a 1.75×10⁶±0.0245×10⁶ M⁻¹ s⁻¹, k_d 1.47×10⁻⁴±0.0505×10⁻⁴ s⁻¹, K_D 8.43×10⁻¹¹±0.407×10⁻¹¹ M). Our results support earlier findings that differences in affinity and kinetics of leptin binding are unlikely to explain mechanistically the phenotypes that have been linked to this common genetic variant. Future studies will seek to elucidate the mechanism by which this mutation influences susceptibility to metabolic, infectious, and malignant pathologies.     

Probability Distribution of Drug-Resistant Mutants in Unselected Populations of Mycobacterium tuberculosis

The fluctuation test shows that Mycobacterium tuberculosis mutates to resistance to isoniazid, streptomycin, ethambutol and rifampin spontaneously and at random. The average mutation rates for the drugs, in the same order, were calculated to be 2.56 × 10⁻⁸, 2.95 × 10⁻⁸, 10⁻⁷, and 2.25 × 10⁻¹⁰ mutation per bacterium per generation. The relatively high mutation rate to ethambutol resistance and the low mutation rate to rifampin resistance were confirmed by analyzing the increase in the proportion of mutants with time in a growing population of the tubercle bacilli. The highest proportions of mutants to be expected in unselected populations of the tubercle bacilli were calculated from the results of fluctuation tests.     

Low base‐substitution mutation rate and predominance of insertion‐deletion events in the acidophilic bacterium Acidobacterium capsulatum

Analyses of spontaneous mutation have shown that total genome‐wide mutation rates are quantitatively similar for most prokaryotic organisms. However, this view is mainly based on organisms that grow best around neutral pH values (6.0–8.0). In particular, the whole‐genome mutation rate has not been determined for an acidophilic organism. Here, we have determined the genome‐wide rate of spontaneous mutation in the acidophilic Acidobacterium capsulatum using a direct and unbiased method: a mutation‐accumulation experiment followed by whole‐genome sequencing. Evaluation of 69 mutation accumulation lines of A. capsulatum after an average of ~2900 cell divisions yielded a base‐substitution mutation rate of 1.22 × 10⁻¹⁰ per site per generation or 4 × 10⁻⁴ per genome per generation, which is significantly lower than the consensus value (2.5−4.6 × 10⁻³) of mesothermophilic (~15–40°C) and neutrophilic (pH 6–8) prokaryotic organisms. However, the insertion‐deletion rate (0.43 × 10⁻¹⁰ per site per generation) is high relative to the base‐substitution mutation rate. Organisms with a similar effective population size and a similar expected effect of genetic drift should have similar mutation rates. Because selection operates on the total mutation rate, it is suggested that the relatively high insertion‐deletion rate may be balanced by a low base‐substitution rate in A. capsulatum, with selection operating on the total mutation rate.     

Y chromosome lineage- and village-specific genes on chromosomes 1p22 and 6q27 control visceral leishmaniasis in Sudan

Familial clustering and ethnic differences suggest that visceral leishmaniasis caused by Leishmania donovani is under genetic control. A recent genome scan provided evidence for a major susceptibility gene on Chromosome 22q12 in the Aringa ethnic group in Sudan. We now report a genome-wide scan using 69 families with 173 affected relatives from two villages occupied by the related Masalit ethnic group. A primary ten-centimorgan scan followed by refined mapping provided evidence for major loci at 1p22 (LOD score 5.65; nominal p = 1.72 × 10⁻⁷; empirical p < 1 × 10⁻⁵; λ_S = 5.1) and 6q27 (LOD score 3.74; nominal p = 1.68 × 10⁻⁵; empirical p < 1 × 10⁻⁴; λ_S = 2.3) that were Y chromosome–lineage and village-specific. Neither village supported a visceral leishmaniasis susceptibility gene on 22q12. The results suggest strong lineage-specific genes due to founder effect and consanguinity in these recently immigrant populations. These chance events in ethnically uniform African populations provide a powerful resource in the search for genes and mechanisms that regulate this complex disease.     

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⁻⁹).     

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.     

Mitogenomic analysis of a 50-generation chicken pedigree reveals a rapid rate of mitochondrial evolution and evidence for paternal mtDNA inheritance

Mitochondrial genomes represent a valuable source of data for evolutionary research, but studies of their short-term evolution have typically been limited to invertebrates, humans and laboratory organisms. Here we present a detailed study of 12 mitochondrial genomes that span a total of 385 transmissions in a well-documented 50-generation pedigree in which two lineages of chickens were selected for low and high juvenile body weight. These data allowed us to test the hypothesis of time-dependent evolutionary rates and the assumption of strict maternal mitochondrial transmission, and to investigate the role of mitochondrial mutations in determining phenotype. The identification of a non-synonymous mutation in ND4L and a synonymous mutation in CYTB, both novel mutations in Gallus, allowed us to estimate a molecular rate of 3.13 × 10⁻⁷ mutations/site/year (95% confidence interval 3.75 × 10⁻⁸–1.12 × 10⁻⁶). This is substantially higher than avian rate estimates based upon fossil calibrations. Ascertaining which of the two novel mutations was present in an additional 49 individuals also revealed an instance of paternal inheritance of mtDNA. Lastly, an association analysis demonstrated that neither of the point mutations was strongly associated with the phenotypic differences between the two selection lines. Together, these observations reveal the highly dynamic nature of mitochondrial evolution over short time periods.     

Leveraging Distant Relatedness to Quantify Human Mutation and Gene-Conversion Rates

The rate at which human genomes mutate is a central biological parameter that has many implications for our ability to understand demographic and evolutionary phenomena. We present a method for inferring mutation and gene-conversion rates by using the number of sequence differences observed in identical-by-descent (IBD) segments together with a reconstructed model of recent population-size history. This approach is robust to, and can quantify, the presence of substantial genotyping error, as validated in coalescent simulations. We applied the method to 498 trio-phased sequenced Dutch individuals and inferred a point mutation rate of 1.66 × 10⁻⁸ per base per generation and a rate of 1.26 × 10⁻⁹ for <20 bp indels. By quantifying how estimates varied as a function of allele frequency, we inferred the probability that a site is involved in non-crossover gene conversion as 5.99 × 10⁻⁶. We found that recombination does not have observable mutagenic effects after gene conversion is accounted for and that local gene-conversion rates reflect recombination rates. We detected a strong enrichment of recent deleterious variation among mismatching variants found within IBD regions and observed summary statistics of local sharing of IBD segments to closely match previously proposed metrics of background selection; however, we found no significant effects of selection on our mutation-rate estimates. We detected no evidence of strong variation of mutation rates in a number of genomic annotations obtained from several recent studies. Our analysis suggests that a mutation-rate estimate higher than that reported by recent pedigree-based studies should be adopted in the context of DNA-based demographic reconstruction.     

Empirical estimates of the mutation rate for an alphabaculovirus

Mutation rates are of key importance for understanding evolutionary processes and predicting their outcomes. Empirical mutation rate estimates are available for a number of RNA viruses, but few are available for DNA viruses, which tend to have larger genomes. Whilst some viruses have very high mutation rates, lower mutation rates are expected for viruses with large genomes to ensure genome integrity. Alphabaculoviruses are insect viruses with large genomes and often have high levels of polymorphism, suggesting high mutation rates despite evidence of proofreading activity by the replication machinery. Here, we report an empirical estimate of the mutation rate per base per strand copying (s/n/r) of Autographa californica multiple nucleopolyhedrovirus (AcMNPV). To avoid biases due to selection, we analyzed mutations that occurred in a stable, non-functional genomic insert after five serial passages in Spodoptera exigua larvae. Our results highlight that viral demography and the stringency of mutation calling affect mutation rate estimates, and that using a population genetic simulation model to make inferences can mitigate the impact of these processes on estimates of mutation rate. We estimated a mutation rate of μ = 1×10⁻⁷ s/n/r when applying the most stringent criteria for mutation calling, and estimates of up to μ = 5×10⁻⁷ s/n/r when relaxing these criteria. The rates at which different classes of mutations accumulate provide good evidence for neutrality of mutations occurring within the inserted region. We therefore present a robust approach for mutation rate estimation for viruses with stable genomes, and strong evidence of a much lower alphabaculovirus mutation rate than supposed based on the high levels of polymorphism observed.     

Steady State Sodium and Rubidium Effluxes in Pisum sativum Roots

Steady state effluxes of potassium and sodium ions were measured on Pisum sativum var. Alaska root segments excised from seedlings which had grown in a nutrient solution containing the major inorganic ions and either ⁸⁶Rb as a tracer for K or ²²Na as a tracer for Na. Fluxes appeared to be from 2 cellular compartments, a small compartment with a high flux rate and a larger compartment with a slow flux rate. Cell wall exchange fluxes are believed to have been negligible. Efflux rates for 11.3% and 88.7% of cellular potassium ions were 6 × 10⁻⁷ and 1.32 × 10⁻⁷ respectively; rates for 33.7% and 66.3% of cellular sodium ions were 1.48 × 10⁻⁷ and 3.83 × 10⁻⁸ respectively, (equivalents per gram fr wt per hr). The sodium flux measurements, with previous measurements of ionic concentrations and transmembrane potentials, support the theory that sodium is transported actively from Pisum roots.     

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.     

The Comparative Genomics of Human Respiratory Syncytial Virus Subgroups A and B: Genetic Variability and Molecular Evolutionary Dynamics

Genomic variation and related evolutionary dynamics of human respiratory syncytial virus (RSV), a common causative agent of severe lower respiratory tract infections, may affect its transmission behavior. RSV evolutionary patterns are likely to be influenced by a precarious interplay between selection favoring variants with higher replicative fitness and variants that evade host immune responses. Studying RSV genetic variation can reveal both the genes and the individual codons within these genes that are most crucial for RSV survival. In this study, we conducted genetic diversity and evolutionary rate analyses on 36 RSV subgroup B (RSV-B) whole-genome sequences. The attachment protein, G, was the most variable protein; accordingly, the G gene had a higher substitution rate than other RSV-B genes. Overall, less genetic variability was found among the available RSV-B genome sequences than among RSV-A genome sequences in a comparable sample. The mean substitution rates of the two subgroups were, however, similar (for subgroup A, 6.47 × 10⁻⁴ substitutions/site/year [95% credible interval {CI 95%}, 5.56 × 10⁻⁴ to 7.38 × 10⁻⁴]; for subgroup B, 7.76 × 10⁻⁴ substitutions/site/year [CI 95%, 6.89 × 10⁻⁴ to 8.58 × 10⁻⁴]), with the time to their most recent common ancestors (TMRCAs) being much lower for RSV-B (19 years) than for RSV-A (46.8 years). The more recent RSV-B TMRCA is apparently the result of a genetic bottleneck that, over longer time scales, is still compatible with neutral population dynamics. Whereas the immunogenic G protein seems to require high substitution rates to ensure immune evasion, strong purifying selection in conserved proteins such as the fusion protein and nucleocapsid protein is likely essential to preserve RSV viability.     

Requirement for Cell Dispersion Prior to Selection of Induced Azaguanine-Resistant Colonies of Chinese Hamster Cells

With V79 Chinese hamster cell cultures treated with a mutagen, the maximum frequency of colonies resistant to 8-azaguanine (AZG) was attained when the cells were dispersed after a suitable expression time before adding the selection medium. V79–4 cells were exposed to 500 µM MMS, 7 µM AFAA, or 10 µM MNNG and allowed to multiply before being reseeded at 4 x 10⁴ cells/60 mm dish and selected with 10 µg/ml AZG. Maximum frequencies of 4 x 10^-5, 4 x 10^-4, and 2.4 x 10^-3 were obtained about 100, 130, and 200 hrs after exposure to MMS, AFAA, and MNNG, respectively. The maximum frequencies following MMS or MNNG treatments were about 10-fold greater than those obtained when induction and selection of AZG-resistant colonies were performed in the same culture dish. The reseeding of treated cells eliminated the possibility of metabolic cooperation within mosaic colonies of wild-type and mutant cells and achieved expression of the induced changes before intercolony crossfeeding reduced the frequency of resistant colonies.—AZG-resistant colonies were selected in medium containing dialyzed fetal bovine serum, and the selection medium was replaced at least twice. Both serum dialysis and selection medium replacement were necessary for consistent achievement of background frequencies of resistant colonies near 10^-6. Reconstruction experiments with AZG-resistant V79 lines showed that the efficiency of recovery of resistant cells in the selection medium was constant over a range of 0–20 colonies observed/dish. A mixed population of V79 and AZG-resistant cells was also correctly analyzed by the procedure used in mutagenesis studies.     

Genetic Effects of Acute Spermatogonial X-Irradiation of the Laboratory Rat

The genetic effects of one generation of spermatogonial X-irradiation in rats, by a single dose of 600r in one experiment and by a fractionated dose of 450r in another, were measured in three generations of their descendants. Estimates of dominant lethal mutation rates—(2 to 3) x 10 ^-4/gamete/r—from litter size differences between irradiated and nonirradiated stock were consistent with previous estimates from rats and mice. Similar consistency was found for estimates of sex-linked recessive mutation rates—(1 to 2) x 10^-4 chromosome/r—from male proportions within strains; however, when measured in crossbreds the proportion of males was higher in the irradiated than in the nonirradiated lines. This inconsistency in results is in keeping with the contradictory results reported for recessive sex-linked lethal mutation rates in mice. The effects used to estimate recessive lethal mutation rates which were unusually high—(2 to 14) x 10^-4/gamete/r—were not significant. Other factors that could have contributed to the observed effects are postulated.     

Rotating Magnetic Field-Assisted Reactor Enhances Mechanisms of Phage Adsorption on Bacterial Cell Surface

Growing interest in bacteriophage research and use, especially as an alternative treatment option for multidrug-resistant bacterial infection, requires rapid development of production methods and strengthening of bacteriophage activities. Bacteriophage adsorption to host cells initiates the process of infection. The rotating magnetic field (RMF) is a promising biotechnological method for process intensification, especially for the intensification of micromixing and mass transfer. This study evaluates the use of RMF to enhance the infection process by influencing bacteriophage adsorption rate. The RMF exposition decreased the t₅₀ and t₇₅ of bacteriophages T4 on Escherichia coli cells and vb_SauM_A phages on Staphylococcus aureus cells. The T4 phage adsorption rate increased from 3.13 × 10⁻⁹ mL × min⁻¹ to 1.64 × 10⁻⁸ mL × min⁻¹. The adsorption rate of vb_SauM_A phages exposed to RMF increased from 4.94 × 10⁻⁹ mL × min⁻¹ to 7.34 × 10⁻⁹ mL × min⁻¹. Additionally, the phage T4 zeta potential changed under RMF from −11.1 ± 0.49 mV to −7.66 ± 0.29 for unexposed and RMF-exposed bacteriophages, respectively.     

Mutation Rates for Enzyme and Morphological Loci in Barley ( HORDEUM VULGARE L.)

Spontaneous mutation rates were estimated by assaying 84,126 seedlings of a highly homozygous barley line (isogenic line 2025) for five enzyme loci. No mutants were observed in 841,260 allele replications. This result excludes, at probability level 0.95, a spontaneous mutation rate larger than 3.56 x 10^-6/locus/gamete/generation for these enzyme loci. Isogenic line 2025 also was scored for mutants at four loci governing morphological variants. No mutants were observed in 3,386,850 allele replications which indicates that the upper bound for the mutation rate for these loci is 8.85 x 10^-7. It was concluded that, even though spontaneous mutation has been important in creating variability in the barley species at the loci scored, the rate is too low to have much affect on the short-term dynamics of barley populations.     

High-Frequency Phage-Mediated Gene Transfer among Escherichia coli Cells, Determined at the Single-Cell Level

Recent whole-genome analysis suggests that lateral gene transfer by bacteriophages has contributed significantly to the genetic diversity of bacteria. To accurately determine the frequency of phage-mediated gene transfer, we employed cycling primed in situ amplification-fluorescent in situ hybridization (CPRINS-FISH) and investigated the movement of the ampicillin resistance gene among Escherichia coli cells mediated by phage at the single-cell level. Phages P1 and T4 and the newly isolated E. coli phage EC10 were used as vectors. The transduction frequencies determined by conventional plating were 3 × 10⁻⁸ to 2 × 10⁻⁶, 1 × 10⁻⁸ to 4 × 10⁻⁸, and <4 × 10⁻⁹ to 4 × 10⁻⁸ per PFU for phages P1, T4, and EC10, respectively. The frequencies of DNA transfer determined by CPRINS-FISH were 7 × 10⁻⁴ to 1 × 10⁻³, 9 × 10⁻⁴ to 3 × 10⁻³, and 5 × 10⁻⁴ to 4 × 10⁻³ for phages P1, T4, and EC10, respectively. Direct viable counting combined with CPRINS-FISH revealed that more than 20% of the cells carrying the transferred gene retained their viabilities. These results revealed that the difference in the number of viable cells carrying the transferred gene and the number of cells capable of growth on the selective medium was 3 to 4 orders of magnitude, indicating that phage-mediated exchange of DNA sequences among bacteria occurs with unexpectedly high frequency.     

Rates of Microbial Metabolism in Deep Coastal Plain Aquifers

Rates of microbial metabolism in deep anaerobic aquifers of the Atlantic coastal plain of South Carolina were investigated by both microbiological and geochemical techniques. Rates of [2-¹⁴C]acetate and [U-¹⁴C]glucose oxidation as well as geochemical evidence indicated that metabolic rates were faster in the sandy sediments composing the aquifers than in the clayey sediments of the confining layers. In the sandy aquifer sediments, estimates of the rates of CO₂ production (millimoles of CO₂ per liter per year) based on the oxidation of [2-¹⁴C] acetate were 9.4 × 10⁻³ to 2.4 × 10⁻¹ for the Black Creek aquifer, 1.1 × 10⁻² for the Middendorf aquifer, and <7 × 10⁻⁵ for the Cape Fear aquifer. These estimates were at least 2 orders of magnitude lower than previously published estimates that were based on the accumulation of CO₂ in laboratory incubations of similar deep subsurface sediments. In contrast, geochemical modeling of groundwater chemistry changes along aquifer flowpaths gave rate estimates that ranged from 10⁻⁴ to 10⁻⁶ mmol of CO₂ per liter per year. The age of these sediments (ca. 80 million years) and their organic carbon content suggest that average rates of CO₂ production could have been no more than 10⁻⁴ mmol per liter per year. Thus, laboratory incubations may greatly overestimate the in situ rates of microbial metabolism in deep subsurface environments. This has important implications for the use of laboratory incubations in attempts to estimate biorestoration capacities of deep aquifers. The rate estimates from geochemical modeling indicate that deep aquifers are among the most oligotrophic aquatic environments in which there is ongoing microbial metabolism.     

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas

De novo mutations are central for evolution, since they provide the raw material for natural selection by regenerating genetic variation. However, studying de novo mutations is challenging and is generally restricted to model species, so we have a limited understanding of the evolution of the mutation rate and spectrum between closely related species. Here, we present a mutation accumulation (MA) experiment to study de novo mutation in the unicellular green alga Chlamydomonas incerta and perform comparative analyses with its closest known relative, Chlamydomonas reinhardtii. Using whole-genome sequencing data, we estimate that the median single nucleotide mutation (SNM) rate in C. incerta is μ = 7.6 × 10⁻¹⁰, and is highly variable between MA lines, ranging from μ = 0.35 × 10⁻¹⁰ to μ = 131.7 × 10⁻¹⁰. The SNM rate is strongly positively correlated with the mutation rate for insertions and deletions between lines (r > 0.97). We infer that the genomic factors associated with variation in the mutation rate are similar to those in C. reinhardtii, allowing for cross-prediction between species. Among these genomic factors, sequence context and complexity are more important than GC content. With the exception of a remarkably high C→T bias, the SNM spectrum differs markedly between the two Chlamydomonas species. Our results suggest that similar genomic and biological characteristics may result in a similar mutation rate in the two species, whereas the SNM spectrum has more freedom to diverge.