Core genes can have higher recombination rates than accessory genes within global microbial populations
Abstract
Recombination is essential to microbial evolution, and is involved in the spread of antibiotic resistance, antigenic variation, and adaptation to the host niche. However, assessing the impact of homologous recombination on accessory genes which are only present in a subset of strains of a given species remains challenging due to their complex phylogenetic relationships. Quantifying homologous recombination for accessory genes (which are important for niche-specific adaptations) in comparison to core genes (which are present in all strains and have essential functions) is critical to understanding how selection acts on variation to shape species diversity and genome structures of bacteria. Here, we apply a computationally efficient, non-phylogenetic approach to measure homologous recombination rates in the core and accessory genome using >100,000 whole genome sequences from Streptococcus pneumoniae and several additional species. By analyzing diverse sets of sequence clusters, we show that core genes often have higher recombination rates than accessory genes, and for some bacterial species the associated effect sizes for these differences are pronounced. In a subset of species, we find that gene frequency and homologous recombination rate are positively correlated. For S. pneumoniae and several additional species, we find that while the recombination rate is higher for the core genome, the mutational divergence is lower, indicating that divergence-based homologous recombination barriers could contribute to differences in recombination rates between the core and accessory genome. Homologous recombination may therefore play a key role in increasing the efficiency of selection in the most conserved parts of the genome.
Data availability
Lists of SRA accession numbers corresponding to the raw reads used to build the multi-sequence alignments analyzed in this manuscript are included as Figure 2 - source data 1 and Figure 3 - source data 1. All SRA files, reference genomes, and complete genome assemblies are available through NCBI. All sequence collections used are listed in Supplementary File 5. For the PubMLST sequence collections, PubMLST was used to identify whole genome sequences (by filtering for strains in the 'Genome Collection' of each species where the sequence length is at least that of the reference genome), then the raw reads were downloaded from NCBI using their SRA numbers. Accession numbers for reference genomes used for each microbial species are also listed in Supplementary File 5.All original code has been deposited at GitHub and is publicly available. Links are given below:- https://github.com/kussell-lab/mcorr- https://github.com/kussell-lab/mcorr-clustering- https://github.com/kussell-lab/ReferenceAlignmentGenerator- https://github.com/kussell-lab/PangenomeAlignmentGenerator
Article and author information
Author details
Funding
National Institutes of Health (R01-GM097356)
- Edo Kussell
Simons Foundation (Simons Foundation Awardee of the Life Sciences Research Foundation)
- Asher Preska Steinberg
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2022, Preska Steinberg et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,767
- views
-
- 531
- downloads
-
- 19
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
- Microbiology and Infectious Disease
The sustained success of Mycobacterium tuberculosis as a pathogen arises from its ability to persist within macrophages for extended periods and its limited responsiveness to antibiotics. Furthermore, the high incidence of resistance to the few available antituberculosis drugs is a significant concern, especially since the driving forces of the emergence of drug resistance are not clear. Drug-resistant strains of Mycobacterium tuberculosis can emerge through de novo mutations, however, mycobacterial mutation rates are low. To unravel the effects of antibiotic pressure on genome stability, we determined the genetic variability, phenotypic tolerance, DNA repair system activation, and dNTP pool upon treatment with current antibiotics using Mycobacterium smegmatis. Whole-genome sequencing revealed no significant increase in mutation rates after prolonged exposure to first-line antibiotics. However, the phenotypic fluctuation assay indicated rapid adaptation to antibiotics mediated by non-genetic factors. The upregulation of DNA repair genes, measured using qPCR, suggests that genomic integrity may be maintained through the activation of specific DNA repair pathways. Our results, indicating that antibiotic exposure does not result in de novo adaptive mutagenesis under laboratory conditions, do not lend support to the model suggesting antibiotic resistance development through drug pressure-induced microevolution.
-
- Microbiology and Infectious Disease
Bacillus velezensis is a species of Bacillus that has been widely investigated because of its broad-spectrum antimicrobial activity. However, most studies on B. velezensis have focused on the biocontrol of plant diseases, with few reports on antagonizing Salmonella Typhimurium infections. In this investigation, it was discovered that B. velezensis HBXN2020, which was isolated from healthy black pigs, possessed strong anti-stress and broad-spectrum antibacterial activity. Importantly, B. velezensis HBXN2020 did not cause any adverse side effects in mice when administered at various doses (1×107, 1×108, and 1×109 CFU) for 14 days. Supplementing B. velezensis HBXN2020 spores, either as a curative or preventive measure, dramatically reduced the levels of S. Typhimurium ATCC14028 in the mice’s feces, ileum, cecum, and colon, as well as the disease activity index (DAI), in a model of infection caused by this pathogen in mice. Additionally, supplementing B. velezensis HBXN2020 spores significantly regulated cytokine levels (Tnfa, Il1b, Il6, and Il10) and maintained the expression of tight junction proteins and mucin protein. Most importantly, adding B. velezensis HBXN2020 spores to the colonic microbiota improved its stability and increased the amount of beneficial bacteria (Lactobacillus and Akkermansia). All together, B. velezensis HBXN2020 can improve intestinal microbiota stability and barrier integrity and reduce inflammation to help treat infection by S. Typhimurium.