Seasonal and episodic lake mixing stimulate differential planktonic bacterial dynamics

Metagenomic analysis reveals gene taxonomic and functional diversity response to microplastics and cadmium in an agricultural soil

Both microplastics (MPs) and heavy metals are common soil pollutants and can interact to generate combined toxicity to soil ecosystems, but their impact on soil microbial communities (e.g., archaea and viruses) remains poorly studied. Here, metagenomic analysis was used to explore the response of soil microbiome in an agricultural soil exposed to MPs [i.e., polyethylene (PE), polystyrene (PS), and polylactic acid (PLA)] and/or Cd. Results showed that MPs had more profound effects on microbial community composition, diversity, and gene abundances when compared to Cd or their combination. Metagenomic analysis indicated that the gene taxonomic diversity and functional diversity of microbial communities varied with MPs type and dose. MPs affected the relative abundance of major microbial phyla and genera, while their coexistence with Cd influenced dominant fungi and viruses. Nitrogen-transforming and pathogenic genera, which were more sensitive to MPs variations, could serve as the indicative taxa for MPs contamination. High-dose PLA treatments (10%, w/w) not only elevated nitrogen metabolism and pathogenic genes, but also enriched copiotrophic microbes from the Proteobacteria phylum. Overall, MPs and Cd showed minimal interactions on soil microbial communities. This study highlights the microbial shifts due to co-occurring MPs and Cd, providing evidence for understanding their environmental risks.

Aqueous Geochemical Controls on the Sestonic Microbial Community in Lakes Michigan and Superior

Despite being the largest freshwater lake system in the world, relatively little is known about the sestonic microbial community structure in the Laurentian Great Lakes. The goal of this research was to better understand this ecosystem using high-throughput sequencing of microbial communities as a function of water depth at six locations in the westernmost Great Lakes of Superior and Michigan. The water column was characterized by gradients in temperature, dissolved oxygen (DO), pH, and other physicochemical parameters with depth. Mean nitrate concentrations were 32 μmol/L, with only slight variation within and between the lakes, and with depth. Mean available phosphorus was 0.07 μmol/L, resulting in relatively large N:P ratios (97:1) indicative of P limitation. Abundances of the phyla Actinobacteria, Bacteroidetes, Cyanobacteria, Thaumarchaeota, and Verrucomicrobia differed significantly among the Lakes. Candidatus Nitrosopumilus was present in greater abundance in Lake Superior compared to Lake Michigan, suggesting the importance of ammonia-oxidating archaea in water column N cycling in Lake Superior. The Shannon diversity index was negatively correlated with pH, temperature, and salinity, and positively correlated with DO, latitude, and N₂ saturation. Results of this study suggest that DO, pH, temperature, and salinity were major drivers shaping the community composition in the Great Lakes.

The impacts of the hydraulic retention effect and typhoon disturbance on the carbon flux in shallow subtropical mountain lakes

A typhoon is extreme weather that flushes terrestrial carbon (C) loads and temporally mixes the entire water columns of lakes in subtropical regions. A C flux varies based on the trophic level associated with the ecological cycle related to hydraulic retention time (residence time). Herein, we sought to clarify how the hydraulic retention time and the disturbance from a typhoon affect the C flux regimes in two subtropical mountain lakes in a humid region of Taiwan with different trophic levels-oligotrophic and mesotrophic. We investigated the meteorological data and vertical profiles of the water temperature, dissolved inorganic carbon (DIC), dissolved organic C (DOC), and chlorophyll a (Chl. a) during the pre-typhoon period (April-July), during the typhoon period (August-November), and the post-typhoon period (December-March) for five years (2009-2010 and 2015-2017). We applied a three-dimensional environmental model (Fantom) to investigate the hydraulic retention effect on the net ecosystem production (NEP) using the residence time in stratified lakes. The results demonstrate that typhoon-induced mixing associated with the hydraulic retention effect plays one of the critical roles in controlling the NEP and C flux in shallow subtropical lakes.

Impacts of Extreme Weather Events on Bacterial Community Composition of a Temperate Humic Lake

Extreme weather events are projected to increase in frequency and intensity as climate change continues. Heterotrophic bacteria play a critical role in lake ecosystems, yet little research has been done to determine how they are affected by such extremes. The purpose of this study was to use high-throughput sequencing to explore the bacterial community composition of a humic oligotrophic lake on the North Atlantic Irish coast and to assess the impacts on composition dynamics related to extreme weather events. Samples for sequencing were collected from Lough Feeagh on a fortnightly basis from April to November 2018. Filtration was used to separate free-living and particle-associated bacterial communities and amplicon sequencing was performed for the 16S rRNA V4 region. Two named storms, six high discharge events, and one drought period occurred during the sampling period. These events had variable, context-dependent effects on bacterial communities in Lough Feeagh. The particle-associated community was found to be more likely to respond to physical changes, such as mixing, while the free-living population responded to changes in nutrient and carbon concentrations. Generally, however, the high stability of the bacterial community observed in Lough Feeagh suggests that the bacterial community is relatively resilient to extreme weather events.

Microdiversity ensures the maintenance of functional microbial communities under changing environmental conditions

Microdiversity can lead to different ecotypes within the same species. These are assumed to provide stability in time and space to those species. However, the role of microdiversity in the stability of whole microbial communities remains underexplored. Understanding the drivers of microbial community stability is necessary to predict community response to future disturbances. Here, we analyzed 16S rRNA gene amplicons from eight different temperate bog lakes at the 97% OTU and amplicon sequence variant (ASV) levels and found ecotypes within the same OTU with different distribution patterns in space and time. We observed that these ecotypes are adapted to different values of environmental factors such as water temperature and oxygen concentration. Our results showed that the existence of several ASVs within a OTU favored its persistence across changing environmental conditions. We propose that microdiversity aids the stability of microbial communities in the face of fluctuations in environmental factors.

The sensitivity and stability of bacterioplankton community structure to wind-wave turbulence in a large, shallow, eutrophic lake

Lakes are strongly influenced by wind-driven wave turbulence. The direct physical effects of turbulence on bacterioplankton community structure however, have not yet been addressed and remains poorly understood. To examine the stability of bacterioplankton communities under turbulent conditions, we simulated conditions in the field to evaluate the responses of the bacterioplankton community to physical forcing in Lake Taihu, using high-throughput sequencing and flow cytometry. A total of 4,520,231 high quality sequence reads and 74,842 OTUs were obtained in all samples with α-proteobacteria, γ-proteobacteria and Actinobacteria being the most dominant taxa. The diversity and structure of bacterioplankton communities varied during the experiment, but were highly similar based on the same time of sampling, suggesting that bacterioplankton communities are insensitive to wind wave turbulence in the lake. This stability could be associated with the traits associated with bacteria. In particular, turbulence favored the growth of bacterioplankton, which enhanced biogeochemical cycling of nutrients in the lake. This study provides a better understanding of bacterioplankton communities in lake ecosystems exposed to natural mixing/disturbances.

Microbial community diversity patterns are related to physical and chemical differences among temperate lakes near Beaver Island, MI

Lakes are dynamic and complex ecosystems that can be influenced by physical, chemical, and biological processes. Additionally, individual lakes are often chemically and physically distinct, even within the same geographic region. Here we show that differences in physicochemical conditions among freshwater lakes located on (and around) the same island, as well as within the water column of each lake, are significantly related to aquatic microbial community diversity. Water samples were collected over time from the surface and bottom-water within four freshwater lakes located around Beaver Island, MI within the Laurentian Great Lakes region. Three of the sampled lakes experienced seasonal lake mixing events, impacting either O₂, pH, temperature, or a combination of the three. Microbial community alpha and beta diversity were assessed and individual microbial taxa were identified via high-throughput sequencing of the 16S rRNA gene. Results demonstrated that physical and chemical variability (temperature, dissolved oxygen, and pH) were significantly related to divergence in the beta diversity of surface and bottom-water microbial communities. Despite its correlation to microbial community structure in unconstrained analyses, constrained analyses demonstrated that dissolved organic carbon (DOC) concentration was not strongly related to microbial community structure among or within lakes. Additionally, several taxa were correlated (either positively or negatively) to environmental variables, which could be related to aerobic and anaerobic metabolisms. This study highlights the measurable relationships between environmental conditions and microbial communities within freshwater temperate lakes around the same island.

A unique assemblage of cosmopolitan freshwater bacteria and higher community diversity differentiate an urbanized estuary from oligotrophic Lake Michigan

Water quality is impacted significantly by urbanization. The delivery of increased nutrient loads to waterways is a primary characteristic of this land use change. Despite the recognized effects of nutrient loading on aquatic systems, the influence of urbanization on the bacterial community composition of these systems is not understood. We used massively-parallel sequencing of bacterial 16S rRNA genes to examine the bacterial assemblages in transect samples spanning the heavily urbanized estuary of Milwaukee, WI to the relatively un-impacted waters of Lake Michigan. With this approach, we found that genera and lineages common to freshwater lake epilimnia were common and abundant in both the high nutrient, urban-impacted waterways, and the low nutrient Lake Michigan. Although the two environments harbored many taxa in common, we identified a significant change in the community assemblage across the urban-influence gradient, and three distinct community features drove this change. First, we found the urban-influenced waterways harbored significantly greater bacterial richness and diversity than Lake Michigan (i.e., taxa augmentation). Second, we identified a shift in the relative abundance among common freshwater lineages, where acI, acTH1, Algoriphagus and LD12, had decreased representation and Limnohabitans, Polynucleobacter, and Rhodobacter had increased representation in the urban estuary. Third, by oligotyping 18 common freshwater genera/lineages, we found that oligotypes (highly resolved sequence clusters) within many of these genera/lineages had opposite preferences for the two environments. With these data, we suggest many of the defined cosmopolitan freshwater genera/lineages contain both oligotroph and more copiotroph species or populations, promoting the idea that within-genus lifestyle specialization, in addition to shifts in the dominance among core taxa and taxa augmentation, drive bacterial community change in urbanized waters.

Temporal patterns of rarity provide a more complete view of microbial diversity

Recently, conditionally rare taxa (CRTs)--those taxa that are typically in very low abundance but occasionally achieve prevalence--were shown to contribute to patterns of microbial diversity because their collective dynamics explained a large proportion of temporal variability in microbial community structure. Here the benefits and challenges of characterizing the presence and interpreting the role of CRTs are further explored, along with questions about CRT ecology. We also introduce a conceptual model for thinking about microbial taxa as dynamic components along the dimensions of occurrence and abundance. Accounting for CRTs in interpretations of microbial ecological dynamics is essential if we are to understand community stability and ecoevolutionary interactions.

The bias associated with amplicon sequencing does not affect the quantitative assessment of bacterial community dynamics

The performance of two sets of primers targeting variable regions of the 16S rRNA gene V1–V3 and V4 was compared in their ability to describe changes of bacterial diversity and temporal turnover in full-scale activated sludge. Duplicate sets of high-throughput amplicon sequencing data of the two 16S rRNA regions shared a collection of core taxa that were observed across a series of twelve monthly samples, although the relative abundance of each taxon was substantially different between regions. A case in point was the changes in the relative abundance of filamentous bacteria Thiothrix, which caused a large effect on diversity indices, but only in the V1–V3 data set. Yet the relative abundance of Thiothrix in the amplicon sequencing data from both regions correlated with the estimation of its abundance determined using fluorescence in situ hybridization. In nonmetric multidimensional analysis samples were distributed along the first ordination axis according to the sequenced region rather than according to sample identities. The dynamics of microbial communities indicated that V1–V3 and the V4 regions of the 16S rRNA gene yielded comparable patterns of: 1) the changes occurring within the communities along fixed time intervals, 2) the slow turnover of activated sludge communities and 3) the rate of species replacement calculated from the taxa–time relationships. The temperature was the only operational variable that showed significant correlation with the composition of bacterial communities over time for the sets of data obtained with both pairs of primers. In conclusion, we show that despite the bias introduced by amplicon sequencing, the variable regions V1–V3 and V4 can be confidently used for the quantitative assessment of bacterial community dynamics, and provide a proper qualitative account of general taxa in the community, especially when the data are obtained over a convenient time window rather than at a single time point.

Assessing the effects of severe rainstorm-induced mixing on a subtropical, subalpine lake

Severe rainstorms cause vertical mixing that modifies the internal dynamics (e.g., internal seiche, thermal structure, and velocity filed) in warm polymictic lakes. Yuan Yang Lake (YYL), a subtropical, subalpine, and seasonally stratified small lake in the north-central region of Taiwan, is normally affected by typhoons accompanied with strong wind and heavy rainfall during the summer and fall. In this study, we used the field data, statistical analysis, spectral analysis, and numerical modeling to investigate severe rainstorm-induced mixing in the lake. Statistical determination of the key meteorological and environmental conditions underlying the observed vertical mixing suggests that the vertical mixing, caused by heat loss during severe rainstorms, was likely larger than wind-induced mixing and that high inflow discharge strongly increased heat loss through advection heat. Spectral analysis revealed that internal seiches at the basin scale occurred under non-rainstorm meteorological conditions and that the internal seiches under the rainstorm were modified on the increase of the internal seiche frequencies. Based upon observed frequencies of the internal seiches, a two-dimensional model was simulated and then appropriate velocity patterns of the internal seiches were determined under non-rainstorm conditions. Moreover, the model implemented with inflow boundary condition was conducted for rainstorm events. The model results showed that the severe rainstorms promoted thermal destratification and changed vertical circulation of the basin-scale, internal seiche motion into riverine flow.

Seasonal variations and resilience of bacterial communities in a sewage polluted urban river

The Zenne River in Brussels (Belgium) and effluents of the two wastewater treatment plants (WWTPs) of Brussels were chosen to assess the impact of disturbance on bacterial community composition (BCC) of an urban river. Organic matters, nutrients load and oxygen concentration fluctuated highly along the river and over time because of WWTPs discharge. Tag pyrosequencing of bacterial 16S rRNA genes revealed the significant effect of seasonality on the richness, the bacterial diversity (Shannon index) and BCC. The major grouping: -winter/fall samples versus spring/summer samples- could be associated with fluctuations of in situ bacterial activities (dissolved and particulate organic carbon biodegradation associated with oxygen consumption and N transformation). BCC of the samples collected upstream from the WWTPs discharge were significantly different from BCC of downstream samples and WWTPs effluents, while no significant difference was found between BCC of WWTPs effluents and the downstream samples as revealed by ANOSIM. Analysis per season showed that allochthonous bacteria brought by WWTPs effluents triggered the changes in community composition, eventually followed by rapid post-disturbance return to the original composition as observed in April (resilience), whereas community composition remained altered after the perturbation by WWTPs effluents in the other seasons.

Fine-scale spatial patterns in bacterial community composition and function within freshwater ponds

The extent to which non-host-associated bacterial communities exhibit small-scale biogeographic patterns in their distribution remains unclear. Our investigation of biogeography in bacterial community composition and function compared samples collected across a smaller spatial scale than most previous studies conducted in freshwater. Using a grid-based sampling design, we abstracted 100+ samples located between 3.5 and 60 m apart within each of three alpine ponds. For every sample, variability in bacterial community composition was monitored using a DNA-fingerprinting methodology (automated ribosomal intergenic spacer analysis) whereas differences in bacterial community function (that is, carbon substrate utilisation patterns) were recorded from Biolog Ecoplates. The exact spatial position and dominant physicochemical conditions (for example, pH and temperature) were simultaneously recorded for each sample location. We assessed spatial differences in bacterial community composition and function within each pond and found that, on average, community composition or function differed significantly when comparing samples located >20 m apart within any pond. Variance partitioning revealed that purely spatial variation accounted for more of the observed variability in both bacterial community composition and function (range: 24-38% and 17-39%) than the combination of purely environmental variation and spatially structured environmental variation (range: 17-32% and 15-20%). Clear spatial patterns in bacterial community composition, but not function were observed within ponds. We therefore suggest that some of the observed variation in bacterial community composition is functionally 'redundant'. We confirm that distinct bacterial communities are present across unexpectedly small spatial scales suggesting that populations separated by distances of >20 m may be dispersal limited, even within the highly continuous environment of lentic water.

Resource availability and spatial heterogeneity control bacterial community response to nutrient enrichment in lakes

The diversity and composition of ecological communities often co-vary with ecosystem productivity. However, the relative importance of productivity, or resource abundance, versus the spatial distribution of resources in shaping those ecological patterns is not well understood, particularly for the bacterial communities that underlie most important ecosystem functions. Increasing ecosystem productivity in lakes has been shown to influence the composition and ecology of bacterial communities, but existing work has only evaluated the effect of increasing resource supply and not heterogeneity in how those resources are distributed. We quantified how bacterial communities varied with the trophic status of lakes and whether community responses differed in surface and deep habitats in response to heterogeneity in nutrient resources. Using ARISA fingerprinting, we found that bacterial communities were more abundant, richer, and more distinct among habitats as lake trophic state and vertical heterogeneity in nutrients increased, and that spatial resource variation produced habitat specific responses of bacteria in response to increased productivity. Furthermore, changes in communities in high nutrient lakes were not produced by turnover in community composition but from additional taxa augmenting core bacterial communities found in lower productivity lakes. These data suggests that bacterial community responses to nutrient enrichment in lakes vary spatially and are likely influenced disproportionately by rare taxa.

Lineage-specific responses of microbial communities to environmental change

A great challenge facing microbial ecology is how to define ecologically relevant taxonomic units. To address this challenge, we investigated how changing the definition of operational taxonomic units (OTUs) influences the perception of ecological patterns in microbial communities as they respond to a dramatic environmental change. We used pyrosequenced tags of the bacterial V2 16S rRNA region, as well as clone libraries constructed from the cytochrome oxidase C gene ccoN, to provide additional taxonomic resolution for the common freshwater genus Polynucleobacter. At the most highly resolved taxonomic scale, we show that distinct genotypes associated with the abundant Polynucleobacter lineages exhibit divergent spatial patterns and dramatic changes over time, while the also abundant Actinobacteria OTUs are highly coherent. This clearly demonstrates that different bacterial lineages demand different taxonomic definitions to capture ecological patterns. Based on the temporal distribution of highly resolved taxa in the hypolimnion, we demonstrate that change in the population structure of a single genotype can provide additional insight into the mechanisms of community-level responses. These results highlight the importance and feasibility of examining ecological change in microbial communities across taxonomic scales while also providing valuable insight into the ecological characteristics of ecologically coherent groups in this system.

Lake microbial communities are resilient after a whole-ecosystem disturbance

Disturbances act as powerful structuring forces on ecosystems. To ask whether environmental microbial communities have capacity to recover after a large disturbance event, we conducted a whole-ecosystem manipulation, during which we imposed an intense disturbance on freshwater microbial communities by artificially mixing a temperate lake during peak summer thermal stratification. We employed environmental sensors and water chemistry analyses to evaluate the physical and chemical responses of the lake, and bar-coded 16S ribosomal RNA gene pyrosequencing and automated ribosomal intergenic spacer analysis (ARISA) to assess the bacterial community responses. The artificial mixing increased mean lake temperature from 14 to 20 °C for seven weeks after mixing ended, and exposed the microorganisms to very different environmental conditions, including increased hypolimnion oxygen and increased epilimnion carbon dioxide concentrations. Though overall ecosystem conditions remained altered (with hypolimnion temperatures elevated from 6 to 20 °C), bacterial communities returned to their pre-manipulation state as some environmental conditions, such as oxygen concentration, recovered. Recovery to pre-disturbance community composition and diversity was observed within 7 (epilimnion) and 11 (hypolimnion) days after mixing. Our results suggest that some microbial communities have capacity to recover after a major disturbance.

Beyond the Venn diagram: the hunt for a core microbiome

Discovering a core microbiome is important for understanding the stable, consistent components across complex microbial assemblages. A core is typically defined as the suite of members shared among microbial consortia from similar habitats, and is represented by the overlapping areas of circles in Venn diagrams, in which each circle contains the membership of the sample or habitats being compared. Ecological insight into core microbiomes can be enriched by 'omics approaches that assess gene expression, thereby extending the concept of the core beyond taxonomically defined membership to community function and behaviour. Parameters defined by traditional ecology theory, such as composition, phylogeny, persistence and connectivity, will also create a more complex portrait of the core microbiome and advance understanding of the role of key microorganisms and functions within and across ecosystems.

Resistance, resilience and recovery: aquatic bacterial dynamics after water column disturbance

For lake microbes, water column mixing acts as a disturbance because it homogenizes thermal and chemical gradients known to define the distributions of microbial taxa. Our first objective was to isolate hypothesized drivers of lake bacterial response to water column mixing. To accomplish this, we designed an enclosure experiment with three treatments to independently test key biogeochemical changes induced by mixing: oxygen addition to the hypolimnion, nutrient addition to the epilimnion, and full water column mixing. We used molecular fingerprinting to observe bacterial community dynamics in the treatment and control enclosures, and in ambient lake water. We found that oxygen and nutrient amendments simulated the physical-chemical water column environment following mixing and resulted in similar bacterial communities to the mixing treatment, affirming that these were important drivers of community change. These results demonstrate that specific environmental changes can replicate broad disturbance effects on microbial communities. Our second objective was to characterize bacterial community stability by quantifying community resistance, recovery and resilience to an episodic disturbance. The communities in the nutrient and oxygen amendments changed quickly (had low resistance), but generally matched the control composition by the 10th day after treatment, exhibiting resilience. These results imply that aquatic bacterial assemblages are generally stable in the face of disturbance.

Our microbial selves: what ecology can teach us

Advances in DNA sequencing have allowed us to characterize microbial communities--including those associated with the human body--at a broader range of spatial and temporal scales than ever before. We can now answer fundamental questions that were previously inaccessible and use well-tested ecological theories to gain insight into changes in the microbiome that are associated with normal development and human disease. Perhaps unsurprisingly, the ecosystems associated with our body follow trends identified in communities at other sites and scales, and thus studies of the microbiome benefit from ecological insight. Here, we assess human microbiome research in the context of ecological principles and models, focusing on diversity, biological drivers of community structure, spatial patterning and temporal dynamics, and suggest key directions for future research that will bring us closer to the goal of building predictive models for personalized medicine.