Planktonic eukaryotes in the Chesapeake Bay: contrasting responses of abundant and rare taxa to estuarine gradients

Phytoplankton are important drivers of aquatic ecosystem function and environmental health. Their community compositions and distributions are directly impacted by environmental processes and human activities, including in the largest estuary in North America, the Chesapeake Bay. It is crucial to uncover how planktonic eukaryotes play fundamental roles as primary producers and trophic links and sustain estuarine ecosystems. In this study, we investigated the detailed community structure and spatiotemporal variations of planktonic eukaryotes in the Chesapeake Bay across space and time for three consecutive years. A clear seasonal and spatial shift of total, abundant, and rare planktonic eukaryotes was evident, and the pattern recurred interannually. Multiple harmful algal species have been identified in the Bay with varied distribution patterns, such as Karlodinium, Heterosigma akashiwo, Protoperidinium sp., etc. Compared to abundant taxa, rare subcommunities were more sensitive to environmental disturbance in terms of richness, diversity, and distribution. The combined effects of temporal variation (13.3%), nutrient availability (10.0%), and spatial gradients (8.8%) structured the distribution of eukaryotic microbial communities in the Bay. Similar spatiotemporal patterns between planktonic prokaryotes and eukaryotes suggest common mechanisms of adjustment, replacement, and species interaction for planktonic microbiomes under strong estuarine gradients. To our best knowledge, this work represents the first systematic study on planktonic eukaryotes in the Bay. A comprehensive view of the distribution of planktonic microbiomes and their interactions with environmental processes is critical in understanding the underlying microbial mechanisms involved in maintaining the stability, function, and environmental health of estuarine ecosystems.

IMPORTANCE

Deep sequencing analysis of planktonic eukaryotes in the Chesapeake Bay reveals high community diversity with many newly recognized phytoplankton taxa. The Chesapeake Bay planktonic eukaryotes show distinct seasonal and spatial variability, with recurring annual patterns of total, abundant, and rare groups. Rare taxa mainly contribute to eukaryotic diversity compared to abundant groups, and they are more sensitive to spatiotemporal variations and environmental filtering. Temporal variations, nutrient availability, and spatial gradients significantly affect the distribution of eukaryotic microbial communities. Similar spatiotemporal patterns in prokaryotes and eukaryotes suggest common mechanisms of adjustment, substitution, and species interactions in planktonic microbiomes under strong estuarine gradients. Interannually recurring patterns demonstrate that diverse eukaryotic taxa have well adapted to the estuarine environment with a long residence time. Further investigations of how human activities impact estuarine planktonic eukaryotes are critical in understanding their essential ecosystem roles and in maintaining environmental safety and public health.

Virus-to-prokaryote ratio in the Salar de Huasco and different ecosystems of the Southern hemisphere and its relationship with physicochemical and biological parameters

The virus-to-prokaryote ratio (VPR) has been used in many ecosystems to study the relationship between viruses and their hosts. While high VPR values indicate a high rate of prokaryotes' cell lysis, low values are interpreted as a decrease in or absence of viral activity. Salar de Huasco is a high-altitude wetland characterized by a rich microbial diversity associated with aquatic sites like springs, ponds, streams and a lagoon with variable physicochemical conditions. Samples from two ponds, Poza Rosada (PR) and Poza Verde (PV), were analyzed by epifluorescence microscopy to determine variability of viral and prokaryotic abundance and to calculate the VPR in a dry season. In addition, to put Salar de Huasco results into perspective, a compilation of research articles on viral and prokaryotic abundance, VPR, and metadata from various Southern hemisphere ecosystems was revised. The ecosystems were grouped into six categories: high-altitude wetlands, Pacific, Atlantic, Indian, and Southern Oceans and Antarctic lakes. Salar de Huasco ponds recorded similar VPR values (an average of 7.4 and 1.7 at PR and PV, respectively), ranging from 3.22 to 15.99 in PR. The VPR variability was associated with VA and chlorophyll a, when considering all data available for this ecosystem. In general, high-altitude wetlands recorded the highest VPR average (53.22 ± 95.09), followed by the Oceans, Southern (21.91 ± 25.72), Atlantic (19.57 ± 15.77) and Indian (13.43 ± 16.12), then Antarctic lakes (11.37 ± 15.82) and the Pacific Ocean (6.34 ± 3.79). Physicochemical variables, i.e., temperature, conductivity, nutrients (nitrate, ammonium, and phosphate) and chlorophyll a as a biological variable, were found to drive the VPR in the ecosystems analyzed. Thus, the viral activity in the Wetland followed similar trends of previous reports based on larger sets of metadata analyses. In total, this study highlights the importance of including viruses as a biological variable to study microbial temporal dynamics in wetlands considering their crucial role in the carbon budgets of these understudied ecosystems in the southern hemisphere.

Greenhouse gases and biogeochemical diel fluctuations in a high-altitude wetland

The landscapes of high-altitude wetland ecosystems are characterized by different kinds of aquatic sites, including ponds holding conspicuous microbial life. Here, we examined a representative pond of the wetland landscape for dynamics of greenhouse gases, and their association with other relevant biogeochemical conditions including diel shifts of microbial communities' structure and activity over two consecutive days. Satellite image analysis indicates that the area of ponds cover 238 of 381.3 Ha (i.e., 62.4%), representing a significant landscape in this wetland. Solar radiation, wind velocity and temperature varied daily and between the days sampled, influencing the biogeochemical dynamics in the pond, shifting the pond reservoir of inorganic versus dissolved organic nitrogen/phosphorus bioavailability, between day 1 and day 2. Day 2 was characterized by high dissolved organic nitrogen/phosphorus and N₂O accumulation. CH₄ presented a positive excess showing maxima at hours of high radiation during both days. The microbial community in the sediment was diverse and enriched in keystone active groups potentially related with GHG recycling including bacteria and archaea, such as Cyanobacteria, Verrucomicrobia, Rhodobacterales and Nanoarchaeaota (Woesearchaeia). Archaea account for the microbial community composition changes between both days and for the secondary productivity in the water measured during day 2. The results indicate that an intense recycling of organic matter occurs in the pond systems and that the activity of the microbial community is correlated with the availability of nutrients. Together, the above results indicate a net sink of CO₂ and N₂O, which has also been reported for other natural and artificial ponds. Overall, our two-day fluctuation study in a representative pond of a high-altitude wetland aquatic landscape indicates the need to explore in more detail the short-term besides the long-term biogeochemical variability in arid ecosystems of the Andes plateau, where wetlands are hotspots of life currently under high anthropogenic pressure.

Active microbiome structure and its association with environmental factors and viruses at different aquatic sites of a high-altitude wetland

Salar de Huasco is a high‐altitude wetland characterized by a highly diverse microbial life adapted to extreme climatic and environmental conditions. Our study aims to determine active microbial community structure changes within different aquatic sites and its relationship with environmental factors and viruses as potential drivers of diversification in different aquatic areas of this ecosystem. In this study, bacteria and archaea composition (16S rRNA subunit pyrolibraries) and picoplankton and viral abundance were determined at ponds, springs and lagoon sites of the wetland during wet and dry seasons (February and July 2012, respectively). In general, mixosaline waters (1,400–51,000 μS/cm) usually found in ponds and lagoon presented higher picoplanktonic abundances compared to freshwater (<800 μS/cm) spring sites, ranging from 1.07 × 10⁵ to 1.83 × 10⁷ cells/ml. Viral abundance and viral to picoplankton ratio (VPR) also presented greater values at ponds compared to spring sites, reaching up to 4.78 × 10⁸ viruses‐like particles and up to 351 for VPR. In general, ponds hold a higher microbial diversity and complexity associated also with the presence of microbial mats compared with water sources or lagoon (Shannon index H′ 2.6–3.9 vs. <2.0). A greater richness of archaea was also detected in ponds characterized by functional groups such as known methanogens and ammonia oxidizers, and uncultured groups. In total, our results indicate that among the different aquatic sites of the wetland, ponds presented a great microbial community diversification associated to a higher top‐down control by viruses which may influence nutrient and greenhouse gases cycling.

Distribution of greenhouse gases in hyper-arid and arid areas of northern Chile and the contribution of the high altitude wetland microbiome (Salar de Huasco, Chile)

Northern Chile harbors different bioclimatic zones including hyper-arid and arid ecosystems and hotspots of microbial life, such as high altitude wetlands, which may contribute differentially to greenhouse gases (GHG) such as carbon dioxide (CO₂), methane (CH₄) and nitrous oxide (N₂O). In this study, we explored ground level GHG distribution and the potential role of a wetland situated at 3800 m.a.s.l, and characterized by high solar radiation < 1600 W m^-2, extreme temperature ranges (-12 to 24 °C) and wind stress (< 17 m s^-1). The water source of the wetland is mainly groundwater springs, which generates streams and ponds surrounded by peatlands. These sites support a rich microbial aquatic life including diverse bacteria and archaea communities, which transiently form more complex structures, such as microbial mats. In this study, GHG were measured in the water and above ground level air at the wetland site and along an elevation gradient in different bioclimatic areas from arid to hyper-arid zones. The microbiome from the water and sediments was described by high-throughput sequencing 16S rRNA and rDNA genes. The results indicate that GHG at ground level were variable along the elevation gradient potentially associated with different bioclimatic zones, reaching high values at the high Andean steppe and variable but lower values in the Atacama Desert and at the wetland. The water areas of the wetland presented high concentrations of CH₄ and CO₂, particularly at the spring areas and in air bubbles below microbial mats. The microbial community was rich (> 40 phyla), including archaea and bacteria potentially active in the different matrices studied (water, sediments and mats). Functional microbial groups associated with GHG recycling were detected at low frequency, i.e., < 2.5% of total sequences. Our results indicate that hyper-arid and arid areas of northern Chile are sites of GHG exchange associated with various bioclimatic zones and particularly in aquatic areas of the wetland where this ecosystem could represent a net sink of N₂O and a source for CH₄ and CO₂.

Bacterial Active Community Cycling in Response to Solar Radiation and Their Influence on Nutrient Changes in a High-Altitude Wetland

Microbial communities inhabiting high-altitude spring ecosystems are subjected to extreme changes in solar irradiance and temperature throughout the diel cycle. Here, using 16S rRNA gene tag pyrosequencing (cDNA) we determined the composition of actively transcribing bacteria from spring waters experimentally exposed through the day (morning, noon, and afternoon) to variable levels of solar radiation and light quality, and evaluated their influence on nutrient recycling. Solar irradiance, temperature, and changes in nutrient dynamics were associated with changes in the active bacterial community structure, predominantly by Cyanobacteria, Verrucomicrobia, Proteobacteria, and 35 other Phyla, including the recently described Candidate Phyla Radiation (e.g., Parcubacteria, Gracilibacteria, OP3, TM6, SR1). Diversity increased at noon, when the highest irradiances were measured (3.3–3.9 H′, 1125 W m^-2) compared to morning and afternoon (0.6–2.8 H′). This shift was associated with a decrease in the contribution to pyrolibraries by Cyanobacteria and an increase of Proteobacteria and other initially low frequently and rare bacteria phyla (< 0.5%) in the pyrolibraries. A potential increase in the activity of Cyanobacteria and other phototrophic groups, e.g., Rhodobacterales, was observed and associated with UVR, suggesting the presence of photo-activated repair mechanisms to resist high levels of solar radiation. In addition, the percentage contribution of cyanobacterial sequences in the afternoon was similar to those recorded in the morning. The shifts in the contribution by Cyanobacteria also influenced the rate of change in nitrate, nitrite, and phosphate, highlighted by a high level of nitrate accumulation during hours of high radiation and temperature associated with nitrifying bacteria activity. We did not detect ammonia or nitrite oxidizing bacteria in situ, but both functional groups (Nitrosomona and Nitrospira) appeared mainly in pyrolibraries generated from dark incubations. In total, our results reveal that both the structure and the diversity of the active bacteria community was extremely dynamic through the day, and showed marked shifts in composition that influenced nutrient recycling, highlighting how abiotic variation affects potential ecosystem functioning.

Detection and phylogenetic analysis of infectious pancreatic necrosis virus in Chile

Infectious pancreatic necrosis virus (IPNV) is the etiological agent of a highly contagious disease that is endemic to salmon farming in Chile and causes great economic losses to the industry. Here we compared different diagnostic methods to detect IPNV in field samples, including 3 real-time reverse transcription PCR (qRT-PCR) assays, cell culture isolation, and indirect fluorescent antibody test (IFAT). Additionally, we performed a phylogenetic analysis to investigate the genogroups prevailing in Chile, as well as their geographic distribution and virulence. The 3 qRT-PCR assays used primers that targeted regions of the VP2 and VP1 genes of the virus and were tested in 46 samples, presenting a fair agreement within their results. All samples were positive for at least 2 of the qRT-PCR assays, 29 were positive for cell culture, and 23 for IFAT, showing less sensitivity for these latter 2 methods. For the phylogenetic analysis, portions of 1180 and 523 bp of the VP2 region of segment A were amplified by RT-PCR, sequenced and compared with sequences from reference strains and from isolates reported by previous studies carried out in Chile. Most of the sequenced isolates belonged to genogroup 5 (European origin), and 5 were classified within genogroup 1 (American origin). Chilean isolates formed clusters within each of the genogroups found, evidencing a clear differentiation from the reference strains. To our knowledge, this is the most extensive study completed for IPNV in Chile, covering isolates from sea- and freshwater salmon farms and showing a high prevalence of this virus in the country.

ULTRASTRUCTURAL CHARACTERIZATION OF THE LYTIC CYCLE OF AN INTRANUCLEAR VIRUS INFECTING THE DIATOM CHAETOCEROS CF. WIGHAMII(BACILLARIOPHYCEAE) FROM CHESAPEAKE BAY, USA(1)

Numerous microalgal species are infected by viruses that have the potential to control phytoplankton dynamics by reducing host populations, preventing bloom formation, or causing the collapse of blooms. Here we describe a virus infecting the diatom Chaetoceros cf. wighamii Brightw. from the Chesapeake Bay. To characterize the morphology and lytic cycle of this virus, we conducted a time-course experiment, sampling every 4 h over 72 h following viral inoculation. In vivo fluorescence began to decline 16 h after inoculation and was reduced to <19% of control cultures by the end of experiment. TEM confirmed infection within the first 8 h of inoculation, as indicated by the presence of virus-like particles (VLP) in the nuclei. VLP were present in two different arrangements: rod-like structures that appeared in cross-section as paracrystalline arrays of hexagonal-shaped profiles measuring 12 ± 2 nm in diameter and uniformly electron-dense hexagonal-shaped particles measuring ∼ 22-28 nm in diameter. Nuclei containing paracrystalline arrays were most prevalent early in the infection cycle, while cells containing VLP increased and then declined toward the end of the cycle. The proportion of nuclei containing both paracrystalline arrays and VLP remained relatively constant. This pattern suggests that rod-like paracrystalline arrays fragmented to produce icosahedral VLP. C. cf. wighamii nuclear inclusion virus (CwNIV) is characterized by a high burst size (averaged 26,400 viruses per infected cell) and fast generation time that could have ecological implications on C. cf. wighamii population control.