Long-Read-Resolved, Ecosystem-Wide Exploration of Nucleotide and Structural Microdiversity of Lake Bacterioplankton Genomes

Two decades of bacterial ecology and evolution in a freshwater lake

Ecology and evolution are considered distinct processes that interact on contemporary time scales in microbiomes. Here, to observe these processes in a natural system, we collected a two-decade, 471-metagenome time series from Lake Mendota (Wisconsin, USA). We assembled 2,855 species-representative genomes and found that genomic change was common and frequent. By tracking strain composition via single nucleotide variants, we identified cyclical seasonal patterns in 80% and decadal shifts in 20% of species. In the dominant freshwater family Nanopelagicaceae, environmental extremes coincided with shifts in strain composition and positive selection of amino acid and nucleic acid metabolism genes. These genes identify organic nitrogen compounds as potential drivers of freshwater responses to global change. Seasonal and long-term strain dynamics could be regarded as ecological processes or, equivalently, as evolutionary change. Rather than as distinct interacting processes, we propose a conceptualization of ecology and evolution as a continuum to better describe change in microbial communities.

Genome-resolved year-round dynamics reveal a broad range of giant virus microdiversity

Giant viruses are crucial for marine ecosystem dynamics because they regulate microeukaryotic community structure, accelerate carbon and nutrient cycles, and drive the evolution of their hosts through co-evolutionary processes. Previously reported long-term observations revealed that these viruses display seasonal fluctuations in abundance. However, the underlying genetic mechanisms driving such dynamics of these viruses remain largely unknown. In this study, we investigated the dynamics of giant viruses using time-series metagenomes from eutrophic coastal seawater samples collected over 20 months. A newly developed computational pipeline generated 1,065 high-quality genomes covering six major giant virus lineages. These genomic data revealed year-round recovery of the viral community structure at the study site and distinct dynamics of viral populations that were classified as persistent (n = 9), seasonal (n = 389), sporadic (n = 318), or others. By profiling the intra-species nucleotide-resolved microdiversity through read mapping, we also identified year-round recovery dynamics at subpopulation level for viruses classified as persistent or seasonal. Our results further indicated that giant viruses with broader niche breadth tended to exhibit higher levels of microdiversity. We argue that greater microdiversity of viruses likely enhances adaptability and thus survival under the virus-host arms race during prolonged interactions with their hosts.IMPORTANCERecent genome-resolved metagenomic surveys have uncovered the vast genomic diversity of giant viruses, which play significant roles in aquatic ecosystems by acting as bloom terminators and influencing biogeochemical cycles. However, the relationship between the ecological dynamics of giant viruses and underlying genetic structures of viral populations remains unresolved. In this study, we performed deep metagenomic sequencing of seawater samples collected across a time-series from a coastal area in Japan. The results revealed a significant positive correlation between microdiversity and temporal persistence of giant virus populations, suggesting that population structure is a crucial factor for adaptation and survival in the interactions with their hosts.

Bacterial ecology and evolution converge on seasonal and decadal scales

Ecology and evolution are often viewed as distinct processes, which interact on contemporary time scales in microbiomes. To observe these processes in a natural system, we collected a two-decade, 471-sample freshwater lake time series, creating the longest metagenome dataset to date. Among 2,855 species-representative genomes, diverse species and strains followed cyclical seasonal patterns, and one in five species experienced decadal shifts in strain composition. The most globally abundant freshwater bacterium had constant species-level abundance, but environmental extremes appeared to trigger a shift in strain composition and positive selection of amino acid and nucleic acid metabolism genes. These genes identify organic nitrogen compounds as potential drivers of freshwater responses to global change. Seasonal and long-term strain dynamics could be regarded as ecological processes or equivalently as evolutionary change. Rather than as distinct processes that interact, we propose a conceptualization where ecology and evolution converge along a continuum to better describe change in diverse microbial communities.

Vaginal microbiomes show ethnic evolutionary dynamics and positive selection of Lactobacillus adhesins driven by a long-term niche-specific process

The vaginal microbiome's composition varies among ethnicities. However, the evolutionary landscape of the vaginal microbiome in the multi-ethnic context remains understudied. We perform a systematic evolutionary analysis of 351 vaginal microbiome samples from 35 multi-ethnic pregnant women, in addition to two validation cohorts, totaling 462 samples from 90 women. Microbiome alpha diversity and community state dynamics show strong ethnic signatures. Lactobacillaceae have a higher ratio of non-synonymous to synonymous polymorphism and lower nucleotide diversity than non-Lactobacillaceae in all ethnicities, with a large repertoire of positively selected genes, including the mucin-binding and cell wall anchor genes. These evolutionary dynamics are driven by the long-term evolutionary process unique to the human vaginal niche. Finally, we propose an evolutionary model reflecting the environmental niches of microbes. Our study reveals the extensive ethnic signatures in vaginal microbial ecology and evolution, highlighting the importance of studying the host-microbiome ecosystem from an evolutionary perspective.

Genome-resolved metagenomics of Venice Lagoon surface sediment bacteria reveals high biosynthetic potential and metabolic plasticity as successful strategies in an impacted environment

Bacteria living in sediments play essential roles in marine ecosystems and deeper insights into the ecology and biogeochemistry of these largely unexplored organisms can be obtained from 'omics' approaches. Here, we characterized metagenome-assembled-genomes (MAGs) from the surface sediment microbes of the Venice Lagoon (northern Adriatic Sea) in distinct sub-basins exposed to various natural and anthropogenic pressures. MAGs were explored for biodiversity, major marine metabolic processes, anthropogenic activity-related functions, adaptations at the microscale, and biosynthetic gene clusters. Starting from 126 MAGs, a non-redundant dataset of 58 was compiled, the majority of which (35) belonged to (Alpha- and Gamma-) Proteobacteria. Within the broad microbial metabolic repertoire (including C, N, and S metabolisms) the potential to live without oxygen emerged as one of the most important features. Mixotrophy was also found as a successful lifestyle. Cluster analysis showed that different MAGs encoded the same metabolic patterns (e.g., C fixation, sulfate oxidation) thus suggesting metabolic redundancy. Antibiotic and toxic compounds resistance genes were coupled, a condition that could promote the spreading of these genetic traits. MAGs showed a high biosynthetic potential related to antimicrobial and biotechnological classes and to organism defense and interactions as well as adaptive strategies for micronutrient uptake and cellular detoxification. Our results highlighted that bacteria living in an impacted environment, such as the surface sediments of the Venice Lagoon, may benefit from metabolic plasticity as well as from the synthesis of a wide array of secondary metabolites, promoting ecosystem resilience and stability toward environmental pressures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00192-z.

Flexible genomic island conservation across freshwater and marine Methylophilaceae

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Large-scale single-virus genomics uncovers hidden diversity of river water viruses and diversified gene profiles

Environmental viruses (primarily bacteriophages) are widely recognized as playing an important role in ecosystem homeostasis through the infection of host cells. However, the majority of environmental viruses are still unknown as their mosaic structure and frequent mutations in their sequences hinder genome construction in current metagenomics. To enable the large-scale acquisition of environmental viral genomes, we developed a new single-viral genome sequencing platform with microfluidic-generated gel beads. Amplification of individual DNA viral genomes in mass-produced gel beads allows high-throughput genome sequencing compared to conventional single-virus genomics. The sequencing analysis of river water samples yielded 1431 diverse viral single-amplified genomes, whereas viral metagenomics recovered 100 viral metagenome-assembled genomes at the comparable sequence depth. The 99.5% of viral single-amplified genomes were determined novel at the species level, most of which could not be recovered by a metagenomic assembly. The large-scale acquisition of diverse viral genomes identified protein clusters commonly detected in different viral strains, allowing the gene transfer to be tracked. Moreover, comparative genomics within the same viral species revealed that the profiles of various methyltransferase subtypes were diverse, suggesting an enhanced escape from host bacterial internal defense mechanisms. Our use of gel bead-based single-virus genomics will contribute to exploring the nature of viruses by accelerating the accumulation of draft genomes of environmental DNA viruses.

Spatiotemporal dynamics of giant viruses within a deep freshwater lake reveal a distinct dark-water community

Giant viruses (GVs) significantly regulate the ecological dynamics of diverse ecosystems. Although metagenomics has expanded our understanding of their diversity and ecological roles played in marine environments, little is known about GVs of freshwater ecosystems. Most previous studies have employed short-read sequencing and therefore resulted in fragmented genomes, hampering accurate assessment of genetic diversity. We sought to bridge this knowledge gap and overcome previous technical limitations. We subjected spatiotemporal (2 depths × 12 months) samples from Lake Biwa to metagenome-assembled genome reconstruction enhanced by long-read metagenomics. This yielded 293 GV metagenome-assembled genomes. Of these, 285 included previously unknown species in five orders of nucleocytoviruses and the first representatives of freshwater mirusviruses, which exhibited marked divergence from marine-derived lineages. The good performance of our long-read metagenomic assembly was demonstrated by the detection of 42 (14.3%) genomes composed of single contigs with completeness values >90%. GVs were partitioned across water depths, with most species specific to either the sunlit epilimnion or the dark hypolimnion. Epilimnion-specific members tended to be transient and exhibit short and intense abundance peaks, in line with the fact that they regulate the surface algal blooms. During the spring bloom, mirusviruses and members of three nucleocytovirus families were among the most abundant viruses. In contrast, hypolimnion-specific ones, including a mirusvirus genome, were typically more persistent in the hypolimnion throughout the water-stratified period, suggesting that they infect hosts specific to the hypolimnion and play previously unexplored ecological roles in dark water microbial ecosystems.

Ubiquitous genome streamlined Acidobacteriota in freshwater environments

Acidobacteriota are abundant in soil, peatlands, and sediments, but their ecology in freshwater environments remains understudied. UBA12189, an Acidobacteriota genus, is an uncultivated, genome-streamlined lineage with a small genome size found in aquatic environments where detailed genomic analyses are lacking. Here, we analyzed 66 MAGs of UBA12189 (including one complete genome) from freshwater lakes and rivers in Europe, North America, and Asia. UBA12189 has small genome sizes (<1.4 Mbp), low GC content, and a highly diverse pangenome. In freshwater lakes, this bacterial lineage is abundant from the surface waters (epilimnion) down to a 300-m depth (hypolimnion). UBA12189 appears to be free-living from CARD-FISH analysis. When compared to other genome-streamlined bacteria such as Nanopelagicales and Methylopumilus, genome reduction has caused UBA12189 to have a more limited metabolic repertoire in carbon, sulfur, and nitrogen metabolisms, limited numbers of membrane transporters, as well as a higher degree of auxotrophy for various amino acids, vitamins, and reduced sulfur. Despite having reduced genomes, UBA12189 encodes proteorhodopsin, complete biosynthesis pathways for heme and vitamin K2, cbb3-type cytochrome c oxidases, and heme-requiring enzymes. These genes may give a selective advantage during the genome streamlining process. We propose the new genus Acidiparvus, with two new species named "A. lacustris" and "A. fluvialis". Acidiparvus is the first described genome-streamlined lineage under the phylum Acidobacteriota, which is a free-living, slow-growing scavenger in freshwater environments.

A genome catalogue of lake bacterial diversity and its drivers at continental scale

Lakes are heterogeneous ecosystems inhabited by a rich microbiome whose genomic diversity is poorly defined. We present a continental-scale study of metagenomes representing 6.5 million km2 of the most lake-rich landscape on Earth. Analysis of 308 Canadian lakes resulted in a metagenome-assembled genome (MAG) catalogue of 1,008 mostly novel bacterial genomospecies. Lake trophic state was a leading driver of taxonomic and functional diversity among MAG assemblages, reflecting the responses of communities profiled by 16S rRNA amplicons and gene-centric metagenomics. Coupling the MAG catalogue with watershed geomatics revealed terrestrial influences of soils and land use on assemblages. Agriculture and human population density were drivers of turnover, indicating detectable anthropogenic imprints on lake bacteria at the continental scale. The sensitivity of bacterial assemblages to human impact reinforces lakes as sentinels of environmental change. Overall, the LakePulse MAG catalogue greatly expands the freshwater genomic landscape, advancing an integrative view of diversity across Earth's microbiomes.

In the right place, at the right time: the integration of bacteria into the Plankton Ecology Group model

BACKGROUND

Planktonic microbial communities have critical impacts on the pelagic food web and water quality status in freshwater ecosystems, yet no general model of bacterial community assembly linked to higher trophic levels and hydrodynamics has been assessed. In this study, we utilized a 2-year survey of planktonic communities from bacteria to zooplankton in three freshwater reservoirs to investigate their spatiotemporal dynamics.

RESULTS

We observed site-specific occurrence and microdiversification of bacteria in lacustrine and riverine environments, as well as in deep hypolimnia. Moreover, we determined recurrent bacterial seasonal patterns driven by both biotic and abiotic conditions, which could be integrated into the well-known Plankton Ecology Group (PEG) model describing primarily the seasonalities of larger plankton groups. Importantly, bacteria with different ecological potentials showed finely coordinated successions affiliated with four seasonal phases, including the spring bloom dominated by fast-growing opportunists, the clear-water phase associated with oligotrophic ultramicrobacteria, the summer phase characterized by phytoplankton bloom-associated bacteria, and the fall/winter phase driven by decay-specialists.

CONCLUSIONS

Our findings elucidate the major principles driving the spatiotemporal microbial community distribution in freshwater ecosystems. We suggest an extension to the original PEG model by integrating new findings on recurrent bacterial seasonal trends. Video Abstract.

Adaptive genetic traits in pelagic freshwater microbes

Pelagic microbes have adopted distinct strategies to inhabit the pelagial of lakes and oceans and can be broadly categorized in two groups: free-living, specialized oligotrophs and patch-associated generalists or copiotrophs. In this review, we aim to identify genomic traits that enable pelagic freshwater microbes to thrive in their habitat. To do so, we discuss the main genetic differences of pelagic marine and freshwater microbes that are both dominated by specialized oligotrophs and the difference to freshwater sediment microbes, where copiotrophs are more prevalent. We phylogenomically analysed a collection of >7700 metagenome-assembled genomes, classified habitat preferences on different taxonomic levels, and compared the metabolic traits of pelagic freshwater, marine, and freshwater sediment microbes. Metabolic differences are mainly associated with transport functions, environmental information processing, components of the electron transport chain, osmoregulation and the isoelectric point of proteins. Several lineages with known habitat transitions (Nitrososphaeria, SAR11, Methylophilaceae, Synechococcales, Flavobacteriaceae, Planctomycetota) and the underlying mechanisms in this process are discussed in this review. Additionally, the distribution, ecology and genomic make-up of the most abundant freshwater prokaryotes are described in details in separate chapters for Actinobacteriota, Bacteroidota, Burkholderiales, Verrucomicrobiota, Chloroflexota, and 'Ca. Patescibacteria'.

Evolutionary ecology of microbial populations inhabiting deep sea sediments associated with cold seeps

Deep sea cold seep sediments host abundant and diverse microbial populations that significantly influence biogeochemical cycles. While numerous studies have revealed their community structure and functional capabilities, little is known about genetic heterogeneity within species. Here, we examine intraspecies diversity patterns of 39 abundant species identified in sediment layers down to 430 cm below the sea floor across six cold seep sites. These populations are grouped as aerobic methane-oxidizing bacteria, anaerobic methanotrophic archaea and sulfate-reducing bacteria. Different evolutionary trajectories are observed at the genomic level among these physiologically and phylogenetically diverse populations, with generally low rates of homologous recombination and strong purifying selection. Functional genes related to methane (pmoA and mcrA) and sulfate (dsrA) metabolisms are under strong purifying selection in most species investigated. These genes differ in evolutionary trajectories across phylogenetic clades but are functionally conserved across sites. Intrapopulation diversification of genomes and their mcrA and dsrA genes is depth-dependent and subject to different selection pressure throughout the sediment column redox zones at different sites. These results highlight the interplay between ecological processes and the evolution of key bacteria and archaea in deep sea cold seep extreme environments, shedding light on microbial adaptation in the subseafloor biosphere.

A Fast and Easy Method to Co-extract DNA and RNA from an Environmental Microbial Sample

We herein propose a fast and easy DNA and RNA co-extraction method for environmental microbial samples. It combines bead beating and phenol-chloroform phase separation followed by the separation and purification of DNA and RNA using the Qiagen AllPrep DNA/RNA mini kit. With a handling time of ~3 h, our method simultaneously extracted high-quality DNA (peak size >10-15‍ ‍kb) and RNA (RNA integrity number >6) from lake bacterioplankton filtered samples. The method is also applicable to low-biomass samples (expected DNA or RNA yield <50‍ ‍ng) and eukaryotic microbial samples, providing an easy option for more versatile eco-genomic applications.