Search

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., Competing Interests: Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

What a strain is and how many strains make up a natural bacterial population remain elusive concepts despite their apparent importance for assessing the role of intra-population diversity in disease emergence or response to environmental perturbations. To advance these concepts, we sequenced 138 randomly selected Salinibacter ruber isolates from two solar salterns and assessed these genomes against companion short-read metagenomes from the same samples. The distribution of genome-aggregate average nucleotide identity (ANI) values among these isolates revealed a bimodal distribution, with four-fold lower occurrence of values between 99.2% and 99.8% relative to ANI >99.8% or <99.2%, revealing a natural "gap" in the sequence space within species. Accordingly, we used this ANI gap to define genomovars and a higher ANI value of >99.99% and shared gene-content >99.0% to define strains. Using these thresholds and extrapolating from how many metagenomic reads each genomovar uniquely recruited, we estimated that -although our 138 isolates represented about 80% of the Sal. ruber population- the total population in one saltern pond is composed of 5,500 to 11,000 genomovars, the great majority of which appear to be rare in-situ. These data also revealed that the most frequently recovered isolate in lab media was often not the most abundant genomovar in-situ, suggesting that cultivation biases are significant, even in cases that cultivation procedures are thought to be robust. The methodology and ANI thresholds outlined here should represent a useful guide for future microdiversity surveys of additional microbial species., (© 2024. The Author(s).)

The specialization-disturbance hypothesis predicts that, in the event of a disturbance, generalists are favored, while specialists are selected against. This hypothesis has not been rigorously tested in microbial systems and it remains unclear to what extent it could explain microbial community succession patterns following perturbations. Previous field observations of Pensacola Beach sands that were impacted by the Deepwater Horizon (DWH) oil spill provided evidence in support of the specialization-disturbance hypothesis. However, ecological drift as well as uncounted environmental fluctuations (e.g., storms) could not be ruled out as confounding factors driving these field results. In this study, the specialization-disturbance hypothesis was tested on beach sands, disturbed by DWH crude oil, ex situ in closed laboratory advective-flow chambers that mimic in situ conditions in saturated beach sediments. The chambers were inoculated with weathered DWH oil and unamended chambers served as controls. The time series of shotgun metagenomic and 16S rRNA gene amplicon sequence data from a two-month long incubation showed that functional diversity significantly increased while taxonomic diversity significantly declined, indicating a decrease in specialist taxa. Thus, results from this laboratory study corroborate field observations, providing verification that the specialization-disturbance hypothesis can explain microbial succession patterns in crude oil impacted beach sands., (© 2021. The Author(s).)

Modeling crude-oil biodegradation in sediments remains a challenge due in part to the lack of appropriate model organisms. Here we report the metagenome-guided isolation of a novel organism that represents a phylogenetically narrow (>97% 16S rRNA gene identity) group of previously uncharacterized, crude-oil degraders. Analysis of available sequence data showed that these organisms are highly abundant in oiled sediments of coastal marine ecosystems across the world, often comprising ~30% of the total community, and virtually absent in pristine sediments or seawater. The isolate genome encodes functional nitrogen fixation and hydrocarbon degradation genes together with putative genes for biosurfactant production that apparently facilitate growth in the typically nitrogen-limited, oiled environment. Comparisons to available genomes revealed that this isolate represents a novel genus within the Gammaproteobacteria, for which we propose the provisional name "Candidatus Macondimonas diazotrophica" gen. nov., sp. nov. "Ca. M. diazotrophica" appears to play a key ecological role in the response to oil spills around the globe and could be a promising model organism for studying ecophysiological responses to oil spills.

Recent diversity studies have revealed that microbial communities of natural environments are dominated by species-like, sequence-discrete populations. However, how stable the sequence and gene-content diversity are within these populations and especially in highly dynamic lotic habitats remain unclear. Here we quantified the dynamics of intra-population diversity in samples spanning two years and five sites in the Kalamas River (Northwest Greece). A significant positive correlation was observed between higher intra-population sequence diversity and longer persistence over time, revealing that more diverse populations tended to represent more autochthonous (vs. allochthonous) community members. Assessment of intra-population gene-content changes caused by strain replacement or gene loss over time revealed different profiles with the majority of populations exhibiting gene-content changes close to 10% of the total genes, while one population exhibited ~21% change. The variable genes were enriched in hypothetical proteins and mobile elements, and thus, were probably functionally neutral or attributable to phage predation. A few notable exceptions to this pattern were also noted such as phototrophy-related proteins in summer vs. winter populations. Taken together, these results revealed that some freshwater genomes are remarkably dynamic, even across short time and spatial scales, and have implications for the bacterial species concept and microbial source tracking.

A fundamental question in microbiology is whether there is continuum of genetic diversity among genomes, or clear species boundaries prevail instead. Whole-genome similarity metrics such as Average Nucleotide Identity (ANI) help address this question by facilitating high resolution taxonomic analysis of thousands of genomes from diverse phylogenetic lineages. To scale to available genomes and beyond, we present FastANI, a new method to estimate ANI using alignment-free approximate sequence mapping. FastANI is accurate for both finished and draft genomes, and is up to three orders of magnitude faster compared to alignment-based approaches. We leverage FastANI to compute pairwise ANI values among all prokaryotic genomes available in the NCBI database. Our results reveal clear genetic discontinuity, with 99.8% of the total 8 billion genome pairs analyzed conforming to >95% intra-species and <83% inter-species ANI values. This discontinuity is manifested with or without the most frequently sequenced species, and is robust to historic additions in the genome databases.

Although petroleum hydrocarbons discharged from the Deepwater Horizon (DWH) blowout were shown to have a pronounced impact on indigenous microbial communities in the Gulf of Mexico, effects on nearshore or coastal ecosystems remain understudied. This study investigated the successional patterns of functional and taxonomic diversity for over 1 year after the DWH oil was deposited on Pensacola Beach sands (FL, USA), using metagenomic and 16S rRNA gene amplicon techniques. Gamma- and Alphaproteobacteria were enriched in oiled sediments, in corroboration of previous studies. In contrast to previous studies, we observed an increase in the functional diversity of the community in response to oil contamination and a functional transition from generalist populations within 4 months after oil came ashore to specialists a year later, when oil was undetectable. At the latter time point, a typical beach community had reestablished that showed little to no evidence of oil hydrocarbon degradation potential, was enriched in archaeal taxa known to be sensitive to xenobiotics, but differed significantly from the community before the oil spill. Further, a clear succession pattern was observed, where early responders to oil contamination, likely degrading aliphatic hydrocarbons, were replaced after 3 months by populations capable of aromatic hydrocarbon decomposition. Collectively, our results advance the understanding of how natural benthic microbial communities respond to crude oil perturbation, supporting the specialization-disturbance hypothesis; that is, the expectation that disturbance favors generalists, while providing (microbial) indicator species and genes for the chemical evolution of oil hydrocarbons during degradation and weathering.

Coral reefs are deteriorating at an alarming rate mainly as a consequence of the emergence of coral diseases. The white plague disease (WPD) is the most prevalent coral disease in the southwestern Caribbean, affecting dozens of coral species. However, the identification of a single causal agent has proved problematic. This suggests more complex etiological scenarios involving alterations in the dynamic interaction between environmental factors, the coral immune system and the symbiotic microbial communities. Here we compare the microbiome of healthy and WPD-affected corals from the two reef-building species Diploria strigosa and Siderastrea siderea collected at the Tayrona National Park in the Caribbean of Colombia. Microbiomes were analyzed by combining culture-dependent methods and pyrosequencing of 16S ribosomal DNA (rDNA) V5-V6 hypervariable regions. A total of 20,410 classifiable 16S rDNA sequences reads were obtained including all samples. No significant differences in operational taxonomic unit diversity were found between healthy and affected tissues; however, a significant increase of Alphaproteobacteria and a concomitant decrease in the Beta- and Gammaproteobacteria was observed in WPD-affected corals of both species. Significant shifts were also observed in the orders Rhizobiales, Caulobacteriales, Burkholderiales, Rhodobacterales, Aleteromonadales and Xanthomonadales, although they were not consistent between the two coral species. These shifts in the microbiome structure of WPD-affected corals suggest a loss of community-mediated growth control mechanisms on bacterial populations specific for each holobiont system.