MOLECULAR DIVERSITY OF MARINE PHYTOPLANKTON COMMUNITIES BASED ON KEY FUNCTIONAL GENES(1)

Differential association of key bacterial groups with diatoms and Phaeocystis spp. during spring blooms in the Southern Ocean

Interactions between phytoplankton and heterotrophic bacteria significantly influence the cycling of organic carbon in the ocean, with many of these interactions occurring at the micrometer scale. We explored potential associations between specific phytoplankton and bacteria in two size fractions, 0.8-3 µm and larger than 3 µm, at three naturally iron-fertilized stations and one high nutrient low chlorophyll station in the Southern Ocean. The composition of phytoplankton and bacterial communities was determined by sequencing the rbcL gene and 16S rRNA gene from DNA and RNA extracts, which represent presence and potential activity, respectively. Diatoms, particularly Thalassiosira, contributed significantly to the DNA sequences in the larger size fractions, while haptophytes were dominant in the smaller size fraction. Correlation analysis between the most abundant phytoplankton and bacterial operational taxonomic units revealed strong correlations between Phaeocystis and picoeukaryotes with SAR11, SAR116, Magnetospira, and Planktomarina. In contrast, most Thalassiosira operational taxonomic units showed the highest correlations with Polaribacter, Sulfitobacteria, Erythrobacter, and Sphingobium, while Fragilariopsis, Haslea, and Thalassionema were correlated with OM60, Fluviicola, and Ulvibacter. Our in-situ observations suggest distinct associations between phytoplankton and bacterial taxa, which could play crucial roles in nutrient cycling in the Southern Ocean.

Spatial and temporal changes in the assembly mechanism and co-occurrence network of the chromophytic phytoplankton communities in coastal ecosystems under anthropogenic influences

As a typical semiclosed coastal sea area in China, the ecological environment of Bohai Bay has been significantly disturbed by human activities. As primary producers, the chromophytic phytoplankton are the basis of ecosystems, especially in coastal ecosystems, and changes in the chromophytic phytoplankton community can affect the stability of the entire ecosystem. In this study, we investigated the effects of the human activity-induced spatial and temporal environmental heterogeneity on the community composition, diversity, assembly mechanisms, and co-occurrence networks of chromophytic phytoplankton in Bohai Bay during the wet season and the dry season. The results showed that in both seasons, there was obvious environmental heterogeneity between the nearshore area and the offshore area, and the nearshore areas were more affected by human disturbance. Although higher diversity was supported by the abundance of nutrients in nearshore areas, co-occurrence network analysis revealed that the chromophytic phytoplankton were less closely connected to each other in nearshore areas than in offshore areas due to chemical oxygen demand (COD), eutrophication index (EI), and dissolved inorganic nitrogen (DIN). The nearshore network was less stable than the offshore co-occurrence network in both seasons, which was related to the concentration of dissolved oxygen and COD. Both stochastic and deterministic processes dominated the assembly of the chromophytic phytoplankton communities, with different importance rankings of stochastic and deterministic processes in the nearshore and offshore areas. Drift dominated the assembly of the communities in nearshore areas, while variable selection dominated the assembly of the communities in offshore areas. DIN, EI, and COD, rather than geographic distance, were the main environmental factors affecting the phylogenetic turnover of the chromophytic phytoplankton. Our study showed that environmental heterogeneity caused by human disturbance had a greater impact on the chromophytic phytoplankton communities in Bohai Bay than natural factors such as temperature and salinity.

Deciphering the diversity and distribution of chromophytic phytoplankton in the Bohai Sea and the Yellow Sea via RuBisCO genes (rbcL)

Present study investigated composition and distribution of chromophytic phytoplankton in the Bohai Sea (BS) and the Yellow Sea (YS) by using rbcL genes. Bacillariophyceae, Haptophyceae and Pelagophyceae were the most abundant phytoplankton groups. Distinct phytoplankton communities were observed in the BS and the YS: offshore stations were dominated by bloom forming genera Thalassiosira and Skeletonema, while brown tide-forming species including Chrysochromulina spp. and Aureococcus anophagefferens were commonly found in the nearshore areas. Redundancy analysis showed that phosphate, temperature and silicic acid play key roles in structuring chromophytic phytoplankton, such as phytoplankton at nearshore stations were affected by nutrient runoff from adjacent rivers (Yellow River). Anthropogenic activities in the Bohai Sea and seasonal circulation of ocean currents may also contribute to shaping chromophytic phytoplankton communities. This study provides data support and foundational observations of chromophytic phytoplankton in the BS and the YS, and their responses to environmental gradients and human activities.

The Enzymology of Ocean Global Change

A small subset of marine microbial enzymes and surface transporters have a disproportionately important influence on the cycling of carbon and nutrients in the global ocean. As a result, they largely determine marine biological productivity and have been the focus of considerable research attention from microbial oceanographers. Like all biological catalysts, the activity of these keystone biomolecules is subject to control by temperature and pH, leaving the crucial ecosystem functions they support potentially vulnerable to anthropogenic environmental change. We summarize and discuss both consensus and conflicting evidence on the effects of sea surface warming and ocean acidification for five of these critical enzymes [carbonic anhydrase, ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), nitrogenase, nitrate reductase, and ammonia monooxygenase] and one important transporter (proteorhodopsin). Finally, we forecast how the responses of these few but essential biocatalysts to ongoing global change processes may ultimately help to shape the microbial communities and biogeochemical cycles of the future greenhouse ocean.

Distribution of Chromophytic Phytoplankton in the Eddy-Induced Upwelling Region of the West Pacific Ocean Revealed Using rbcL Genes

Marine chromophytic phytoplankton are a diverse group of algae and contribute significantly to the total oceanic primary production. However, the spatial distribution of chromophytic phytoplankton is understudied in the West Pacific Ocean (WPO). In this study, we have investigated the community structure and spatial distribution of chromophytic phytoplankton using RuBisCO genes (Form ID rbcL). Our results showed that Haptophyceae, Pelagophyceae, Cyanophyceae, Xanthophyceae, and Bacillariophyceae were the dominant groups. Further, chromophytic phytoplankton can be distinguished between upwelling and non-upwelling zones of the WPO. Surface and 75 m depths of a non-upwelling area were dominated by Prochlorococcus strains, whereas chromophytic phytoplankton were homogenously distributed at the surface layer in the upwelling zone. Meanwhile, Pelagomonas-like sequences were dominant at DCM (75 m) and 150 m depths of the upwelling zone. Non-metric multidimensional scaling (NMDS) analysis did not differentiate between chromophytic phytoplankton in the upwelling and non-upwelling areas, however, it showed clear trends of them at different depths. Further, redundancy analysis (RDA) showed the influence of physicochemical parameters on the distribution of chromophytic phytoplankton. Along with phosphate (p < 0.01), temperature and other dissolved nutrients were important in driving community structure. The upwelling zone was impacted by a decrease in temperature, salinity, and re-supplement of nutrients, where Pelagomonas-like sequences outnumbered other chromophytic groups presented.

The spatiotemporal contribution of the phytoplankton community and environmental variables to the carbon sequestration potential in an urban river

The phytoplankton (internal driving forces) and environmental variables that affect complex biochemical reactions (external driving forces) play an important role in regulating photosynthetic carbon fixation. Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) exists in various phytoplankton species and is an important enzyme in the photosynthetic process. To investigate the phytoplankton composition (internal driving forces), we selected the functional gene of the Rubisco large subunit (rbcL) as the target gene for this study. Phytoplankton gross primary productivity was measured using light and dark biological oxygen demand bottles to assess the carbon sequestration potential. The fundamental environmental indicators were determined to analyze the mechanisms that drive the carbon fixation process. The correlation results indicated that green algae were only controlled by nitrate, and that diatoms were positively correlated with phosphate. The cluster analysis results demonstrated that nitrite was the major driver controlling phytoplankton primary productivity. During the wet seasons (spring and summer), the contribution of the planktonic community respiration to the carbon sequestration potential was higher than net primary productivity (NPP), followed by dissolved organic carbon and nitrate. During the dry season (autumn), NPP, total nitrogen, and nitrite ranked highest in terms of carbon sequestration potential. The contributions of green algae and diatoms to the carbon sequestration potential were temporally higher than those of cyanobacteria. The maximum carbon sequestration potential occurred during autumn because of diatom production and the function of phosphate, whereas the minimum carbon sequestration potential occurred in summer. Spatially, the upstream carbon sequestration potential was higher compared with downstream because of the effect (contribution) of cyanobacteria (Phormidium), diatoms (Surirella solea and Thalassiosira pseudonana), and environmental variable (nitrite). These findings provide a better understanding of the underlying mechanisms of phytoplankton productivity and the influences of environmental variables on carbon sequestration in urban river ecosystems.

Diversity and Spatial Distribution of Chromophytic Phytoplankton in the Bay of Bengal Revealed by RuBisCO Genes (rbcL)

Phytoplankton are the basis of primary production and play important roles in regulating energy export in marine ecosystems. Compared to other regions, chromophytic phytoplankton are considerably understudied in the Bay of Bengal (BOB). Here, we investigated community structure and spatial distribution of chromophytic phytoplankton in the BOB by using RuBisCO genes (Form ID rbcL). High throughput sequencing of rbcL genes revealed that diatoms, cyanobacteria (Cyanophyceae), Pelagophyceae, Haptophyceae, Chrysophyceae, Eustigamatophyceae, Xanthophyceae, Cryptophyceae, Dictyochophyceae, and Pinguiophyceae were the most abundant groups recovered in the BOB. Abundances and distribution of diatoms and Pelagophyceae were further verified using quantitative PCR analyses which showed the dominance of these groups near the Equator region (p < 0.01) where upwelling was likely the source of nutrients. Further, redundancy analysis demonstrated that temperature was an important environmental driver in structuring distributions of Cyanophyceae and dominant chromophytic phytoplankton. Morphological identification and quantification confirmed the dominance of diatoms, and also detected other cyanobacteria and dinoflagellates that were missing in our molecular characterizations. Pearson’s correlations of these morphologically identified phytoplankton with environmental gradients also indicated that nutrients and temperature were key variables shaping community structure. Combination of molecular characterization and morphological identification provided a comprehensive overview of chromophytic phytoplankton. This is the first molecular study of chromophytic phytoplankton accomplished in the BOB, and our results highlight a combination of molecular analysis targeting rbcL genes and microscopic detection in examining phytoplankton composition and diversity.

Nitrate Consumers in Arctic Marine Eukaryotic Communities: Comparative Diversities of 18S rRNA, 18S rRNA Genes, and Nitrate Reductase Genes

For photosynthetic microbial eukaryotes, the rate-limiting step in NO₃ ^- assimilation is its reduction to nitrite (NO₂ ^-), which is catalyzed by assimilatory nitrate reductase (NR). Oceanic productivity is primarily limited by available nitrogen and, although nitrate is the most abundant form of available nitrogen in oceanic waters, little is known about the identity of microbial eukaryotes that take up nitrate. This lack of knowledge is especially severe for ice-covered seas that are being profoundly affected by climate change. To address this, we examined the distribution and diversity of NR genes in the Arctic region by way of clone libraries and data mining of available metagenomes (total of 4.24 billion reads). We directly compared NR clone phylogenies with the V4 region of the 18S rRNA gene (DNA pool) and 18S rRNA (RNA pool) at two ice-influenced stations in the Canada Basin (Beaufort Sea). The communities from the two nucleic acid templates were similar at the level of major groups, and species identified by way of NR gene phylogeny and microscopy were a subset of the 18S results. Most NR genes from arctic clone libraries matched diatoms and chromist nanoflagellates, including novel clades, while the NR genes in arctic eukaryote metagenomes were dominated by chlorophyte NR, in keeping with the ubiquitous occurrence of Mamiellophyceae in the Arctic Ocean. Overall, these data suggest that a dynamic and mixed eukaryotic community utilizes nitrate across the Arctic region, and they show the potential utility of NR as a tool to identify ongoing changes in arctic photosynthetic communities.IMPORTANCE To better understand the diversity of primary producers in the Arctic Ocean, we targeted a nitrogen cycle gene, NR, which is required for phytoplankton to assimilate nitrate into organic forms of nitrogen macromolecules. We compared this to the more detailed taxonomy from ice-influenced stations using a general taxonomic gene (18S rRNA). NR genes were ubiquitous and could be classified as belonging to diatoms, dinoflagellates, other flagellates, chlorophytes, and unknown microbial eukaryotes, suggesting novel diversity of both species and metabolism in arctic phytoplankton.

Response of Spring Diatoms to CO2 Availability in the Western North Pacific as Determined by Next-Generation Sequencing

Next-generation sequencing (NGS) technologies have enabled us to determine phytoplankton community compositions at high resolution. However, few studies have adopted this approach to assess the responses of natural phytoplankton communities to environmental change. Here, we report the impact of different CO₂ levels on spring diatoms in the Oyashio region of the western North Pacific as estimated by NGS of the diatom-specific rbcL gene (DNA), which encodes the large subunit of RubisCO. We also examined the abundance and composition of rbcL transcripts (cDNA) in diatoms to assess their physiological responses to changing CO₂ levels. A short-term (3-day) incubation experiment was carried out on-deck using surface Oyashio waters under different pCO₂ levels (180, 350, 750, and 1000 μatm) in May 2011. During the incubation, the transcript abundance of the diatom-specific rbcL gene decreased with an increase in seawater pCO₂ levels. These results suggest that CO₂ fixation capacity of diatoms decreased rapidly under elevated CO₂ levels. In the high CO₂ treatments (750 and 1000 μatm), diversity of diatom-specific rbcL gene and its transcripts decreased relative to the control treatment (350 μatm), as well as contributions of Chaetocerataceae, Thalassiosiraceae, and Fragilariaceae to the total population, but the contributions of Bacillariaceae increased. In the low CO₂ treatment, contributions of Bacillariaceae also increased together with other eukaryotes. These suggest that changes in CO₂ levels can alter the community composition of spring diatoms in the Oyashio region. Overall, the NGS technology provided us a deeper understanding of the response of diatoms to changes in CO₂ levels in terms of their community composition, diversity, and photosynthetic physiology.

A comprehensive framework for functional diversity patterns of marine chromophytic phytoplankton using rbcL phylogeny

Marine chromophytes are taxonomically diverse group of algae and contribute approximately half of the total oceanic primary production. To understand the global patterns of functional diversity of chromophytic phytoplankton, robust bioinformatics and statistical analyses including deep phylogeny based on 2476 form ID rbcL gene sequences representing seven ecologically significant oceanographic ecoregions were undertaken. In addition, 12 form ID rbcL clone libraries were generated and analyzed (148 sequences) from Sundarbans Biosphere Reserve representing the world’s largest mangrove ecosystem as part of this study. Global phylogenetic analyses recovered 11 major clades of chromophytic phytoplankton in varying proportions with several novel rbcL sequences in each of the seven targeted ecoregions. Majority of OTUs was found to be exclusive to each ecoregion, whereas some were shared by two or more ecoregions based on beta-diversity analysis. Present phylogenetic and bioinformatics analyses provide a strong statistical support for the hypothesis that different oceanographic regimes harbor distinct and coherent groups of chromophytic phytoplankton. It has been also shown as part of this study that varying natural selection pressure on form ID rbcL gene under different environmental conditions could lead to functional differences and overall fitness of chromophytic phytoplankton populations.

Analysis of diversity of chromophytic phytoplankton in a mangrove ecosystem using rbcL gene sequencing

Phytoplankton forms the basis of primary production in mangrove environments. The phylogeny and diversity based on the amplification and sequencing of rbcL, the large subunit encoding the key enzyme ribulose-1, 5-bisphosphate carboxylase/oxygenase was investigated for improved understanding of the community structure and temporal trends of chromophytic eukaryotic phytoplankton assemblages in Sundarbans, the world's largest continuous mangrove. Diatoms (Bacillariophyceae) were by far the most frequently detected group in clone libraries (485 out of 525 clones), consistent with their importance as a major bloom-forming group. Other major chromophytic algal groups including Cryptophyceae, Haptophyceae, Pelagophyceae, Eustigmatophyceae, and Raphidophyceae which are important component of the assemblages were detected for the first time from Sundarbans based on rbcL approach. Many of the sequences from Sundarbans rbcL clone libraries showed identity with key bloom forming diatom genera namely Thalassiosira, Skeletonema and Nitzschia. Similarly, several rbcL sequences which were diatom-like were also detected highlighting the need to explore diatom communities from the study area. Some of the rbcL sequences detected from Sundarbans were ubiquitous in distribution showing 100% identities with uncultured rbcL sequences targeted previously from the Gulf of Mexico and California upwelling system that are geographically separated from study area. Novel rbcL lineages were also detected highlighting the need to culture and sequence phytoplankton from the ecoregion. Principal component analysis revealed that nitrate is an important variable that is associated with observed variation in phytoplankton assemblages (operational taxonomic units). This study applied molecular tools to highlight the ecological significance of diatoms, in addition to other chromophytic algal groups in Sundarbans.

Temporal Succession of Phytoplankton Assemblages in a Tidal Creek System of the Sundarbans Mangroves: An Integrated Approach

Sundarbans, the world's largest mangrove ecosystem, is unique and biologically diverse. A study was undertaken to track temporal succession of phytoplankton assemblages at the generic level (≥10 µm) encompassing 31 weeks of sampling (June 2010–May 2011) in Sundarbans based on microscopy and hydrological measurements. As part of this study, amplification and sequencing of type ID rbcL subunit of RuBisCO enzyme were also applied to infer chromophytic algal groups (≤10 µm size) from one of the study points. We report the presence of 43 genera of Bacillariophyta, in addition to other phytoplankton groups, based on microscopy. Phytoplankton cell abundance, which was highest in winter and spring, ranged between 300 and 27,500 cells/L during this study. Cell biovolume varied between winter of 2010 (90–35281.04 µm3) and spring-summer of 2011 (52–33962.24 µm3). Winter supported large chain forming diatoms, while spring supported small sized diatoms, followed by other algal groups in summer. The clone library approach showed dominance of Bacillariophyta-like sequences, in addition to Cryptophyta-, Haptophyta-, Pelagophyta-, and Eustigmatophyta-like sequences which were detected for the first time highlighting their importance in mangrove ecosystem. This study clearly shows that a combination of microscopy and molecular tools can improve understanding of phytoplankton assemblages in mangrove environments.

Linking phytoplankton community composition to seasonal changes in f-ratio

Seasonal changes in nitrogen assimilation have been studied in the western English Channel by sampling at approximately weekly intervals for 12 months. Nitrate concentrations showed strong seasonal variations. Available nitrogen in the winter was dominated by nitrate but this was close to limit of detection from May to September, after the spring phytoplankton bloom. The (15)N uptake experiments showed that nitrate was the nitrogen source for the spring phytoplankton bloom but regenerated nitrogen supported phytoplankton productivity throughout the summer. The average annual f-ratio was 0.35, which demonstrated the importance of ammonia regeneration in this dynamic temperate region. Nitrogen uptake rate measurements were related to the phytoplankton responsible by assessing the relative abundance of nitrate reductase (NR) genes and the expression of NR among eukaryotic phytoplankton. Strong signals were detected from NR sequences that are not associated with known phylotypes or cultures. NR sequences from the diatom Phaeodactylum tricornutum were highly represented in gene abundance and expression, and were significantly correlated with f-ratio. The results demonstrate that analysis of functional genes provides additional information, and may be able to give better indications of which phytoplankton species are responsible for the observed seasonal changes in f-ratio than microscopic phytoplankton identification.

Application of rbcL based molecular diversity analysis to algae in wastewater treatment plants

The molecular diversity of algae in the final clarifier or denitrification filter outfall from three wastewater treatment plants (WWTPs) with activated sludge based treatment was analyzed using the rbcL gene as a phylogenetic marker. The rbcL gene encodes the large subunit of the CO(2) fixing enzyme RuBisCO. Among algae identified at the WWTPs were diatoms, green algae, cyanobacteria, Eustigmatophyceae, and unknown heterokonts. A high level of diversity was observed within WWTPs with 19-24 unique rbcL sequences detected at each plant. Algae composition also varied between treatment plants. Our results show that the rbcL gene can be used as a phylogenetic marker for algae diversity analysis in a wastewater treatment context.