A Systematics for Discovering the Fundamental Units of Bacterial Diversity

Emergence of potentially disinfection-resistant, naturalized Escherichia coli populations across food- and water-associated engineered environments

The Escherichia coli species is comprised of several 'ecotypes' inhabiting a wide range of host and natural environmental niches. Recent studies have suggested that novel naturalized ecotypes have emerged across wastewater treatment plants and meat processing facilities. Phylogenetic and multilocus sequence typing analyses clustered naturalized wastewater and meat plant E. coli strains into two main monophyletic clusters corresponding to the ST635 and ST399 sequence types, with several serotypes identified by serotyping, potentially representing distinct lineages that have naturalized across wastewater treatment plants and meat processing facilities. This evidence, taken alongside ecotype prediction analyses that distinguished the naturalized strains from their host-associated counterparts, suggests these strains may collectively represent a novel ecotype that has recently emerged across food- and water-associated engineered environments. Interestingly, pan-genomic analyses revealed that the naturalized strains exhibited an abundance of biofilm formation, defense, and disinfection-related stress resistance genes, but lacked various virulence and colonization genes, indicating that their naturalization has come at the cost of fitness in the original host 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.

Microevolutionary patterns in ecotypes of the symbiotic cyanobacterium UCYN-A revealed from a Northwest Atlantic coastal time series

UCYN-A is a globally important nitrogen-fixing symbiotic microbe often found in colder regions and coastal areas where nitrogen fixation has been overlooked. We present a 3-year coastal Northwest Atlantic time series of UCYN-A by integrating oceanographic data with weekly nifH and16S rRNA gene sequencing and quantitative PCR assays for UCYN-A ecotypes. High UCYN-A relative abundances dominated by A1 to A4 ecotypes reoccurred annually in the coastal Northwest Atlantic. Although UCYN-A was detected every summer/fall, the ability to observe separate ecotypes may be highly dependent on sampling time given intense interannual and weekly variability of ecotype-specific occurrences. Additionally, much of UCYN-A's rarer diversity was populated by short-lived neutral mutational variants, therefore providing insight into UCYN-A's microevolutionary patterns. For instance, rare ASVs exhibited community composition restructuring annually, while also sharing a common connection to a dominant ASV within each ecotype. Our study provides additional perspectives for interpreting UCYN-A intraspecific diversity and underscores the need for high-resolution datasets when deciphering spatiotemporal ecologies within UCYN-A.

Prevalence of Viral Frequency-Dependent Infection in Coastal Marine Prokaryotes Revealed Using Monthly Time Series Virome Analysis

Viruses infecting marine prokaryotes have a large impact on the diversity and dynamics of their hosts. Model systems suggest that viral infection is frequency dependent and constrained by the virus-host encounter rate. However, it is unclear whether frequency-dependent infection is pervasive among the abundant prokaryotic populations with different temporal dynamics. To address this question, we performed a comparison of prokaryotic and viral communities using 16S rRNA amplicon and virome sequencing based on samples collected monthly for 2 years at a Japanese coastal site, Osaka Bay. Concurrent seasonal shifts observed in prokaryotic and viral community dynamics indicated that the abundance of viruses correlated with that of their predicted host phyla (or classes). Cooccurrence network analysis between abundant prokaryotes and viruses revealed 6,423 cooccurring pairs, suggesting a tight coupling of host and viral abundances and their "one-to-many" correspondence. Although stable dominant species, such as SAR11, showed few cooccurring viruses, a fast succession of their viruses suggests that viruses infecting these populations changed continuously. Our results suggest that frequency-dependent viral infection prevails in coastal marine prokaryotes regardless of host taxa and temporal dynamics. IMPORTANCE There is little room for doubt that viral infection is prevalent among abundant marine prokaryotes regardless of their taxa or growth strategy. However, comprehensive evaluations of viral infections in natural prokaryotic communities are still technically difficult. In this study, we examined viral infection in abundant prokaryotes by monitoring the monthly dynamics of prokaryotic and viral communities at a eutrophic coastal site, Osaka Bay. We compared the community dynamics of viruses with those of their putative hosts based on genome-based in silico host prediction. We observed frequent cooccurrence among the predicted virus-host pairs, suggesting that viral infection is prevalent in abundant prokaryotes regardless of their taxa or temporal dynamics. This likely indicates that frequent lysis of the abundant prokaryotes via viral infection has a considerable contribution to the biogeochemical cycling and maintenance of prokaryotic community diversity.

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.

Genomic Insights into the Radiation-Resistant Capability of Sphingomonas qomolangmaensis S5-59T and Sphingomonas glaciei S8-45T, Two Novel Bacteria from the North Slope of Mount Everest

Mount Everest provides natural advantages to finding radiation-resistant extremophiles that are functionally mechanistic and possess commercial significance. (1) Background: Two bacterial strains, designated S5-59T and S8-45T, were isolated from moraine samples collected from the north slope of Mount Everest at altitudes of 5700m and 5100m above sea level. (2) Methods: The present study investigated the polyphasic features and genomic characteristics of S5-59^T and S8-45^T. (3) Results: The major fatty acids and the predominant respiratory menaquinone of S5-59^T and S8-45^T were summed as feature 3 (comprising C16:1 ω6c and/or C16:1 ω7c) and ubiquinone-10 (Q-10). Phylogenetic analyses based on 16S rRNA sequences and average nucleotide identity values among these two strains and their reference type strains were below the species demarcation thresholds of 98.65% and 95%. Strains S5-59^T and S8-45^T harbored great radiation resistance. The genomic analyses showed that DNA damage repair genes, such as mutL, mutS, radA, radC, recF, recN, etc., were present in the S5-59^T and S8-45^T strains. Additionally, strain S5-59^T possessed more genes related to DNA protection proteins. The pan-genome analysis and horizontal gene transfers revealed that strains of Sphingomonas had a consistently homologous genetic evolutionary radiation resistance. Moreover, enzymatic antioxidative proteins also served critical roles in converting ROS into harmless molecules that resulted in resistance to radiation. Further, pigments and carotenoids such as zeaxanthin and alkylresorcinols of the non-enzymatic antioxidative system were also predicted to protect them from radiation. (4) Conclusions: Type strains S5-59^T (=JCM 35564T =GDMCC 1.3193T) and S8-45^T (=JCM 34749T =GDMCC 1.2715T) represent two novel species of the genus Sphingomonas with the proposed name Sphingomonas qomolangmaensis sp. nov. and Sphingomonas glaciei sp. nov. The type strains, S5-59^T and S8-45^T, were assessed in a deeply genomic study of their radiation-resistant mechanisms and this thus resulted in a further understanding of their greater potential application for the development of anti-radiation protective drugs.

Exploring the mobilome and resistome of Enterococcus faecium in a One Health context across two continents

Enterococcus faecium is a ubiquitous opportunistic pathogen that is exhibiting increasing levels of antimicrobial resistance (AMR). Many of the genes that confer resistance and pathogenic functions are localized on mobile genetic elements (MGEs), which facilitate their transfer between lineages. Here, features including resistance determinants, virulence factors and MGEs were profiled in a set of 1273 E. faecium genomes from two disparate geographic locations (in the UK and Canada) from a range of agricultural, clinical and associated habitats. Neither lineages of E. faecium, type A and B, nor MGEs are constrained by geographic proximity, but our results show evidence of a strong association of many profiled genes and MGEs with habitat. Many features were associated with a group of clinical and municipal wastewater genomes that are likely forming a new human-associated ecotype within type A. The evolutionary dynamics of E. faecium make it a highly versatile emerging pathogen, and its ability to acquire, transmit and lose features presents a high risk for the emergence of new pathogenic variants and novel resistance combinations. This study provides a workflow for MGE-centric surveillance of AMR in Enterococcus that can be adapted to other pathogens.

Genomic insights revealed the environmental adaptability of Planococcus halotolerans Y50 isolated from petroleum-contaminated soil on the Qinghai-Tibet Plateau

The Tibetan Plateau niche provides unprecedented opportunities to find microbes that are functional and commercial significance. The present study investigated the physiological and genomic characteristics of Planococcus halotolerans Y50 that was isolated from a petroleum-contaminated soil sample from the Qinghai-Tibet Plateau, and it displayed psychrotolerant, antiradiation, and oil-degraded characteristics. Whole genome sequencing indicated that strain Y50 has a 3.52 Mb genome and 44.7% G + C content, and it possesses 3377 CDSs. The presence of a wide range of UV damage repair genes uvrX and uvsE, DNA repair genes radA and recN, superoxide dismutase, peroxiredoxin and dioxygenase genes provided the genomic basis for the adaptation of the plateau environment polluted by petroleum. Related experiments also verified that the Y50 strain could degrade n-alkanes from C₁₁-C₂₃, and approximately 30% of the total petroleum at 25 °C within 7 days. Meanwhile, strain Y50 could withstand 5 × 10³ J/m² UVC and 10 KGy gamma ray radiation, and it had strong antioxidant and high radical scavengers for superoxide anion, hydroxyl radical and DPPH. In addition, pan-genome analysis and horizontal gene transfers revealed that strains with different niches have obtained various genes through horizontal gene transfer in the process of evolution, and the more similar their geographical locations, the more similar their members are genetically and ecologically. In conclusion, P. halotolerans Y50 possesses high potential of applications in the bioremediation of alpine hydrocarbons contaminated environment.

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.

Genetic and Ecological Diversity of Escherichia coli and Cryptic Escherichia Clades in Subtropical Aquatic Environments

Escherichia coli not only inhabit the large intestines of human and warm-blooded animals but could also persist in the external environment. However, current knowledge was largely based on host-associated strains. Moreover, cryptic Escherichia clades that were often misidentified as E. coli by conventional diagnostic methods were discovered. Failure to distinguish them from E. coli sensu stricto could lead to inaccurate conclusions about the population genetics of E. coli. Based on seven housekeeping genes, we determine the genetic and ecological diversity of E. coli and cryptic clades as they occupy aquatic habitats with different characteristics and human impact levels in subtropical Hong Kong. Contrary to previous reports, clade II was the most abundant cryptic lineage co-isolated with E. coli, being especially abundant in relatively pristine subtropical aquatic environments. The phylogenetically distinct cryptic clades and E. coli showed limited recombination and significant genetic divergence. Analyses indicated that these clade II strains were ecologically differentiated from typical E. coli; some may even represent novel environmental Escherichia clades that were closely related to the original clade II strains of fecal origins. E. coli of diverse origins exhibited clonality amidst divergent genotypes STs, echoing other studies in that recombination in housekeeping genes was insufficient to disrupt phylogenetic signals of the largely clonal E. coli. Notably, environmental E. coli were less diverse than fecal isolates despite contributing many new alleles and STs. Finally, we demonstrated that human activities influenced the distribution of E. coli and clade II in a small aquatic continuum. Moving from relatively pristine sites toward areas with higher human disturbance, the abundance of clade II isolates and new E. coli genotypes reduces, while E. coli bearing class I integrons and belonging to CCs of public health concern accumulates. Altogether, this work revealed the new genetic diversity of E. coli and cryptic clades embedded in selected subtropical aquatic habitats, especially relatively pristine sites, which will aid a more thorough understanding of the extent of their genetic and functional variations in relation to diverse habitats with varied conditions.

Fine-scale genetic diversity and putative ecotypes of oxymonad protists coinhabiting the hindgut of Reticulitermes speratus

The hindgut of lower termites is generally coinhabited by multiple morphologically identifiable protist species. However, it is unclear how many protist species truly coexist in this miniaturized environment, and moreover, it is difficult to define the fundamental unit of protist diversity. Species delineation of termite gut protists has therefore been guided without a theory-based concept of species. Here, we focused on the hindgut of the termite Reticulitermes speratus, where 10 or 11 morphologically distinct oxymonad cell types, that is, morphospecies, coexist. We elucidated the phylogenetic structure of all co-occurring oxymonads and addressed whether their diversity can be explained by the "ecotype" hypothesis. Oxymonad-specific 18S rRNA gene amplicon sequencing analyses of whole-gut samples, combined with single-cell 18S rRNA sequencing of the oxymonad morphospecies, identified 210 one-nucleotide-level variants. The phylogenetic analysis of these variants revealed the presence of microdiverse clusters typically within 1% sequence divergence. Each known oxymonad morphospecies comprised one to several monophyletic or paraphyletic microdiverse clusters. Using these sequence data sets, we conducted computational simulation to predict the rates of ecotype formation and periodic selection, and to demarcate putative ecotypes. Our simulations suggested that the oxymonad genetic divergence is constrained primarily by strong selection, in spite of limited population size and possible bottlenecks during intergenerational transmission. A total of 33 oxymonad ecotypes were predicted, and most of the putative ecotypes were consistently detected among different colonies and host individuals. These findings provide a possible theoretical basis for species diversity and underlying mechanisms of coexistence of termite gut protists.

Genotyping and multivariate regression trees reveal ecological diversification within the Microcystis aeruginosa complex along a wide environmental gradient

Addressing the ecological and evolutionary processes underlying biodiversity patterns is essential to identify the mechanisms shaping community structure and function. In bacteria, the formation of new ecologically distinct populations (ecotypes) is proposed as one of the main drivers of diversification. New ecotypes arise when mutations in key functional genes or acquisition of new metabolic pathways by horizontal gene transfer allow the population to exploit new resources, permitting their coexistence with the parental population. We previously reported the presence of microcystin-producing organisms of the Microcystis aeruginosa complex (toxic MAC) through an 800 km environmental gradient ranging from freshwater to estuarine-marine waters in South America. We hypothesize that the success of toxic MAC in such a gradient is due to the existence of very closely related populations that are ecologically distinct (ecotypes), each specialized to a specific arrangement of environmental variables. Here, we analyzed toxic MAC genetic diversity through qPCR and high-resolution melting analysis (HRMA) of a functional gene (mcyJ, microcystin synthetase cluster). We explored the variability of the mcyJ gene along the environmental gradient by multivariate classification and regression trees (mCART). Six groups of mcyJ genotypes were distinguished and associated with different combinations of water temperature, conductivity and turbidity. We propose that each mcyJ variant associated to a defined environmental condition is an ecotype (or species) whose relative abundances vary according to their fitness in the local environment. This mechanism would explain the success of toxic MAC in such a wide array of environmental conditions. Importance Organisms of the Microcystis aeruginosa Complex form harmful algal blooms (HABs) in nutrient-rich water bodies worldwide. MAC HABs are difficult to manage owing to the production of potent toxins (microcystins) that resist water treatment. Besides, the role of microcystins in the ecology of MAC organisms is still elusive, meaning that the environmental conditions driving the toxicity of the bloom are not clear. Furthermore, the lack of coherence between morphology-based and genomic-based species classification makes it difficult to draw sound conclusions about when and where each member species of the MAC will dominate the bloom. Here, we propose that the diversification process and success of toxic MAC in a wide range of waterbodies involves the generation of ecotypes, each specialized in a particular niche, whose relative abundance varies according to its fitness in the local environment. This knowledge can improve the generation of accurate prediction models of MAC growth and toxicity, helping to prevent human and animal intoxication.

Ecological and Evolutionary responses to Antibiotic Treatment in the Human Gut Microbiota

The potential for antibiotics to affect the ecology and evolution of the human gut microbiota is well recognised and has wide-ranging implications for host health. Here, we review the findings of key studies that surveyed the human gut microbiota during antibiotic treatment. We find several broad patterns including the loss of diversity, disturbance of community composition, suppression of bacteria in the Actinobacteria phylum, amplification of bacteria in the Bacteroidetes phylum, and promotion of antibiotic resistance. Such changes to the microbiota were often, but not always, recovered following the end of treatment. However, many studies reported unique and/or contradictory results, which highlights our inability to meaningfully predict or explain the effects of antibiotic treatment on the human gut microbiome. This problem arises from variation between existing studies in three major categories: differences in dose, class and combinations of antibiotic treatments used; differences in demographics, lifestyles, and locations of subjects; and differences in measurements, analyses and reporting styles used by researchers. To overcome this, we suggest two integrated approaches: (i) a top-down approach focused on building predictive models through large sample sizes, deep metagenomic sequencing, and effective collaboration; and (ii) a bottom-up reductionist approach focused on testing hypotheses using model systems.

'Species' without species

Biological science uses multiple species concepts. Order can be brought to this diversity if we recognize two key features. First, any given species concept is likely to have a patchwork structure, generated by repeated application of the concept to new domains. We illustrate this by showing how two species concepts (biological and ecological) have been modified from their initial eukaryotic applications to apply to prokaryotes. Second, both within and between patches, distinct species concepts may interact and hybridize. We thus defend a semantic picture of the species concept as a collection of interacting patchwork structures. Thus, although not all uses of the term pick out the same kind of unit in nature, the diversity of uses reflects something more than mere polysemy. We suggest that the emphasis on the use of species to pick out natural units is itself problematic, because that is not the term's sole function. In particular, species concepts are used to manage inquiry into processes of speciation, even when these processes do not produce clearly delimited species.

Negative frequency-dependent selection and asymmetrical transformation stabilise multi-strain bacterial population structures

Streptococcus pneumoniae can be divided into many strains, each a distinct set of isolates sharing similar core and accessory genomes, which co-circulate within the same hosts. Previous analyses suggested the short-term vaccine-associated dynamics of S. pneumoniae strains may be mediated through multi-locus negative frequency-dependent selection (NFDS), which maintains accessory loci at equilibrium frequencies. Long-term simulations demonstrated NFDS stabilised clonally-evolving multi-strain populations through preventing the loss of variation through drift, based on polymorphism frequencies, pairwise genetic distances and phylogenies. However, allowing symmetrical recombination between isolates evolving under multi-locus NFDS generated unstructured populations of diverse genotypes. Replication of the observed data improved when multi-locus NFDS was combined with recombination that was instead asymmetrical, favouring deletion of accessory loci over insertion. This combination separated populations into strains through outbreeding depression, resulting from recombinants with reduced accessory genomes having lower fitness than their parental genotypes. Although simplistic modelling of recombination likely limited these simulations' ability to maintain some properties of genomic data as accurately as those lacking recombination, the combination of asymmetrical recombination and multi-locus NFDS could restore multi-strain population structures from randomised initial populations. As many bacteria inhibit insertions into their chromosomes, this combination may commonly underlie the co-existence of strains within a niche.

Long-run bacteria-phage coexistence dynamics under natural habitat conditions in an environmental biotechnology system

Bacterial viruses are widespread and abundant across natural and engineered habitats. They influence ecosystem functioning through interactions with their hosts. Laboratory studies of phage-host pairs have advanced our understanding of phenotypic and genetic diversification in bacteria and phages. However, the dynamics of phage-host interactions have been seldom recorded in complex natural environments. We conducted an observational metagenomic study of the dynamics of interaction between Gordonia and their phages using a three-year data series of samples collected from a full-scale wastewater treatment plant. The aim was to obtain a comprehensive picture of the coevolution dynamics in naturally evolving populations at relatively high time resolution. Coevolution was followed by monitoring changes over time in the CRISPR loci of Gordonia metagenome-assembled genome, and reciprocal changes in the viral genome. Genome-wide analysis indicated low strain variability of Gordonia, and almost clonal conservation of the trailer end of the CRISPR loci. Incorporation of newer spacers gave rise to multiple coexisting bacterial populations. The host population carrying a shorter CRISPR locus that contain only ancestral spacers, which has not acquired newer spacers against the coexisting phages, accounted for more than half of the total host abundance in the majority of samples. Phages genome co-evolved by introducing directional changes, with no preference for mutations within the protospacer and PAM regions. Metagenomic reconstruction of time-resolved variants of host and viral genomes revealed how the complexity at the population level has important consequences for bacteria-phage coexistence.

A fine-scale map of genome-wide recombination in divergent Escherichia coli population

Recombination is one of the most important molecular mechanisms of prokaryotic genome evolution, but its exact roles are still in debate. Here we try to infer genome-wide recombination within a species, utilizing a dataset of 149 complete genomes of Escherichia coli from diverse animal hosts and geographic origins, including 45 in-house sequenced with the single-molecular real-time platform. Two major clades identified based on physiological, clinical and ecological characteristics form distinct genetic lineages based on scarcity of interclade gene exchanges. By defining gene-based syntenies for genomic segments within and between the two clades, we build a fine-scale recombination map for this representative global E. coli population. The map suggests extensive within-clade recombination that often breaks physical linkages among individual genes but seldom interrupts the structure of genome organizational frameworks as well as primary metabolic portfolios supported by the framework integrity, possibly due to strong natural selection for both physiological compatibility and ecological fitness. In contrast, the between-clade recombination declines drastically when phylogenetic distance increases to the extent where a 10-fold reduction can be observed, establishing a firm genetic barrier between clades. Our empirical data suggest a critical role for such recombination events in the early stage of speciation where recombination rate is associated with phylogenetic distance in addition to sequence and gene variations. The extensive intraclade recombination binds sister strains into a quasisexual group and optimizes genes or alleles to streamline physiological activities, whereas the sharply declined interclade recombination split the population into clades adaptive to divergent ecological niches.

Genetic Variation and Preliminary Indications of Divergent Niche Adaptation in Cryptic Clade II of Escherichia

The evolution, habitat, and lifestyle of the cryptic clade II of Escherichia, which were first recovered at low frequency from non-human hosts and later from external environments, were poorly understood. Here, the genomes of selected strains were analyzed for preliminary indications of ecological differentiation within their population. We adopted the delta bitscore metrics to detect functional divergence of their orthologous genes and trained a random forest classifier to differentiate the genomes according to habitats (gastrointestinal vs external environment). Model was built with inclusion of other Escherichia genomes previously demonstrated to have exhibited genomic traits of adaptation to one of the habitats. Overall, gene degradation was more prominent in the gastrointestinal strains. The trained model correctly classified the genomes, identifying a set of predictor genes that were informative of habitat association. Functional divergence in many of these genes were reflective of ecological divergence. Accuracy of the trained model was confirmed by its correct prediction of the habitats of an independent set of strains with known habitat association. In summary, the cryptic clade II of Escherichia displayed genomic signatures that are consistent with divergent adaptation to gastrointestinal and external environments.

The evolution of bacterial pathogens in the Anthropocene

Humankind has become a primary driver of global environmental and climate change. The extent of planetary change is such that it has been proposed to classify the current geological age as the 'Anthropocene'. Anthropogenic environmental degradation presents numerous threats to human health and wellbeing, including an increased risk of infectious disease. This review focuses on how processes such as pollution, climate change and human-mediated dispersal could affect the evolution of bacterial pathogens. Effects of environmental change on the 'big five' of evolution: mutation rate, recombination (horizontal gene transfer), migration, selection and drift are discussed. Microplastic pollution is used as a case study to highlight the combined effects of some of these processes on the evolutionary diversification of human pathogens. Although the evidence is still incomplete, a picture is emerging that environmental pathogens could evolve at increased rates in the Anthropocene, with potential consequences for human infection.

Towards Genomic Criteria for Delineating Fungal Species

The discussion of fungal species delineation has yet to reach a consensus, despite the advancements in technology, which helped modernise traditional approaches. In particular, the phylogenetic species concept was one of the tools that has been used with considerable success across the fungal kingdom. The fast rise of fungal genomics provides an unprecedented opportunity to expand measuring the relatedness of fungal strains to the level of whole genomes. However, the use of genomic information in taxonomy has only just begun, and few methodological guidelines have been suggested so far. Here, a simple approach of computationally measuring genomic distances and their use as a standard for species delineation is investigated. A fixed threshold genomic distance calculated by the quick and easy-to-use tools Mash and Dashing proved to be an unexpectedly widely applicable and robust criterion for determining whether two genomes belong to the same or to different species. The accuracy of species delineation in an uncurated dataset of GenBank fungal genomes was close to 90%—and exceeded 90% with minimal curation. As expected, the discriminative power of this approach was lower at higher taxonomic ranks, but still significantly larger than zero. Simple instructions for calculation of a genomic distance between two genomes and species similarity thresholds at different k-mer sizes are suggested. The calculation of genomic distance is identified as a powerful approach for delineating fungal species and is proposed—not as the only criterion—but as an additional tool in the versatile toolbox of fungal taxonomy.

The Evolutionary Success of the Marine Bacterium SAR11 Analyzed through a Metagenomic Perspective

As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

The SAR11 clade of Alphaproteobacteria is the most abundant group of planktonic cells in the near-surface epipelagic waters of the ocean, but the mechanisms underlying its exceptional success have not been fully elucidated. Here, we applied a metagenomic approach to explore microdiversity patterns by measuring the accumulation of synonymous and nonsynonymous mutations as well as homologous recombination in populations of SAR11 from different aquatic habitats (marine epipelagic, bathypelagic, and surface freshwater). The patterns of mutation accumulation and recombination were compared to those of other groups of representative marine microbes with multiple ecological strategies that share the same marine habitat, namely, Cyanobacteria (Prochlorococcus and Synechococcus), Archaea (“Candidatus Nitrosopelagicus” and Marine Group II Thalassoarchaea), and some heterotrophic marine bacteria (Alteromonas and Erythrobacter). SAR11 populations showed widespread recombination among distantly related members, preventing divergence leading to a genetically stable population. Moreover, their high intrapopulation sequence diversity with an enrichment in synonymous replacements supports the idea of a very ancient divergence and the coexistence of multiple different clones. However, other microbes analyzed seem to follow different evolutionary dynamics where processes of diversification driven by geographic and ecological instability produce a higher number of nonsynonymous replacements and lower intrapopulation sequence diversity. Together, these data shed light on some of the evolutionary and ecological processes that lead to the large genomic diversity in SAR11. Furthermore, this approach can be applied to other similar microbes that are difficult to culture in the laboratory, but abundant in nature, to investigate the underlying dynamics of their genomic evolution.

IMPORTANCE As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

Stabilization of extensive fine-scale diversity by ecologically driven spatiotemporal chaos

It has recently become apparent that the diversity of microbial life extends far below the species level to the finest scales of genetic differences. Remarkably, extensive fine-scale diversity can coexist spatially. How is this diversity stable on long timescales, despite selective or ecological differences and other evolutionary processes? Most work has focused on stable coexistence or assumed ecological neutrality. We present an alternative: extensive diversity maintained by ecologically driven spatiotemporal chaos, with no assumptions about niches or other specialist differences between strains. We study generalized Lotka-Volterra models with antisymmetric correlations in the interactions inspired by multiple pathogen strains infecting multiple host strains. Generally, these exhibit chaos with increasingly wild population fluctuations driving extinctions. But the simplest spatial structure, many identical islands with migration between them, stabilizes a diverse chaotic state. Some strains (subspecies) go globally extinct, but many persist for times exponentially long in the number of islands. All persistent strains have episodic local blooms to high abundance, crucial for their persistence as, for many, their average population growth rate is negative. Snapshots of the abundance distribution show a power law at intermediate abundances that is essentially indistinguishable from the neutral theory of ecology. But the dynamics of the large populations are much faster than birth-death fluctuations. We argue that this spatiotemporally chaotic "phase" should exist in a wide range of models, and that even in rapidly mixed systems, longer-lived spores could similarly stabilize a diverse chaotic phase.

Potential Enterotoxicity of Phylogenetically Diverse Bacillus cereus Sensu Lato Soil Isolates from Different Geographical Locations

Bacillus cereus sensu lato comprises Gram-positive spore-forming bacteria producing toxins associated with foodborne diseases. Three pore-forming enterotoxins, nonhemolytic enterotoxin (Nhe), hemolysin BL (Hbl), and cytotoxin K (CytK), are considered the primary factors in B. cereus sensu lato diarrhea. The aim of this study was to determine the potential risk of enterotoxicity among soil B. cereus sensu lato isolates representing diverse phylogroups and originated from different geographic locations with various climates (Burkina Faso, Kenya, Argentina, Kazakhstan, and Poland). While nheA- and hblA-positive isolates were present among all B. cereus sensu lato populations and distributed across all phylogenetic groups, cytK-2-positive strains predominated in geographic regions with an arid hot climate (Africa) and clustered together on a phylogenetic tree mainly within mesophilic groups III and IV. The highest in vitro cytotoxicity to Caco-2 and HeLa cells was demonstrated by the strains clustered within phylogroups II and IV. Overall, our results suggest that B. cereus sensu lato pathogenicity is a comprehensive process conditioned by many intracellular factors and diverse environmental conditions.IMPORTANCE This research offers a new route for a wider understanding of the dependency between pathogenicity and phylogeny of a natural bacterial population, specifically within Bacillus cereus sensu lato, that is widely distributed around the world and easily transferred into food products. Our study indicates differences in the phylogenetic and geographical distributions of potential enterotoxigenic B. cereus sensu lato strains. Hence, these bacilli possess a risk for human health, and rapid testing methods for their identification are greatly needed. In particular, the detection of the CytK enterotoxin should be a supporting strategy for the identification of pathogenic B. cereus sensu lato.

Selection Is a Significant Driver of Gene Gain and Loss in the Pangenome of the Bacterial Genus Sulfurovum in Geographically Distinct Deep-Sea Hydrothermal Vents

Microbial genomes have highly variable gene content, and the evolutionary history of microbial populations is shaped by gene gain and loss mediated by horizontal gene transfer and selection. To evaluate the influence of selection on gene content variation in hydrothermal vent microbial populations, we examined 22 metagenome-assembled genomes (MAGs) (70 to 97% complete) from the ubiquitous vent Epsilonbacteraeota genus Sulfurovum that were recovered from two deep-sea hydrothermal vent regions, Axial Seamount in the northeastern Pacific Ocean (13 MAGs) and the Mid-Cayman Rise in the Caribbean Sea (9 MAGs). Genes involved in housekeeping functions were highly conserved across Sulfurovum lineages. However, genes involved in environment-specific functions, and in particular phosphate regulation, were found mostly in Sulfurovum genomes from the Mid-Cayman Rise in the low-phosphate Atlantic Ocean environment, suggesting that nutrient limitation is an important selective pressure for these bacteria. Furthermore, genes that were rare within the pangenome were more likely to undergo positive selection than genes that were highly conserved in the pangenome, and they also appeared to have experienced gene-specific sweeps. Our results suggest that selection is a significant driver of gene gain and loss for dominant microbial lineages in hydrothermal vents and highlight the importance of factors like nutrient limitation in driving microbial adaptation and evolution.IMPORTANCE Microbes can alter their gene content through the gain and loss of genes. However, there is some debate as to whether natural selection or neutral processes play a stronger role in molding the gene content of microbial genomes. In this study, we examined variation in gene content for the Epsilonbacteraeota genus Sulfurovum from deep-sea hydrothermal vents, which are dynamic habitats known for extensive horizontal gene transfer within microbial populations. Our results show that natural selection is a strong driver of Sulfurovum gene content and that nutrient limitation in particular has shaped the Sulfurovum genome, leading to differences in gene content between ocean basins. Our results also suggest that recently acquired genes undergo stronger selection than genes that were acquired in the more distant past. Overall, our results highlight the importance of natural selection in driving the evolution of microbial populations in these dynamic habitats.

Biogeography of American Northwest Hot Spring A/B-Lineage Synechococcus Populations

Previous analyses have shown how diversity among unicellular cyanobacteria inhabiting island-like hot springs is structured relative to physical separation and physiochemical differences among springs, especially at local to regional scales. However, these studies have been limited by the low resolution provided by the molecular markers surveyed. We analyzed large datasets obtained by high-throughput sequencing of a segment of the photosynthesis gene psaA from samples collected in hot springs from geothermal basins in Yellowstone National Park, Montana, and Oregon, all known from previous studies to contain populations of A/B'-lineage Synechococcus. The fraction of identical sequences was greater among springs separated by <50 km than among springs separated by >50 km, and springs separated by >800 km shared sequence variants only rarely. Phylogenetic analyses provided evidence for endemic lineages that could be related to geographic isolation and/or geochemical differences on regional scales. Ecotype Simulation 2 was used to predict putative ecotypes (ecologically distinct populations), and their membership, and canonical correspondence analysis was used to examine the geographical and geochemical bases for variation in their distribution. Across the range of Oregon and Yellowstone, geographical separation explained the largest percentage of the differences in distribution of ecotypes (9.5% correlated to longitude; 9.4% to latitude), with geochemical differences explaining the largest percentage of the remaining differences in distribution (7.4-9.3% correlated to magnesium, sulfate, and sulfide). Among samples within the Greater Yellowstone Ecosystem, geochemical differences significantly explained the distribution of ecotypes (6.5-9.3% correlated to magnesium, boron, sulfate, silicon dioxide, chloride, and pH). Nevertheless, differences in the abundance and membership of ecotypes in Yellowstone springs with similar chemistry suggested that allopatry may be involved even at local scales. Synechococcus populations have diverged both by physical isolation and physiochemical differences, and populations on surprisingly local scales have been evolving independently.

Advancing the understanding of treponemal disease in the past and present

Syphilis was perceived to be a new disease in Europe in the late 15th century, igniting a debate about its origin that continues today in anthropological, historical, and medical circles. We move beyond this age-old debate using an interdisciplinary approach that tackles broader questions to advance the understanding of treponemal infection (syphilis, yaws, bejel, and pinta). How did the causative organism(s) and humans co-evolve? How did the related diseases caused by Treponema pallidum emerge in different parts of the world and affect people across both time and space? How are T. pallidum subspecies related to the treponeme causing pinta? The current state of scholarship in specific areas is reviewed with recommendations made to stimulate future work. Understanding treponemal biology, genetic relationships, epidemiology, and clinical manifestations is crucial for vaccine development today and for investigating the distribution of infection in both modern and past populations. Paleopathologists must improve diagnostic criteria and use a standard approach for recording skeletal lesions on archaeological human remains. Adequate contextualization of cultural and environmental conditions is necessary, including site dating and justification for any corrections made for marine or freshwater reservoir effects. Biogeochemical analyses may assess aquatic contributions to diet, physiological changes arising from treponemal disease and its treatments (e.g., mercury), or residential mobility of those affected. Shifting the focus from point of origin to investigating who is affected (e.g., by age/sex or socioeconomic status) and disease distribution (e.g., coastal/ inland, rural/urban) will advance our understanding of the treponemal disease and its impact on people through time.

Comparative Genomics Analysis of Lactobacillus ruminis from Different Niches

Lactobacillus ruminis is a commensal motile lactic acid bacterium living in the intestinal tract of humans and animals. Although a few genomes of L. ruminis were published, most of them were animal derived. To explore the genetic diversity and potential niche-specific adaptation changes of L. ruminis, in the current work, draft genomes of 81 L. ruminis strains isolated from human, bovine, piglet, and other animals were sequenced, and comparative genomic analysis was performed. The genome size and GC content of L. ruminis on average were 2.16 Mb and 43.65%, respectively. Both the origin and the sampling distance of these strains had a great influence on the phylogenetic relationship. For carbohydrate utilization, the human-derived L. ruminis strains had a higher consistency in the utilization of carbon source compared to the animal-derived strains. L. ruminis mainly increased the competitiveness of niches by producing class II bacteriocins. The type of clustered regularly interspaced short palindromic repeats /CRISPR-associated (CRISPR/Cas) system presented in L. ruminis was mainly subtype IIA. The diversity of CRISPR/Cas locus depended on the high denaturation of spacer number and sequence, although cas1 protein was relatively conservative. The genetic differences in those newly sequenced L. ruminis strains highlighted the gene gains and losses attributed to niche adaptations.

Single-amino acid variants reveal evolutionary processes that shape the biogeography of a global SAR11 subclade

Members of the SAR11 order Pelagibacterales dominate the surface oceans. Their extensive diversity challenges emerging operational boundaries defined for microbial 'species' and complicates efforts of population genetics to study their evolution. Here, we employed single-amino acid variants (SAAVs) to investigate ecological and evolutionary forces that maintain the genomic heterogeneity within ubiquitous SAR11 populations we accessed through metagenomic read recruitment using a single isolate genome. Integrating amino acid and protein biochemistry with metagenomics revealed that systematic purifying selection against deleterious variants governs non-synonymous variation among very closely related populations of SAR11. SAAVs partitioned metagenomes into two main groups matching large-scale oceanic current temperatures, and six finer proteotypes that connect distant oceanic regions. These findings suggest that environmentally-mediated selection plays a critical role in the journey of cosmopolitan surface ocean microbial populations, and the idea 'everything is everywhere but the environment selects' has credence even at the finest resolutions.

Adaptation of host transmission cycle during Clostridium difficile speciation

Bacterial speciation is a fundamental evolutionary process characterized by diverging genotypic and phenotypic properties. However, the selective forces that affect genetic adaptations and how they relate to the biological changes that underpin the formation of a new bacterial species remain poorly understood. Here, we show that the spore-forming, healthcare-associated enteropathogen Clostridium difficile is actively undergoing speciation. Through large-scale genomic analysis of 906 strains, we demonstrate that the ongoing speciation process is linked to positive selection on core genes in the newly forming species that are involved in sporulation and the metabolism of simple dietary sugars. Functional validation shows that the new C. difficile produces spores that are more resistant and have increased sporulation and host colonization capacity when glucose or fructose is available for metabolism. Thus, we report the formation of an emerging C. difficile species, selected for metabolizing simple dietary sugars and producing high levels of resistant spores, that is adapted for healthcare-mediated transmission.

Ecological Divergence with Gene Flow in a Thermophilic Cyanobacterium

How ecological diversity is maintained and distributed within populations is a longstanding question in microbial ecology. In the thermophilic cyanobacterium Synechococcus B', high observed levels of recombination are predicted to maintain ecological variation despite the simultaneous action of diverse selective pressures on different regions of the genome. To investigate ecological diversity in these bacteria, we directly isolated laboratory strains of Synechococcus B' from samples collected along the thermal gradients of two geothermal environments in Yellowstone National Park. Extensive recombination was evident for a multi-locus sequence data set, and, consequently, our sample did not exhibit the sequence clustering expected for distinct ecotypes evolving by periodic clonal selection. Evidence for local selective sweeps at specific loci suggests that sweeps may be common but that recombination is effective for maintaining diversity of unlinked genomic regions. Thermal performance for strain growth was positively associated with the temperature of the environment, indicating that Synechococcus B' populations consist of locally adapted ecological specialists that occupy specific thermal niches. Because this ecological differentiation is observed despite the absence of dispersal barriers among sites, we conclude that these bacteria may freely exchange much of the genome but that barriers to gene flow exist for loci under direct temperature selection.

Delineation of ecologically distinct units of marine Bacteroidetes in the Northwestern Mediterranean Sea

Bacteroidetes is one of the dominant phyla of ocean bacterioplankton, yet its diversity and population structure is poorly understood. To advance in the delineation of ecologically meaningful units within this group, we constructed near full-length 16S rRNA gene clone libraries from contrasting marine environments in the NW Mediterranean. Based on phylogeny and the associated ecological variables (depth and season), 24 different Bacteroidetes clades were delineated. By considering their relative abundance (from iTag amplicon sequencing studies), we described the distribution patterns of each of these clades, delimiting them as Ecologically Significant Taxonomic Units (ESTUs). Spatially, there was almost no overlap among ESTUs at different depths. In deep waters there was predominance of Owenweeksia, Leeuwenhoekiella, Muricauda-related genera, and some depth-associated ESTUs within the NS5 and NS2b marine clades. Seasonally, multi-annual dynamics of recurring ESTUs were present with dominance of some ESTUs within the NS4, NS5 and NS2b marine clades along most of the year, but with variable relative frequencies between months. A drastic change towards the predominance of Formosa-related ESTUs and one ESTU from the NS5 marine clade was typically present after the spring bloom. Even though there are no isolates available for these ESTUs to determine their physiology, correlation analyses identified the environmental preference of some of them. Overall, our results suggest that there is a high degree of niche specialisation within these closely related clades. This work constitutes a step forward in disentangling the ecology of marine Bacteroidetes, which are essential players in organic matter processing in the oceans.

The role of recombination, niche-specific gene pools and flexible genomes in the ecological speciation of bacteria

Bacteria diversify into genetic clusters analogous to those observed in sexual eukaryotes, but the definition of bacterial species is an ongoing problem. Recent work has focused on adaptation to distinct ecological niches as the main driver of clustering, but there remains debate about the role of recombination in that process. One view is that homologous recombination occurs too rarely for gene flow to constrain divergent selection. Another view is that homologous recombination is frequent enough in many bacterial populations that barriers to gene flow are needed to permit divergence. Niche-specific gene pools have been proposed as a general mechanism to limit gene flow. We use theoretical models to evaluate additional hypotheses that evolving genetic architecture, specifically the effect sizes of genes and gene gain and loss, can limit gene flow between diverging populations. Our model predicts that (a) in the presence of gene flow and recombination, ecological divergence is concentrated in few loci of large effect and (b) high rates of gene flow plus recombination promote gene loss and favor the evolution of niche-specific genes. The results show that changing genetic architecture and gene loss can facilitate ecological divergence, even without niche-specific gene pools. We discuss these results in the context of recent studies of sympatric divergence in microbes.

Population Genomics Insights into Adaptive Evolution and Ecological Differentiation in Streptomycetes

Deciphering the genomic variation that represents microevolutionary processes toward species divergence is key to understanding microbial speciation, which has long been under debate. Streptomycetes are filamentous bacteria that are ubiquitous in nature and the richest source of antibiotics; however, their speciation processes remain unknown. To tackle this issue, we performed a comprehensive population genomics analysis on Streptomyces albidoflavus residing in different habitats and with a worldwide distribution and identified and characterized the foundational changes within the species. We detected three well-defined phylogenomic clades, of which clades I and III mainly contained free-living (soil/marine) and insect-associated strains, respectively, and clade II had a mixed origin. By performing genome-wide association studies (GWAS), we identified a number of genetic variants associated with free-living or entomic (denoting or relating to insects) habitats in both the accessory and core genomes. These variants contributed collectively to the population structure and had annotated or confirmed functions that likely facilitate differential adaptation of the species. In addition, we detected higher levels of homologous recombination within each clade and in the free-living group than within the whole species and in the entomic group. A subset of the insect-associated strains (clade III) showed a relatively independent evolutionary trajectory with more symbiosis-favorable genes but little genetic interchange with the other lineages. Our results demonstrate that ecological adaptation promotes genetic differentiation in S. albidoflavus, suggesting a model of ecological speciation with gene flow in streptomycetes.IMPORTANCE Species are the fundamental units of ecology and evolution, and speciation leads to the astounding diversity of life on Earth. Studying speciation is thus of great significance to understand, protect, and exploit biodiversity, but it is a challenge in the microbial world. In this study, using population genomics, we placed Streptomyces albidoflavus strains in a spectrum of speciation and showed that the genetic differences between phylogenomic clusters evolved mainly by environmental selection and gene-specific sweeps. These findings highlight the role of ecology in structuring recombining bacterial species, making a step toward a deeper understanding of microbial speciation. Our results also raise concerns of an underrated microbial diversity at the intraspecies level, which can be utilized for mining of ecologically relevant natural products.

Historical Factors Associated With Past Environments Influence the Biogeography of Thermophilic Endospores in Arctic Marine Sediments

Selection by the local, contemporary environment plays a prominent role in shaping the biogeography of microbes. However, the importance of historical factors in microbial biogeography is more debatable. Historical factors include past ecological and evolutionary circumstances that may have influenced present-day microbial diversity, such as dispersal and past environmental conditions. Diverse thermophilic sulfate-reducing Desulfotomaculum are present as dormant endospores in marine sediments worldwide where temperatures are too low to support their growth. Therefore, they are dispersed to here from elsewhere, presumably a hot, anoxic habitat. While dispersal through ocean currents must influence their distribution in cold marine sediments, it is not clear whether even earlier historical factors, related to the source habitat where these organisms were once active, also have an effect. We investigated whether these historical factors may have influenced the diversity and distribution of thermophilic endospores by comparing their diversity in 10 Arctic fjord surface sediments. Although community composition varied spatially, clear biogeographic patterns were only evident at a high level of taxonomic resolution (>97% sequence similarity of the 16S rRNA gene) achieved with oligotyping. In particular, the diversity and distribution of oligotypes differed for the two most prominent OTUs (defined using a standard 97% similarity cutoff). One OTU was dominated by a single ubiquitous oligotype, while the other OTU consisted of ten more spatially localized oligotypes that decreased in compositional similarity with geographic distance. These patterns are consistent with differences in historical factors that occurred when and where the taxa were once active, prior to sporulation. Further, the influence of history on biogeographic patterns was only revealed by analyzing microdiversity within OTUs, suggesting that populations within standard OTU-level groupings do not necessarily share a common ecological and evolutionary history.

Probing the ecological and evolutionary history of a thermophilic cyanobacterial population via statistical properties of its microdiversity

Despite extensive DNA sequencing data derived from natural microbial communities, it remains a major challenge to identify the key evolutionary and ecological forces that shape microbial populations. We have focused on the extensive microdiversity of the cyanobacterium Synechococcus sp., which is a dominant member of the dense phototrophic biofilms in the hot springs of Yellowstone National Park. From deep amplicon sequencing of many loci and statistical analyses of these data, we showed previously that the population has undergone an unexpectedly high degree of homologous recombination, unlinking synonymous SNP-pair correlations even on intragenic length scales. Here, we analyze the genic amino acid diversity, which provides new evidence of selection and insights into the evolutionary history of the population. Surprisingly, some features of the data, including the spectrum of distances between genic-alleles, appear consistent with primarily asexual neutral drift. Yet the non-synonymous site frequency spectrum has too large an excess of low-frequency polymorphisms to result from negative selection on deleterious mutations given the distribution of coalescent times that we infer. And our previous analyses showed that the population is not asexual. Taken together, these apparently contradictory data suggest that selection, epistasis, and hitchhiking all play essential roles in generating and stabilizing the diversity. We discuss these as well as potential roles of ecological niches at genomic and genic levels. From quantitative properties of the diversity and comparative genomic data, we infer aspects of the history and inter-spring dispersal of the meta-population since it was established in the Yellowstone Caldera. Our investigations illustrate the need for combining multiple types of sequencing data and quantitative statistical analyses to develop an understanding of microdiversity in natural microbial populations.

Temperature modulates Fischerella thermalis ecotypes in Porcelana Hot Spring

In the Porcelana Hot Spring (Northern Patagonia), true-branching cyanobacteria are the dominant primary producers in microbial mats, and they are mainly responsible for carbon and nitrogen fixation. However, little is known about their metabolic and genomic adaptations at high temperatures. Therefore, in this study, a total of 81 Fischerella thermalis strains (also known as Mastigocladus laminosus) were isolated from mat samples in a thermal gradient between 61-46°C. The complementary use of proteomic comparisons from these strains, and comparative genomics of F. thermalis pangenomes, suggested that at least two different ecotypes were present within these populations. MALDI-TOF MS analysis separated the strains into three clusters; two with strains obtained from mats within the upper temperature range (61 and 54°C), and a third obtained from mats within the lower temperature range (51 and 46°C). Both groups possessed different but synonymous nifH alleles. The main proteomic differences were associated with the abundance of photosynthesis-related proteins. Three F. thermalis metagenome assembled genomes (MAGs) were described from 66, 58 and 48°C metagenomes. These pangenomes indicated a divergence of orthologous genes and a high abundance of exclusive genes at 66°C. These results improved the current understanding of thermal adaptation of F. thermalis and the evolution of these thermophilic cyanobacterial species.

Transmission in the Origins of Bacterial Diversity, From Ecotypes to Phyla

Any two lineages, no matter how distant they are now, began their divergence as one population splitting into two lineages that could coexist indefinitely. The rate of origin of higher-level taxa is therefore the product of the rate of speciation times the probability that two new species coexist long enough to reach a particular level of divergence. Here I have explored these two parameters of disparification in bacteria. Owing to low recombination rates, sexual isolation is not a necessary milestone of bacterial speciation. Rather, irreversible and indefinite divergence begins with ecological diversification, that is, transmission of a bacterial lineage to a new ecological niche, possibly to a new microhabitat but at least to new resources. Several algorithms use sequence data from a taxon of focus to identify phylogenetic groups likely to bear the dynamic properties of species. Identifying these newly divergent lineages allows us to characterize the genetic bases of speciation, as well as the ecological dimensions upon which new species diverge. Speciation appears to be least frequent when a given lineage has few new resources it can adopt, as exemplified by photoautotrophs, C1 heterotrophs, and obligately intracellular pathogens; speciation is likely most rapid for generalist heterotrophs. The genetic basis of ecological divergence may determine whether ecological divergence is irreversible and whether lineages will diverge indefinitely into the future. Long-term coexistence is most likely when newly divergent lineages utilize at least some resources not shared with the other and when the resources themselves will coexist into the remote future.

Trajectories and Drivers of Genome Evolution in Surface-Associated Marine Phaeobacter

The extent of genome divergence and the evolutionary events leading to speciation of marine bacteria have mostly been studied for (locally) abundant, free-living groups. The genus Phaeobacter is found on different marine surfaces, seems to occupy geographically disjunct habitats, and is involved in different biotic interactions, and was therefore targeted in the present study. The analysis of the chromosomes of 32 closely related but geographically spread Phaeobacter strains revealed an exceptionally large, highly syntenic core genome. The flexible gene pool is constantly but slightly expanding across all Phaeobacter lineages. The horizontally transferred genes mostly originated from bacteria of the Roseobacter group and horizontal transfer most likely was mediated by gene transfer agents. No evidence for geographic isolation and habitat specificity of the different phylogenomic Phaeobacter clades was detected based on the sources of isolation. In contrast, the functional gene repertoire and physiological traits of different phylogenomic Phaeobacter clades were sufficiently distinct to suggest an adaptation to an associated lifestyle with algae, to additional nutrient sources, or toxic heavy metals. Our study reveals that the evolutionary trajectories of surface-associated marine bacteria can differ significantly from free-living marine bacteria or marine generalists.

Clustering of Pan- and Core-genome of Lactobacillus provides Novel Evolutionary Insights for Differentiation

Background

Bacterial taxonomy aims to classify bacteria based on true evolutionary events and relies on a polyphasic approach that includes phenotypic, genotypic and chemotaxonomic analyses. Until now, complete genomes are largely ignored in taxonomy. The genus Lactobacillus consists of 173 species and many genomes are available to study taxonomy and evolutionary events.

Results

We analyzed and clustered 98 completely sequenced genomes of the genus Lactobacillus and 234 draft genomes of 5 different Lactobacillus species, i.e. L. reuteri, L. delbrueckii, L. plantarum, L. rhamnosus and L. helveticus. The core-genome of the genus Lactobacillus contains 266 genes and the pan-genome 20′800 genes. Clustering of the Lactobacillus pan- and core-genome resulted in two highly similar trees. This shows that evolutionary history is traceable in the core-genome and that clustering of the core-genome is sufficient to explore relationships. Clustering of core- and pan-genomes at species’ level resulted in similar trees as well. Detailed analyses of the core-genomes showed that the functional class “genetic information processing” is conserved in the core-genome but that “signaling and cellular processes” is not. The latter class encodes functions that are involved in environmental interactions. Evolution of lactobacilli seems therefore directed by the environment. The type species L. delbrueckii was analyzed in detail and its pan-genome based tree contained two major clades whose members contained different genes yet identical functions. In addition, evidence for horizontal gene transfer between strains of L. delbrueckii, L. plantarum, and L. rhamnosus, and between species of the genus Lactobacillus is presented. Our data provide evidence for evolution of some lactobacilli according to a parapatric-like model for species differentiation.

Conclusions

Core-genome trees are useful to detect evolutionary relationships in lactobacilli and might be useful in taxonomic analyses. Lactobacillus’ evolution is directed by the environment and HGT.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4601-5) contains supplementary material, which is available to authorized users.

Phylogenomic analysis of the species of the Mycobacterium tuberculosis complex demonstrates that Mycobacterium africanum, Mycobacterium bovis, Mycobacterium caprae, Mycobacterium microti and Mycobacterium pinnipedii are later heterotypic synonyms of Mycobacterium tuberculosis

The species within the Mycobacterium tuberculosis Complex (MTBC) have undergone numerous taxonomic and nomenclatural changes, leaving the true structure of the MTBC in doubt. We used next-generation sequencing (NGS), digital DNA-DNA hybridization (dDDH), and average nucleotide identity (ANI) to investigate the relationship between these species. The type strains of Mycobacterium africanum, Mycobacterium bovis, Mycobacterium caprae, Mycobacterium microti and Mycobacterium pinnipedii were sequenced via NGS. Pairwise dDDH and ANI comparisons between these, previously sequenced MTBC type strain genomes (including 'Mycobacterium canettii', 'Mycobacterium mungi' and 'Mycobacterium orygis') and M. tuberculosis H37Rv^T were performed. Further, all available genome sequences in GenBank for species in or putatively in the MTBC were compared to H37Rv^T. Pairwise results indicated that all of the type strains of the species are extremely closely related to each other (dDDH: 91.2-99.2 %, ANI: 99.21-99.92 %), greatly exceeding the respective species delineation thresholds, thus indicating that they belong to the same species. Results from the GenBank genomes indicate that all the strains examined are within the circumscription of H37Rv^T (dDDH: 83.5-100 %). We, therefore, formally propose a union of the species of the MTBC as M. tuberculosis. M. africanum, M. bovis, M. caprae, M. microti and M. pinnipedii are reclassified as later heterotypic synonyms of M. tuberculosis. 'M. canettii', 'M. mungi', and 'M. orygis' are classified as strains of the species M. tuberculosis. We further recommend use of the infrasubspecific term 'variant' ('var.') and infrasubspecific designations that generally retain the historical nomenclature associated with the groups or otherwise convey such characteristics, e.g. M. tuberculosis var. bovis.

Genomic divergence and cohesion in a species of pelagic freshwater bacteria

Background

In many prokaryotic genera a clustered phylogeny is observed, akin to the occurrence of species in sexually reproducing organisms. For some taxa, homologous recombination has been invoked as the underlying mechanism providing genomic cohesion among conspecific individuals. Whether this mechanism is applicable to prokaryotes in freshwaters with low habitat connectivity – i.e. elevated geographic barriers to gene flow – is unclear. To investigate further we studied genomic trends within the globally abundant PnecC cluster (genus Polynucleobacter, Betaproteobacteria) and analyzed homologous recombination within the affiliated species P. asymbioticus.

Results

Comparisons among 20 PnecC genomes revealed a clearly discontinuous distribution of nucleotide sequence similarities. Among the nine conspecific individuals (P. asymbioticus) all average nucleotide identity (ANI) values were greater than 97%, whereas all other comparisons exhibited ANI values lower than 85%. The reconstruction of recombination and mutation events for the P. asymbioticus core genomes yielded an r/m ratio of 7.4, which is clearly above estimated thresholds for recombination to act as a cohesive force. Hotspots of recombination were found to be located in the flanking regions of genomic islands. Even between geographically separated habitats a high flux of recombination was evident. While a biogeographic population structure was suggested from MLST data targeting rather conserved loci, such a structure was barely visible when whole genome data was considered. However, both MLST and whole genome data showed evidence of differentiation between two lineages of P. asymbioticus. The ratios of non-synonymous to synonymous substitution rates as well as growth rates in transplantation experiments suggested that this divergence was not selectively neutral.

Conclusions

The high extent of homologous recombination among P. asymbioticus bacteria can act as a cohesive force that effectively counteracts genetic divergence. At least on a regional scale, homologous recombination can act across geographically separated ecosystems and therefore plays an important role in the evolution and consistency of bacterial freshwater species. A species model akin to the biological species concept may be applicable for P. asymbioticus. Nonetheless, two genetically distinct lineages have emerged and further research may clarify if their divergence has been initiated by reinforced geographical barriers or has been evolving in sympatry.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4199-z) contains supplementary material, which is available to authorized users.

Theoretical microbial ecology without species

Ecosystems are commonly conceptualized as networks of interacting species. However, partitioning natural diversity of organisms into discrete units is notoriously problematic and mounting experimental evidence raises the intriguing question whether this perspective is appropriate for the microbial world. Here an alternative formalism is proposed that does not require postulating the existence of species as fundamental ecological variables and provides a naturally hierarchical description of community dynamics. This formalism allows approaching the species problem from the opposite direction. While the classical models treat a world of imperfectly clustered organism types as a perturbation around well-clustered species, the presented approach allows gradually adding structure to a fully disordered background. The relevance of this theoretical construct for describing highly diverse natural ecosystems is discussed.

A null model for microbial diversification

Whether prokaryotes (Bacteria and Archaea) are naturally organized into phenotypically and genetically cohesive units comparable to animal or plant species remains contested, frustrating attempts to estimate how many such units there might be, or to identify the ecological roles they play. Analyses of gene sequences in various closely related prokaryotic groups reveal that sequence diversity is typically organized into distinct clusters, and processes such as periodic selection and extensive recombination are understood to be drivers of cluster formation ("speciation"). However, observed patterns are rarely compared with those obtainable with simple null models of diversification under stochastic lineage birth and death and random genetic drift. Via a combination of simulations and analyses of core and phylogenetic marker genes, we show that patterns of diversity for the genera Escherichia, Neisseria, and Borrelia are generally indistinguishable from patterns arising under a null model. We suggest that caution should thus be taken in interpreting observed clustering as a result of selective evolutionary forces. Unknown forces do, however, appear to play a role in Helicobacter pylori, and some individual genes in all groups fail to conform to the null model. Taken together, we recommend the presented birth-death model as a null hypothesis in prokaryotic speciation studies. It is only when the real data are statistically different from the expectations under the null model that some speciation process should be invoked.