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.

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.

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.

Host genotype structures the microbiome of a globally dispersed marine phytoplankton.

Phytoplankton support complex bacterial microbiomes that rely on phytoplankton-derived extracellular compounds and perform functions necessary for algal growth. Recent work has revealed sophisticated interactions and exchanges of molecules between specific phytoplankton-bacteria pairs, but the role of host genotype in regulating those interactions is unknown. Here, we show how phytoplankton microbiomes are shaped by intraspecific genetic variation in the host using global environmental isolates of the model phytoplankton host Thalassiosira rotula and a laboratory common garden experiment. A set of 81 environmental T. rotula genotypes from three ocean basins and eight genetically distinct populations did not reveal a core microbiome. While no single bacterial phylotype was shared across all genotypes, we found strong genotypic influence of T. rotula, with microbiomes associating more strongly with host genetic population than with environmental factors. The microbiome association with host genetic population persisted across different ocean basins, suggesting that microbiomes may be associated with host populations for decades. To isolate the impact of host genotype on microbiomes, a common garden experiment using eight genotypes from three distinct host populations again found that host genotype influenced microbial community composition, suggesting that a process we describe as genotypic filtering, analogous to environmental filtering, shapes phytoplankton microbiomes. In both the environmental and laboratory studies, microbiome variation between genotypes suggests that other factors influenced microbiome composition but did not swamp the dominant signal of host genetic background. The long-term association of microbiomes with specific host genotypes reveals a possible mechanism explaining the evolution and maintenance of complex phytoplankton-bacteria chemical exchanges.

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.

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.

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.

Evidence for environmental and ecological selection in a microbe with no geographic limits to gene flow.

The ability for organisms to disperse throughout their environment is thought to strongly influence population structure and thus evolution of diversity within species. A decades-long debate surrounds processes that generate and support high microbial diversity, particularly in the ocean. The debate concerns whether diversification occurs primarily through geographic partitioning (where distance limits gene flow) or through environmental selection, and remains unresolved due to lack of empirical data. Here we show that gene flow in a diatom, an ecologically important eukaryotic microbe, is not limited by global-scale geographic distance. Instead, environmental and ecological selection likely play a more significant role than dispersal in generating and maintaining diversity. We detected significantly diverged populations (FST > 0.130) and discovered temporal genetic variability at a single site that was on par with spatial genetic variability observed over distances of 15,000 km. Relatedness among populations was decoupled from geographic distance across the global ocean and instead, correlated significantly with water temperature and whole-community chlorophyll a Correlations with temperature point to the importance of environmental selection in structuring populations. Correlations with whole-community chlorophyll a, a proxy for autotrophic biomass, suggest that ecological selection via interactions with other plankton may generate and maintain population genetic structure in marine microbes despite global-scale dispersal. Here, we provide empirical evidence for global gene flow in a marine eukaryotic microbe, suggesting that everything holds the potential to be everywhere, with environmental and ecological selection rather than geography or dispersal dictating the structure and evolution of diversity over space and time.

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.

Metatranscriptome analyses indicate resource partitioning between diatoms in the field.

Diverse communities of marine phytoplankton carry out half of global primary production. The vast diversity of the phytoplankton has long perplexed ecologists because these organisms coexist in an isotropic environment while competing for the same basic resources (e.g., inorganic nutrients). Differential niche partitioning of resources is one hypothesis to explain this "paradox of the plankton," but it is difficult to quantify and track variation in phytoplankton metabolism in situ. Here, we use quantitative metatranscriptome analyses to examine pathways of nitrogen (N) and phosphorus (P) metabolism in diatoms that cooccur regularly in an estuary on the east coast of the United States (Narragansett Bay). Expression of known N and P metabolic pathways varied between diatoms, indicating apparent differences in resource utilization capacity that may prevent direct competition. Nutrient amendment incubations skewed N/P ratios, elucidating nutrient-responsive patterns of expression and facilitating a quantitative comparison between diatoms. The resource-responsive (RR) gene sets deviated in composition from the metabolic profile of the organism, being enriched in genes associated with N and P metabolism. Expression of the RR gene set varied over time and differed significantly between diatoms, resulting in opposite transcriptional responses to the same environment. Apparent differences in metabolic capacity and the expression of that capacity in the environment suggest that diatom-specific resource partitioning was occurring in Narragansett Bay. This high-resolution approach highlights the molecular underpinnings of diatom resource utilization and how cooccurring diatoms adjust their cellular physiology to partition their niche space.

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.

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.

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

Diatoms generate nearly half of marine primary production and are comprised of a diverse array of species that are often morphologically cryptic or difficult to identify using light microscopy. Here, species composition and realized thermal niches of species in the diatom genus Thalassiosira were examined at the site of the Narragansett Bay (NBay) Long-Term Plankton Time Series using a combination of light microscopy (LM), high-throughput sequencing (HTS) of the 18S rDNA V4 region and historical records. Thalassiosira species were identified over 6 years using a combination of LM and DNA sequences. Sixteen Thalassiosira taxa were identified using HTS: nine were newly identified in NBay. Several newly identified species have small cell diameters and are difficult to identify using LM. However, they appeared frequently and thus may play a significant ecological role in NBay, particularly since their realized niches suggest they are eurythermal and able to tolerate the >25 °C temperature range of NBay. Four distinct species assemblages that grouped by season were best explained by surface water temperature. When compared to historical records, we found that the cold-water species Thalassiosira nordenskioeldii has decreased in persistence over time, suggesting that increasing surface water temperature has influenced the ecology of phytoplankton in NBay.

Metapopulation Structure in the Planktonic Diatom Ditylum brightwellii (Bacillariophyceae).

Approximately 200,000 diatom species are thought to exist and yet the underlying processes of speciation in diatoms are unknown. Because genetic subdivision within species can reveal potential speciation mechanisms, we examined genetic differentiation and patterns of gene flow among four populations of the diatom Ditylum brightwellii. Single-cell isolates were examined at two microsatellite markers and two rDNA loci (18S and internal transcribed spacer region I (ITSI)). Among isolates, rDNA sequences varied by 0.08+/-0.04% (18S) and 0.7+/-0.3% (ITSI) and there were no compensatory base pair changes in the predicted ITSI secondary structure, all suggesting that a single species was represented. Two numerically dominant ITSI sequence types were detected and their distribution among isolates from genetically distinct populations was significantly different. Two populations shared ITSI sequence type 1 and two shared ITSI sequence type 2, indicating differences in relatedness among populations. The signature of unequal gene flow among populations suggested that D. brightwellii exhibited a metapopulation structure: the species was subdivided into populations of populations. The identification of metapopulations suggests a possible mechanism of speciation through reduced levels of gene flow, providing newly evolved taxa with a large repository of genetic and physiological diversity and perhaps significant adaptive potential.

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.

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.