The archaeome in metaorganism research, with a focus on marine models and their bacteria-archaea interactions

Metaorganism research contributes substantially to our understanding of the interaction between microbes and their hosts, as well as their co-evolution. Most research is currently focused on the bacterial community, while archaea often remain at the sidelines of metaorganism-related research. Here, we describe the archaeome of a total of eleven classical and emerging multicellular model organisms across the phylogenetic tree of life. To determine the microbial community composition of each host, we utilized a combination of archaea and bacteria-specific 16S rRNA gene amplicons. Members of the two prokaryotic domains were described regarding their community composition, diversity, and richness in each multicellular host. Moreover, association with specific hosts and possible interaction partners between the bacterial and archaeal communities were determined for the marine models. Our data show that the archaeome in marine hosts predominantly consists of Nitrosopumilaceae and Nanoarchaeota, which represent keystone taxa among the porifera. The presence of an archaeome in the terrestrial hosts varies substantially. With respect to abundant archaeal taxa, they harbor a higher proportion of methanoarchaea over the aquatic environment. We find that the archaeal community is much less diverse than its bacterial counterpart. Archaeal amplicon sequence variants are usually host-specific, suggesting adaptation through co-evolution with the host. While bacterial richness was higher in the aquatic than the terrestrial hosts, a significant difference in diversity and richness between these groups could not be observed in the archaeal dataset. Our data show a large proportion of unclassifiable archaeal taxa, highlighting the need for improved cultivation efforts and expanded databases.

Seagrass genomes reveal ancient polyploidy and adaptations to the marine environment

We present chromosome-level genome assemblies from representative species of three independently evolved seagrass lineages: Posidonia oceanica, Cymodocea nodosa, Thalassia testudinum and Zostera marina. We also include a draft genome of Potamogeton acutifolius, belonging to a freshwater sister lineage to Zosteraceae. All seagrass species share an ancient whole-genome triplication, while additional whole-genome duplications were uncovered for C. nodosa, Z. marina and P. acutifolius. Comparative analysis of selected gene families suggests that the transition from submerged-freshwater to submerged-marine environments mainly involved fine-tuning of multiple processes (such as osmoregulation, salinity, light capture, carbon acquisition and temperature) that all had to happen in parallel, probably explaining why adaptation to a marine lifestyle has been exceedingly rare. Major gene losses related to stomata, volatiles, defence and lignification are probably a consequence of the return to the sea rather than the cause of it. These new genomes will accelerate functional studies and solutions, as continuing losses of the 'savannahs of the sea' are of major concern in times of climate change and loss of biodiversity.

Origin matters: Using a local reference genome improves measures in population genomics

Genome sequencing enables answering fundamental questions about the genetic basis of adaptation, population structure and epigenetic mechanisms. Yet, we usually need a suitable reference genome for mapping population-level resequencing data. In some model systems, multiple reference genomes are available, giving the challenging task of determining which reference genome best suits the data. Here, we compared the use of two different reference genomes for the three-spined stickleback (Gasterosteus aculeatus), one novel genome derived from a European gynogenetic individual and the published reference genome of a North American individual. Specifically, we investigated the impact of using a local reference versus one generated from a distinct lineage on several common population genomics analyses. Through mapping genome resequencing data of 60 sticklebacks from across Europe and North America, we demonstrate that genetic distance among samples and the reference genomes impacts downstream analyses. Using a local reference genome increased mapping efficiency and genotyping accuracy, effectively retaining more and better data. Despite comparable distributions of the metrics generated across the genome using SNP data (i.e. π, Tajima's D and FST ), window-based statistics using different references resulted in different outlier genes and enriched gene functions. A marker-based analysis of DNA methylation distributions had a comparably high overlap in outlier genes and functions, yet with distinct differences depending on the reference genome. Overall, our results highlight how using a local reference genome decreases reference bias to increase confidence in downstream analyses of the data. Such results have significant implications in all reference-genome-based population genomic analyses.

An evolutionary epigenetic clock in plants

Molecular clocks are the basis for dating the divergence between lineages over macroevolutionary timescales (~105 to 108 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes display a clocklike behavior. This "epimutation clock" is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation clocks recapitulate known topologies and branching times of intraspecies phylogenetic trees in the self-fertilizing plant Arabidopsis thaliana and the clonal seagrass Zostera marina, which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

Worldwide Population Genomics Reveal Long-Term Stability of the Mitochondrial Genome Architecture in a Keystone Marine Plant

Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

Population genetic differentiation of the ubiquitous brooding coral Pocillopora acuta along Phuket Island reefs in the Andaman Sea, Thailand

BACKGROUND

The widespread Indo-Pacific coral species Pocillopora acuta Lamarck, 1816 displays varying levels of asexual versus sexual reproduction, with strong repercussions on genetic diversity, connectivity and genetic structuring within and among populations. For many geographic regions, baseline information on genetic diversity is still lacking, particularly in the Andaman Sea. The region suffered a massive heat-induced bleaching event in 2010 with high coral cover loss of branching coral species such as P. acuta. A subsequent bleaching in 2016, however, revealed a mild bleaching response in pocilloporids compared to other coral taxa in the region, suggesting that rare, heat tolerant genotypes had been selected by the 2010 bleaching event. In order to test whether this potential 'evolutionary rescue' event has led to a low genetic diversity, we conducted a population genetic survey covering a total of nine different P. acuta populations (336 individuals) along a 50 km coastal stretch around Phuket Island, Thailand. We used six microsatellite markers to assess genotypic diversity and to determine the prevalent mode of reproduction (i.e. sexual or asexual recruitment).

RESULTS

In contrast to other Indian Ocean P. acuta populations, the majority of corals in this study adopted a sexual reproduction mode (75% across all populations). At the same time, substantial regional gene flow was observed around Phuket Island with strong genetic differentiation as indicated by three genetic clusters that were separated by only a few kilometers. Patterns of isolation by distance over 0.7 - 40 km suggest small-scale genetic barriers, such as changing currents throughout each monsoonal season, potentially contributing to locally restricted dispersal of P. acuta larvae.

CONCLUSIONS

The occurrence of distinct genetic clusters within short coastal stretches suggests that the 2010 bleaching event has not led to extreme genetic impoverishment. While more in-depth genomic analyses are necessary to investigate changes in genetic diversity following extreme bleaching events, our results will help guide conservation efforts to maintain genetic diversity of a coral species that likely will be dominant in future, warmer Andaman Sea reefs.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Detecting clonemate pairs in multicellular diploid clonal species based on a shared heterozygosity index

Clonal reproduction, the formation of nearly identical individuals via mitosis in the absence of genetic recombination, is a very common reproductive mode across plants, fungi and animals. To detect clonal genetic structure, genetic similarity indices based on shared alleles are widely used, such as the Jaccard index, or identity by state. Here we propose a new pairwise genetic similarity index, the SH index, based on segregating genetic marker loci (typically single nucleotide polymorphisms) that are identically heterozygous for pairs of samples (N_SH ). To test our method, we analyse two old seagrass clones (Posidonia australis, estimated to be around 8500 years old; Zostera marina, >750 years old) along with two young Z. marina clones of known age (17 years old). We show that focusing on shared heterozygosity amplifies the power to distinguish sample pairs belonging to different clones compared to methods focusing on all shared alleles. Our proposed workflow can successfully detect clonemates at a location dominated by a single clone. When the collected samples involve two or more clones, the SH index shows a clear gap between clonemate pairs and interclone sample pairs. Ideally N_SH should be on the order of approximately ≥3000, a number easily achievable via restriction-site associated DNA (RAD) sequencing or whole-genome resequencing. Another potential application of the SH index is to detect possible parent-descendant pairs under selfing. Our proposed workflow takes advantage of the availability of the larger number of genetic markers in the genomic era, and improves the ability to distinguish clonemates from nonclonemates in multicellular diploid clonal species.

A trait-based framework for seagrass ecology: Trends and prospects

In the last three decades, quantitative approaches that rely on organism traits instead of taxonomy have advanced different fields of ecological research through establishing the mechanistic links between environmental drivers, functional traits, and ecosystem functions. A research subfield where trait-based approaches have been frequently used but poorly synthesized is the ecology of seagrasses; marine angiosperms that colonized the ocean 100M YA and today make up productive yet threatened coastal ecosystems globally. Here, we compiled a comprehensive trait-based response-effect framework (TBF) which builds on previous concepts and ideas, including the use of traits for the study of community assembly processes, from dispersal and response to abiotic and biotic factors, to ecosystem function and service provision. We then apply this framework to the global seagrass literature, using a systematic review to identify the strengths, gaps, and opportunities of the field. Seagrass trait research has mostly focused on the effect of environmental drivers on traits, i.e., "environmental filtering" (72%), whereas links between traits and functions are less common (26.9%). Despite the richness of trait-based data available, concepts related to TBFs are rare in the seagrass literature (15% of studies), including the relative importance of neutral and niche assembly processes, or the influence of trait dominance or complementarity in ecosystem function provision. These knowledge gaps indicate ample potential for further research, highlighting the need to understand the links between the unique traits of seagrasses and the ecosystem services they provide.

Improved chromosome-level genome assembly and annotation of the seagrass, Zostera marina (eelgrass)

Background: Seagrasses (Alismatales) are the only fully marine angiosperms.  Zostera marina (eelgrass) plays a crucial role in the functioning of coastal marine ecosystems and global carbon sequestration. It is the most widely studied seagrass and has become a marine model system for exploring adaptation under rapid climate change. The original draft genome (v.1.0) of the seagrass  Z. marina (L.) was based on a combination of Illumina mate-pair libraries and fosmid-ends. A total of 25.55 Gb of Illumina and 0.14 Gb of Sanger sequence was obtained representing 47.7× genomic coverage. The assembly resulted in ~2000 unordered scaffolds (L50 of 486 Kb), a final genome assembly size of 203MB, 20,450 protein coding genes and 63% TE content. Here, we present an upgraded chromosome-scale genome assembly and compare v.1.0 and the new v.3.1, reconfirming previous results from Olsen et al. (2016), as well as pointing out new findings.  

Methods: The same high molecular weight DNA used in the original sequencing of the Finnish clone was used. A high-quality reference genome was assembled with the MECAT assembly pipeline combining PacBio long-read sequencing and Hi-C scaffolding. 

Results: In total, 75.97 Gb PacBio data was produced. The final assembly comprises six pseudo-chromosomes and 304 unanchored scaffolds with a total length of 260.5Mb and an N50 of 34.6 MB, showing high contiguity and few gaps (~0.5%). 21,483 protein-encoding genes are annotated in this assembly, of which 20,665 (96.2%) obtained at least one functional assignment based on similarity to known proteins. 

Conclusions: As an important marine angiosperm, the improved  Z. marina genome assembly will further assist evolutionary, ecological, and comparative genomics at the chromosome level. The new genome assembly will further our understanding into the structural and physiological adaptations from land to marine life.

Phenotypic plasticity under rapid global changes: The intrinsic force for future seagrasses survival

Coastal oceans are particularly affected by rapid and extreme environmental changes with dramatic consequences for the entire ecosystem. Seagrasses are key ecosystem engineering or foundation species supporting diverse and productive ecosystems along the coastline that are particularly susceptible to fast environmental changes. In this context, the analysis of phenotypic plasticity could reveal important insights into seagrasses persistence, as it represents an individual property that allows species' phenotypes to accommodate and react to fast environmental changes and stress. Many studies have provided different definitions of plasticity and related processes (acclimation and adaptation) resulting in a variety of associated terminology. Here, we review different ways to define phenotypic plasticity with particular reference to seagrass responses to single and multiple stressors. We relate plasticity to the shape of reaction norms, resulting from genotype by environment interactions, and examine its role in the presence of environmental shifts. The potential role of genetic and epigenetic changes in underlying seagrasses plasticity in face of environmental changes is also discussed. Different approaches aimed to assess local acclimation and adaptation in seagrasses are explored, explaining strengths and weaknesses based on the main results obtained from the most recent literature. We conclude that the implemented experimental approaches, whether performed with controlled or field experiments, provide new insights to explore the basis of plasticity in seagrasses. However, an improvement of molecular analysis and the application of multi-factorial experiments are required to better explore genetic and epigenetic adjustments to rapid environmental shifts. These considerations revealed the potential for selecting the best phenotypes to promote assisted evolution with fundamental implications on restoration and preservation efforts.

Deep-sea predator niche segregation revealed by combined cetacean biologging and eDNA analysis of cephalopod prey

Fundamental insight on predator-prey dynamics in the deep sea is hampered by a lack of combined data on hunting behavior and prey spectra. Deep-sea niche segregation may evolve when predators target specific prey communities, but this hypothesis remains untested. We combined environmental DNA (eDNA) metabarcoding with biologging to assess cephalopod community composition in the deep-sea foraging habitat of two top predator cetaceans. Risso's dolphin and Cuvier's beaked whale selectively targeted distinct epi/meso- and bathypelagic foraging zones, holding eDNA of 39 cephalopod taxa, including 22 known prey. Contrary to expectation, extensive taxonomic overlap in prey spectra between foraging zones indicated that predator niche segregation was not driven by prey community composition alone. Instead, intraspecific prey spectrum differences may drive differentiation for hunting fewer, more calorific, mature cephalopods in deeper waters. The novel combination of methods presented here holds great promise to disclose elusive deep-sea predator-prey systems, aiding in their protection.

A novel metabarcoding primer pair for environmental DNA analysis of Cephalopoda (Mollusca) targeting the nuclear 18S rRNA region

Cephalopods are pivotal components of marine food webs, but biodiversity studies are hampered by challenges to sample these agile marine molluscs. Metabarcoding of environmental DNA (eDNA) is a potentially powerful technique to study oceanic cephalopod biodiversity and distribution but has not been applied thus far. We present a novel universal primer pair for metabarcoding cephalopods from eDNA, Ceph18S (Forward: 5'-CGC GGC GCT ACA TAT TAG AC-3', Reverse: 5'-GCA CTT AAC CGA CCG TCG AC-3'). The primer pair targets the hypervariable region V2 of the nuclear 18S rRNA gene and amplifies a relatively short target sequence of approximately 200 bp in order to allow the amplification of degraded DNA. In silico tests on a reference database and empirical tests on DNA extracts from cephalopod tissue estimate that 44-66% of cephalopod species, corresponding to about 310-460 species, can be amplified and identified with this primer pair. A multi-marker approach with the novel Ceph18S and two previously published cephalopod mitochondrial 16S rRNA primer sets targeting the same region (Jarman et al. 2006 Mol. Ecol. Notes. 6, 268-271; Peters et al. 2015 Mar. Ecol. 36, 1428-1439) is estimated to amplify and identify 89% of all cephalopod species, of which an estimated 19% can only be identified by Ceph18S. All sequences obtained with Ceph18S were submitted to GenBank, resulting in new 18S rRNA sequences for 13 cephalopod taxa.

From ecosystems to socio-economic benefits: A systematic review of coastal ecosystem services in the Baltic Sea

Seagrass meadows, algal forests and mussel beds are widely regarded as foundation species that support communities providing valuable ecosystem services in many coastal regions; however, quantitative evidence of the relationship is scarce. Using the Baltic Sea as a case study, a region of significant socio-economic importance in the northern hemisphere, we systematically synthesized the primary literature and summarized the current knowledge on ecosystem services derived from seagrass, macroalgae, and mussels (see animated video summary of the manuscript: Video abstract). We found 1740 individual ecosystem service records (ESR), 61% of which were related to macroalgae, 26% to mussel beds and 13% to seagrass meadows. The most frequently reported ecosystem services were raw material (533 ESR), habitat provision (262 ESR) and regulation of pollutants (215 ESR). Toxins (356 ESR) and nutrients (302 ESR) were the most well-documented pressures to services provided by coastal ecosystems. Next, we assessed the current state of knowledge as well as knowledge transfer of ecosystem services to policies through natural, social, human and economic dimensions, using a systematic scoring tool, the Eco-GAME matrix. We found good quantitative information about how ecosystems generated the service but almost no knowledge of how they translate into socio-economic benefits (8 out of 657 papers, 1.2%). While we are aware that research on Baltic Sea socio-economic benefits does exist, the link with ecosystems providing the service is mostly missing. To close this knowledge gap, we need a better analytical framework that is capable of directly linking existing quantitative information about ecosystem service generation with human benefit.

Modeling eelgrass spatial response to nutrient abatement measures in a changing climate

For many coastal areas including the Baltic Sea, ambitious nutrient abatement goals have been set to curb eutrophication, but benefits of such measures were normally not studied in light of anticipated climate change. To project the likely responses of nutrient abatement on eelgrass (Zostera marina), we coupled a species distribution model with a biogeochemical model, obtaining future water turbidity, and a wave model for predicting the future hydrodynamics in the coastal area. Using this, eelgrass distribution was modeled for different combinations of nutrient scenarios and future wind fields. We are the first to demonstrate that while under a business as usual scenario overall eelgrass area will not recover, nutrient reductions that fulfill the Helsinki Commission's Baltic Sea Action Plan (BSAP) are likely to lead to a substantial areal expansion of eelgrass coverage, primarily at the current distribution's lower depth limits, thereby overcompensating losses in shallow areas caused by a stormier climate.

Differences in the microbiota of native and non-indigenous gelatinous zooplankton organisms in a low saline environment

The translocation of non-indigenous species (NIS) around the world, especially in marine systems, is increasingly being recognized as a matter of concern. Species translocations have been shown to lead to wide ranging changes in food web structure and functioning. In addition to the direct effects of NIS, they could facilitate the accumulation or translocation of bacteria as part of their microbiomes. The Baltic Sea harbours many non-indigenous species, with most recent detection of the jellyfish Blackfordia virginica and the comb jelly Mnemiopsis leidyi in the low saline southwestern Baltic Sea. In this study, we used a multidisciplinary approach and investigated three gelatinous zooplankton species that co-occur in the same environment and feed on similar zooplankton food sources but show different histories of origin. The aim was to conduct a comparative microbiome analysis of indigenous and non-indigenous gelatinous zooplankton species in the low-saline southwestern Baltic Sea. Next-generation 16S rRNA marker gene sequencing of the V1/V2 region was employed to study the bacterial microbiome compositions. All tested species showed significant differences in their microbiome compositions (one way ANOSIM, R = 1, P < 0.008) with dissimilarities ranging from 85 to 92%. The indigenous jellyfish Aurelia aurita showed the highest bacterial operational taxonomic unit (OTU) richness. The overall differentiation between microbiomes was driven by eight indicator OTUs, which included Mycoplasma and Vibrio species. These bacteria can be problematic, as they include known pathogenic strains that are relevant to human health and aquaculture activities. Our results suggest that the impact assessment of NIS should consider potential pathogenic bacteria, enriched in the environment due to invasion, as potential risks to aquaculture activities.

Cultivable microbiota associated with Aurelia aurita and Mnemiopsis leidyi

The associated microbiota of marine invertebrates plays an important role to the host in relation to fitness, health, and homeostasis. Cooperative and competitive interactions between bacteria, due to release of, for example, antibacterial substances and quorum sensing (QS)/quorum quenching (QQ) molecules, ultimately affect the establishment and dynamics of the associated microbial community. Aiming to address interspecies competition of cultivable microbes associated with emerging model species of the basal animal phyla Cnidaria (Aurelia aurita) and Ctenophora (Mnemiopsis leidyi), we performed a classical isolation approach. Overall, 84 bacteria were isolated from A. aurita medusae and polyps, 64 bacteria from M. leidyi, and 83 bacteria from ambient seawater, followed by taxonomically classification by 16S rRNA gene analysis. The results show that A. aurita and M. leidyi harbor a cultivable core microbiome consisting of typical marine ubiquitous bacteria also found in the ambient seawater. However, several bacteria were restricted to one host suggesting host‐specific microbial community patterns. Interbacterial interactions were assessed by (a) a growth inhibition assay and (b) QS interference screening assay. Out of 231 isolates, 4 bacterial isolates inhibited growth of 17 isolates on agar plates. Moreover, 121 of the 231 isolates showed QS‐interfering activities. They interfered with the acyl‐homoserine lactone (AHL)‐based communication, of which 21 showed simultaneous interference with autoinducer 2. Overall, this study provides insights into the cultivable part of the microbiota associated with two environmentally important marine non‐model organisms and into interbacterial interactions, which are most likely considerably involved in shaping a healthy and resilient microbiota.

Somatic genetic drift and multilevel selection in a clonal seagrass

All multicellular organisms are genetic mosaics owing to somatic mutations. The accumulation of somatic genetic variation in clonal species undergoing asexual (or clonal) reproduction may lead to phenotypic heterogeneity among autonomous modules (termed ramets). However, the abundance and dynamics of somatic genetic variation under clonal reproduction remain poorly understood. Here we show that branching events in a seagrass (Zostera marina) clone or genet lead to population bottlenecks of tissue that result in the evolution of genetically differentiated ramets in a process of somatic genetic drift. By studying inter-ramet somatic genetic variation, we uncovered thousands of single nucleotide polymorphisms that segregated among ramets. Ultra-deep resequencing of single ramets revealed that the strength of purifying selection on mosaic genetic variation was greater within than among ramets. Our study provides evidence for multiple levels of selection during the evolution of seagrass genets. Somatic genetic drift during clonal propagation leads to the emergence of genetically unique modules that constitute an elementary level of selection and individuality in long-lived clonal species.

Genome-Wide Genotype-Expression Relationships Reveal Both Copy Number and Single Nucleotide Differentiation Contribute to Differential Gene Expression between Stickleback Ecotypes

Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.

Two different epigenetic information channels in wild three-spined sticklebacks are involved in salinity adaptation

Epigenetic inheritance has been proposed to contribute to adaptation and acclimation via two information channels: (i) inducible epigenetic marks that enable transgenerational plasticity and (ii) noninducible epigenetic marks resulting from random epimutations shaped by selection. We studied both postulated channels by sequencing methylomes and genomes of Baltic three-spined sticklebacks (Gasterosteus aculeatus) along a salinity cline. Wild populations differing in salinity tolerance revealed differential methylation (pop-DMS) at genes enriched for osmoregulatory processes. A two-generation experiment demonstrated that 62% of these pop-DMS were noninducible by salinity manipulation, suggesting that they are the result of either direct selection or associated genomic divergence at cis- or trans-regulatory sites. Two-thirds of the remaining inducible pop-DMS increased in similarity to patterns detected in wild populations from corresponding salinities. The level of similarity accentuated over consecutive generations, indicating a mechanism of transgenerational plasticity. While we can attribute natural DNA methylation patterns to the two information channels, their interplay with genomic variation in salinity adaptation is still unresolved.

Comparative analysis of amplicon and metagenomic sequencing methods reveals key features in the evolution of animal metaorganisms

BACKGROUND

The interplay between hosts and their associated microbiome is now recognized as a fundamental basis of the ecology, evolution, and development of both players. These interdependencies inspired a new view of multicellular organisms as "metaorganisms." The goal of the Collaborative Research Center "Origin and Function of Metaorganisms" is to understand why and how microbial communities form long-term associations with hosts from diverse taxonomic groups, ranging from sponges to humans in addition to plants.

METHODS

In order to optimize the choice of analysis procedures, which may differ according to the host organism and question at hand, we systematically compared the two main technical approaches for profiling microbial communities, 16S rRNA gene amplicon and metagenomic shotgun sequencing across our panel of ten host taxa. This includes two commonly used 16S rRNA gene regions and two amplification procedures, thus totaling five different microbial profiles per host sample.

CONCLUSION

While 16S rRNA gene-based analyses are subject to much skepticism, we demonstrate that many aspects of bacterial community characterization are consistent across methods. The resulting insight facilitates the selection of appropriate methods across a wide range of host taxa. Overall, we recommend single- over multi-step amplification procedures, and although exceptions and trade-offs exist, the V3 V4 over the V1 V2 region of the 16S rRNA gene. Finally, by contrasting taxonomic and functional profiles and performing phylogenetic analysis, we provide important and novel insight into broad evolutionary patterns among metaorganisms, whereby the transition of animals from an aquatic to a terrestrial habitat marks a major event in the evolution of host-associated microbial composition.

Ecological-economic sustainability of the Baltic cod fisheries under ocean warming and acidification

Human-induced climate change such as ocean warming and acidification, threatens marine ecosystems and associated fisheries. In the Western Baltic cod stock socio-ecological links are particularly important, with many relying on cod for their livelihoods. A series of recent experiments revealed that cod populations are negatively affected by climate change, but an ecological-economic assessment of the combined effects, and advice on optimal adaptive management are still missing. For Western Baltic cod, the increase in larval mortality due to ocean acidification has experimentally been quantified. Time-series analysis allows calculating the temperature effect on recruitment. Here, we include both processes in a stock-recruitment relationship, which is part of an ecological-economic optimization model. The goal was to quantify the effects of climate change on the triple bottom line (ecological, economic, social) of the Western Baltic cod fishery. Ocean warming has an overall negative effect on cod recruitment in the Baltic. Optimal management would react by lowering fishing mortality with increasing temperature, to create a buffer against climate change impacts. The negative effects cannot be fully compensated, but even at 3 °C warming above the 2014 level, a reduced but viable fishery would be possible. However, when accounting for combined effects of ocean warming and acidification, even optimal fisheries management cannot adapt to changes beyond a warming of +1.5° above the current level. Our results highlight the need for multi-factorial climate change research, in order to provide the best available, most realistic, and precautionary advice for conservation of exploited species as well as their connected socio-economic systems.

Genetic analyses reveal complex dynamics within a marine fish management area

Genetic data have great potential for improving fisheries management by identifying the fundamental management units-that is, the biological populations-and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

Divergent responses of Atlantic cod to ocean acidification and food limitation

In order to understand the effect of global change on marine fishes, it is imperative to quantify the effects on fundamental parameters such as survival and growth. Larval survival and recruitment of the Atlantic cod (Gadus morhua) were found to be heavily impaired by end-of-century levels of ocean acidification. Here, we analysed larval growth among 35-36 days old surviving larvae, along with organ development and ossification of the skeleton. We combined CO₂ treatments (ambient: 503 µatm, elevated: 1,179 µatm) with food availability in order to evaluate the effect of energy limitation in addition to the ocean acidification stressor. As expected, larval size (as a proxy for growth) and skeletogenesis were positively affected by high food availability. We found significant interactions between acidification and food availability. Larvae fed ad libitum showed little difference in growth and skeletogenesis due to the CO₂ treatment. Larvae under energy limitation were significantly larger and had further developed skeletal structures in the elevated CO₂ treatment compared to the ambient CO₂ treatment. However, the elevated CO₂ group revealed impairments in critically important organs, such as the liver, and had comparatively smaller functional gills indicating a mismatch between size and function. It is therefore likely that individual larvae that had survived acidification treatments will suffer from impairments later during ontogeny. Our study highlights important allocation trade-off between growth and organ development, which is critically important to interpret acidification effects on early life stages of fish.

Transgenerational plasticity and selection shape the adaptive potential of sticklebacks to salinity change

In marine climate change research, salinity shifts have been widely overlooked. While widespread desalination effects are expected in higher latitudes, salinity is predicted to increase closer to the equator. We took advantage of the steep salinity gradient of the Baltic Sea as a space-for-time design to address effects of salinity change on populations. Additionally, genetic diversity, a prerequisite for adaptive responses, is reduced in Baltic compared to Atlantic populations. On the one hand, adaptive transgenerational plasticity (TGP) might buffer the effects of environmental change, which may be of particular importance under reduced genetic variation. On the other hand, physiological trade-offs due to environmental stress may hamper parental provisioning to offspring thereby intensifying the impact of climate change across generations (nonadaptive TGP). Here, we studied both hypothesis of adaptive and nonadaptive TGP in the three-spined stickleback (Gasterosteus aculeatus) fish model along the strong salinity gradient of the Baltic Sea in a space-for-time experiment. Each population tolerated desalination well, which was not altered by parental exposure to low salinity. Despite a common marine ancestor, populations locally adapted to low salinity lost their ability to cope with fully marine conditions, resulting in lower survival and reduced relative fitness. Negative transgenerational effects were evident in early life stages, but disappeared after selection via mortality occurred during the first 12-30 days posthatch. Modeling various strengths of selection, we showed that nonadaptive transgenerational plasticity accelerated evolution by increasing directional selection within the offspring generation. Qualitatively, when genetic diversity is large, we predict that such effects will facilitate rapid adaptation and population persistence, while below a certain threshold populations suffer a higher risk of local extinction. Overall, our results suggest that transgenerational plasticity and selection are not independent of each other and thereby highlight a current gap in TGP studies.

Effects of parental acclimation and energy limitation in response to high CO2 exposure in Atlantic cod

Ocean acidification (OA), the dissolution of excess anthropogenic carbon dioxide in ocean waters, is a potential stressor to many marine fish species. Whether species have the potential to acclimate and adapt to changes in the seawater carbonate chemistry is still largely unanswered. Simulation experiments across several generations are challenging for large commercially exploited species because of their long generation times. For Atlantic cod (Gadus morhua), we present first data on the effects of parental acclimation to elevated aquatic CO₂ on larval survival, a fundamental parameter determining population recruitment. The parental generation in this study was exposed to either ambient or elevated aquatic CO₂ levels simulating end-of-century OA levels (~1100 µatm CO₂) for six weeks prior to spawning. Upon fully reciprocal exposure of the F1 generation, we quantified larval survival, combined with two larval feeding regimes in order to investigate the potential effect of energy limitation. We found a significant reduction in larval survival at elevated CO₂ that was partly compensated by parental acclimation to the same CO₂ exposure. Such compensation was only observed in the treatment with high food availability. This complex 3-way interaction indicates that surplus metabolic resources need to be available to allow a transgenerational alleviation response to ocean acidification.

The Baltic Sea as a time machine for the future coastal ocean

Coastal global oceans are expected to undergo drastic changes driven by climate change and increasing anthropogenic pressures in coming decades. Predicting specific future conditions and assessing the best management strategies to maintain ecosystem integrity and sustainable resource use are difficult, because of multiple interacting pressures, uncertain projections, and a lack of test cases for management. We argue that the Baltic Sea can serve as a time machine to study consequences and mitigation of future coastal perturbations, due to its unique combination of an early history of multistressor disturbance and ecosystem deterioration and early implementation of cross-border environmental management to address these problems. The Baltic Sea also stands out in providing a strong scientific foundation and accessibility to long-term data series that provide a unique opportunity to assess the efficacy of management actions to address the breakdown of ecosystem functions. Trend reversals such as the return of top predators, recovering fish stocks, and reduced input of nutrient and harmful substances could be achieved only by implementing an international, cooperative governance structure transcending its complex multistate policy setting, with integrated management of watershed and sea. The Baltic Sea also demonstrates how rapidly progressing global pressures, particularly warming of Baltic waters and the surrounding catchment area, can offset the efficacy of current management approaches. This situation calls for management that is (i) conservative to provide a buffer against regionally unmanageable global perturbations, (ii) adaptive to react to new management challenges, and, ultimately, (iii) multisectorial and integrative to address conflicts associated with economic trade-offs.

Carrying Capacity and Colonization Dynamics of Curvibacter in the Hydra Host Habitat

Most eukaryotic species are colonized by a microbial community – the microbiota – that is acquired during early life stages and is critical to host development and health. Much research has focused on the microbiota biodiversity during the host life, however, empirical data on the basic ecological principles that govern microbiota assembly is lacking. Here we quantify the contribution of colonizer order, arrival time and colonization history to microbiota assembly on a host. We established the freshwater polyp Hydra vulgaris and its dominant colonizer Curvibacter as a model system that enables the visualization and quantification of colonizer population size at the single cell resolution, in vivo, in real time. We estimate the carrying capacity of a single Hydra polyp as 2 × 10⁵Curvibacter cells, which is robust among individuals and time. Colonization experiments reveal a clear priority effect of first colonizers that depends on arrival time and colonization history. First arriving colonizers achieve a numerical advantage over secondary colonizers within a short time lag of 24 h. Furthermore, colonizers primed for the Hydra habitat achieve a numerical advantage in the absence of a time lag. These results follow the theoretical expectations for any bacterial habitat with a finite carrying capacity. Thus, Hydra colonization and succession processes are largely determined by the habitat occupancy over time and Curvibacter colonization history. Our experiments provide empirical data on the basic steps of host-associated microbiota establishment – the colonization stage. The presented approach supplies a framework for studying habitat characteristics and colonization dynamics within the host–microbe setting.

Formation and mosaicity of coccolith segment calcite of the marine algae Emiliania huxleyi

Coccolithophores belong to the most abundant calcium carbonate mineralizing organisms. Coccolithophore biomineralization is a complex and highly regulated process, resulting in a product that strongly differs in its intricate morphology from the abiogenically produced mineral equivalent. Moreover, unlike extracellularly formed biological carbonate hard tissues, coccolith calcite is neither a hybrid composite, nor is it distinguished by a hierarchical microstructure. This is remarkable as the key to optimizing crystalline biomaterials for mechanical strength and toughness lies in the composite nature of the biological hard tissue and the utilization of specific microstructures. To obtain insight into the pathway of biomineralization of Emiliania huxleyi coccoliths, we examine intracrystalline nanostructural features of the coccolith calcite in combination with cell ultrastructural observations related to the formation of the calcite in the coccolith vesicle within the cell. With TEM diffraction and annular dark-field imaging, we prove the presence of planar imperfections in the calcite crystals such as planar mosaic block boundaries. As only minor misorientations occur, we attribute them to dislocation networks creating small-angle boundaries. Intracrystalline occluded biopolymers are not observed. Hence, in E. huxleyi calcite mosaicity is not caused by occluded biopolymers, as it is the case in extracellularly formed hard tissues of marine invertebrates, but by planar defects and dislocations which are typical for crystals formed by classical ion-by-ion growth mechanisms. Using cryo-preparation techniques for SEM and TEM, we found that the membrane of the coccolith vesicle and the outer membrane of the nuclear envelope are in tight proximity, with a well-controlled constant gap of ~4 nm between them. We describe this conspicuous connection as a not yet described interorganelle junction, the "nuclear envelope junction". The narrow gap of this junction likely facilitates transport of Ca²⁺ ions from the nuclear envelope to the coccolith vesicle. On the basis of our observations, we propose that formation of the coccolith utilizes the nuclear envelope-endoplasmic reticulum Ca²⁺ -store of the cell for the transport of Ca²⁺ ions from the external medium to the coccolith vesicle and that E. huxleyi calcite forms by ion-by-ion growth rather than by a nanoparticle accretion mechanism.

Inter- and intraspecific phenotypic plasticity of three phytoplankton species in response to ocean acidification

Phenotypic plasticity describes the phenotypic adjustment of the same genotype to different environmental conditions and is best described by a reaction norm. We focus on the effect of ocean acidification on inter- and intraspecific reaction norms of three globally important phytoplankton species (Emiliania huxleyi, Gephyrocapsa oceanica and Chaetoceros affinis). Despite significant differences in growth rates between the species, they all showed a high potential for phenotypic buffering (similar growth rates between ambient and high CO₂ conditions). Only three coccolithophore genotypes showed a reduced growth in high CO₂. Diverging responses to high CO₂ of single coccolithophore genotypes compared with the respective mean species responses, however, raise the question of whether an extrapolation to the population level is possible from single-genotype experiments. We therefore compared the mean response of all tested genotypes with a total species response comprising the same genotypes, which was not significantly different in the coccolithophores. Assessing species reaction norms to different environmental conditions on short time scale in a genotype-mix could thus reduce sampling effort while increasing predictive power.

Microbial contributions to the persistence of coral reefs

On contemplating the adaptive capacity of reef organisms to a rapidly changing environment, the microbiome offers significant and greatly unrecognised potential. Microbial symbionts contribute to the physiology, development, immunity and behaviour of their hosts, and can respond very rapidly to changing environmental conditions, providing a powerful mechanism for acclimatisation and also possibly rapid evolution of coral reef holobionts. Environmentally acquired fluctuations in the microbiome can have significant functional consequences for the holobiont phenotype upon which selection can act. Environmentally induced changes in microbial abundance may be analogous to host gene duplication, symbiont switching / shuffling as a result of environmental change can either remove or introduce raw genetic material into the holobiont; and horizontal gene transfer can facilitate rapid evolution within microbial strains. Vertical transmission of symbionts is a key feature of many reef holobionts and this would enable environmentally acquired microbial traits to be faithfully passed to future generations, ultimately facilitating microbiome-mediated transgenerational acclimatisation (MMTA) and potentially even adaptation of reef species in a rapidly changing climate. In this commentary, we highlight the capacity and mechanisms for MMTA in reef species, propose a modified Price equation as a framework for assessing MMTA and recommend future areas of research to better understand how microorganisms contribute to the transgenerational acclimatisation of reef organisms, which is essential if we are to reliably predict the consequences of global change for reef ecosystems.

Widespread introgression in deep-sea hydrothermal vent mussels

Background

The analysis of hybrid zones is crucial for gaining a mechanistic understanding of the process of speciation and the maintenance of species boundaries. Hybrid zones have been studied intensively in terrestrial and shallow-water ecosystems, but very little is known about their occurrence in deep-sea environments. Here we used diagnostic, single nucleotide polymorphisms in combination with one mitochondrial gene to re-examine prior hypotheses about a contact zone involving deep-sea hydrothermal vent mussels, Bathymodiolus azoricus and B. puteoserpentis, living along the Mid-Atlantic Ridge.

Results

Admixture was found to be asymmetric with respect to the parental species, while introgression was more widespread geographically than previously recognized. Admixed individuals with a majority of alleles from one of the parental species were most frequent in habitats corresponding to that species. Mussels found at a geographically intermediate vent field constituted a genetically mixed population that showed no evidence for hybrid incompatibilities, a finding that does not support a previously inferred tension zone model.

Conclusions

Our analyses indicate that B. azoricus and B. puteoserpentis hybridize introgressively across a large geographic area without evidence for general hybrid incompatibilities. While these findings shed new light onto the genetic structure of this hybrid zone, many aspects about its nature still remain obscure. Our study sets a baseline for further research that should primarily focus on the acquisition of additional mussel samples and environmental data, a detailed exploration of vent areas and hidden populations as well as genomic analyses in both mussel hosts and their bacterial symbionts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0862-2) contains supplementary material, which is available to authorized users.

Ocean Acidification Effects on Atlantic Cod Larval Survival and Recruitment to the Fished Population

How fisheries will be impacted by climate change is far from understood. While some fish populations may be able to escape global warming via range shifts, they cannot escape ocean acidification (OA), an inevitable consequence of the dissolution of anthropogenic carbon dioxide (CO₂) emissions in marine waters. How ocean acidification affects population dynamics of commercially important fish species is critical for adapting management practices of exploited fish populations. Ocean acidification has been shown to impair fish larvae’s sensory abilities, affect the morphology of otoliths, cause tissue damage and cause behavioural changes. Here, we obtain first experimental mortality estimates for Atlantic cod larvae under OA and incorporate these effects into recruitment models. End-of-century levels of ocean acidification (~1100 μatm according to the IPCC RCP 8.5) resulted in a doubling of daily mortality rates compared to present-day CO₂ concentrations during the first 25 days post hatching (dph), a critical phase for population recruitment. These results were consistent under different feeding regimes, stocking densities and in two cod populations (Western Baltic and Barents Sea stock). When mortality data were included into Ricker-type stock-recruitment models, recruitment was reduced to an average of 8 and 24% of current recruitment for the two populations, respectively. Our results highlight the importance of including vulnerable early life stages when addressing effects of climate change on fish stocks.

Biophysical and Population Genetic Models Predict the Presence of "Phantom" Stepping Stones Connecting Mid-Atlantic Ridge Vent Ecosystems

Deep-sea hydrothermal vents are patchily distributed ecosystems inhabited by specialized animal populations that are textbook meta-populations. Many vent-associated species have free-swimming, dispersive larvae that can establish connections between remote populations. However, connectivity patterns among hydrothermal vents are still poorly understood because the deep sea is undersampled, the molecular tools used to date are of limited resolution, and larval dispersal is difficult to measure directly. A better knowledge of connectivity is urgently needed to develop sound environmental management plans for deep-sea mining. Here, we investigated larval dispersal and contemporary connectivity of ecologically important vent mussels (Bathymodiolus spp.) from the Mid-Atlantic Ridge by using high-resolution ocean modeling and population genetic methods. Even when assuming a long pelagic larval duration, our physical model of larval drift suggested that arrival at localities more than 150 km from the source site is unlikely and that dispersal between populations requires intermediate habitats ("phantom" stepping stones). Dispersal patterns showed strong spatiotemporal variability, making predictions of population connectivity challenging. The assumption that mussel populations are only connected via additional stepping stones was supported by contemporary migration rates based on neutral genetic markers. Analyses of population structure confirmed the presence of two southern and two hybridizing northern mussel lineages that exhibited a substantial, though incomplete, genetic differentiation. Our study provides insights into how vent animals can disperse between widely separated vent habitats and shows that recolonization of perturbed vent sites will be subject to chance events, unless connectivity is explicitly considered in the selection of conservation areas.

Long-term dynamics of adaptive evolution in a globally important phytoplankton species to ocean acidification

Marine phytoplankton may adapt to ocean change, such as acidification or warming, because of their large population sizes and short generation times. Long-term adaptation to novel environments is a dynamic process, and phenotypic change can take place thousands of generations after exposure to novel conditions. We conducted a long-term evolution experiment (4 years = 2100 generations), starting with a single clone of the abundant and widespread coccolithophore Emiliania huxleyi exposed to three different CO2 levels simulating ocean acidification (OA). Growth rates as a proxy for Darwinian fitness increased only moderately under both levels of OA [+3.4% and +4.8%, respectively, at 1100 and 2200 μatm partial pressure of CO2 (Pco2)] relative to control treatments (ambient CO2, 400 μatm). Long-term adaptation to OA was complex, and initial phenotypic responses of ecologically important traits were later reverted. The biogeochemically important trait of calcification, in particular, that had initially been restored within the first year of evolution was later reduced to levels lower than the performance of nonadapted populations under OA. Calcification was not constitutively lost but returned to control treatment levels when high CO2-adapted isolates were transferred back to present-day control CO2 conditions. Selection under elevated CO2 exacerbated a general decrease of cell sizes under long-term laboratory evolution. Our results show that phytoplankton may evolve complex phenotypic plasticity that can affect biogeochemically important traits, such as calcification. Adaptive evolution may play out over longer time scales (>1 year) in an unforeseen way under future ocean conditions that cannot be predicted from initial adaptation responses.

Swift thermal reaction norm evolution in a key marine phytoplankton species

Temperature has a profound effect on the species composition and physiology of marine phytoplankton, a polyphyletic group of microbes responsible for half of global primary production. Here, we ask whether and how thermal reaction norms in a key calcifying species, the coccolithophore Emiliania huxleyi, change as a result of 2.5 years of experimental evolution to a temperature ≈2°C below its upper thermal limit. Replicate experimental populations derived from a single genotype isolated from Norwegian coastal waters were grown at two temperatures for 2.5 years before assessing thermal responses at 6 temperatures ranging from 15 to 26°C, with pCO ₂ (400/1100/2200 μatm) as a fully factorial additional factor. The two selection temperatures (15°/26.3°C) led to a marked divergence of thermal reaction norms. Optimal growth temperatures were 0.7°C higher in experimental populations selected at 26.3°C than those selected at 15.0°C. An additional negative effect of high pCO ₂ on maximal growth rate (8% decrease relative to lowest level) was observed. Finally, the maximum persistence temperature (T_max) differed by 1-3°C between experimental treatments, as a result of an interaction between pCO ₂ and the temperature selection. Taken together, we demonstrate that several attributes of thermal reaction norms in phytoplankton may change faster than the predicted progression of ocean warming.

Transcriptome profiling of immune tissues reveals habitat-specific gene expression between lake and river sticklebacks

The observation of habitat-specific phenotypes suggests the action of natural selection. The three-spined stickleback (Gasterosteus aculeatus) has repeatedly colonized and adapted to diverse freshwater habitats across the northern hemisphere since the last glaciation, while giving rise to recurring phenotypes associated with specific habitats. Parapatric lake and river populations of sticklebacks harbour distinct parasite communities, a factor proposed to contribute to adaptive differentiation between these ecotypes. However, little is known about the transcriptional response to the distinct parasite pressure of those fish in a natural setting. Here, we sampled wild-caught sticklebacks across four geographical locations from lake and river habitats differing in their parasite load. We compared gene expression profiles between lake and river populations using 77 whole-transcriptome libraries from two immune-relevant tissues, the head kidney and the spleen. Differential expression analyses revealed 139 genes with habitat-specific expression patterns across the sampled population pairs. Among the 139 differentially expressed genes, eight are annotated with an immune function and 42 have been identified as differentially expressed in previous experimental studies in which fish have been immune challenged. Together, these findings reinforce the hypothesis that parasites contribute to adaptation of sticklebacks in lake and river habitats.

Dietary β-glucan (MacroGard®) enhances survival of first feeding turbot (Scophthalmus maximus) larvae by altering immunity, metabolism and microbiota

Reflecting the natural biology of mass spawning fish aquaculture production of fish larvae is often hampered by high and unpredictable mortality rates. The present study aimed to enhance larval performance and immunity via the oral administration of an immunomodulator, β-glucan (MacroGard(®)) in turbot (Scophthalmus maximus). Rotifers (Brachionus plicatilis) were incubated with or without yeast β-1,3/1,6-glucan in form of MacroGard(®) at a concentration of 0.5 g/L. Rotifers were fed to first feeding turbot larvae once a day. From day 13 dph onwards all tanks were additionally fed untreated Artemia sp. nauplii (1 nauplius ml/L). Daily mortality was monitored and larvae were sampled at 11 and 24 dph for expression of 30 genes, microbiota analysis, trypsin activity and size measurements. Along with the feeding of β-glucan daily mortality was significantly reduced by ca. 15% and an alteration of the larval microbiota was observed. At 11 dph gene expression of trypsin and chymotrypsin was elevated in the MacroGard(®) fed fish, which resulted in heightened tryptic enzyme activity. No effect on genes encoding antioxidative proteins was observed, whilst the immune response was clearly modulated by β-glucan. At 11 dph complement component c3 was elevated whilst cytokines, antimicrobial peptides, toll like receptor 3 and heat shock protein 70 were not affected. At the later time point (24 dph) an anti-inflammatory effect in form of a down-regulation of hsp 70, tnf-α and il-1β was observed. We conclude that the administration of MacroGard(®) induced an immunomodulatory response and could be used as an effective measure to increase survival in rearing of turbot.

It is the economy, stupid! Projecting the fate of fish populations using ecological-economic modeling

Four marine fish species are among the most important on the world market: cod, salmon, tuna, and sea bass. While the supply of North American and European markets for two of these species - Atlantic salmon and European sea bass - mainly comes from fish farming, Atlantic cod and tunas are mainly caught from wild stocks. We address the question what will be the status of these wild stocks in the midterm future, in the year 2048, to be specific. Whereas the effects of climate change and ecological driving forces on fish stocks have already gained much attention, our prime interest is in studying the effects of changing economic drivers, as well as the impact of variable management effectiveness. Using a process-based ecological-economic multispecies optimization model, we assess the future stock status under different scenarios of change. We simulate (i) technological progress in fishing, (ii) increasing demand for fish, and (iii) increasing supply of farmed fish, as well as the interplay of these driving forces under different scenarios of (limited) fishery management effectiveness. We find that economic change has a substantial effect on fish populations. Increasing aquaculture production can dampen the fishing pressure on wild stocks, but this effect is likely to be overwhelmed by increasing demand and technological progress, both increasing fishing pressure. The only solution to avoid collapse of the majority of stocks is institutional change to improve management effectiveness significantly above the current state. We conclude that full recognition of economic drivers of change will be needed to successfully develop an integrated ecosystem management and to sustain the wild fish stocks until 2048 and beyond.

Immunity comes first: the effect of parasite genotypes on adaptive immunity and immunization in three-spined sticklebacks

Adaptive immunity in vertebrates can confer increased resistance against invading pathogens upon re-infection. But how specific parasite genotypes affect the temporal transition from innate to adaptive immunity under continual exposure to parasites is poorly understood. Here, we investigated the effects of homologous and heterologous exposures of genetically distinct parasite lineages of the eye fluke Diplostomum pseudospathaceum on gene expression patterns of adaptive immunity in sticklebacks (Gasterosteus aculeatus). Observable differences in gene expression were largely attributable to final exposures while there was no transcription pattern characteristic for a general response to repeated infections with D. pseudospathaceum. None of the final exposure treatments was able to erase the distinct expression patterns resulting from a heterologous pre-exposed fish. Interestingly, heterologous final exposures showed similarities between different treatment groups subjected to homologous pre-exposure. The observed pattern was supported by parasite infection rates and suggests that host immunization was optimized towards an adaptive immune response that favored effectiveness against parasite diversity over specificity.

Experimental assessment of critical anthropogenic sediment burial in eelgrass Zostera marina

Seagrass meadows, one of the world's most important and productive coastal habitats, are threatened by a range of anthropogenic actions. Burial of seagrass plants due to coastal activities is one important anthropogenic pressure leading to the decline of local populations. In our study, we assessed the response of eelgrass Zostera marina to sediment burial from physiological, morphological, and population parameters. In a full factorial field experiment, burial level (5-20cm) and burial duration (4-16weeks) were manipulated. Negative effects were visible even at the lowest burial level (5cm) and shortest duration (4weeks), with increasing effects over time and burial level. Buried seagrasses showed higher shoot mortality, delayed growth and flowering and lower carbohydrate storage. The observed effects will likely have an impact on next year's survival of buried plants. Our results have implications for the management of this important coastal plant.

Genomics of divergence along a continuum of parapatric population differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.

A variety of evolutionary forces influence the genomic landscape of divergence during ecological speciation. Here we characterize the evolution of genomic divergence patterns based on 60 fully sequenced three-spined stickleback genomes, contrasting lake and river populations that differ in parasite abundance. Our comparison of the size and abundance of divergent regions in the genomes across a continuum of population differentiation suggests that selection and the hitchhiking effect on neutral sites mainly contributes to the observed heterogeneous patterns of genomic divergence. Additional divergent regions of the genome can be explained by a local reduction of gene flow. Our description of genomic divergence patterns across a continuum of population differentiation combined with an analysis of molecular signatures of evolution highlights how adaptation shapes the differentiation of sticklebacks in freshwater habitats.

Gene expression changes in the coccolithophore Emiliania huxleyi after 500 generations of selection to ocean acidification

Coccolithophores are unicellular marine algae that produce biogenic calcite scales and substantially contribute to marine primary production and carbon export to the deep ocean. Ongoing ocean acidification particularly impairs calcifying organisms, mostly resulting in decreased growth and calcification. Recent studies revealed that the immediate physiological response in the coccolithophore Emiliania huxleyi to ocean acidification may be partially compensated by evolutionary adaptation, yet the underlying molecular mechanisms are currently unknown. Here, we report on the expression levels of 10 candidate genes putatively relevant to pH regulation, carbon transport, calcification and photosynthesis in E. huxleyi populations short-term exposed to ocean acidification conditions after acclimation (physiological response) and after 500 generations of high CO2 adaptation (adaptive response). The physiological response revealed downregulation of candidate genes, well reflecting the concomitant decrease of growth and calcification. In the adaptive response, putative pH regulation and carbon transport genes were up-regulated, matching partial restoration of growth and calcification in high CO2-adapted populations. Adaptation to ocean acidification in E. huxleyi likely involved improved cellular pH regulation, presumably indirectly affecting calcification. Adaptive evolution may thus have the potential to partially restore cellular pH regulatory capacity and thereby mitigate adverse effects of ocean acidification.

Specific gene expression responses to parasite genotypes reveal redundancy of innate immunity in vertebrates

Vertebrate innate immunity is the first line of defense against an invading pathogen and has long been assumed to be largely unspecific with respect to parasite/pathogen species. However, recent phenotypic evidence suggests that immunogenetic variation, i.e. allelic variability in genes associated with the immune system, results in host-parasite genotype-by-genotype interactions and thus specific innate immune responses. Immunogenetic variation is common in all vertebrate taxa and this reflects an effective immunological function in complex environments. However, the underlying variability in host gene expression patterns as response of innate immunity to within-species genetic diversity of macroparasites in vertebrates is unknown. We hypothesized that intra-specific variation among parasite genotypes must be reflected in host gene expression patterns. Here we used high-throughput RNA-sequencing to examine the effect of parasite genotypes on gene expression patterns of a vertebrate host, the three-spined stickleback (Gasterosteus aculeatus). By infecting naïve fish with distinct trematode genotypes of the species Diplostomum pseudospathaceum we show that gene activity of innate immunity in three-spined sticklebacks depended on the identity of an infecting macroparasite genotype. In addition to a suite of genes indicative for a general response against the trematode we also find parasite-strain specific gene expression, in particular in the complement system genes, despite similar infection rates of single clone treatments. The observed discrepancy between infection rates and gene expression indicates the presence of alternative pathways which execute similar functions. This suggests that the innate immune system can induce redundant responses specific to parasite genotypes.

Current European Labyrinthula zosterae are not virulent and modulate seagrass (Zostera marina) defense gene expression

Pro- and eukaryotic microbes associated with multi-cellular organisms are receiving increasing attention as a driving factor in ecosystems. Endophytes in plants can change host performance by altering nutrient uptake, secondary metabolite production or defense mechanisms. Recent studies detected widespread prevalence of Labyrinthula zosterae in European Zostera marina meadows, a protist that allegedly caused a massive amphi-Atlantic seagrass die-off event in the 1930's, while showing only limited virulence today. As a limiting factor for pathogenicity, we investigated genotype × genotype interactions of host and pathogen from different regions (10-100 km-scale) through reciprocal infection. Although the endophyte rapidly infected Z. marina, we found little evidence that Z. marina was negatively impacted by L. zosterae. Instead Z. marina showed enhanced leaf growth and kept endophyte abundance low. Moreover, we found almost no interaction of protist × eelgrass-origin on different parameters of L. zosterae virulence/Z. marina performance, and also no increase in mortality after experimental infection. In a target gene approach, we identified a significant down-regulation in the expression of 6/11 genes from the defense cascade of Z. marina after real-time quantitative PCR, revealing strong immune modulation of the host's defense by a potential parasite for the first time in a marine plant. Nevertheless, one gene involved in phenol synthesis was strongly up-regulated, indicating that Z. marina plants were probably able to control the level of infection. There was no change in expression in a general stress indicator gene (HSP70). Mean L. zosterae abundances decreased below 10% after 16 days of experimental runtime. We conclude that under non-stress conditions L. zosterae infection in the study region is not associated with substantial virulence.