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.

Factors shaping community patterns of protists and bacteria on a European scale.

Factors shaping community patterns of microorganisms are controversially discussed. Physical and chemical factors certainly limit the survival of individual taxa and maintenance of diversity. In recent years, a contribution of geographic distance and dispersal barriers to distribution patterns of protists and bacteria has been demonstrated. Organismic interactions such as competition, predation and mutualism further modify community structure and maintenance of distinct taxa. Here, we address the relative importance of these different factors in shaping protists and bacterial communities on a European scale using high-throughput sequencing data obtained from lentic freshwater ecosystems. We show that community patterns of protists are similar to those of bacteria. Our results indicate that cross-domain organismic factors are important variables with a higher influence on protists as compared with bacteria. Abiotic physical and chemical factors also contributed significantly to community patterns. The contribution of these latter factors was higher for bacteria, which may reflect a stronger biogeochemical coupling. The contribution of geographical distance was similar for both microbial groups.

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

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

Differential associations of five riverine organism groups with multiple stressors.

The decline of river and stream biodiversity results from multiple simultaneous occuring stressors, yet few studies explore responses explore responses across various taxonomic groups at the same locations. In this study, we address this shortcoming by using a coherent data set to study the association of nine commonly occurring stressors (five chemical, one morphological and three hydraulic) with five taxonomic groups (bacteria, fungi, diatoms, macro-invertebrates and fish). According to studies on single taxonomic groups, we hypothesise that gradients of chemical stressors structure community composition of all taxonomic groups, while gradients of hydraulic and morphological stressors are mainly related to larger organisms such as benthic macro-invertebrates and fish. Organisms were sampled over two years at 20 sites in two catchments: a recently restored urban lowland catchment (Boye) and a moderately disturbed rural mountainous catchment (Kinzig). Dissimilarity matrices were computed for each taxonomic group within a catchment. Taxonomic dissimilarities between sites were linked to stressor dissimilarities using multivariable Generalized Linear Mixed Models. Stressor gradients were longer in the Boye, but did in contrast to the Kinzig not cover low stress intensities. Accordingly, responses of the taxonomic groups were stronger in the Kinzig catchment than in the recently restored Boye catchment. The discrepancy between catchments underlines that associations to stressors strongly depend on which part of the stressor gradient is covered in a catchment. All taxonomic groups were related to conductivity. Bacteria, fungi and macro-invertebrates change with dissolved oxygen, and bacteria and fungi with total nitrogen. Morphological and hydraulic stressors had minor correlations with bacteria, fungi and diatoms, while macro-invertebrates were strongly related to fine sediment and discharge, and fish to high flow peaks. The results partly support our hypotheses about the differential associations of the different taxonomic groups with the stressors.

A multi-omics study on quantifying antimicrobial resistance in European freshwater lakes.

The surveillance of wastewater for the Covid-19 virus during this unprecedented pandemic and mapped to the distribution and magnitude of the infected in the population near real-time exemplifies the importance of tracking rapidly changing trends of pathogens or public health problems at a large scale. The rising trends of antimicrobial resistance (AMR) with multidrug-resistant pathogens from the environmental water have similarly gained much attention in recent years. Wastewater-based epidemiology from water samples has shown that a wide range of AMR-related genes is frequently detected. Albeit sewage is treated before release and thus, the abundance of pathogens should be significantly reduced or even pathogen-free, several studies indicated the contrary. Pathogens are still measurable in the released water, ultimately entering freshwaters, such as rivers and lakes. Furthermore, socio-economic and environmental factors, such as chemical industries and animal farming nearby, impact the presence of AMR. Many bacterial species from the environment are intrinsically resistant and also contribute to the resistome of freshwater lakes. This study collected the most extensive standardized freshwater data set from hundreds of European lakes and conducted a comprehensive multi-omics analysis on antimicrobial resistance from these freshwater lakes. Our research shows that genes encoding for AMR against tetracyclines, cephalosporins, and quinolones were commonly identified, while for some, such as sulfonamides, resistance was less frequently present. We provide an estimation of the characteristic resistance of AMR in European lakes, which can be used as a comprehensive resistome dataset to facilitate and monitor temporal changes in the development of AMR in European freshwater lakes.

A Genomic Perspective Across Earth's Microbiomes Reveals That Genome Size in Archaea and Bacteria Is Linked to Ecosystem Type and Trophic Strategy.

Our view of genome size in Archaea and Bacteria has remained skewed as the data has been dominated by genomes of microorganisms that have been cultivated under laboratory settings. However, the continuous effort to catalog Earth's microbiomes, specifically propelled by recent extensive work on uncultivated microorganisms, provides an opportunity to revise our perspective on genome size distribution. We present a meta-analysis that includes 26,101 representative genomes from 3 published genomic databases; metagenomic assembled genomes (MAGs) from GEMs and stratfreshDB, and isolates from GTDB. Aquatic and host-associated microbial genomes present on average the smallest estimated genome sizes (3.1 and 3.0 Mbp, respectively). These are followed by terrestrial microbial genomes (average 3.7 Mbp), and genomes from isolated microorganisms (average 4.3 Mbp). On the one hand, aquatic and host-associated ecosystems present smaller genomes sizes in genera of phyla with genome sizes above 3 Mbp. On the other hand, estimated genome size in phyla with genomes under 3 Mbp showed no difference between ecosystems. Moreover, we observed that when using 95% average nucleotide identity (ANI) as an estimator for genetic units, only 3% of MAGs cluster together with genomes from isolated microorganisms. Although there are potential methodological limitations when assembling and binning MAGs, we found that in genome clusters containing both environmental MAGs and isolate genomes, MAGs were estimated only an average 3.7% smaller than isolate genomes. Even when assembly and binning methods introduce biases, estimated genome size of MAGs and isolates are very similar. Finally, to better understand the ecological drivers of genome size, we discuss on the known and the overlooked factors that influence genome size in different ecosystems, phylogenetic groups, and trophic strategies.

Ecological Differentiation in Two Major Freshwater Bacterial Taxa Along Environmental Gradients.

Polynucleobacter (Burkholderiaceae, Betaproteobacteria) and Limnohabitans (Comamonadaceae, Betaproteobacteria) are abundant freshwater bacteria comprising large genetic and taxonomic diversities, with species adapted to physico-chemically distinct types of freshwater systems. The relative importance of environmental drivers, i.e., physico-chemistry, presence of microeukaryotes and geographic position for the diversity and prevalence has not been investigated for both taxa before. Here, we present the first pan-European study on this topic, comprising 255 freshwater lakes. We investigated Limnohabitans and Polynucleobacter using an amplicon sequencing approach of partial 16S rRNA genes along environmental gradients. We show that physico-chemical factors had the greatest impact on both genera. Analyses on environmental gradients revealed an exceptionally broad ecological spectrum of operational taxonomic units (OTUs). Despite the coarse resolution of the genetic marker, we found OTUs with contrasting environmental preferences within Polynucleobacter and Limnohabitans subclusters. Such an ecological differentiation has been characterized for PnecC and LimC before but was so far unknown for less well studied subclusters such as PnecA and PnecB. Richness and abundance of OTUs are geographically clustered, suggesting that geographic diversity patterns are attributable to region-specific physico-chemical characteristics (e.g., pH and temperature) rather than latitudinal gradients or lake sizes.

Cryptophyta as major bacterivores in freshwater summer plankton.

Small bacterivorous eukaryotes play a cardinal role in aquatic food webs and their taxonomic classification is currently a hot topic in aquatic microbial ecology. Despite increasing interest in their diversity, core questions regarding predator-prey specificity remain largely unanswered, e.g., which heterotrophic nanoflagellates (HNFs) are the main bacterivores in freshwaters and which prokaryotes support the growth of small HNFs. To answer these questions, we fed natural communities of HNFs from Rímov reservoir (Czech Republic) with five different bacterial strains of the ubiquitous betaproteobacterial genera Polynucleobacter and Limnohabitans. We combined amplicon sequencing and catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH) targeting eukaryotic 18 S rRNA genes to track specific responses of the natural HNF community to prey amendments. While amplicon sequencing provided valuable qualitative data and a basis for designing specific probes, the number of reads was insufficient to accurately quantify certain eukaryotic groups. We also applied a double-hybridization technique that allows simultaneous phylogenetic identification of both predator and prey. Our results show that community composition of HNFs is strongly dependent upon prey type. Surprisingly, Cryptophyta were the most abundant bacterivores, although this phylum has been so far assumed to be mainly autotrophic. Moreover, the growth of a small lineage of Cryptophyta (CRY1 clade) was strongly stimulated by one Limnohabitans strain in our experiment. Thus, our study is the first report that colorless Cryptophyta are major bacterivores in summer plankton samples and can play a key role in the carbon transfer from prokaryotes to higher trophic levels.

Responses of stream microbes to multiple anthropogenic stressors in a mesocosm study.

Stream ecosystems are affected by multiple anthropogenic stressors worldwide. Even though effects of many single stressors are comparatively well studied, the effects of multiple stressors are difficult to predict. In particular bacteria and protists, which are responsible for the majority of ecosystem respiration and element flows, are infrequently studied with respect to multiple stressors responses. We conducted a stream mesocosm experiment to characterize the responses of single and multiple stressors on microbiota. Two functionally important stream habitats, leaf litter and benthic phototrophic rock biofilms, were exposed to three stressors in a full factorial design: fine sediment deposition, increased chloride concentration (salinization) and reduced flow velocity. We analyzed the microbial composition in the two habitat types of the mesocosms using an amplicon sequencing approach. Community analysis on different taxonomic levels as well as principle component analyses (PCoAs) based on realtive abundances of operational taxonomic units (OTUs) showed treatment specific shifts in the eukaryotic biofilm community. Analysis of variance (ANOVA) revealed that Bacillariophyta responded positively salinity and sediment increase, while the relative read abundance of chlorophyte taxa decreased. The combined effects of multiple stressors were mainly antagonistic. Therefore, the community composition in multiply stressed environments resembled the composition of the unstressed control community in terms of OTU occurrence and relative abundances.