Metabarcoding Reveals Temporal Patterns of Community Composition and Realized Thermal Niches of Thalassiosira Spp. (Bacillariophyceae) from the Narragansett Bay Long-Term Plankton Time Series

Image-derived indicators of phytoplankton community responses to Pseudo-nitzschia blooms

Phytoplankton populations in the natural environment interact with each other. Despite rising global concern with Pseudo-nitzschia blooms, which can produce the potent neurotoxin domoic acid, we still do not fully understand how other phytoplankton genera respond to the presence of Pseudo-nitzschia. Here, we used a 4-year high-resolution imaging dataset for 9 commonly found phytoplankton genera in Narragansett Bay, alongside environmental data, to identify potential interactions between phytoplankton genera and their response to elevated Pseudo-nitzschia abundance. Our results indicate that Pseudo-nitzschia tends to bloom either concurrently with or right after other phytoplankton genera. Such bloom periods coincide with higher water temperatures and lower salinity. Pseudo-nitzschia image abundance tends to increase the most from March-May and peaks during May-Jun, whereas the image-derived biovolume and width of Pseudo-nitzschia chains increase the most during Jan-Feb. For most phytoplankton genera, their relationship with Pseudo-nitzschia abundance is noticeably different from their relationship with Pseudo-nitzschia image features. Despite the complexity in the phytoplankton community, our analysis suggests several ecological indicators that may be used to determine the risk of harmful algal blooms.

Surface currents shape protist community structure across the Indo-Pacific

Biogeographic structure in marine protist communities is shaped by a combination of dispersal potential and environmental selection. High-throughput sequencing and global sampling efforts have helped better resolve the composition and functions of these communities in the world's oceans using both molecular and visual methods. However, molecular barcoding data are critically lacking across the Indo-Pacific, a region widely considered the epicenter of marine biodiversity. To fill this gap, we characterized protist communities in four sampling regions across Indonesia that represent the latitudinal, longitudinal, and human population gradients of the region: Lombok, Wakatobi, Misool, and Waigeo. We show high spatial structuring in marine protist communities across Indonesia, and biotic factors appear to play little role in driving this observed structure. Our results appear to be driven by abiotic factors linked to surface current patterns across the Indo-Pacific as a result of: (1) a choke point in circulation at the Indonesian Throughflow leading to low diatom diversity in Lombok, Wakatobi, and Misool; (2) an increase in nutrient availability at the edge of the Halmahera Eddy in Waigeo, leading to an increase in diatom diversity; and/or (3) seasonal variations in protist communities in line with shifts in velocity of the Indonesian Throughflow. Overall, our results highlight the importance of abiotic factors in shaping protist communities on broad geographic scales over biotic, top-down pressures, such as grazing from higher trophic levels.

Metabarcoding reveals potentially mixotrophic flagellates and picophytoplankton as key groups of phytoplankton in the Elbe estuary

In estuaries, phytoplankton are faced with strong environmental forcing (e.g. high turbidity, salinity gradients). Taxa that appear under such conditions may play a critical role in maintaining food webs and biological carbon pumping, but knowledge about estuarine biota remains limited. This is also the case in the Elbe estuary where the lower 70 km of the water body are largely unexplored. In the present study, we investigated the phytoplankton composition in the Elbe estuary via metabarcoding. Our aim was to identify key taxa in the unmonitored reaches of this ecosystem and compare our results from the monitored area with available microscopy data. Phytoplankton communities followed distinct seasonal and spatial patterns. Community composition was similar across methods. Contributions of key classes and genera were correlated to each other (p < 0.05) when obtained from reads and biovolume (R2 = 0.59 and 0.33, respectively). Centric diatoms (e.g. Stephanodiscus) were the dominant group - comprising on average 55 % of the reads and 66-69 % of the biovolume. However, results from metabarcoding imply that microscopy underestimates the prevalence of picophytoplankton and flagellates with a potential for mixotrophy (e.g. cryptophytes). This might be due to their small size and sensitivity to fixation agents. We argue that mixotrophic flagellates are ecologically relevant in the mid to lower estuary, where, e.g., high turbidity render living conditions rather unfavorable, and skills such as phagotrophy provide fundamental advantages. Nevertheless, further findings - e.g. important taxa missing from the metabarcoding dataset - emphasize potential limitations of this method and quantitative biases can result from varying numbers of gene copies in different taxa. Further research should address these methodological issues but also shed light on the causal relationship of taxa with the environmental conditions, also with respect to active mixotrophic behavior.

Long-term declines in chlorophyll a and variable phenology revealed by a 60-year estuarine plankton time series

Long-term ecological time series provide a unique perspective on the emergent properties of ecosystems. In aquatic systems, phytoplankton form the base of the food web and their biomass, measured as the concentration of the photosynthetic pigment chlorophyll a (chl a), is an indicator of ecosystem quality. We analyzed temporal trends in chl a from the Long-Term Plankton Time Series in Narragansett Bay, Rhode Island, USA, a temperate estuary experiencing long-term warming and changing anthropogenic nutrient inputs. Dynamic linear models were used to impute and model environmental variables (1959 to 2019) and chl a concentrations (1968 to 2019). A long-term chl a decrease was observed with an average decline in the cumulative annual chl a concentration of 49% and a marked decline of 57% in winter-spring bloom magnitude. The long-term decline in chl a concentration was directly and indirectly associated with multiple environmental factors that are impacted by climate change (e.g., warming temperatures, water column stratification, reduced nutrient concentrations) indicating the importance of accounting for regional climate change effects in ecosystem-based management. Analysis of seasonal phenology revealed that the winter-spring bloom occurred earlier, at a rate of 4.9 ± 2.8 d decade-1. Finally, the high degree of temporal variation in phytoplankton biomass observed in Narragansett Bay appears common among estuaries, coasts, and open oceans. The commonality among these marine ecosystems highlights the need to maintain a robust set of phytoplankton time series in the coming decades to improve signal-to-noise ratios and identify trends in these highly variable environments.

Rich diversity and active spatial-temporal dynamics of Thalassiosira species revealed by time-series metabarcoding analysis

Thalassiosira is a species-rich genus in Bacillariophyta that not only contributes positively as primary producer, but also poses negative impacts on ecosystems by causing harmful algal blooms. Although taxonomical studies have identified a large number of Thalassiosira species, however, the composition of Thalassiosira species and their geographical distribution in marine ecosystems were not well understood due primarily to the lack of resolution of morphology-based approaches used previously in ecological expeditions. In this study, we systematically analyzed the composition and spatial-temporal dynamic distributions of Thalassiosira in the model marine ecosystem Jiaozhou Bay by applying metabarcoding analysis. Through analyzing samples collected monthly from 12 sampling sites, 14 Thalassiosira species were identified, including five species that were not previously reported in Jiaozhou Bay, demonstrating the resolution and effectiveness of metabarcoding analysis in ecological research. Many Thalassiosira species showed prominent temporal preferences in Jiaozhou Bay, with some displaying spring-winter preference represented by Thalassiosira tenera, while others displaying summer-autumn preference represented by Thalassiosira lundiana and Thalassiosira minuscula, indicating that the temperature is an important driving factor in the temporal dynamics. The application of metabarcoding analysis, equipped with appropriate molecular markers with high resolution and high specificity and databases of reference molecular marker sequences for potential all Thalassiosira species, will revolutionize ecological research of Thalassiosira species in Jiaozhou Bay and other marine ecosystems.

Oligotrophic waters of the Northwest Atlantic support taxonomically diverse diatom communities that are distinct from coastal waters

Diatoms are important components of the marine food web and one of the most species-rich groups of phytoplankton. The diversity and composition of diatoms in eutrophic nearshore habitats have been well documented due to the outsized influence of diatoms on coastal ecosystem functioning. In contrast, patterns of both diatom diversity and community composition in offshore oligotrophic regions where diatom biomass is low have been poorly resolved. To compare the diatom diversity and community composition in oligotrophic and eutrophic waters, diatom communities were sampled along a 1,250 km transect from the oligotrophic Sargasso Sea to the coastal waters of the northeast US shelf. Diatom community composition was determined by amplifying and sequencing the 18S rDNA V4 region. Of the 301 amplicon sequence variants (ASVs) identified along the transect, the majority (70%) were sampled exclusively from oligotrophic waters of the Gulf Stream and Sargasso Sea and included the genera Bacteriastrum, Haslea, Hemiaulus, Pseudo-nitzschia, and Nitzschia. Diatom ASV richness did not vary along the transect, indicating that the oligotrophic Sargasso Sea and Gulf Stream are occupied by a diverse diatom community. Although ASV richness was similar between oligotrophic and coastal waters, diatom community composition in these regions differed significantly and was correlated with temperature and phosphate, two environmental variables known to influence diatom metabolism and geographic distribution. In sum, oligotrophic waters of the western North Atlantic harbor diverse diatom assemblages that are distinct from coastal regions, and these open ocean diatoms warrant additional study, as they may play critical roles in oligotrophic ecosystems.

Dual thermal ecotypes coexist within a nearly genetically identical population of the unicellular marine cyanobacterium Synechococcus

The extent and ecological significance of intraspecific functional diversity within marine microbial populations is still poorly understood, and it remains unclear if such strain-level microdiversity will affect fitness and persistence in a rapidly changing ocean environment. In this study, we cultured 11 sympatric strains of the ubiquitous marine picocyanobacterium Synechococcus isolated from a Narragansett Bay (RI) phytoplankton community thermal selection experiment. Thermal performance curves revealed selection at cool and warm temperatures had subdivided the initial population into thermotypes with pronounced differences in maximum growth temperatures. Curiously, the genomes of all 11 isolates were almost identical (average nucleotide identities of >99.99%, with >99% of the genome aligning) and no differences in gene content or single nucleotide variants were associated with either cool or warm temperature phenotypes. Despite a very high level of genomic similarity, sequenced epigenomes for two strains showed differences in methylation on genes associated with photosynthesis. These corresponded to measured differences in photophysiology, suggesting a potential pathway for future mechanistic research into thermal microdiversity. Our study demonstrates that present-day marine microbial populations can harbor cryptic but environmentally relevant thermotypes which may increase their resilience to future rising temperatures.

Methodological 'lenses' influence the characterization of phytoplankton dynamics in a coastal upwelling ecosystem

New technologies enable the opportunity to improve our monitoring and understanding of marine phytoplankton communities. However, careful consideration for how different methodological approaches, or 'lenses', influence our interpretation of phytoplankton ecology is important, particularly when drawing conclusions about change over time or space. Using both high-throughput 18S rRNA gene sequencing and microscopy, we explored how phytoplankton community structure varied over the course of a year within a nearshore semi-enclosed coastal embayment along the Central Coast of California. The seasonal shift in the relative community dominance (i.e., diatom vs. dinoflagellate dominance) was captured in the microscopy results but not effectively captured in the molecular-based findings. However, the molecular approach explained more of the variability in composition across seasons relative to the microscopy approach. Temporal dynamics of specific bloom-forming taxa also differed between the molecular and microscopy results. Overall, the observed differences between the molecular- and microscopy-derived characterization of phytoplankton dynamics suggest that the approaches are best suited to answer different research questions. Moreover, the approaches complement each other for a more comprehensive perspective of a coastal phytoplankton ecosystem. Therefore, identifying the biases of each approach within natural communities is necessary to effectively and accurately characterize phytoplankton communities.

Emerging harmful algal blooms caused by distinct seasonal assemblages of a toxic diatom

Diatoms in the Pseudo-nitzschia genus produce the neurotoxin domoic acid. Domoic acid bioaccumulates in shellfish, causing illness in humans and marine animals upon ingestion. In 2017, high domoic acid levels in shellfish meat closed shellfish harvest in Narragansett Bay, Rhode Island for the first and only time in history, although abundant Pseudo-nitzschia have been observed for over 60 years. To investigate whether an environmental factor altered endemic Pseudo-nitzschia physiology or new domoic acid-producing strain(s) were introduced to Narragansett Bay, we conducted weekly sampling from 2017 to 2019 and compared closure samples. Plankton-associated domoic acid was quantified by LC-MS/MS and Pseudo-nitzschia spp. were identified using a taxonomically improved high-throughput rDNA sequencing approach. Comparison with environmental data revealed a detailed understanding of domoic acid dynamics and seasonal multi-species assemblages. Plankton-associated domoic acid was low throughout 2017-2019, but recurred in fall and early summer maxima. Fall domoic acid maxima contained known toxic species as well as a novel Pseudo-nitzschia genotype. Summer domoic acid maxima included fewer species but also known toxin producers. Most 2017 closure samples contained the particularly concerning toxic species, P. australis, which also appeared infrequently during 2017-2019. Recurring Pseudo-nitzschia assemblages were driven by seasonal temperature changes, and plankton-associated domoic acid correlated with low dissolved inorganic nitrogen. Thus, the Narragansett Bay closures were likely caused by both resident assemblages that become toxic depending on nutrient status as well as the episodic introductions of toxic species from oceanographic and climatic shifts.

Co-occurrence networks reveal the central role of temperature in structuring the plankton community of the Thau Lagoon

To identify the environmental factors that drive plankton community composition and structure in coastal waters, a shallow northwestern Mediterranean lagoon was monitored from winter to spring in two contrasting years. The campaign was based on high-frequency recordings of hydrological and meteorological parameters and weekly samplings of nutrients and the plankton community. The collected data allowed the construction of correlation networks, which revealed that water temperature was the most important factor governing community composition, structure and succession at different trophic levels, suggesting its ubiquitous food web control. Temperature favoured phytoplanktonic flagellates (Cryptophyceae, Chrysophyceae, and Chlorophyceae) and ciliates during winter and early spring. In contrast, it favoured Bacillariophyceae, dinoflagellates, phytoplankton < 6 µm and aloricate Choreotrichida during spring. The secondary factors were light, which influenced phytoplankton, and wind, which may regulate turbidity and the nutrient supply from land or sediment, thus affecting benthic species such as Nitzschia sp. and Uronema sp. or salinity-tolerant species such as Prorocentrum sp. The central role of temperature in structuring the co-occurrence network suggests that future global warming could deeply modify plankton communities in shallow coastal zones, affecting whole-food web functioning.

Grazing Pressure Is Independent of Prey Size in a Generalist Herbivorous Protist: Insights from Experimental Temperature Gradients

Grazing by herbivorous protists contributes to structuring plankton communities through its effect on the growth, biomass, and competitiveness of prey organisms and also impacts the transfer of primary production towards higher trophic levels. Previous evidence shows that heterotrophic processes (grazing rates, g) are more sensitive to temperature than autotrophic ones (phytoplankton growth rates, μ) and also that small cells tend to be more heavily predated than larger ones; however, it remains unresolved how the interplay between changes in temperature and cell size modulates grazing pressure (i.e., g:μ ratio). We addressed this problem by conducting an experiment with four phytoplankton populations, from pico- to microphytoplankton, over a 12 °C gradient and in the presence/absence of a generalist herbivorous protist, Oxyrrhis marina. We found that highest g rates coincided with highest μ rates, which corresponded to intermediate cell sizes. There were no significant differences in either μ or g between the smallest and largest cell sizes considered. The g:μ ratio was largely independent of cell size and C:N ratios, and its thermal dependence was low although species-specific differences were large. We suggest that the similar g:μ found could be the consequence that the energetic demand imposed by rising temperatures would be a more important issue than the mechanical constriction to ingestion derived from prey cell size. Despite the difficulty of quantifying μ and g in natural planktonic communities, we suggest that the g:μ ratio is a key response variable to evaluate thermal sensitivity of food webs because it gives a more integrative view of trophic functioning than both rates separately.

Exploring biogeochemical and ecological redundancy in phytoplankton communities in the global ocean

Climate-change-induced alterations of oceanic conditions will lead to the ecological niches of some marine phytoplankton species disappearing, at least regionally. How will such losses affect the ecosystem and the coupled biogeochemical cycles? Here, we couch this question in terms of ecological redundancy (will other species be able to fill the ecological roles of the extinct species) and biogeochemical redundancy (can other species replace their biogeochemical roles). Prior laboratory and field studies point to a spectrum in the degree of redundancy. We use a global three-dimensional computer model with diverse planktonic communities to explore these questions further. The model includes 35 phytoplankton types that differ in size, biogeochemical function and trophic strategy. We run two series of experiments in which single phytoplankton types are either partially or fully eliminated. The niches of the targeted types were not completely reoccupied, often with a reduction in the transfer of matter from autotrophs to heterotrophs. Primary production was often decreased, but sometimes increased due to reduction in grazing pressure. Complex trophic interactions (such as a decrease in the stocks of a predator's grazer) led to unexpected reshuffling of the community structure. Alterations in resource utilization may cause impacts beyond the regions where the type went extinct. Our results suggest a lack of redundancy, especially in the 'knock on' effects on higher trophic levels. Redundancy appeared lowest for types on the edges of trait space (e.g. smallest) or with unique competitive strategies. Though highly idealized, our modelling findings suggest that the results from laboratory or field studies often do not adequately capture the ramifications of functional redundancy. The modelled, often counterintuitive, responses-via complex food web interactions and bottom-up versus top-down controls-indicate that changes in planktonic community will be key determinants of future ocean global change ecology and biogeochemistry.