Spontaneous Mutation Rate Research Articles

The Challenging Life of Mutators: How Pseudomonas aeruginosa Survives between Persistence and Evolution in Cystic Fibrosis Lung.

Cystic fibrosis (CF) is a life-threatening genetic disease characterised by chronic lung infections sustained by opportunistic pathogens such as Pseudomonas aeruginosa. During the chronic long-lasting lung infections, P. aeruginosa adapts to the host environment. Hypermutability, mainly due to defects in the DNA repair system, resulting in an increased spontaneous mutation rate, represents a way to boost the rapid adaptation frequently encountered in CF P. aeruginosa isolates. We selected 609 isolates from 51 patients with CF chronically colonised by P. aeruginosa to study, by full-length genome sequencing, the longitudinal evolution of the bacterium. We recovered at least one hypermutable (mutator) isolate in 57% of patients. By combining genomic information and phenotypic analyses, we followed the evolutionary pathways of the P. aeruginosa mutator strains, identifying their contribution to multi-drug resistance and the emergence of new sub-lineages. By implementing patient clinical data, we observed that mutators preferentially follow a specific evolutionary trajectory in patients with a negative clinical outcome and that maintenance antibiotic polytherapy, based on alternating molecules, apparently reduces the occurrence of hypermutability. Finally, we draw attention to the possibility that modulator-induced changes in the pulmonary environment may be associated with the onset of hypermutability.

Assessment of the In Vitro Phosphatidylinositol Glycan Class A (PIG-A) Gene Mutation Assay Using Human TK6 and Mouse Hepa1c1c7 Cell Lines.

Gene mutations linked to diseases like cancer may be caused by exposure to environmental chemicals. The X-linked phosphatidylinositol glycan class A (PIG-A) gene, required for glycosylphosphatidylinositol (GPI) anchor biosynthesis, is a key target locus for in vitro genetic toxicity assays. Various organisms and cell lines may respond differently to genotoxic agents. Here, we compared the mutagenic potential of directly genotoxic ethyl methane sulfonate (EMS) to metabolically activated pro-mutagenic polycyclic aromatic hydrocarbons (PAHs). The two classes of mutagens were compared in an in vitro PIG-A gene mutation test using the metabolically active murine hepatoma Hepa1c1c7 cell line and the human TK6 cell line, which has limited metabolic capability. Determination of cell viability is required for quantifying mutagenicity. Two common cell viability tests, the MTT assay and propidium iodide (PI) staining measured by flow cytometry, were evaluated. The MTT assay overestimated cell viability in adherent cells at high benzo[a]pyrene (B[a]P) exposure concentrations, so PI-based cytotoxicity was used in calculations. The spontaneous mutation rates for TK6 and Hepa1c1c7 cells were 1.87 and 1.57 per million cells per cell cycle, respectively. TK6 cells exposed to 600 µM and 800 µM EMS showed significantly higher mutation frequencies (36 and 47 per million cells per cell cycle, respectively). Exposure to the pro-mutagen benzo[a]pyrene (B[a]P, 10 µM) did not increase mutation frequency in TK6 cells. In Hepa1c1c7 cells, mutation frequencies varied across exposure groups (50, 50, 29, and 81 per million cells per cell cycle when exposed to 10 µM B[a]P, 5-methylcholanthrene (5-MC), chrysene, or 16,000 µM EMS, respectively). We demonstrate that the choice of cytotoxicity assay and cell line can determine the outcome of the Pig-A mutagenesis assay when assessing a specific mutagen.

Meiosis I causes a high spontaneous mutation rate in a multicellular red alga (Pyropia yezoensis) with a complex life cycle

Mutations are the origin of genetic diversity and are fundamental parameters needed to understand the molecular evolution of species. Estimations of mutation rates have been conducted for many diverse taxa, although rates in several major eukaryotic lineages remain unexplored. Here, the first estimation is reported of the spontaneous mutation rate for the multicellular eukaryote red alga, Pyropia yezoensis, which exhibits a complex life cycle. An estimated mutation rate of 2.97 × 10−8 (95 % CI: 2.16 × 10−8–3.99 × 10−8) per site per generation was generated for the primary life cycle, the sexual cycle, which is the highest sexual mutation rate among published sexual plants. Combined with tetrad analysis, meiosis I was identified as the primary period responsible for the high mutation rate during the complex life cycle of P. yezoensis. This result provides direct evidence for the “meiosis is mutagenic” hypothesis for multicellular organisms. The accurate estimate of the mutation rate of P. yezoensis also informs several immediate applications. Based on the above estimate, the effective population size (Ne) of P. yezoensis was estimated at about 19,000, with extensive haploid phases and asexual reproduction through monospores possibly leading to linked selection that may reduce the genome-wide genetic diversity of P. yezoensis and consequently influence Ne estimation. Lastly, P. yezoensis was estimated to have diverged from P. haitanensis about 4.2 Ma, representing a more recent date than estimates from fossil-calibrated phylogenies. These findings provide valuable new information for understanding the evolution of red algae, in addition to the underlying mechanism of mutations.

Insights into macrolide resistance in Arcobacter butzleri: potential resistance mechanisms and impact on bacterial fitness and virulence.

Macrolides are recommended for treating the emerging enteropathogen Arcobacter butzleri; nonetheless, this bacterium often exhibits highly variable resistance rates, and the mechanisms behind this resistance phenotype remain largely unexplored. To understand the phenotypic and genotypic consequences associated with the acquisition of erythromycin resistance in A. butzleri, as well as the effects on the fitness of this species. Resistant strains resulting from spontaneous mutations and adaptive laboratory evolution under increasing erythromycin concentrations were examined regarding their cross-resistance and collateral susceptibility profiles. Genetic causes of phenotypic antibiotic resistance were analysed by sequencing and bioinformatics, with functional correlation through ethidium bromide accumulation assays. Growth profiles in the presence and absence of erythromycin, motility and biofilm formation abilities were assessed to detect potential changes in fitness and virulence. Clones from spontaneous mutation rate evolution demonstrated decreased susceptibility to erythromycin and other classes of antibiotics, associated with mutations in the transcriptional repressor areR, causing overexpression of the AreABC efflux pump. In turn, WGS analysis of the evolved strain showed additional mutations in the ribosomal proteins L4 and L22 and in the areR gene. Furthermore, the acquisition of macrolide resistance altered A. butzleri virulence and entailed a high biological cost. The findings of this study have proved that efflux activity contributes synergistically with mutations in the ribosomal proteins L4 and L22 to A. butzleri's high-level macrolide resistance. The results further suggest an impact on the bacterial physiology and virulence, with the increased fitness cost justifying the low worldwide prevalence of high-level resistant circulating strains.

Dual genetic level modification engineering accelerate genome evolution of Corynebacterium glutamicum.

High spontaneous mutation rate is crucial for obtaining ideal phenotype and exploring the relationship between genes and phenotype. How to break the genetic stability of organisms and increase the mutation frequency has become a research hotspot. Here, we present a practical and controllable evolutionary tool (oMut-Cgts) based on dual genetic level modification engineering for Corynebacterium glutamicum. Firstly, the modification engineering of transcription and replication levels based on RNA polymerase α subunit and DNA helicase Cgl0854 as the 'dock' of cytidine deaminase (pmCDA1) significantly increased the mutation rate, proving that the localization of pmCDA1 around transient ssDNA is necessary for genome mutation. Then, the combined modification and optimization of engineering at dual genetic level achieved 1.02×104-fold increased mutation rate. The genome sequencing revealed that the oMut-Cgts perform uniform and efficient C:G→T:A transitions on a genome-wide scale. Furthermore, oMut-Cgts-mediated rapid evolution of C. glutamicum with stress (acid, oxidative and ethanol) tolerance proved that the tool has powerful functions in multi-dimensional biological engineering (rapid phenotype evolution, gene function mining and protein evolution). The strategies for rapid genome evolution provided in this study are expected to be applicable to a variety of applications in all prokaryotic cells.

Identification of the mutual gliding locus as a factor for gut colonization in non-native bee hosts using the ARTP mutagenesis

BackgroundThe gut microbiota and their hosts profoundly affect each other’s physiology and evolution. Identifying host-selected traits is crucial to understanding the processes that govern the evolving interactions between animals and symbiotic microbes. Current experimental approaches mainly focus on the model bacteria, like hypermutating Escherichia coli or the evolutionary changes of wild stains by host transmissions. A method called atmospheric and room temperature plasma (ARTP) may overcome the bottleneck of low spontaneous mutation rates while maintaining mild conditions for the gut bacteria.ResultsWe established an experimental symbiotic system with gnotobiotic bee models to unravel the molecular mechanisms promoting host colonization. By in vivo serial passage, we tracked the genetic changes of ARTP-treated Snodgrassella strains from Bombus terrestris in the non-native honeybee host. We observed that passaged isolates showing genetic changes in the mutual gliding locus have a competitive advantage in the non-native host. Specifically, alleles in the orphan mglB, the GTPase activating protein, promoted colonization potentially by altering the type IV pili-dependent motility of the cells. Finally, competition assays confirmed that the mutations out-competed the ancestral strain in the non-native honeybee gut but not in the native host.ConclusionsUsing the ARTP mutagenesis to generate a mutation library of gut symbionts, we explored the potential genetic mechanisms for improved gut colonization in non-native hosts. Our findings demonstrate the implication of the cell mutual-gliding motility in host association and provide an experimental system for future study on host-microbe interactions.ArN8Xjm3N6Dubzus4DCSQmVideo

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

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.

Human-mouse comparison of the multistage nature of radiation carcinogenesis in a mathematical model.

Mouse models are vital for assessing risk from environmental carcinogens, including ionizing radiation, yet the interspecies difference in the dose response precludes direct application of experimental evidence to humans. Herein, we take a mathematical approach to delineate the mechanism underlying the human-mouse difference in radiation-related cancer risk. We used a multistage carcinogenesis model assuming a mutational action of radiation to analyze previous data on cancer mortality in the Japanese atomic bomb survivors and in lifespan mouse experiments. Theoretically, the model predicted that exposure will chronologically shift the age-related increase in cancer risk forward by a period corresponding to the time in which the spontaneous mutational process generates the same mutational burden as that the exposure generates. This model appropriately fitted both human and mouse data and suggested a linear dose response for the time shift. The effect per dose decreased with increasing age at exposure similarly between humans and mice on a per-lifespan basis (0.72- and 0.71-fold, respectively, for every tenth lifetime). The time shift per dose was larger by two orders of magnitude in humans (7.8 and 0.046 years per Gy for humans and mice, respectively, when exposed at ~35% of their lifetime). The difference was mostly explained by the two orders of magnitude difference in spontaneous somatic mutation rates between the species plus the species-independent radiation-induced mutation rate. Thus, the findings delineate the mechanism underlying the interspecies difference in radiation-associated cancer mortality and may lead to the use of experimental evidence for risk prediction in humans.

Unveiling elevated spontaneous mutation rates in Phyllostachys edulis (Moso Bamboo) through whole genome sequencing (WGS) and investigating the impact of atmospheric and room temperature plasma (ARTP) induced mutagenesis

Moso bamboo, recognized for its wide distribution and economic importance, encounters challenges in varietal enhancement owing to its protracted sexual reproduction cycle. This study employed whole-genome resequencing to identify spontaneous mutations in Moso bamboo and investigated mutagenesis using atmospheric and room temperature plasma (ARTP). Through the sequencing results, we identified a population of flowering bamboo as an asexual breeding line. Notably, the flowering Moso bamboo population, exclusively derived from asexual reproduction, exhibited a high spontaneous mutation rate (4.54 × 10−4 to 1.15 × 10−3/bp) during sexual reproduction, considering parental and cross-pollination effects. Genetic disparities between the offspring and parents exhibited a bimodal distribution, indicating a substantial cross-pollination rate. ARTP mutagenesis increased structural variations in offspring, whereas changes in SNPs and INDELs were less pronounced. Sanger sequencing validated the gene subset, providing a foundation for the investigation of spontaneous mutation rates via whole-genome sequencing. These insights, particularly from mutagenized offspring sequencing, can contribute to Moso bamboo breeding strategies.

Spontaneous mutations and mutational responses to penicillin treatment in the bacterial pathogen Streptococcus pneumoniae D39

Bacteria with functional DNA repair systems are expected to have low mutation rates due to strong natural selection for genomic stability. However, our study of the wild-type Streptococcus pneumoniae D39, a pathogen responsible for many common diseases, revealed a high spontaneous mutation rate of 0.02 per genome per cell division in mutation-accumulation (MA) lines. This rate is orders of magnitude higher than that of other non-mutator bacteria and is characterized by a high mutation bias in the A/T direction. The high mutation rate may have resulted from a reduction in the overall efficiency of selection, conferred by the tiny effective population size in nature. In line with this, S. pneumoniae D39 also exhibited the lowest DNA mismatch-repair (MMR) efficiency among bacteria. Treatment with the antibiotic penicillin did not elevate the mutation rate, as penicillin did not induce DNA damage and S. pneumoniae lacks a stress response pathway. Our findings suggested that the MA results are applicable to within-host scenarios and provide insights into pathogen evolution.

Bacillus subtilis stress-associated mutagenesis and developmental DNA repair.

SUMMARYThe metabolic conditions that prevail during bacterial growth have evolved with the faithful operation of repair systems that recognize and eliminate DNA lesions caused by intracellular and exogenous agents. This idea is supported by the low rate of spontaneous mutations (10-9) that occur in replicating cells, maintaining genome integrity. In contrast, when growth and/or replication cease, bacteria frequently process DNA lesions in an error-prone manner. DNA repairs provide cells with the tools needed for maintaining homeostasis during stressful conditions and depend on the developmental context in which repair events occur. Thus, different physiological scenarios can be anticipated. In nutritionally stressed bacteria, different components of the base excision repair pathway may process damaged DNA in an error-prone approach, promoting genetic variability. Interestingly, suppressing the mismatch repair machinery and activating specific DNA glycosylases promote stationary-phase mutations. Current evidence also suggests that in resting cells, coupling repair processes to actively transcribed genes may promote multiple genetic transactions that are advantageous for stressed cells. DNA repair during sporulation is of interest as a model to understand how transcriptional processes influence the formation of mutations in conditions where replication is halted. Current reports indicate that transcriptional coupling repair-dependent and -independent processes operate in differentiating cells to process spontaneous and induced DNA damage and that error-prone synthesis of DNA is involved in these events. These and other noncanonical ways of DNA repair that contribute to mutagenesis, survival, and evolution are reviewed in this manuscript.

Correction of non-random mutational biases along a linear bacterial chromosome by the mismatch repair endonuclease NucS.

The linear chromosome of Streptomyces exhibits a highly compartmentalized structure with a conserved central region flanked by variable arms. As double strand break (DSB) repair mechanisms play a crucial role in shaping the genome plasticity of Streptomyces, we investigated the role of EndoMS/NucS, a recently characterized endonuclease involved in a non-canonical mismatch repair (MMR) mechanism in archaea and actinobacteria, that singularly corrects mismatches by creating a DSB. We showed that Streptomyces mutants lacking NucS display a marked colonial phenotype and a drastic increase in spontaneous mutation rate. In vitro biochemical assays revealed that NucS cooperates with the replication clamp to efficiently cleave G/T, G/G and T/T mismatched DNA by producing DSBs. These findings are consistent with the transition-shifted mutational spectrum observed in the mutant strains and reveal that NucS-dependent MMR specific task is to eliminate G/T mismatches generated by the DNA polymerase during replication. Interestingly, our data unveil a crescent-shaped distribution of the transition frequency from the replication origin towards the chromosomal ends, shedding light on a possible link between NucS-mediated DSBs and Streptomyces genome evolution.

Horizontal gene transfer facilitates the molecular reverse-evolution of antibiotic sensitivity in experimental populations of H. pylori.

Reversing the evolution of traits harmful to humans, such as antimicrobial resistance, is a key ambition of applied evolutionary biology. A major impediment to reverse evolution is the relatively low spontaneous mutation rates that revert evolved genotypes back to their ancestral state. However, the repeated re-introduction of ancestral alleles by horizontal gene transfer (HGT) could make reverse evolution likely. Here we evolve populations of an antibiotic-resistant strain of Helicobacter pylori in growth conditions without antibiotics while introducing an ancestral antibiotic-sensitive allele by HGT. We evaluate reverse evolution using DNA sequencing and find that HGT facilitates the molecular reverse evolution of the antibiotic resistance allele, and that selection for high rates of HGT drives the evolution of increased HGT rates in low-HGT treatment populations. Finally, we use a theoretical model and carry out simulations to infer how the fitness costs of antibiotic resistance, rates of HGT and effects of genetic drift interact to determine the probability and predictability of reverse evolution.

Parthenogenomics: Insights on mutation rates and nucleotide diversity in parthenogenetic Panagrolaimus nematodes.

Asexual reproduction is assumed to lead to the accumulation of deleterious mutations, and reduced heterozygosity due to the absence of recombination. Panagrolaimid nematode species display different modes of reproduction. Sexual reproduction with distinct males and females, asexual reproduction through parthenogenesis in the genus Panagrolaimus, and hermaphroditism in Propanagrolaimus. Here, we compared genomic features of free-living nematodes in populations and species isolated from geographically distant regions to study diversity, and genome-wide differentiation under different modes of reproduction. We firstly estimated genome-wide spontaneous mutation rates in a triploid parthenogenetic Panagrolaimus, and a diploid hermaphroditic Propanagrolaimus via long-term mutation accumulation lines. Secondly, we calculated population genetic parameters including nucleotide diversity, and fixation index (F ST) between populations of asexually and sexually reproducing nematodes. Thirdly, we used phylogenetic network methods on sexually and asexually reproducing Panagrolaimus populations to understand evolutionary relationships between them. The estimated mutation rate was slightly lower for the asexual population, as expected for taxa with this reproductive mode. Natural polyploid asexual populations revealed higher nucleotide diversity. Despite their common ancestor, a gene network revealed a high level of genetic differentiation among asexual populations. The elevated heterozygosity found in the triploid parthenogens could be explained by the third genome copy. Given their tendentially lower mutation rates it can be hypothesized that this is part of the mechanism to evade Muller's ratchet. Our findings in parthenogenetic triploid nematode populations seem to challenge common expectations of evolution under asexuality.

Signification and Application of Mutator and Antimutator Phenotype-Induced Genetic Variations in Evolutionary Adaptation and Cancer Therapeutics.

Mutations present a dichotomy in their implications for cellular processes. They primarily arise from DNA replication errors or damage repair processes induced by environmental challenges. Cumulative mutations underlie genetic variations and drive evolution, yet also contribute to degenerative diseases such as cancer and aging. The mutator phenotype elucidates the heightened mutation rates observed in malignant tumors. Evolutionary adaptation, analogous to bacterial and eukaryotic systems, manifests through mutator phenotypes during changing environmental conditions, highlighting the delicate balance between advantageous mutations and their potentially detrimental consequences. Leveraging the genetic tractability of Saccharomyces cerevisiae offers unique insights into mutator phenotypes and genome instability akin to human cancers. Innovative reporter assays in yeast model organisms enable the detection of diverse genome alterations, aiding a comprehensive analysis of mutator phenotypes. Despite significant advancements, our understanding of the intricate mechanisms governing spontaneous mutation rates and preserving genetic integrity remains incomplete. This review outlines various cellular pathways affecting mutation rates and explores the role of mutator genes and mutation-derived phenotypes, particularly prevalent in malignant tumor cells. An in-depth comprehension of mutator and antimutator activities in yeast and higher eukaryotes holds promise for effective cancer control strategies.

The function of BCL11B in base excision repair contributes to its dual role as an oncogene and a haplo-insufficient tumor suppressor gene.

Genetic studies in mice and human cancers established BCL11B as a haploinsufficient tumor suppressor gene. Paradoxically, BCL11B is overexpressed in some human cancers where its knockdown is synthetic lethal. We identified the BCL11B protein in a proximity-dependent biotinylation screen performed with the DNA glycosylase NTHL1. In vitro DNA repair assays demonstrated that both BCL11B and a small recombinant BCL11B213-560 protein lacking transcription regulation potential can stimulate the enzymatic activities of two base excision repair (BER) enzymes: NTHL1 and Pol β. In cells, BCL11B is rapidly recruited to sites of DNA damage caused by laser microirradiation. BCL11B knockdown delays, whereas ectopic expression of BCL11B213-560 accelerates, the repair of oxidative DNA damage. Inactivation of one BCL11B allele in TK6 lymphoblastoid cells causes an increase in spontaneous and radiation-induced mutation rates. In turn, ectopic expression of BCL11B213-560 cooperates with the RAS oncogene in cell transformation by reducing DNA damage and cellular senescence. These findings indicate that BCL11B functions as a BER accessory factor, safeguarding normal cells from acquiring mutations. Paradoxically, it also enables the survival of cancer cells that would otherwise undergo senescence or apoptosis due to oxidative DNA damage resulting from the elevated production of reactive oxygen species.

DRAGON: Harnessing the power of DNA repair for accelerating genome evolution in Corynebacterium glutamicum

Hypermutation is a robust phenotype characterized by high elevation of spontaneous mutation rates, which has been shown to facilitate rapid adaptation to the stressful environments by hitchhiking with favorable mutations. Accumulating evidence argues that deficient DNA repair can give rise to hypermutation events in bacteria. Here, we provided a comprehensive survey of DNA repair systems to identify promising targets ensuring high DNA fidelity in Corynebacterium glutamicum. Four effective DNA repair factors, including nucS, tag, xpb, and dinP, were found to be strongly associated with the occurrence of hypermutable phenotypes, and these targets were then engineered to establish a CRISPRi-based all-in-one plasmid system for genome mutagenesis. On the basis of these findings, we presented a novel evolutionary engineering method named “DNA repair-assisted genome evolution (DRAGON)”. As a proof-of-concept, DRAGON strategy was successfully applied to facilitate rapid acquisition of microbial robustness in C. glutamicum, such as increased tolerances towards kanamycin, acidic pH and high L-serine, showing its promise and potential for rapid strain improvement. Overall, our study will offer new insights into the understanding of DNA repair and evolutionary adaptation in C. glutamicum.

Vibrio vulnificus mutation rate: an in vitro approach.

Vibrio vulnificus is a multi-host pathogenic species currently subdivided into five phylogenetic lineages (L) plus one pathovar with the ability to infect fish due to a transmissible virulence plasmid. This plasmid (or a fragment of it) has been transmitted between lineages within the species, contributing to the evolution of V. vulnificus. This study aimed to provide an experimental approximation to the V. vulnificus mutation rate by determining spontaneous mutation rates from bacterial cultures of representants of the different lineages by whole-genome sequencing. To this purpose, synonymous SNP differences, i.e., spontaneous mutation not subjected to the evolutive forces, between initial and final culture after serial growth were evaluated and used for mutation rate calculation.

Mutation rates and adaptive variation among the clinically dominant clusters of Mycobacterium abscessus

Mycobacterium abscessus (Mab) is a multidrug-resistant pathogen increasingly responsible for severe pulmonary infections. Analysis of whole-genome sequences (WGS) of Mab demonstrates dense genetic clustering of clinical isolates collected from disparate geographic locations. This has been interpreted as supporting patient-to-patient transmission, but epidemiological studies have contradicted this interpretation. Here, we present evidence for a slowing of the Mab molecular clock rate coincident with the emergence of phylogenetic clusters. We performed phylogenetic inference using publicly available WGS from 483 Mab patient isolates. We implement a subsampling approach in combination with coalescent analysis to estimate the molecular clock rate along the long internal branches of the tree, indicating a faster long-term molecular clock rate compared to branches within phylogenetic clusters. We used ancestry simulation to predict the effects of clock rate variation on phylogenetic clustering and found that the degree of clustering in the observed phylogeny is more easily explained by a clock rate slowdown than by transmission. We also find that phylogenetic clusters are enriched in mutations affecting DNA repair machinery and report that clustered isolates have lower spontaneous mutation rates in vitro. We propose that Mab adaptation to the host environment through variation in DNA repair genes affects the organism's mutation rate and that this manifests as phylogenetic clustering. These results challenge the model that phylogenetic clustering in Mab is explained by person-to-person transmission and inform our understanding of transmission inference in emerging, facultative pathogens.