Metabarcoding reveals unique microbial mat communities and evidence of biogeographic influence in low-oxygen, high-sulfur sinkholes and springs.

High-sulfur, low-oxygen environments formed by underwater sinkholes and springs create unique habitats populated by microbial mat communities. To explore the diversity and biogeography of these mats, samples were collected from three sites in Alpena, Michigan, one site in Monroe, Michigan, and one site in Palm Coast, Florida. Our study investigated previously undescribed eukaryotic diversity in these habitats and further explored their bacterial communities. Mat samples and water parameters were collected from sulfur spring sites during the spring, summer, and fall of 2022. Cyanobacteria and diatoms were cultured from mat subsamples to create a culture-based DNA reference library. Remaining mat samples were used for metabarcoding of the 16S and rbcL regions to explore bacterial and diatom diversity, respectively. Analyses of water chemistry, alpha diversity, and beta diversity articulated a range of high-sulfur, low-oxygen habitats, each with distinct microbial communities. Conductivity, pH, dissolved oxygen, temperature, sulfate, and chloride had significant influences on community composition but did not describe the differences between communities well. Chloride concentration had the strongest correlation with microbial community structure. Mantel tests revealed that biogeography contributed to differences between communities as well. Our results provide novel information on microbial mat composition and present evidence that both local conditions and biogeography influence these unique communities.

Unbiased analyses of ITS folding motifs in a taxonomically confusing lineage: Anagnostidinema visiae sp. nov. (cyanobacteria).

Cyanobacteria are diverse prokaryotic, photosynthetic organisms present in nearly every known ecosystem. Recent investigations around the world have recovered vast amounts of novel biodiversity in seldom sampled habitats. One phylogenetically significant character, the secondary folding structures of the 16S-23S ITS rDNA region, has allowed an unprecedented capacity to erect new species. However, two questions arise: Is this feature as informative as is proposed, and how do we best employ these features? Submerged sinkholes with oxygen-poor, sulfur-rich ground water in Lake Huron (USA) contain microbial mats dominated by both oxygenic and anoxygenic cyanobacteria. We sought to document some of this unique cyanobacterial diversity. Using culture-based investigations, we recovered 45 strains, of which 23 were analyzed employing 16S-23S rDNA sequences, ITS folding patterns, ecology, and morphology. With scant morphological discontinuities and nebulous 16S rDNA gene sequence divergence, ITS folding patterns were effective at articulating cryptic biodiversity. However, we would have missed these features had we not folded all the available motifs from the strains, including those with highly similar 16S rDNA gene sequences. If we had relied solely on morphological or 16S rDNA gene data, then we might well have missed the diversity of Anagnostidinema. Thus, in order to avoid conformation basis, which is potentially common when employing ITS structures, we advocate clustering strains based on ITS rDNA region patterns independently and comparing them back to 16S rDNA gene phylogenies. Using a total evidence approach, we erected a new taxon according to the International Code of Nomenclature for Algae, Fungi, and Plants: Anagnostidinema visiae.

Omics-Inferred Partitioning and Expression of Diverse Biogeochemical Functions in a Low-O2 Cyanobacterial Mat Community.

Cyanobacterial mats profoundly influenced Earth's biological and geochemical evolution and still play important ecological roles in the modern world. However, the biogeochemical functioning of cyanobacterial mats under persistent low-O2 conditions, which dominated their evolutionary history, is not well understood. To investigate how different metabolic and biogeochemical functions are partitioned among community members, we conducted metagenomics and metatranscriptomics on cyanobacterial mats in the low-O2, sulfidic Middle Island sinkhole (MIS) in Lake Huron. Metagenomic assembly and binning yielded 144 draft metagenome assembled genomes, including 61 of medium quality or better, and the dominant cyanobacteria and numerous Proteobacteria involved in sulfur cycling. Strains of a Phormidium autumnale-like cyanobacterium dominated the metagenome and metatranscriptome. Transcripts for the photosynthetic reaction core genes psaA and psbA were abundant in both day and night. Multiple types of psbA genes were expressed from each cyanobacterium, and the dominant psbA transcripts were from an atypical microaerobic type of D1 protein from Phormidium. Further, cyanobacterial transcripts for photosystem I genes were more abundant than those for photosystem II, and two types of Phormidium sulfide quinone reductase were recovered, consistent with anoxygenic photosynthesis via photosystem I in the presence of sulfide. Transcripts indicate active sulfur oxidation and reduction within the cyanobacterial mat, predominately by Gammaproteobacteria and Deltaproteobacteria, respectively. Overall, these genomic and transcriptomic results link specific microbial groups to metabolic processes that underpin primary production and biogeochemical cycling in a low-O2 cyanobacterial mat and suggest mechanisms for tightly coupled cycling of oxygen and sulfur compounds in the mat ecosystem. IMPORTANCE Cyanobacterial mats are dense communities of microorganisms that contain photosynthetic cyanobacteria along with a host of other bacterial species that play important yet still poorly understood roles in this ecosystem. Although such cyanobacterial mats were critical agents of Earth's biological and chemical evolution through geological time, little is known about how they function under the low-oxygen conditions that characterized most of their natural history. Here, we performed sequencing of the DNA and RNA of modern cyanobacterial mat communities under low-oxygen and sulfur-rich conditions from the Middle Island sinkhole in Lake Huron. The results reveal the organisms and metabolic pathways that are responsible for both oxygen-producing and non-oxygen-producing photosynthesis as well as interconversions of sulfur that likely shape how much O2 is produced in such ecosystems. These findings indicate tight metabolic reactions between community members that help to explain the limited the amount of O2 produced in cyanobacterial mat ecosystems.

Microhabitats are associated with diversity-productivity relationships in freshwater bacterial communities.

Eukaryotic communities commonly display a positive relationship between biodiversity and ecosystem function (BEF) but the results have been mixed when assessed in bacterial communities. Habitat heterogeneity, a factor in eukaryotic BEFs, may explain these variable observations but it has not been thoroughly evaluated in bacterial communities. Here, we examined the impact of habitat on the relationship between diversity assessed based on the (phylogenetic) Hill diversity metrics and heterotrophic productivity. We sampled co-occurring free-living (more homogenous) and particle-associated (more heterogeneous) bacterial habitats in a freshwater, estuarine lake over three seasons: spring, summer and fall. There was a strong, positive, linear relationship between particle-associated bacterial richness and heterotrophic productivity that strengthened when considering dominant taxa. There were no observable BEF trends in free-living bacterial communities for any diversity metric. Biodiversity, richness and Inverse Simpson's index, were the best predictors of particle-associated production whereas pH was the best predictor of free-living production. Our findings show that heterotrophic productivity is positively correlated with the effective number of taxa and that BEF relationships are associated with microhabitats. These results add to the understanding of the highly distinct contributions to diversity and functioning contributed by bacteria in free-living and particle-associated habitats.

Randomized Lasso Links Microbial Taxa with Aquatic Functional Groups Inferred from Flow Cytometry.

High-nucleic-acid (HNA) and low-nucleic-acid (LNA) bacteria are two operational groups identified by flow cytometry (FCM) in aquatic systems. A number of reports have shown that HNA cell density correlates strongly with heterotrophic production, while LNA cell density does not. However, which taxa are specifically associated with these groups, and by extension, productivity has remained elusive. Here, we addressed this knowledge gap by using a machine learning-based variable selection approach that integrated FCM and 16S rRNA gene sequencing data collected from 14 freshwater lakes spanning a broad range in physicochemical conditions. There was a strong association between bacterial heterotrophic production and HNA absolute cell abundances (R 2 = 0.65), but not with the more abundant LNA cells. This solidifies findings, mainly from marine systems, that HNA and LNA bacteria could be considered separate functional groups, the former contributing a disproportionately large share of carbon cycling. Taxa selected by the models could predict HNA and LNA absolute cell abundances at all taxonomic levels. Selected operational taxonomic units (OTUs) ranged from low to high relative abundance and were mostly lake system specific (89.5% to 99.2%). A subset of selected OTUs was associated with both LNA and HNA groups (12.5% to 33.3%), suggesting either phenotypic plasticity or within-OTU genetic and physiological heterogeneity. These findings may lead to the identification of system-specific putative ecological indicators for heterotrophic productivity. Generally, our approach allows for the association of OTUs with specific functional groups in diverse ecosystems in order to improve our understanding of (microbial) biodiversity-ecosystem functioning relationships.IMPORTANCE A major goal in microbial ecology is to understand how microbial community structure influences ecosystem functioning. Various methods to directly associate bacterial taxa to functional groups in the environment are being developed. In this study, we applied machine learning methods to relate taxonomic data obtained from marker gene surveys to functional groups identified by flow cytometry. This allowed us to identify the taxa that are associated with heterotrophic productivity in freshwater lakes and indicated that the key contributors were highly system specific, regularly rare members of the community, and that some could possibly switch between being low and high contributors. Our approach provides a promising framework to identify taxa that contribute to ecosystem functioning and can be further developed to explore microbial contributions beyond heterotrophic production.

Correction: Verrucomicrobia are prevalent in north-temperate freshwater lakes and display class-level preferences between lake habitats.

[This corrects the article DOI: 10.1371/journal.pone.0195112.].

Verrucomicrobia are prevalent in north-temperate freshwater lakes and display class-level preferences between lake habitats.

The bacterial phylum Verrucomicrobia was formally described two decades ago and originally believed to be a minor member of many ecosystems; however, it is now recognized as ubiquitous and abundant in both soil and aquatic systems. Nevertheless, knowledge of the drivers of its relative abundance and within-phylum habitat preferences remains sparse, especially in lake systems. Here, we documented the distribution of Verrucomicrobia in 12 inland lakes in Southeastern Michigan, a Laurentian Great Lake (Lake Michigan), and a freshwater estuary, which span a gradient in lake sizes, depths, residence times, and trophic states. A wide range of physical and geochemical parameters was covered by sampling seasonally from the surface and bottom of each lake, and by separating samples into particle-associated and free-living fractions. On average, Verrucomicrobia was the 4th most abundant phylum (range 1.7-41.7%). Fraction, season, station, and depth explained up to 70% of the variance in Verrucomicrobia community composition and preference for these habitats was phylogenetically conserved at the class-level. When relative abundance was linearly modeled against environmental data, Verrucomicrobia and non-Verrucomicrobia bacterial community composition correlated to similar quantitative environmental parameters, although there were lake system-dependent differences and > 55% of the variance remained unexplained. A majority of the phylum exhibited preference for the particle-associated fraction and two classes (Opitutae and Verrucomicrobiae) were identified to be more abundant during the spring season. This study highlights the high relative abundance of Verrucomicrobia in north temperate lake systems and expands insights into drivers of within-phylum habitat preferences of the Verrucomicrobia.

Predominant Non-additive Effects of Multiple Stressors on Autotroph C:N:P Ratios Propagate in Freshwater and Marine Food Webs.

A continuing challenge for scientists is to understand how multiple interactive stressor factors affect biological interactions, and subsequently, ecosystems-in ways not easily predicted by single factor studies. In this review, we have compiled and analyzed available research on how multiple stressor pairs composed of temperature (T), light (L), ultraviolet radiation (UVR), nutrients (Nut), carbon dioxide (CO2), dissolved organic carbon (DOC), and salinity (S) impact the stoichiometry of autotrophs which in turn shapes the nature of their ecological interactions within lower trophic levels in streams, lakes and oceans. Our analysis from 66 studies with 320 observations of 11 stressor pairs, demonstrated that non-additive responses predominate across aquatic ecosystems and their net interactive effect depends on the stressor pair at play. Across systems, there was a prevalence of antagonism in freshwater (60-67% vs. 47% in marine systems) compared to marine systems where synergism was more common (49% vs. 33-40% in freshwaters). While the lack of data impeded comparisons among all of the paired stressors, we found pronounced system differences for the L × Nut interactions. For this interaction, our data for C:P and N:P is consistent with the initial hypothesis that the interaction was primarily synergistic in the oceans, but not for C:N. Our study found a wide range of variability in the net effects of the interactions in freshwater systems, with some observations supporting antagonism, and others synergism. Our results suggest that the nature of the stressor pairs interactions on C:N:P ratios regulates the "continuum" commensalistic-competitive-predatory relationship between algae and bacteria and the food chain efficiency at the algae-herbivore interface. Overall, the scarce number of studies with even more fewer replications in each study that are available for freshwater systems have prevented a more detailed, insightful analysis. Our findings highlighting the preponderance of antagonistic and synergistic effects of stressor interactions in aquatic ecosystems-effects that play key roles in the functioning of feedback loops in the biosphere-also stress the need for further studies evaluating the interactive effects of multiple stressors in a rapidly changing world facing a confluence of tipping points.

Novel Large Sulfur Bacteria in the Metagenomes of Groundwater-Fed Chemosynthetic Microbial Mats in the Lake Huron Basin.

Little is known about large sulfur bacteria (LSB) that inhabit sulfidic groundwater seeps in large lakes. To examine how geochemically relevant microbial metabolisms are partitioned among community members, we conducted metagenomic analysis of a chemosynthetic microbial mat in the Isolated Sinkhole, which is in a deep, aphotic environment of Lake Huron. For comparison, we also analyzed a white mat in an artesian fountain that is fed by groundwater similar to Isolated Sinkhole, but that sits in shallow water and is exposed to sunlight. De novo assembly and binning of metagenomic data from these two communities yielded near complete genomes and revealed representatives of two families of LSB. The Isolated Sinkhole community was dominated by novel members of the Beggiatoaceae that are phylogenetically intermediate between known freshwater and marine groups. Several of these Beggiatoaceae had 16S rRNA genes that contained introns previously observed only in marine taxa. The Alpena fountain was dominated by populations closely related to Thiothrix lacustris and an SM1 euryarchaeon known to live symbiotically with Thiothrix spp. The SM1 genomic bin contained evidence of H2-based lithoautotrophy. Genomic bins of both the Thiothrix and Beggiatoaceae contained genes for sulfur oxidation via the rDsr pathway, H2 oxidation via Ni-Fe hydrogenases, and the use of O2 and nitrate as electron acceptors. Mats at both sites also contained Deltaproteobacteria with genes for dissimilatory sulfate reduction (sat, apr, and dsr) and hydrogen oxidation (Ni-Fe hydrogenases). Overall, the microbial mats at the two sites held low-diversity microbial communities, displayed evidence of coupled sulfur cycling, and did not differ largely in their metabolic potentials, despite the environmental differences. These results show that groundwater-fed communities in an artesian fountain and in submerged sinkholes of Lake Huron are a rich source of novel LSB, associated heterotrophic and sulfate-reducing bacteria, and archaea.

Ecological and genetic interactions between cyanobacteria and viruses in a low-oxygen mat community inferred through metagenomics and metatranscriptomics.

Metagenomic and metatranscriptomic sequencing was conducted on cyanobacterial mats of the Middle Island Sinkhole (MIS), Lake Huron. Metagenomic data from 14 samples collected over 5 years were used to reconstruct genomes of two genotypes of a novel virus, designated PhV1 type A and PhV1 type B. Both viral genotypes encode and express nblA, a gene involved in degrading phycobilisomes, which are complexes of pigmented proteins that harvest light for photosynthesis. Phylogenetic analysis indicated that the viral-encoded nblA is derived from the host cyanobacterium, Phormidium MIS-PhA. The cyanobacterial host also has two complete CRISPR (clustered regularly interspaced short palindromic repeats) systems that serve as defence mechanisms for bacteria and archaea against viruses and plasmids. One 45 bp CRISPR spacer from Phormidium had 100% nucleotide identity to PhV1 type B, but this region was absent from PhV1 type A. Transcripts from PhV1 and the Phormidium CRISPR loci were detected in all six metatranscriptomic data sets (three during the day and three at night), indicating that both are transcriptionally active in the environment. These results reveal ecological and genetic interactions between viruses and cyanobacteria at MIS, highlighting the value of parallel analysis of viruses and hosts in understanding ecological interactions in natural communities.

Seeking sunlight: rapid phototactic motility of filamentous mat-forming cyanobacteria optimize photosynthesis and enhance carbon burial in Lake Huron's submerged sinkholes.

We studied the motility of filamentous mat-forming cyanobacteria consisting primarily of Oscillatoria-like cells growing under low-light, low-oxygen, and high-sulfur conditions in Lake Huron's submerged sinkholes using in situ observations, in vitro measurements and time-lapse microscopy. Gliding movement of the cyanobacterial trichomes (100-10,000 µm long filaments, composed of cells ∼10 µm wide and ∼3 µm tall) revealed individual as well as group-coordinated motility. When placed in a petri dish and dispersed in ground water from the sinkhole, filaments re-aggregated into defined colonies within minutes, then dispersed again. Speed of individual filaments increased with temperature from ∼50 µm min(-1) or ∼15 body lengths min(-1) at 10°C to ∼215 µm min(-1) or ∼70 body lengths min(-1) at 35°C - rates that are rapid relative to non-flagellated/ciliated microbes. Filaments exhibited precise and coordinated positive phototaxis toward pinpoints of light and congregated under the light of foil cutouts. Such light-responsive clusters showed an increase in photosynthetic yield - suggesting phototactic motility aids in light acquisition as well as photosynthesis. Once light source was removed, filaments slowly spread out evenly and re-aggregated, demonstrating coordinated movement through inter-filament communication regardless of light. Pebbles and pieces of broken shells placed upon intact mat were quickly covered by vertically motile filaments within hours and became fully buried in the anoxic sediments over 3-4 diurnal cycles - likely facilitating the preservation of falling debris. Coordinated horizontal and vertical filament motility optimize mat cohesion and dynamics, photosynthetic efficiency and sedimentary carbon burial in modern-day sinkhole habitats that resemble the shallow seas in Earth's early history. Analogous cyanobacterial motility may have played a key role in the oxygenation of the planet by optimizing photosynthesis while favoring carbon burial.

Systematically variable planktonic carbon metabolism along a land-to-lake gradient in a Great Lakes coastal zone.

During the summers of 2002-2013, we measured rates of carbon metabolism in surface waters of six sites across a land-to-lake gradient from the upstream end of drowned river-mouth Muskegon Lake (ML) (freshwater estuary) to 19 km offshore in Lake Michigan (LM) (a Great Lake). Despite considerable inter-year variability, the average rates of gross production (GP), respiration (R) and net production (NP) across ML (604 ± 58, 222 ± 22 and 381 ± 52 µg C L-1 day-1, respectively) decreased steeply in the furthest offshore LM site (22 ± 3, 55 ± 17 and -33 ± 15 µg C L-1day-1, respectively). Along this land-to-lake gradient, GP decreased by 96 ± 1%, whereas R only decreased by 75 ± 9%, variably influencing the carbon balance along this coastal zone. All ML sites were consistently net autotrophic (mean GP:R = 2.7), while the furthest offshore LM site was net heterotrophic (mean GP:R = 0.4). Our study suggests that pelagic waters of this Great Lakes coastal estuary are net carbon sinks that transition into net carbon sources offshore. Reactive and dynamic estuarine coastal zones everywhere may contribute similarly to regional and global carbon cycles.

Benthic bacterial diversity in submerged sinkhole ecosystems.

Physicochemical characterization, automated ribosomal intergenic spacer analysis (ARISA) community profiling, and 16S rRNA gene sequencing approaches were used to study bacterial communities inhabiting submerged Lake Huron sinkholes inundated with hypoxic, sulfate-rich groundwater. Photosynthetic cyanobacterial mats on the sediment surface were dominated by Phormidium autumnale, while deeper, organically rich sediments contained diverse and active bacterial communities.

Temperature affects stoichiometry and biochemical composition of Escherichia coli.

Temperature is a master variable controlling biochemical processes in organisms, and its effects are manifested on many organizational levels in organisms and ecosystems. We examined the effects of temperature on the biochemical composition and stoichiometry of a model heterotrophic bacterium, Escherichia coli K-12, held at constant growth rate in chemostats. Increasing temperature led to increased cellular organic carbon (C) and organic nitrogen (N) with decreased phosphorus (P) content, leading to increased C/P and N/P biomass ratios. P content was related to cellular RNA, which is P-rich (9-10% by weight) and nonnucleic acid P (presumably composed of mostly phospholipids, intracellular phosphate, and polyphosphate). These results indicate that E. coli allocates an increased proportion of its P cell quota toward assembly (ribosomes) at low temperatures and an increasing proportion toward resource acquisition machinery (membranes) at higher temperatures. If these results are relevant to the behavior of prokaryotic heterotrophs in natural settings (the gut, soils, lakes, oceans, etc.), it suggests greater nutrient regeneration and less microbial nutrient retention as temperatures increase.