Temporal genomics help in deciphering neutral and adaptive patterns in the contemporary evolution of kelp populations.

The impact of climate change on populations will be contingent upon their contemporary adaptive evolution. In this study, we investigated the contemporary evolution of four populations of the cold-water kelp Laminaria digitata by analysing their spatial and temporal genomic variation using ddRAD-sequencing. These populations were sampled from the center to the southern margin of its north-eastern Atlantic distribution at two-time points, spanning at least two generations. Through genome scans for local adaptation at a single time point, we identified candidate loci that showed clinal variation correlated with changes in sea surface temperature (SST) along latitudinal gradients. This finding suggests that SST may drive the adaptive response of these kelp populations, although factors such as species' demographic history should also be considered. Additionally, we performed a simulation approach to distinguish the effect of selection from genetic drift in allele frequency changes over time. This enabled the detection of loci in the southernmost population that exhibited temporal differentiation beyond what would be expected from genetic drift alone: these are candidate loci which could have evolved under selection over time. In contrast, we did not detect any outlier locus based on temporal differentiation in the population from the North Sea, which also displayed low and decreasing levels of genetic diversity. The diverse evolutionary scenarios observed among populations can be attributed to variations in the prevalence of selection relative to genetic drift across different environments. Therefore, our study highlights the potential of temporal genomics to offer valuable insights into the contemporary evolution of marine foundation species facing climate change.

Red macroalgae in the genomic era.

Rhodophyta (or red algae) are a diverse and species-rich group that forms one of three major lineages in the Archaeplastida, a eukaryotic supergroup whose plastids arose from a single primary endosymbiosis. Red algae are united by several features, such as relatively small intron-poor genomes and a lack of cytoskeletal structures associated with motility like flagella and centrioles, as well as a highly efficient photosynthetic capacity. Multicellular red algae (or macroalgae) are one of the earliest diverging eukaryotic lineages to have evolved complex multicellularity, yet despite their ecological, evolutionary, and commercial importance, they have remained a largely understudied group of organisms. Considering the increasing availability of red algal genome sequences, we present a broad overview of fundamental aspects of red macroalgal biology and posit on how this is expected to accelerate research in many domains of red algal biology in the coming years.

The Rhodoexplorer Platform for Red Algal Genomics and Whole-Genome Assemblies for Several Gracilaria Species.

Macroalgal (seaweed) genomic resources are generally lacking as compared with other eukaryotic taxa, and this is particularly true in the red algae (Rhodophyta). Understanding red algal genomes is critical to understanding eukaryotic evolution given that red algal genes are spread across eukaryotic lineages from secondary endosymbiosis and red algae diverged early in the Archaeplastids. The Gracilariales is a highly diverse and widely distributed order including species that can serve as ecosystem engineers in intertidal habitats and several notorious introduced species. The genus Gracilaria is cultivated worldwide, in part for its production of agar and other bioactive compounds with downstream pharmaceutical and industrial applications. This genus is also emerging as a model for algal evolutionary ecology. Here, we report new whole-genome assemblies for two species (Gracilaria chilensis and Gracilaria gracilis), a draft genome assembly of Gracilaria caudata, and genome annotation of the previously published Gracilaria vermiculophylla genome. To facilitate accessibility and comparative analysis, we integrated these data in a newly created web-based portal dedicated to red algal genomics (https://rhodoexplorer.sb-roscoff.fr). These genomes will provide a resource for understanding algal biology and, more broadly, eukaryotic evolution.

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.

Evolution of life cycles and reproductive traits: Insights from the brown algae.

A vast diversity of types of life cycles exists in nature, and several theories have been advanced to explain how this diversity has evolved and how each type of life cycle is retained over evolutionary time. Here, we exploited the diversity of life cycles and reproductive traits of the brown algae (Phaeophyceae) to test several hypotheses on the evolution of life cycles. We investigated the evolutionary dynamics of four life-history traits: life cycle, sexual system, level of gamete dimorphism and gamete parthenogenetic capacity. We assigned states to up to 77 representative species of the taxonomic diversity of the brown algal group, in a multi-gene phylogeny. We used maximum likelihood and Bayesian analyses of correlated evolution, while taking the phylogeny into account, to test for correlations between traits and to investigate the chronological sequence of trait acquisition. Our analyses are consistent with the prediction that diploid growth evolves when sexual reproduction is preferred over asexual reproduction, possibly because it allows the complementation of deleterious mutations. We also found that haploid sex determination is ancestral in relation to diploid sex determination. However, our results could not address whether increased zygotic and diploid growth are associated with increased sexual dimorphism. Our analyses suggest that in the brown algae, isogamous species evolved from anisogamous ancestors, contrary to the commonly reported pattern where evolution proceeds from isogamy to anisogamy.

Exploring the Genetic Consequences of Clonality in Haplodiplontic Taxa.

Partially clonality is an incredibly common reproductive mode found across all the major eukaryotic lineages. Yet, population genetic theory is based on exclusive sexuality or exclusive asexuality, and partial clonality is often ignored. This is particularly true in haplodiplontic eukaryotes, including algae, ferns, mosses, and fungi, where somatic development occurs in both the haploid and diploid stages. Haplodiplontic life cycles are predicted to be correlated with asexuality, but tests of this prediction are rare. Moreover, there are unique consequences of having long-lived haploid and diploid stages in the same life cycle. For example, clonal processes uncouple the life cycle such that the repetition of the diploid stage via clonality leads to the loss of the haploid stage. Here, we surveyed the literature to find studies that had genotyped both haploid and diploid stages and recalculated population genetic summary metrics for seven red algae, one green alga, three brown algae, and three mosses. We compared these data to recent simulations that explicitly addressed the population genetic consequences of partial clonality in haplodiplontic life cycles. Not only was partial clonality found to act as a homogenizing force, but the combined effects of proportion of haploids, rate of clonality, and the relative strength of mutation versus genetic drift impacts the distributions of population genetic indices. We found remarkably similar patterns across commonly used population genetic metrics between our empirical and recent theoretical expectations. To facilitate future studies, we provide some recommendations for sampling and analyzing population genetic parameters for haplodiplontic taxa.

Seascape Genomics of the Sugar Kelp Saccharina latissima along the North Eastern Atlantic Latitudinal Gradient.

Temperature is one of the most important range-limiting factors for many seaweeds. Driven by the recent climatic changes, rapid northward shifts of species' distribution ranges can potentially modify the phylogeographic signature of Last Glacial Maximum. We explored this question in detail in the cold-tolerant kelp species Saccharina latissima, using microsatellites and double digest restriction site-associated DNA sequencing ( ddRAD-seq) derived single nucleotide polymorphisms (SNPs) to analyze the genetic diversity and structure in 11 sites spanning the entire European Atlantic latitudinal range of this species. In addition, we checked for statistical correlation between genetic marker allele frequencies and three environmental proxies (sea surface temperature, salinity, and water turbidity). Our findings revealed that genetic diversity was significantly higher for the northernmost locality (Spitsbergen) compared to the southern ones (Northern Iberia), which we discuss in light of the current state of knowledge on phylogeography of S. latissima and the potential influence of the recent climatic changes on the population structure of this species. Seven SNPs and 12 microsatellite alleles were found to be significantly associated with at least one of the three environmental variables. We speculate on the putative adaptive functions of the genes associated with the outlier markers and the importance of these markers for successful conservation and aquaculture strategies for S. latissima in this age of rapid global change.

After a catastrophe, a little bit of sex is better than nothing: Genetic consequences of a major earthquake on asexual and sexual populations.

Catastrophic events can have profound effects on the demography of a population and consequently on genetic diversity. The dynamics of postcatastrophic recovery and the role of sexual versus asexual reproduction in buffering the effects of massive perturbations remain poorly understood, in part because the opportunity to document genetic diversity before and after such events is rare. Six natural (purely sexual) and seven cultivated (mainly clonal due to farming practices) populations of the red alga Agarophyton chilense were surveyed along the Chilean coast before, in the days after and 2 years after the 8.8 magnitude earthquake in 2010. The genetic diversity of sexual populations appeared sensitive to this massive perturbation, notably through the loss of rare alleles immediately after the earthquake. By 2012, the levels of diversity returned to those observed before the catastrophe, probably due to migration. In contrast, enhanced rates of clonality in cultivated populations conferred a surprising ability to buffer the instantaneous loss of diversity. After the earthquake, farmers increased the already high rate of clonality to maintain the few surviving beds, but most of them collapsed rapidly. Contrasting fates between sexual and clonal populations suggest that betting on strict clonality to sustain production is risky, probably because this extreme strategy hampered adaptation to the brutal environmental perturbation induced by the catastrophe.

Species specific gene expression dynamics during harmful algal blooms.

Harmful algal blooms are caused by specific members of microbial communities. Understanding the dynamics of these events requires comparing the strategies developed by the problematic species to cope with environmental fluctuations to the ones developed by the other members of the community. During three consecutive years, the meta-transcriptome of micro-eukaryote communities was sequenced during blooms of the toxic dinoflagellate Alexandrium minutum. The dataset was analyzed to investigate species specific gene expression dynamics. Major shifts in gene expression were explained by the succession of different species within the community. Although expression patterns were strongly correlated with fluctuation of the abiotic environment, and more specifically with nutrient concentration, transcripts specifically involved in nutrient uptake and metabolism did not display extensive changes in gene expression. Compared to the other members of the community, A. minutum displayed a very specific expression pattern, with lower expression of photosynthesis transcripts and central metabolism genes (TCA cycle, glucose metabolism, glycolysis…) and contrasting expression pattern of ion transporters across environmental conditions. These results suggest the importance of mixotrophy, cell motility and cell-to-cell interactions during A. minutum blooms.

Dual influence of terrestrial and marine historical processes on the phylogeography of the Brazilian intertidal red alga Gracilaria caudata.

In this study, we explored how past terrestrial and marine climate changes have interacted to shape the phylogeographic patterns of the intertidal red seaweed Gracilaria caudata, an economically important species exploited for agar production in the Brazilian north-east. Seven sites were sampled along the north-east tropical and south-east sub-tropical Brazilian coast. The genetic diversity and structure of G. caudata was inferred using a combination of mitochondrial (COI and cox2-3), chloroplast (rbcL) and 15 nuclear microsatellite markers. A remarkable congruence between nuclear, mitochondrial and chloroplast data revealed clear separation between the north-east (from 03° S to 08° S) and the south-east (from 20° S to 23° S) coast of Brazil. These two clades differ in their demographic histories, with signatures of recent demographic expansions in the north-east and divergent populations in the south-east, suggesting the maintenance of several refugia during the last glacial maximum due to sea-level rise and fall. The Bahia region (around 12° S) occupies an intermediate position between both clades. Microsatellites and mtDNA markers showed additional levels of genetic structure within each sampled site located south of Bahia. The separation between the two main groups in G. caudata is likely recent, probably occurring during the Quaternary glacial cycles. The genetic breaks are concordant with (i) those separating terrestrial refugia, (ii) major river outflows and (iii) frontiers between tropical and subtropical regions. Taken together with previously published eco-physiological studies that showed differences in the physiological performance of the strains from distinct locations, these results suggest that the divergent clades in G. caudata correspond to distinct ecotypes in the process of incipient speciation and thus should be considered for the management policy of this commercially important species.

Combining niche shift and population genetic analyses predicts rapid phenotypic evolution during invasion.

The rapid evolution of non-native species can facilitate invasion success, but recent reviews indicate that such microevolution rarely yields expansion of the climatic niche in the introduced habitats. However, because some invasions originate from a geographically restricted portion of the native species range and its climatic niche, it is possible that the frequency, direction, and magnitude of phenotypic evolution during invasion have been underestimated. We explored the utility of niche shift analyses in the red seaweed Gracilaria vermiculophylla, which expanded its range from the northeastern coastline of Japan to North America, Europe, and northwestern Africa within the last 100 years. A genetically informed climatic niche shift analysis indicates that native source populations occur in colder and highly seasonal habitats, while most non-native populations typically occur in warmer, less seasonal habitats. This climatic niche expansion predicts that non-native populations evolved greater tolerance for elevated heat conditions relative to native source populations. We assayed 935 field-collected and 325 common-garden thalli from 40 locations, and as predicted, non-native populations had greater tolerance for ecologically relevant extreme heat (40°C) than did Japanese source populations. Non-native populations also had greater tolerance for cold and low-salinity stresses relative to source populations. The importance of local adaptation to warm temperatures during invasion was reinforced by evolution of parallel clines: Populations from warmer, lower-latitude estuaries had greater heat tolerance than did populations from colder, higher-latitude estuaries in both Japan and eastern North America. We conclude that rapid evolution plays an important role in facilitating the invasion success of this and perhaps other non-native marine species. Genetically informed ecological niche analyses readily generate clear predictions of phenotypic shifts during invasions and may help to resolve debate over the frequency of niche conservatism versus rapid adaptation during invasion.

ClonEstiMate, a Bayesian method for quantifying rates of clonality of populations genotyped at two-time steps.

Partial clonality is commonly used in eukaryotes and has large consequences for their evolution and ecology. Assessing accurately the relative importance of clonal vs. sexual reproduction matters for studying and managing such species. Here, we proposed a Bayesian approach, ClonEstiMate, to infer rates of clonality c from populations sampled twice over a short time interval, ideally one generation time. The method relies on the likelihood of the transitions between genotype frequencies of ancestral and descendent populations, using an extended Wright-Fisher model explicitly integrating reproductive modes. Our model provides posterior probability distribution of inferred c, given the assumed rates of mutation, as well as inbreeding and selfing when occurring. Tested under various conditions, this model provided accurate inferences of c, especially when the amount of information was modest, that is low sample sizes, few loci, low polymorphism and strong linkage disequilibrium. Inferences remained robust when mutation models and rates were misinformed. However, the method was sensitive to moderate frequencies of null alleles and when the time interval between required samplings exceeding two generations. Misinformed rates on mating modes (inbreeding and selfing) also resulted in biased inferences. Our method was tested on eleven data sets covering five partially clonal species, for which the extent of clonality was formerly deciphered. It delivered highly consistent results with previous information on the biology of those species. ClonEstiMate represents a powerful tool for detecting and inferring clonality in finite populations, genotyped with SNPs or microsatellites. It is freely available at https://www6.rennes.inra.fr/igepp_eng/Productions/Software.

Evidence for parasite-mediated selection during short-lasting toxic algal blooms.

Parasites play a role in the control of transient algal blooms, but it is not known whether parasite-mediated selection results in coevolution of the host and the parasites over this short time span. We investigated the presence of coevolution between the toxic dinoflagellate Alexandrium minutum and two naturally occurring endoparasites during blooms lasting a month in two river estuaries, using cross-inoculation experiments across time and space. Higher parasite abundance was associated with a large daily reduction in relative A. minutum abundances, demonstrating strong parasite-mediated selection. There was genetic variability in infectivity in both parasite species, and in resistance in the host. We found no evidence for coevolution in one estuary; however, in the other estuary, we found high genetic diversity in the two parasite species, fluctuations in infectivity and suggestion that the two parasites are well adapted to their host, as in 'Red Queen' dynamics. Thus, coevolution is possible over the short time span of a bloom, but geographically variable, and may feedback on community dynamics.

Patterns of genetic diversity of the cryptogenic red alga Polysiphonia morrowii (Ceramiales, Rhodophyta) suggest multiple origins of the Atlantic populations.

The red alga Polysiphonia morrowii, native to the North Pacific (Northeast Asia), has recently been reported worldwide. To determine the origin of the French and Argentine populations of this introduced species, we compared samples from these two areas with samples collected in Korea and at Hakodate, Japan, the type locality of the species. Combined analyses of chloroplastic (rbcL) and mitochondrial (cox1) DNA revealed that the French and Argentine populations are closely related and differ substantially from the Korean and Japanese populations. The genetic structure of P. morrowii populations from South Atlantic and North Atlantic, which showed high haplotype diversity compared with populations from the North Pacific, suggested the occurrence of multiple introduction events from areas outside of the so-called native regions. Although similar, the French and Argentine populations are not genetically identical. Thus, the genetic structure of these two introduced areas may have been modified by cryptic and recurrent introduction events directly from Asia or from other introduced areas that act as introduction relays. In addition, the large number of private cytoplasmic types identified in the two introduced regions strongly suggests that local populations of P. morrowii existed before the recent detection of these invasions. Our results suggest that the most likely scenario is that the source population(s) of the French and Argentine populations was not located only in the North Pacific and/or that P. morrowii is a cryptogenic species.

Evolutionary processes and cellular functions underlying divergence in Alexandrium minutum.

Understanding divergence in the highly dispersive and seemingly homogeneous pelagic environment for organisms living as free drifters in the water column remains a challenge. Here, we analysed the transcriptome-wide mRNA sequences, as well as the morphology of 18 strains of Alexandrium minutum, a dinoflagellate responsible for harmful algal blooms worldwide, to investigate the functional bases of a divergence event. Analysis of the joint site frequency spectrum (JSFS) pointed towards an ancestral divergence in complete isolations followed by a secondary contact resulting in gene flow between the two diverging groups, but heterogeneous across sites. The sites displaying fixed SNPs were associated with a highly restricted gene flow and a strong overrepresentation of nonsynonymous polymorphism, suggesting the importance of selective pressures as drivers of the divergence. The most divergent transcripts were homologs to genes involved in calcium/potassium fluxes across the membrane, calcium transduction signal and saxitoxin production. The implication of these results in terms of ecological divergence and build-up of reproductive isolation is discussed. Dinoflagellates are especially difficult to study in the field at the ecological level due to their small size and the dynamic nature of their natural environment, but also at the genomic level due to their huge and complex genome and the absence of closely related model organism. This study illustrates the possibility to identify the traits of primary importance in ecology and evolution starting from high-throughput sequencing data, even for such organisms.

Evolution and maintenance of haploid-diploid life cycles in natural populations: The case of the marine brown alga Ectocarpus.

The evolutionary stability of haploid-diploid life cycles is still controversial. Mathematical models indicate that niche differences between ploidy phases may be a necessary condition for the evolution and maintenance of these life cycles. Nevertheless, experimental support for this prediction remains elusive. In the present work, we explored this hypothesis in natural populations of the brown alga Ectocarpus. Consistent with the life cycle described in culture, Ectocarpus crouaniorum in NW France and E. siliculosus in SW Italy exhibited an alternation between haploid gametophytes and diploid sporophytes. Our field data invalidated, however, the long-standing view of an isomorphic alternation of generations. Gametophytes and sporophytes displayed marked differences in size and, conforming to theoretical predictions, occupied different spatiotemporal niches. Gametophytes were found almost exclusively on the alga Scytosiphon lomentaria during spring whereas sporophytes were present year-round on abiotic substrata. Paradoxically, E. siliculosus in NW France exhibited similar habitat usage despite the absence of alternation of ploidy phases. Diploid sporophytes grew both epilithically and epiphytically, and this mainly asexual population gained the same ecological advantage postulated for haploid-diploid populations. Consequently, an ecological interpretation of the niche differences between haploid and diploid individuals does not seem to satisfactorily explain the evolution of the Ectocarpus life cycle.

Species are hypotheses: avoid connectivity assessments based on pillars of sand.

Connectivity among populations determines the dynamics and evolution of populations, and its assessment is essential in ecology in general and in conservation biology in particular. The robust basis of any ecological study is the accurate delimitation of evolutionary units, such as populations, metapopulations and species. Yet a disconnect still persists between the work of taxonomists describing species as working hypotheses and the use of species delimitation by molecular ecologists interested in describing patterns of gene flow. This problem is particularly acute in the marine environment where the inventory of biodiversity is relatively delayed, while for the past two decades, molecular studies have shown a high prevalence of cryptic species. In this study, we illustrate, based on marine case studies, how the failure to recognize boundaries of evolutionary-relevant unit leads to heavily biased estimates of connectivity. We review the conceptual framework within which species delimitation can be formalized as falsifiable hypotheses and show how connectivity studies can feed integrative taxonomic work and vice versa. Finally, we suggest strategies for spatial, temporal and phylogenetic sampling to reduce the probability of inadequately delimiting evolutionary units when engaging in connectivity studies.

Molecular evidence for the coexistence of two sibling species in Pylaiella littoralis (Ectocarpales, Phaeophyceae) along the Brittany coast.

The great phenotypic variability and the lack of diagnostic characters in the genus Pylaiella render the systematic study of this genus problematic. In this study, we investigated the diversity of Pylaiella littoralis along the Brittany (France) coast using a DNA barcoding multilocus approach with mitochondrial (cox1, nad1, and atp9) and chloroplastic (rbcL and atpB) markers associated with a population genetics approach using 10 microsatellite markers. In addition, spatio-temporal sampling was conducted along the Brittany coast. We sampled 140 individuals from four sites located between Saint-Malo and Concarneau (380 km) from April to October. Mitochondrial sequence data revealed the occurrence of two sibling species, with a minimum of 2.4% divergence between them. Microsatellite genotypic data congruently revealed two well-supported clusters matching the two mitochondrial clades of Pylaiella. Although gene flow is limited between species, occurrence of genetic admixtures in some populations suggested that reproductive isolation is not complete. Our study highlighted the complementarity of barcoding and population genetics approaches to shed light on the evolutionary processes that lead to speciation.

Tracing the trans-pacific evolutionary history of a domesticated Seaweed (Gracilaria chilensis) with archaeological and genetic data.

The history of a domesticated marine macroalga is studied using archaeological, phylogeographic and population genetic tools. Phylogeographic and population genetic analyses demonstrated that the cultivated red alga Gracilaria chilensis colonised the Chilean coast from New Zealand. Combining archaeological observations with phylogeographic data provided evidence that exchanges between New Zealand and Chile have occurred at least before the Holocene, likely at the end of the Last Glacial Maximum (LGM) and we suggest that migration probably occurred via rafting. Furthermore, the remarkably low microsatellite diversity found in the Chilean populations compared to those in New Zealand is consistent with a recent genetic bottleneck as a result of over-exploitation of natural populations and/or the process of domestication. Therefore, the aquaculture of this seaweed, based essentially on clonal propagation, is occurring from genetically depressed populations and may be driving the species to an extinction vortex in Chile.

Alteration of sexual reproduction and genetic diversity in the kelp species Laminaria digitata at the southern limit of its range.

Adaptation to marginal habitats at species range-limits has often been associated with parthenogenetic reproduction in terrestrial animals and plants. Laboratory observations have shown that brown algae exhibit a high propensity for parthenogenesis by various mechanisms. The kelp Laminaria digitata is an important component of the ecosystem in Northern European rocky intertidal habitats. We studied four L. digitata populations for the effects of marginality on genetic diversity and sexual reproduction. Two populations were marginal: One (Locquirec, in Northern Brittany) was well within the geographic range, but was genetically isolated from other populations by large stretches of sandy beaches. Another population was at the range limits of the species (Quiberon, in Southern Brittany) and was exposed to much higher seasonal temperature changes. Microsatellite analyses confirmed that these populations showed decreased genetic and allelic diversity, consistent with marginality and genetic isolation. Sporophytes from both marginal populations showed greatly diminished spore-production compared to central populations, but only the southern-limit population (Quiberon) showed a high propensity for producing unreduced (2N) spores. Unreduced 2N spores formed phenotypically normal gametophytes with nuclear area consistent with ≥2N DNA contents, and microsatellite studies suggested these were produced at least in part by automixis. However, despite this being the dominant path of spore production in Quiberon sporophyte individuals, the genetic evidence indicated the population was maintained mostly by sexual reproduction. Thus, although spore production and development showed the expected tendency of geographical parthenogenesis in marginal populations, this appeared to be a consequence of maladaptation, rather than an adaptation to, life in a marginal habitat.