A new era of mutation rate analyses: Concepts and methods.

The mutation rate is a pivotal biological characteristic, intricately governed by natural selection and historically garnering considerable attention. Recent advances in high-throughput sequencing and analytical methodologies have profoundly transformed our understanding in this domain, ushering in an unprecedented era of mutation rate research. This paper aims to provide a comprehensive overview of the key concepts and methodologies frequently employed in the study of mutation rates. It examines various types of mutations, explores the evolutionary dynamics and associated theories, and synthesizes both classical and contemporary hypotheses. Furthermore, this review comprehensively explores recent advances in understanding germline and somatic mutations in animals and offers an overview of experimental methodologies, mutational patterns, molecular mechanisms, and driving forces influencing variations in mutation rates across species and tissues. Finally, it proposes several potential research directions and pressing questions for future investigations.

Coevolution of longevity and female germline maintenance.

An often-overlooked aspect of life-history optimization is the allocation of resources to protect the germline and secure safe transmission of genetic information. While failure to do so renders significant fitness consequences in future generations, germline maintenance comes with substantial costs. Thus, germline allocation should trade off with other life-history decisions and be optimized in accordance with an organism's reproductive schedule. Here, we tested this hypothesis by studying germline maintenance in lines of seed beetle, selected for early (E) or late (L) reproduction for 350 and 240 generations, respectively. Female animals provide maintenance and screening of male gametes in their reproductive tract and oocytes. Here, we reveal the ability of young and aged E- and L-females to provide this form of germline maintenance by mating them to males with ejaculates with artificially elevated levels of protein and DNA damage. We find that germline maintenance in E-females peaks at young age and then declines, while the opposite is true for L-females, in accordance with the age of reproduction in the respective regime. These findings identify the central role of allocation to secure germline integrity in life-history evolution and highlight how females can play a crucial role in mitigating the effects of male germline decisions on mutation rate and offspring quality.

The rate and nature of mitochondrial DNA mutations in human pedigrees.

We examined the rate and nature of mitochondrial DNA (mtDNA) mutations in humans using sequence data from 64,806 contemporary Icelanders from 2,548 matrilines. Based on 116,663 mother-child transmissions, 8,199 mutations were detected, providing robust rate estimates by nucleotide type, functional impact, position, and different alleles at the same position. We thoroughly document the true extent of hypermutability in mtDNA, mainly affecting the control region but also some coding-region variants. The results reveal the impact of negative selection on viable deleterious mutations, including rapidly mutating disease-associated 3243A>G and 1555A>G and pre-natal selection that most likely occurs during the development of oocytes. Finally, we show that the fate of new mutations is determined by a drastic germline bottleneck, amounting to an average of 3 mtDNA units effectively transmitted from mother to child.

Plasma ctDNA enhances the tissue-based detection of oncodriver mutations in colorectal cancer.

PURPOSE: The advent of circulating tumor DNA (ctDNA) technology has provided a convenient and noninvasive means to continuously monitor cancer genomic data, facilitating personalized cancer treatment. This study aimed to evaluate the supplementary benefits of plasma ctDNA alongside traditional tissue-based next-generation sequencing (NGS) in identifying targetable mutations and tumor mutational burden (TMB) in colorectal cancers (CRC). METHODS: Our study involved 76 CRC patients, collecting both tissue and plasma samples for NGS. We assessed the concordance of gene mutational status between ctDNA and tissue, focusing on actionable genes such as KRAS, NRAS, PIK3CA, BRAF, and ERBB2. Logistic regression analysis was used to explore variables associated with discordance and positive mutation rates. RESULTS: In total, 26 cancer-related genes were identified. The most common variants in tumor tissues and plasma samples were in APC (57.9% vs 19.7%), TP53 (55.3% vs 22.4%) and KRAS (47.4% vs 43.4%). Tissue and ctDNA showed an overall concordance of 73.53% in detecting actionable gene mutations. Notably, plasma ctDNA improved detection for certain genes and gene pools. Variables significantly associated with discordance included gender and peritoneal metastases. TMB analysis revealed a higher detection rate in tissues compared to plasma, but combining both increased detection. CONCLUSIONS: Our study highlights the importance of analyzing both tissue and plasma for detecting actionable mutations in CRC, with plasma ctDNA offering added value. Discordance is associated with gender and peritoneal metastases, and TMB analysis can benefit from a combination of tissue and plasma data. This approach provides valuable insights for personalized CRC treatment.

The Tradeoffs Between Persistence and Mutation Rates at Sub-Inhibitory Antibiotic Concentrations in Staphylococcus aureus.

The rational design of the antibiotic treatment of bacterial infections employs these drugs to reach concentrations that exceed the minimum needed to prevent the replication of the target bacteria. However, within a treated patient, spatial and physiological heterogeneity promotes antibiotic gradients such that the concentration of antibiotics at specific sites is below the minimum needed to inhibit bacterial growth. Here, we investigate the effects of sub-inhibitory antibiotic concentrations on three parameters central to bacterial infection and the success of antibiotic treatment, using in vitro experiments with Staphylococcus aureus and mathematical-computer simulation models. Our results, using drugs of six different classes, demonstrate that exposure to sub-inhibitory antibiotic concentrations not only alters the dynamics of bacterial growth but also increases the mutation rate to antibiotic resistance and decreases the rate of production of persister cells thereby reducing the persistence level. Understanding this trade-off between mutation rates and persistence levels resulting from sub-inhibitory antibiotic exposure is crucial for optimizing, and mitigating the failure of, antibiotic therapy.

Environmental and genetic influence on the rate and spectrum of spontaneous mutations in Escherichia coli.

Spontaneous mutations are the ultimate source of novel genetic variation on which evolution operates. Although mutation rate is often discussed as a single parameter in evolution, it comprises multiple distinct types of changes at the level of DNA. Moreover, the rates of these distinct changes can be independently influenced by genomic background and environmental conditions. Using fluctuation tests, we characterized the spectrum of spontaneous mutations in Escherichia coli grown in low and high glucose environments. These conditions are known to affect the rate of spontaneous mutation in wild-type MG1655, but not in a ΔluxS deletant strain - a gene with roles in both quorum sensing and the recycling of methylation products used in E. coli's DNA repair process. We find an increase in AT>GC transitions in the low glucose environment, suggesting that processes relating to the production or repair of this mutation could drive the response of overall mutation rate to glucose concentration. Interestingly, this increase in AT>GC transitions is maintained by the glucose non-responsive ΔluxS deletant. Instead, an elevated rate of GC>TA transversions, more common in a high glucose environment, leads to a net non-responsiveness of overall mutation rate for this strain. Our results show how relatively subtle changes, such as the concentration of a carbon substrate or loss of a regulatory gene, can substantially influence the amount and nature of genetic variation available to selection.

Mutational analysis of 23 autosomal short tandem repeats based on trio paternity testing in the Korean population.

This study aimed to estimate A-STR mutation rates in 2,317 Korean parent-child trios by examining 20 Combined DNA Index System (CODIS) core loci (D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S51, D21S11, CSF1PO, FGA, TH01, TPOX, vWA, D1S1656, D2S441, D2S1338, D10S1248, D12S391, D19S433, and D22S1045) and three non-CODIS loci (Penta E, Penta D, and SE33). Locus-specific mutation rate estimates varied from 0.00 to 8.63 × 10-3 per generation, with an average mutation rate of 1.62 × 10-3 (95 % CI, 1.39-1.88 × 10-3). We also combined data from previous studies to obtain comprehensive genetic values for the Korean population, and the average mutation rate was 1.59 × 10-3 (95 % CI, 1.38-1.82 × 10-3). Single-step mutations (95.69 %) and double-step mutations (3.35 %) were observed in the mutation pattern analysis, and cases expected to have multi-step mutations (0.96 %) were also observed. Large-sized alleles exhibited more loss mutations than gain mutations, and paternal mutations (62.68 %) were more frequently observed than maternal mutations (19.62 %). The calculated values and features of the 23 A-STRs explored in this study are expected to play a crucial role in establishing criteria for forensic genetic interpretation.

What are we missing in teaching the Luria-Delbrück experiment?

The importance of teaching the Luria-Delbrück experiment to biology students is increasingly recognized by educators, and improved pedagogical methods for teaching the classic experiment have been proposed and tested in the classroom. However, there are still obstacles that impede the proper teaching of the classic experiment. This note proposes two strategies to further improve the teaching of the classic experiment. The first strategy is to be frank with an inherent limitation of the classic experiment, and instructors should explain from a logical point of view why the classic experiment cannot be used to refute the possibility of directed mutation. The second strategy is to emphasize the pioneering work of Delbrück on developing the mutant distribution that enables researchers to estimate microbial mutation rates using data generated by fluctuation experiments, and instructors should shift their attention to the overlooked essential role of the mutant distribution.

Clonal expansion in normal tissues.

Cancer originates from a single ancestral cell that acquires a driver mutation, which confers a growth or survival advantage, followed by the acquisition of additional driver mutations by descendant cells. Recently, it has become evident that somatic cell mutations accumulate in normal tissues with aging and exposure to environmental factors, such as alcohol, smoking, and UV rays, increases the mutation rate. Clones harboring driver mutations expand with age, leading to tissue remodeling. Lineage analysis of myeloproliferative neoplasms and der(1;16)-positive breast cancer revealed that driver mutations were acquired early in our lives and that the development of cancer takes decades, unveiling the previously unknown early process of cancer development. Evidence that clonal hematopoiesis affects various diseases, including nonneoplastic diseases, highlights the potential role of the identification and functional analysis of mutated clones in unraveling unknown pathologies. In this review, we summarize the recent updates on clonal expansion in normal tissues and the natural history of cancer revealed through lineage analysis of noncancerous and cancerous tissues.

Low mutation rate of spontaneous mutants enables detection of causative genes by comparing whole genome sequences.

In the early 1900s, mutation breeding to select varieties with desirable traits using spontaneous mutation was actively conducted around the world, including Japan. In rice, the number of fixed mutations per generation was estimated to be 1.38-2.25. Although this low mutation rate was a major problem for breeding in those days, in the modern era with the development of next-generation sequencing (NGS) technology, it was conversely considered to be an advantage for efficient gene identification. In this paper, we proposed an in silico approach using NGS to compare the whole genome sequence of a spontaneous mutant with that of a closely related strain with a nearly identical genome, to find polymorphisms that differ between them, and to identify the causal gene by predicting the functional variation of the gene caused by the polymorphism. Using this approach, we found four causal genes for the dwarf mutation, the round shape grain mutation and the awnless mutation. Three of these genes were the same as those previously reported, but one was a novel gene involved in awn formation. The novel gene was isolated from Bozu-Aikoku, a mutant of Aikoku with the awnless trait, in which nine polymorphisms were predicted to alter gene function by their whole-genome comparison. Based on the information on gene function and tissue-specific expression patterns of these candidate genes, Os03g0115700/LOC_Os03g02460, annotated as a short-chain dehydrogenase/reductase SDR family protein, is most likely to be involved in the awnless mutation. Indeed, complementation tests by transformation showed that it is involved in awn formation. Thus, this method is an effective way to accelerate genome breeding of various crop species by enabling the identification of useful genes that can be used for crop breeding with minimal effort for NGS analysis.

Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

Considerations on the application of a mutation model for Y-STR interpretation.

If Y-STR profiling is to be more effective in criminal casework, the methods used to evaluate evidential weight require improvement. Many forensic scientists assign an evidential weight by estimating the number of times a Y-STR profile obtained from a questioned sample has been observed in YHRD datasets. More sophisticated models have been suggested but not yet implemented into routine casework, e.g. Andersen & Balding [1]. Mutation is inherent to STR meiosis (or inheritance) and is encountered in practice. We evaluated a mutation model that can be incorporated into a method for assigning evidential weight to Y-STR profiles, an essential part of bringing any method into practice. Since an important part of implementation to casework is communication, the article is written in an accessible format for practitioners as well as statisticians. The mutation component within the MUTEA model by Willems et al. [2] incorporates the potential for multistep mutations and a tendency for alleles to revert towards a central length, reflecting observed mutation data, e.g. [3]. We have estimated the parameters in this model and in a simplified symmetric version of this model, using sequence data from father/son pairs [4] and deep-rooted pedigrees [5]. Both datasets contain multistep mutations, which may have an effect on models based on simulations [1]. We introduce Beta-Binomial and Beta-Geometric conjugate analyses for estimating rate and step parameters for the mutation models presented here, which require only summations and multiplications. We proved mathematically that the parameters can be estimated independently. We show the importance of reporting the variability of the parameters and not only a point estimate. The parameters can be easily incorporated into statistical models, and updated sequentially as more data becomes available. We recommend fuller publication of data to enable the development and evaluation of a wider range of mutation models.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus.

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Manufacturing DNA in E. coli yields higher-fidelity DNA than in vitro enzymatic synthesis.

Biotechnologies such as gene therapy have brought DNA vectors to the forefront of pharmaceuticals. The quality of starting material plays a pivotal role in determining final product quality. Here, we examined the fidelity of DNA replication using enzymatic methods (in vitro) compared to plasmid DNA produced in vivo in E. coli. Next-generation sequencing approaches rely on in vitro polymerases, which have inherent limitations in sensitivity. To address this challenge, we introduce a novel assay based on loss-of-function (LOF) mutations in the conditionally toxic sacB gene. Our findings show that DNA production in E. coli results in significantly fewer LOF mutations (80- to 3,000-fold less) compared to enzymatic DNA replication methods such as polymerase chain reaction (PCR) and rolling circle amplification (RCA). These results suggest that using DNA produced by PCR or RCA may introduce a substantial number of mutation impurities, potentially affecting the quality and yield of final pharmaceutical products. Our study underscores that DNA synthesized in vitro has a significantly higher mutation rate than DNA produced traditionally in E. coli. Therefore, utilizing in vitro enzymatically produced DNA in biotechnology and biomanufacturing may entail considerable fidelity-related risks, while using DNA starting material derived from E. coli substantially mitigates this risk.

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history.

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

Chemical manipulation of m1A mediates its detection in human tRNA.

N 1-methyl adenosine (m1A) is a widespread RNA modification present in tRNA, rRNA, and mRNA. m1A modification sites in tRNAs are evolutionarily conserved and its formation on tRNA is catalyzed by methyltransferase TRMT61A and TRMT6 complex. m1A promotes translation initiation and elongation. Due to its positive charge under physiological conditions, m1A can notably modulate RNA structure. It also blocks Watson-Crick-Franklin base-pairing and causes mutation and truncation during reverse transcription. Several misincorporation-based high-throughput sequencing methods have been developed to sequence m1A. In this study, we introduce a reduction-based m1A sequencing (red-m1A-seq). We report that NaBH4 reduction of m1A can improve the mutation and readthrough rates using commercially available RT enzymes to give a better positive signature, while alkaline-catalyzed Dimroth rearrangement can efficiently convert m1A to m6A to provide good controls, allowing the detection of m1A with higher sensitivity and accuracy. We applied red-m1A-seq to sequence human small RNA, and we not only detected all the previously reported tRNA m1A sites, but also new m1A sites in mt-tRNAAsn-GTT and 5.8S rRNA.

Temporal challenges in detecting balancing selection from population genomic data.

The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice.

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Emergence of a multilocus mutator genotype in mutator Escherichia coli experimental populations under repeated lethal selection.

Mutator alleles, which confer increased mutation rates, are known to spontaneously emerge and "hitchhike" to fixation in evolving asexual populations. Theory predicts that in an evolving asexual mutator population, a second mutator allele may spontaneously arise and hitchhike to fixation. Here, we describe an empirical test of the hypothesis of repeated hitchhiking. The starting population was a clonal strain of mutL-Escherichia coli whose mutation rate was 100-fold higher than wild type. We exposed the mutL- strain to a series of three antibiotics in increasing order of selective strength: fosfomycin, rifampicin, and streptomycin. Two independent replicates of the experiment were performed. As predicted, elevated mutation rates and enrichment for multilocus mutators (which bear more than one mutator allele) were observed in the end point populations of both experiments. DNA sequencing revealed an identical spontaneous 1-bp insertion in the mutator gene mutT in both end point populations. In the multilocus mutators, the causal relationship between the mutT- mutations and the increase in mutation rate was supported with mutT+ plasmid complementation tests. Surprisingly, when the experiment was repeated with the antibiotics deployed in decreasing order of selective strength, enrichment for multilocus mutators was not observed. Our data support the likelihood that the mutT- mutations rose to fixation in both populations, consistent with the hypothesis of repeated mutator hitchhiking. The escalation of mutation rates in asexual populations is relevant to multiple biological scenarios, including antibiotic resistance, host-pathogen interactions, and carcinogenesis.

Effects of parental age and polymer composition on short tandem repeat de novo mutation rates.

Short tandem repeats (STRs) are hotspots of genomic variability in the human germline because of their high mutation rates, which have long been attributed largely to polymerase slippage during DNA replication. This model suggests that STR mutation rates should scale linearly with a father's age, as progenitor cells continually divide after puberty. In contrast, it suggests that STR mutation rates should not scale with a mother's age at her child's conception, since oocytes spend a mother's reproductive years arrested in meiosis II and undergo a fixed number of cell divisions that are independent of the age at ovulation. Yet, mirroring recent findings, we find that STR mutation rates covary with paternal and maternal age, implying that some STR mutations are caused by DNA damage in quiescent cells rather than polymerase slippage in replicating progenitor cells. These results echo the recent finding that DNA damage in oocytes is a significant source of de novo single nucleotide variants and corroborate evidence of STR expansion in postmitotic cells. However, we find that the maternal age effect is not confined to known hotspots of oocyte mutagenesis, nor are postzygotic mutations likely to contribute significantly. STR nucleotide composition demonstrates divergent effects on de novo mutation (DNM) rates between sexes. Unlike the paternal lineage, maternally derived DNMs at A/T STRs display a significantly greater association with maternal age than DNMs at G/C-containing STRs. These observations may suggest the mechanism and developmental timing of certain STR mutations and contradict prior attribution of replication slippage as the primary mechanism of STR mutagenesis.