SAG-RAD: A Method for Single-Cell Population Genomics of Unicellular Eukaryotes.

Sequencing of reduced representation libraries enables genotyping of many individuals for population genomic studies. However, high amounts of DNA are required, and the method cannot be applied directly on single cells, preventing its use on most microbes. We developed and implemented the analysis of single amplified genomes followed by restriction-site-associated DNA sequencing to bypass labor-intensive culturing and to avoid culturing bias in population genomic studies of unicellular eukaryotes. This method thus opens the way for addressing important questions about the genetic diversity, gene flow, adaptation, dispersal, and biogeography of hitherto unexplored species.

Strain-specific metabarcoding reveals rapid evolution of copper tolerance in populations of the coastal diatom Skeletonema marinoi.

Phytoplankton have short generation times, flexible reproduction strategies, large population sizes and high standing genetic diversity, traits that should facilitate rapid evolution under directional selection. We quantified local adaptation of copper tolerance in a population of the diatom Skeletonema marinoi from a mining-exposed inlet in the Baltic Sea and in a non-exposed population 100 km away. We hypothesized that mining pollution has driven evolution of elevated copper tolerance in the impacted population of S. marinoi. Assays of 58 strains originating from sediment resting stages revealed no difference in the average tolerance to copper between the two populations. However, variation within populations was greater at the mining site, with three strains displaying hyper-tolerant phenotypes. In an artificial evolution experiment, we used a novel intraspecific metabarcoding locus to track selection and quantify fitness of all 58 strains during co-cultivation in one control and one toxic copper treatment. As expected, the hyper-tolerant strains enabled rapid evolution of copper tolerance in the mining-exposed population through selection on available strain diversity. Within 42 days, in each experimental replicate a single strain dominated (30%-99% abundance) but different strains dominated the different treatments. The reference population developed tolerance beyond expectations primarily due to slowly developing plastic response in one strain, suggesting that different modes of copper tolerance are present in the two populations. Our findings provide novel empirical evidence that standing genetic diversity of phytoplankton resting stage allows populations to evolve rapidly (20-50 generations) and flexibly on timescales relevant for seasonal bloom progressions.

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

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

Genome-wide single nucleotide polymorphism markers reveal population structure and dispersal direction of an expanding nuisance algal bloom species.

Species invasion and range expansion are currently under scrutiny due to increasing anthropogenic impact on the natural environment. This is also true for harmful algal blooms, which have been reported to have increased in frequency. However, this research is challenging due to the ephemeral nature, small size and mostly low concentrations of microalgae in the environment. One such species is the nuisance microalga Gonyostomum semen (Raphidophyceae), which has increased in occurrence in northern Europe in recent decades. The question of whether the species has expanded its habitat range or if it was already present in the lakes but was too rare to be detected remains unanswered. The aim of the present study was to determine the genetic structure and dispersal pathways of G. semen using RAD (restriction-site-associated DNA) tag sequencing. For G. semen, which has a huge genome (32 Gbp), we faced particular challenges, but were nevertheless able to recover over 1000 single nucleotide polymorphisms at high coverage. Our data revealed a distinct population genetic structure, demonstrating a divide of western and eastern populations that probably represent different lineages. Despite significant genetic differentiation among lakes, we found only limited isolation-by-distance. While we had expected a pattern of recent expansion northwards, the data demonstrated gene flow from the northeast/east towards the southwest/west. This genetic signature suggests that the observed gene flow may be due to dispersal by autumn migratory birds, which act as dispersal vectors of resistant resting propagules that form at the end of the G. semen blooms.

High Iron Requirements for Growth in the Nuisance Alga Gonyostomum semen (Raphidophyceae).

The bloom-forming freshwater alga Gonyostomum semen is associated with acidic, mesotrophic brown water lakes in boreal regions. However, researchers have been unable to conclusively link G. semen abundance and bloom formation to typical brown water lake traits, that is, high water color and DOC (dissolved organic carbon) concentrations. Iron is a main driver of water color in boreal lakes, and a recent study of lake monitoring data indicated a connection between lakes with high G. semen abundance and iron concentrations >200 µg · L-1 . Thus, iron may be the missing link in explaining G. semen abundance and growth dynamics. We experimentally assessed the effects of different iron concentrations above or below 200 µg · L-1 on the growth of G. semen batch monocultures. Iron concentrations <200 µg · L-1 limited G. semen growth, while iron concentrations >200 µg · L-1 did not. Moreover, the iron concentration of the medium required for growth was higher than for other common phytoplankton (Microcystis botrys and Chlamydomonas sp.) included in the experiment. These results indicate that G. semen requires high levels of iron in the lake environment. Consequently, this and previous findings using lake monitoring data support the hypothesis that high concentrations of iron favor the formation of high-density G. semen blooms in boreal brown water lakes. As lakes get browner in a changing climate, monitoring iron levels could be a potential tool to identify lakes at risk for G. semen blooms, especially among lakes that provide ecosystem services to society.

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.

Contrasting prevalence of selection and drift in the community structuring of bacteria and microbial eukaryotes.

Whether or not communities of microbial eukaryotes are structured in the same way as bacteria is a general and poorly explored question in ecology. Here, we investigated this question in a set of planktonic lake microbiotas in Eastern Antarctica that represent a natural community ecology experiment. Most of the analysed lakes emerged from the sea during the last 6000 years, giving rise to waterbodies that originally contained marine microbiotas and that subsequently evolved into habitats ranging from freshwater to hypersaline. We show that habitat diversification has promoted selection driven by the salinity gradient in bacterial communities (explaining ∼ 72% of taxa turnover), while microeukaryotic counterparts were predominantly structured by ecological drift (∼72% of the turnover). Nevertheless, we also detected a number of microeukaryotes with specific responses to salinity, indicating that albeit minor, selection has had a role in the structuring of specific members of their communities. In sum, we conclude that microeukaryotes and bacteria inhabiting the same communities can be structured predominantly by different processes. This should be considered in future studies aiming to understand the mechanisms that shape microbial assemblages.

Delineating closely related dinoflagellate lineages using phylotranscriptomics.

Recently radiated dinoflagellates Apocalathium aciculiferum (collected in Lake Erken, Sweden), Apocalathium malmogiense (Baltic Sea) and Apocalathium aff. malmogiense (Highway Lake, Antarctica) represent a lineage with an unresolved phylogeny. We determined their phylogenetic relationships using phylotranscriptomics based on 792 amino acid sequences. Our results showed that A. aciculiferum diverged from the other two closely related lineages, consistent with their different morphologies in cell size, relative cell length and presence of spines. We hypothesized that A. aff. malmogiense and A. malmogiense, which inhabit different hemispheres, are evolutionarily more closely related because they diverged from a marine common ancestor, adapting to a wide salinity range, while A. aciculiferum colonized a freshwater habitat, by acquiring adaptations to this environment, in particular, salinity intolerance. We show that phylotranscriptomics can resolve the phylogeny of recently diverged protists. This has broad relevance, given that many phytoplankton species are morphologically very similar, and single genes sometimes lack the information to determine species' relationships.

The potential for dispersal of microalgal resting cysts by migratory birds.

Most microalgal species are geographically widespread, but little is known about how they are dispersed. One potential mechanism for long-distance dispersal is through birds, which may transport cells internally (endozoochory) and deposit them during, or in-between, their migratory stopovers. We hypothesize that dinoflagellates, in particular resting stages, can tolerate bird digestion; that bird temperature, acidity, and retention time negatively affect dinoflagellate viability; and that recovered cysts can germinate after passage through the birds' gut, contributing to species-specific dispersal of the dinoflagellates across scales. Tolerance of two dinoflagellate species (Peridiniopsis borgei, a warm-water species and Apocalathium malmogiense, a cold-water species) to Mallard gut passage was investigated using in vitro experiments simulating the gizzard and caeca conditions. The effect of in vitro digestion and retention time on cell integrity, cell viability, and germination capacity of the dinoflagellate species was examined targeting both their vegetative and resting stages. Resting stages (cysts) of both species were able to survive simulated bird gut passage, even if their survival rate and germination were negatively affected by exposure to acidic condition and bird internal temperature. Cysts of A. malmogiense were more sensitive than P. borgei to treatments and to the presence of digestive enzymes. Vegetative cells did not survive conditions of bird internal temperature and formed pellicle cysts when exposed to gizzard-like acid conditions. We show that dinoflagellate resting cysts serve as dispersal propagules through migratory birds. Assuming a retention time of viable cysts of 2-12 h to duck stomach conditions, cysts could be dispersed 150-800 km and beyond.

Evidence of concurrent local adaptation and high phenotypic plasticity in a polar microeukaryote.

Here we investigated whether there is evidence of local adaptation in strains of an ancestrally marine dinoflagellate to the lacustrine environment they now inhabit (optimal genotypes) and/or if they have evolved phenotypic plasticity (a range of phenotypes). Eleven strains of Polarella glacialis were isolated and cultured from three different environments: the polar seas, a hyposaline and a hypersaline Antarctic lake. Local adaptation was tested by comparing growth rates of lacustrine and marine strains at their own and reciprocal site conditions. To determine phenotypic plasticity, we measured the reaction norm for salinity. We found evidence of both, limited local adaptation and higher phenotypic plasticity in lacustrine strains when compared with marine ancestors. At extreme high salinities, local lake strains outperformed other strains, and at extreme low salinities, strains from the hyposaline lake outperformed all other strains. The data suggest that lake populations may have evolved higher phenotypic plasticity in the lake habitats compared with the sea, presumably due to the high temporal variability in salinity in the lacustrine systems. Moreover, the interval of salinity tolerance differed between strains from the hyposaline and hypersaline lakes, indicating local adaptation promoted by different salinity.

Freshwater protists do not go with the flow: population structure in Gonyostomum semen independent of connectivity among lakes.

Many recent studies have found genetically differentiated populations in microorganisms despite potentially high dispersal. We designed a study to specifically examine the importance of physical dispersal barriers, i.e. geographic distance and lack of hydrological connectivity, in restricting gene flow and enhancing divergence in limnic microorganisms. We focused on the nuisance microalga Gonyostomum semen, which has recently expanded in Northern Europe and differentiated into genetically distinct populations. G. semen was sampled from six lakes distributed in two adjacent watersheds, which thereby comprised, both connected and non-connected lakes. The individual isolates were genotyped by amplified fragment length polymorphism. Several lake populations were differentiated from each other, but connectivity within watersheds could not explain the observed population genetic pattern. However, isolation by distance was moderate and might limit the gene flow among distant populations. In addition, we found low, but significant linkage disequilibrium, which indicates regular sexual recombination in this species, despite its high degree of asexual reproduction. Therefore, we conclude that the genetic properties of microalgae with occasional sexual reproduction essentially mirror regularly recombining species. Furthermore, the data indicated bottlenecks supporting the hypothesized recent range expansion of this species.

Priority effects in a planktonic bloom-forming marine diatom.

Priority effects occur when a species or genotype with earlier arrival has an advantage such that its relative abundance in the community or population is increased compared with later-arriving species. Few studies have dealt with this concept in the context of within-species competition. Skeletonema marinoi is a marine diatom that shows a high degree of genetic differentiation between populations over small geographical distances. To test whether historical events such as priority effects may have been important in inducing these patterns of population differentiation, we performed microcosm experiments with successive inoculation of different S. marinoi strains. Our results show that even in the absence of a numerical advantage, significant priority effects were evident. We propose that priority effects may be an important mechanism in initiating population genetic differentiation.

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.

Pigment composition and photoacclimation as keys to the ecological success of Gonyostomum semen (Raphidophyceae, Stramenopiles).

Aquatic habitats are usually structured by light attenuation with depth resulting in different microalgal communities, each one adapted to a certain light regime by their specific pigment composition. Several taxa contain pigments restricted to one phylogenetic group, making them useful as marker pigments in phytoplankton community studies. The nuisance and invasive freshwater microalga Gonyostomum semen (Raphidophyceae) is mainly found in brown water lakes with sharp vertical gradients in light intensity and color. However, its pigment composition and potential photoadaptations have not been comprehensively studied. We analyzed the photopigment composition of 12 genetically different strains of G. semen by high performance liquid chromatography after acclimation to different light conditions. We confirmed the pigments chl a, chl c1c2, diadinoxanthin, trans-neoxanthin, cis-neoxanthin, α and ß carotene, which have already been reported for G. semen. In addition, we identified, for the first time, the pigments violaxan-thin, zeaxanthin, and alloxanthin in this species. Alloxanthin has never been observed in raphidophytes before, suggesting differences in evolutionary plastid acquisition between freshwater lineages and the well-described marine species. The amount of total chl a per cell generally decreased with increasing light intensity. In contrast, the increasing ratios of the prominent pigments diadinoxanthin and alloxanthin per chl a with light intensity suggest photoprotective functions. In addition, we found significant variation in cell-specific pigment concentration among strains, grouped by lake of origin, which might correspond to genetic differences between strains and populations.

Single-cell genomics of a bloom-forming phytoplankton species reveals population genetic structure across continents.

The study of microbial diversity over time and space is fundamental to the understanding of their ecology and evolution. The underlying processes driving these patterns are not fully resolved but can be studied using population genetic approaches. Here we investigated the population genetic structure of Gonyostomum semen, a bloom-forming phytoplankton species, across two continents. The species appears to be expanding in Europe, whereas similar trends are not observed in the USA. Our aim was to investigate if populations of Gonyostomum semen in Europe and in the USA are genetically differentiated, if there is population genetic structure within the continents, and what the potential drivers of differentiation are. To this end, we used a novel method based on single-amplified genomes combined with Restriction-site Associated DNA sequencing that allows de novo genotyping of natural single-cell isolates without the need for culturing. We amplified over 900 single-cell genomes from 25 lake populations across Europe and the USA and identified two distinct population clusters, one in Europe and another in the USA. Low genetic diversity in European populations supports the hypothesized recent expansion of Gonyostomum semen on this continent. Geographic population structure within each continent was associated with differences in environmental variables that may have led to ecological divergence of population clusters. Overall, our results show that single-amplified genomes combined with Restriction-site Associated DNA sequencing can be used to analyze microalgal population structure and differentiation based on single-cell isolates from natural, uncultured samples.

Population genomic analyses reveal that salinity and geographic isolation drive diversification in a free-living protist.

Protists make up the vast diversity of eukaryotic life and play a critical role in biogeochemical cycling and in food webs. Because of their small size, cryptic life cycles, and large population sizes, our understanding of speciation in these organisms is very limited. We performed population genomic analyses on 153 strains isolated from eight populations of the recently radiated dinoflagellate genus Apocalathium, to explore the drivers and mechanisms of speciation processes. Species of this genus inhabit both freshwater and saline habitats, lakes and seas, and are found in cold temperate environments across the world. RAD sequencing analyses revealed that the populations were overall highly differentiated, but morphological similarity was not congruent with genetic similarity. While geographic isolation was to some extent coupled to genetic distance, this pattern was not consistent. Instead, we found evidence that the environment, specifically salinity, is a major factor in driving ecological speciation in Apocalathium. While saline populations were unique in loci coupled to genes involved in osmoregulation, freshwater populations appear to lack these. Our study highlights that adaptation to freshwater through loss of osmoregulatory genes may be an important speciation mechanism in free-living aquatic protists.

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.

Sterols of the green-pigmented, freshwater raphidophyte, Gonyostomum semen, from Scandinavian lakes.

Sterols are a class of membrane-reinforcing, ringed lipids which have a long history of examination in algae as a means of deriving chemotaxonomic relationships and as potential lipidic biomarkers. The Raphidophyceae represent a class of harmful, bloom-forming, marine and freshwater algae. To date, there have been four published examinations of their sterol composition, focusing primarily on brown-pigmented, marine species within the genera, Chattonella, Fibrocapsa, and Heterosigma. Lacking in these examinations has been the species Gonyostomum semen Ehrenb., which is a green-pigmented, freshwater raphidophyte with a worldwide distribution. The goal of this study was to examine the sterol composition of this nuisance alga, determine the potential of using its sterol profile as a biomarker, and finally to determine if there is any intraspecific variability between isolates. We have examined 21 isolates of G. semen from a number of Scandinavian lakes, and all were found to produce two major sterols, 24-ethylcholesta-5,22E-dien-3ß-ol and 24-ethylcholest-5-en-3ß-ol, and 24-methylcholest-5-en-3ß-ol as a minor sterol; the presence of 24-ethylcholesta-5,22E-dien-3ß-ol differentiates G. semen from brown-pigmented, marine raphidophytes which generally lack it. The results of this study indicate that isolates of G. semen from geographically separate lakes across Finland and Scandinavia have the same sterol biosynthetic pathway, and that there is no evolutionary divergence between the isolates with regard to sterol composition. The sterols of G. semen are not considered to be useful biomarkers for this particular organism because they are commonly found in other algae and plants.

Local adaptation through countergradient selection in northern populations of Skeletonema marinoi.

Marine microorganisms have the potential to disperse widely with few obvious barriers to gene flow. However, among microalgae, several studies have demonstrated that species can be highly genetically structured with limited gene flow among populations, despite hydrographic connectivity. Ecological differentiation and local adaptation have been suggested as drivers of such population structure. Here we tested whether multiple strains from two genetically distinct Baltic Sea populations of the diatom Skeletonema marinoi showed evidence of local adaptation to their local environments: the estuarine Bothnian Sea and the marine Kattegat Sea. We performed reciprocal transplants of multiple strains between culture media based on water from the respective environments, and we also allowed competition between strains of estuarine and marine origin in both salinities. When grown alone, both marine and estuarine strains performed best in the high-salinity environment, and estuarine strains always grew faster than marine strains. This result suggests local adaptation through countergradient selection, that is, genetic effects counteract environmental effects. However, the higher growth rate of the estuarine strains appears to have a cost in the marine environment and when strains were allowed to compete, marine strains performed better than estuarine strains in the marine environment. Thus, other traits are likely to also affect fitness. We provide evidence that tolerance to pH could be involved and that estuarine strains that are adapted to a more fluctuating pH continue growing at higher pH than marine strains.

Single-cell adaptations shape evolutionary transitions to multicellularity in green algae.

The evolution of multicellular life has played a pivotal role in shaping biological diversity. However, we know surprisingly little about the natural environmental conditions that favour the formation of multicellular groups. Here we experimentally examine how key environmental factors (predation, nitrogen and water turbulence) combine to influence multicellular group formation in 35 wild unicellular green algae strains (19 Chlorophyta species). All environmental factors induced the formation of multicellular groups (more than four cells), but there was no evidence this was adaptive, as multicellularity (% cells in groups) was not related to population growth rate under any condition. Instead, population growth was related to extracellular matrix (ECM) around single cells and palmelloid formation, a unicellular life-cycle stage where two to four cells are retained within a mother-cell wall after mitosis. ECM production increased with nitrogen levels resulting in more cells being in palmelloids and higher rates of multicellular group formation. Examining the distribution of 332 algae species across 478 lakes monitored over 55 years, showed that ECM and nitrogen availability also predicted patterns of obligate multicellularity in nature. Our results highlight that adaptations of unicellular organisms to cope with environmental challenges may be key to understanding evolutionary routes to multicellular life.