The role of intraspecific variation in the ecological and evolutionary success of diatoms in changing environments

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

Allele-dependent expression and functionality of lipid enzyme phospholipid:diacylglycerol acyltransferase affect diatom carbon storage and growth

In the acyl-CoA-independent pathway of triacylglycerol (TAG) synthesis unique to plants, fungi, and algae, TAG formation is catalyzed by the enzyme phospholipid:diacylglycerol acyltransferase (PDAT). The unique PDAT gene of the model diatom Phaeodactylum tricornutum strain CCMP2561 boasts 47 single nucleotide variants within protein coding regions of the alleles. To deepen our understanding of TAG synthesis, we observed the allele-specific expression of PDAT by the analysis of 87 published RNA-sequencing (RNA-seq) data and experimental validation. The transcription of one of the two PDAT alleles, Allele 2, could be specifically induced by decreasing nitrogen concentrations. Overexpression of Allele 2 in P. tricornutum substantially enhanced the accumulation of TAG by 44% to 74% under nutrient stress; however, overexpression of Allele 1 resulted in little increase of TAG accumulation. Interestingly, a more serious growth inhibition was observed in the PDAT Allele 1 overexpression strains compared with Allele 2 counterparts. Heterologous expression in yeast (Saccharomyces cerevisiae) showed that enzymes encoded by PDAT Allele 2 but not Allele 1 had TAG biosynthetic activity, and 7 N-terminal and 3 C-terminal amino acid variants between the 2 allele-encoded proteins substantially affected enzymatic activity. P. tricornutum PDAT, localized in the innermost chloroplast membrane, used monogalactosyldiacylglycerol and phosphatidylcholine as acyl donors as demonstrated by the increase of the 2 lipids in PDAT knockout lines, which indicated a common origin in evolution with green algal PDATs. Our study reveals unequal roles among allele-encoded PDATs in mediating carbon storage and growth in response to nitrogen stress and suggests an unsuspected strategy toward lipid and biomass improvement for biotechnological purposes.

Temperature-Related Short-Term Succession Events of Bacterial Phylotypes in Potter Cove, Antarctica

In recent years, our understanding of the roles of bacterial communities in the Antarctic Ocean has substantially improved. It became evident that Antarctic marine bacteria are metabolically versatile, and even closely related strains may differ in their functionality and, therefore, affect the ecosystem differently. Nevertheless, most studies have been focused on entire bacterial communities, with little attention given to individual taxonomic groups. Antarctic waters are strongly influenced by climate change; thus, it is crucial to understand how changes in environmental conditions, such as changes in water temperature and salinity fluctuations, affect bacterial species in this important area. In this study, we show that an increase in water temperature of 1 °C was enough to alter bacterial communities on a short-term temporal scale. We further show the high intraspecific diversity of Antarctic bacteria and, subsequently, rapid intra-species succession events most likely driven by various temperature-adapted phylotypes. Our results reveal pronounced changes in microbial communities in the Antarctic Ocean driven by a single strong temperature anomaly. This suggests that long-term warming may have profound effects on bacterial community composition and presumably functionality in light of continuous and future climate change.

Diatom endemism and taxonomic turnover: Assessment in high-altitude alpine lakes covering a large geographical range

Diatoms are widely used as ecological indicators and show various degrees of endemism. Many studies that support the idea of endemic species integrate several climate zones, a variety of ecosystem types, and often focus on a global scale. Here, we investigated whether endemism could be detected when considering a homogeneous type of ecosystem in a single climate zone. We sampled stone biofilms at 40-50 cm depth in high-altitude lakes in the Alpine climate zone. A total of 149 samples were obtained from the French and Georgian mountains, two areas separated by ∼3000 km. Using Amplicon Sequence Variants derived from DNA metabarcoding, we assessed taxonomic turnover and Zeta-diversity (a measure of endemism). We ran haplotype networks and phylogenetic tests to measure geographical signal in the phylogenies of dominant taxa. The French and Georgian communities shared 51 % of species. Species that were not shared across both regions were mostly rare, and often not characteristic of lakes but of neighboring habitats instead. In contrast, at the sub-species level, 87 % of the genotypes showed restricted distributions. Whereas endemism was the rule at sub-species level, most species were shared across both French and Georgian lakes, suggesting that geographic barriers strongly limited dispersal at the sub-species level but not species level. Dominant species hosted higher levels of sub-specific diversity than rare species. In contrast to global-scale studies, we did not find any significant geographical structuring in the phylogeny of the investigated species. This could indicate ongoing dispersal at a frequency fast enough to prevent allopatric divergence, yet slow enough to prevent sharing most haplotypes between France and Georgia. These results have implications for biomonitoring: depending on the taxonomic level chosen, robust generic tools (species level) or tools dedicated to a region able to discriminate fine pressures differences (sub-species level) may be developed.

How to improve scaling from traits to ecosystem processes

Scaling approaches in ecology assume that traits are the main attributes by which organisms influence ecosystem functioning. However, several recent empirical papers have found only weak links between traits and ecosystem functioning, questioning the usefulness of trait-based ecology (TBE). We argue that these studies often suffer from one or more widespread misconceptions. Specifically, these studies often (i) conflict with the conceptual foundations of TBE, (ii) lack theory- or hypothesis-driven selection and use of traits, (iii) tend to ignore intraspecific variation, and (iv) use experimental or study designs that are not well suited to make strong tests of TBE assumptions. Addressing these aspects could significantly improve our ability to scale from traits to ecosystem functioning.

Winners and Losers of Atlantification: The Degree of Ocean Warming Affects the Structure of Arctic Microbial Communities

Arctic microbial communities (i.e., protists and bacteria) are increasingly subjected to an intrusion of new species via Atlantification and an uncertain degree of ocean warming. As species differ in adaptive traits, these oceanic conditions may lead to compositional changes with functional implications for the ecosystem. In June 2021, we incubated water from the western Fram Strait at three temperatures (2 °C, 6 °C, and 9 °C), mimicking the current and potential future properties of the Arctic Ocean. Our results show that increasing the temperature to 6 °C only minorly affects the community, while an increase to 9 °C significantly lowers the diversity and shifts the composition. A higher relative abundance of large hetero- and mixotrophic protists was observed at 2 °C and 6 °C compared to a higher abundance of intermediate-sized temperate diatoms at 9 °C. The compositional differences at 9 °C led to a higher chlorophyll a:POC ratio, but the C:N ratio remained similar. Our results contradict the common assumption that smaller organisms and heterotrophs are favored under warming and strongly indicate a thermal limit between 6 °C and 9 °C for many Arctic species. Consequently, the magnitude of temperature increase is a crucial factor for microbial community reorganization and the ensuing ecological consequences in the future Arctic Ocean.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Geographical distribution, molecular and toxin diversity of the dinoflagellate species Gambierdiscus honu in the Pacific region

An increase in cases of ciguatera poisoning (CP) and expansion of the causative species in the South Pacific region highlight the need for baseline data on toxic microalgal species to help identify new areas of risk and manage known hot spots. Gambierdiscus honu is a toxin producing and potential CP causing dinoflagellate species, first described in 2017. Currently no high-resolution geographical distribution, intraspecific genetic variation or toxin production diversity data is available for G. honu. This research aimed to further characterize G. honu by investigating its distribution using species-specific real-time polymerase chain reaction assays at 25 sites in an area spanning ∼8000 km of the Coral Sea/Pacific Ocean, and assessing intraspecific genetic variation, toxicity and toxin production of isolated strains. Assessment of genetic variation of the partial rRNA operon of isolates demonstrated no significant intraspecific population structure, in addition to a lack of adherence to isolation by distance (IBD) model of evolution. The detected distribution of G. honu in the Pacific region was within the expected tropical to temperate latitudinal ranges of 10° to -30° and extended from Australia to French Polynesia. In the lipophilic fractions, the neuroblastoma cell-based assay (CBA-N2a) showed no ciguatoxin (CTX)-like activity for nine of the 10 isolates, and an atypical pattern for CAWD233 isolate which showed cytotoxic activity in OV- and OV+ conditions. In the same way, liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis confirmed no Pacific-CTXs (CTX-3B, CTX-3C, CTX-4A, CTX-4B) were produced by the ten strains. The CBA-N2a assessment of the hydrophilic fractions showed moderate to high cytotoxicity in both OV- and OV+ condition for all the strains showing a cytotoxic profile similar to that of gambierone. Indeed, this study is the first to show the cytotoxic activity of gambierone on mouse neuroblastoma cells while no cytotoxicity was observed when 44-MG was analysed at the same concentrations using the CBA-N2a. Analysis of the hydrophilic via LC-MS/MS confirmed production of gambierone in all isolates, ranging from 2.1 to 38.1 pg/cell, with 44-methylgambierone (44-MG) also produced by eight of the isolates, ranging from 0.3 to 42.9 pg/cell. No maitotoxin-1 was detected in any of the isolates. Classification of the G. honu strains according to the quantities of gambierone produced aligned with the classification of their cytotoxicity using the CBA-N2a. Finally, no maitotoxin-1 (MTX) was detected in any of the isolates. This study shows G. honu is widely distributed within the Pacific region with no significant intraspecific population structure present. This aligns with the view of microalgal populations as global metapopulations, however more in-depth assessment with other genetic markers could detect further structure. Toxicity diversity across 10 isolates assessed did not display any geographical patterns.

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.

Diversity strengthens competing teams

How does the composition of a collection of individuals affect its outcome in competition with other collections of individuals? Assuming that individuals can be different, we develop a model to interpolate between individual-level interactions and collective-level consequences. Rooted in theoretical mathematics, the model is not constrained to any specific context. Potential applications include research, education, sports, politics, ecology, agriculture, algorithms and finance. Our first main contribution is a game theoretic model that interpolates between the internal composition of an ensemble of individuals and the repercussions for the ensemble as a whole in competition with others. The second main contribution is the rigorous identification of all equilibrium points and strategies. These equilibria suggest a mechanistic underpinning for biological and physical systems to tend towards increasing diversity due to the strength it imparts to the system in competition with others.

Strain-specific transcriptional responses overshadow salinity effects in a marine diatom sampled along the Baltic Sea salinity cline

The salinity gradient separating marine and freshwater environments represents a major ecological divide for microbiota, yet the mechanisms by which marine microbes have adapted to and ultimately diversified in freshwater environments are poorly understood. Here, we take advantage of a natural evolutionary experiment: the colonization of the brackish Baltic Sea by the ancestrally marine diatom Skeletonema marinoi. To understand how diatoms respond to low salinity, we characterized transcriptomic responses of acclimated S. marinoi grown in a common garden. Our experiment included eight strains from source populations spanning the Baltic Sea salinity cline. Gene expression analysis revealed that low salinities induced changes in the cellular metabolism of S. marinoi, including upregulation of photosynthesis and storage compound biosynthesis, increased nutrient demand, and a complex response to oxidative stress. However, the strain effect overshadowed the salinity effect, as strains differed significantly in their response, both regarding the strength and the strategy (direction of gene expression) of their response. The high degree of intraspecific variation in gene expression observed here highlights an important but often overlooked source of biological variation associated with how diatoms respond to environmental change.

Mycelial Growth-promoting Potential of Extracellular Metabolites of Paraburkholderia spp. Isolated from Rhizopogon roseolus Sporocarp

This study aimed to investigate the effect of potential metabolite(s) produced by Paraburkholderia spp. isolated from the Rhizopogon roseolus (shouro mushroom) sporocarp on the mycelial growth of R. roseolus. For this purpose, we selected two molecularly identified bacteria: P. fungorum GIB024 and P. caledonica KN1. Direct confrontation assay at three different distances, a pour plate method that sampled bacterial spent broth either with and without agitation at 25 °C, and an indirect confrontation assay was carried out in order to assess the R. roseolus growth-promoting ability of Paraburkholderia spp. These assessments were carried out in a 1:5 diluted Melin-Norkran-modified medium with glucose (hs-dMMN) and without glucose (ls-dMMN). GIB024 promoted the growth of R. roseolus in ls-dMMN in short distance, whereas KN1 inhibited the growth of the fungus in that condition. In hs-dMMN, both bacteria have neutral or slightly promotion effect toward R. roseolus. We determined from the spent broth analysis that Paraburkholderia spp. that grew axenically under static conditions had a more pronounced mycelial growth-promoting effect on R. roseolus than under agitation conditions. We also found that high concentration of spent broth resulted in a decrease in mycelial growth-promoting ability. Volatile metabolite(s) produced by both bacteria did not promote the mycelial growth of R. roseolus. In conclusion, Paraburkholderia spp. exhibited a species- and nutrient (sugar)-dependent ability to promote the mycelial growth of R. roseolus, and the bacterial soluble metabolite(s) play a crucial role in their growth-promoting ability.

Evolution of Phytoplankton as Estimated from Genetic Diversity

Phytoplankton are photosynthetic, single-celled organisms producing almost half of all oxygen on Earth and play a central role as prey for higher organisms, making them irreplaceable in the marine food web. As Global Change proceeds, imposing rapidly intensifying selection pressures, phytoplankton are forced to undergo evolution, local extinction, or redistribution, with potentially cascading effects throughout the marine ecosystem. Recent results from the field of population genetics display high levels of standing genetic diversity in natural phytoplankton populations, providing ample ‘evolutionary options’ and implying high adaptive potential to changing conditions. This potential for adaptive evolution is realized in several studies of experimental evolution, even though most of these studies investigate the evolution of only single strains. This, however, shows that phytoplankton not only evolve from standing genetic diversity, but also rely on de novo mutations. Recent global sampling campaigns show that the immense intraspecific diversity of phytoplankton in the marine ecosystem has been significantly underestimated, meaning we are only studying a minor portion of the relevant variability in the context of Global Change and evolution. An increased understanding of genomic diversity is primarily hampered by the low number of ecologically representative reference genomes of eukaryotic phytoplankton and the functional annotation of these. However, emerging technologies relying on metagenome and transcriptome data may offer a more realistic understanding of phytoplankton diversity.

Influence of land use and occupation on water quality and on the diatom community of reservoirs in northeast Brazil

Historically, the reduction in water quality has been associated with changes in land use, especially changes in urbanization and expansion of agricultural frontiers. The objective of this study was to understand the characteristics of the diatomological community and the use and occupation of the soil around three reservoirs (Arara, Rio do Canto and Araçagi) and natural lagoon (Lagoa do Paó) (Northeast Brazil). Land use was determined using high spatial resolution satellite images. Abiotic variables of water and total phosphorus, organic matter and diatoms in surface sediments were analyzed. Abiotic variables and assemblages of diatoms were analyzed using Principal Component Analysis (PCA). Four types of land use have been identified, associated with anthropogenic change and intensive use of the landscape. The PCA separated the reservoirs by trophic status based on the total phosphorus in the water. 24 taxa of diatoms distributed in 12 genera and 10 families were identified in the sediment. Our results indicate the occurrence of planktonic, epiphytic and benthic diatoms adapted to different trophic gradients

Sex-specific thermal tolerance limits in the ditch shrimp Palaemon varians: Eco-evolutionary implications under a warming ocean

As global temperatures continue to rise due to global change, marine heatwaves are also becoming more frequent and intense, impacting marine biodiversity patterns worldwide. Organisms inhabiting shallow water environments, such as the commercially relevant ditch shrimp Palaemon varians, are expected to be the most affected by rising temperatures. Thus, addressing species' thermal ecology and climate extinction-risk is crucial to foster climate-smart conservation strategies for shallow water ecosystems. Here, we estimated sex-specific upper thermal tolerance limits for P. varians via the Critical Thermal Maximum method (CTmax), using loss of equilibrium as endpoint. We further calculated thermal safety margins for males and females and tested for correlations between upper thermal limits and shrimps' body size. To determine sex-biased variation in P. varians' traits (CTmax, weight and length), we compared trait variation between females and males through the coefficient of variation ratio (lnCVR). Females displayed an average CTmax value 1.8% lower than males (CTmax_females = 37.0 °C vs CTmax_males = 37.7 °C). This finding may be related to the larger body size exhibited by females (156% heavier and 39% larger than males), as both length and weight had a significant effect on CTmax. The high energetic investment of females in offspring may also contribute to the differences recorded in thermal tolerance. Overall, organisms with a smaller body-size displayed a greater tolerance to elevated temperature, thus suggesting that smaller individuals may be positively selected in warmer environments. This selection may result in a reduction of size-at-maturity and shifts in sex ratio, given the sexual dimorphism in body size of shrimps. The thermal safety margin of P. varians was narrow (∼2.2 °C for males and ∼1.5 °C for females), revealing the vulnerability of this species to ocean warming and heatwaves.

Host genotype structures the microbiome of a globally dispersed marine phytoplankton

Microscale interactions between marine phytoplankton and their bacterial microbiomes can influence ecosystem functioning and global biogeochemical cycling through complex exchanges of metabolites and sophisticated ecological processes. Previous investigation of the phytoplankton microbiome has not focused on the role of a host’s underlying genetic background. Through examination of a single phytoplankton species’ microbiome across the global ocean, we found that host genotype strongly influenced microbiome community composition, with associations that potentially persist across generations and ocean basins but assemble rapidly (within days). The long-term association of microbiomes with host genetic background could explain the evolution and maintenance of intricate phytoplankton–bacteria interactions.

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.

Seasonal genotype dynamics of a marine dinoflagellate: Pelagic populations are homogeneous and as diverse as benthic seed banks

Genetic diversity is the basis for evolutionary adaptation and selection under changing environmental conditions. Phytoplankton populations are genotypically diverse, can become genetically differentiated within small spatiotemporal scales and many species form resting stages. Resting stage accumulations in sediments (seed banks) are expected to serve as reservoirs for genetic information, but so far their role in maintaining phytoplankton diversity and in evolution has remained unclear. In this study we used the toxic dinoflagellate Alexandrium ostenfeldii (Dinophyceae) as a model organism to investigate if (i) the benthic seed bank is more diverse than the pelagic population and (ii) the pelagic population is seasonally differentiated. Resting stages (benthic) and plankton (pelagic) samples were collected at a coastal bloom site in the Baltic Sea, followed by cell isolation and genotyping using microsatellite markers (MS) and restriction site associated DNA sequencing (RAD). High clonal diversity (98%–100%) combined with intermediate to low gene diversity (0.58–0.03, depending on the marker) was found. Surprisingly, the benthic and pelagic fractions of the population were equally diverse, and the pelagic fraction was temporally homogeneous, despite seasonal fluctuation of environmental selection pressures. The results of this study suggest that continuous benthic–pelagic coupling, combined with frequent sexual reproduction, as indicated by persistent linkage equilibrium, prevent the dominance of single clonal lineages in a dynamic environment. Both processes harmonize the pelagic with the benthic population and thus prevent seasonal population differentiation. At the same time, frequent sexual reproduction and benthic–pelagic coupling maintain high clonal diversity in both habitats.

Multivariate trait analysis reveals diatom plasticity constrained to a reduced set of biological axes

Trait-based approaches to phytoplankton ecology have gained traction in recent decades as phenotypic traits are incorporated into ecological and biogeochemical models. Here, we use high-throughput phenotyping to explore both intra- and interspecific constraints on trait combinations that are expressed in the cosmopolitan marine diatom genus Thalassiosira. We demonstrate that within Thalassiosira, phenotypic diversity cannot be predicted from genotypic diversity, and moreover, plasticity can create highly divergent phenotypes that are incongruent with taxonomic grouping. Significantly, multivariate phenotypes can be represented in reduced dimensional space using principal component analysis with 77.7% of the variance captured by two orthogonal axes, here termed a 'trait-scape'. Furthermore, this trait-scape can be recovered with a reduced set of traits. Plastic responses to the new environments expanded phenotypic trait values and the trait-scape, however, the overall pattern of response to the new environments was similar between strains and many trait correlations remained constant. These findings demonstrate that trait-scapes can be used to reveal common constraints on multi-trait plasticity in phytoplankton with divergent underlying phenotypes. Understanding how to integrate trait correlational constraints and trade-offs into theoretical frameworks like biogeochemical models will be critical to predict how microbial responses to environmental change will impact elemental cycling now and into the future.

A High-Throughput Assay for Quantifying Phenotypic Traits of Microalgae

High-throughput methods for phenotyping microalgae are in demand across a variety of research and commercial purposes. Many microalgae can be readily cultivated in multi-well plates for experimental studies which can reduce overall costs, while measuring traits from low volume samples can reduce handling. Here we develop a high-throughput quantitative phenotypic assay (QPA) that can be used to phenotype microalgae grown in multi-well plates. The QPA integrates 10 low-volume, relatively high-throughput trait measurements (growth rate, cell size, granularity, chlorophyll a, neutral lipid content, silicification, reactive oxygen species accumulation, and photophysiology parameters: ETRmax, Ik, and alpha) into one workflow. We demonstrate the utility of the QPA on Thalassiosira spp., a cosmopolitan marine diatom, phenotyping six strains in a standard nutrient rich environment (f/2 media) using the full 10-trait assay. The multivariate phenotypes of strains can be simplified into two dimensions using principal component analysis, generating a trait-scape. We determine that traits show a consistent pattern when grown in small volume compared to more typical large volumes. The QPA can thus be used for quantifying traits across different growth environments without requiring exhaustive large-scale culturing experiments, which facilitates experiments on trait plasticity. We confirm that this assay can be used to phenotype newly isolated diatom strains within 4 weeks of isolation. The QPA described here is highly amenable to customisation for other traits or unicellular taxa and provides a framework for designing high-throughput experiments. This method will have applications in experimental evolution, modelling, and for commercial applications where screening of phytoplankton traits is of high importance.

Marine microalgae for outdoor biomass production-A laboratory study simulating seasonal light and temperature for the west coast of Sweden

At Nordic latitudes, year-round outdoor cultivation of microalgae is debatable due to seasonal variations in productivity. Shall the same species/strains be used throughout the year, or shall seasonal-adapted ones be used? To elucidate this, a laboratory study was performed where two out of 167 marine microalgal strains were selected for intended cultivation at the west coast of Sweden. The two local strains belong to Nannochloropsis granulata (Ng) and Skeletonema marinoi (Sm142). They were cultivated in photobioreactors and compared in conditions simulating variations in light and temperature of a year divided into three growth seasons (spring, summer and winter). The strains grew similarly well in summer (and also in spring), but Ng produced more biomass (0.225 vs. 0.066 g DW L^-1 day^-1 ) which was more energy rich (25.0 vs. 16.6 MJ kg^-1 DW). In winter, Sm142 grew faster and produced more biomass (0.017 vs. 0.007 g DW L^-1 day^-1 ), having similar energy to the other seasons. The higher energy of the Ng biomass is attributed to a higher lipid content (40 vs. 16% in summer). The biomass of both strains was richest in proteins (65%) in spring. In all seasons, Sm142 was more effective in removing phosphorus from the cultivation medium (6.58 vs. 4.14 mg L^-1 day^-1 in summer), whereas Ng was more effective in removing nitrogen only in summer (55.0 vs. 30.8 mg L^-1 day^-1 ). Our results suggest that, depending on the purpose, either the same or different local species can be cultivated, and are relevant when designing outdoor studies.

Restoration, conservation and phytoplankton hysteresis

Phytoplankton growth depends not only upon external factors that are not strongly altered by the presence of phytoplankton, such as temperature, but also upon factors that are strongly influenced by activity of phytoplankton, including photosynthetically active radiation, and the availability of the macronutrients carbon, nitrogen, phosphorus and, for some, silicate. Since phytoplankton therefore modify, and to an extent create, their own habitats, established phytoplankton communities can show resistance and resilience to change, including managed changes in nutrient regimes. Phytoplankton blooms and community structures can be predicted from the overall biogeochemical setting and inputs, but restorations may be influenced by the physiological responses of established phytoplankton taxa to nutrient inputs, temperature, second-order changes in illumination and nutrient recycling. In this review we discuss the contributions of phytoplankton ecophysiology to biogeochemical hysteresis and possible effects on community composition in the face of management, conservation or remediation plans.

Biodiversity of marine microbes is safeguarded by phenotypic heterogeneity in ecological traits

Why, contrary to theoretical predictions, do marine microbe communities harbor tremendous phenotypic heterogeneity? How can so many marine microbe species competing in the same niche coexist? We discovered a unifying explanation for both phenomena by investigating a non-cooperative game that interpolates between individual-level competitions and species-level outcomes. We identified all equilibrium strategies of the game. These strategies represent the probability distribution of competitive abilities (e.g. traits) and are characterized by maximal phenotypic heterogeneity. They are also neutral towards each other in the sense that an unlimited number of species can co-exist while competing according to the equilibrium strategies. Whereas prior theory predicts that natural selection would minimize trait variation around an optimum value, here we obtained a mathematical proof that species with maximally variable traits are those that endure. This discrepancy may reflect a disparity between predictions from models developed for larger organisms in contrast to our microbe-centric model. Rigorous mathematics proves that phenotypic heterogeneity is itself a mechanistic underpinning of microbial diversity. This discovery has fundamental ramifications for microbial ecology and may represent an adaptive reservoir sheltering biodiversity in changing environmental conditions.

Resolving cryptic species complexes in marine protists: phylogenetic haplotype networks meet global DNA metabarcoding datasets

Marine protists have traditionally been assumed to be lowly diverse and cosmopolitan. Yet, several recent studies have shown that many protist species actually consist of cryptic complexes of species whose members are often restricted to particular biogeographic regions. Nonetheless, detection of cryptic species is usually hampered by sampling coverage and application of methods (e.g. phylogenetic trees) that are not well suited to identify relatively recent divergence and ongoing gene flow. In this paper, we show how these issues can be overcome by inferring phylogenetic haplotype networks from global metabarcoding datasets. We use the Chaetoceros curvisetus (Bacillariophyta) species complex as study case. Using two complementary metabarcoding datasets (Ocean Sampling Day and Tara Oceans), we equally resolve the cryptic complex in terms of number of inferred species. We detect new hypothetical species in both datasets. Gene flow between most of species is absent, but no barcoding gap exists. Some species have restricted distribution patterns whereas others are widely distributed. Closely related taxa occupy contrasting biogeographic regions, suggesting that geographic and ecological differentiation drive speciation. In conclusion, we show the potential of the analysis of metabarcoding data with evolutionary approaches for systematic and phylogeographic studies of marine protists.

Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms

Diatoms, an evolutionarily successful group of microalgae, display high levels of intraspecific genetic variability in natural populations. However, the contribution of various mechanisms generating such diversity is unknown. Here we estimated the genetic micro-diversity within a natural diatom population and mapped the genomic changes arising within clonally propagated diatom cell cultures. Through quantification of haplotype diversity by next-generation sequencing and amplicon re-sequencing of selected loci, we documented a rapid accumulation of multiple haplotypes accompanied by the appearance of novel protein variants in cell cultures initiated from a single founder cell. Comparison of the genomic changes between mother and daughter cells revealed copy number variation and copy-neutral loss of heterozygosity leading to the fixation of alleles within individual daughter cells. The loss of heterozygosity can be accomplished by recombination between homologous chromosomes. To test this hypothesis, we established an endogenous readout system and estimated that the frequency of interhomolog mitotic recombination was under standard growth conditions 4.2 events per 100 cell divisions. This frequency is increased under environmental stress conditions, including treatment with hydrogen peroxide and cadmium. These data demonstrate that copy number variation and mitotic recombination between homologous chromosomes underlie clonal variability in diatom populations. We discuss the potential adaptive evolutionary benefits of the plastic response in the interhomolog mitotic recombination rate, and we propose that this may have contributed to the ecological success of diatoms.

Environmental impacts on single-cell variation within a ubiquitous diatom: The role of growth rate

Morphological and physiological characteristics of phytoplankton cells are highly sensitive to changes in environmental conditions and, in turn, influence the dynamics of phytoplankton populations and communities. To cope with environmental change, trait variability and phenotypic plasticity may play an important role. Since global change comprises simultaneous changes in abiotic parameters, we assessed the impact of multiple drivers on functional traits of the diatom Thalassiosira (Conticribra) weissflogii by manipulating concurrently temperature, pCO₂, and dissolved nitrogen:phosphorus (N:P) ratio. We tested three scenarios: ambient (ambient temperature and atmospheric pCO₂; 16 N:P ratio), moderate future scenario (+1.5°C and 800 ppm CO₂; 25 N:P ratio), and more severe future scenario (+3°C and 1000 ppm CO₂; 25 N:P ratio). We applied flow cytometry to measure on single-cell levels to investigate trait variability and phenotypic plasticity within one strain of diatoms. Growth rates differed significantly between the treatments and were strongly correlated with cell size and cellular chlorophyll a content. We observed a negative correlation of growth rate with chlorophyll a variability among single strain populations and a negative correlation with the phenotypic plasticity of cell size, i.e. when growth rates were higher, the cell size cell-to-cell variability within cultures was lower. Additionally, the phenotypic plasticity in cell size was lower under the global change scenarios. Overall, our study shows that multiple traits are interlinked and driven by growth rate and that this interconnection may partly be shaped by environmental factors.

Revealing environmentally driven population dynamics of an Arctic diatom using a novel microsatellite PoolSeq barcoding approach

Ecological stability under environmental change is determined by both interspecific and intraspecific processes. Particularly for planktonic microorganisms, it is challenging to follow intraspecific dynamics over space and time. We propose a new method, microsatellite PoolSeq barcoding (MPB), for tracing allele frequency changes in protist populations. We successfully applied this method to experimental community incubations and field samples of the diatom Thalassiosira hyalina from the Arctic, a rapidly changing ecosystem. Validation of the method found compelling accuracy in comparison with established genotyping approaches within different diversity contexts. In experimental and environmental samples, we show that MPB can detect meaningful patterns of population dynamics, resolving allelic stability and shifts within a key diatom species in response to experimental treatments as well as different bloom phases and years. Through our novel MPB approach, we produced a large dataset of populations at different time-points and locations with comparably little effort. Results like this can add insights into the roles of selection and plasticity in natural protist populations under stable experimental but also variable field conditions. Especially for organisms where genotype sampling remains challenging, MPB holds great potential to efficiently resolve eco-evolutionary dynamics and to assess the mechanisms and limits of resilience to environmental stressors.

Phenotypic trait variability as an indication of adaptive capacity in a cosmopolitan marine diatom

Determining the adaptive capacity of marine phytoplankton is important in predicting changes in phytoplankton responses to ocean warming. Phytoplankton may consist of high levels of standing phenotypic and genetic variability, the basis of rapid evolution; however, few studies have quantified trait variability within and amongst closely related diatom species. Using 35 clonal cultures of the ubiquitous marine diatom Leptocylindrus isolated from six locations, spanning 2000 km of the south-eastern Australian coastline, we found evidence of significant intraspecific morphological and metabolic trait variability, which for 8 of 9 traits (growth rate, biovolume, C:N, silica deposition, silica incorporation rate, chl-a, and photosynthetic efficiency under dark adapted, growth irradiance, and high-light adaptation) were greater within a species than between species. Moreover, only two traits revealed a latitudinal trend with strains isolated from lower latitudes showing significantly higher silicification rates and protein:lipid content compared to their higher latitude counterparts. These data mirror recent studies on diatom intraspecific genetic diversity, which has found comparable levels of genetic diversity at a single site to those thousands of kilometres apart, and provide evidence of a functional role of diatom diversity that will allow for rapid adaptation via ecological selection on standing variation in response to changing conditions.

Lake-depth related pattern of genetic and morphological diatom diversity in boreal Lake Bolshoe Toko, Eastern Siberia

Large, old and heterogenous lake systems are valuable sources of biodiversity. The analysis of current spatial variability within such lakes increases our understanding of the origin and establishment of biodiversity. The environmental sensitivity and the high taxonomic richness of diatoms make them ideal organisms to investigate intra-lake variability. We investigated modern intra-lake diatom diversity in the large and old sub-arctic Lake Bolshoe Toko in Siberia. Our study uses diatom-specific metabarcoding, applying a short rbcL marker combined with next-generation sequencing and morphological identification to analyse the diatom diversity in modern sediment samples of 17 intra-lake sites. We analysed abundance-based compositional taxonomic diversity and generic phylogenetic diversity to investigate the relationship of diatom diversity changes with water depth. The two approaches show differences in taxonomic identification and alpha diversity, revealing a generally higher diversity with the genetic approach. With respect to beta diversity and ordination analyses, both approaches result in similar patterns. Water depth or related lake environmental conditions are significant factors influencing intra-lake diatom patterns, showing many significant negative correlations between alpha and beta diversity and water depth. Further, one near-shore and two lagoon lake sites characterized by low (0-10m) and medium (10-30m) water depth are unusual with unique taxonomic compositions. At deeper (>30m) water sites we identified strongest phylogenetic clustering in Aulacoseira, but generally much less in Staurosira, which supports that water depth is a strong environmental filter on the Aulacoseira communities. Our study demonstrates the utility of combining analyses of genetic and morphological as well as phylogenetic diversity to decipher compositional and generic phylogenetic patterns, which are relevant in understanding intra-lake heterogeneity as a source of biodiversity in the sub-arctic glacial Lake Bolshoe Toko.

A genomics approach reveals the global genetic polymorphism, structure, and functional diversity of ten accessions of the marine model diatom Phaeodactylum tricornutum

Diatoms emerged in the Mesozoic period and presently constitute one of the main primary producers in the world's ocean and are of a major economic importance. In the current study, using whole genome sequencing of ten accessions of the model diatom Phaeodactylum tricornutum, sampled at broad geospatial and temporal scales, we draw a comprehensive landscape of the genomic diversity within the species. We describe strong genetic subdivisions of the accessions into four genetic clades (A-D) with constituent populations of each clade possessing a conserved genetic and functional makeup, likely a consequence of the limited dispersal of P. tricornutum in the open ocean. We further suggest dominance of asexual reproduction across all the populations, as implied by high linkage disequilibrium. Finally, we show limited yet compelling signatures of genetic and functional convergence inducing changes in the selection pressure on many genes and metabolic pathways. We propose these findings to have significant implications for understanding the genetic structure of diatom populations in nature and provide a framework to assess the genomic underpinnings of their ecological success and impact on aquatic ecosystems where they play a major role. Our work provides valuable resources for functional genomics and for exploiting the biotechnological potential of this model diatom species.

Digest: Fossils, evolutionary models, and diatoms

Adaptation to new environments by organisms has been and is of much interest in scientific research. Does the diatom lineage Stephanodiscus niagarae/yellowstonensis follow models compatible with adaptation to a new environment? Voje (2019) shows that in this lineage, none of the current models tested-decelerated evolution, Ornstein-Uhlenbeck, and random walk-fully explained how multiple traits-valve diameter, the number of costae, and spines per valve-evolved. More accurate models are needed to understand the evolutionary history of these diatoms.

Differential Influence of Life Cycle on Growth and Toxin Production of three Pseudo-nitzschia Species (Bacillariophyceae)

We used a multistrain approach to study the intra- and interspecific variability of the growth rates of three Pseudo-nitzschia species - P. australis, P. fraudulenta, and P. pungens - and of their domoic acid (DA) production. We carried out mating and batch experiments to investigate the respective effects of strain age and cell size, and thus the influence of their life cycle on the physiology of these species. The cell size - life cycle relationship was characteristic of each species. The influence of age and cell size on the intraspecific variability of growth rates suggests that these characteristics should be considered cautiously for the strains used in physiological studies on Pseudo-nitzschia species. The results from all three species do not support the hypothesis of a decrease in DA production with time since isolation from natural populations. In P. australis, the cellular DA content was rather a function of cell size. More particularly, cells at the gametangia stage of their life cycle contained up to six times more DA than smaller or larger cells incapable of sexual reproduction. These findings reveal a link between P. australis life cycle and cell toxicity. This suggest that life cycle dynamics in Pseudo-nitzschia natural populations may influence bloom toxicity.

Company matters: The presence of other genotypes alters traits and intraspecific selection in an Arctic diatom under climate change

Arctic phytoplankton and their response to future conditions shape one of the most rapidly changing ecosystems on the planet. We tested how much the phenotypic responses of strains from the same Arctic diatom population diverge and whether the physiology and intraspecific composition of multistrain populations differs from expectations based on single strain traits. To this end, we conducted incubation experiments with the diatom Thalassiosira hyalina under present-day and future temperature and pCO2 treatments. Six fresh isolates from the same Svalbard population were incubated as mono- and multistrain cultures. For the first time, we were able to closely follow intraspecific selection within an artificial population using microsatellites and allele-specific quantitative PCR. Our results showed not only that there is substantial variation in how strains of the same species cope with the tested environments but also that changes in genotype composition, production rates, and cellular quotas in the multistrain cultures are not predictable from monoculture performance. Nevertheless, the physiological responses as well as strain composition of the artificial populations were highly reproducible within each environment. Interestingly, we only detected significant strain sorting in those populations exposed to the future treatment. This study illustrates that the genetic composition of populations can change on very short timescales through selection from the intraspecific standing stock, indicating the potential for rapid population level adaptation to climate change. We further show that individuals adjust their phenotype not only in response to their physicochemical but also to their biological surroundings. Such intraspecific interactions need to be understood in order to realistically predict ecosystem responses to global change.

Intraspecific Variation in Protists: Clues for Microevolution from Poteriospumella lacustris (Chrysophyceae)

Species delimitation in protists is still a challenge, attributable to the fact that protists are small, difficult to observe and many taxa are poor in morphological characters, whereas most current phylogenetic approaches only use few marker genes to measure genetic diversity. To address this problem, we assess genome-level divergence and microevolution in strains of the protist Poteriospumella lacustris, one of the first free-living, nonmodel organisms to study genome-wide intraspecific variation. Poteriospumella lacustris is a freshwater protist belonging to the Chrysophyceae with an assumed worldwide distribution. We examined three strains from different geographic regions (New Zealand, China, and Austria) by sequencing their genomes with the Illumina and PacBio platforms. The assembled genomes were small with 49-55 Mb but gene-rich with 16,000-19,000 genes, of which ∼8,000 genes could be assigned to functional categories. At least 68% of these genes were shared by all three species. Genetic variation occurred predominantly in genes presumably involved in ecological niche adaptation. Most surprisingly, we detected differences in genome ploidy between the strains (diploidy, triploidy, and tetraploidy). In analyzing intraspecific variation, several mechanisms of diversification were identified including SNPs, change of ploidy and genome size reduction.

Critical light-related gene expression varies in two different strains of the dinoflagellate Karlodinium veneficum in response to the light spectrum and light intensity

The toxic dinoflagellate Karlodinium veneficum is widely distributed in cosmopolitan estuaries and is responsible for massive fish mortality worldwide. Intraspecific biodiversity is important for the spread to various habitats, interspecific competition to dominate a population, and bloom formation and density maintenance. Strategies for light adaptation may help determine the ecological niches of different ecotypes. However, the mechanism of phenotypic biodiversity is still unclear. In this study, intraspecific differences in genetic regulatory mechanisms in response to varied light intensities and qualities were comparatively researched on two different strains isolated from coastal areas of the East China Sea, namely, GM2 and GM3. In GM2, the expression of genes in the Calvin cycle, namely, rbcL and SBPase, and a light-related gene that correlated with cellular motility, rhodopsin, were significantly inhibited under high light intensities. Thus, this strain was adapted to low light. In contrast, the gene expression levels were promoted by high light conditions in GM3. These upregulated genes in the GM3 strain probably compensated for the negative effects on the maximum quantum yields of PSII (Fv/Fm) under high light stress, which inhibited both strains, enabling GM3 to maintain a constant growth rate. Thus, this strain was adapted to high light. Compared with white light, monochromatic blue light had negative effects on Fv/Fm and the relative electron transfer rate (ETR) in both strains. Under blue light, gene expression levels of rbcL and SBPase in GM2 were inhibited; in contrast, the levels of these genes, especially rbcL, were promoted in GM3. rbcL was significantly upregulated in the blue light groups. Monochromatic red light promoted rhodopsin gene expression in the two strains in a similar manner. These intraspecific diverse responses to light play important roles in the motor characteristics, diel vertical migration, interspecific relationships and photosynthetic or phagotrophic activities of K. veneficum and can determine the population distribution, population maintenance and bloom formation.

Phytoplankton adaptation in ecosystem models

We compare two different approaches to model adaptation of phytoplankton through trait value changes. Both consider mutation and selection (MuSe) but differ with respect to the underlying conceptual framework. The first one (MuSe-IBM) explicitly considers a population of individuals that are subject to random mutation during cell division. The second is a deterministic multi-compartment model (MuSe-MCM) that considers numerous genotypes of the population and where mutations are treated as a transfer of biomass between neighboring genotypes (i.e., a diffusion of characteristics in trait space). Focusing on the adaptation of optimal temperature, we show model results for different scenarios: a sudden change in environmental temperature, a seasonal variation and high frequency fluctuations. In addition, we investigate the effect of different shapes of thermal reaction norms as well as the role of alternating growth and resting phases on the adaptation process. For all cases, the differences between MuSe-IBM and MuSe-MCM are found to be negligible. Both models produce a number of well-known and plausible features. While the IBM has the advantage of including more mechanistic (i.e., probabilistic) processes, the MCM is much less computationally demanding and therefore suitable for implementation in three-dimensional ecosystem models.

Influence of sudden salinity variation on the physiology and domoic acid production by two strains of Pseudo-nitzschia australis

Several coastal countries including France have experienced serious and increasing problems related to Pseudo-nitzschia toxic blooms. These toxic blooms occur in estuarine and coastal waters potentially subject to fluctuations in salinity. In this study, we document for the first time the viability, growth, photosynthetic efficiency, and toxin production of two strains of Pseudo-nitzschia australis grown under conditions with sudden salinity changes. Following salinity variation, the two strains survived over a restricted salinity range of 30-35, with favorable physiological responses, as the growth, effective quantum yield and toxin content were high compared to the other conditions. In addition, high cellular quotas of domoic acid (DA) were observed at a salinity of 40 for the strain IFR-PAU-16.1 in comparison with the other strain IFR-PAU-16.2 where the cell DA content was directly released into the medium. On the other hand, the osmotic stress imposed at lower salinities, 20 and 10, resulted in cell lysis and a sudden DA leakage in the medium. Intra-specific variability was observed in growth and toxin production, with the strain IFR-PAU-16.1 apparently able to withstand higher salinities than the strain IFR-PAU-16.2. On the whole, DA does not appear to act as an osmolyte in response to sudden salinity changes. Since most of the shellfish harvesting areas of bivalve molluscs in France are located in areas where the salinity generally varies between 30 and 35, Pseudo-nitzschia australis blooms might potentially impact public health and commercial shellfish resources in these places.

Diatom diversity through HTS-metabarcoding in coastal European seas

Diatoms constitute a diverse lineage of unicellular organisms abundant and ecologically important in aquatic ecosystems. Compared to other protists, their biology and taxonomy are well-studied, offering the opportunity to combine traditional approaches and new technologies. We examined a dataset of diatom 18S rRNA- and rDNA- (V4 region) reads from different plankton size-fractions and sediments from six European coastal marine sites, with the aim of identifying peculiarities and commonalities with respect to the whole protistan community. Almost all metabarcodes (99.6%) were assigned to known genera (121) and species (236), the most abundant of which were those already known from classic studies and coincided with those seen in light microscopy. rDNA and rRNA showed comparable patterns for the dominant taxa, but rRNA revealed a much higher diversity particularly in the sediment communities. Peculiar to diatoms is a tight bentho-pelagic coupling, with many benthic or planktonic species colonizing both water column and sediments and the dominance of planktonic species in both habitats. Overall metabarcoding results reflected the marked specificity of diatoms compared to other protistan groups in terms of morphological and ecological characteristics, at the same time confirming their great potential in the description of protist communities.

Growth of a common planktonic diatom quantified using solid medium culturing

The ability to grow on solid culture medium is a pre-requisite for a successful microbial genetic model organism. Skeletonema marinoi, a bloom-forming, planktonic marine microalga, is widely used in ecological, evolutionary and population genetics studies. We have tested and confirmed the ability of this common organism to grow on solid culture medium (agar) under experimentally manipulated conditions. We established a protocol for quantifying growth characteristics - length of lag phase, growth rate, maximum biomass yield - on agar medium. The procedure was tested under experimental treatments and the resulting growth changes correlated with those observed in standard liquid culture. The ability to grow on solid medium broadens the use of S. marinoi as a molecular model, where agar is routinely used for various purposes (growth, selection, storage); and the possibility to quantify colony growth opens the way for high throughput, automated, or semi-automated phenotyping solutions.