Mutation and Recombination Rates Vary Across Bacterial Chromosome

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.

Ecological and evolutionary mechanisms driving within-patient emergence of antimicrobial resistance

The ecological and evolutionary mechanisms of antimicrobial resistance (AMR) emergence within patients and how these vary across bacterial infections are poorly understood. Increasingly widespread use of pathogen genome sequencing in the clinic enables a deeper understanding of these processes. In this Review, we explore the clinical evidence to support four major mechanisms of within-patient AMR emergence in bacteria: spontaneous resistance mutations; in situ horizontal gene transfer of resistance genes; selection of pre-existing resistance; and immigration of resistant lineages. Within-patient AMR emergence occurs across a wide range of host niches and bacterial species, but the importance of each mechanism varies between bacterial species and infection sites within the body. We identify potential drivers of such differences and discuss how ecological and evolutionary analysis could be embedded within clinical trials of antimicrobials, which are powerful but underused tools for understanding why these mechanisms vary between pathogens, infections and individuals. Ultimately, improving understanding of how host niche, bacterial species and antibiotic mode of action combine to govern the ecological and evolutionary mechanism of AMR emergence in patients will enable more predictive and personalized diagnosis and antimicrobial therapies.

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.

Mutation bias and adaptation in bacteria

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

Pangenome-level analysis of nucleoid-associated proteins in the Acidithiobacillia class: insights into their functional roles in mobile genetic elements biology

Mobile genetic elements (MGEs) are relevant agents in bacterial adaptation and evolutionary diversification. Stable appropriation of these DNA elements depends on host factors, among which are the nucleoid-associated proteins (NAPs). NAPs are highly abundant proteins that bind and bend DNA, altering its topology and folding, thus affecting all known cellular DNA processes from replication to expression. Even though NAP coding genes are found in most prokaryotic genomes, their functions in host chromosome biology and xenogeneic silencing are only known for a few NAP families. Less is known about the occurrence, abundance, and roles of MGE-encoded NAPs in foreign elements establishment and mobility. In this study, we used a combination of comparative genomics and phylogenetic strategies to gain insights into the diversity, distribution, and functional roles of NAPs within the class Acidithiobacillia with a special focus on their role in MGE biology. Acidithiobacillia class members are aerobic, chemolithoautotrophic, acidophilic sulfur-oxidizers, encompassing substantial genotypic diversity attributable to MGEs. Our search for NAP protein families (PFs) in more than 90 genomes of the different species that conform the class, revealed the presence of 1,197 proteins pertaining to 12 different NAP families, with differential occurrence and conservation across species. Pangenome-level analysis revealed 6 core NAP PFs that were highly conserved across the class, some of which also existed as variant forms of scattered occurrence, in addition to NAPs of taxa-restricted distribution. Core NAPs identified are reckoned as essential based on the conservation of genomic context and phylogenetic signals. In turn, various highly diversified NAPs pertaining to the flexible gene complement of the class, were found to be encoded in known plasmids or, larger integrated MGEs or, present in genomic loci associated with MGE-hallmark genes, pointing to their role in the stabilization/maintenance of these elements in strains and species with larger genomes. Both core and flexible NAPs identified proved valuable as markers, the former accurately recapitulating the phylogeny of the class, and the later, as seed in the bioinformatic identification of novel episomal and integrated mobile elements.

Results of the 2020 Genomic Proficiency Test for the network of European Union Reference Laboratory for Antimicrobial Resistance assessing whole-genome-sequencing capacities

The global surveillance and outbreak investigation of antimicrobial resistance (AMR) is amidst a paradigm shift from traditional biology to bioinformatics. This is due to developments in whole-genome-sequencing (WGS) technologies, bioinformatics tools, and reduced costs. The increased use of WGS is accompanied by challenges such as standardization, quality control (QC), and data sharing. Thus, there is global need for inter-laboratory WGS proficiency test (PT) schemes to evaluate laboratories' capacity to produce reliable genomic data. Here, we present the results of the first iteration of the Genomic PT (GPT) organized by the Global Capacity Building Group at the Technical University of Denmark in 2020. Participating laboratories sequenced two isolates and corresponding DNA of Salmonella enterica, Escherichia coli and Campylobacter coli, using WGS methodologies routinely employed at their laboratories. The participants' ability to obtain consistently good-quality WGS data was assessed based on several QC WGS metrics. A total of 21 laboratories from 21 European countries submitted WGS and meta-data. Most delivered high-quality sequence data with only two laboratories identified as overall underperforming. The QC metrics, N50 and number of contigs, were identified as good indicators for high-sequencing quality. We propose QC thresholds for N50 greater than 20 000 and 25 000 for Campylobacter coli and Escherichia coli, respectively, and number of contigs >200 bp greater than 225, 265 and 100 for Salmonella enterica, Escherichia coli and Campylobacter coli, respectively. The GPT2020 results confirm the importance of systematic QC procedures, ensuring the submission of reliable WGS data for surveillance and outbreak investigation to meet the requirements of the paradigm shift in methodology.

Signatures of selection in core and accessory genomes indicate different ecological drivers of diversification among Bacillus cereus clades

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan‐genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well‐studied Bacillus cereus sensu lato group, we identified clade‐specific core genes (present in all clade members) and then used clade‐specific allelic diversity to identify genes under purifying and diversifying selection. Clade‐specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade‐specific core or accessory genomes. Genes under purifying selection were often clade‐specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade‐specific selection, allowing ready identification of strains and their ecological niche.

Natural Chromosome-Chromid Fusion across rRNA Operons in a Burkholderiaceae Bacterium

Chromids (secondary chromosomes) in bacterial genomes that are present in addition to the main chromosome appear to be evolutionarily conserved in some specific bacterial groups. In rare cases among these groups, a small number of strains from Rhizobiales and Vibrionales were shown to possess a naturally fused single chromosome that was reported to have been generated through intragenomic homologous recombination between repeated sequences on the chromosome and chromid. Similar examples have never been reported in the family Burkholderiaceae, a well-documented group that conserves chromids. Here, an in-depth genomic characterization was performed on a Burkholderiaceae bacterium that was isolated from a soil bacterial consortium maintained on diesel fuel and mutagenic benzo[a]pyrene. This organism, Cupriavidus necator strain KK10, was revealed to carry a single chromosome with unexpectedly large size (>6.6 Mb), and results of comparative genomics with the genome of C. necator N-1T indicated that the single chromosome of KK10 was generated through fusion of the prototypical chromosome and chromid at the rRNA operons. This fusion hypothetically occurred through homologous recombination with a crossover between repeated rRNA operons on the chromosome and chromid. Some metabolic functions that were likely expressed from genes on the prototypical chromid region were indicated to be retained. If this phenomenon-the bacterial chromosome-chromid fusion across the rRNA operons through homologous recombination-occurs universally in prokaryotes, the multiple rRNA operons in bacterial genomes may not only contribute to the robustness of ribosome function, but also provide more opportunities for genomic rearrangements through frequent recombination. IMPORTANCE A bacterial chromosome that was naturally fused with the secondary chromosome, or "chromid," and presented as an unexpectedly large single replicon was discovered in the genome of Cupriavidus necator strain KK10, a biotechnologically useful member of the family Burkholderiaceae. Although Burkholderiaceae is a well-documented group that conserves chromids in their genomes, this chromosomal fusion event has not been previously reported for this family. This fusion has hypothetically occurred through intragenomic homologous recombination between repeated rRNA operons and, if so, provides novel insight into the potential of multiple rRNA operons in bacterial genomes to lead to chromosome-chromid fusion. The harsh conditions under which strain KK10 was maintained-a genotoxic hydrocarbon-enriched milieu-may have provided this genotype with a niche in which to survive.

Evaluation of Bacterial Diversity and Evolutionary Dynamics of Gut Bifidobacterium longum Isolates Obtained from Older Individuals in Hubei Province, China

Bifidobacterium longum predominates in the human gut throughout the life span, from birth to old age, and could alter the intestinal microbial population and immune function in the elderly. We investigated the intestinal bacterial diversity in the elderly, and further evaluated the genetic diversity and population structure of B. longum. The results revealed a distinct difference in gut bacterial populations between the elderly from Xiangyang and its neighboring region, Enshi city. A total of 62 bifidobacterial strains were isolated, 30 of which were found to be B. longum. The multilocus sequence typing (MLST) analysis also revealed that 437 B. longum isolates from diverse regions worldwide, including the 30 isolated in this study, could be classified into 341 sequence types (STs). They could be further clustered into 10 clonal complexes and 127 singleton STs, indicating a highly genetic diversity among B. longum isolates. Two putative clone complexes (CCs) containing the isolates from Xiangyang were found to be geographically specific, and a 213-bp recombination fragment was detected. Phylogenetic trees divided these 437 isolates into three lineages, corresponding to the three subspecies of B. longum. It is noteworthy that two isolates from the elderly were identified to be B. longum subsp. suis, while the others were B. longum subsp. longum. Together, our study characterized the intestinal bacterial diversity and evolution of B. longum in the elderly, and it could contribute to further studies on the genotyping and discrimination of B. longum.

IMPORTANCEBifidobacterium longum are common inhabitants of the human gut throughout the life span, and have been associated with health-promoting effects, yet little is known about the genotype profile and evolution of these isolates. Our study showed that there was significant difference in gut bacterial community and abundance of B. longum between the elderly from two neighboring cities. Furthermore, the possible geographically specific STs, CCs, and intraspecies recombination fragment were found among the B. longum isolates from elderly.

Special Issue: Role of Bacterial Chromatin in Environmental Sensing, Adaptation and Evolution

A typical bacterial cell is micron-sized and contains a genome several million base pairs in length [...].

Spatiotemporal Manipulation of the Mismatch Repair System of Pseudomonas putida Accelerates Phenotype Emergence

The development of complex phenotypes in industrially relevant bacteria is a major goal of metabolic engineering, which encompasses the implementation of both rational and random approaches. In the latter case, several tools have been developed toward increasing mutation frequencies, yet the precise control of mutagenesis processes in cell factories continues to represent a significant technical challenge. Pseudomonas species are endowed with one of the most efficient DNA mismatch repair (MMR) systems found in the bacterial domain. Here, we investigated if the endogenous MMR system could be manipulated as a general strategy to artificially alter mutation rates in Pseudomonas species. To bestow a conditional mutator phenotype in the platform bacterium Pseudomonas putida, we constructed inducible mutator devices to modulate the expression of the dominant-negative mutL^E36K allele. Regulatable overexpression of mutL^E36K in a broad-host-range, easy-to-cure plasmid format resulted in a transitory inhibition of the MMR machinery, leading to a significant increase (up to 438-fold) in DNA mutation frequencies and a heritable fixation of mutations in the genome. Following such an accelerated mutagenesis-followed by selection approach, three phenotypes were successfully evolved: resistance to antibiotics streptomycin and rifampicin (either individually or combined) and reversion of a synthetic uracil auxotrophy. Thus, these mutator devices could be applied to accelerate the evolution of metabolic pathways in long-term evolutionary experiments, alternating cycles of (inducible) mutagenesis coupled to selection schemes toward the desired phenotype(s).

Xenogeneic Silencing and Bacterial Genome Evolution: Mechanisms for DNA Recognition Imply Multifaceted Roles of Xenogeneic Silencers

Horizontal gene transfer (HGT) is a major driving force for bacterial evolution. To avoid the deleterious effects due to the unregulated expression of newly acquired foreign genes, bacteria have evolved specific proteins named xenogeneic silencers to recognize foreign DNA sequences and suppress their transcription. As there is considerable diversity in genomic base compositions among bacteria, how xenogeneic silencers distinguish self- from nonself DNA in different bacteria remains poorly understood. This review summarizes the progress in studying the DNA binding preferences and the underlying molecular mechanisms of known xenogeneic silencer families, represented by H-NS of Escherichia coli, Lsr2 of Mycobacterium, MvaT of Pseudomonas, and Rok of Bacillus. Comparative analyses of the published data indicate that the differences in DNA recognition mechanisms enable these xenogeneic silencers to have clear characteristics in DNA sequence preferences, which are further correlated with different host genomic features. These correlations provide insights into the mechanisms of how these xenogeneic silencers selectively target foreign DNA in different genomic backgrounds. Furthermore, it is revealed that the genomic AT contents of bacterial species with the same xenogeneic silencer family proteins are distributed in a limited range and are generally lower than those species without any known xenogeneic silencers in the same phylum/class/genus, indicating that xenogeneic silencers have multifaceted roles on bacterial genome evolution. In addition to regulating horizontal gene transfer, xenogeneic silencers also act as a selective force against the GC to AT mutational bias found in bacterial genomes and help the host genomic AT contents maintained at relatively low levels.

Pseudouridines of tRNA Anticodon Stem-Loop Have Unexpected Role in Mutagenesis in Pseudomonas sp

Pseudouridines are known to be important for optimal translation. In this study we demonstrate an unexpected link between pseudouridylation of tRNA and mutation frequency in Pseudomonas species. We observed that the lack of pseudouridylation activity of pseudouridine synthases TruA or RluA elevates the mutation frequency in Pseudomonas putida 3 to 5-fold. The absence of TruA but not RluA elevates mutation frequency also in Pseudomonas aeruginosa. Based on the results of genetic studies and analysis of proteome data, the mutagenic effect of the pseudouridylation deficiency cannot be ascribed to the involvement of error-prone DNA polymerases or malfunctioning of DNA repair pathways. In addition, although the deficiency in TruA-dependent pseudouridylation made P. putida cells more sensitive to antimicrobial compounds that may cause oxidative stress and DNA damage, cultivation of bacteria in the presence of reactive oxygen species (ROS)-scavenging compounds did not eliminate the mutator phenotype. Thus, the elevated mutation frequency in the absence of tRNA pseudouridylation could be the result of a more specific response or, alternatively, of a cumulative effect of several small effects disturbing distinct cellular functions, which remain undetected when studied independently. This work suggests that pseudouridines link the translation machinery to mutation frequency.

Nucleoid Associated Proteins: The Small Organizers That Help to Cope With Stress

The bacterial chromosome must be efficiently compacted to fit inside the small and crowded cell while remaining accessible for the protein complexes involved in replication, transcription, and DNA repair. The dynamic organization of the nucleoid is a consequence of both intracellular factors (i.e., simultaneously occurring cell processes) and extracellular factors (e.g., environmental conditions, stress agents). Recent studies have revealed that the bacterial chromosome undergoes profound topological changes under stress. Among the many DNA-binding proteins that shape the bacterial chromosome structure in response to various signals, NAPs (nucleoid associated proteins) are the most abundant. These small, basic proteins bind DNA with low specificity and can influence chromosome organization under changing environmental conditions (i.e., by coating the chromosome in response to stress) or regulate the transcription of specific genes (e.g., those involved in virulence).