Alexander Jueterbock

Calculations of population differentiation based on GST and D: forget GST but not all of statistics!

G ST‐values and its relatives (FST) belong to the most used parameters to define genetic differences between populations. Originally, they were developed for allozymes with very low number of alleles. Using highly polymorphic microsatellite markers it was often puzzling that GST‐values were very low but statistically significant. In their papers, Jost (2008) and Hedrick (2005) explained that GST‐values do not show genetic differentiation, and Jost suggested calculating D‐values instead. Theoretical mathematical considerations are often difficult to follow; therefore, we chose an applied approach comparing two artificial populations with different number of alleles at equal frequencies and known genetic divergence. Our results show that even for more than one allele per population GST‐values do not calculate population differentiation correctly; in contrast, D‐values do reflect the genetic differentiation indicating that data based on GST‐values need to be re‐evaluated. In our approach, statistical evaluations remained similar. We provide information about the impact of different sample sizes on D‐values in relation to number of alleles and genetic divergence.

The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea

Seagrasses colonized the sea on at least three independent occasions to form the basis of one of the most productive and widespread coastal ecosystems on the planet. Here we report the genome of Zostera marina (L.), the first, to our knowledge, marine angiosperm to be fully sequenced. This reveals unique insights into the genomic losses and gains involved in achieving the structural and physiological adaptations required for its marine lifestyle, arguably the most severe habitat shift ever accomplished by flowering plants. Key angiosperm innovations that were lost include the entire repertoire of stomatal genes, genes involved in the synthesis of terpenoids and ethylene signalling, and genes for ultraviolet protection and phytochromes for far-red sensing. Seagrasses have also regained functions enabling them to adjust to full salinity. Their cell walls contain all of the polysaccharides typical of land plants, but also contain polyanionic, low-methylated pectins and sulfated galactans, a feature shared with the cell walls of all macroalgae and that is important for ion homoeostasis, nutrient uptake and O2/CO2 exchange through leaf epidermal cells. The Z. marina genome resource will markedly advance a wide range of functional ecological studies from adaptation of marine ecosystems under climate warming, to unravelling the mechanisms of osmoregulation under high salinities that may further inform our understanding of the evolution of salt tolerance in crop plants.

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.

Improving transferability of introduced species' distribution models: new tools to forecast the spread of a highly invasive seaweed.

The utility of species distribution models for applications in invasion and global change biology is critically dependent on their transferability between regions or points in time, respectively. We introduce two methods that aim to improve the transferability of presence-only models: density-based occurrence thinning and performance-based predictor selection. We evaluate the effect of these methods along with the impact of the choice of model complexity and geographic background on the transferability of a species distribution model between geographic regions. Our multifactorial experiment focuses on the notorious invasive seaweed Caulerpacylindracea (previously Caulerpa racemosa var. cylindracea ) and uses Maxent, a commonly used presence-only modeling technique. We show that model transferability is markedly improved by appropriate predictor selection, with occurrence thinning, model complexity and background choice having relatively minor effects. The data shows that, if available, occurrence records from the native and invaded regions should be combined as this leads to models with high predictive power while reducing the sensitivity to choices made in the modeling process. The inferred distribution model of Caulerpacylindracea shows the potential for this species to further spread along the coasts of Western Europe, western Africa and the south coast of Australia.

Genomic divergence between the migratory and stationary ecotypes of Atlantic cod

Atlantic cod displays a range of phenotypic and genotypic variations, which includes the differentiation into coastal stationary and offshore migratory types of cod that co‐occur in several parts of its distribution range and are often sympatric on the spawning grounds. Differentiation of these ecotypes may involve both historical separation and adaptation to ecologically distinct environments, the genetic basis of which is now beginning to be unravelled. Genomic analyses based on recent sequencing advances are able to document genomic divergence in more detail and may facilitate the exploration of causes and consequences of genome‐wide patterns. We examined genomic divergence between the stationary and migratory types of cod in the Northeast Atlantic, using next‐generation sequencing of pooled DNA from each of two population samples. Sequence data was mapped to the published cod genome sequence, arranged in more than 6000 scaffolds (611 Mb). We identified 25 divergent scaffolds (26 Mb) with a higher than average gene density, against a backdrop of overall moderate genomic differentiation. Previous findings of localized genomic divergence in three linkage groups were confirmed, including a large (15 Mb) genomic region, which seems to be uniquely involved in the divergence of migratory and stationary cod. The results of the pooled sequencing approach support and extend recent findings based on single‐nucleotide polymorphism markers and suggest a high degree of reproductive isolation between stationary and migratory cod in the North‐east Atlantic.

Phylogeographic differentiation versus transcriptomic adaptation to warm temperatures in Zostera marina, a globally important seagrass

Populations distributed across a broad thermal cline are instrumental in addressing adaptation to increasing temperatures under global warming. Using a space‐for‐time substitution design, we tested for parallel adaptation to warm temperatures along two independent thermal clines in Zostera marina, the most widely distributed seagrass in the temperate Northern Hemisphere. A North–South pair of populations was sampled along the European and North American coasts and exposed to a simulated heatwave in a common‐garden mesocosm. Transcriptomic responses under control, heat stress and recovery were recorded in 99 RNAseq libraries with ~13 000 uniquely annotated, expressed genes. We corrected for phylogenetic differentiation among populations to discriminate neutral from adaptive differentiation. The two southern populations recovered faster from heat stress and showed parallel transcriptomic differentiation, as compared with northern populations. Among 2389 differentially expressed genes, 21 exceeded neutral expectations and were likely involved in parallel adaptation to warm temperatures. However, the strongest differentiation following phylogenetic correction was between the three Atlantic populations and the Mediterranean population with 128 of 4711 differentially expressed genes exceeding neutral expectations. Although adaptation to warm temperatures is expected to reduce sensitivity to heatwaves, the continued resistance of seagrass to further anthropogenic stresses may be impaired by heat‐induced downregulation of genes related to photosynthesis, pathogen defence and stress tolerance.

Variation in thermal stress response in two populations of the brown seaweed, Fucus distichus, from the Arctic and subarctic intertidal

It is unclear whether intertidal organisms are ‘preadapted’ to cope with the increase of temperature and temperature variability or if they are currently at their thermal tolerance limits. To address the dichotomy, we focused on an important ecosystem engineer of the Arctic intertidal rocky shores, the seaweed Fucus distichus and investigated thermal stress responses of two populations from different temperature regimes (Svalbard and Kirkenes, Norway). Thermal stress responses at 20°C, 24°C and 28°C were assessed by measuring photosynthetic performance and expression of heat shock protein (HSP) genes (shsp, hsp90 and hsp70). We detected population-specific responses between the two populations of F. distichus, as the Svalbard population revealed a smaller decrease in photosynthesis performance but a greater activation of molecular defence mechanisms (indicated by a wider repertoire of HSP genes and their stronger upregulation) compared with the Kirkenes population. Although the temperatures used in our study exceed temperatures encountered by F. distichus at the study sites, we believe response to these temperatures may serve as a proxy for the species’ potential to respond to climate-related stresses.

Pre-zygotic isolation in the macroalgal genus Fucus from four contact zones spanning 100-10 000 years: A tale of reinforcement?

Hybrid zones provide an ideal natural experiment to study the selective forces driving evolution of reproductive barriers and speciation. If hybrid offspring are less fit than the parental species, pre-zygotic isolating barriers can evolve and strengthen in response to selection against the hybrids (reinforcement). Four contact zones between the intertidal macroalgae Fucus serratus (Fs) and Fucus distichus (Fd), characterized by varying times of sympatry and order of species introduction provide an opportunity to investigate reinforcement. We examined patterns of hybridization and reproductive isolation between Fs and Fd in: (i) northern Norway (consisting of two natural sites, 10 000 years old), (ii) the Kattegat near Denmark (Fd introduced, nineteenth century) and (iii) Iceland (Fs introduced, nineteenth century). Using 10 microsatellites and chloroplast DNA, we showed that hybridization and introgression decreased with increasing duration of sympatry. The two younger contact zones revealed 13 and 24% hybrids and several F1 individuals, in contrast to the older contact zone with 2–3% hybrids and an absence of F1s. Cross-fertilization experiments revealed that the reduction in hybridization in the oldest zone is consistent with increased gametic incompatibility.

Decadal stability in genetic variation and structure in the intertidal seaweed Fucus serratus (Heterokontophyta: Fucaceae)

BackgroundThe spatial distribution of genetic diversity and structure has important implications for conservation as it reveals a species’ strong and weak points with regard to stability and evolutionary capacity. Temporal genetic stability is rarely tested in marine species other than commercially important fishes, but is crucial for the utility of temporal snapshots in conservation management. High and stable diversity can help to mitigate the predicted northward range shift of seaweeds under the impact of climate change. Given the key ecological role of fucoid seaweeds along rocky shores, the positive effect of genetic diversity may reach beyond the species level to stabilize the entire intertidal ecosystem along the temperate North Atlantic. In this study, we estimated the effective population size, as well as temporal changes in genetic structure and diversity of the seaweed F. serratus using 22 microsatellite markers. Samples were taken across latitudes and a range of temperature regimes at seven locations with decadal sampling (2000 and 2010).ResultsAcross latitudes, genetic structure and diversity remained stable over 5–10 generations. Stable small-scale structure enhanced regional diversity throughout the species’ range. In accordance with its biogeographic history, effective population size and diversity peaked in the species’ mid-range in Brittany (France), and declined towards its leading and trailing edge to the north and south. At the species’ southern edge, multi-locus-heterozygosity displayed a strong decline from 1999 to 2010.ConclusionTemporally stable genetic structure over small spatial scales is a potential driver for local adaptation and species radiation in the genus Fucus. Survival and adaptation of the low-diversity leading edge of F. serratus may be enhanced by regional gene flow and ‘surfing’ of favorable mutations or impaired by the accumulation of deleterious mutations. Our results have clear implications for the conservation of F. serratus at its genetically unique southern edge in Northwest Iberia, where increasing temperatures are likely the major cause for the decline not only of F. serratus, but also other intertidal and subtidal macroalgae. We expect that F. serratus will disappear from Northwest Iberia by 2100 if genetic rescue is not induced by the influx of genetic variation from Brittany.

Table S4 (Excel): Annotations of mapped reads.

ZosmaID: gene ID based on the Zostera marina genome annotation v2.1 from the ORCAE database, GenBank Accession: LFYR00000000); mRNAID: the mapped sequence ID; name: the gene description; annotation.source: the source of annotation (inference from homology or top BLAST hit); and GO.terms: the associated Gene Ontology terms.