Metapopulation Structure in the Planktonic Diatom Ditylum brightwellii (Bacillariophyceae)

Whole-genome scanning reveals environmental selection mechanisms that shape diversity in populations of the epipelagic diatom Chaetoceros

Diatoms form a diverse and abundant group of photosynthetic protists that are essential players in marine ecosystems. However, the microevolutionary structure of their populations remains poorly understood, particularly in polar regions. Exploring how closely related diatoms adapt to different environments is essential given their short generation times, which may allow rapid adaptations, and their prevalence in marine regions dramatically impacted by climate change, such as the Arctic and Southern Oceans. Here, we address genetic diversity patterns in Chaetoceros, the most abundant diatom genus and one of the most diverse, using 11 metagenome-assembled genomes (MAGs) reconstructed from Tara Oceans metagenomes. Genome-resolved metagenomics on these MAGs confirmed a prevalent distribution of Chaetoceros in the Arctic Ocean with lower dispersal in the Pacific and Southern Oceans as well as in the Mediterranean Sea. Single-nucleotide variants identified within the different MAG populations allowed us to draw a landscape of Chaetoceros genetic diversity and revealed an elevated genetic structure in some Arctic Ocean populations. Gene flow patterns of closely related Chaetoceros populations seemed to correlate with distinct abiotic factors rather than with geographic distance. We found clear positive selection of genes involved in nutrient availability responses, in particular for iron (e.g., ISIP2a, flavodoxin), silicate, and phosphate (e.g., polyamine synthase), that were further supported by analysis of Chaetoceros transcriptomes. Altogether, these results highlight the importance of environmental selection in shaping diatom diversity patterns and provide new insights into their metapopulation genomics through the integration of metagenomic and environmental data.

The anchoring effect-long-term dormancy and genetic population structure

Understanding the genetic structure of populations is key to revealing past and present demographic and evolutionary processes in a species. In the past decade high genetic differentiation has been observed in many microbial species challenging the previous view of cosmopolitan distribution. Populations have displayed high genetic differentiation, even at small spatial scales, despite apparent high dispersal. Numerous species of microalgae have a life-history strategy that includes a long-term resting stage, which can accumulate in sediments and serve as refuge during adverse conditions. It is presently unclear how these seed banks affect the genetic structure of populations in aquatic environments. Here we provide a conceptual framework, using a simple model, to show that long-term resting stages have an anchoring effect on populations leading to increased genetic diversity and population differentiation in the presence of gene flow. The outcome that species with resting stages have a higher degree of genetic differentiation compared to species without, is supported by empirical data obtained from a systematic literature review. With this work we propose that seed banks in aquatic microalgae play an important role in the contradicting patterns of gene flow, and ultimately the adaptive potential and population dynamics in species with long-term resting stages.

Marteilia refringens and Marteilia pararefringens sp. nov. are distinct parasites of bivalves and have different European distributions

Marteilia refringens causes marteiliosis in oysters, mussels and other bivalve molluscs. This parasite previously comprised two species, M. refringens and Marteilia maurini, which were synonymized in 2007 and subsequently referred to as M. refringens ‘O-type’ and ‘M-type’. O-type has caused mass mortalities of the flat oyster Ostrea edulis. We used high throughput sequencing and histology to intensively screen flat oysters and mussels (Mytilus edulis) from the UK, Sweden and Norway for infection by both types and to generate multi-gene datasets to clarify their genetic distinctiveness. Mussels from the UK, Norway and Sweden were more frequently polymerase chain reaction (PCR)-positive for M-type (75/849) than oysters (11/542). We did not detect O-type in any northern European samples, and no histology-confirmed Marteilia-infected oysters were found in the UK, Norway and Sweden, even where co-habiting mussels were infected by the M-type. The two genetic lineages within ‘M. refringens’ are robustly distinguishable at species level. We therefore formally define them as separate species: M. refringens (previously O-type) and Marteilia pararefringens sp. nov. (M-type). We designed and tested new Marteilia-specific PCR primers amplifying from the 3’ end of the 18S rRNA gene through to the 5.8S gene, which specifically amplified the target region from both tissue and environmental samples.

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

Intraspecific variation in diatoms has been shown to play a key role in species' responses to several important environmental factors such as light, salinity, temperature and nutrients. Furthermore, modelling efforts indicate that this variation within species extends bloom periods, and likely provides sufficient variability in competitive interactions between species under hydrographically variable conditions. The intraspecific variation most likely corresponds to optimal fitness in temporary microhabitats and may help to explain the paradox of the plankton. Here, we examine the implications of intraspecific variation for the ecology and success of diatoms in general and emphasize the potential implications for our understanding of carbon metabolism in these important organisms. Additionally, data from palaeoecological studies have the potential for evaluating genetic variation through past climate changes, going thousands of years back in time. We suggest pathways for future research including the adoption of multiple strains of individual species into studies of diatom carbon metabolism, to refine our understanding of the variation within and between species, and the inclusion of experimental evolution as a tool for understanding potential evolutionary responses of diatom carbon metabolism to climate change.This article is part of the themed issue 'The peculiar carbon metabolism in diatoms'.

Morphological diversity and phylogeny of the diatom genus Entomoneis (Bacillariophyta) in marine plankton: six new species from the Adriatic Sea

The diatom genus Entomoneis is known from the benthos and plankton of marine, brackish, and freshwaters. Entomoneis includes diatoms with a bilobate keel elevated above the valve surface, a sigmoid canal raphe, and numerous girdle bands. Owing mostly to the scarcity of molecular data for a diverse set of species, the phylogeny of Entomoneis has not been investigated in depth. The few previous studies that included Entomoneis were focused on broader questions and the available data were from a small number of either unidentified Entomoneis or well-known species (e.g., E. paludosa). Since the first description of new species combining both molecular and morphological characters (E. tenera), we have continued to cultivate and investigate Entomoneis in the plankton of the Adriatic Sea. Combined multigene phylogeny (SSU rDNA sequences, rbcL, and psbC genes) and morphological observations (LM, SEM and TEM) revealed six new Entomoneis species supported by phylogenetic and morphological data: E. pusilla, E. gracilis, E. vilicicii, E. infula, E. adriatica, and E. umbratica. The most important morphological features for species delineation were cell shape, the degree and mode of torsion, valve apices, the appearance and structure of the transition between keel and valve body, the ultrastructure and the shape of the girdle bands, and the arrangement and density of perforations along the valve and valvocopulae. Our results highlight the underappreciated diversity of Entomoneis and call for a more in-depth morphological and molecular investigation of this genus especially in planktonic habitats.

Phylogeography of the freshwater raphidophyte Gonyostomum semen confirms a recent expansion in northern Europe by a single haplotype

Gonyostmum semen is a freshwater raphidophyte that has increased in occurrence and abundance in several countries in northern Europe since the 1980s. More recently, the species has expanded rapidly also in north-eastern Europe, and it is frequently referred to as invasive. To better understand the species history, we have explored the phylogeography of G. semen using strains from northern Europe, United States, and Japan. Three regions of the ribosomal RNA gene (small subunit [SSU], internal transcribed spacer [ITS] and large subunit [LSU]) and one mitochondrial DNA marker (cox1) were analyzed. The SSU and partial LSU sequences were identical in all strains, confirming that they belong to the same species. The ITS region differentiated the American from the other strains, but showed high intra-strain variability. In contrast, the mitochondrial marker cox1 showed distinct differences between the European, American, and Japanese strains. Interestingly, only one cox1 haplotype was detected in European strains. The overall low diversity and weak geographic structure within northern European strains supported the hypothesis of a recent invasion of new lakes by G. semen. Our data also show that the invasive northern European lineage is genetically distinct from the lineages from the other continents. Finally, we concluded that the mitochondrial cox1 was the most useful marker in determining large-scale biogeographic patterns in this species.

Intraspecific facilitation by allelochemical mediated grazing protection within a toxigenic dinoflagellate population

Dinoflagellates are a major cause of harmful algal blooms (HABs), with consequences for coastal marine ecosystem functioning and services. Alexandrium fundyense (previously Alexandrium tamarense) is one of the most abundant and widespread toxigenic species in the temperate Northern and Southern Hemisphere and produces paralytic shellfish poisoning toxins as well as lytic allelochemical substances. These bioactive compounds may support the success of A. fundyense and its ability to form blooms. Here we investigate the impact of grazing on monoclonal and mixed set-ups of highly (Alex2) and moderately (Alex4) allelochemically active A. fundyense strains and a non-allelochemically active conspecific (Alex5) by the heterotrophic dinoflagellate Polykrikos kofoidii. While Alex4 and particularly Alex5 were strongly grazed by P. kofoidii when offered alone, both strains grew well in the mixed assemblages (Alex4 + Alex5 and Alex2 + Alex5). Hence, the allelochemical active strains facilitated growth of the non-active strain by protecting the population as a whole against grazing. Based on our results, we argue that facilitation among clonal lineages within a species may partly explain the high genotypic and phenotypic diversity of Alexandrium populations. Populations of Alexandrium may comprise multiple cooperative traits that act in concert with intraspecific facilitation, and hence promote the success of this notorious HAB species.

Hidden diversity in the freshwater planktonic diatom Asterionella formosa

Many freshwater and marine algal species are described as having cosmopolitan distributions. Whether these widely distributed morphologically similar algae also share a similar gene pool remains often unclear. In the context of island biogeography theory, stronger spatial isolation deemed typical of freshwater lakes should restrict gene flow and lead to higher genetic differentiation among lakes. Using nine microsatellite loci, we investigate the genetic diversity of a widely distributed freshwater planktonic diatom, Asterionella formosa, across different lakes in Switzerland and the Netherlands. We applied a hierarchical spatial sampling design to determine the geographical scale at which populations are structured. A subset of the isolates was additionally analysed using amplified fragment length polymorphism (AFLP) markers. Our results revealed complex and unexpected population structure in A. formosa with evidence for both restricted and moderate to high gene flow at the same time. Different genetic markers (microsatellites and AFLPs) analysed with a variety of multivariate methods consistently revealed that genetic differentiation within lakes was much stronger than among lakes, indicating the presence of cryptic species within A. formosa. We conclude that the hidden diversity found in this study is expected to have implications for the further use of A. formosa in biogeographical, conservation and ecological studies. Further research using species-level phylogenetic markers is necessary to place the observed differentiation in an evolutionary context of speciation.

Temporal changes in population structure of a marine planktonic diatom

A prevailing question in phytoplankton research addresses changes of genetic diversity in the face of huge population sizes and apparently unlimited dispersal capabilities. We investigated population genetic structure of the pennate planktonic marine diatom Pseudo-nitzschia multistriata at the LTER station MareChiara in the Gulf of Naples (Italy) over four consecutive years and explored possible changes over seasons and from year to year. A total of 525 strains were genotyped using seven microsatellite markers, for a genotypic diversity of 75.05%, comparable to that found in other Pseudo-nitzschia species. Evidence from Bayesian clustering analysis (BA) identified two genetically distinct clusters, here interpreted as populations, and several strains that could not be assigned with ≥90% probability to either population, here interpreted as putative hybrids. Principal Component Analysis (PCA) recovered these two clusters in distinct clouds with most of the putative hybrids located in-between. Relative proportions of the two populations and the putative hybrids remained similar within years, but changed radically between 2008 and 2009 and between 2010 and 2011, when the 2008-population apparently became the dominant one again. Strains from the two populations are inter-fertile, and so is their offspring. Inclusion of genotypes of parental strains and their offspring shows that the majority of the latter could not be assigned to any of the two parental populations. Therefore, field strains classified by BA as the putative hybrids could be biological hybrids. We hypothesize that P. multistriata population dynamics in the Gulf of Naples follows a meta-population-like model, including establishment of populations by cell inocula at the beginning of each growth season and remixing and dispersal governed by moving and mildly turbulent water masses.

Exploration of nuclear DNA markers for population structure assessment in the desmid Micrasterias rotata (Zygnematophyceae, Streptophyta)

Freshwater green microalgae are diverse and widely distributed across the globe, yet the population structuring of these organisms is poorly understood. We assessed the degree of genetic diversity and differentiation of the desmid species, Micrasterias rotata. First, we compared the sequences of four nuclear regions (actin, gapC1, gapC2, and oee1) in 25 strains and selected the gapC1 and actin regions as the most appropriate markers for population structure assessment in this species. Population genetic structure was subsequently analyzed, based on seven populations from the Czech Republic and Ireland. Hudson's Snn statistics indicated that nearest-neighbor sequences occurred significantly more frequently within geographical populations than within the wider panmictic population. Moreover, Irish populations consistently showed higher genetic diversity than the Czech samples. These results are in accordance with the unbalanced distribution of alleles in many land plant species; however, the large genetic diversity in M. rotata differs from levels of genetic diversity found in most land plants.

High evolutionary potential of marine zooplankton

Open ocean zooplankton often have been viewed as slowly evolving species that have limited capacity to respond adaptively to changing ocean conditions. Hence, attention has focused on the ecological responses of zooplankton to current global change, including range shifts and changing phenology. Here, we argue that zooplankton also are well poised for evolutionary responses to global change. We present theoretical arguments that suggest plankton species may respond rapidly to selection on mildly beneficial mutations due to exceptionally large population size, and consider the circumstantial evidence that supports our inference that selection may be particularly important for these species. We also review all primary population genetic studies of open ocean zooplankton and show that genetic isolation can be achieved at the scale of gyre systems in open ocean habitats (100s to 1000s of km). Furthermore, population genetic structure often varies across planktonic taxa, and appears to be linked to the particular ecological requirements of the organism. In combination, these characteristics should facilitate adaptive evolution to distinct oceanographic habitats in the plankton. We conclude that marine zooplankton may be capable of rapid evolutionary as well as ecological responses to changing ocean conditions, and discuss the implications of this view. We further suggest two priority areas for future research to test our hypothesis of high evolutionary potential in open ocean zooplankton, which will require (1) assessing how pervasive selection is in driving population divergence and (2) rigorously quantifying the spatial and temporal scales of population differentiation in the open ocean. Recent attention has focused on the ecological responses of open ocean zooplankton to current global change, including range shifts and changing phenology. Here, we argue that marine zooplankton also are well poised for evolutionary responses to global change.

Different transcriptional responses of heat shock protein 70/90 in the marine diatom Ditylum brightwellii exposed to metal compounds and endocrine-disrupting chemicals

Environmental hazard assessments using diatoms have been well documented; however, their molecular toxicology has not been sufficiently studied. In this study, we characterized heat shock protein (HSP) 70/90 of the diatom Ditylum brightwellii (Db) and evaluated their transcriptional profiles in response to various environmental stresses (e.g., thermal shocks and metal and non-metal pollutants). Putative DbHSP70 (658aa, 71.7 kDa) and DbHSP90 (707aa, 80.2 kDa) proteins had conserved HSP family motifs but different C-terminus motifs, that is, "EEVD" in DbHSP70 and "MEEVD" in DbHSP90. Phylogenetic analyses of both proteins showed that D. brightwellii was well clustered with other diatoms. Real-time PCR analysis showed that thermal stress considerably upregulated DbHSP70 and DbHSP90. As for chemical pollutants, DbHSP70 greatly responded to CuSO4 and NiSO4 exposure, but not CuCl2 or NiCl2. However, DbHSP90 was significantly upregulated by all the metal compounds tested (CuSO4, NiSO4, CuCl2, and NiCl2). Strikingly, the expression of both genes was not induced by the organic pollutants tested, such as endocrine-disrupting chemicals. These data suggest that DbHSP70 and DbHSP90 are differentially involved in the defense response against various environmental stressors. Moreover, metal toxicity may be specifically affected by the conjugated anion in the metal compounds (e.g., SO4(2-) and Cl(-)).

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.

Genetic diversity within and genetic differentiation between blooms of a microalgal species

The field of genetic diversity in protists, particularly phytoplankton, is under expansion. However, little is known regarding variation in genetic diversity within populations over time. The aim of our study was to investigate intrapopulation genetic diversity and genetic differentiation in the freshwater bloom-forming microalga Gonyostomum semen (Raphidophyceae). The study covered a 2-year period including all phases of the bloom. Amplified fragment length polymorphism (AFLP) was used to determine the genetic structure and diversity of the population. Our results showed a significant differentiation between samples collected during the two blooms from consecutive years. Also, an increase of gene diversity and a loss of differentiation among sampling dates were observed over time within a single bloom. The latter observations may reflect the continuous germination of cysts from the sediment. The life cycle characteristics of G. semen, particularly reproduction and recruitment, most likely explain a high proportion of the observed variation. This study highlights the importance of the life cycle for the intraspecific genetic diversity of microbial species, which alternates between sexual and asexual reproduction.

Genome size differentiates co-occurring populations of the planktonic diatom Ditylum brightwellii (Bacillariophyta)

BACKGROUND

Diatoms are one of the most species-rich groups of eukaryotic microbes known. Diatoms are also the only group of eukaryotic micro-algae with a diplontic life history, suggesting that the ancestral diatom switched to a life history dominated by a duplicated genome. A key mechanism of speciation among diatoms could be a propensity for additional stable genome duplications. Across eukaryotic taxa, genome size is directly correlated to cell size and inversely correlated to physiological rates. Differences in relative genome size, cell size, and acclimated growth rates were analyzed in isolates of the diatom Ditylum brightwellii. Ditylum brightwellii consists of two main populations with identical 18s rDNA sequences; one population is distributed globally at temperate latitudes and the second appears to be localized to the Pacific Northwest coast of the USA. These two populations co-occur within the Puget Sound estuary of WA, USA, although their peak abundances differ depending on local conditions.

RESULTS

All isolates from the more regionally-localized population (population 2) possessed 1.94 +/- 0.74 times the amount of DNA, grew more slowly, and were generally larger than isolates from the more globally distributed population (population 1). The ITS1 sequences, cell sizes, and genome sizes of isolates from New Zealand were the same as population 1 isolates from Puget Sound, but their growth rates were within the range of the slower-growing population 2 isolates. Importantly, the observed genome size difference between isolates from the two populations was stable regardless of time in culture or the changes in cell size that accompany the diatom life history.

CONCLUSIONS

The observed two-fold difference in genome size between the D. brightwellii populations suggests that whole genome duplication occurred within cells of population 1 ultimately giving rise to population 2 cells. The apparent regional localization of population 2 is consistent with a recent divergence between the populations, which are likely cryptic species. Genome size variation is known to occur in other diatom genera; we hypothesize that genome duplication may be an active and important mechanism of genetic and physiological diversification and speciation in diatoms.