Phytoplankton thermal trait parameterization alters community structure and biogeochemical processes in a modeled ocean

Phytoplankton exhibit diverse physiological responses to temperature which influence their fitness in the environment and consequently alter their community structure. Here, we explored the sensitivity of phytoplankton community structure to thermal response parameterization in a modelled marine phytoplankton community. Using published empirical data, we evaluated the maximum thermal growth rates (μmax ) and temperature coefficients (Q10 ; the rate at which growth scales with temperature) of six key Phytoplankton Functional Types (PFTs): coccolithophores, cyanobacteria, diatoms, diazotrophs, dinoflagellates, and green algae. Following three well-documented methods, PFTs were either assumed to have (1) the same μmax and the same Q10 (as in to Eppley, 1972), (2) a unique μmax but the same Q10 (similar to Kremer et al., 2017), or (3) a unique μmax and a unique Q10 (following Anderson et al., 2021). These trait values were then implemented within the Massachusetts Institute of Technology biogeochemistry and ecosystem model (called Darwin) for each PFT under a control and climate change scenario. Our results suggest that applying a μmax and Q10 universally across PFTs (as in Eppley, 1972) leads to unrealistic phytoplankton communities, which lack diatoms globally. Additionally, we find that accounting for differences in the Q10 between PFTs can significantly impact each PFT's competitive ability, especially at high latitudes, leading to altered modeled phytoplankton community structures in our control and climate change simulations. This then impacts estimates of biogeochemical processes, with, for example, estimates of export production varying by ~10% in the Southern Ocean depending on the parameterization. Our results indicate that the diversity of thermal response traits in phytoplankton not only shape community composition in the historical and future, warmer ocean, but that these traits have significant feedbacks on global biogeochemical cycles.

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.

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.

Thermal trait variation may buffer Southern Ocean phytoplankton from anthropogenic warming

Despite the potential of standing genetic variation to rescue communities and shape future adaptation to climate change, high levels of uncertainty are associated with intraspecific trait variation in marine phytoplankton. Recent model intercomparisons have pointed to an urgent need to reduce uncertainty in the projected responses of marine ecosystems to climate change, including Southern Ocean (SO) surface waters, which are among the most rapidly warming habitats on Earth. Because SO phytoplankton growth responses to warming sea surface temperature (SST) are poorly constrained, we developed a high-throughput growth assay to simultaneously examine inter- and intra-specific thermal trait variation in a group of 43 taxonomically diverse and biogeochemically important SO phytoplankton called diatoms. We found significant differential growth performance among species across thermal traits, including optimum and maximum tolerated growth temperatures. Within species, coefficients of variation ranged from 3% to 48% among strains for those same key thermal traits. Using SO SST projections for 2100, we predicted biogeographic ranges that differed by up to 97% between the least and most tolerant strains for each species, illustrating the role that strain-specific differences in temperature response can play in shaping predictions of future phytoplankton biogeography. Our findings revealed the presence and scale of thermal trait variation in SO phytoplankton and suggest these communities may already harbour the thermal trait diversity required to withstand projected 21st-century SST change in the SO even under severe climate forcing scenarios.

Identification and Expression Analyses of the Nitrate Transporter Gene (NRT2) Family Among Skeletonema species (Bacillariophyceae)

High-affinity nitrate transporters are considered to be the major transporter system for nitrate uptake in diatoms. In the diatom genus Skeletonema, three forms of genes encoding high-affinity nitrate transporters (NRT2) were newly identified from transcriptomes generated as part of the marine microbial eukaryote transcriptome sequencing project. To examine the expression of each form of NRT2 under different nitrogen environments, laboratory experiments were conducted under nitrate-sufficient, ammonium-sufficient, and nitrate-limited conditions using three ecologically important Skeletonema species: S. dohrnii, S. menzelii, and S. marinoi. Primers were developed for each NRT2 form and species and Q-RT-PCR was performed. For each NRT2 form, the three Skeletonema species had similar transcriptional patterns. The transcript levels of NRT2:1 were significantly elevated under nitrogen-limited conditions, but strongly repressed in the presence of ammonium. The transcript levels of NRT2:2 were also repressed by ammonium, but increased 5- to 10-fold under nitrate-sufficient and nitrogen-limited conditions. Finally, the transcript levels of NRT2:3 did not vary significantly under various nitrogen conditions, and behaved more like a constitutively expressed gene. Based on the observed transcript variation among NRT2 forms, we propose a revised model describing nitrate uptake kinetics regulated by multiple forms of nitrate transporter genes in response to various nitrogen conditions in Skeletonema. The differential NRT2 transcriptional responses among species suggest that species-specific adaptive strategies exist within this genus to cope with environmental changes.

Navigating in a sea of genes

Randomly generated microbial genomes produce realistic nutrient gradients in a model ocean

pCO2 effects on species composition and growth of an estuarine phytoplankton community

The effects of ongoing changes in ocean carbonate chemistry on plankton ecology have important implications for food webs and biogeochemical cycling. However, conflicting results have emerged regarding species-specific responses to pCO2 enrichment and thus community responses have been difficult to predict. To assess community level effects (e.g., production) of altered carbonate chemistry, studies are needed that capitalize on the benefits of controlled experiments but also retain features of intact ecosystems that may exacerbate or ameliorate the effects observed in single-species or single cohort experiments. We performed incubations of natural plankton communities from Narragansett Bay, RI, USA in winter at ambient bay temperatures (5-13 °C), light and nutrient concentrations under three levels of controlled and constant CO2 concentrations, simulating past, present and future conditions at mean pCO2 levels of 224, 361, and 724 μatm respectively. Samples for carbonate analysis, chlorophyll a, plankton size-abundance, and plankton species composition were collected daily and phytoplankton growth rates in three different size fractions (<5, 5-20, and >20 μm) were measured at the end of the 7-day incubation period. Community composition changed during the incubation period with major increases in relative diatom abundance, which were similar across pCO2 treatments. At the end of the experiment, 24-hr growth responses to pCO2 levels varied as a function of cell size. The smallest size fraction (<5 μm) grew faster at the elevated pCO2 level. In contrast, the 5-20 μm size fraction grew fastest in the Present treatment and there were no significant differences in growth rate among treatments in the > 20 μm size fraction. Cell size distribution shifted toward smaller cells in both the Past and Future treatments but remained unchanged in the Present treatment. Similarity in Past and Future treatments for cell size distribution and growth rate (5-20 μm size fraction) illustrate non-monotonic effects of increasing pCO2 on ecological indicators and may be related to opposing physiological effects of high CO2 and low pH both within and among species. Interaction of these effects with other factors (e.g., nutrients, light, temperature, grazing, initial species composition) may explain variability among published studies. The absence of clear treatment-specific effects at the community level suggest that extrapolation of species-specific responses or experiments with only present day and future pCO2 treatments levels would produce misleading predictions of ocean acidification impacts on plankton production.

Persistent Intra-Specific Variation in Genetic and Behavioral Traits in the Raphidophyte, Heterosigma akashiwo

Motility is a key trait that phytoplankton utilize to navigate the heterogeneous marine environment. Quantifying both intra- and inter-specific variability in trait distributions is key to utilizing traits to distinguish groups of organisms and assess their ecological function. Because examinations of intra-specific variability are rare, here we measured three-dimensional movement behaviors and distribution patterns of seven genetically distinct strains of the ichthyotoxic raphidophyte, Heterosigma akashiwo. Strains were collected from different ocean basins but geographic distance between isolates was a poor predictor of genetic relatedness among strains. Observed behaviors were significantly different among all strains examined, with swimming speed and turning rate ranging from 33–115 μm s^-1 and 41–110° s^-1, respectively. Movement behaviors were consistent over at least 12 h, and in one case identical when measured several years apart. Movement behaviors were not associated with a specific cell size, carbon content, genetic relatedness, or geographic distance. These strain-specific behaviors resulted in algal populations that had distinct vertical distributions in the experimental tank. This study demonstrates that the traits of genetic identity and motility can provide resolution to distinguish strains of species, where variations in size or biomass are insufficient characteristics.

Experimental evolution gone wild

Because of their large population sizes and rapid cell division rates, marine microbes have, or can generate, ample variation to fuel evolution over a few weeks or months, and subsequently have the potential to evolve in response to global change. Here we measure evolution in the marine diatom Skeletonema marinoi evolved in a natural plankton community in CO2-enriched mesocosms deployed in situ. Mesocosm enclosures are typically used to study how the species composition and biogeochemistry of marine communities respond to environmental shifts, but have not been used for experimental evolution to date. Using this approach, we detect a large evolutionary response to CO2 enrichment in a focal marine diatom, where population growth rate increased by 1.3-fold in high CO2-evolved lineages. This study opens an exciting new possibility of carrying out in situ evolution experiments to understand how marine microbial communities evolve in response to environmental change.

Evolutionary potential of marine phytoplankton under ocean acidification

Marine phytoplankton have many obvious characters, such as rapid cell division rates and large population sizes, that give them the capacity to evolve in response to global change on timescales of weeks, months or decades. However, few studies directly investigate if this adaptive potential is likely to be realized. Because of this, evidence of to whether and how marine phytoplankton may evolve in response to global change is sparse. Here, we review studies that help predict evolutionary responses to global change in marine phytoplankton. We find limited support from experimental evolution that some taxa of marine phytoplankton may adapt to ocean acidification, and strong indications from studies of variation and structure in natural populations that selection on standing genetic variation is likely. Furthermore, we highlight the large body of literature on plastic responses to ocean acidification available, and evolutionary theory that may be used to link plastic and evolutionary responses. Because of the taxonomic breadth spanned by marine phytoplankton, and the diversity of roles they fill in ocean ecosystems and biogeochemical cycles, we stress the necessity of treating taxa or functional groups individually.

Marine phytoplankton temperature versus growth responses from polar to tropical waters--outcome of a scientific community-wide study

"It takes a village to finish (marine) science these days" Paraphrased from Curtis Huttenhower (the Human Microbiome project) The rapidity and complexity of climate change and its potential effects on ocean biota are challenging how ocean scientists conduct research. One way in which we can begin to better tackle these challenges is to conduct community-wide scientific studies. This study provides physiological datasets fundamental to understanding functional responses of phytoplankton growth rates to temperature. While physiological experiments are not new, our experiments were conducted in many laboratories using agreed upon protocols and 25 strains of eukaryotic and prokaryotic phytoplankton isolated across a wide range of marine environments from polar to tropical, and from nearshore waters to the open ocean. This community-wide approach provides both comprehensive and internally consistent datasets produced over considerably shorter time scales than conventional individual and often uncoordinated lab efforts. Such datasets can be used to parameterise global ocean model projections of environmental change and to provide initial insights into the magnitude of regional biogeographic change in ocean biota in the coming decades. Here, we compare our datasets with a compilation of literature data on phytoplankton growth responses to temperature. A comparison with prior published data suggests that the optimal temperatures of individual species and, to a lesser degree, thermal niches were similar across studies. However, a comparison of the maximum growth rate across studies revealed significant departures between this and previously collected datasets, which may be due to differences in the cultured isolates, temporal changes in the clonal isolates in cultures, and/or differences in culture conditions. Such methodological differences mean that using particular trait measurements from the prior literature might introduce unknown errors and bias into modelling projections. Using our community-wide approach we can reduce such protocol-driven variability in culture studies, and can begin to address more complex issues such as the effect of multiple environmental drivers on ocean biota.

Identifying reference genes with stable expression from high throughput sequence data

Genes that are constitutively expressed across multiple environmental stimuli are crucial to quantifying differentially expressed genes, particularly when employing quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) assays. However, the identification of these potential reference genes in non-model organisms is challenging and is often guided by expression patterns in distantly related organisms. Here, transcriptome datasets from the diatom Thalassiosira pseudonana grown under replete, phosphorus-limited, iron-limited, and phosphorus and iron co-limited nutrient regimes were analyzed through literature-based searches for homologous reference genes, k-means clustering, and analysis of sequence counts (ASC) to identify putative reference genes. A total of 9759 genes were identified and screened for stable expression. Literature-based searches surveyed 18 generally accepted reference genes, revealing 101 homologs in T. pseudonana with variable expression and a wide range of mean tags per million. k-means analysis parsed the whole transcriptome into 15 clusters. The two most stable clusters contained 709 genes, but still had distinct patterns in expression. ASC analyses identified 179 genes that were stably expressed (posterior probability < 0.1 for 1.25 fold change). Genes known to have a stable expression pattern across the test treatments, like actin, were identified in this pool of 179 candidate genes. ASC can be employed on data without biological replicates and was more robust than the k-means approach in isolating genes with stable expression. The intersection of the genes identified through ASC with commonly used reference genes from the literature suggests that actin and ubiquitin ligase may be useful reference genes for T. pseudonana and potentially other diatoms. With the wealth of transcriptome sequence data becoming available, ASC can be easily applied to transcriptome datasets from other phytoplankton to identify reference genes.

Molecular subdivision of the marine diatom Thalassiosira rotula in relation to geographic distribution, genome size, and physiology

Background

Marine phytoplankton drift passively with currents, have high dispersal potentials and can be comprised of morphologically cryptic species. To examine molecular subdivision in the marine diatom Thalassiosira rotula, variations in rDNA sequence, genome size, and growth rate were examined among isolates collected from the Atlantic and Pacific Ocean basins. Analyses of rDNA included T. gravida because morphological studies have argued that T. rotula and T. gravida are conspecific.

Results

Culture collection isolates of T. gravida and T. rotula diverged by 7.0 ± 0.3% at the ITS1 and by 0.8 ± 0.03% at the 28S. Within T. rotula, field and culture collection isolates were subdivided into three lineages that diverged by 0.6 ± 0.3% at the ITS1 and 0% at the 28S. The predicted ITS1 secondary structure revealed no compensatory base pair changes among lineages. Differences in genome size were observed among isolates, but were not correlated with ITS1 lineages. Maximum acclimated growth rates of isolates revealed genotype by environment effects, but these were also not correlated with ITS1 lineages. In contrast, intra-individual variation in the multi-copy ITS1 revealed no evidence of recombination amongst lineages, and molecular clock estimates indicated that lineages diverged 0.68 Mya. The three lineages exhibited different geographic distributions and, with one exception, each field sample was dominated by a single lineage.

Conclusions

The degree of inter- and intra-specific divergence between T. gravida and T. rotula suggests they should continue to be treated as separate species. The phylogenetic distinction of the three closely-related T. rotula lineages was unclear. On the one hand, the lineages showed no physiological differences, no consistent genome size differences and no significant changes in the ITS1 secondary structure, suggesting there are no barriers to interbreeding among lineages. In contrast, analysis of intra-individual variation in the multicopy ITS1 as well as molecular clock estimates of divergence suggest these lineages have not interbred for significant periods of time. Given the current data, these lineages should be considered a single species. Furthermore, these T. rotula lineages may be ecologically relevant, given their differential abundance over large spatial scales.

The transcriptome and proteome of the diatom Thalassiosira pseudonana reveal a diverse phosphorus stress response

Phosphorus (P) is a critical driver of phytoplankton growth and ecosystem function in the ocean. Diatoms are an abundant class of marine phytoplankton that are responsible for significant amounts of primary production. With the control they exert on the oceanic carbon cycle, there have been a number of studies focused on how diatoms respond to limiting macro and micronutrients such as iron and nitrogen. However, diatom physiological responses to P deficiency are poorly understood. Here, we couple deep sequencing of transcript tags and quantitative proteomics to analyze the diatom Thalassiosira pseudonana grown under P-replete and P-deficient conditions. A total of 318 transcripts were differentially regulated with a false discovery rate of <0.05, and a total of 136 proteins were differentially abundant (p<0.05). Significant changes in the abundance of transcripts and proteins were observed and coordinated for multiple biochemical pathways, including glycolysis and translation. Patterns in transcript and protein abundance were also linked to physiological changes in cellular P distributions, and enzyme activities. These data demonstrate that diatom P deficiency results in changes in cellular P allocation through polyphosphate production, increased P transport, a switch to utilization of dissolved organic P through increased production of metalloenzymes, and a remodeling of the cell surface through production of sulfolipids. Together, these findings reveal that T. pseudonana has evolved a sophisticated response to P deficiency involving multiple biochemical strategies that are likely critical to its ability to respond to variations in environmental P availability.

Learning to read the oceans genomics of marine phytoplankton

The phytoplankton are key members of marine ecosystems, generating about half of global primary productivity, supporting valuable fisheries and regulating global biogeochemical cycles. Marine phytoplankton are phylogenetically diverse and are comprised of both prokaryotic and eukaryotic species. In the last decade, new insights have been gained into the ecology and evolution of these important organisms through whole genome sequencing projects and more recently, through both transcriptomics and targeted metagenomics approaches. Sequenced genomes of cyanobacteria are generally small, ranging in size from 1.8 to 9 million base pairs (Mbp). Eukaryotic genomes, in general, have a much larger size range and those that have been sequenced range from 12 to 57 Mbp. Whole genome sequencing projects have revealed key features of the evolutionary history of marine phytoplankton, their varied responses to environmental stress, their ability to scavenge and store nutrients and their unique ability to form elaborate cellular coverings. We have begun to learn how to read the 'language' of marine phytoplankton, as written in their DNA. Here, we review the ecological and evolutionary insights gained from whole genome sequencing projects, illustrate how these genomes are yielding information on marine natural products and informing nanotechnology as well as make suggestions for future directions in the field of marine phytoplankton genomics.

Empirical bayes analysis of sequencing-based transcriptional profiling without replicates

Background

Recent technological advancements have made high throughput sequencing an increasingly popular approach for transcriptome analysis. Advantages of sequencing-based transcriptional profiling over microarrays have been reported, including lower technical variability. However, advances in technology do not remove biological variation between replicates and this variation is often neglected in many analyses.

Results

We propose an empirical Bayes method, titled Analysis of Sequence Counts (ASC), to detect differential expression based on sequencing technology. ASC borrows information across sequences to establish prior distribution of sample variation, so that biological variation can be accounted for even when replicates are not available. Compared to current approaches that simply tests for equality of proportions in two samples, ASC is less biased towards highly expressed sequences and can identify more genes with a greater log fold change at lower overall abundance.

Conclusions

ASC unifies the biological and statistical significance of differential expression by estimating the posterior mean of log fold change and estimating false discovery rates based on the posterior mean. The implementation in R is available at http://www.stat.brown.edu/Zwu/research.aspx.

Maintenance of clonal diversity during a spring bloom of the centric diatom Ditylum brightwellii

Maintenance of genetic diversity in eukaryotic microbes reflects a synergism between reproductive mode (asexual vs. sexual) and environmental conditions. We determined clonal diversity in field samples of the planktonic marine diatom, Ditylum brightwellii, during a bloom, when cell number increased by seven-fold because of rapid asexual division. The genotypes at three microsatellite loci were determined for 607 individual cell lines isolated during the 11 days of sampling. Genetic diversity remained high during the bloom and 87% of the cells sampled each day were genetically distinct. Sixty-nine clonal lineages were sampled two or more times during the bloom, and two clones were sampled seven times. Based on the frequency of resampled clonal lineages, capture-recapture statistics were used to determine that at least 2400 genetically distinct clonal lineages comprised the bloom population. No significant differences in microsatellite allele frequencies were observed among daily samples indicating that the bloom was comprised of a single population. No sexual stages were observed, although linkage equilibrium at two loci, high levels of allelic and genotypic diversity, and heterozygote deficiencies were all indicative of past sexual reproduction events. At the height of the bloom, a windstorm diluted cell numbers by 51% and coincided with a change in the frequency distribution of some resampled lineages. The extensive clonal diversity generated through past sexual reproduction events coupled with frequent environmental changes appear to prevent individual clonal lineages from becoming numerically dominant, maintaining genetic diversity and the adaptive potential of the population.

The Genome of the Diatom Thalassiosira Pseudonana: Ecology, Evolution, and Metabolism

Diatoms are unicellular algae with plastids acquired by secondary endosymbiosis. They are responsible for approximately 20% of global carbon fixation. We report the 34 million-base pair draft nuclear genome of the marine diatom Thalassiosira pseudonana and its 129 thousand-base pair plastid and 44 thousand-base pair mitochondrial genomes. Sequence and optical restriction mapping revealed 24 diploid nuclear chromosomes. We identified novel genes for silicic acid transport and formation of silica-based cell walls, high-affinity iron uptake, biosynthetic enzymes for several types of polyunsaturated fatty acids, use of a range of nitrogenous compounds, and a complete urea cycle, all attributes that allow diatoms to prosper in aquatic environments.