Genomic scans detect signatures of selection along a salinity gradient in populations of the intertidal seaweed Fucus serratus on a 12 km scale

Local thermal adaptation and limited gene flow constrain future climate responses of a marine ecosystem engineer

Rising ocean temperatures and extreme temperature events have precipitated declines and local extinctions in many marine species globally, but patterns of loss are often uneven across species ranges for reasons that are poorly understood. Knowledge of the extent of local adaptation and gene flow may explain such patterns and help predict future trajectories under scenarios of climate change. We test the extent to which local differentiation in thermal tolerance is influenced by gene flow and local adaptation using a widely distributed intertidal seaweed (Hormosira banksii) from temperate Australia. Population surveys across ~2,000 km of the species range revealed strong genetic structuring at regional and local scales (global F ST = 0.243) reflecting extremely limited gene flow, while common garden experiments (14-day exposures to 15, 18, 21°C) revealed strong site differences in early development and mortality in response to elevated temperature. Embryos from many sites spanning a longitudinal thermal gradient showed suppressed development and increased mortality to elevated water temperatures, but populations originating from warmer and more variable thermal environments tended to be less susceptible to warming. Notably, there was significant local-scale variation in the thermal responses of embryos within regions which was corroborated by the finding of small-scale genetic differences. We expect the observed genetic and phenotypic differentiation to lead to uneven responses to warming sea surface temperatures in this important marine foundation species. The study highlights the challenges of predicting species responses to thermal stress and the importance of management strategies that incorporate evolutionary potential for "climate-proofing" marine ecosystems.

Historical isolation and contemporary gene flow drive population diversity of the brown alga Sargassum thunbergii along the coast of China

BACKGROUND

Long-term survival in isolated marginal seas of the China coast during the late Pleistocene ice ages is widely believed to be an important historical factor contributing to population genetic structure in coastal marine species. Whether or not contemporary factors (e.g. long-distance dispersal via coastal currents) continue to shape diversity gradients in marine organisms with high dispersal capability remains poorly understood. Our aim was to explore how historical and contemporary factors influenced the genetic diversity and distribution of the brown alga Sargassum thunbergii, which can drift on surface water, leading to long-distance dispersal.

RESULTS

We used 11 microsatellites and the plastid RuBisCo spacer to evaluate the genetic diversity of 22 Sargassum thunbergii populations sampled along the China coast. Population structure and differentiation was inferred based on genotype clustering and pairwise F ST and allele-frequency analyses. Integrated genetic analyses revealed two genetic clusters in S. thunbergii that dominated in the Yellow-Bohai Sea (YBS) and East China Sea (ECS) respectively. Higher levels of genetic diversity and variation were detected among populations in the YBS than in the ECS. Bayesian coalescent theory was used to estimate contemporary and historical gene flow. High levels of contemporary gene flow were detected from the YBS (north) to the ECS (south), whereas low levels of historical gene flow occurred between the two regions.

CONCLUSIONS

Our results suggest that the deep genetic divergence in S. thunbergii along the China coast may result from long-term geographic isolation during glacial periods. The dispersal of S. thunbergii driven by coastal currents may facilitate the admixture between southern and northern regimes. Our findings exemplify how both historical and contemporary forces are needed to understand phylogeographical patterns in coastal marine species with long-distance dispersal.

High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga Ectocarpus

Deciphering the genetic architecture of adaptation of brown algae to environmental stresses such as temperature and salinity is of evolutionary as well as of practical interest. The filamentous brown alga Ectocarpus sp. is a model for the brown algae and its genome has been sequenced. As sessile organisms, brown algae need to be capable of resisting the various abiotic stressors that act in the intertidal zone (e.g. osmotic pressure, temperature, salinity, UV radiation) and previous studies have shown that an important proportion of the expressed genes is regulated in response to hyposaline, hypersaline or oxidative stress conditions. Using the double digest RAD sequencing method, we constructed a dense genetic map with 3,588 SNP markers and identified 39 QTLs for growth-related traits and their plasticity under different temperature and salinity conditions (tolerance to high temperature and low salinity). GO enrichment tests within QTL intervals highlighted membrane transport processes such as ion transporters. Our study represents a significant step towards deciphering the genetic basis of adaptation of Ectocarpus sp. to stress conditions and provides a substantial resource to the increasing list of tools generated for the species.

Environmental versus geographical determinants of genetic structure in two subalpine conifers

Alpine ecosystems are facing rapid human-induced environmental changes, and so more knowledge about tree adaptive potential is needed. This study investigated the relative role of isolation by distance (IBD) versus isolation by adaptation (IBA) in explaining population genetic structure in Abies alba and Larix decidua, based on 231 and 233 single nucleotide polymorphisms (SNPs) sampled across 36 and 22 natural populations, respectively, in the Alps and Apennines. Genetic structure was investigated for both geographical and environmental groups, using analysis of molecular variance (AMOVA). For each species, nine environmental groups were defined using climate variables selected from a multiple factor analysis. Complementary methods were applied to identify outliers based on these groups, and to test for IBD versus IBA. AMOVA showed weak but significant genetic structure for both species, with higher values in L. decidua. Among the potential outliers detected, up to two loci were found for geographical groups and up to seven for environmental groups. A stronger effect of IBD than IBA was found in both species; nevertheless, once spatial effects had been removed, temperature and soil in A. alba, and precipitation in both species, were relevant factors explaining genetic structure. Based on our findings, in the Alpine region, genetic structure seems to be affected by both geographical isolation and environmental gradients, creating opportunities for local adaptation.

High and distinct range-edge genetic diversity despite local bottlenecks

The genetic consequences of living on the edge of distributional ranges have been the subject of a largely unresolved debate. Populations occurring along persistent low latitude ranges (rear-edge) are expected to retain high and unique genetic diversity. In contrast, currently less favourable environmental conditions limiting population size at such range-edges may have caused genetic erosion that prevails over past historical effects, with potential consequences on reducing future adaptive capacity. The present study provides an empirical test of whether population declines towards a peripheral range might be reflected on decreasing diversity and increasing population isolation and differentiation. We compare population genetic differentiation and diversity with trends in abundance along a latitudinal gradient towards the peripheral distribution range of Saccorhiza polyschides, a large brown seaweed that is the main structural species of kelp forests in SW Europe. Signatures of recent bottleneck events were also evaluated to determine whether the recently recorded distributional shifts had a negative influence on effective population size. Our findings show decreasing population density and increasing spatial fragmentation and local extinctions towards the southern edge. Genetic data revealed two well supported groups with a central contact zone. As predicted, higher differentiation and signs of bottlenecks were found at the southern edge region. However, a decrease in genetic diversity associated with this pattern was not verified. Surprisingly, genetic diversity increased towards the edge despite bottlenecks and much lower densities, suggesting that extinctions and recolonizations have not strongly reduced diversity or that diversity might have been even higher there in the past, a process of shifting genetic baselines.

Testing the consistency of connectivity patterns for a widely dispersing marine species

Connectivity is widely recognized as an important component in developing effective management and conservation strategies. Although managers are generally most interested in demographic, rather than genetic connectivity, new analytic approaches are able to provide estimates of both demographic and genetic connectivity measures from genetic data. Combining such genetic data with mathematical models represents a powerful approach for accurately determining patterns of population connectivity. Here, we use microsatellite markers to investigate the genetic population structure of the New Zealand Rock Lobster, Jasus edwardsii, which has one of the longest known larval durations of all marine species (>2 years), a very large geographic range (>5500 km), and has been the subject of extensive dispersal modeling. Despite earlier mitochondrial DNA studies finding homogeneous genetic structure, the mathematical model suggests that there are source-sink dynamics for this species. We found evidence of genetic structure in J. edwardsii populations with three distinct genetic groups across New Zealand and a further Australian group; these groups and patterns of gene flow were generally congruent with the earlier mathematical model. Of particular interest was the consistent identification of a self-recruiting population/region from both modeling and genetic approaches. Although there is the potential for selection and harvesting to influence the patterns we observed, we believe oceanographic processes are most likely responsible for the genetic structure observed in J. edwardsii. Our results, using a species at the extreme end of the dispersal spectrum, demonstrate that source-sink population dynamics may still exist for such species.

Climate change impact on seaweed meadow distribution in the North Atlantic rocky intertidal

The North-Atlantic has warmed faster than all other ocean basins and climate change scenarios predict sea surface temperature isotherms to shift up to 600 km northwards by the end of the 21st century. The pole-ward shift has already begun for many temperate seaweed species that are important intertidal foundation species. We asked the question: Where will climate change have the greatest impact on three foundational, macroalgal species that occur along North-Atlantic shores: Fucus serratus, Fucus vesiculosus, and Ascophyllum nodosum? To predict distributional changes of these key species under three IPCC (Intergovernmental Panel on Climate Change) climate change scenarios (A2, A1B, and B1) over the coming two centuries, we generated Ecological Niche Models with the program MAXENT. Model predictions suggest that these three species will shift northwards as an assemblage or "unit" and that phytogeographic changes will be most pronounced in the southern Arctic and the southern temperate provinces. Our models predict that Arctic shores in Canada, Greenland, and Spitsbergen will become suitable for all three species by 2100. Shores south of 45° North will become unsuitable for at least two of the three focal species on both the Northwest- and Northeast-Atlantic coasts by 2200. If these foundational species are unable to adapt to the rising temperatures, they will lose their centers of genetic diversity and their loss will trigger an unpredictable shift in the North-Atlantic intertidal ecosystem.

Ecological and evolutionary genomics of marine photosynthetic organisms

Environmental (ecological) genomics aims to understand the genetic basis of relationships between organisms and their abiotic and biotic environments. It is a rapidly progressing field of research largely due to recent advances in the speed and volume of genomic data being produced by next generation sequencing (NGS) technologies. Building on information generated by NGS-based approaches, functional genomic methodologies are being applied to identify and characterize genes and gene systems of both environmental and evolutionary relevance. Marine photosynthetic organisms (MPOs) were poorly represented amongst the early genomic models, but this situation is changing rapidly. Here we provide an overview of the recent advances in the application of ecological genomic approaches to both prokaryotic and eukaryotic MPOs. We describe how these approaches are being used to explore the biology and ecology of marine cyanobacteria and algae, particularly with regard to their functions in a broad range of marine ecosystems. Specifically, we review the ecological and evolutionary insights gained from whole genome and transcriptome sequencing projects applied to MPOs and illustrate how their genomes are yielding information on the specific features of these organisms.

Range-edge genetic diversity: locally poor extant southern patches maintain a regionally diverse hotspot in the seagrass Zostera marina

Refugial populations at the rear edge are predicted to contain higher genetic diversity than those resulting from expansion, such as in post-glacial recolonizations. However, peripheral populations are also predicted to have decreased diversity compared to the centre of a species' distribution. We aim to test these predictions by comparing genetic diversity in populations at the limits of distribution of the seagrass Zostera marina, with populations in the species' previously described central diversity 'hotspot'. Zostera marina populations show decreased allelic richness, heterozygosity and genotypic richness in both the 'rear' edge and the 'leading' edge compared to the diversity 'hotspot' in the North Sea/Baltic region. However, when populations are pooled, genetic diversity at the southern range is as high as in the North Sea/Baltic region while the 'leading edge' remains low in genetic diversity. The decreased genetic diversity in these southern Iberian populations compared to more central populations is possibly the effect of drift because of small effective population size, as a result of reduced habitat, low sexual reproduction and low gene flow. However, when considering the whole southern edge of distribution rather than per population, diversity is as high as in the central 'hotspot' in the North Sea/Baltic region. We conclude that diversity patterns assessed per population can mask the real regional richness that is typical of rear edge populations, which have played a key role in the species biogeographical history and as marginal diversity hotspots have very high conservation value.

In silico mining and characterization of 12 EST-SSRs for the invasive slipper limpet Crepidula fornicata

In silico mining of an expressed sequence tags (ESTs) library was found to be efficient at isolating simple sequence repeats (SSRs) loci in the non-indigenous marine mollusc Crepidula fornicata. Twelve SSR loci were developed for routine genotyping. Cross-species amplification to 8 other Crepidula species showed that the 12 loci are highly specific for C. fornicata. Mendelian inheritance was shown for 11 of them (1 being monomorphic in the analyzed offspring array). The genetic diversity for 88 adults was found to be variable across the 12 loci (2-40 alleles, expected heterozygosity between 0.023 and 0.898) with a high overall exclusion probability of 0.99. The degree of genetic polymorphism found here is similar to that shown for 7 anonymous SSRs previously developed and here used on the same samples. This set of 12 specific loci is relevant to perform reliable population and relatedness analyses in Crepidula fornicata.

Stress ecology in fucus: abiotic, biotic and genetic interactions

Stress regimes defined as the synchronous or sequential action of abiotic and biotic stresses determine the performance and distribution of species. The natural patterns of stress to which species are more or less well adapted have recently started to shift and alter under the influence of global change. This was the motivation to review our knowledge on the stress ecology of a benthic key player, the macroalgal genus Fucus. We first provide a comprehensive review of the genus as an ecological model including what is currently known about the major lineages of Fucus species with respect to hybridization, ecotypic differentiation and speciation; as well as life history, population structure and geographic distribution. We then review our current understanding of both extrinsic (abiotic/biotic) and intrinsic (genetic) stress(es) on Fucus species and how they interact with each other. It is concluded that (i) interactive stress effects appear to be equally distributed over additive, antagonistic and synergistic categories at the level of single experiments, but are predominantly additive when averaged over all studies in a meta-analysis of 41 experiments; (ii) juvenile and adult responses to stress frequently differ and (iii) several species or particular populations of Fucus may be relatively unaffected by climate change as a consequence of pre-adapted ecotypes that collectively express wide physiological tolerences. Future research on Fucus should (i) include additional species, (ii) include marginal populations as models for responses to environmental stress; (iii) assess a wider range of stress combinations, including their temporal fluctuations; (iv) better differentiate between stress sensitivity of juvenile versus adult stages; (v) include a functional genomic component in order to better integrate Fucus' ecological and evolutionary responses to stress regimes and (vi) utilize a multivariate modelling approach in order to develop and understand interaction networks.