Ecogenomics sheds light on diverse lifestyle strategies in freshwater CPR

Genome streamlining in Parcubacteria transitioning from soil to groundwater

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-resolved metagenomics reveals the effect of nutrient availability on bacterial genomic properties across 44 European freshwater lakes

Understanding intricate microbial interactions in the environment is crucial. This is especially true for the relationships between nutrients and bacteria, as phosphorus, nitrogen and organic carbon availability are known to influence bacterial population dynamics. It has been suggested that low nutrient conditions prompt the evolutionary process of genome streamlining. This process helps conserve scarce nutrients and allows for proliferation. Genome streamlining is associated with genomic properties such as %GC content, genes encoding sigma factors, percent coding regions, gene redundancy, and functional shifts in processes like cell motility and ATP binding cassette transporters, among others. The current study aims to unveil the impact of nutrition on the genome size, %GC content, and functional properties of pelagic freshwater bacteria. We do this at finer taxonomic resolutions for many metagenomically characterized communities. Our study confirms the interplay of trophic level and genomic properties. It also highlights that different nutrient types, particularly phosphorus and nitrogen, impact these properties differently. We observed a covariation of functional traits with genome size. Larger genomes exhibit enriched pathways for motility, environmental interaction, and regulatory genes. ABC transporter genes reflect the availability of nutrients in the environment, with small genomes presumably relying more on metabolites from other organisms. We also discuss the distinct strategies different phyla adopt to adapt to oligotrophic environments. The findings contribute to our understanding of genomic adaptations within complex microbial communities.

Global abundance patterns, diversity, and ecology of Patescibacteria in wastewater treatment plants

BACKGROUND

Microorganisms are responsible for nutrient removal and resource recovery in wastewater treatment plants (WWTPs), and their diversity is often studied by 16S rRNA gene amplicon sequencing. However, this approach underestimates the abundance and diversity of Patescibacteria due to the low coverage of commonly used PCR primers for this highly divergent bacterial phylum. Therefore, our current understanding of the global diversity, distribution, and ecological role of Patescibacteria in WWTPs is very incomplete. This is particularly relevant as Patescibacteria are considered to be associated with microbial host cells and can therefore influence the abundance and temporal variability of other microbial groups that are important for WWTP functioning.

RESULTS

Here, we evaluated the in silico coverage of widely used 16S rRNA gene-targeted primer pairs and redesigned a primer pair targeting the V4 region of bacterial and archaeal 16S rRNA genes to expand its coverage for Patescibacteria. We then experimentally evaluated and compared the performance of the original and modified V4-targeted primers on 565 WWTP samples from the MiDAS global sample collection. Using the modified primer pair, the percentage of ASVs classified as Patescibacteria increased from 5.9 to 23.8%, and the number of detected patescibacterial genera increased from 560 to 1576, while the detected diversity of the remaining microbial community remained similar. Due to this significantly improved coverage of Patescibacteria, we identified 23 core genera of Patescibacteria in WWTPs and described the global distribution pattern of these unusual microbes in these systems. Finally, correlation network analysis revealed potential host organisms that might be associated with Patescibacteria in WWTPs. Interestingly, strong indications were found for an association between Patescibacteria of the Saccharimonadia and globally abundant polyphosphate-accumulating organisms of the genus Ca. Phosphoribacter.

CONCLUSIONS

Our study (i) provides an improved 16S rRNA gene V4 region-targeted amplicon primer pair inclusive of Patescibacteria with little impact on the detection of other taxa, (ii) reveals the diversity and distribution patterns of Patescibacteria in WWTPs on a global scale, and (iii) provides new insights into the ecological role and potential hosts of Patescibacteria in WWTPs. Video Abstract.

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.

Enrichment of different taxa of the enigmatic candidate phyla radiation bacteria using a novel picolitre droplet technique

The candidate phyla radiation (CPR) represents a distinct monophyletic clade and constitutes a major portion of the tree of life. Extensive efforts have focused on deciphering the functional diversity of its members, primarily using sequencing-based techniques. However, cultivation success remains scarce, presenting a significant challenge, particularly in CPR-dominated groundwater microbiomes characterized by low biomass. Here, we employ an advanced high-throughput droplet microfluidics technique to enrich CPR taxa from groundwater. Utilizing a low-volume filtration approach, we successfully harvested a microbiome resembling the original groundwater microbial community. We assessed CPR enrichment in droplet and aqueous bulk cultivation for 30 days using a novel CPR-specific primer to rapidly track the CPR fraction through the cultivation attempts. The combination of soil extract and microbial-derived necromass provided the most supportive conditions for CPR enrichment. Employing these supplemented conditions, droplet cultivation proved superior to bulk cultivation, resulting in up to a 13-fold CPR enrichment compared to a 1- to 2-fold increase in bulk cultivation. Amplicon sequencing revealed 10 significantly enriched CPR orders. The highest enrichment in CPRs was observed for some unknown members of the Parcubacteria order, Cand. Jorgensenbacteria, and unclassified UBA9983. Furthermore, we identified co-enriched putative host taxa, which may guide more targeted CPR isolation approaches in subsequent investigations.

First shotgun metagenomics study of Juan de Fuca deep-sea sediments reveals distinct microbial communities above, within, between, and below sulfate methane transition zones

The marine deep subsurface is home to a vast microbial ecosystem, affecting biogeochemical cycles on a global scale. One of the better-studied deep biospheres is the Juan de Fuca (JdF) Ridge, where hydrothermal fluid introduces oxidants into the sediment from below, resulting in two sulfate methane transition zones (SMTZs). In this study, we present the first shotgun metagenomics study of unamplified DNA from sediment samples from different depths in this stratified environment. Bioinformatic analyses showed a shift from a heterotrophic, Chloroflexota-dominated community above the upper SMTZ to a chemolithoautotrophic Proteobacteria-dominated community below the secondary SMTZ. The reintroduction of sulfate likely enables respiration and boosts active cells that oxidize acetate, iron, and complex carbohydrates to degrade dead biomass in this low-abundance, low-diversity environment. In addition, analyses showed many proteins of unknown function as well as novel metagenome-assembled genomes (MAGs). The study provides new insights into microbial communities in this habitat, enabled by an improved DNA extraction protocol that allows a less biased view of taxonomic composition and metabolic activities, as well as uncovering novel taxa. Our approach presents the first successful attempt at unamplified shotgun sequencing samples from beyond 50 meters below the seafloor and opens new ways for capturing the true diversity and functional potential of deep-sea sediments.

Influence of soil nutrients on the presence and distribution of CPR bacteria in a long-term crop rotation experiment

Bacteria affiliated with the Candidate Phyla Radiation (CPR) are a hyper-diverse group of ultra-small bacteria with versatile yet sparse metabolisms. However, most insights into this group come from a surprisingly small number of environments, and recovery of CPR bacteria from soils has been hindered due to their extremely low abundance within complex microbial assemblages. In this study we enriched soil samples from 14 different soil fertility treatments for ultra-small (<0.45 μm) bacteria in order to study rare soil CPR. 42 samples were sequenced, enabling the reconstruction of 27 quality CPR metagenome-assembled genomes (MAGs) further classified as Parcubacteria/Paceibacteria, Saccharibacteria/Saccharimonadia and ABY1, in addition to representative genomes from Gemmatimonadetes, Dependentiae and Chlamydae phyla. These genomes were fully annotated and used to reconstruct the CPR community across all 14 plots. Additionally, for five of these plots, the entire microbiota was reconstructed using 16S amplification, showing that specific soil CPR may form symbiotic relationships with a varied and circumstantial range of hosts. Cullars CPR had a prevalence of enzymes predicted to degrade plant-derived carbohydrates, which suggests they have a role in plant biomass degradation. Parcubacteria appear to be more apt at microfauna necromass degradation. Cullars Saccharibacteria and a Parcubacteria group were shown to carry a possible aerotolerance mechanism coupled with potential for aerobic respiration, which appear to be a unique adaptation to the oxic soil environment. Reconstruction of CPR communities across treatment plots showed that they were not impacted by changes in nutrient levels or microbiota composition, being only impacted by extreme conditions, causing some CPR to dominate the community. These findings corroborate the understanding that soil-dwelling CPR bacteria have a very broad symbiont range and have metabolic capabilities associated to soil environments which allows them to scavenge resources and form resilient communities. The contributions of these microbial dark matter species to soil ecology and plant interactions will be of significant interest in future studies.

Improvement of carbon source composition reduces antibiotic resistance genes in the ectopic fermentation system

Effectively reduce antibiotic resistance genes (ARGs) in ectopic fermentation system (EFS) is essential for practical production. In this study, three experiments were performed to explore how to remove ARGs in EFS effectively. Results demonstrated that ARGs were easily enriched in rice-husk-sawdust padding; simultaneous addition of laccase and cellulase suppressed the ARGs, mainly by increasing soluble carbohydrate concentration and promoting humic acid concentration; addition of corn stalks into rice-husk-sawdust decreased the abundance of ARGs by improving the carbon source structure and enhancing cellulase activity. In conclusion, the present study provides a guidance to reduce the threat of ARGs in EFS, which paved a potential pathway to safely utilize manure resources.

Using Genomics to Decipher the Enigmatic Properties and Survival Adaptation of Candidate Phyla Radiation

Microbial ecology is a critical field for understanding the composition, diversity, and functions of microorganisms in various environmental and health-related processes. The discovery of Candidate Phyla Radiation (CPR) through culture-independent methods has introduced a new division of microbes characterized by a symbiotic/parasitic lifestyle, small cell size, and small genome. Despite being poorly understood, CPRs have garnered significant attention in recent years due to their widespread detection in a variety of environmental and clinical samples. These microorganisms have been found to exhibit a high degree of genetic diversity compared to other microbes. Several studies have shed light on their potential importance in global biogeochemical cycles and their impact on various human activities. In this review, we provide a systematic overview of the discovery of CPRs. We then focus on describing how the genomic characteristics of CPRs have helped them interact with and adapt to other microbes in different ecological niches. Future works should focus on discovering the metabolic capacities of CPRs and, if possible, isolating them to obtain a better understanding of these microorganisms.

Baseline metagenome-assembled genome (MAG) data of Sikkim hot springs from Indian Himalayan geothermal belt (IHGB) showcasing its potential CAZymes, and sulfur-nitrogen metabolic activity

Here we present the construction and characterization of metagenome assembled genomes (MAGs) from two hot springs residing in the vicinity of Indian Himalayan Geothermal Belt (IHGB). A total of 78 and 7 taxonomic bins were obtained for Old Yume Samdong (OYS) and New Yume Samdong (NYS) hot springs respectively. After passing all the criteria only 21 and 4 MAGs were further studied based on the successful prediction of their 16 S rRNA. Various databases were used such as GTDB, Kaiju, EzTaxon, BLAST XY Plot and NCBI BLAST to get the taxonomic classification of various 16 S rRNA predicted MAGs. The bacterial genomes found were from both thermophilic and mesophilic bacteria among which Proteobacteria, Chloroflexi, Bacteroidetes and Firmicutes were the abundant phyla. However, in case of OYS, two genomes belonged to archaeal Methanobacterium and Methanocaldococcus. Functional characterization revealed the richness of CAZymes such as Glycosyl Transferase (GT) (56.7%), Glycoside Hydrolase (GH) (37.4%), Carbohydrate Esterase family (CE) (8.2%), and Polysaccharide Lyase (PL) (1.9%). There were negligible antibiotic resistance genes in the MAGs however, a significant heavy metal tolerance gene was found in the MAGs. Thus, it may be assumed that there is no coexistence of antibiotic and heavy metal resistance genes in these hot spring microbiomes. Since the selected hot springs possess good sulfur content thus, we also checked the presence of genes for sulfur and nitrogen metabolism. It was found that MAGs from both the hot springs possess significant number of genes related to sulfur and nitrogen metabolism.

Genome-centric metagenomics reveals the host-driven dynamics and ecological role of CPR bacteria in an activated sludge system

BACKGROUND

Candidate phyla radiation (CPR) constitutes highly diverse bacteria with small cell sizes and are likely obligate intracellular symbionts. Given their distribution and complex associations with bacterial hosts, genetic and biological features of CPR bacteria in low-nutrient environments have received increasing attention. However, CPR bacteria in wastewater treatment systems remain poorly understood. We utilized genome-centric metagenomics to answer how CPR communities shift over 11 years and what kind of ecological roles they act in an activated sludge system.

RESULTS

We found that approximately 9% (135) of the 1,526 non-redundant bacterial and archaeal metagenome-assembled genomes were affiliated with CPR. CPR bacteria were consistently abundant with a relative abundance of up to 7.5% in the studied activated sludge system. The observed striking fluctuations in CPR community compositions and the limited metabolic and biosynthetic capabilities in CPR bacteria collectively revealed the nature that CPR dynamics may be directly determined by the available hosts. Similarity-based network analysis further confirmed the broad bacterial hosts of CPR lineages. The proteome contents of activated sludge-associated CPR had a higher similarity to those of environmental-associated CPR than to those of human-associated ones. Comparative genomic analysis observed significant enrichment of genes for oxygen stress resistance in activated sludge-associated CPR bacteria. Furthermore, genes for carbon cycling and horizontal gene transfer were extensively identified in activated sludge-associated CPR genomes.

CONCLUSIONS

These findings highlight the presence of specific host interactions among CPR lineages in activated sludge systems. Despite the lack of key metabolic pathways, these small, yet abundant bacteria may have significant involvements in biogeochemical cycling and bacterial evolution in activated sludge systems. 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'.

A novel and diverse group of Candidatus Patescibacteria from bathypelagic Lake Baikal revealed through long-read metagenomics

BACKGROUND

Lake Baikal, the world's deepest freshwater lake, contains important numbers of Candidatus Patescibacteria (formerly CPR) in its deepest reaches. However, previously obtained CPR metagenome-assembled genomes recruited very poorly indicating the potential of other groups being present. Here, we have applied for the first time a long-read (PacBio CCS) metagenomic approach to analyze in depth the Ca. Patescibacteria living in the bathypelagic water column of Lake Baikal at 1600 m.

RESULTS

The retrieval of nearly complete 16S rRNA genes before assembly has allowed us to detect the presence of a novel and a likely endemic group of Ca. Patescibacteria inhabiting bathypelagic Lake Baikal. This novel group seems to possess extremely high intra-clade diversity, precluding complete genomes' assembly. However, read binning and scaffolding indicate that these microbes are similar to other Ca. Patescibacteria (i.e. parasites or symbionts), although they seem to carry more anabolic pathways, likely reflecting the extremely oligotrophic habitat they inhabit. The novel bins have not been found anywhere, but one of the groups appears in small amounts in an oligotrophic and deep alpine Lake Thun. We propose this novel group be named Baikalibacteria.

CONCLUSION

The recovery of 16S rRNA genes via long-read metagenomics plus the use of long-read binning to uncover highly diverse "hidden" groups of prokaryotes are key strategies to move forward in ecogenomic microbiology. The novel group possesses enormous intraclade diversity akin to what happens with Ca. Patescibacteria at the interclade level, which is remarkable in an environment that has changed little in the last 25 million years.

Single-amplified genomes reveal most streamlined free-living marine bacteria

Evolutionary adaptations of prokaryotes to the environment sometimes result in genome reduction. Our knowledge of this phenomenon among free-living bacteria remains scarce. We address the dynamics and limits of genome reduction by examining one of the most abundant bacteria in the ocean, the SAR86 clade. Despite its abundance, comparative genomics has been limited by the absence of pure cultures and the poor representation in metagenome-assembled genomes. We co-assembled multiple previously available single-amplified genomes to obtain the first complete genomes from members of the four families. All families showed a convergent evolutionary trajectory with characteristic features of streamlined genomes, most pronounced in the TMED112 family. This family has a genome size of ca. 1 Mb and only 1 bp as median intergenic distance, exceeding values found in other abundant microbes such as SAR11, OM43 and Prochlorococcus. This genomic simplification led to a reduction in the biosynthesis of essential molecules, DNA repair-related genes, and the ability to sense and respond to environmental factors, which could suggest an evolutionary dependence on other co-occurring microbes for survival (Black Queen hypothesis). Therefore, these reconstructed genomes within the SAR86 clade provide new insights into the limits of genome reduction in free-living marine bacteria.

Into the darkness of the microbial dark matter in situ activities through expression profiles of Patescibacteria populations

Patescibacteria form a highly diverse and widespread superphylum of uncultured microorganisms representing a third of the global microbial diversity. Most of our knowledge on Patescibacteria putative physiology relies on metagenomic mining and metagenome-assembled genomes, but the in situ activities and the ecophysiology of these microorganisms have been rarely explored, leaving the role of Patescibacteria in ecosystems elusive. Using a genome-centric metatranscriptomic approach, we analyzed the diel and seasonal gene transcription profiles of 18 Patescibacteria populations in brackish microbial mats to test whether our understanding of Patescibacteria metabolism allows the extrapolation of their in situ activities. Although our results revealed a circadian cycle in Patescibacteria activities, a strong streamlined genetic expression characterized the Patescibacteria populations. This result has a major consequence for the extrapolation of their physiology and environmental function since most transcribed genes were uncharacterized, indicating that the ecophysiology of Patescibacteria cannot be yet reliably predicted from genomic data.

Occurrence, Diversity, and Genomes of "Candidatus Patescibacteria" along the Early Diagenesis of Marine Sediments

The phylum "Candidatus Patescibacteria" (or Candidate Phyla Radiation [CPR]) accounts for roughly one-quarter of microbial diversity on Earth, but the presence and diversity of these bacteria in marine sediments have been rarely charted. Here, we investigate the abundance, diversity, and metabolic capacities of CPR bacteria in three sediment sites (Mohns Ridge, North Pond, and Costa Rica Margin) with samples covering a wide range of redox zones formed during the early diagenesis of organic matter. Through metagenome sequencing, we found that all investigated sediment horizons contain "Ca. Patescibacteria" (0.4 to 28% of the total communities), which are affiliated with the classes "Ca. Paceibacteria," "Ca. Gracilibacteria," "Ca. Microgenomatia," "Ca. Saccharimonadia," "Ca. ABY1," and "Ca. WWE3." However, only a subset of the diversity of marine sediment "Ca. Patescibacteria," especially the classes "Ca. Paceibacteria" and "Ca. Gracilibacteria," can be captured by 16S rRNA gene amplicon sequencing with commonly used universal primers. We recovered 11 metagenome-assembled genomes (MAGs) of CPR from these sediments, most of which are novel at the family or genus level in the "Ca. Paceibacteria" class and are missed by the amplicon sequencing. While individual MAGs are confined to specific anoxic niches, the lack of capacities to utilize the prevailing terminal electron acceptors indicates that they may not be directly selected by the local redox conditions. These CPR bacteria lack essential biosynthesis pathways and may use a truncated glycolysis pathway to conserve energy as fermentative organotrophs. Our findings suggest that marine sediments harbor some novel yet widespread CPR bacteria during the early diagenesis of organic matter, which needs to be considered in population dynamics assessments in this vast environment. IMPORTANCE Ultrasmall-celled "Ca. Patescibacteria" have been estimated to account for one-quarter of the total microbial diversity on Earth, the parasitic lifestyle of which may exert a profound control on the overall microbial population size of the local ecosystems. However, their diversity and metabolic functions in marine sediments, one of the largest yet understudied ecosystems on Earth, remain virtually uncharacterized. By applying cultivation-independent approaches to a range of sediment redox zones, we reveal that "Ca. Patescibacteria" members are rare but widespread regardless of the prevailing geochemical conditions. These bacteria are affiliated with novel branches of "Ca. Patescibacteria" and have been largely missed in marker gene-based surveys. They do not have respiration capacity but may conserve energy by fermenting organic compounds from their episymbiotic hosts. Our findings suggest that these novel "Ca. Patescibacteria" are among the previously overlooked microbes in diverse marine sediments.

Widespread but Poorly Understood Bacteria: Candidate Phyla Radiation

Candidate Phyla Radiation (CPR) bacteria is a bacterial division composed mainly of candidate phyla bacteria with ultra-small cell sizes, streamlined genomes, and limited metabolic capacity, which are generally considered to survive in a parasitic or symbiotic manner. Despite their wide distribution and rich diversity, CPR bacteria have received little attention until recent years, and are therefore poorly understood. This review systematically summarizes the history of CPR research, the parasitic/symbiotic lifestyle, and the ecological distribution and unique metabolic features of CPR bacteria, hoping to provide guidance for future ecological and physiological research on CPR bacteria.

The combined impact of low temperatures and shifting phosphorus availability on the competitive ability of cyanobacteria

In freshwater systems, cyanobacteria are strong competitors under enhanced temperature and eutrophic conditions. Understanding their adaptive and evolutionary potential to multiple environmental states allows us to accurately predict their response to future conditions. To better understand if the combined impacts of temperature and nutrient limitation could suppress the cyanobacterial blooms, a single strain of Microcystis aeruginosa was inoculated into natural phytoplankton communities with different nutrient conditions: oligotrophic, eutrophic and eutrophic with the addition of bentophos. We found that the use of the bentophos treatment causes significant differences in prokaryotic and eukaryotic communities. This resulted in reduced biodiversity among the eukaryotes and a decline in cyanobacterial abundance suggesting phosphorus limitation had a strong impact on the community structure. The low temperature during the experiment lead to the disappearance of M. aeruginosa in all treatments and gave other phytoplankton groups a competitive advantage leading to the dominance of the eukaryotic families that have diverse morphologies and nutritional modes. These results show cyanobacteria have a reduced competitive advantage under certain temperature and nutrient limiting conditions and therefore, controlling phosphorus concentrations could be a possible mitigation strategy for managing harmful cyanobacterial blooms in a future warmer climate.

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

Reconstruction of metagenome-assembled genomes (MAGs) has become a fundamental approach in microbial ecology. However, a MAG is hardly complete and overlooks genomic microdiversity because metagenomic assembly fails to resolve microvariants among closely related genotypes. Aiming at understanding the universal factors that drive or constrain prokaryotic genome diversification, we performed an ecosystem-wide high-resolution metagenomic exploration of microdiversity by combining spatiotemporal (2 depths × 12 months) sampling from a pelagic freshwater system, high-quality MAG reconstruction using long- and short-read metagenomic sequences, and profiling of single nucleotide variants (SNVs) and structural variants (SVs) through mapping of short and long reads to the MAGs, respectively. We reconstructed 575 MAGs, including 29 circular assemblies, providing high-quality reference genomes of freshwater bacterioplankton. Read mapping against these MAGs identified 100 to 101,781 SNVs/Mb and 0 to 305 insertions, 0 to 467 deletions, 0 to 41 duplications, and 0 to 6 inversions for each MAG. Nonsynonymous SNVs were accumulated in genes potentially involved in cell surface structural modification to evade phage recognition. Most (80.2%) deletions overlapped with a gene coding region, and genes of prokaryotic defense systems were most frequently (>8% of the genes) overlapped with a deletion. Some such deletions exhibited a monthly shift in their allele frequency, suggesting a rapid turnover of genotypes in response to phage predation. MAGs with extremely low microdiversity were either rare or opportunistic bloomers, suggesting that population persistency is key to their genomic diversification. The results concluded that prokaryotic genomic diversification is driven primarily by viral load and constrained by a population bottleneck. IMPORTANCE Identifying intraspecies genomic diversity (microdiversity) is crucial to understanding microbial ecology and evolution. However, microdiversity among environmental assemblages is not well investigated, because most microbes are difficult to culture. In this study, we performed cultivation-independent exploration of bacterial genomic microdiversity in a lake ecosystem using a combination of short- and long-read metagenomic analyses. The results revealed the broad spectrum of genomic microdiversity among the diverse bacterial species in the ecosystem, which has been overlooked by conventional approaches. Our ecosystem-wide exploration further allowed comparative analysis among the genomes and genes and revealed factors behind microbial genomic diversification, namely, that diversification is driven primarily by resistance against viral infection and constrained by the population size.