Primer Design for an Accurate View of Picocyanobacterial Community Structure by Using High-Throughput Sequencing

Ecophysiological analysis reveals distinct environmental preferences in closely related Baltic Sea picocyanobacteria

Cluster 5 picocyanobacteria significantly contribute to primary productivity in aquatic ecosystems. Estuarine populations are highly diverse and consist of many co-occurring strains, but their physiology remains largely understudied. In this study, we characterized 17 novel estuarine picocyanobacterial strains. Phylogenetic analysis of the 16S rRNA and pigment genes (cpcB and cpeBA) uncovered multiple estuarine and freshwater-related clusters and pigment types. Assays with five representative strains (three phycocyanin rich and two phycoerythrin rich) under temperature (10-30°C), light (10-190 μmol photons m-2  s-1 ), and salinity (2-14 PSU) gradients revealed distinct growth optima and tolerance, indicating that genetic variability was accompanied by physiological diversity. Adaptability to environmental conditions was associated with differential pigment content and photosynthetic performance. Amplicon sequence variants at a coastal and an offshore station linked population dynamics with phylogenetic clusters, supporting that strains isolated in this study represent key ecotypes within the Baltic Sea picocyanobacterial community. The functional diversity found within strains with the same pigment type suggests that understanding estuarine picocyanobacterial ecology requires analysis beyond the phycocyanin and phycoerythrin divide. This new knowledge of the environmental preferences in estuarine picocyanobacteria is important for understanding and evaluating productivity in current and future ecosystems.

Seasonal dynamics in picocyanobacterial abundance and clade composition at coastal and offshore stations in the Baltic Sea

Picocyanobacteria (< 2 µm in diameter) are significant contributors to total phytoplankton biomass. Due to the high diversity within this group, their seasonal dynamics and relationship with environmental parameters, especially in brackish waters, are largely unknown. In this study, the abundance and community composition of phycoerythrin rich picocyanobacteria (PE-SYN) and phycocyanin rich picocyanobacteria (PC-SYN) were monitored at a coastal (K-station) and at an offshore station (LMO; ~ 10 km from land) in the Baltic Sea over three years (2018–2020). Cell abundances of picocyanobacteria correlated positively to temperature and negatively to nitrate (NO₃) concentration. While PE-SYN abundance correlated to the presence of nitrogen fixers, PC-SYN abundance was linked to stratification/shallow waters. The picocyanobacterial targeted amplicon sequencing revealed an unprecedented diversity of 2169 picocyanobacterial amplicons sequence variants (ASVs). A unique assemblage of distinct picocyanobacterial clades across seasons was identified. Clade A/B dominated the picocyanobacterial community, except during summer when low NO_3, high phosphate (PO₄) concentrations and warm temperatures promoted S5.2 dominance. This study, providing multiyear data, links picocyanobacterial populations to environmental parameters. The difference in the response of the two functional groups and clades underscore the need for further high-resolution studies to understand their role in the ecosystem.

Seasonal and Spatial Variations in Synechococcus Abundance and Diversity Throughout the Gullmar Fjord, Swedish Skagerrak

The picophytoplankton Synechococcus is a globally abundant autotroph that contributes significantly to primary production in the oceans and coastal areas. These cyanobacteria constitute a diverse genus of organisms that have developed independent niche spaces throughout aquatic environments. Here, we use the 16S V3–V4 rRNA gene region and flow cytometry to explore the diversity of Synechococcus within the picophytoplankton community in the Gullmar Fjord, on the west coast of Sweden. We conducted a station-based 1-year time series and two transect studies of the fjord. Our analysis revealed that within the large number of Synechococcus amplicon sequence variants (ASVs; 239 in total), prevalent ASVs phylogenetically clustered with clade representatives in both marine subcluster 5.1 and 5.2. The near-surface composition of ASVs shifted from spring to summer, when a 5.1 subcluster dominated community developed along with elevated Synechococcus abundances up to 9.3 × 104 cells ml–1. This seasonal dominance by subcluster 5.1 was observed over the length of the fjord (25 km), where shifts in community composition were associated with increasing depth. Unexpectedly, the community shift was not associated with changes in salinity. Synechococcus abundance dynamics also differed from that of the photosynthetic picoeukaryote community. These results highlight how seasonal variations in environmental conditions influence the dynamics of Synechococcus clades in a high latitude threshold fjord.

Novel primers for 16S rRNA gene-based archaeal and bacterial community analysis in oceanic trench sediments

High-throughput sequencing of the 16S ribosomal RNA (16S rRNA) gene has been successfully applied to explore the microbial structure and dynamics in various environments. The distinctive microbial communities in oceanic trench sediments are expected because of the extremely high pressure and V-shape topology that caused the isolation from the other marine sediments. However, they have only been primarily targeted using 'universal' primers that provide variable performances for different environments. It is necessary to design specific primers to improve the detection resolution of unique microbial groups in oceanic trenches. Here, we designed one pair of bacterial and two pairs of archaeal specific primers based on 16S rRNA gene full-length sequences that truly come from trench sediment and tested their performances in 30 oceanic trench sediment samples. An in silico analysis showed that the V3-V4 hypervariable region was the most informative and representative for oceanic trench microbial groups. Compared with the 'universal' primers, 46 bacterial families were only detected by newly designed primer B344F/B749R, and eight archaeal families were only detected by the newly designed primer A306F/A713R which covered the one or two orders of magnitude more ASVs (amplicon sequence variants) (1,470,216) in the tested total 30 samples. Moreover, A306F/A713R had the largest number of observed ASVs suggesting its better performance in discovering more archaeal species which were easily ignored in universal primer-based experiments for oceanic trench sediments. The novel primers designed in this research could be a better option to access the unique microbial communities in extreme oceanic trench sediments.Key points• Defining V3-V4 as the most adequate hypervariable region for archaea and bacteria from oceanic trench sediments.• Three sets of bacterial and archaeal primers appear validity and advantage in revealing the real trench microbial communities.• The novel primers provide a better option to specifically detect the unique microbial communities in extreme oceanic trench sediments.

Comparing microscopy and DNA metabarcoding techniques for identifying cyanobacteria assemblages across hundreds of lakes

Accurately identifying the species present in an ecosystem is vital to lake managers and successful bioassessment programs. This is particularly important when monitoring cyanobacteria, as numerous taxa produce toxins and can have major negative impacts on aquatic ecosystems. Increasingly, DNA-based techniques such as metabarcoding are being used for measuring aquatic biodiversity, as they could accelerate processing time, decrease costs and reduce some of the biases associated with traditional light microscopy. Despite the continuing use of traditional microscopy and the growing use of DNA metabarcoding to identify cyanobacteria assemblages, methodological comparisons between the two approaches have rarely been reported from a wide suite of lake types. Here, we compare planktonic cyanobacteria assemblages generated by inverted light microscopy and DNA metabarcoding from a 379-lake dataset spanning a longitudinal and trophic gradient. We found moderate levels of congruence between methods at the broadest taxonomic levels (i.e., Order, RV=0.40, p < 0.0001). This comparison revealed distinct cyanobacteria communities from lakes of different trophic states, with Microcystis, Aphanizomenon and Dolichospermum dominating with both methods in eutrophic and hypereutrophic sites. This finding supports the use of either method when monitoring eutrophication in lake surface waters. The biggest difference between the two methods was the detection of picocyanobacteria, which are typically underestimated by light microscopy. This reveals that the communities generated by each method currently are complementary as opposed to identical and promotes a combined-method strategy when monitoring a range of trophic systems. For example, microscopy can provide measures of cyanobacteria biomass, which are critical data in managing lakes. Going forward, we believe that molecular genetic methods will be increasingly adopted as reference databases are routinely updated with more representative sequences and will improve as cyanobacteria taxonomy is resolved with the increase in available genetic information.

Increases in Picocyanobacteria Abundance in Agriculturally Eutrophic Pampean Lakes Inferred from Historical Records of Secchi Depth and Chlorophyll-a

Phytoplankton size structure has profound consequences on food-web organization and energy transfer. Presently, picocyanobacteria (size < 2 µm) represent a major fraction of the autotrophic plankton of Pampean lakes. Glyphosate is known to stimulate the development of picocyanobacteria capable of degrading the herbicide. Due to the worldwide adoption of glyphosate-resistant crops, herbicide usage has increased sharply since the mid-1990s. Unfortunately, there are very few studies (none for the Pampa region) reporting picocyanobacteria abundance before 2000. The proliferation of µm sized particles should decrease Secchi disc depth (ZSD). Therefore ZSD, conditional to chlorophyll-a, may serve as an indicator of picocyanobacteria abundance. We use generalized additive models (GAMs) to analyze a “validation” dataset consisting of 82 records of ZSD, chlorophyll-a, and picocyanobacteria abundance from two Pampean lakes surveys (2009 and 2015). In support of the hypothesis, ZSD was negatively related to picocyanobacteria after accounting for the effect of chlorophyll-a. We then fitted a “historical” dataset using hierarchical GAMs to compare ZSD conditional to chlorophyll-a, before and after 2000. We estimated that ZSD levels during 2000–2021 were, on average, only about half as deep as those during 1980–1999. We conclude that the adoption of glyphosate-resistant crops has stimulated outbreaks of picocyanobacteria populations, resulting in lower water transparency.

The dynamics of picocyanobacteria from a hypereutrophic shallow lake is affected by light-climate and small-bodied zooplankton: a 10-year cytometric time-series analysis

In aquatic systems, an interplay between bottom-up and top-down processes determines the dynamic of picocyanobacteria (Pcy) abundance and community structure. Here, we analyzed a 10-year time series (sampled fortnightly) from a hypereutrophic turbid shallow lake located within the Pampa Region of South America, generating the first long-term record of freshwater Pcy from the Southern Hemisphere. We used a cytometric approach to study Pcy community, and focused on its relations with nutrient and light conditions (bottom-up) and potential grazers (top-down). A novel Pcy abundance seasonality with winter maximums was observed for years with relatively stable hydrological levels, related with decreased abundance of seasonal rotifers during colder seasons. Pcy showed lower abundance and higher cytometric alpha diversity during summer, probably due to a strong predation exerted by rotifers. In turn, a direct effect of the non-seasonal small cladocerans Bosmina spp. decreased Pcy abundance and induced a shift from single-cell Pcy into aggregated forms. This structuring effect of Bosmina spp. was further confirmed by Pcy cytometric (dis)similarity analyses from the time series and in situ experimental data. Remarkably, Pcy showed acclimatization to underwater light variations, resembling the relevance of light in this turbid system.