Population genomics of early events in the ecological differentiation of bacteria

Taxonomic reassessment of genomes from a divergent population of Streptococcus suis by average nucleotide identity analysis

Streptococcus spp., including the emerging zoonotic pathogen S. suis, represent a dominant bacterial population in the porcine nasopharynx. Species identification within the Streptococcus genus remains challenging and frequently ambiguous, resulting in indistinct species boundaries. By employing comparative genomic analyses, a previous study categorized S. suis into a central population and divergent lineages, based on the single nucleotide polymorphisms (SNPs) within core genes and the presence or absence of accessory genes, indicating evolutionary divergence. The divergent lineages were designated as the "out population" in this study for clarity. The 16S rRNA gene sequences of seven putative novel Streptococcus strains isolated from the throats of healthy pigs in China exhibited 100 % similarity to that of the 684_17B strain of S. suis, which clustered in the out population. This study established a threshold average nucleotide identity (ANI) value of 93.17 % for the identification of authentic S. suis. All the 645 genomes from the out population fell below this threshold, indicating that they did not belong to S. suis. Further taxonomic assessment of the 645 genomes from the out population revealed 18 clusters based on pairwise ANI comparisons, using a 92.33 % threshold, determined by pairwise ANI comparisons among the 2422 genomes from the central population of S. suis. These clusters were identified as 12 novel Streptococcus spp. (Streptococcus sp. nov-1-12) and six known Streptococcus spp. by ANI comparisons with type or reference genomes of 98 known Streptococcus spp. The study provides a useful framework for the identification of authentic S. suis and the determination of Streptococcus sp. nov.

Seepage area of the cold seep exhibits strong homogeneous selection on prokaryotic community assembly and supports high depth variability of both archaeal and bacterial communities

Haima Cold Seep is an active cold seep system in the South China Sea with a thriving biological community that plays an important role in the global carbon cycle. Despite its ecological importance, we know little about the variability of its microbial communities and depth-related structure. In this study, we comprehensively investigated the microbial depth variability and its underlying assembly process in the seepage and non-seepage areas of Haima Cold Seep. We found significant inter-group differences in archaeal (ACC) and bacterial (BCC) other than eukaryotic (ECC) community compositions, between the seepage and non-seepage areas. The seepage area had significantly more archaeal and bacterial 16S rRNA gene copies than the non-seepage area and showed greater fluctuation with depth. Anaerobic methanotrophic archaea (primarily, ANME-3) and sulfate-reducing bacteria (primarily, SEEP-SRB1) were predominant in the seepage area, suggesting that ANME-3 might work in conjunction with sulfate-reducing microorganisms to support the carbon cycle in the cold seep environment. Moreover, there were more lineages found only in specific depth ranges, supporting higher depth variability of both ACC and BCC in the seepage area than in the non-seepage area. The greater depth variability of ACC and BCC in the seepage area appeared to be primarily driven by stronger homogeneous selection imposed by environmental factors (e.g., ammonium). By contrast, eukaryotic community assembly was influenced by random processes (primarily drift) and exhibited no depth dependence. Our findings may help broaden our understanding of the cold seep ecosystem and thus provide clues for its resource utilization.IMPORTANCEMarine cold seeps are characterized by the discharge of hydrocarbons and reducing fluids. Rising geological fluids in cold seeps may act as physical transport vectors for deep biosphere microorganisms from the subsurface environment to the surface environment, and thus may influence the depth patterns of microbial community assembly. Despite the ecological importance of microbial communities in cold seeps, we have limited knowledge about their responses to environmental changes along sediment depth and the underlying processes driving these responses. Our study showed that compared with non-seepage area, seepage area exhibited stronger homogeneous selection on prokaryotic community assembly, had more depth-related specialized microorganisms, and supported higher depth variability of both archaeal and bacterial communities. Our findings may provide a theoretical basis for protection and resource utilization of the cold seep ecosystem.

Environmental Heterogeneity Drives Ecological Differentiation in Vibrio Populations Across Subtropical Marine Habitats

Elucidating how environmental gradients structure bacterial communities remains fundamental to microbial ecology. We investigated Vibrio population dynamics across contrasting subtropical marine environments in Hong Kong over a year period. Using an integrated approach combining cultivation techniques with molecular analyses of Hsp60 and 16S rRNA genes, we characterised the population structure between a coastal site (Clear Water Bay) and an estuarine site (Deep Bay). The estuarine environment consistently harboured higher Vibrio abundances (104-107 copies/mL) compared to coastal waters (102-104 copies/mL), with significantly greater phylogenetic diversity. Multivariate analyses revealed salinity as the primary driver of community differentiation between sites, while temperature governed seasonal succession patterns. Phylogenetic analysis of 1521 Vibrio isolates identified three distinct ecological groups corresponding to specific temperature-salinity niches, with evidence of habitat-specific thermal adaptations among closely related strains. Experimental characterisation of thermal performance curves confirmed physiological differentiation between warm- and cool-temperature adapted strains despite high genetic similarity (> 97% Hsp60 gene sequence identity). Several abundant species detected via amplicon sequencing (including V. navarrensis and V. mimicus) displayed site-specific ecotypes but remained uncultivated, highlighting methodological constraints in community characterisation. Our findings demonstrate how environmental heterogeneity drives fine-scale ecological differentiation in Vibrio populations, providing insights into mechanisms of bacterial adaptation in dynamic marine environments.

Context matters: assessing the impacts of genomic background and ecology on microbial biosynthetic gene cluster evolution

Encoded within many microbial genomes, biosynthetic gene clusters (BGCs) underlie the synthesis of various secondary metabolites that often mediate ecologically important functions. Several studies and bioinformatics methods developed over the past decade have advanced our understanding of both microbial pangenomes and BGC evolution. In this minireview, we first highlight challenges in broad evolutionary analysis of BGCs, including delineation of BGC boundaries and clustering of BGCs across genomes. We further summarize key findings from microbial comparative genomics studies on BGC conservation across taxa and habitats and discuss the potential fitness effects of BGCs in different settings. Afterward, recent research showing the importance of genomic context on the production of secondary metabolites and the evolution of BGCs is highlighted. These studies draw parallels to recent, broader, investigations on gene-to-gene associations within microbial pangenomes. Finally, we describe mechanisms by which microbial pangenomes and BGCs evolve, ranging from the acquisition or origination of entire BGCs to micro-evolutionary trends of individual biosynthetic genes. An outlook on how expansions in the biosynthetic capabilities of some taxa might support theories that open pangenomes are the result of adaptive evolution is also discussed. We conclude with remarks about how future work leveraging longitudinal metagenomics across diverse ecosystems is likely to significantly improve our understanding on the evolution of microbial genomes and BGCs.

Using core genome and machine learning for serovar prediction in Salmonella enterica subspecies I strains

This study presents a dual investigation of Salmonella enterica subspecies I, focusing on serovar prediction and core genome characteristics. We utilized two large genomic datasets (panX and NCBI Pathogen Detection) to test machine learning methods for predicting Salmonella serovars based on genomic differences. Among the four tested algorithms, the Random Forest model demonstrated higher performance, achieving 90.3% accuracy with the panX dataset and 95.3% with the NCBI dataset, particularly effective when trained on >50% of available data. When combined with hierarchical clustering validation, our approach achieved 100% prediction accuracy on the simulated data. Parallel analysis of panX core genome characteristics revealed that pathogenicity-related genes (including sseA, invA, mgtC, phoP, phoQ, and sitA) exhibited similar phylogenetic topologies distinct from the core genome phylogenetic tree, suggesting shared evolutionary histories. Notably, all identified core antibiotic resistance genes and virulence factors showed evidence of negative selection, indicating their essential role in bacterial survival. This study not only presents a promising machine learning-based alternative for Salmonella serovar classification, particularly valuable when analyzing newly identified serovars alongside known reference strains but also provides insights into the evolutionary dynamics of core virulence-associated genes, contributing to our understanding of Salmonella genomic architecture and pathogenicity.

Relating ecological diversity to genetic discontinuity across bacterial species

BACKGROUND

Genetic discontinuity represents abrupt breaks in genomic identity among species. Advances in genome sequencing have enhanced our ability to track and characterize genetic discontinuity in bacterial populations. However, exploring the degree to which bacterial diversity exists as a continuum or sorted into discrete and readily defined species remains a challenge in microbial ecology. Here, we aim to quantify the genetic discontinuity ( δ ) and investigate how this metric is related to ecology.

RESULTS

We harness a dataset comprising 210,129 genomes to systematically explore genetic discontinuity patterns across several distantly related species, finding clear breakpoints which vary depending on the taxa in question. By delving into pangenome characteristics, we uncover a significant association between pangenome saturation and genetic discontinuity. Closed pangenomes are associated with more pronounced breaks, exemplified by Mycobacterium tuberculosis. Additionally, through a machine learning approach, we detect key features such as gene conservation patterns and functional annotations that significantly impact genetic discontinuity prediction.

CONCLUSIONS

Our study clarifies bacterial genetic patterns and their ecological impacts, enhancing the delineation of species boundaries and deepening our understanding of microbial diversity.

Coevolution of marine phytoplankton and Alteromonas bacteria in response to pCO2 and coculture

As a result of human activity, Earth's atmosphere and climate are changing at an unprecedented pace. Models based on short-term experiments predict major changes will occur in marine phytoplankton communities in the future ocean, but rarely consider how evolution or interactions with other microbes may influence these changes. Here, we experimentally evolved several phytoplankton in coculture with a heterotrophic bacterium, Alteromonas sp. EZ55, under either present-day or predicted future pCO2 conditions. Growth rates of phytoplankton generally increased over time under both conditions, but only Thalassiosira oceanica had evidence of a growth rate tradeoff in the ancestral environment after evolution at elevated pCO2. The growth defects observed in ancestral Prochlorococcus cultures at elevated pCO2 and in axenic culture were diminished after evolution, possibly due to regulatory mutations in antioxidant genes. Except for Prochlorococcus, mutational profiles suggested phytoplankton experienced primarily purifying selection, but most Alteromonas lineages showed evidence of directional selection, where evolution appeared to favor a metabolic switch between growth on small organic acids with cyanobacteria versus catabolism of more complex carbon substrates with eukaryotic phytoplankton. Evolved Alteromonas were also poorer "helpers" for Prochlorococcus, consistent with that interaction being a competitive Black Queen process rather than a true mutualism. This work provides new insights on how phytoplankton will respond to increased pCO2 and on the evolutionary mechanisms governing phytoplankton:bacteria interactions. It also clearly demonstrates that both evolution and interspecies interactions must be considered to predict future marine biogeochemistry.

How Big Is Big? The Effective Population Size of Marine Bacteria

Genome-reduced bacteria constitute most of the cells in surface-ocean bacterioplankton communities. Their extremely large census population sizes (N c) have been unfoundedly translated to huge effective population sizes (N e)-the size of an ideal population carrying as much neutral genetic diversity as the actual population. As N e scales inversely with the strength of genetic drift, constraining the magnitude of N e is key to evaluating whether natural selection can overcome the power of genetic drift to drive evolutionary events. Determining the N e of extant species requires measuring the genomic mutation rate, a challenging step for most genome-reduced bacterioplankton lineages. Results for genome-reduced Prochlorococcus and CHUG are surprising-their N e values are an order of magnitude lower than those of less abundant lineages carrying large genomes, such as Ruegeria and Vibrio. As bacterioplankton genome reduction commonly occurred in the distant past, appreciating their population genetic mechanisms requires constraining their ancient N e values by other methods.

Diversity and biogeography of the bacterial microbiome in glacier-fed streams

The rapid melting of mountain glaciers and the vanishing of their streams is emblematic of climate change1,2. Glacier-fed streams (GFSs) are cold, oligotrophic and unstable ecosystems in which life is dominated by microbial biofilms2,3. However, current knowledge on the GFS microbiome is scarce4,5, precluding an understanding of its response to glacier shrinkage. Here, by leveraging metabarcoding and metagenomics, we provide a comprehensive survey of bacteria in the benthic microbiome across 152 GFSs draining the Earth's major mountain ranges. We find that the GFS bacterial microbiome is taxonomically and functionally distinct from other cryospheric microbiomes. GFS bacteria are diverse, with more than half being specific to a given mountain range, some unique to single GFSs and a few cosmopolitan and abundant. We show how geographic isolation and environmental selection shape their biogeography, which is characterized by distinct compositional patterns between mountain ranges and hemispheres. Phylogenetic analyses furthermore uncovered microdiverse clades resulting from environmental selection, probably promoting functional resilience and contributing to GFS bacterial biodiversity and biogeography. Climate-induced glacier shrinkage puts this unique microbiome at risk. Our study provides a global reference for future climate-change microbiology studies on the vanishing GFS ecosystem.

Co-evolution and Gene Transfers Drive Speciation Patterns in Host-Associated Bacteria

Microbial communities that maintain symbiotic relationships with animals evolve by adapting to the specific environmental niche provided by their host, yet understanding their patterns of speciation remains challenging. Whether bacterial speciation occurs primarily through allopatric or sympatric processes remains an open question. In addition, patterns of DNA transfers, which are pervasive in bacteria, are more constrained in a closed host-gut system. Eusocial bees have co-evolved with their specialized microbiota for over 85 million years, constituting a simple and valuable system to study the complex dynamics of host-associated microbial interactions. Here, we studied the patterns of speciation and evolution of seven specialized gut bacteria from three clades of eusocial bee species: western honey bees, eastern honey bees, and bumblebees. We conducted genomic analyses to infer species delineation relative to the patterns of homologous recombination (HR), and horizontal gene transfer (HGT). The studied bacteria presented various modes of evolution and speciation relative to their hosts, but some trends were consistent across all of them. We observed a clear interruption of HR between bacteria inhabiting different bee hosts, which is consistent with a mechanism of allopatric speciation, but we also identified interruptions of HR within hosts, suggesting recent or ongoing sympatric speciation. In contrast to HR, we observed that HGT events were not constrained by species borders. Overall, our findings show that in host-associated bacterial populations, patterns of HR and HGT have different impacts on speciation patterns, which are driven by both allopatric and sympatric speciation processes.

Pervasive selective sweeps across human gut microbiomes

The human gut microbiome is composed of a highly diverse consortia of species which are continually evolving within and across hosts. The ability to identify adaptations common to many human gut microbiomes would not only reveal shared selection pressures across hosts, but also key drivers of functional differentiation of the microbiome that may affect community structure and host traits. However, to date there has not been a systematic scan for adaptations that have spread across human gut microbiomes. Here, we develop a novel selection scan statistic named the integrated Linkage Disequilibrium Score (iLDS) that can detect the spread of adaptive haplotypes across host microbiomes via migration and horizontal gene transfer. Specifically, iLDS leverages signals of hitchhiking of deleterious variants with the beneficial variant. Application of the statistic to ~30 of the most prevalent commensal gut species from 24 populations around the world revealed more than 300 selective sweeps across species. We find an enrichment for selective sweeps at loci involved in carbohydrate metabolism-potentially indicative of adaptation to features of host diet-and we find that the targets of selection significantly differ between Westernized and non-Westernized populations. Underscoring the potential role of diet in driving selection, we find a selective sweep absent from non-Westernized populations but ubiquitous in Westernized populations at a locus known to be involved in the metabolism of maltodextrin, a synthetic starch that has recently become a widespread component of Western diets. In summary, we demonstrate that selective sweeps across host microbiomes are a common feature of the evolution of the human gut microbiome, and that targets of selection may be strongly impacted by host diet.

Microbial species and intraspecies units exist and are maintained by ecological cohesiveness coupled to high homologous recombination

Recent genomic analyses have revealed that microbial communities are predominantly composed of persistent, sequence-discrete species and intraspecies units (genomovars), but the mechanisms that create and maintain these units remain unclear. By analyzing closely-related isolate genomes from the same or related samples and identifying recent recombination events using a novel bioinformatics methodology, we show that high ecological cohesiveness coupled to frequent-enough and unbiased (i.e., not selection-driven) horizontal gene flow, mediated by homologous recombination, often underlie these diversity patterns. Ecological cohesiveness was inferred based on greater similarity in temporal abundance patterns of genomes of the same vs. different units, and recombination was shown to affect all sizable segments of the genome (i.e., be genome-wide) and have two times or greater impact on sequence evolution than point mutations. These results were observed in both Salinibacter ruber, an environmental halophilic organism, and Escherichia coli, the model gut-associated organism and an opportunistic pathogen, indicating that they may be more broadly applicable to the microbial world. Therefore, our results represent a departure compared to previous models of microbial speciation that invoke either ecology or recombination, but not necessarily their synergistic effect, and answer an important question for microbiology: what a species and a subspecies are.

Case report of disseminated borrelial lymphocytoma with isolation of Borrelia burgdorferi sensu stricto in chronic lymphatic leukemia stage Binet A-an 11 year follow up

We report a rare manifestation of cutaneous borreliosis in a patient with pre-existing malignant lymphoproliferative disease, in particular chronic lymphocytic B cell leukemia (B-CLL). The patient's cutaneous lesions were initially diagnosed histologically as leukemia cutis. Distribution pattern of the skin lesions were in typical localizations for borrelial lymphocytoma. Borrelia burgdorferi sensu stricto was isolated and cultured from two sites (ear, mammilla). Antibiotic therapy improved the cutaneous lesions and the general condition of the patient. However, a second round of antibiotic therapy was required to resolve the lesions. At eleven years of follow-up the patient's skin was clear and she still had a stable condition of B-CLL without chemotherapy. In conclusion, the patient suffered from Lyme borreliosis (Borrelia lymphocytoma) and the cutaneous symptoms were aggravated by the underlying condition of chronic B-CLL condition.

Origin of biogeographically distinct ecotypes during laboratory evolution

Resource partitioning is central to the incredible productivity of microbial communities, including gigatons in annual methane emissions through syntrophic interactions. Previous work revealed how a sulfate reducer (Desulfovibrio vulgaris, Dv) and a methanogen (Methanococcus maripaludis, Mm) underwent evolutionary diversification in a planktonic context, improving stability, cooperativity, and productivity within 300-1000 generations. Here, we show that mutations in just 15 Dv and 7 Mm genes within a minimal assemblage of this evolved community gave rise to co-existing ecotypes that were spatially enriched within a few days of culturing in a fluidized bed reactor. The spatially segregated communities partitioned resources in the simulated subsurface environment, with greater lactate utilization by attached Dv but partial utilization of resulting H2 by low affinity hydrogenases of Mm in the same phase. The unutilized H2 was scavenged by high affinity hydrogenases of planktonic Mm, producing copious amounts of methane. Our findings show how a few mutations can drive resource partitioning amongst niche-differentiated ecotypes, whose interplay synergistically improves productivity of the entire mutualistic community.

Barley farmland harbors a highly homogeneous soil bacterial community compared to wild ecosystems in the Qinghai-Xizang Plateau

Introduction

Understanding patterns and processes of microbial biogeography in soils is important for monitoring ecological responses to human activities, particularly in ecologically vulnerable areas such as the Qinghai-Xizang Plateau. Highland barley is the staple food of local people and has mainly been cultivated along the Yarlung Zangbo River valley in Xizang.

Methods

Here we investigated soil bacterial communities from 33 sampling sites of highland barley farmland in this region and compared them to those from wild ecosystems including alpine tundra, meadow, forest, and swamp. Additionally, the effects of environmental factors on bacterial communities, as well as the relative importance of stochastic and deterministic processes in shaping the beta diversity of soil bacterial communities in alpine ecosystems were assessed.

Results

In contrast to soils of wild ecosystems, these farmland samples harbored a highly homogeneous bacterial community without significant correlations with geographic, elevation, and edaphic distances. Discriminant bacterial taxa identified for farmland samples belong to Acidobacteria, with Acidobacteria Gp4 as the dominant clade. Although Acidobacteria were the most abundant members in all ecosystems, characterized bacterial taxa of meadow and forest were members of other phyla such as Proteobacteria and Verrucomicrobia. pH and organic matter were major edaphic attributes shaping these observed patterns across ecosystems. Null model analyses revealed that the deterministic assembly was dominant in bacterial communities in highland barley farmland and tundra soils, whereas stochastic assembly also contributed a large fraction to the assembly of bacterial communities in forest, meadow and swamp soils.

Discussion

These findings provide an insight into the consequences of human activities and agricultural intensification on taxonomic homogenization of soil bacterial communities in the Qinghai-Xizang Plateau.

Genome streamlining in Parcubacteria transitioning from soil to groundwater

BACKGROUND

To better understand the influence of habitat on the genetic content of bacteria, with a focus on members of Candidate Phyla Radiation (CPR) bacteria, we studied the effects of transitioning from soil via seepage waters to groundwater on genomic composition of ultra-small Parcubacteria, the dominating CPR class in seepage waters, using genome resolved metagenomics.

RESULTS

Bacterial metagenome-assembled genomes (MAGs), (318 total, 32 of Parcubacteria) were generated from seepage waters and compared directly to groundwater counterparts. The estimated average genome sizes of members of major phyla Proteobacteria, Bacteroidota and Cand. Patescibacteria (Candidate Phyla Radiation - CPR bacteria) were significantly higher in soil-seepage water as compared to their groundwater counterparts. Seepage water Parcubacteria (Paceibacteria) exhibited 1.18-fold greater mean genome size and 2-fold lower mean proportion of pseudogenes than those in groundwater. Bacteroidota and Proteobacteria also showed a similar trend of reduced genomes in groundwater compared to seepage. While exploring gene loss and adaptive gains in closely related CPR lineages in groundwater, we identified a membrane protein, and a lipoglycopeptide resistance gene unique to a seepage Parcubacterium genome. A nitrite reductase gene was also identified and was unique to the groundwater Parcubacteria genomes, likely acquired from other planktonic microbes via horizontal gene transfer.

CONCLUSIONS

Overall, our data suggest that bacteria in seepage waters, including ultra-small Parcubacteria, have significantly larger genomes and higher metabolic enrichment than their groundwater counterparts, highlighting possible genome streamlining of the latter in response to habitat selection in an oligotrophic environment.

Adaptations to nitrogen availability drive ecological divergence of chemosynthetic symbionts

Bacterial symbionts, with their shorter generation times and capacity for horizontal gene transfer (HGT), play a critical role in allowing marine organisms to cope with environmental change. The closure of the Isthmus of Panama created distinct environmental conditions in the Tropical Eastern Pacific (TEP) and Caribbean, offering a "natural experiment" for studying how closely related animals evolve and adapt under environmental change. However, the role of bacterial symbionts in this process is often overlooked. We sequenced the genomes of endosymbiotic bacteria in two sets of sister species of chemosymbiotic bivalves from the genera Codakia and Ctena (family Lucinidae) collected on either side of the Isthmus, to investigate how differing environmental conditions have influenced the selection of symbionts and their metabolic capabilities. The lucinid sister species hosted different Candidatus Thiodiazotropha symbionts and only those from the Caribbean had the genetic potential for nitrogen fixation, while those from the TEP did not. Interestingly, this nitrogen-fixing ability did not correspond to symbiont phylogeny, suggesting convergent evolution of nitrogen fixation potential under nutrient-poor conditions. Reconstructing the evolutionary history of the nifHDKT operon by including other lucinid symbiont genomes from around the world further revealed that the last common ancestor (LCA) of Ca. Thiodiazotropha lacked nif genes, and populations in oligotrophic habitats later re-acquired the nif operon through HGT from the Sedimenticola symbiont lineage. Our study suggests that HGT of the nif operon has facilitated niche diversification of the globally distributed Ca. Thiodiazotropha endolucinida species clade. It highlights the importance of nitrogen availability in driving the ecological diversification of chemosynthetic symbiont species and the role that bacterial symbionts may play in the adaptation of marine organisms to changing environmental conditions.

Genomic language model predicts protein co-regulation and function

Deciphering the relationship between a gene and its genomic context is fundamental to understanding and engineering biological systems. Machine learning has shown promise in learning latent relationships underlying the sequence-structure-function paradigm from massive protein sequence datasets. However, to date, limited attempts have been made in extending this continuum to include higher order genomic context information. Evolutionary processes dictate the specificity of genomic contexts in which a gene is found across phylogenetic distances, and these emergent genomic patterns can be leveraged to uncover functional relationships between gene products. Here, we train a genomic language model (gLM) on millions of metagenomic scaffolds to learn the latent functional and regulatory relationships between genes. gLM learns contextualized protein embeddings that capture the genomic context as well as the protein sequence itself, and encode biologically meaningful and functionally relevant information (e.g. enzymatic function, taxonomy). Our analysis of the attention patterns demonstrates that gLM is learning co-regulated functional modules (i.e. operons). Our findings illustrate that gLM's unsupervised deep learning of the metagenomic corpus is an effective and promising approach to encode functional semantics and regulatory syntax of genes in their genomic contexts and uncover complex relationships between genes in a genomic region.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

Microbial evolution-An under-appreciated driver of soil carbon cycling

Although substantial advances in predicting the ecological impacts of global change have been made, predictions of the evolutionary impacts have lagged behind. In soil ecosystems, microbes act as the primary energetic drivers of carbon cycling; however, microbes are also capable of evolving on timescales comparable to rates of global change. Given the importance of soil ecosystems in global carbon cycling, we assess the potential impact of microbial evolution on carbon-climate feedbacks in this system. We begin by reviewing the current state of knowledge concerning microbial evolution in response to global change and its specific effect on soil carbon dynamics. Through this integration, we synthesize a roadmap detailing how to integrate microbial evolution into ecosystem biogeochemical models. Specifically, we highlight the importance of microscale mechanistic soil carbon models, including choosing an appropriate evolutionary model (e.g., adaptive dynamics, quantitative genetics), validating model predictions with 'omics' and experimental data, scaling microbial adaptations to ecosystem level processes, and validating with ecosystem-scale measurements. The proposed steps will require significant investment of scientific resources and might require 10-20 years to be fully implemented. However, through the application of multi-scale integrated approaches, we will advance the integration of microbial evolution into predictive understanding of ecosystems, providing clarity on its role and impact within the broader context of environmental change.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Horizontal gene transfer is predicted to overcome the diversity limit of competing microbial species

Natural microbial ecosystems harbor substantial diversity of competing species. Explaining such diversity is challenging, because in classic theories it is extremely infeasible for a large community of competing species to stably coexist in homogeneous environments. One important aspect mostly overlooked in these theories, however, is that microbes commonly share genetic materials with their neighbors through horizontal gene transfer (HGT), which enables the dynamic change of species growth rates due to the fitness effects of the mobile genetic elements (MGEs). Here, we establish a framework of species competition by accounting for the dynamic gene flow among competing microbes. Combining theoretical derivation and numerical simulations, we show that in many conditions HGT can surprisingly overcome the biodiversity limit predicted by the classic model and allow the coexistence of many competitors, by enabling dynamic neutrality of competing species. In contrast with the static neutrality proposed by previous theories, the diversity maintained by HGT is highly stable against random perturbations of microbial fitness. Our work highlights the importance of considering gene flow when addressing fundamental ecological questions in the world of microbes and has broad implications for the design and engineering of complex microbial consortia.

An ANI gap within bacterial species that advances the definitions of intra-species units

IMPORTANCE

Bacterial strains and clonal complexes are two cornerstone concepts for microbiology that remain loosely defined, which confuses communication and research. Here we identify a natural gap in genome sequence comparisons among isolate genomes of all well-sequenced species that has gone unnoticed so far and could be used to more accurately and precisely define these and related concepts compared to current methods. These findings advance the molecular toolbox for accurately delineating and following the important units of diversity within prokaryotic species and thus should greatly facilitate future epidemiological and micro-diversity studies across clinical and environmental settings.

Turtle species and ecology drive carapace microbiome diversity in three seasonally interconnected wetland habitats

Different species of freshwater turtles exhibit primary behaviours ranging from aerial basking to benthic bottom-walking, cycle between wet and dry conditions at different time intervals, and undertake short-distance overland movements between aquatic habitats. These behaviours in turn may impact the accumulation of microbes on external shell surfaces of turtles and provide novel niches for differentiation of microbial communities. We assessed microbial diversity using 16S and 18S rRNA metabarcoding on carapace surfaces of six species of freshwater turtles residing in three adjacent and seasonally interconnected wetland habitats in southeast Oklahoma (United States). Communities were highly diverse, with nearly 4200 prokaryotic and 500 micro-eukaryotic amplicon sequence variants recovered, and included taxa previously reported as common or differentially abundant on turtle shells. The 16S rRNA alpha diversity tended to be highest for two species of benthic turtles, while 18S rRNA alpha diversity was highest for two basking and one shallow-water benthic species. Beta diversity of communities was more strongly differentiated by turtle species than by collection site, and ordination patterns were largely reflective of turtle species' primary habits (i.e. benthic, basking, or benthic-basking). Our data support that freshwater turtles could play a role in microbial ecology and evolution in freshwater habitats and warrant additional exploration including in areas with high native turtle diversity and inter-habitat turtle movements.

Geographic population structure and distinct intra-population dynamics of globally abundant freshwater bacteria

Implications of geographic separation and temporal dynamics on the evolution of free-living bacterial species are widely unclear. However, the vast amount of metagenome sequencing data generated during the last decades from various habitats around the world provides an unprecedented opportunity for such investigations. Here, we exploited publicly available and new freshwater metagenomes in combination with the genomes of abundant freshwater bacteria to reveal geographic and temporal population structure. We focused on species that were detected across broad geographic ranges at high enough sequence coverage for meaningful population genomic analyses, associated with the predominant freshwater taxa acI, LD12, Polynucleobacter, and Candidatus Methylopumilus. Despite the broad geographic ranges, each species appeared as a sequence-discrete cluster, in contrast to abundant marine taxa, for which continuous diversity structures were reported on a global scale. Population differentiation increased significantly with spatial distance in all species, but notable dispersal barriers (e.g. oceanic) were not apparent. Yet, the different species showed contrasting rates of geographic divergence and strikingly different intra-population dynamics in time series within individual habitats. The change in an LD12 population over 7 years was minor (FST = 0.04) compared to differentiation between lakes, whereas a Polynucleobacter population displayed strong changes within merely 2 months (FST up to 0.54), similar in scale to differentiation between populations separated by thousands of kilometers. The slowly and steadily evolving LD12 population showed high strain diversity, whereas the dynamic Polynucleobacter population exhibited alternating clonal expansions of mostly two strains only. Based on the contrasting population structures, we propose distinct models of speciation.

Sequence-discrete species for prokaryotes and other microbes: A historical perspective and pending issues

Whether prokaryotes, and other microorganisms, form distinct clusters that can be recognized as species remains an issue of paramount theoretical as well as practical consequence in identifying, regulating, and communicating about these organisms. In the past decade, comparisons of thousands of genomes of isolates and hundreds of metagenomes have shown that prokaryotic diversity may be predominantly organized in such sequence-discrete clusters, albeit organisms of intermediate relatedness between the identified clusters are also frequently found. Accumulating evidence suggests, however, that the latter "intermediate" organisms show enough ecological and/or functional distinctiveness to be considered different species. Notably, the area of discontinuity between clusters often-but not always-appears to be around 85%-95% genome-average nucleotide identity, consistently among different taxa. More recent studies have revealed remarkably similar diversity patterns for viruses and microbial eukaryotes as well. This high consistency across taxa implies a specific mechanistic process that underlies the maintenance of the clusters. The underlying mechanism may be a substantial reduction in the efficiency of homologous recombination, which mediates (successful) horizontal gene transfer, around 95% nucleotide identity. Deviations from the 95% threshold (e.g., species showing lower intraspecies diversity) may be caused by ecological differentiation that imposes barriers to otherwise frequent gene transfer. While this hypothesis that clusters are driven by ecological differentiation coupled to recombination frequency (i.e., higher recombination within vs. between groups) is appealing, the supporting evidence remains anecdotal. The data needed to rigorously test the hypothesis toward advancing the species concept are also outlined.

Vibrio cyclitrophicus population-specific biofilm formation and epibiotic growth on marine copepods

Vibrio spp. form a part of the microbiome of copepods-an abundant component of marine mesozooplankton. The biological mechanisms of the Vibrio-copepod association are largely unknown. In this study we compared biofilm formation of V. cyclitrophicus isolated from copepods (L-strains related to other particle-associated strains) and closely related strains originating from seawater (S-strains), and visualized and quantified their attachment and growth on copepods. The S- and L-strains formed similar biofilms in the presence of complete sea salts, suggesting previously unknown biofilm mechanisms in the S-strains. No biofilms formed if sodium chloride was present as the only salt but added calcium significantly enhanced biofilms in the L-strains. GFP-L-strain cells attached to live copepods at higher numbers than the S-strains, suggesting distinct mechanisms, potentially including calcium, support their colonization of copepods. The cells grew on live copepods after attachment, demonstrating that copepods sustain epibiotic V. cyclitrophicus growth in situ. The results demonstrate that in spite of their 99.1% average nucleotide identity, these V. cyclitrophicus strains have a differential capacity to colonize marine copepods. The introduced V. cyclitrophicus-A. tonsa model could be informative in future studies on Vibrio-copepod association.

Data-driven identification of major axes of functional variation in bacteria

The discovery of major axes of correlated functional variation among species and habitats has revealed the fundamental trade-offs structuring both functional and taxonomic diversity in eukaryotes such as plants. Whether such functional axes exist in the bacterial realm and whether they could explain bacterial taxonomic turnover among ecosystems remains unknown. Here, we use a data-driven approach to leverage global genomic and metagenomic datasets to reveal the existence of major axes of functional variation explaining both evolutionary differentiation within Bacteria and their ecological sorting across diverse habitats. We show that metagenomic variation among bacterial communities from various ecosystems is structured along a few axes of correlated functional pathways. Similar clusters of traits explained phylogenetic trait variation among >16,000 bacterial genomes, suggesting that functional turnover among bacterial communities from distinct habitats does not only result from the differential filtering of similar functions among communities, but also from phylogenetic correlations among these functions. Concordantly, functional pathways associated with trait clusters that were most important for defining functional turnover among bacterial communities were also those that had the highest phylogenetic signal in the bacterial genomic phylogeny. This study overall underlines the important role of evolutionary history in shaping contemporary distributions of bacteria across ecosystems.

Microdiversity sustains the distribution of rhizosphere-associated bacterial species from the root surface to the bulk soil region in maize crop fields

Over the years, the microbial community of maize (Zea mays) rhizosphere has been extensively studied; however, the role of microdiversity sustain rhizosphere-associated microbial species distribution from root surface to bulk soil in mature maize is still unclear. Although operational taxonomic units (OTUs) have been used to classify species, amplicon sequence variants (ASVs) have been shown to be effective in representing microdiversity within OTUs at a finer genetic scale. Therefore, the aim of this study was to examine the role of microdiversity in influencing the distribution of rhizosphere-associated microbial species across environmental gradients from root surface to bulk soil at the OTU and ASV levels. Here, the microbial community structures of bulk, loosely bound, and tightly bound soil samples from maize rhizosphere were examined at OTU and ASV levels. The results showed that OTU and ASV methods exhibited similar microbial community structures in rhizosphere. Additionally, different ecotypes with varying distributions and habitat preferences were observed within the same bacterial OTU at the ASV level, indicating a rich bacterial microdiversity. In contrast, the fungal community exhibited low microdiversity, with no significant relationship between fungal microdiversity and persistence and variability. Moreover, the ecotypes observed within the bacterial OTUs were found to be positively or negatively associated with environmental factors, such as soil organic carbon (SOC), NO3 --N, NH4 +-N contents, and pH. Overall, the results showed that the rich microdiversity could sustain the distribution of rhizosphere-associated bacterial species across environmental gradients from root surface to bulk soil. Further genetic analyses of rhizosphere-associated bacterial species could have considerable implications for potential mediation of microdiversity for sustainable crop production.

Metacommunity structure preserves genome diversity in the presence of gene-specific selective sweeps under moderate rates of horizontal gene transfer

The horizontal transfer of genes is fundamental for the eco-evolutionary dynamics of microbial communities, such as oceanic plankton, soil, and the human microbiome. In the case of an acquired beneficial gene, classic population genetics would predict a genome-wide selective sweep, whereby the genome spreads clonally within the community and together with the beneficial gene, removing genome diversity. Instead, several sources of metagenomic data show the existence of "gene-specific sweeps", whereby a beneficial gene spreads across a bacterial community, maintaining genome diversity. Several hypotheses have been proposed to explain this process, including the decreasing gene flow between ecologically distant populations, frequency-dependent selection from linked deleterious allelles, and very high rates of horizontal gene transfer. Here, we propose an additional possible scenario grounded in eco-evolutionary principles. Specifically, we show by a mathematical model and simulations that a metacommunity where species can occupy multiple patches, acting together with a realistic (moderate) HGT rate, helps maintain genome diversity. Assuming a scenario of patches dominated by single species, our model predicts that diversity only decreases moderately upon the arrival of a new beneficial gene, and that losses in diversity can be quickly restored. We explore the generic behaviour of diversity as a function of three key parameters, frequency of insertion of new beneficial genes, migration rates and horizontal transfer rates.Our results provides a testable explanation for how diversity can be maintained by gene-specific sweeps even in the absence of high horizontal gene transfer rates.

Multi-host infection and phylogenetically diverse lineages shape the recombination and gene pool dynamics of Staphylococcus aureus

BACKGROUND

Staphylococcus aureus can infect and adapt to multiple host species. However, our understanding of the genetic and evolutionary drivers of its generalist lifestyle remains inadequate. This is particularly important when considering local populations of S. aureus, where close physical proximity between bacterial lineages and between host species may facilitate frequent and repeated interactions between them. Here, we aim to elucidate the genomic differences between human- and animal-derived S. aureus from 437 isolates sampled from disease cases in the northeast region of the United States.

RESULTS

Multi-locus sequence typing revealed the existence of 75 previously recognized sequence types (ST). Our population genomic analyses revealed heterogeneity in the accessory genome content of three dominant S. aureus lineages (ST5, ST8, ST30). Genes related to antimicrobial resistance, virulence, and plasmid types were differentially distributed among isolates according to host (human versus non-human) and among the three major STs. Across the entire population, we identified a total of 1,912 recombination events that occurred in 765 genes. The frequency and impact of homologous recombination were comparable between human- and animal-derived isolates. Low-frequency STs were major donors of recombined DNA, regardless of the identity of their host. The most frequently recombined genes (clfB, aroA, sraP) function in host infection and virulence, which were also frequently shared between the rare lineages.

CONCLUSIONS

Taken together, these results show that frequent but variable patterns of recombination among co-circulating S. aureus lineages, including the low-frequency lineages, that traverse host barriers shape the structure of local gene pool and the reservoir of host-associated genetic variants. Our study provides important insights to the genetic and evolutionary factors that contribute to the ability of S. aureus to colonize and cause disease in multiple host species. Our study highlights the importance of continuous surveillance of S. aureus circulating in different ecological host niches and the need to systematically sample from them. These findings will inform development of effective measures to control S. aureus colonization, infection, and transmission across the One Health continuum.

Comparative genomics unveils extensive genomic variation between populations of Listeria species in natural and food-associated environments

Comprehending bacterial genomic variation linked to distinct environments can yield novel insights into mechanisms underlying differential adaptation and transmission of microbes across environments. Gaining such insights is particularly crucial for pathogens as it benefits public health surveillance. However, the understanding of bacterial genomic variation is limited by a scarcity of investigations in genomic variation coupled with different ecological contexts. To address this limitation, we focused on Listeria, an important bacterial genus for food safety that includes the human pathogen L. monocytogenes, and analyzed a large-scale genomic dataset collected by us from natural and food-associated environments across the United States. Through comparative genomics analyses on 449 isolates from the soil and 390 isolates from agricultural water and produce processing facilities representing L. monocytogenes, L. seeligeri, L. innocua, and L. welshimeri, we find that the genomic profiles strongly differ by environments within each species. This is supported by the environment-associated subclades and differential presence of plasmids, stress islands, and accessory genes involved in cell envelope biogenesis and carbohydrate transport and metabolism. Core genomes of Listeria species are also strongly associated with environments and can accurately predict isolation sources at the lineage level in L. monocytogenes using machine learning. We find that the large environment-associated genomic variation in Listeria appears to be jointly driven by soil property, climate, land use, and accompanying bacterial species, chiefly representing Actinobacteria and Proteobacteria. Collectively, our data suggest that populations of Listeria species have genetically adapted to different environments, which may limit their transmission from natural to food-associated environments.

Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population

Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.

Revisiting the Design of the Long-Term Evolution Experiment with Escherichia coli

The long-term evolution experiment (LTEE) with Escherichia coli began in 1988 and it continues to this day, with its 12 populations having recently reached 75,000 generations of evolution in a simple, well-controlled environment. The LTEE was designed to explore open-ended questions about the dynamics and repeatability of phenotypic and genetic evolution. Here I discuss various aspects of the LTEE's experimental design that have enabled its stability and success, including the choices of the culture regime, growth medium, ancestral strain, and statistical replication. I also discuss some of the challenges associated with a long-running project, such as handling procedural errors (e.g., cross-contamination) and managing the expanding collection of frozen samples. The simplicity of the experimental design and procedures have supported the long-term stability of the LTEE. That stability-along with the inherent creativity of the evolutionary process and the emergence of new genomic technologies-provides a platform that has allowed talented students and collaborators to pose questions, collect data, and make discoveries that go far beyond anything I could have imagined at the start of the LTEE.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Tracking contemporary microbial evolution in a changing ocean

Ocean microbes are fundamental for the functioning of the Earth system. Yet, our understanding of how they are reacting to global change in terms of evolution is limited. Microbes typically grow in large populations and reproduce quickly, which may allow them to rapidly adapt to environmental stressors compared to larger organisms. However, genetic evidence of contemporary evolution in wild microbes is scarce. We must begin coordinated efforts to establish new microbial time-series and explore novel tools, experiments, and data to fill this knowledge gap. The development of coordinated microbial 'genomic' observatories will provide the unprecedented opportunity to track contemporary microbial evolution in the ocean and explore the role of evolution in enabling wild microbes to respond to global change.

Differential carbon utilization enables co-existence of recently speciated Campylobacteraceae in the cow rumen epithelial microbiome

The activities of different microbes in the cow rumen have been shown to modulate the host's ability to utilize plant biomass, while the host-rumen interface has received little attention. As datasets collected worldwide have pointed to Campylobacteraceae as particularly abundant members of the rumen epithelial microbiome, we targeted this group in a subset of seven cows with meta- and isolate genome analysis. We show that the dominant Campylobacteraceae lineage has recently speciated into two populations that were structured by genome-wide selective sweeps followed by population-specific gene import and recombination. These processes led to differences in gene expression and enzyme domain composition that correspond to the ability to utilize acetate, the main carbon source for the host, at the cost of inhibition by propionate. This trade-off in competitive ability further manifests itself in differential dynamics of the two populations in vivo. By exploring population-level adaptations that otherwise remain cryptic in culture-independent analyses, our results highlight how recent evolutionary dynamics can shape key functional roles in the rumen microbiome.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Genome-Centric Dynamics Shape the Diversity of Oral Bacterial Populations

Two major viewpoints have been put forward for how microbial populations change, differing in whether adaptation is driven principally by gene-centric or genome-centric processes. Longitudinal sampling at microbially relevant timescales, i.e., days to weeks, is critical for distinguishing these mechanisms. Because of its significance for both microbial ecology and human health and its accessibility and high level of curation, we used the oral microbiota to study bacterial intrapopulation genome dynamics. Metagenomes were generated by shotgun sequencing of total community DNA from the healthy tongues of 17 volunteers at four to seven time points obtained over intervals of days to weeks. We obtained 390 high-quality metagenome-assembled genomes (MAGs) defining population genomes from 55 genera. The vast majority of genes in each MAG were tightly linked over the 2-week sampling window, indicating that the majority of the population's genomes were temporally stable at the MAG level. MAG-defined populations were composed of up to 5 strains, as determined by single-nucleotide-variant frequencies. Although most were stable over time, individual strains carrying over 100 distinct genes that rose from low abundance to dominance in a population over a period of days were detected. These results indicate a genome-wide as opposed to a gene-level process of population change. We infer that genome-wide selection of ecotypes is the dominant mode of adaptation in the oral populations over short timescales. IMPORTANCE The oral microbiome represents a microbial community of critical relevance to human health. Recent studies have documented the diversity and dynamics of different bacteria to reveal a rich, stable ecosystem characterized by strain-level dynamics. However, bacterial populations and their genomes are neither monolithic nor static; their genomes are constantly evolving to lose, gain, or alter their functional potential. To better understand how microbial genomes change in complex communities, we used culture-independent approaches to reconstruct the genomes (MAGs) for bacterial populations that approximated different species, in 17 healthy donors' mouths over a 2-week window. Our results underscored the importance of strain-level dynamics, which agrees with and expands on the conclusions of previous research. Altogether, these observations reveal patterns of genomic dynamics among strains of oral bacteria occurring over a matter of days.

Differences in genetic flux in invasive Streptococcus pneumoniae associated with bacteraemia and meningitis

Background

Genetic flux, a crucial process of pneumococcal evolution, is an essential aspect of bacterial physiology during human pathogenesis. However, the role of these genetic changes and the selective forces that drive them is not fully understood. Elucidating the underlying selective forces that determine the magnitude and direction (gene gain or loss) of gene transfer is important for better understanding the pathogenesis process, and may also highlight potential therapeutic and diagnostic targets.

Methods

Here, we leveraged data from high throughput genome sequencing and robust probabilistic models to discover the magnitude and likely direction of genetic flux events, but not the source, in 209 multi-lineage invasive pneumococcal genomes generated from blood (n = 147) and CSF (n = 62) isolates, associated with bacteremia and meningitis respectively. The Gain and Loss Mapping Engine (GLOOME) was used to infer gene gain and loss more accurately by taking into account differences in rates of gene gain and loss among gene families, as well as independent evolution within and across lineages.

Results

Our results show the likely extent and direction of gene fluctuations at different niche, during pneumococcal pathogenesis, highlighting that evolutionary dynamics are important for tissue-specific host invasion and survival.

Conclusion

These findings improve insights on evolutionary dynamics during invasive pneumococcal disease, and highlight potential diagnostic and therapeutic targets.

Gene flow and introgression are pervasive forces shaping the evolution of bacterial species

BACKGROUND

Although originally thought to evolve clonally, studies have revealed that most bacteria exchange DNA. However, it remains unclear to what extent gene flow shapes the evolution of bacterial genomes and maintains the cohesion of species.

RESULTS

Here, we analyze the patterns of gene flow within and between >2600 bacterial species. Our results show that fewer than 10% of bacterial species are truly clonal, indicating that purely asexual species are rare in nature. We further demonstrate that the taxonomic criterion of ~95% genome sequence identity routinely used to define bacterial species does not accurately represent a level of divergence that imposes an effective barrier to gene flow across bacterial species. Interruption of gene flow can occur at various sequence identities across lineages, generally from 90 to 98% genome identity. This likely explains why a ~95% genome sequence identity threshold has empirically been judged as a good approximation to define bacterial species. Our results support a universal mechanism where the availability of identical genomic DNA segments required to initiate homologous recombination is the primary determinant of gene flow and species boundaries in bacteria. We show that these barriers of gene flow remain porous since many distinct species maintain some level of gene flow, similar to introgression in sexual organisms.

CONCLUSIONS

Overall, bacterial evolution and speciation are likely shaped by similar forces driving the evolution of sexual organisms. Our findings support a model where the interruption of gene flow-although not necessarily the initial cause of speciation-leads to the establishment of permanent and irreversible species borders.

Response of soil viral communities to land use changes

Soil viruses remain understudied when compared to virus found in aquatic ecosystems. Here, we investigate the ecological patterns of soil viral communities across various land use types encompassing forest, agricultural, and urban soil in Xiamen, China. We recovered 59,626 viral operational taxonomic units (vOTUs) via size-fractioned viromic approach with additional mitomycin C treatment to induce virus release from bacterial fraction. Our results show that viral communities are significantly different amongst the land use types considered. A microdiversity analysis indicates that selection act on soil vOTUs, resulting in disparities between land use associated viral communities. Soil pH is one of the major determinants of viral community structure, associated with changes of in-silico predicted host compositions of soil vOTUs. Habitat disturbance and variation of soil moisture potentially contribute to the dynamics of putative lysogenic vOTUs. These findings provide mechanistic understandings of the ecology and evolution of soil viral communities in changing environments.

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.

Incipient Sympatric Speciation and Evolution of Soil Bacteria Revealed by Metagenomic and Structured Non-Coding RNAs Analysis

Simple Summary

The microevolutionary dynamics of soil bacteria under microclimatic differences are largely unexplored in contrast to our improving knowledge of their vast diversity. In this study, we performed a comparative metagenomic analysis of two sharply divergent rocks and soil types at the Evolution Plateau (EP) in eastern Upper Galilee, Israel. We have identified the significant differences in bacterial taxonomic diversity, functions, and patterns of RNA-based gene regulation between the bacteria from two different soil types. Furthermore, we have identified several species with a significant genetic divergence of the same species between the two soil types, highlighting the soil bacteria’s incipient sympatric speciation.

Abstract

Soil bacteria respond rapidly to changes in new environmental conditions. For adaptation to the new environment, they could mutate their genome, which impacts the alternation of the functional and regulatory landscape. Sometimes, these genetic and ecological changes may drive the bacterial evolution and sympatric speciation. Although sympatric speciation has been controversial since Darwin suggested it in 1859, there are several strong theoretical or empirical evidences to support it. Sympatric speciation associated with soil bacteria remains largely unexplored. Here, we provide potential evidence of sympatric speciation of soil bacteria by comparison of metagenomics from two sharply contrasting abutting divergence rock and soil types (Senonian chalk and its rendzina soil, and abutting Pleistocene basalt rock and basalt soil). We identified several bacterial species with significant genetic differences in the same species between the two soil types and ecologies. We show that the bacterial community composition has significantly diverged between the two soils; correspondingly, their functions were differentiated in order to adapt to the local ecological stresses. The ecologies, such as water availability and pH value, shaped the adaptation and speciation of soil bacteria revealed by the clear-cut genetic divergence. Furthermore, by a novel analysis scheme of riboswitches, we highlight significant differences in structured non-coding RNAs between the soil bacteria from two divergence soil types, which could be an important driver for functional adaptation. Our study provides new insight into the evolutionary divergence and incipient sympatric speciation of soil bacteria under microclimatic ecological differences.

Coastal Transient Niches Shape the Microdiversity Pattern of a Bacterioplankton Population with Reduced Genomes

Globally dominant marine bacterioplankton lineages are often limited in metabolic versatility, owing to their extensive genome reductions, and thus cannot take advantage of transient nutrient patches. It is therefore perplexing how the nutrient-poor bulk seawater sustains the pelagic streamlined lineages, each containing numerous populations. Here, we sequenced the genomes of 33 isolates of the recently discovered CHUG lineage (~2.6 Mbp), which have some of the smallest genomes in the globally abundant Roseobacter group (commonly over 4 Mbp). These genome-reduced bacteria were isolated from a transient habitat: seawater surrounding the brown alga, Sargassum hemiphyllum. Population genomic analyses showed that: (i) these isolates, despite sharing identical 16S rRNA genes, were differentiated into several genetically isolated populations through successive speciation events; (ii) only the first speciation event led to the genetic separation of both core and accessory genomes; and (iii) populations resulting from this event are differentiated at many loci involved in carbon utilization and oxygen respiration, corroborated by BiOLOG phenotype microarray assays and oxygen uptake kinetics experiments, respectively. These differentiated traits match well with the dynamic nature of the macroalgal seawater, in which the quantity and quality of carbon sources and the concentration of oxygen likely vary spatially and temporally, though other habitats, like fresh organic aggregates, cannot be ruled out. Our study implies that transient habitats in the overall nutrient-poor ocean can shape the microdiversity and population structure of genome-reduced bacterioplankton lineages.

Differentiated Evolutionary Strategies of Genetic Diversification in Atlantic and Pacific Thaumarchaeal Populations

Some marine microbes are seemingly “ubiquitous,” thriving across a wide range of environmental conditions. While the increased depth in metagenomic sequencing has led to a growing body of research on within-population heterogeneity in environmental microbial populations, there have been fewer systematic comparisons and characterizations of population-level genetic diversity over broader expanses of time and space. Here, we investigated the factors that govern the diversification of ubiquitous microbial taxa found within and between ocean basins. Specifically, we use mapped metagenomic paired reads to examine the genetic diversity of ammonia-oxidizing archaeal (“Candidatus Nitrosopelagicus brevis”) populations in the Pacific (Hawaii Ocean Time-series [HOT]) and Atlantic (Bermuda Atlantic Time Series [BATS]) Oceans sampled over 2 years. We observed higher nucleotide diversity in “Ca. N. brevis” at HOT, driven by a higher rate of homologous recombination. In contrast, “Ca. N. brevis” at BATS featured a more open pangenome with a larger set of genes that were specific to BATS, suggesting a history of dynamic gene gain and loss events. Furthermore, we identified highly differentiated genes that were regulatory in function, some of which exhibited evidence of recent selective sweeps. These findings indicate that different modes of genetic diversification likely incur specific adaptive advantages depending on the selective pressures that they are under. Within-population diversity generated by the environment-specific strategies of genetic diversification is likely key to the ecological success of “Ca. N. brevis.”

IMPORTANCE Ammonia-oxidizing archaea (AOA) are one of the most abundant chemolithoautotrophic microbes in the marine water column and are major contributors to global carbon and nitrogen cycling. Despite their ecological importance and geographical pervasiveness, there have been limited systematic comparisons and characterizations of their population-level genetic diversity over time and space. Here, we use metagenomic time series from two ocean observatories to address the fundamental questions of how abiotic and biotic factors shape the population-level genetic diversity and how natural microbial populations adapt across diverse habitats. We show that the marine AOA “Candidatus Nitrosopelagicus brevis” in different ocean basins exhibits distinct modes of genetic diversification in response to their selective regimes shaped by nutrient availability and patterns of environmental fluctuations. Our findings specific to “Ca. N. brevis” have broader implications, particularly in understanding the population-level responses to the changing climate and predicting its impact on biogeochemical cycles.

Diversity and ecological footprint of Global Ocean RNA viruses

DNA viruses are increasingly recognized as influencing marine microbes and microbe-mediated biogeochemical cycling. However, little is known about global marine RNA virus diversity, ecology, and ecosystem roles. In this study, we uncover patterns and predictors of marine RNA virus community- and "species"-level diversity and contextualize their ecological impacts from pole to pole. Our analyses revealed four ecological zones, latitudinal and depth diversity patterns, and environmental correlates for RNA viruses. Our findings only partially parallel those of cosampled plankton and show unexpectedly high polar ecological interactions. The influence of RNA viruses on ecosystems appears to be large, as predicted hosts are ecologically important. Moreover, the occurrence of auxiliary metabolic genes indicates that RNA viruses cause reprogramming of diverse host metabolisms, including photosynthesis and carbon cycling, and that RNA virus abundances predict ocean carbon export.

A network approach to decipher the dynamics of Lysobacteraceae plasmid gene sharing

Plasmids provide an efficient vehicle for gene sharing among bacterial populations, playing a key role in bacterial evolution. Network approaches are particularly suitable to represent multipartite relationships and are useful tools to characterize plasmid-mediated gene sharing events. The Lysobacteraceae bacterial family gathers plant commensal, plant pathogenic and opportunistic human pathogens for which plasmid mediated adaptation was reported. We searched for homologues of plasmid gene sequences from this family in all the diversity of available bacterial genome sequences and built a network of plasmid gene sharing from the results. While plasmid genes are openly shared between the bacteria of the Lysobacteraceae family, taxonomy strongly defined the boundaries of these exchanges, that only barely reached other families. Most inferred plasmid gene sharing events involved a few genes only, and evidence of full plasmid transfers were restricted to taxonomically close taxon. We detected multiple plasmid-chromosome gene transfers, among which the otherwise known sharing of a heavy metal resistance transposon. In the network, bacterial lifestyles shaped sub-structures of isolates colonizing specific ecological niches and harboring specific types of resistance genes. Genes associated to pathogenicity or antibiotic and metal resistance were among those that most importantly structured the network, highlighting the imprints of human-mediated selective pressure on pathogenic populations. A massive sequencing effort on environmental Lysobacteraceae is therefore required to refine our understanding on how this reservoir fuels the emergence and the spread of genes amongst this family and its potential impact on plant, animal and human health.

Horizontal gene transfer and adaptive evolution in bacteria

Horizontal gene transfer (HGT) is arguably the most conspicuous feature of bacterial evolution. Evidence for HGT is found in most bacterial genomes. Although HGT can considerably alter bacterial genomes, not all transfer events may be biologically significant and may instead represent the outcome of an incessant evolutionary process that only occasionally has a beneficial purpose. When adaptive transfers occur, HGT and positive selection may result in specific, detectable signatures in genomes, such as gene-specific sweeps or increased transfer rates for genes that are ecologically relevant. In this Review, we first discuss the various mechanisms whereby HGT occurs, how the genetic signatures shape patterns of genomic variation and the distinct bioinformatic algorithms developed to detect these patterns. We then discuss the evolutionary theory behind HGT and positive selection in bacteria, and discuss the approaches developed over the past decade to detect transferred DNA that may be involved in adaptation to new environments.

A comparative whole-genome approach identifies bacterial traits for marine microbial interactions

Microbial interactions shape the structure and function of microbial communities with profound consequences for biogeochemical cycles and ecosystem health. Yet, most interaction mechanisms are studied only in model systems and their prevalence is unknown. To systematically explore the functional and interaction potential of sequenced marine bacteria, we developed a trait-based approach, and applied it to 473 complete genomes (248 genera), representing a substantial fraction of marine microbial communities. We identified genome functional clusters (GFCs) which group bacterial taxa with common ecology and life history. Most GFCs revealed unique combinations of interaction traits, including the production of siderophores (10% of genomes), phytohormones (3-8%) and different B vitamins (57-70%). Specific GFCs, comprising Alpha- and Gammaproteobacteria, displayed more interaction traits than expected by chance, and are thus predicted to preferentially interact synergistically and/or antagonistically with bacteria and phytoplankton. Linked trait clusters (LTCs) identify traits that may have evolved to act together (e.g., secretion systems, nitrogen metabolism regulation and B vitamin transporters), providing testable hypotheses for complex mechanisms of microbial interactions. Our approach translates multidimensional genomic information into an atlas of marine bacteria and their putative functions, relevant for understanding the fundamental rules that govern community assembly and dynamics.