Impact of Homologous Recombination on the Evolution of Prokaryotic Core Genomes

A marker-free genetic manipulation method for Glaesserella parasuis strains developed by alternately culturing transformants at 37°C and 30°C

BACKGROUND

Glaesserella parasuis (G. parasuis) is the causative agent of Glässer's disease, which causes significant economic losses in the swine industry. However, research on the pathogenesis of G. parasuis has been hampered by the lack of a simple and efficient marker-free knockout system.

RESULTS

In this study, a marker-free knockout system was developed for G. parasuis using a temperature-sensitive vector. By alternating the incubation of transformants at 30°C and 37°C, we optimized the screening process for this system. The system was successfully applied to knockout the KanR cassette from JS0135ΔnanH::KanR, achieving a knockout efficiency of 90% in the final round of screening. To confirm that temperature variation was a key factor, we proceeded with knocking out the nanH and apd genes in the CF7066 strain. The knockout efficiency reached up to 100%, with the shortest screening time being only four days. The knockout of the nanH gene resulted in a significant reduction in the growth vitality of the strains, while the knockout of the apd gene led to an approximate 56% improvement in the adhesion rate. Additionally, we observed that the expression of recombinant genes in transformants was higher at 30℃ than at 37℃, with the recC gene being upregulated approximately 7-fold. In contrast, there was almost no difference in the expression of recombinant genes between 30℃ and 37℃ in the wild-type strains. This discrepancy was likely due to an elevated copy number of target plasmids at 30℃, which may have resulted in the enhanced expression of recombinant genes.

CONCLUSIONS

In conclusion, this newly developed gene knockout system for G. parasuis presents a valuable tool for advancing research on this organism.

Comparative genomic analysis identifies potential adaptive variation in Mycoplasma ovipneumoniae

Mycoplasma ovipneumoniae is associated with respiratory disease in wild and domestic Caprinae globally, with wide variation in disease outcomes within and between host species. To gain insight into phylogenetic structure and mechanisms of pathogenicity for this bacterial species, we compared M. ovipneumoniae genomes for 99 samples from 6 countries (Australia, Bosnia and Herzegovina, Brazil, China, France and USA) and 4 host species (domestic sheep, domestic goats, bighorn sheep and caribou). Core genome sequences of M. ovipneumoniae assemblies from domestic sheep and goats fell into two well-supported phylogenetic clades that are divergent enough to be considered different bacterial species, consistent with each of these two clades having an evolutionary origin in separate host species. Genome assemblies from bighorn sheep and caribou also fell within these two clades, indicating multiple spillover events, most commonly from domestic sheep. Pangenome analysis indicated a high percentage (91.4 %) of accessory genes (i.e. genes found only in a subset of assemblies) compared to core genes (i.e. genes found in all assemblies), potentially indicating a propensity for this pathogen to adapt to within-host conditions. In addition, many genes related to carbon metabolism, which is a virulence factor for Mycoplasmas, showed evidence for homologous recombination, a potential signature of adaptation. The presence or absence of annotated genes was very similar between sheep and goat clades, with only two annotated genes significantly clade-associated. However, three M. ovipneumoniae genome assemblies from asymptomatic caribou in Alaska formed a highly divergent subclade within the sheep clade that lacked 23 annotated genes compared to other assemblies, and many of these genes had functions related to carbon metabolism. Overall, our results suggest that adaptation of M. ovipneumoniae has involved evolution of carbon metabolism pathways and virulence mechanisms related to those pathways. The genes involved in these pathways, along with other genes identified as potentially involved in virulence in this study, are potential targets for future investigation into a possible genomic basis for the high variation observed in disease outcomes within and between wild and domestic host species.

Genome-wide identification and expression analysis of the BZR gene family in Zanthoxylum armatum DC and functional analysis of ZaBZR1 in drought tolerance

MAIN CONCLUSION

In this study, six ZaBZRs were identified in Zanthoxylum armatum DC, and all the ZaBZRs were upregulated by abscisic acid (ABA) and drought. Overexpression of ZaBZR1 enhanced the drought tolerance of transgenic Nicotiana benthamian. Brassinosteroids (BRs) are a pivotal class of sterol hormones in plants that play a crucial role in plant growth and development. BZR (brassinazole resistant) is a crucial transcription factor in the signal transduction pathway of BRs. However, the BZR gene family members have not yet been identified in Zanthoxylum armatum DC. In this study, six members of the ZaBZR family were identified by bioinformatic methods. All six ZaBZRs exhibited multiple phosphorylation sites. Phylogenetic and collinearity analyses revealed a closest relationship between ZaBZRs and ZbBZRs located on the B subgenomes. Expression analysis revealed tissue-specific expression patterns of ZaBZRs in Z. armatum, and their promoter regions contained cis-acting elements associated with hormone response and stress induction. Additionally, all six ZaBZRs showed upregulation upon treatment after abscisic acid (ABA) and polyethylene glycol (PEG), indicating their participation in drought response. Subsequently, we conducted an extensive investigation of ZaBZR1. ZaBZR1 showed the highest expression in the root, followed by the stem and terminal bud. Subcellular localization analysis revealed that ZaBZR1 is present in the cytoplasm and nucleus. Overexpression of ZaBZR1 in transgenic Nicotiana benthamiana improved seed germination rate and root growth under drought conditions, reducing water loss rates compared to wild-type plants. Furthermore, ZaBZR1 increased proline content (PRO) and decreased malondialdehyde content (MDA), indicating improved tolerance to drought-induced oxidative stress. The transgenic plants also showed a reduced accumulation of reactive oxygen species. Importantly, ZaBZR1 up-regulated the expression of drought-related genes such as NbP5CS1, NbDREB2A, and NbWRKY44. These findings highlight the potential of ZaBZR1 as a candidate gene for enhancing drought resistance in transgenic N. benthamiana and provide insight into the function of ZaBZRs in Z. armatum.

Homologous recombination and gene-specific selection co-shape the vertical nucleotide diversity of mangrove sediment microbial populations

Mangrove sediments host a diverse array of microbial populations and are characterized by high heterogeneity along their vertical depths. However, the genetic diversity within these populations is largely unknown, hindering our understanding of their adaptive evolution across the sediment depths. To elucidate their genetic diversity, we utilized metagenome sequencing to identify 16 high-frequency microbial populations comprised of two archaea and 14 bacteria from mangrove sediment cores (0-100 cm, with 10 depths) in Qi'ao Island, China. Our analysis of the genome-wide genetic variation revealed extensive nucleotide diversity in the microbial populations. The genes involved in the transport and the energy metabolism displayed a high nucleotide diversity (HND; 0.0045-0.0195; an indicator of shared minor alleles with the microbial populations). By tracking the processes of homologous recombination, we found that each microbial population was subjected to different purification selection levels at different depths (44.12% genes). This selection resulted in significant differences in synonymous/non-synonymous mutation ratio between 0-20 and 20-100 cm layers, indicating the adaptive evolutionary process of microbial populations. Furthermore, our assessment of differentiation in the allele frequencies between these two layers showed that the functional genes involved in the metabolic processes of amino acids or cofactors were highly differential in more than half of them. Together, we showed that the nucleotide diversity of microbial populations was shaped by homologous recombination and gene-specific selection, finally resulting in the stratified differentiation occurring between 0-20 and 20-100 cm. These results enhance our cognition of the microbial adaptation mechanisms to vertical environmental changes during the sedimentation process of coastal blue carbon ecosystems.

Homologous Recombination Shapes the Architecture and Evolution of Bacterial Genomes

Homologous recombination is a key evolutionary force that varies considerably across bacterial species. However, how the landscape of homologous recombination varies across genes and within individual genomes has only been studied in a few species. Here, we used Approximate Bayesian Computation to estimate the recombination rate along the genomes of 145 bacterial species. Our results show that homologous recombination varies greatly along bacterial genomes and shapes many aspects of genome architecture and evolution. The genomic landscape of recombination presents several key signatures: rates are highest near the origin of replication in most species, patterns of recombination generally appear symmetrical in both replichores (i.e. replicational halves of circular chromosomes) and most species have genomic hotpots of recombination. Furthermore, many closely related species share conserved landscapes of recombination across orthologs indicating that recombination landscapes are conserved over significant evolutionary distances. We show evidence that recombination drives the evolution of GC-content through increasing the effectiveness of selection and not through biased gene conversion, thereby contributing to an ongoing debate. Finally, we demonstrate that the rate of recombination varies across gene function and that many hotspots of recombination are associated with adaptive and mobile regions often encoding genes involved in pathogenicity.

Evolution of homologous recombination rates across bacteria

Bacteria are nonsexual organisms but are capable of exchanging DNA at diverse degrees through homologous recombination. Intriguingly, the rates of recombination vary immensely across lineages where some species have been described as purely clonal and others as "quasi-sexual." However, estimating recombination rates has proven a difficult endeavor and estimates often vary substantially across studies. It is unclear whether these variations reflect natural variations across populations or are due to differences in methodologies. Consequently, the impact of recombination on bacterial evolution has not been extensively evaluated and the evolution of recombination rate-as a trait-remains to be accurately described. Here, we developed an approach based on Approximate Bayesian Computation that integrates multiple signals of recombination to estimate recombination rates. We inferred the rate of recombination of 162 bacterial species and one archaeon and tested the robustness of our approach. Our results confirm that recombination rates vary drastically across bacteria; however, we found that recombination rate-as a trait-is conserved in several lineages but evolves rapidly in others. Although some traits are thought to be associated with recombination rate (e.g., GC-content), we found no clear association between genomic or phenotypic traits and recombination rate. Overall, our results provide an overview of recombination rate, its evolution, and its impact on bacterial evolution.

Impact of Homologous Recombination on Core Genome Evolution and Host Adaptation of Pectobacterium parmentieri

Homologous recombination is a major force mechanism driving bacterial evolution, host adaptability, and acquisition of novel virulence traits. Pectobacterium parmentieri is a plant bacterial pathogen distributed worldwide, primarily affecting potatoes, by causing soft rot and blackleg diseases. The goal of this investigation was to understand the impact of homologous recombination on the genomic evolution of P. parmentieri. Analysis of P. parmentieri genomes using Roary revealed a dynamic pan-genome with 3,742 core genes and over 55% accessory genome variability. Bayesian population structure analysis identified 7 lineages, indicating species heterogeneity. ClonalFrameML analysis displayed 5,125 recombination events, with the lineage 4 exhibiting the highest events. fastGEAR analysis identified 486 ancestral and 941 recent recombination events ranging from 43 bp to 119 kb and 36 bp to 13.96 kb, respectively, suggesting ongoing adaptation. Notably, 11% (412 genes) of the core genome underwent recent recombination, with lineage 1 as the main donor. The prevalence of recent recombination (double compared to ancient) events implies continuous adaptation, possibly driven by global potato trade. Recombination events were found in genes involved in vital cellular processes (DNA replication, DNA repair, RNA processing, homeostasis, and metabolism), pathogenicity determinants (type secretion systems, cell-wall degrading enzymes, iron scavengers, lipopolysaccharides (LPS), flagellum, etc.), antimicrobial compounds (phenazine and colicin) and even CRISPR-Cas genes. Overall, these results emphasize the potential role of homologous recombination in P. parmentieri's evolutionary dynamics, influencing host colonization, pathogenicity, adaptive immunity, and ecological fitness.

Diversity, Distribution, and Chromosomal Rearrangements of TRIP1 Repeat Sequences in Escherichia coli

The bacterial genome contains numerous repeated sequences that greatly affect its genomic plasticity. The Escherichia coli K-12 genome contains three copies of the TRIP1 repeat sequence (TRIP1a, TRIP1b, and TRIP1c). However, the diversity, distribution, and role of the TRIP1 repeat sequence in the E. coli genome are still unclear. In this study, after screening 6725 E. coli genomes, the TRIP1 repeat was found in the majority of E. coli strains (96%: 6454/6725). The copy number and direction of the TRIP1 repeat sequence varied in each genome. Overall, 2449 genomes (36%: 2449/6725) had three copies of TRIP1 (TRIP1a, TRIP1b, and TRIP1c), which is the same as E. coli K-12. Five types of TRIP1 repeats, including two new types (TRIP1d and TRIP1e), are identified in E. coli genomes, located in 4703, 3529, 5741, 1565, and 232 genomes, respectively. Each type of TRIP1 repeat is localized to a specific locus on the chromosome. TRIP1 repeats can cause intra-chromosomal rearrangements. A total of 156 rearrangement events were identified, of which 88% (137/156) were between TRIP1a and TRIP1c. These findings have important implications for future research on TRIP1 repeats.

Genome-based taxonomy of Burkholderia sensu lato: Distinguishing closely related species

The taxonomy of Burkholderia sensu lato (s.l.) has been revisited using genome-based tools, which have helped differentiate closely related species. Many species from this group are indistinguishable through phenotypic traits and 16S rRNA gene sequence analysis. Furthermore, they also exhibit whole-genome Average Nucleotide Identity (ANI) values in the twilight zone for species circumscription (95-96%), which may impair their correct classification. In this work, we provided an updated Burkholderia s.l. taxonomy focusing on closely related species and give other recommendations for those developing genome-based taxonomy studies. We showed that a combination of ANI and digital DNA-DNA hybridization (dDDH) applying the universal cutoff values of 95% and 70%, respectively, successfully discriminates Burkholderia s.l. species. Using genome metrics with this pragmatic criterion, we demonstrated that i) Paraburkholderia insulsa should be considered a later heterotypic synonym of Paraburkholderia fungorum; ii) Paraburkholderia steynii differs from P. terrae by harboring symbiotic genes; iii) some Paraburkholderia are indeed different species based on dDDH values, albeit sharing ANI values close to 95%; iv) some Burkholderia s.l. indeed represent new species from the genomic viewpoint; iv) some genome sequences should be evaluated with care due to quality concerns.

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Ecological differences among hydrothermal vent symbioses may drive contrasting patterns of symbiont population differentiation

The intra-host composition of horizontally transmitted microbial symbionts can vary across host populations due to interactive effects of host genetics, environmental, and geographic factors. While adaptation to local habitat conditions can drive geographic subdivision of symbiont strains, it is unknown how differences in ecological characteristics among host-symbiont associations influence the genomic structure of symbiont populations. To address this question, we sequenced metagenomes of different populations of the deep-sea mussel Bathymodiolus septemdierum, which are common at Western Pacific deep-sea hydrothermal vents and show characteristic patterns of niche partitioning with sympatric gastropod symbioses. Bathymodiolus septemdierum lives in close symbiotic relationship with sulfur-oxidizing chemosynthetic bacteria but supplements its symbiotrophic diet through filter-feeding, enabling it to occupy ecological niches with little exposure to geochemical reductants. Our analyses indicate that symbiont populations associated with B. septemdierum show structuring by geographic location, but that the dominant symbiont strain is uncorrelated with vent site. These patterns are in contrast to co-occurring Alviniconcha and Ifremeria gastropod symbioses that exhibit greater symbiont nutritional dependence and occupy habitats with higher spatial variability in environmental conditions. Our results suggest that relative habitat homogeneity combined with sufficient symbiont dispersal and genomic mixing might promote persistence of similar symbiont strains across geographic locations, while mixotrophy might decrease selective pressures on the host to affiliate with locally adapted symbiont strains. Overall, these data contribute to our understanding of the potential mechanisms influencing symbiont population structure across a spectrum of marine microbial symbioses that occupy contrasting ecological niches. IMPORTANCE Beneficial relationships between animals and microbial organisms (symbionts) are ubiquitous in nature. In the ocean, microbial symbionts are typically acquired from the environment and their composition across geographic locations is often shaped by adaptation to local habitat conditions. However, it is currently unknown how generalizable these patterns are across symbiotic systems that have contrasting ecological characteristics. To address this question, we compared symbiont population structure between deep-sea hydrothermal vent mussels and co-occurring but ecologically distinct snail species. Our analyses show that mussel symbiont populations are less partitioned by geography and do not demonstrate evidence for environmental adaptation. We posit that the mussel's mixotrophic feeding mode may lower its need to affiliate with locally adapted symbiont strains, while microhabitat stability and symbiont genomic mixing likely favors persistence of symbiont strains across geographic locations. Altogether, these findings further our understanding of the mechanisms shaping symbiont population structure in marine environmentally transmitted symbioses.

The Multifaceted Role of Homologous Recombination in a Fastidious Bacterial Plant Pathogen

Homologous recombination plays a key function in the evolution of bacterial genomes. Within Xylella fastidiosa, an emerging plant pathogen with increasing host and geographic ranges, it has been suggested that homologous recombination facilitates host switching, speciation, and the development of virulence. We used 340 whole-genome sequences to study the relationship between inter- and intrasubspecific homologous recombination, random mutation, and natural selection across individual X. fastidiosa genes. Individual gene orthologs were identified and aligned, and a maximum likelihood (ML) gene tree was generated. Each gene alignment and tree pair were then used to calculate gene-wide and branch-specific r/m values (relative effect of recombination to mutation), gene-wide and branch-site nonsynonymous over synonymous substitution rates (dN/dS values; episodic selection), and branch length (as a proxy for mutation rate). The relationships between these variables were evaluated at the global level (i.e., for all genes among and within a subspecies), among specific functional classes (i.e., COGs), and between pangenome components (i.e., accessory versus core genes). Our analysis showed that r/m varied widely among genes as well as across X. fastidiosa subspecies. While r/m and dN/dS values were positively correlated in some instances (e.g., core genes in X. fastidiosa subsp. fastidiosa and both core and accessory genes in X. fastidiosa subsp. multiplex), low correlation coefficients suggested no clear biological significance. Overall, our results indicate that, in addition to its adaptive role in certain genes, homologous recombination acts as a homogenizing and a neutral force across phylogenetic clades, gene functional groups, and pangenome components. IMPORTANCE There is ample evidence that homologous recombination occurs frequently in the economically important plant pathogen Xylella fastidiosa. Homologous recombination has been known to occur among sympatric subspecies and is associated with host-switching events and virulence-linked genes. As a consequence, is it generally assumed that recombinant events in X. fastidiosa are adaptive. This mindset influences expectations of how homologous recombination acts as an evolutionary force as well as how management strategies for X. fastidiosa diseases are determined. Yet, homologous recombination plays roles beyond that of a source for diversification and adaptation. Homologous recombination can act as a DNA repair mechanism, as a means to facilitate nucleotide compositional change, as a homogenization mechanism within populations, or even as a neutral force. Here, we provide a first assessment of long-held beliefs regarding the general role of recombination in adaptation for X. fastidiosa. We evaluate gene-specific variations in homologous recombination rate across three X. fastidiosa subspecies and its relationship to other evolutionary forces (e.g., natural selection, mutation, etc.). These data were used to assess the role of homologous recombination in X. fastidiosa evolution.

Leaky barriers to gene sharing between locally co-existing coagulase-negative Staphylococcus species

Coagulase-negative Staphylococcus (CoNS) are opportunistic pathogens implicated in many human and animal infections. The evolutionary history of CoNS remains obscure because of the historical lack of recognition for their clinical importance and poor taxonomic sampling. Here, we sequenced the genomes of 191 CoNS isolates representing 15 species sampled from diseased animals diagnosed in a veterinary diagnostic laboratory. We found that CoNS are important reservoirs of diverse phages, plasmids and mobilizable genes encoding antimicrobial resistance, heavy metal resistance, and virulence. Frequent exchange of DNA between certain donor-recipient partners suggests that specific lineages act as hubs of gene sharing. We also detected frequent recombination between CoNS regardless of their animal host species, indicating that ecological barriers to horizontal gene transfer can be surmounted in co-circulating lineages. Our findings reveal frequent but structured patterns of transfer that exist within and between CoNS species, which are driven by their overlapping ecology and geographical proximity.

Despite Shared Geography, Campylobacter Isolated from Surface Water Are Genetically Distinct from Campylobacter Isolated from Chickens

We tested the hypothesis that Campylobacter isolated from chicken ceca and river water in an overlapping geographic area would share genetic information. Isolates of C. jejuni from chicken ceca were collected from a commercial slaughter plant and isolates of C. jejuni were also collected from rivers and creeks in the same watershed. Isolates were subjected to whole-genome sequencing and the data were used for core genome multilocus sequence typing (cgMLST). Cluster analysis showed that there were four distinct subpopulations, two from chickens and two from water. Calculation of fixation statistic (Fst) showed that all four subpopulations were significantly distinct. Greater than 90% of the loci were differentiated by subpopulation. Only two genes showed clear differentiation of both chicken subpopulations from both water subpopulations. Sequence fragments of the CJIE4 bacteriophage family were found frequently in the main chicken subpopulation and the water outgroup subpopulation but were sparsely found in the main water population and not at all in the chicken outgroup. CRISPR spacers that targeted the phage sequences were common in the main water subpopulation, only once in the main chicken subpopulation, and not at all in the chicken or water outgroups. Restriction enzyme genes also showed a biased distribution. These data suggest that there is little transfer of C. jejuni genetic material between chickens and nearby river water. Campylobacter differentiation according to these two sources does not show clear evidence of evolutionary selection; the differentiation is probably due to geospatial isolation, genetic drift, and the action of CRISPRs and restriction enzymes. IMPORTANCE Campylobacter jejuni causes gastroenteritis in humans, and chickens and environmental water are leading sources of infection. We tested the hypothesis that Campylobacter isolated from chicken ceca and river water in an overlapping geographic area would share genetic information. Isolates of Campylobacter were collected from water and chicken sources in the same watershed and their genomes were sequenced and analyzed. Four distinct subpopulations were found. There was no evidence of sharing genetic material between the subpopulations. Phage profiles, CRISPR profiles and restriction systems differed by subpopulation.

High genomic differentiation and limited gene flow indicate recent cryptic speciation within the genus Laspinema (cyanobacteria)

The sympatric occurrence of closely related lineages displaying conserved morphological and ecological traits is often characteristic of free-living microbes. Gene flow, recombination, selection, and mutations govern the genetic variability between these cryptic lineages and drive their differentiation. However, sequencing conservative molecular markers (e.g., 16S rRNA) coupled with insufficient population-level sampling hindered the study of intra-species genetic diversity and speciation in cyanobacteria. We used phylogenomics and a population genomic approach to investigate the extent of local genomic diversity and the mechanisms underlying sympatric speciation of Laspinema thermale. We found two cryptic lineages of Laspinema. The lineages were highly genetically diverse, with recombination occurring more frequently within than between them. That suggests the existence of a barrier to gene flow, which further maintains divergence. Genomic regions of high population differentiation harbored genes associated with possible adaptations to high/low light conditions and stress stimuli, although with a weak diversifying selection. Overall, the diversification of Laspinema species might have been affected by both genomic and ecological processes.

Core genes can have higher recombination rates than accessory genes within global microbial populations

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.

Fast and accurate metagenotyping of the human gut microbiome with GT-Pro

Single nucleotide polymorphisms (SNPs) in metagenomics are used to quantify population structure, track strains and identify genetic determinants of microbial phenotypes. However, existing alignment-based approaches for metagenomic SNP detection require high-performance computing and enough read coverage to distinguish SNPs from sequencing errors. To address these issues, we developed the GenoTyper for Prokaryotes (GT-Pro), a suite of methods to catalog SNPs from genomes and use unique k-mers to rapidly genotype these SNPs from metagenomes. Compared to methods that use read alignment, GT-Pro is more accurate and two orders of magnitude faster. Using high-quality genomes, we constructed a catalog of 104 million SNPs in 909 human gut species and used unique k-mers targeting this catalog to characterize the global population structure of gut microbes from 7,459 samples. GT-Pro enables fast and memory-efficient metagenotyping of millions of SNPs on a personal computer.

Mosaic Evolution of Beta-Barrel-Porin-Encoding Genes in Escherichia coli

Bacterial porin-encoding genes are often found under positive selection. Local recombination has also been identified in a few of them to facilitate bacterial rapid adaptation, although it remains unknown whether it is a common evolutionary mechanism for the porins or outer membrane proteins in Gram-negative bacteria. In this study, we investigated the beta-barrel (β-barrel) porin-encoding genes in Escherichia coli that were reported under positive Darwinian selection. Besides fhuA that was found with ingenic local recombination previously, we identified four other genes, i.e., lamB, ompA, ompC, and ompF, all showing the similar mosaic evolution patterns. Comparative analysis of the protein sequences disclosed a list of highly variable regions in each family, which are mostly located in the convex of extracellular loops and coinciding with the binding sites of bacteriophages. For each of the porin families, mosaic recombination leads to unique combinations of the variable regions with different sequence patterns, generating diverse protein groups. Structural modeling indicated a conserved global topology among the different porins, with the extracellular surface varying a lot due to individual or combinatorial variable regions. The conservation of global tertiary structure would ensure the channel activity, while the wide diversity of variable regions may represent selection to avoid the invasion of phages, antibiotics or immune surveillance factors. Our study identified multiple bacterial porin genes with mosaic evolution. We hypothesize that this could be generalized strategy for outer membrane proteins to both maintain normal life processes and evade the attack of unfavored factors rapidly. IMPORTANCE Microevolution studies can disclose more elaborate evolutionary mechanisms of genes, appearing especially important for genes with multifaceted function such as those encoding outer membrane proteins. However, in most cases, the gene is considered as a whole unit, and the evolutionary patterns are disclosed. Here, we report that multiple bacterial porin proteins follow mosaic evolution, with local ingenic recombination combined with spontaneous mutations based on positive Darwinian selection, and conservation for most structural regions. This could represent a common mechanism for bacterial outer membrane proteins. The variable regions within each porin family showed large coincidence with the binding sites of bacteriophages, antibiotics, and immune factors and therefore would represent effective targets for the development of new antibacterial agents or vaccines.

Evolutionary Ecology of Natural Comammox Nitrospira Populations

Microbes commonly exist in diverse and complex communities where species interact, and their genomic repertoires evolve over time. Our understanding of species interaction and evolution has increased during the last decades, but most studies of evolutionary dynamics are based on single species in isolation or in experimental systems composed of few interacting species. Here, we use the microbial ecosystem found in groundwater-fed sand filter as a model to avoid this limitation. In these open systems, diverse microbial communities experience relatively stable conditions, and the coupling between chemical and biological processes is generally well defined. Metagenomic analysis of 12 sand filters communities revealed systematic co-occurrence of at least five comammox Nitrospira species, likely promoted by low ammonium concentrations. These Nitrospira species showed intrapopulation sequence diversity, although possible clonal expansion was detected in a few abundant local comammox populations. Nitrospira species showed low homologous recombination and strong purifying selection, the latter process being especially strong in genes essential in energy metabolism. Positive selection was detected for genes related to resistance to foreign DNA and phages. We found that, compared to other habitats, groundwater-fed sand filters impose strong purifying selection and low recombination on comammox Nitrospira populations. These results suggest that evolutionary processes are more affected by habitat type than by species identity. Together, this study improves our understanding of species interaction and evolution in complex microbial communities and sheds light on the environmental dependency of evolutionary processes. IMPORTANCE Microbial species interact with each other and their environment (ecological processes) and undergo changes in their genomic repertoire over time (evolutionary processes). How these two classes of processes interact is largely unknown, especially for complex communities, as most studies of microbial evolutionary dynamics consider single species in isolation or a few interacting species in simplified experimental systems. In this study, these limitations are circumvented by examining the microbial communities found in stable and well-described groundwater-fed sand filters. Combining metagenomics and strain-level analyses, we identified the microbial interactions and evolutionary processes affecting comammox Nitrospira, a recently discovered bacterial type capable of performing the whole nitrification process. We found that abundant and co-occurrent Nitrospira populations in groundwater-fed sand filters are characterized by low recombination and strong purifying selection. In addition, by comparing these observations with those obtained from Nitrospira species inhabiting other environments, we revealed that evolutionary processes are more affected by habitat type than by species identity.

OUP accepted manuscript

Nowadays, the growing human population exacerbates the need for sustainable resources. Inspiration and achievements in nutrient production or human/animal health might emanate from microorganisms and their adaptive strategies. Here, we exemplify the benefits of lactic acid bacteria (LAB) for numerous biotechnological applications and showcase their natural transformability as a fast and robust method to hereditarily influence their phenotype/traits in fundamental and applied research contexts. We described the biogenesis of the transformation machinery and we analyzed the genome of hundreds of LAB strains exploitable for human needs to predict their transformation capabilities. Finally, we provide a stepwise rational path to stimulate and optimize natural transformation with standard and synthetic biology techniques. A comprehensive understanding of the molecular mechanisms driving natural transformation will facilitate and accelerate the improvement of bacteria with properties that serve broad societal interests.

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 [...].

Pangenome inventory of Burkholderia sensu lato, Burkholderia sensu stricto, and the Burkholderia cepacia complex reveals the uniqueness of Burkholderia catarinensis

Here the pangenome analysis of Burkholderia sensu lato (s.l.) was performed for the first time, together with an updated analysis of the pangenome of Burkholderia sensu stricto, and Burkholderia cepacia complex (Bcc) focusing on the Bcc B. catarinensis specific features of its re-sequenced genome. The pangenome of Burkholderia s.l., Burkholderia s.s., and of the Bcc are open, composed of more than 96% of accessory genes, and more than 62% of unknown genes. Functional annotations showed that secondary metabolism genes belong to the variable portion of genomes, which might explain their production of several compounds with varied bioactivities. Taken together, this work shows the great variability and uniqueness of these genomes and reveals an underexplored unknown potential in poorly characterized genes. Regarding B. catarinensis 89^T, its genome harbors genes related to hydrolases production and plant growth promotion. This draft genome will be valuable for further investigation of its biotechnological potentials.

Bacterial evolution during human infection: Adapt and live or adapt and die

Microbes are constantly evolving. Laboratory studies of bacterial evolution increase our understanding of evolutionary dynamics, identify adaptive changes, and answer important questions that impact human health. During bacterial infections in humans, however, the evolutionary parameters acting on infecting populations are likely to be much more complex than those that can be tested in the laboratory. Nonetheless, human infections can be thought of as naturally occurring in vivo bacterial evolution experiments, which can teach us about antibiotic resistance, pathogenesis, and transmission. Here, we review recent advances in the study of within-host bacterial evolution during human infection and discuss practical considerations for conducting such studies. We focus on 2 possible outcomes for de novo adaptive mutations, which we have termed "adapt-and-live" and "adapt-and-die." In the adapt-and-live scenario, a mutation is long lived, enabling its transmission on to other individuals, or the establishment of chronic infection. In the adapt-and-die scenario, a mutation is rapidly extinguished, either because it carries a substantial fitness cost, it arises within tissues that block transmission to new hosts, it is outcompeted by more fit clones, or the infection resolves. Adapt-and-die mutations can provide rich information about selection pressures in vivo, yet they can easily elude detection because they are short lived, may be more difficult to sample, or could be maladaptive in the long term. Understanding how bacteria adapt under each of these scenarios can reveal new insights about the basic biology of pathogenic microbes and could aid in the design of new translational approaches to combat bacterial infections.

Global Patterns of Recombination across Human Viruses

Viral recombination is a major evolutionary mechanism driving adaptation processes, such as the ability of host-switching. Understanding global patterns of recombination could help to identify underlying mechanisms and to evaluate the potential risks of rapid adaptation. Conventional approaches (e.g., those based on linkage disequilibrium) are computationally demanding or even intractable when sequence alignments include hundreds of sequences, common in viral data sets. We present a comprehensive analysis of recombination across 30 genomic alignments from viruses infecting humans. In order to scale the analysis and avoid the computational limitations of conventional approaches, we apply newly developed topological data analysis methods able to infer recombination rates for large data sets. We show that viruses, such as ZEBOV and MARV, consistently displayed low levels of recombination, whereas high levels of recombination were observed in Sarbecoviruses, HBV, HEV, Rhinovirus A, and HIV. We observe that recombination is more common in positive single-stranded RNA viruses than in negatively single-stranded RNA ones. Interestingly, the comparison across multiple viruses suggests an inverse correlation between genome length and recombination rate. Positional analyses of recombination breakpoints along viral genomes, combined with our approach, detected at least 39 nonuniform patterns of recombination (i.e., cold or hotspots) in 18 viral groups. Among these, noteworthy hotspots are found in MERS-CoV and Sarbecoviruses (at spike, Nucleocapsid and ORF8). In summary, we have developed a fast pipeline to measure recombination that, combined with other approaches, has allowed us to find both common and lineage-specific patterns of recombination among viruses with potential relevance in viral adaptation.

High Rates of Genome Rearrangements and Pathogenicity of Shigella spp

Shigella are pathogens originating within the Escherichia lineage but frequently classified as a separate genus. Shigella genomes contain numerous insertion sequences (ISs) that lead to pseudogenisation of affected genes and an increase of non-homologous recombination. Here, we study 414 genomes of E. coli and Shigella strains to assess the contribution of genomic rearrangements to Shigella evolution. We found that Shigella experienced exceptionally high rates of intragenomic rearrangements and had a decreased rate of homologous recombination compared to pathogenic and non-pathogenic E. coli. The high rearrangement rate resulted in independent disruption of syntenic regions and parallel rearrangements in different Shigella lineages. Specifically, we identified two types of chromosomally encoded E3 ubiquitin-protein ligases acquired independently by all Shigella strains that also showed a high level of sequence conservation in the promoter and further in the 5′-intergenic region. In the only available enteroinvasive E. coli (EIEC) strain, which is a pathogenic E. coli with a phenotype intermediate between Shigella and non-pathogenic E. coli, we found a rate of genome rearrangements comparable to those in other E. coli and no functional copies of the two Shigella-specific E3 ubiquitin ligases. These data indicate that the accumulation of ISs influenced many aspects of genome evolution and played an important role in the evolution of intracellular pathogens. Our research demonstrates the power of comparative genomics-based on synteny block composition and an important role of non-coding regions in the evolution of genomic islands.

The Origin of Niches and Species in the Bacterial World

Niches are spaces for the biological units of selection, from cells to complex communities. In a broad sense, "species" are biological units of individuation. Niches do not exist without individual organisms, and every organism has a niche. We use "niche" in the Hutchinsonian sense as an abstraction of a multidimensional environmental space characterized by a variety of conditions, both biotic and abiotic, whose quantitative ranges determine the positive or negative growth rates of the microbial individual, typically a species, but also parts of the communities of species contained in this space. Microbial organisms ("species") constantly diversify, and such diversification (radiation) depends on the possibility of opening up unexploited or insufficiently exploited niches. Niche exploitation frequently implies "niche construction," as the colonized niche evolves with time, giving rise to new potential subniches, thereby influencing the selection of a series of new variants in the progeny. The evolution of niches and organisms is the result of reciprocal interacting processes that form a single unified process. Centrifugal microbial diversification expands the limits of the species' niches while a centripetal or cohesive process occurs simultaneously, mediated by horizontal gene transfers and recombinatorial events, condensing all of the information recovered during the diversifying specialization into "novel organisms" (possible future species), thereby creating a more complex niche, where the selfishness of the new organism(s) establishes a "homeostatic power" limiting the niche's variation. Once the niche's full carrying capacity has been reached, reproductive isolation occurs, as no foreign organisms can outcompete the established population/community, thereby facilitating speciation. In the case of individualization-speciation of the microbiota, its contribution to the animal' gut structure is a type of "niche construction," the result of crosstalk between the niche (host) and microorganism(s). Lastly, there is a parallelism between the hierarchy of niches and that of microbial individuals. The increasing anthropogenic effects on the biosphere (such as globalization) might reduce the diversity of niches and bacterial individuals, with the potential emergence of highly transmissible multispecialists (which are eventually deleterious) resulting from the homogenization of the microbiosphere, a possibility that should be explored and prevented.

Evidence of gene nucleotide composition favoring replication and growth in a fastidious plant pathogen

Nucleotide composition (GC content) varies across bacteria species, genome regions, and specific genes. In Xylella fastidiosa, a vector-borne fastidious plant pathogen infecting multiple crops, GC content ranges between ∼51-52%; however, these values were gathered using limited genomic data. We evaluated GC content variations across X. fastidiosa subspecies fastidiosa (N = 194), subsp. pauca (N = 107), and subsp. multiplex (N = 39). Genomes were classified based on plant host and geographic origin; individual genes within each genome were classified based on gene function, strand, length, ortholog group, Core vs. Accessory, and Recombinant vs. Non-recombinant. GC content was calculated for each gene within each evaluated genome. The effects of genome and gene level variables were evaluated with a mixed effect ANOVA, and the marginal-GC content was calculated for each gene. Also, the correlation between gene-specific GC content vs. natural selection (dN/dS) and recombination/mutation (r/m) was estimated. Our analyses show that intra-genomic changes in nucleotide composition in X. fastidiosa are small and influenced by multiple variables. Higher AT-richness is observed in genes involved in replication and translation, and genes in the leading strand. In addition, we observed a negative correlation between high-AT and dN/dS in subsp. pauca. The relationship between recombination and GC content varied between core and accessory genes. We hypothesize that distinct evolutionary forces and energetic constraints both drive and limit these small variations in nucleotide composition.

Comparative genomics reveals broad genetic diversity, extensive recombination and nascent ecological adaptation in Micrococcus luteus

Background

Micrococcus luteus is a group of actinobacteria that is widely used in biotechnology and is being thought as an emerging nosocomial pathogen. With one of the smallest genomes of free-living actinobacteria, it is found in a wide range of environments, but intraspecies genetic diversity and adaptation strategies to various environments remain unclear. Here, comparative genomics, phylogenomics, and genome-wide association studies were used to investigate the genomic diversity, evolutionary history, and the potential ecological differentiation of the species.

Results

High-quality genomes of 66 M. luteus strains were downloaded from the NCBI GenBank database and core and pan-genome analysis revealed a considerable intraspecies heterogeneity. Phylogenomic analysis, gene content comparison, and average nucleotide identity calculation consistently indicated that the species has diverged into three well-differentiated clades. Population structure analysis further suggested the existence of an unknown ancestor or the fourth, yet unsampled, clade. Reconstruction of gene gain/loss events along the evolutionary history revealed both early events that contributed to the inter-clade divergence and recent events leading to the intra-clade diversity. We also found convincing evidence that recombination has played a key role in the evolutionary process of the species, with upto two-thirds of the core genes having been affected by recombination. Furthermore, distribution of mammal-associated strains (including pathogens) on the phylogenetic tree suggested that the last common ancestor had a free-living lifestyle, and a few recently diverged lineages have developed a mammal-associated lifestyle separately. Consistently, genome-wide association analysis revealed that mammal-associated strains from different lineages shared genes functionally relevant to the host-associated lifestyle, indicating a recent ecological adaption to the new host-associated habitats.

Conclusions

These results revealed high intraspecies genomic diversity of M. luteus and highlighted that gene gain/loss events and extensive recombination events played key roles in the genome evolution. Our study also indicated that, as a free-living species, some lineages have recently developed or are developing a mammal-associated lifestyle. This study provides insights into the mechanisms that drive the genome evolution and adaption to various environments of a bacterial species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07432-5.

Mechanisms That Shape Microbial Pangenomes

Analyses of multiple whole-genome sequences from the same species have revealed that differences in gene content can be substantial, particularly in prokaryotes. Such variation has led to the recognition of pangenomes, the complete set of genes present in a species - consisting of core genes, present in all individuals, and accessory genes whose presence is variable. Questions now arise about how pangenomes originate and evolve. We describe how gene content variation can arise as a result of the combination of several processes, including random drift, selection, gain/loss balance, and the influence of ecological and epistatic interactions. We believe that identifying the contributions of these processes to pangenomes will need novel theoretical approaches and empirical data.

Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Genomic data reveal single-nucleotide polymorphisms (SNPs) that may carry information about the evolutionary history of bacteria. However, it remains unclear what inferences about selection can be made from genomic SNP data. Bacterial species are often sampled during epidemic outbreaks or within hosts during the course of chronic infections. SNPs obtained from genomic analysis of these data are not necessarily fixed. Treating them as fixed during analysis by using measures such as the ratio of nonsynonymous to synonymous evolutionary changes (dN/dS) may lead to incorrect inferences about the strength and direction of selection. In this study, we consider data from a range of whole-genome sequencing studies of bacterial pathogens and explore patterns of nonsynonymous variation to assess whether evidence of selection can be identified by investigating SNP counts alone across multiple WGS studies. We visualize these SNP data in ways that highlight their relationship to neutral baseline expectations. These neutral expectations are based on a simple model of mutation, from which we simulate SNP accumulation to investigate how SNP counts are distributed under alternative assumptions about positive and negative selection. We compare these patterns with empirical SNP data and illustrate the general difficulty of detecting positive selection from SNP data. Finally, we consider whether SNP counts observed at the between-host population level differ from those observed at the within-host level and find some evidence that suggests that dynamics across these two scales are driven by different underlying processes.IMPORTANCE Identifying selection from SNP data obtained from whole-genome sequencing studies is challenging. Some current measures used to identify and quantify selection acting on genomes rely on fixed differences; thus, these are inappropriate for SNP data where variants are not fixed. With the increase in whole-genome sequencing studies, it is important to consider SNP data in the context of evolutionary processes. How SNPs are counted and analyzed can help in understanding mutation accumulation and trajectories of strains. We developed a tool for identifying possible evidence of selection and for comparative analysis with other SNP data. We propose a model that provides a rule-of-thumb guideline and two new visualization techniques that can be used to interpret and compare SNP data. We quantify the expected proportion of nonsynonymous SNPs in coding regions under neutrality and demonstrate its use in identifying evidence of positive and negative selection from simulations and empirical data.

Evidence of homologous recombination as a driver of diversity in Brachyspira pilosicoli

The enteric, pathogenic spirochaete Brachyspira pilosicoli colonizes and infects a variety of birds and mammals, including humans. However, there is a paucity of genomic data available for this organism. This study introduces 12 newly sequenced draft genome assemblies, boosting the cohort of examined isolates by fourfold and cataloguing the intraspecific genomic diversity of the organism more comprehensively. We used several in silico techniques to define a core genome of 1751 genes and qualitatively and quantitatively examined the intraspecific species boundary using phylogenetic analysis and average nucleotide identity, before contextualizing this diversity against other members of the genus Brachyspira. Our study revealed that an additional isolate that was unable to be species typed against any other Brachyspira lacked putative virulence factors present in all other isolates. Finally, we quantified that homologous recombination has as great an effect on the evolution of the core genome of the B. pilosicoli as random mutation (r/m=1.02). Comparative genomics has informed Brachyspira diversity, population structure, host specificity and virulence. The data presented here can be used to contribute to developing advanced screening methods, diagnostic assays and prophylactic vaccines against this zoonotic pathogen.

Impact of homologous recombination on core genome phylogenies

BACKGROUND

Core genome phylogenies are widely used to build the evolutionary history of individual prokaryote species. By using hundreds or thousands of shared genes, these approaches are the gold standard to reconstruct the relationships of large sets of strains. However, there is growing evidence that bacterial strains exchange DNA through homologous recombination at rates that vary widely across prokaryote species, indicating that core genome phylogenies might not be able to reconstruct true phylogenies when recombination rate is high. Few attempts have been made to evaluate the robustness of core genome phylogenies to recombination, but some analyses suggest that reconstructed trees are not always accurate.

RESULTS

In this study, we tested the robustness of core genome phylogenies to various levels of recombination rates. By analyzing simulated and empirical data, we observed that core genome phylogenies are relatively robust to recombination rates; nevertheless, our results suggest that many reconstructed trees are not completely accurate even when bootstrap supports are high. We found that some core genome phylogenies are highly robust to recombination whereas others are strongly impacted by it, and we identified that the robustness of core genome phylogenies to recombination is highly linked to the levels of selective pressures acting on a species. Stronger selective pressures lead to less accurate tree reconstructions, presumably because selective pressures more strongly bias the routes of DNA transfers, thereby causing phylogenetic artifacts.

CONCLUSIONS

Overall, these results have important implications for the application of core genome phylogenies in prokaryotes.

Distribution of RecBCD and AddAB recombination-associated genes among bacteria in 33 phyla

Homologous recombination plays key roles in fundamental processes such as recovery from DNA damage and in bacterial horizontal gene transfer, yet there are still open questions about the distribution of recognized components of recombination machinery among bacteria and archaea. RecBCD helicase-nuclease plays a central role in recombination among Gammaproteobacteria like Escherichia coli; while bacteria in other phyla, like the Firmicute Bacillus subtilis, use the related AddAB complex. The activity of at least some of these complexes is controlled by short DNA sequences called crossover hotspot instigator (Chi) sites. When RecBCD or AddAB complexes encounter an autologous Chi site during unwinding, they introduce a nick such that ssDNA with a free end is available to invade another duplex. If homologous DNA is present, RecA-dependent homologous recombination is promoted; if not (or if no autologous Chi site is present) the RecBCD/AddAB complex eventually degrades the DNA. We examined the distribution of recBCD and addAB genes among bacteria, and sought ways to distinguish them unambiguously. We examined bacterial species among 33 phyla, finding some unexpected distribution patterns. RecBCD and addAB are less conserved than recA, with the orthologous recB and addA genes more conserved than the recC or addB genes. We were able to classify RecB vs. AddA and RecC vs. AddB in some bacteria where this had not previously been done. We used logo analysis to identify sequence segments that are conserved, but differ between the RecBC and AddAB proteins, to help future differentiation between members of these two families.

Experimental Evolution of Bacillus subtilis Reveals the Evolutionary Dynamics of Horizontal Gene Transfer and Suggests Adaptive and Neutral Effects

Tracing evolutionary processes that lead to fixation of genomic variation in wild bacterial populations is a prime challenge in molecular evolution. In particular, the relative contribution of horizontal gene transfer (HGT) vs.de novo mutations during adaptation to a new environment is poorly understood. To gain a better understanding of the dynamics of HGT and its effect on adaptation, we subjected several populations of competent Bacillus subtilis to a serial dilution evolution on a high-salt-containing medium, either with or without foreign DNA from diverse pre-adapted or naturally salt tolerant species. Following 504 generations of evolution, all populations improved growth yield on the medium. Sequencing of evolved populations revealed extensive acquisition of foreign DNA from close Bacillus donors but not from more remote donors. HGT occurred in bursts, whereby a single bacterial cell appears to have acquired dozens of fragments at once. In the largest burst, close to 2% of the genome has been replaced by HGT. Acquired segments tend to be clustered in integration hotspots. Other than HGT, genomes also acquired spontaneous mutations. Many of these mutations occurred within, and seem to alter, the sequence of flagellar proteins. Finally, we show that, while some HGT fragments could be neutral, others are adaptive and accelerate evolution.

Efficient rational modification of non-ribosomal peptides by adenylation domain substitution

Non-ribosomal peptide synthetase (NRPS) enzymes form modular assembly-lines, wherein each module governs the incorporation of a specific monomer into a short peptide product. Modules are comprised of one or more key domains, including adenylation (A) domains, which recognise and activate the monomer substrate; condensation (C) domains, which catalyse amide bond formation; and thiolation (T) domains, which shuttle reaction intermediates between catalytic domains. This arrangement offers prospects for rational peptide modification via substitution of substrate-specifying domains. For over 20 years, it has been considered that C domains play key roles in proof-reading the substrate; a presumption that has greatly complicated rational NRPS redesign. Here we present evidence from both directed and natural evolution studies that any substrate-specifying role for C domains is likely to be the exception rather than the rule, and that novel non-ribosomal peptides can be generated by substitution of A domains alone. We identify permissive A domain recombination boundaries and show that these allow us to efficiently generate modified pyoverdine peptides at high yields. We further demonstrate the transferability of our approach in the PheATE-ProCAT model system originally used to infer C domain substrate specificity, generating modified dipeptide products at yields that are inconsistent with the prevailing dogma.

Diversity within species: interpreting strains in microbiomes

Studying within-species variation has traditionally been limited to culturable bacterial isolates and low-resolution microbial community fingerprinting. Metagenomic sequencing and technical advances have enabled culture-free, high-resolution strain and subspecies analyses at high throughput and in complex environments. This holds great scientific promise but has also led to an overwhelming number of methods and terms to describe infraspecific variation. This Review aims to clarify these advances by focusing on the diversity within bacterial and archaeal species in the context of microbiomics. We cover foundational microevolutionary concepts relevant to population genetics and summarize how within-species variation can be studied and stratified directly within microbial communities with a focus on metagenomics. Finally, we describe how common applications of within-species variation can be achieved using metagenomic data. We aim to guide the selection of appropriate terms and analytical approaches to facilitate researchers in benefiting from the increasing availability of large, high-resolution microbiome genetic sequencing data.

Development and Application of a Core Genome Multilocus Sequence Typing Scheme for the Health Care-Associated Pathogen Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen that frequently causes health care-associated infections (HAIs). Due to its metabolic diversity and ability to form biofilms, this Gram-negative nonfermenting bacterium can persist in the health care environment, which can lead to prolonged HAI outbreaks. We describe the creation of a core genome multilocus sequence typing (cgMLST) scheme to provide a stable platform for the rapid comparison of P. aeruginosa isolates using whole-genome sequencing (WGS) data. We used a diverse set of 58 complete P. aeruginosa genomes to curate a set of 4,440 core genes found in each isolate, representing ∼64% of the average genome size. We then expanded the alleles for each gene using 1,991 contig-level genome sequences. The scheme was used to analyze genomes from four historical HAI outbreaks to compare the phylogenies generated using cgMLST to those of other means (traditional MLST, pulsed-field gel electrophoresis [PFGE], and single-nucleotide variant [SNV] analysis). The cgMLST scheme provides sufficient resolution for analyzing individual outbreaks, as well as the stability for comparisons across a variety of isolates encountered in surveillance studies, making it a valuable tool for the rapid analysis of P. aeruginosa genomes.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland meadow

Soil microbial diversity is often studied from the perspective of community composition, but less is known about genetic heterogeneity within species. The relative impacts of clonal interference, gene-specific selection, and recombination in many abundant but rarely cultivated soil microbes remain unknown. Here we track genome-wide population genetic variation for 19 highly abundant bacterial species sampled from across a grassland meadow. Genomic inferences about population structure are made using the millions of sequencing reads that are assembled de novo into consensus genomes from metagenomes, as each read pair describes a short genomic sequence from a cell in each population. Genomic nucleotide identity of assembled genomes was significantly associated with local geography for over half of the populations studied, and for a majority of populations within-sample nucleotide diversity could often be as high as meadow-wide nucleotide diversity. Genes involved in metabolite biosynthesis and extracellular transport were characterized by elevated nucleotide diversity in multiple species. Microbial populations displayed varying degrees of homologous recombination and recombinant variants were often detected at 7-36% of loci genome-wide. Within multiple populations we identified genes with unusually high spatial differentiation of alleles, fewer recombinant events, elevated ratios of nonsynonymous to synonymous variants, and lower nucleotide diversity, suggesting recent selective sweeps for gene variants. Taken together, these results indicate that recombination and gene-specific selection commonly shape genetic variation in several understudied soil bacterial lineages.

Genotypic and phenotypic analyses reveal distinct population structures and ecotypes for sugar beet-associated Pseudomonas in Oxford and Auckland

Fluorescent pseudomonads represent one of the largest groups of bacteria inhabiting the surfaces of plants, but their genetic composition in planta is poorly understood. Here, we examined the population structure and diversity of fluorescent pseudomonads isolated from sugar beet grown at two geographic locations (Oxford, United Kingdom and Auckland, New Zealand). To seek evidence for niche adaptation, bacteria were sampled from three types of leaves (immature, mature, and senescent) and then characterized using a combination of genotypic and phenotypic analysis. We first performed multilocus sequence analysis (MLSA) of three housekeeping genes (gapA, gltA, and acnB) in a total of 152 isolates (96 from Oxford, 56 from Auckland). The concatenated sequences were grouped into 81 sequence types and 22 distinct operational taxonomic units (OTUs). Significant levels of recombination were detected, particularly for the Oxford isolates (rate of recombination to mutation (r/m) = 5.23 for the whole population). Subsequent ancestral analysis performed in STRUCTURE found evidence of six ancestral populations, and their distributions significantly differed between Oxford and Auckland. Next, their ability to grow on 95 carbon sources was assessed using the Biolog™ GN2 microtiter plates. A distance matrix was generated from the raw growth data (A 660) and subjected to multidimensional scaling (MDS) analysis. There was a significant correlation between substrate utilization profiles and MLSA genotypes. Both phenotypic and genotypic analyses indicated presence of a geographic structure for strains from Oxford and Auckland. Significant differences were also detected for MLSA genotypes between strains isolated from immature versus mature/senescent leaves. The fluorescent pseudomonads thus showed an ecotypic population structure, suggestive of adaptation to both geographic conditions and local plant niches.

Characterization of Burkholderia cepacia Complex Core Genome and the Underlying Recombination and Positive Selection

Recombination and positive selection are two key factors that play a vital role in pathogenic microorganisms’ population adaptation and diversification. The Burkholderia cepacia complex (Bcc) represents bacterial species with high similarity, which can cause severe infections among cases suffering from the chronic granulomatous disorder and cystic fibrosis (CF). At present, no genome-wide study has been carried out focusing on investigating the core genome of Bcc associated with the two evolutionary forces. The general characteristics of the core genome of Bcc species remain scarce as well. In this study, we explored the core orthologous genes of 116 Bcc strains using comparative genomic analysis and studied the two adaptive evolutionary forces: recombination and positive selection. We estimated 1005 orthogroups consisting entirely of single copy genes. These single copy orthologous genes in some Cluster of Orthologous Groups (COG) categories showed significant differences in the comparison of several evolutionary properties, and the encoding proteins were relatively simple and compact. Our findings showed that 5.8% of the core orthologous genes strongly supported recombination; in the meantime, 1.1% supported positive selection. We found that genes involved in protein synthesis as well as material transport and metabolism are favored by selection pressure. More importantly, homologous recombination contributed more genetic variation to a large number of genes and largely maintained the genetic cohesion in Bcc. This high level of recombination between Bcc species blurs their taxonomic boundaries, which leads Bcc species to be difficult or impossible to distinguish phenotypically and genotypically.

Impacts of local population history and ecology on the evolution of a globally dispersed pathogen

BACKGROUND

Pathogens with a global distribution face diverse biotic and abiotic conditions across populations. Moreover, the ecological and evolutionary history of each population is unique. Xylella fastidiosa is a xylem-dwelling bacterium infecting multiple plant hosts, often with detrimental effects. As a group, X. fastidiosa is divided into distinct subspecies with allopatric historical distributions and patterns of multiple introductions from numerous source populations. The capacity of X. fastidiosa to successfully colonize and cause disease in naïve plant hosts varies among subspecies, and potentially, among populations. Within Central America (i.e. Costa Rica) two X. fastidiosa subspecies coexist: the native subsp. fastidiosa and the introduced subsp. pauca. Using whole genome sequences, the patterns of gene gain/loss, genomic introgression, and genetic diversity were characterized within Costa Rica and contrasted to other X. fastidiosa populations.

RESULTS

Within Costa Rica, accessory and core genome analyses showed a highly malleable genome with numerous intra- and inter-subspecific gain/loss events. Likewise, variable levels of inter-subspecific introgression were found within and between both coexisting subspecies; nonetheless, the direction of donor/recipient subspecies to the recombinant segments varied. Some strains appeared to recombine more frequently than others; however, no group of genes or gene functions were overrepresented within recombinant segments. Finally, the patterns of genetic diversity of subsp. fastidiosa in Costa Rica were consistent with those of other native populations (i.e. subsp. pauca in Brazil).

CONCLUSIONS

Overall, this study shows the importance of characterizing local evolutionary and ecological history in the context of world-wide pathogen distribution.

Telomeric and sub-telomeric regions undergo rapid turnover within a Streptomyces population

Genome dynamics was investigated within natural populations of the soil bacterium Streptomyces. The exploration of a set of closely related strains isolated from micro-habitats of a forest soil exhibited a strong diversity of the terminal structures of the linear chromosome, i.e. terminal inverted repeats (TIRs). Large insertions, deletions and translocations could be observed along with evidence of transfer events between strains. In addition, the telomere and its cognate terminal protein complexes required for terminal replication and chromosome maintenance, were shown to be variable within the population probably reflecting telomere exchanges between the chromosome and other linear replicons (i.e., plasmids). Considering the close genetic relatedness of the strains, these data suggest that the terminal regions are prone to a high turnover due to a high recombination associated with extensive horizontal gene transfer.

Disentangling the impact of environmental and phylogenetic constraints on prokaryotic within-species diversity

Microbial organisms inhabit virtually all environments and encompass a vast biological diversity. The pangenome concept aims to facilitate an understanding of diversity within defined phylogenetic groups. Hence, pangenomes are increasingly used to characterize the strain diversity of prokaryotic species. To understand the interdependence of pangenome features (such as the number of core and accessory genes) and to study the impact of environmental and phylogenetic constraints on the evolution of conspecific strains, we computed pangenomes for 155 phylogenetically diverse species (from ten phyla) using 7,000 high-quality genomes to each of which the respective habitats were assigned. Species habitat ubiquity was associated with several pangenome features. In particular, core-genome size was more important for ubiquity than accessory genome size. In general, environmental preferences had a stronger impact on pangenome evolution than phylogenetic inertia. Environmental preferences explained up to 49% of the variance for pangenome features, compared with 18% by phylogenetic inertia. This observation was robust when the dataset was extended to 10,100 species (59 phyla). The importance of environmental preferences was further accentuated by convergent evolution of pangenome features in a given habitat type across different phylogenetic clades. For example, the soil environment promotes expansion of pangenome size, while host-associated habitats lead to its reduction. Taken together, we explored the global principles of pangenome evolution, quantified the influence of habitat, and phylogenetic inertia on the evolution of pangenomes and identified criteria governing species ubiquity and habitat specificity.

CRISPR-Cas9-mediated genome editing in vancomycin-resistant Enterococcus faecium

The Gram-positive bacterium Enterococcus faecium is becoming increasingly prevalent as a cause of hospital-acquired, antibiotic-resistant infections. A fundamental part of research into E. faecium biology relies on the ability to generate targeted mutants but this process is currently labour-intensive and time-consuming, taking 4 to 5 weeks per mutant. In this report, we describe a method relying on the high recombination rates of E. faecium and the application of the Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-Cas9 genome editing tool to more efficiently generate targeted mutants in the E. faecium chromosome. Using this tool and the multi-drug resistant clinical E. faecium strain E745, we generated a deletion mutant in the lacL gene, which encodes the large subunit of the E. faeciumβ-galactosidase. Blue/white screening using 5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside (X-gal) could be used to distinguish between the wild-type and lacL deletion mutant. We also inserted two copies of gfp into the intrinsic E. faecium macrolide resistance gene msrC to generate stable green fluorescent cells. We conclude that CRISPR-Cas9 can be used to generate targeted genome modifications in E. faecium in 3 weeks, with limited hands-on time. This method can potentially be implemented in other Gram-positive bacteria with high intrinsic recombination rates.

The Bacterial Amyloid-Like Hfq Promotes In Vitro DNA Alignment

The Hfq protein is reported to be involved in environmental adaptation and virulence of several bacteria. In Gram-negative bacteria, Hfq mediates the interaction between regulatory noncoding RNAs and their target mRNAs. Besides these RNA-related functions, Hfq is also associated with DNA and is a part of the bacterial chromatin. Its precise role in DNA structuration is, however, unclear and whether Hfq plays a direct role in DNA-related processes such as replication or recombination is controversial. In previous works, we showed that Escherichia coli Hfq, or more precisely its amyloid-like C-terminal region (CTR), induces DNA compaction into a condensed form. In this paper, we evidence a new property for Hfq; precisely we show that its CTR influences double helix structure and base tilting, resulting in a strong local alignment of nucleoprotein Hfq:DNA fibers. The significance of this alignment is discussed in terms of chromatin structuration and possible functional consequences on evolutionary processes and adaptation to environment.

Population Genetics in the Human Microbiome

While the human microbiome's structure and function have been extensively studied, its within-species genetic diversity is less well understood. However, genetic mutations in the microbiome can confer biomedically relevant traits, such as the ability to extract nutrients from food, metabolize drugs, evade antibiotics, and communicate with the host immune system. The population genetic processes by which these traits evolve are complex, in part due to interacting ecological and evolutionary forces in the microbiome. Advances in metagenomic sequencing, coupled with bioinformatics tools and population genetic models, facilitate quantification of microbiome genetic variation and inferences about how this diversity arises, evolves, and correlates with traits of both microbes and hosts. In this review, we explore the population genetic forces (mutation, recombination, drift, and selection) that shape microbiome genetic diversity within and between hosts, as well as efforts towards predictive models that leverage microbiome genetics.

Linking high GC content to the repair of double strand breaks in prokaryotic genomes

Genomic GC content varies widely among microbes for reasons unknown. While mutation bias partially explains this variation, prokaryotes near-universally have a higher GC content than predicted solely by this bias. Debate surrounds the relative importance of the remaining explanations of selection versus biased gene conversion favoring GC alleles. Some environments (e.g. soils) are associated with a high genomic GC content of their inhabitants, which implies that either high GC content is a selective adaptation to particular habitats, or that certain habitats favor increased rates of gene conversion. Here, we report a novel association between the presence of the non-homologous end joining DNA double-strand break repair pathway and GC content; this observation suggests that DNA damage may be a fundamental driver of GC content, leading in part to the many environmental patterns observed to-date. We discuss potential mechanisms accounting for the observed association, and provide preliminary evidence that sites experiencing higher rates of double-strand breaks are under selection for increased GC content relative to the genomic background.

The overall nucleotide composition of an organism’s genome varies greatly between species. Previous work has identified certain environmental factors (e.g., oxygen availability) associated with the relative number of GC bases as opposed to AT bases in the genomes of species. Many of these environments that are associated with high GC content are also associated with relatively high rates of DNA damage. We show that organisms possessing the non-homologous end-joining DNA repair pathway, which is one mechanism to repair DNA double-strand breaks, have an elevated GC content relative to expectation. We also show that certain sites on the genome that are particularly susceptible to double strand breaks have an elevated GC content. This leads us to suggest that an important underlying driver of variability in nucleotide composition across environments is the rate of DNA damage (specifically double-strand breaks) to which an organism living in each environment is exposed.

Janthinobacterium CG23_2: Comparative Genome Analysis Reveals Enhanced Environmental Sensing and Transcriptional Regulation for Adaptation to Life in an Antarctic Supraglacial Stream

As many bacteria detected in Antarctic environments are neither true psychrophiles nor endemic species, their proliferation in spite of environmental extremes gives rise to genome adaptations. Janthinobacterium sp. CG23_2 is a bacterial isolate from the Cotton Glacier stream, Antarctica. To understand how Janthinobacterium sp. CG23_2 has adapted to its environment, we investigated its genomic traits in comparison to genomes of 35 published Janthinobacterium species. While we hypothesized that genome shrinkage and specialization to narrow ecological niches would be energetically favorable for dwelling in an ephemeral Antarctic stream, the genome of Janthinobacterium sp. CG23_2 was on average 1.7 ± 0.6 Mb larger and predicted 1411 ± 499 more coding sequences compared to the other Janthinobacterium spp. Putatively identified horizontal gene transfer events contributed 0.92 Mb to the genome size expansion of Janthinobacterium sp. CG23_2. Genes with high copy numbers in the species-specific accessory genome of Janthinobacterium sp. CG23_2 were associated with environmental sensing, locomotion, response and transcriptional regulation, stress response, and mobile elements-functional categories which also showed molecular adaptation to cold. Our data suggest that genome plasticity and the abundant complementary genes for sensing and responding to the extracellular environment supported the adaptation of Janthinobacterium sp. CG23_2 to this extreme environment.