Limits to gene flow in a cosmopolitan marine planktonic diatom

Distinct biogeographical patterns in snail gastrointestinal tract bacterial communities compared with sediment and water

The factors that influence the distribution of bacterial community composition are not well understood. The role of geographical patterns, which suggest limited dispersal, is still a topic of debate. Bacteria associated with hosts face unique dispersal challenges as they often rely on their hosts, which provide specific environments for their symbionts. In this study, we examined the effect of biogeographic distances on the bacterial diversity and composition of bacterial communities in the gastrointestinal tract of Ampullaceana balthica. We compared the effects on the host-associated bacterial community to those on bacterial communities in water and sediment. This comparison was made using 16S ribosomal RNA gene sequencing. We found that the bacterial communities we sampled in Estonia, Denmark, and Northern Germany varied between water, sediment, and the gastrointestinal tract. They also varied between countries within each substrate. This indicates that the type of substrate is a dominant factor in determining bacterial community composition. We separately analyzed the turnover rates of water, sediment, and gastrointestinal bacterial communities over increasing geographic distances. We observed that the turnover rate was lower for gastrointestinal bacterial communities compared to water bacterial communities. This implies that the composition of gastrointestinal bacteria remains relatively stable over distances, while water bacterial communities exhibit greater variability. However, the gastrointestinal tract had the lowest percentage of country-specific amplicon sequence variants, suggesting bacterial colonization from local bacterial communities. Since the overlap between the water and gastrointestinal tract was highest, it appears that the gastrointestinal bacterial community is colonized by the water bacterial community. Our study confirmed that biogeographical patterns in host-associated communities differ from those in water and sediment bacterial communities. These host-associated communities consist of numerous facultative symbionts derived from the water bacterial community.

Multiannual patterns of genetic structure and mating type ratios highlight the complex bloom dynamics of a marine planktonic diatom

Understanding the genetic structure of populations and the processes responsible for its spatial and temporal dynamics is vital for assessing species' adaptability and survival in changing environments. We investigate the genetic fingerprinting of blooming populations of the marine diatom Pseudo-nitzschia multistriata in the Gulf of Naples (Mediterranean Sea) from 2008 to 2020. Strains were genotyped using microsatellite fingerprinting and natural samples were also analysed with Microsatellite Pool-seq Barcoding based on Illumina sequencing of microsatellite loci. Both approaches revealed a clonal expansion event in 2013 and a more stable genetic structure during 2017-2020 compared to previous years. The identification of a mating type (MT) determination gene allowed to assign MT to strains isolated over the years. MTs were generally at equilibrium with two notable exceptions, including the clonal bloom of 2013. The populations exhibited linkage equilibrium in most blooms, indicating that sexual reproduction leads to genetic homogenization. Our findings show that P. multistriata blooms exhibit a dynamic genetic and demographic composition over time, most probably determined by deeper-layer cell inocula. Occasional clonal expansions and MT imbalances can potentially affect the persistence and ecological success of planktonic diatoms.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Rapid diversification of a free-living protist is driven by adaptation to climate and habitat

Microbial eukaryotes (protists) have major functional roles in aquatic ecosystems, including the biogeochemical cycling of elements as well as occupying various roles in the food web. Despite their importance for ecosystem function, the factors that drive diversification in protists are not known. Here, we aimed to identify the factors that drive differentiation and, subsequently, speciation in a free-living protist, Synura petersenii (Chrysophyceae). We sampled five different geographic areas and utilized population genomics and quantitative trait analyses. Habitat and climate were the major drivers of diversification on the local geographical scale, while geography played a role over longer distances. In addition to conductivity and temperature, precipitation was one of the most important environmental drivers of differentiation. Our results imply that flushing episodes (floods) drive microalgal adaptation to different niches, highlighting the potential for rapid diversification in protists.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Rapid diversification underlying the global dominance of a cosmopolitan phytoplankton

Marine phytoplankton play important roles in the global ecosystem, with a limited number of cosmopolitan keystone species driving their biomass. Recent studies have revealed that many of these phytoplankton are complexes composed of sibling species, but little is known about the evolutionary processes underlying their formation. Gephyrocapsa huxleyi, a widely distributed and abundant unicellular marine planktonic algae, produces calcified scales (coccoliths), thereby significantly affects global biogeochemical cycles via sequestration of inorganic carbon. This species is composed of morphotypes defined by differing degrees of coccolith calcification, the evolutionary ecology of which remains unclear. Here, we report an integrated morphological, ecological and genomic survey across globally distributed G. huxleyi strains to reconstruct evolutionary relationships between morphotypes in relation to their habitats. While G. huxleyi has been considered a single cosmopolitan species, our analyses demonstrate that it has evolved to comprise at least three distinct species, which led us to formally revise the taxonomy of the G. huxleyi complex. Moreover, the first speciation event occurred before the onset of the last interglacial period (~140 ka), while the second followed during this interglacial. Then, further rapid diversifications occurred during the most recent ice-sheet expansion of the last glacial period and established morphotypes as dominant populations across environmental clines. These results suggest that glacial-cycle dynamics contributed to the isolation of ocean basins and the segregations of oceans fronts as extrinsic drivers of micro-evolutionary radiations in extant marine phytoplankton.

Abiotic and biotic factors controlling sexual reproduction in populations of Pseudo-nitzschia pungens (Bacillariophyceae)

Pseudo-nitzschia pungens is a widely distributed marine pennate diatom. Hybrid zones, regions in which two different genotypes may interbreed, are important areas for speciation and ecology, and have been reported across the globe for this species. However, sexual reproduction between differing clades in the natural environment is yet to be observed and is difficult to predict. Here we carried out experiments using two mono-clonal cultures of P. pungens from different genotypes to measure the frequency and timing of sexual reproduction across varying biotic (growth phases and cell activity potential) and abiotic conditions (nutrients, light, turbulence). We found the mating rates and number of zygotes gradually decreased from exponential to late stationary growth phases. The maximum zygote abundance observed was 1,390 cells mL^-1 and the maximum mating rate was 7.1%, both which occurred during the exponential growth phase. Conversely, only 9 cells mL^-1 and a maximum mating rate of 0.1% was observed during the late stationary phase. We also found the higher the relative potential cell activity (rPCA) in parent cells, as determined by the concentration of chlorophyll a per cell and the ratio of colony formation during parent cultivations, revealed higher mating rates. Furthermore, sexual events were reduced under nutrient enrichment conditions, and mating pairs and zygotes were not formed under aphotic (dark) or shaking culture conditions (150 rpm). In order to understand the sexual reproduction of Pseudo-nitzschia in the natural environment, our results highlight that it is most likely the combination of both biotic (growth phase, Chl. a content) and abiotic factors (nutrients, light, turbulence) that will determine the successful union of intraspecific populations of P. pungens in any given region.

Biogeographical distribution of Hartaetosiga strains (Choanoflagellatea, Craspedida, Salpingoecidae) including morphological and transcriptomic data from a transect across the Atlantic Ocean

The genus Hartaetosiga Carr, Richter and Nitsche, 2017 comprised up to now only three species, H. gracilis (Kent) Carr, Richter, Nitsche, 2017, H. balthica (Wylezich and Karpov) Carr, Richter and Nitsche, 2017 and H. minima (Wylezich and Karpov) Carr, Richter and Nitsche, 2017. Based on distinct molecular data these species were relocated from the strictly freshwater genus Codosiga (Ehrenberg) Bütschli, 1878 to a new genus comprising brackish and marine species. During the cruise MSM82/2 across the Atlantic Ocean in 2019, surface water samples were taken from 15 stations along a transect ranging from 35°S to 23°N. We were able to isolate and cultivate 14 strains of the genus Hartaetosiga. Morphometric data showed no distinct morphological traits allowing for a species delineation, indicating a cryptic species complex within the genus. Based on cultivation, morphological data, and molecular analyses, we recorded H. gracilis for the first time from off-shelf waters of the Atlantic Ocean and could describe a new species, H. australis n. sp. This new species was recorded from sampling stations in the Southern Hemisphere only, which may indicate a potential biogeographic restriction likely caused by the Equatorial Counter Current (ECC), dividing the northern and southern surface waters.

Environmental filtering and mass effect are two important processes driving lake benthic diatoms: Results of a DNA metabarcoding study in a large lake

Environmental filtering is often found to dominate assembly rules in diatoms. These microalgae are diverse, especially at subspecies level, and tend to exhibit a niche phylogenetic conservatism. Therefore, other rules, such as competition or mass effects, should be detectable when environmental gradients and dispersal barriers are limited. We used metabarcoding to analyse benthic littoral diatom communities in 153 sites in a large lake (Geneva) exhibiting weak geographical barriers and weak environmental gradients outside river estuaries. We assessed assembly rules using variance partitioning, phylogenetic and source tracking analyses. No phylogenetic over-dispersion of communities, indicative of exclusive competition, was detected. Instead, we found these communities to be phylogenetically over-clustered, indicating environmental filtering, which was even stronger near river estuaries where environmental gradients are stronger. Finally, using a Bayesian method (SourceTracker), we found that rivers flowing into the lake bring communities that settle, especially in sites close to estuaries. Rivers with the highest discharges are primarily responsible for immigration, explaining 27% of lake composition. Therefore, despite favourable conditions to observe other rules, our results support that diatom communities are prominently assembled by environmental filtering and immigration processes, in particular from rivers. However, this does not exclude that other assembly rules may be at play at a finer spatial, temporal and/or phylogenetic scale.

Assessing the relative importance of stressors to the benthic index, M-AMBI: An example from U.S. estuaries

M-AMBI, a multivariate benthic index, has been used by European and American (U.S.) authorities to assess estuarine and coastal health and has been used in scientific studies throughout the world. It has been shown to be related to multiple pressures and stressors, but the relative importance of individual stressors within a multiple stressor context has not generally been assessed. In this study, we assembled data collected between 1999 and 2015 by the U.S. Environmental Protection Agency using consistent methods. These data included sediment and water quality measures and benthic invertebrate data which were used to calculate M-AMBI. We further assembled watersheds for all US estuaries with benthic data and calculated land use metrics. Random forest (RF) was used to identify those variables most strongly related to M-AMBI. Because RF is a compilation of multiple, nonlinear models, we then assessed which of these variables had a direct relationship with M-AMBI. The resulting variables were then assessed using RF to identify the subsets of variables that produced an effective and parsimonious model. This process was conducted at the national and ecoregional scale and the variables identified as being most important to predict M-AMBI were compared with literature reports of ecological patterns in a given area. At the national scale, better condition was correlated with clearer waters, lower amounts of agriculture in the watershed, and lower carbon and metal concentrations in estuarine sediments. Other stressors were identified as being important at the ecoregional scale, although sediment metal concentrations and watershed agriculture were identified as being important in most ecoregions. Our results suggest that this technique is useful to identify the most important variables impacting M-AMBI at broad spatial scales, even when the percentage of sites in Bad or Poor condition is low. This technique also provides an initial identification of important stressors that can be used to target more intensive local studies.

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

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.

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.

Reduced representation sequencing accurately quantifies relative abundance and reveals population-level variation in Pseudo-nitzschia spp

Certain species within the genus Pseudo-nitzschia are able to produce the neurotoxin domoic acid (DA), which can cause illness in humans, mass-mortality of marine animals, and closure of commercial and recreational shellfisheries during toxic events. Understanding and forecasting blooms of these harmful species is a primary management goal. However, accurately predicting the onset and severity of bloom events remains difficult, in part because the underlying drivers of bloom formation have not been fully resolved. Furthermore, Pseudo-nitzschia species often co-occur, and recent work suggests that the genetic composition of a Pseudo-nitzschia bloom may be a better predictor of toxicity than prevailing environmental conditions. We developed a novel next-generation sequencing assay using restriction site-associated DNA (2b-RAD) genotyping and applied it to mock Pseudo-nitzschia communities generated by mixing cultures of different species in known abundances. On average, 94% of the variance in observed species abundance was explained by the expected abundance. In addition, the false positive rate was low (0.45% on average) and unrelated to read depth, and false negatives were never observed. Application of this method to environmental DNA samples collected during natural Pseudo-nitzschia spp. bloom events in Southern California revealed that increases in DA were associated with increases in the relative abundance of P. australis. Although the absolute correlation across time-points was weak, an independent species fingerprinting assay (Automated Ribosomal Intergenic Spacer Analysis) supported this and identified other potentially toxic species. Finally, we assessed population-level genomic variation by mining SNPs from the environmental 2bRAD dataset. Consistent shifts in allele frequencies in P. pungens and P. subpacifica were detected between high and low DA years, suggesting that different intraspecific variants may be associated with prevailing environmental conditions or the presence of DA. Taken together, this method presents a potentially cost-effective and high-throughput approach for studies aiming to evaluate both population and species dynamics in mixed samples.

Phenological segregation suggests speciation by time in the planktonic diatom Pseudo‐nitzschia allochrona sp. nov

The processes leading to the emergence of new species are poorly understood in marine plankton, where weak physical barriers and homogeneous environmental conditions limit spatial and ecological segregation. Here, we combine molecular and ecological information from a long‐term time series and propose Pseudo‐nitzschia allochrona, a new cryptic planktonic diatom, as a possible case of speciation by temporal segregation. The new species differs in several genetic markers (18S, 28S and ITS rDNA fragments and rbcL) from its closest relatives, which are morphologically very similar or identical, and is reproductively isolated from its sibling species P. arenysensis. Data from a long‐term plankton time series show P. allochrona invariably occurring in summer–autumn in the Gulf of Naples, where its closely related species P. arenysensis, P. delicatissima, and P. dolorosa are instead found in winter–spring. Temperature and nutrients are the main factors associated with the occurrence of P. allochrona, which could have evolved in sympatry by switching its phenology and occupying a new ecological niche. This case of possible speciation by time shows the relevance of combining ecological time series with molecular information to shed light on the eco‐evolutionary dynamics of marine microorganisms.

The biogeographic differentiation of algal microbiomes in the upper ocean from pole to pole

Eukaryotic phytoplankton are responsible for at least 20% of annual global carbon fixation. Their diversity and activity are shaped by interactions with prokaryotes as part of complex microbiomes. Although differences in their local species diversity have been estimated, we still have a limited understanding of environmental conditions responsible for compositional differences between local species communities on a large scale from pole to pole. Here, we show, based on pole-to-pole phytoplankton metatranscriptomes and microbial rDNA sequencing, that environmental differences between polar and non-polar upper oceans most strongly impact the large-scale spatial pattern of biodiversity and gene activity in algal microbiomes. The geographic differentiation of co-occurring microbes in algal microbiomes can be well explained by the latitudinal temperature gradient and associated break points in their beta diversity, with an average breakpoint at 14 °C ± 4.3, separating cold and warm upper oceans. As global warming impacts upper ocean temperatures, we project that break points of beta diversity move markedly pole-wards. Hence, abrupt regime shifts in algal microbiomes could be caused by anthropogenic climate change.

Patterns of protist distribution and diversification in alpine lakes across Europe

Biogeography in Europe is known to be crucially influenced by the large mountain ranges serving as biogeographical islands for cold-adapted taxa and geographical barriers for warm-adapted taxa. While biogeographical patterns are well-known for plants and animals in Europe, we here investigated diversity and distribution patterns of protist freshwater communities on a European scale (256 lakes) in the light of the well-studied post-glacial distribution patterns of macroorganisms. Thus, our study compared 43 alpine protist communities of lakes located in the Alps, Carpathians, Pyrenees, and the Sierra Nevada with that of surrounding lowland lakes. We verified altitudinal diversity gradients of freshwater protists with decreasing richness and diversity across altitudes similar to those observed for plants and animals. Alpine specialists and generalists could be identified differing significantly in richness and diversity, but hardly in occurrence and proportions of major taxonomic groups. High proportions of region-specific alpine specialists indicate an increased occurrence of distinct lineages within each mountain range and thus, suggested either separated glacial refugia or post-glacial diversification within mountain ranges. However, a few alpine specialists were shared between mountain ranges suggesting a post-glacial recolonization from a common lowland pool. Our results identified generalists with wide distribution ranges and putatively wide tolerance ranges toward environmental conditions as main drivers of protist diversification (specification) in alpine lakes, while there was hardly any diversification in alpine specialists.

The importance of type species and their correct identification: A key example from tintinnid ciliates (Alveolata, Ciliophora, Spirotricha)

Types and the corresponding rules in the International Code of Zoological Nomenclature are crucial for taxonomy and are meant to provide nomenclatural stability. In the case of neotypification, especially diligent taxonomic work is required to retain continuity. In the present communication, we first outline the main principles of typification and neotypification. We then discuss a critical case, using a current example from the marine planktonic tintinnid genus Tintinnopsis Stein, 1867 (Alveolata, Ciliophora). This diverse and ubiquitous genus is nonmonophyletic, but its revision and the erection of new related genera is currently prevented by the uncertain affiliation of its type species.

New alleles in the mating type determination region of West Atlantic strains of Pseudo-nitzschia multistriata

The cosmopolitan, species-rich diatom genus Pseudo-nitzschia represents a good system for the study of speciation, evolution and diversity. Understanding elements linked to population dynamics and life cycle regulation for these species is of particular importance in view of their ability to produce the toxin domoic acid and cause harmful blooms. Pseudo-nitzschia multistriata, one of the toxic species that represents a model for the study of life cycle related questions, is the only diatom for which a sex determination mechanism has been described. Populations in the Gulf of Naples (Mediterranean Sea), can share four different allelic variants (A, M, B, N) of the mating type determination region, and one of them (A) is responsible for the determination of the mating type + (MT+), defined by the MT+ restricted expression of the gene MRP3. Here, we analysed the sex determination genomic region in three new strains isolated from the Gulf of Mexico and compared it to the alleles previously described in the Mediterranean strains. We first show that these geographically distant strains of P. multistriata belong to different populations but can interbreed. Next, we show that the two populations share an overall similar structure of the genomic locus although differences can be seen in the polymorphic regions upstream of MRP3. In strain P4-C1, we amplified and sequenced an allele (M) identical to one of those previously characterized in the Mediterranean strains. In the other two strains, P4-C2 and P4-C5, we identified three new alleles, which we named A2, B2 and N2. P4-C2 and P4-C5 are heterozygous and share the common allele A2 linked to the monoallelic expression of the MT+ specific sex determining gene MRP3. Our results expand information on the global distribution of P. multistriata and on the level of conservation of the sex determination region in different populations. The definition of the extent of intra- and inter-specific conservation of this region would be a relevant addition to our understanding of Pseudo-nitzschia diversity and evolution.

Oceanic dispersal barriers in a holoplanktonic gastropod

Pteropods, a group of holoplanktonic gastropods, are regarded as bioindicators of the effects of ocean acidification on open ocean ecosystems, because their thin aragonitic shells are susceptible to dissolution. While there have been recent efforts to address their capacity for physiological acclimation, it is also important to gain predictive understanding of their ability to adapt to future ocean conditions. However, little is known about the levels of genetic variation and large-scale population structuring of pteropods, key characteristics enabling local adaptation. We examined the spatial distribution of genetic diversity in the mitochondrial cytochrome c oxidase I (COI) and nuclear 28S gene fragments, as well as shell shape variation, across a latitudinal transect in the Atlantic Ocean (35°N-36°S) for the pteropod Limacina bulimoides. We observed high levels of genetic variability (COI π = 0.034, 28S π = 0.0021) and strong spatial structuring (COI ΦST  = 0.230, 28S ΦST  = 0.255) across this transect. Based on the congruence of mitochondrial and nuclear differentiation, as well as differences in shell shape, we identified a primary dispersal barrier in the southern Atlantic subtropical gyre (15-18°S). This barrier is maintained despite the presence of expatriates, a gyral current system, and in the absence of any distinct oceanographic gradients in this region, suggesting that reproductive isolation between these populations must be strong. A secondary dispersal barrier supported only by 28S pairwise ΦST comparisons was identified in the equatorial upwelling region (between 15°N and 4°S), which is concordant with barriers observed in other zooplankton species. Both oceanic dispersal barriers were congruent with regions of low abundance reported for a similar basin-scale transect that was sampled 2 years later. Our finding supports the hypothesis that low abundance indicates areas of suboptimal habitat that result in barriers to gene flow in widely distributed zooplankton species. Such species may in fact consist of several populations or (sub)species that are adapted to local environmental conditions, limiting their potential for adaptive responses to ocean changes. Future analyses of genome-wide diversity in pteropods could provide further insight into the strength, formation and maintenance of oceanic dispersal barriers.

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.

Morphological Variability of Pseudo-nitzschia pungens Clade I (Bacillariophyceae) in the Northwestern Adriatic Sea

Pseudo-nitzschia pungens is a common component of the phytoplankton community of the northern Adriatic Sea. In this study, an in-depth morphological analysis of P. pungens was carried out in both cultured strains isolated in different periods and field samples, revealing a surprisingly wide variability in a number of details, with both the gross morphology and ultrastructural levels deviating from the nominal P. pungens. Colonies showed an overlap (from one-third to one-sixth) and a transapical axis (rarely reaching 3 µm), strongly differing from the original description of the species. Moreover, valves may be either symmetrical or slightly asymmetrical, with striae almost always biseriate but sometimes uniseriate or triseriate. Poroids' morphology in cingular bands was characterized by a wide variability (square, circular, or rectangular poroids without or with up to two hymen sectors), with several combination of them, even within the same cingular band. Phylogenetic analyses based on ITS rDNA showed that the P. pungens of the northern Adriatic Sea belonged to clade I. Domoic acid was not detected.

Functional significance of phylogeographic structure in a toxic benthic marine microbial eukaryote over a latitudinal gradient along the East Australian Current

Genetic diversity in marine microbial eukaryotic populations (protists) drives their ecological success by enabling diverse phenotypes to respond rapidly to changing environmental conditions. Despite enormous population sizes and lack of barriers to gene flow, genetic differentiation that is associated with geographic distance, currents, and environmental gradients has been reported from planktonic protists. However, for benthic protists, which have reduced dispersal opportunities, phylogeography and its phenotypic significance are little known. In recent years, the East Australian Current (EAC) has intensified its southward flow, associated with the tropicalization of temperate waters. Benthic harmful algal species have been increasingly found in south-eastern Australia. Yet little is known about the potential of these species to adapt or extend their range in relation to changing conditions. Here, we examine genetic diversity and functional niche divergence in a toxic benthic dinoflagellate, Ostreopsis cf. siamensis, along a 1,500 km north-south gradient in southeastern Australia. Sixty-eight strains were established from eight sampling sites. The study revealed long-standing genetic diversity among strains established from the northern-most sites, along with large phenotypic variation in observed physiological traits such as growth rates, cell volume, production of palytoxin-like compounds, and photophysiological parameters. Strains from the southern populations were more uniform in both genetic and functional traits, and have possibly colonized their habitats more recently. Our study reports significant genetic and functional trait variability in a benthic harmful algal species, indicative of high adaptability, and a possible climate-driven range extension. The observed high trait variation may facilitate development of harmful algal blooms under dynamic coastal environmental conditions.

Sexual reproduction and genetic polymorphism within the cosmopolitan marine diatom Pseudo-nitzschia pungens

Different clades belonging to the cosmopolitan marine diatom Pseudo-nitzschia pungens appear to be present in different oceanic environments, however, a ‘hybrid zone’, where populations of different clades interbreed, has also been reported. Many studies have investigated the sexual reproduction of P. pungens, focused on morphology and life cycle, rather than the role of sexual reproduction in mixing the genomes of their parents. We carried out crossing experiments to determine the sexual compatibility/incompatibility between different clades of P. pungens, and examined the genetic polymorphism in the ITS2 region. Sexual reproduction did not occur only between clades II and III under any of experimental temperature conditions. Four offspring strains were established between clade I and III successfully. Strains established from offspring were found interbreed with other offspring strains as well as viable with their parental strains. We confirmed the hybrid sequence patterns between clades I and III and found novel sequence types including polymorphic single nucleotide polymorphisms (SNPs) in the offspring strains. Our results implicate that gene exchange and mixing between different clades are still possible, and that sexual reproduction is a significant ecological strategy to maintain the genetic diversity within this diatom species.

Global radiation in a rare biosphere soil diatom

Soil micro-organisms drive the global carbon and nutrient cycles that underlie essential ecosystem functions. Yet, we are only beginning to grasp the drivers of terrestrial microbial diversity and biogeography, which presents a substantial barrier to understanding community dynamics and ecosystem functioning. This is especially true for soil protists, which despite their functional significance have received comparatively less interest than their bacterial counterparts. Here, we investigate the diversification of Pinnularia borealis, a rare biosphere soil diatom species complex, using a global sampling of >800 strains. We document unprecedented high levels of species-diversity, reflecting a global radiation since the Eocene/Oligocene global cooling. Our analyses suggest diversification was largely driven by colonization of novel geographic areas and subsequent evolution in isolation. These results illuminate our understanding of how protist diversity, biogeographical patterns, and members of the rare biosphere are generated, and suggest allopatric speciation to be a powerful mechanism for diversification of micro-organisms.

It is generally thought many microbes, owing to their ubiquity and dispersal capability, lack biogeographic structuring and clear speciation patterns compared to macroorganisms. However, Pinseel et al. demonstrate multiple cycles of colonization and diversification in Pinnularia borealis, a rare biosphere soil diatom.

Disentangling the mechanisms shaping the surface ocean microbiota

BACKGROUND

The ocean microbiota modulates global biogeochemical cycles and changes in its configuration may have large-scale consequences. Yet, the underlying ecological mechanisms structuring it are unclear. Here, we investigate how fundamental ecological mechanisms (selection, dispersal and ecological drift) shape the smallest members of the tropical and subtropical surface-ocean microbiota: prokaryotes and minute eukaryotes (picoeukaryotes). Furthermore, we investigate the agents exerting abiotic selection on this assemblage as well as the spatial patterns emerging from the action of ecological mechanisms. To explore this, we analysed the composition of surface-ocean prokaryotic and picoeukaryotic communities using DNA-sequence data (16S- and 18S-rRNA genes) collected during the circumglobal expeditions Malaspina-2010 and TARA-Oceans.

RESULTS

We found that the two main components of the tropical and subtropical surface-ocean microbiota, prokaryotes and picoeukaryotes, appear to be structured by different ecological mechanisms. Picoeukaryotic communities were predominantly structured by dispersal-limitation, while prokaryotic counterparts appeared to be shaped by the combined action of dispersal-limitation, selection and drift. Temperature-driven selection appeared as a major factor, out of a few selected factors, influencing species co-occurrence networks in prokaryotes but not in picoeukaryotes, indicating that association patterns may contribute to understand ocean microbiota structure and response to selection. Other measured abiotic variables seemed to have limited selective effects on community structure in the tropical and subtropical ocean. Picoeukaryotes displayed a higher spatial differentiation between communities and a higher distance decay when compared to prokaryotes, consistent with a scenario of higher dispersal limitation in the former after considering environmental heterogeneity. Lastly, random dynamics or drift seemed to have a more important role in structuring prokaryotic communities than picoeukaryotic counterparts.

CONCLUSIONS

The differential action of ecological mechanisms seems to cause contrasting biogeography, in the tropical and subtropical ocean, among the smallest surface plankton, prokaryotes and picoeukaryotes. This suggests that the idiosyncrasy of the main constituents of the ocean microbiota should be considered in order to understand its current and future configuration, which is especially relevant in a context of global change, where the reaction of surface ocean plankton to temperature increase is still unclear. Video Abstract.

Ecological time series and integrative taxonomy unveil seasonality and diversity of the toxic diatom Pseudo-nitzschia H. Peragallo in the northern Adriatic Sea

Pseudo-nitzschia H. Peragallo (1900) is a globally distributed genus of pennate diatoms that are important components of phytoplankton communities worldwide. Some members of the genus produce the neurotoxin domoic acid, so regular monitoring is in place. However, the identification of toxic members in routine samplings remains problematic. In this study, the diversity and seasonal occurrence of Pseudo-nitzschia species were investigated in the Gulf of Trieste, a shallow gulf in the northern Adriatic Sea. We used time series data from 2005 to 2018 to describe the seasonal and inter-annual occurrence of the genus in the area and its contribution to the phytoplankton community. On average, the genus accounted for about 15 % of total diatom abundance and peaked in spring and autumn, with occasional outbreaks during summer and large inter-annual fluctuations. Increased water temperature and decreased salinity positively affected the presence of some members of the genus, while strong effects could be masked by an unsuitable definition of the species complexes used for monitoring purposes. Therefore, combining morphological (TEM) and molecular analyses by sequencing the ITS, 28S and rbcL markers, eight species were identified from 83 isolated monoclonal strains: P. calliantha, P. fraudulenta, P. delicatissima, P. galaxiae, P. mannii, P. multistriata, P. pungens and P. subfraudulenta. A genetic comparison between the isolated strains and other strains in the Mediterranean was carried out and rbcL was inspected as a potential barcode marker in respect to our results. This is the first study in the Gulf of Trieste on Pseudo-nitzschia time series from a long-term ecological research (LTER) site coupled with molecular data. We show that meaningful ecological conclusions can be drawn by applying integrative methodology, as opposed to the approach that only considers species complexes. The results of this work will provide guidance for further monitoring efforts as well as research activities, including population genetics and genomics, associated with seasonal distribution and toxicity profiles.

Acclimation of the Marine Diatom Pseudo-nitzschia australis to Different Salinity Conditions: Effects on Growth, Photosynthetic Activity, and Domoic Acid Content1

Toxic Pseudo-nitzschia australis strains isolated from French coastal waters were studied to investigate their capacity to adapt to different salinities. Their acclimation to different salinity conditions (10, 20, 30, 35, and 40) was studied on growth, photosynthetic capacity, cell biovolume, and domoic acid (DA) content. The strains showed an ability to acclimate to a salinity range from 20 to 40, with optimal growth rates between salinities 30 and 40. The highest cell biovolume was observed at the lowest salinity 20 and was associated with the lowest growth rate. Salinity did not affect the photosynthetic activity; F_v /F_m values and the pigment contents remained high with no significant difference among salinities. An enhanced production of zeaxanthin was, however, observed in the late stationary and decline phases in all cultures except for those acclimated to salinity 20. In terms of cellular toxin content, DA concentrations were 2 to 3-fold higher at the lowest salinity (20) than at the other salinities and were combined with a low amount of dissolved DA. The fact that P. australis accumulate more DA per cell in less saline waters, illustrates that climate-related changes in salinity may affect Pseudo-nitzschia physiology through direct effects on growth, physiology, and toxin content.

Genetic diversity and novel lineages in the cosmopolitan copepod Pleuromamma abdominalis in the Southeast Pacific

Across boundary currents, zooplankton are subject to strong oceanographic gradients and hence strong selective pressures. How such gradients interact with the speciation process of pelagic organisms is still poorly understood in the open ocean realm. Here we report on genetic diversity within the pelagic copepod Pleuromamma abdominalis in the poorly known Southeast Pacific region, with samples spanning an ocean gradient from coastal upwelling to the oligotrophic South Pacific Subtropical Gyre. We assessed variation in fragments of the mitochondrial (mt) genes cytochrome c oxidase subunit I (COI) and Cytochrome b as well as in the nuclear internal transcribed spacer (ITS) region and 28 S rRNA. Phylogenetic analyses revealed the presence of 8 divergent lineages occurring across the gradient with genetic distances in the range of 0.036 and 0.44 (mt genes), and GMYC species delimitation methods support their inference as distinct (undescribed) species. Genetic lineages occurring across the zonal gradient showed strong genetic structuring, with the presence of at least two new lineages within the coastal upwelling zone, revealing an unexpectedly high level of endemism within the Humboldt Current System. Multivariate analyses found strong correlation between genetic variation and surface chlorophyll-a and salinity, suggesting an important role for hydrographic gradients in maintaining genetic diversity. However, the presence of cryptic lineages within the upwelling zone cannot be easily accounted for by environmental heterogeneity and poses challenging questions for understanding the speciation process for oceanic zooplankton.

Speciation in protists: Spatial and ecological divergence processes cause rapid species diversification in a freshwater chrysophyte

Although eukaryotic microorganisms are extremely numerous, diverse and essential to global ecosystem functioning, they are largely understudied by evolutionary biologists compared to multicellular macroscopic organisms. In particular, very little is known about the speciation mechanisms which may give rise to the diversity of microscopic eukaryotes. It was postulated that the enormous population sizes and ubiquitous distribution of these organisms could lead to a lack of population differentiation and therefore very low speciation rates. However, such assumptions have traditionally been based on morphospecies, which may not accurately reflect the true diversity, missing cryptic taxa. In this study, we aim to articulate the major diversification mechanisms leading to the contemporary molecular diversity by using a colonial freshwater flagellate, Synura sphagnicola, as an example. Phylogenetic analysis of five sequenced loci showed that S. sphagnicola differentiated into two morphologically distinct lineages approximately 15.4 million years ago, which further diverged into several evolutionarily recent haplotypes during the late Pleistocene. The most recent haplotypes are ecologically and biogeographically much more differentiated than the old lineages, presumably because of their persistent differentiation after the allopatric speciation events. Our study shows that in microbial eukaryotes, species diversification via the colonization of new geographical regions or ecological resources occurs much more readily than was previously thought. Consequently, divergence times of microorganisms in some lineages may be equivalent to the estimated times of speciation in plants and animals.

Multilocus Sequence Typing (MLST) and Random Polymorphic DNA (RAPD) Comparisons of Geographic Isolates of Neoparamoeba perurans, the Causative Agent of Amoebic Gill Disease

Neoparamoba perurans, is the aetiological agent of amoebic gill disease (AGD), a disease that affects farmed Atlantic salmon worldwide. Multilocus sequence typing (MLST) and Random Amplified Polymorphic DNA (RAPD) are PCR-based typing methods that allow for the highly reproducible genetic analysis of population structure within microbial species. To the best of our knowledge, this study represents the first use of these typing methods applied to N. perurans with the objective of distinguishing geographical isolates. These analyses were applied to a total of 16 isolates from Australia, Canada, Ireland, Scotland, Norway, and the USA. All the samples from Australia came from farm sites on the island state of Tasmania. Genetic polymorphism among isolates was more evident from the RAPD analysis compared to the MLST that used conserved housekeeping genes. Both techniques consistently identified that isolates of N. perurans from Tasmania, Australia were more similar to each other than to the isolates from other countries. While genetic differences were identified between geographical isolates, a BURST analysis provided no evidence of a founder genotype. This suggests that emerging outbreaks of AGD are not due to rapid translocation of this important salmonid pathogen from the same area.

Contributions of single-cell genomics to our understanding of planktonic marine archaea

Single-cell genomics has transformed many fields of biology, marine microbiology included. Here, we consider the impact of single-cell genomics on a specific group of marine microbes-the planktonic marine archaea. Despite single-cell enabled discoveries of novel metabolic function in the marine thaumarchaea, population-level investigations are hindered by an overall lower than expected recovery of thaumarchaea in single-cell studies. Metagenome-assembled genomes have so far been a more useful method for accessing genome-resolved insights into the Marine Group II euryarchaea. Future progress in the application of single-cell genomics to archaeal biology in the ocean would benefit from more targeted sorting approaches, and a more systematic investigation of potential biases against archaea in single-cell workflows including cell lysis, genome amplification and genome screening. This article is part of a discussion meeting issue 'Single cell ecology'.

Repeated species radiations in the recent evolution of the key marine phytoplankton lineage Gephyrocapsa

Phytoplankton account for nearly half of global primary productivity and strongly affect the global carbon cycle, yet little is known about the forces that drive the evolution of these keystone microscopic organisms. Here we combine morphometric data from the fossil record of the ubiquitous coccolithophore genus Gephyrocapsa with genomic analyses of extant species to assess the genetic processes underlying Pleistocene palaeontological patterns. We demonstrate that all modern diversity in Gephyrocapsa (including Emiliania huxleyi) originated in a rapid species radiation during the last 0.6 Ma, coincident with the latest of the three pulses of Gephyrocapsa diversification and extinction documented in the fossil record. Our evolutionary genetic analyses indicate that new species in this genus have formed in sympatry or parapatry, with occasional hybridisation between species. This sheds light on the mode of speciation during evolutionary radiation of marine phytoplankton and provides a model of how new plankton species form.

The phytoplankton Gephyrocapsa have gone through repeated macroevolutionary shifts in size. Here, Bendif et al. combine fossil and genomic data to show the latest shift was coincident with a species radiation and suggest that previous shifts have also resulted from cycles of radiation and extinction.

Thermal Niche Differentiation in the Benthic Diatom Cylindrotheca closterium (Bacillariophyceae) Complex

Coastal waters are expected to undergo severe warming in the coming decades. Very little is known about how diatoms, the dominant primary producers in these habitats, will cope with these changes. We investigated the thermal niche of Cylindrotheca closterium, a widespread benthic marine diatom, using 24 strains collected over a wide latitudinal gradient. A multi-marker phylogeny in combination with a species delimitation approach shows that C. closterium represents a (pseudo)cryptic species complex, and this is reflected in distinct growth response patterns in terms of optimum growth temperature, maximum growth rate, and thermal niche width. Strains from the same clade displayed a similar thermal response, suggesting niche conservation between closely related strains. Due to their lower maximum growth rate and smaller thermal niche width, we expect the polar species to be particularly sensitive to warming, and, in the absence of adaptation, to be replaced with species from lower latitudes.

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.

Annotated 18S and 28S rDNA reference sequences of taxa in the planktonic diatom family Chaetocerotaceae

The species-rich diatom family Chaetocerotaceae is common in the coastal marine phytoplankton worldwide where it is responsible for a substantial part of the primary production. Despite its relevance for the global cycling of carbon and silica, many species are still described only morphologically, and numerous specimens do not fit any described taxa. Nowadays, studies to assess plankton biodiversity deploy high throughput sequencing metabarcoding of the 18S rDNA V4 region, but to translate the gathered metabarcodes into biologically meaningful taxa, there is a need for reference barcodes. However, 18S reference barcodes for this important family are still relatively scarce. We provide 18S rDNA and partial 28S rDNA reference sequences of 443 morphologically characterized chaetocerotacean strains. We gathered 164 of the 216 18S sequences and 244 of the 413 28S sequences of strains from the Gulf of Naples, Atlantic France, and Chile. Inferred phylogenies showed 84 terminal taxa in seven principal clades. Two of these clades included terminal taxa whose rDNA sequences contained spliceosomal and Group IC1 introns. Regarding the commonly used metabarcode markers in planktonic diversity studies, all terminal taxa can be discriminated with the 18S V4 hypervariable region; its primers fit their targets in all but two species, and the V4-tree topology is similar to that of the 18S. Hence V4-metabarcodes of unknown Chaetocerotaceae are assignable to the family. Regarding the V9 hypervariable region, most terminal taxa can be discriminated, but several contain introns in their primer targets. Moreover, poor phylogenetic resolution of the V9 region affects placement of metabarcodes of putative but unknown chaetocerotacean taxa, and hence, uncertainty in taxonomic assignment, even of higher taxa.

Novel Widespread Marine Oomycetes Parasitising Diatoms, Including the Toxic Genus Pseudo-nitzschia: Genetic, Morphological, and Ecological Characterisation

Parasites are key drivers of phytoplankton bloom dynamics and related aquatic ecosystem processes. Yet, the dearth of morphological and molecular information hinders the assessment of their diversity and ecological role. Using single-cell techniques, we characterise morphologically and molecularly, intracellular parasitoids infecting four potentially toxin-producing Pseudo-nitzschia and one Melosira species on the North Atlantic coast. These sequences define two, morphologically indistinguishable clades within the phylum Oomycota, related to the genera of algal parasites Anisolpidium and Olpidiopsis and the diatom parasitoid species Miracula helgolandica. Our morphological data are insufficient to attribute either clade to the still unsequenced genus Ectrogella; hence it is proposed to name the clades OOM_1 and OOM_2. A screening of global databases of the barcode regions V4 and V9 of the 18S rDNA demonstrate the presence of these parasitoids beyond the North Atlantic coastal region. During a biweekly metabarcoding survey (Concarneau Bay, France), reads associated with one sequenced parasitoid coincided with the decline of Cerataulina pelagica bloom, whilst the other parasitoids co-occurred at low abundance with Pseudo-nitzschia. Our data highlight a complex and unexplored diversity of the oomycete parasitoids of diatoms and calls for the investigation of their phenology, evolution, and potential contribution in controlling their host spatial-temporal dynamics.

Pseudo-nitzschia, Nitzschia, and domoic acid: New research since 2011

Some diatoms of the genera Pseudo-nitzschia and Nitzschia produce the neurotoxin domoic acid (DA), a compound that caused amnesic shellfish poisoning (ASP) in humans just over 30 years ago (December 1987) in eastern Canada. This review covers new information since two previous reviews in 2012. Nitzschia bizertensis was subsequently discovered to be toxigenic in Tunisian waters. The known distribution of N. navis-varingica has expanded from Vietnam to Malaysia, Indonesia, the Philippines and Australia. Furthermore, 15 new species (and one new variety) of Pseudo-nitzschia have been discovered, bringing the total to 52. Seven new species were found to produce DA, bringing the total of toxigenic species to 26. We list all Pseudo-nitzschia species, their ability to produce DA, and show their global distribution. A consequence of the extended distribution and increased number of toxigenic species worldwide is that DA is now found more pervasively in the food web, contaminating new marine organisms (especially marine mammals), affecting their physiology and disrupting ecosystems. Recent findings highlight how zooplankton grazers can induce DA production in Pseudo-nitzschia and how bacteria interact with Pseudo-nitzschia. Since 2012, new discoveries have been reported on physiological controls of Pseudo-nitzschia growth and DA production, its sexual reproduction, and infection by an oomycete parasitoid. Many advances are the result of applying molecular approaches to discovering new species, and to understanding the population genetic structure of Pseudo-nitzschia and mechanisms used to cope with iron limitation. The availability of genomes from three Pseudo-nitzschia species, coupled with a comparative transcriptomic approach, has allowed advances in our understanding of the sexual reproduction of Pseudo-nitzschia, its signaling pathways, its interactions with bacteria, and genes involved in iron and vitamin B₁₂ and B₇ metabolism. Although there have been no new confirmed cases of ASP since 1987 because of monitoring efforts, new blooms have occurred. A massive toxic Pseudo-nitzschia bloom affected the entire west coast of North America during 2015-2016, and was linked to a 'warm blob' of ocean water. Other smaller toxic blooms occurred in the Gulf of Mexico and east coast of North America. Knowledge gaps remain, including how and why DA and its isomers are produced, the world distribution of potentially toxigenic Nitzschia species, the prevalence of DA isomers, and molecular markers to discriminate between toxigenic and non-toxigenic species and to discover sexually reproducing populations in the field.

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.

Sexual ancestors generated an obligate asexual and globally dispersed clone within the model diatom species Thalassiosira pseudonana

Sexual reproduction roots the eukaryotic tree of life, although its loss occurs across diverse taxa. Asexual reproduction and clonal lineages persist in these taxa despite theoretical arguments suggesting that individual clones should be evolutionarily short-lived due to limited phenotypic diversity. Here, we present quantitative evidence that an obligate asexual lineage emerged from a sexual population of the marine diatom Thalassiosira pseudonana and rapidly expanded throughout the world’s oceans. Whole genome comparisons identified two lineages with characteristics expected of sexually reproducing strains in Hardy-Weinberg equilibrium. A third lineage displays genomic signatures for the functional loss of sexual reproduction followed by a recent global colonization by a single ancestral genotype. Extant members of this lineage are genetically differentiated and phenotypically plastic, potentially allowing for rapid adaptation when they are challenged by natural selection. Such mechanisms may be expected to generate new clones within marginal populations of additional unicellular species, facilitating the exploration and colonization of novel environments, aided by exponential growth and ease of dispersal.

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'.

A planktonic diatom displays genetic structure over small spatial scales

Marine planktonic microalgae have potentially global dispersal, yet reduced gene flow has been confirmed repeatedly for several species. Over larger distances (>200 km) geographic isolation and restricted oceanographic connectivity have been recognized as instrumental in driving population divergence. Here we investigated whether similar patterns, that is, structured populations governed by geographic isolation and/or oceanographic connectivity, can be observed at smaller (6-152 km) geographic scales. To test this we established 425 clonal cultures of the planktonic diatom Skeletonema marinoi collected from 11 locations in the Archipelago Sea (northern Baltic Sea). The region is characterized by a complex topography, entailing several mixing regions of which four were included in the sampling area. Using eight microsatellite markers and conventional F-statistics, significant genetic differentiation was observed between several sites. Moreover, Bayesian cluster analysis revealed the co-occurrence of two genetic groups spread throughout the area. However, geographic isolation and oceanographic connectivity could not explain the genetic patterns observed. Our data reveal hierarchical genetic structuring whereby despite high dispersal potential, significantly diverged populations have developed over small spatial scales. Our results suggest that biological characteristics and historical events may be more important in generating barriers to gene flow than physical barriers at small spatial scales.

Revealing the distinct habitat ranges and hybrid zone of genetic sub-populations within Pseudo-nitzschia pungens (Bacillariophyceae) in the West Pacific area

Genetic sub-populations (clades) of cosmopolitan marine diatom Pseudo-nitzschia pungens might have distinct habitats, and their hybrid zone is suspected in higher latitude area of the West Pacific area, however, it is still unrevealed because of technical difficulties and lack of evidences in natural environments. The aim of this study is to investigate the habitat characteristics of each clade of P. pungens on geographical distribution with the habitat temperature ranges of each clade and to reveal their hybrid zone in the West Pacific area. We employed the 137 number of nucleotide sequences of P. pungens and its sampling data (spatial and temporal scale) originated from the West Pacific area, and used field application of qPCR assay for intra-specific level of P. pungens. Only two genotypes, clade I and III, were identified in the West Pacific area. Clade I was distributed from 39 to 32.3°N, and clade III were from 1.4 to 34.4°N. The estimated habitat temperature for the clade I and clade III ranges were 8.1-26.9 °C and 24.2-31.2 °C, respectively. The latitudinal distributions and temperature ranges of each clade were significantly different. The qPCR assay employed, and results suggested that the hybrid zone for clade I and III has been observed in the southern Korean coasts, and clade III might be introduced from the Southern Pacific area. The cell abundances of clade III were strongly related with the higher seawater temperature and warm current force. This study has defined distinct habitat characteristics of genetically different sub-populations of P. pungens, and revealed its hybrid zone in natural environment for the first time. We also provided strong evidences about dispersion of the population of clade III to higher latitude in the West Pacific area.

Sandwich hybridization probes for the detection of Pseudo-nitzschia (Bacillariophyceae) species: An update to existing probes and a description of new probes

New sandwich hybridization assay (SHA) probes for detecting Pseudo-nitzschia species (P. arenysensis, P. fraudulenta, P. hasleana, P. pungens) are presented, along with updated cross-reactivity information on historical probes (SHA and FISH; fluorescence in situ hybridization) targeting P. australis and P. multiseries. Pseudo-nitzschia species are a cosmopolitan group of diatoms that produce varying levels of domoic acid (DA), a neurotoxin that can accumulate in finfish and shellfish and transfer throughout the food web. Consumption of infected food sources can lead to illness in humans (amnesic shellfish poisoning; ASP) and marine wildlife (domoic acid poisoning; DAP). The threat of human illness, along with economic loss from fishery closures has resulted in the implementation of monitoring protocols and intensive ecological studies. SHA probes have been instrumental in some of these efforts, as the technique performs well in complex heterogeneous sample matrices and has been adapted to benchtop and deployable (Environmental Sample Processor) platforms. The expanded probe set will enhance future efforts towards understanding spatial, temporal and successional patterns in species during bloom and non-bloom periods.