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

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.     

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.     

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.     

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.     

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.     

Calibrating the Human Mutation Rate via Ancestral Recombination Density in Diploid Genomes

The human mutation rate is an essential parameter for studying the evolution of our species, interpreting present-day genetic variation, and understanding the incidence of genetic disease. Nevertheless, our current estimates of the rate are uncertain. Most notably, recent approaches based on counting de novo mutations in family pedigrees have yielded significantly smaller values than classical methods based on sequence divergence. Here, we propose a new method that uses the fine-scale human recombination map to calibrate the rate of accumulation of mutations. By comparing local heterozygosity levels in diploid genomes to the genetic distance scale over which these levels change, we are able to estimate a long-term mutation rate averaged over hundreds or thousands of generations. We infer a rate of 1.61 ± 0.13 × 10⁻⁸ mutations per base per generation, which falls in between phylogenetic and pedigree-based estimates, and we suggest possible mechanisms to reconcile our estimate with previous studies. Our results support intermediate-age divergences among human populations and between humans and other great apes.     

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.     

Rates of transposition in Escherichia coli

The evolutionary role of transposable elements (TEs) is still highly controversial. Two key parameters, the transposition rate (u and w, for replicative and non-replicative transposition) and the excision rate (e) are fundamental to understanding their evolution and maintenance in populations. We have estimated u, w and e for six families of TEs (including eight members: IS1, IS2, IS3, IS4, IS5, IS30, IS150 and IS186) in Escherichia coli, using a mutation accumulation (MA) experiment. In this experiment, mutations accumulate essentially at the rate at which they appear, during a period of 80 500 (1610 generations × 50 lines) generations, and spontaneous transposition events can be detected. This differs from other experiments in which insertions accumulated under strong selective pressure or over a limited genomic target. We therefore provide new estimates for the spontaneous rates of transposition and excision in E. coli. We observed 25 transposition and three excision events in 50 MA lines, leading to overall rate estimates of u ∼ 1.15 × 10^–5, w ∼ 4 × 10⁻⁸ and e ∼ 1.08 × 10⁻⁶ (per element, per generation). Furthermore, extensive variation between elements was found, consistent with previous knowledge of the mechanisms and regulation of transposition for the different elements.     

Genetic variation at the FADS1-FADS2 gene locus influences delta-5 desaturase activity and LC-PUFA proportions after fish oil supplement

Delta-5 and delta-6 desaturases (D5D and D6D) are key enzymes in endogenous synthesis of long-chain PUFAs. In this sample of healthy subjects (n = 310), genotypes of single nucleotide polymorphisms (SNPs) rs174537, rs174561, and rs3834458 in the FADS1-FADS2 gene cluster were strongly associated with proportions of LC-PUFAs and desaturase activities estimated in plasma and erythrocytes. Minor allele carriage associated with decreased activities of D5D (FADS1) (5.84 × 10⁻¹⁹ ≤ P ≤ 4.5 × 10⁻¹⁸) and D6D (FADS2) (6.05 × 10⁻⁸ ≤ P ≤ 4.20 × 10⁻⁷) was accompanied by increased substrate and decreased product proportions (0.05 ≤ P ≤ 2.49 × 10⁻¹⁶). The significance of haplotype association with D5D activity (P = 2.19 × 10⁻¹⁷) was comparable to that of single SNPs, but haplotype association with D6D activity (P = 3.39 × 10⁻²⁸) was much stronger. In a randomized controlled dietary intervention, increasing eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3) intake significantly increased D5D (P = 4.0 × 10⁻⁹) and decreased D6D activity (P = 9.16 × 10⁻⁶) after doses of 0.45, 0.9, and 1.8 g/day for six months. Interaction of rs174537 genotype with treatment was a determinant of D5D activity estimated in plasma (P = 0.05). In conclusion, different sites at the FADS1-FADS2 locus appear to influence D5D and D6D activity, and rs174537 genotype interacts with dietary EPA+DHA to modulate D5D.     

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.     

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.     

Quantifying the Diversification of Hepatitis C Virus (HCV) during Primary Infection: Estimates of the In Vivo Mutation Rate

Hepatitis C virus (HCV) is present in the host with multiple variants generated by its error prone RNA-dependent RNA polymerase. Little is known about the initial viral diversification and the viral life cycle processes that influence diversity. We studied the diversification of HCV during acute infection in 17 plasma donors, with frequent sampling early in infection. To analyze these data, we developed a new stochastic model of the HCV life cycle. We found that the accumulation of mutations is surprisingly slow: at 30 days, the viral population on average is still 46% identical to its transmitted viral genome. Fitting the model to the sequence data, we estimate the median in vivo viral mutation rate is 2.5×10⁻⁵ mutations per nucleotide per genome replication (range 1.6–6.2×10⁻⁵), about 5-fold lower than previous estimates. To confirm these results we analyzed the frequency of stop codons (N = 10) among all possible non-sense mutation targets (M = 898,335), and found a mutation rate of 2.8–3.2×10⁻⁵, consistent with the estimate from the dynamical model. The slow accumulation of mutations is consistent with slow turnover of infected cells and replication complexes within infected cells. This slow turnover is also inferred from the viral load kinetics. Our estimated mutation rate, which is similar to that of other RNA viruses (e.g., HIV and influenza), is also compatible with the accumulation of substitutions seen in HCV at the population level. Our model identifies the relevant processes (long-lived cells and slow turnover of replication complexes) and parameters involved in determining the rate of HCV diversification.     

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.     

Characteristics of the Inhibition of Potato (Solanum tuberosum) Invertase by an Endogenous Proteinaceous Inhibitor in Potatoes

Effect of several parameters on inhibition of potato (Solanum tuberosum) invertase by its endogenous proteinaceous inhibitor was determined using homogeneous preparations of both proteins. The inhibitor and invertase formed an inactive complex with an observed association rate constant at pH 4.70 and 37°C of 8.82 × 10² per molar per second and a dissociation rate constant of 3.3 × 10⁻³ per minute. The inhibitor appeared to bind to invertase in more than one step. Initial interaction (measured by loss of invertase activity) was rapid, relatively weak, readily reversible (K_i of 2 × 10⁻⁶ molar) and noncompetitive with substrate at pH 4.70. Initial interaction was probably followed by isomerization to a tighter (K_i of 6.23 × 10⁻⁸ molar) complex, which dissociated slowly with a half-time of 3.5 hour. Interaction between enzyme and inhibitor appeared to be of ionic character and essentially pH independent between pH 3.5 and 7.4.     

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.     

Phenotypic mutation rates and the abundance of abnormal proteins in yeast

Phenotypic mutations are errors that occur during protein synthesis. These errors lead to amino acid substitutions that give rise to abnormal proteins. Experiments suggest that such errors are quite common. We present a model to study the effect of phenotypic mutation rates on the amount of abnormal proteins in a cell. In our model, genes are regulated to synthesize a certain number of functional proteins. During this process, depending on the phenotypic mutation rate, abnormal proteins are generated. We use data on protein length and abundance in Saccharomyces cerevisiae to parametrize our model. We calculate that for small phenotypic mutation rates most abnormal proteins originate from highly expressed genes that are on average nearly twice as large as the average yeast protein. For phenotypic mutation rates much above 5 × 10⁻⁴, the error-free synthesis of large proteins is nearly impossible and lowly expressed, very large proteins contribute more and more to the amount of abnormal proteins in a cell. This fact leads to a steep increase of the amount of abnormal proteins for phenotypic mutation rates above 5 × 10⁻⁴. Simulations show that this property leads to an upper limit for the phenotypic mutation rate of approximately 2 × 10⁻³ even if the costs for abnormal proteins are extremely low. We also consider the adaptation of individual proteins. Individual genes/proteins can decrease their phenotypic mutation rate by using preferred codons or by increasing their robustness against amino acid substitutions. We discuss the similarities and differences between the two mechanisms and show that they can only slow down but not prevent the rapid increase of the amount of abnormal proteins. Our work allows us to estimate the phenotypic mutation rate based on data on the fraction of abnormal proteins. For S. cerevisiae, we predict that the value for the phenotypic mutation rate is between 2 × 10⁻⁴ and 6 × 10⁻⁴.     

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.     

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.     

Pedigree-based and phylogenetic methods support surprising patterns of mutation rate and spectrum in the gray mouse lemur

Mutations are the raw material on which evolution acts, and knowledge of their frequency and genomic distribution is crucial for understanding how evolution operates at both long and short timescales. At present, the rate and spectrum of de novo mutations have been directly characterized in relatively few lineages. Our study provides the first direct mutation-rate estimate for a strepsirrhine (i.e., the lemurs and lorises), which comprises nearly half of the primate clade. Using high-coverage linked-read sequencing for a focal quartet of gray mouse lemurs (Microcebus murinus), we estimated the mutation rate to be among the highest calculated for a mammal at 1.52 × 10^–8 (95% credible interval: 1.28 × 10⁻⁸–1.78 × 10⁻⁸) mutations/site/generation. Further, we found an unexpectedly low count of paternal mutations, and only a modest overrepresentation of mutations at CpG sites. Despite the surprising nature of these results, we found both the rate and spectrum to be robust to the manipulation of a wide range of computational filtering criteria. We also sequenced a technical replicate to estimate a false-negative and false-positive rate for our data and show that any point estimate of a de novo mutation rate should be considered with a large degree of uncertainty. For validation, we conducted an independent analysis of context-dependent substitution types for gray mouse lemur and five additional primate species for which de novo mutation rates have also been estimated. These comparisons revealed general consistency of the mutation spectrum between the pedigree-based and the substitution-rate analyses for all species compared.Subject terms: Evolution, Molecular evolution