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

The metabolic network response and tolerance mechanism of Thalassia hemprichii under high sulfide based on widely targeted metabolome and transcriptome

Costal eutrophication leads to increased sulfide levels in sediments, which has been identified as a major cause of the global decline in seagrass beds. The seagrass Thalassia hemprichii, a dominant tropical species in the Indo-Pacific, is facing a potential threat from sulfide, which can be easily reduced from sulfate in porewater under the influence of global climate change and eutrophication. However, its metabolic response and tolerance mechanisms to high sulfide remain unclear. Thus, the current study investigated the physiological responses and programmed metabolic networks of T. hemprichii through a three-week mesocosm experiment, integrating physiology, stable isotope, widely targeted metabolomics, transcriptomics, and microbial diversity assessments. High sulfide reduced the sediment microbial diversity, while increased sediment sulfate reduced bacterial abundance and δ34S. The exposure to sulfide enhanced root δ34S while decreased leaf δ34S in T. hemprichii. High sulfide was shown to inhibit photosynthesis via damaging PSII, which further reduced ATP production. In response, abundant up-regulated differentially expressed genes in energy metabolism, especially in oxidative phosphorylation, were activated to compensate high energy requirement. High sulfide also promoted autophagy by overexpressing the genes related to phagocytosis and phagolysosome. Meanwhile, metabolomic profiling revealed that the contents of many primary metabolites, such as carbohydrates and amino acids, were reduced in both leaves and roots, likely to provide more energy and synthesize stress-responsive secondary metabolites. Genes related to nitrate reduction and transportation were up-regulated to promote N uptake for sulfide detoxification. High sulfide levels specifically enhanced thiamine in roots, while increased jasmonic acid and flavonoid levels in leaves. The distinct differences in metabolism between roots and leaves might be related to sulfide levels and the growth-defense trade-off. Collectively, our work highlights the specific mechanisms underlying the response and tolerance of T. hemprichii to high sulfide, providing new insights into seagrass strategies for resisting sulfide.

Decline of seagrass (Posidonia oceanica) production over two decades in the face of warming of the Eastern Mediterranean Sea

The response of Posidonia oceanica meadows to global warming of the Eastern Mediterranean Sea, where the increase in sea surface temperature (SST) is particularly severe, is poorly investigated. Here, we reconstructed the long-term P. oceanica production in 60 meadows along the Greek Seas over two decades (1997-2018), using lepidochronology. We determined the effect of warming on production by reconstructing the annual and maximum (i.e. August) SST, considering the role of other production drivers related to water quality (i.e. Chla, suspended particulate matter, Secchi depth). Grand mean (±SE) production across all sites and the study period was 48 ± 1.1 mg DW per shoot yr-1 . Production over the last two decades followed a trajectory of decrease, which was related to the concurrent increase in annual SST and SSTaug . Annual SST > 20°C and SSTaug  > 26.5°C was related to production decline (GAMM, P < 0.05), while the rest of the tested factors did not help explain the production pattern. Our results indicate a persistent and increasing threat for Eastern Mediterranean meadows, drawing attention to management authorities, highlighting the necessity of reducing local impacts to enhance the resilience of seagrass meadows to global change threats.

Signs of local adaptation by genetic selection and isolation promoted by extreme temperature and salinity in the Mediterranean seagrass Posidonia oceanica

Adaptation to local conditions is known to occur in seagrasses; however, knowledge of the genetic basis underlying this phenomenon remains scarce. Here, we analysed Posidonia oceanica from six sites within and around the Stagnone di Marsala, a semi-enclosed coastal lagoon where salinity and temperature exceed the generally described tolerance thresholds of the species. Sea surface temperatures (SSTs) were measured and plant samples were collected for the assessment of morphology, flowering rate and for screening genome-wide polymorphisms using double digest restriction-site-associated DNA sequencing. Results demonstrated more extreme SSTs and salinity levels inside the lagoon than the outer lagoon regions. Morphological results showed significantly fewer and shorter leaves and reduced rhizome growth of P. oceanica from the inner lagoon and past flowering events were recorded only for a meadow farthest away from the lagoon. Using an array of 51,329 single nucleotide polymorphisms, we revealed a clear genetic structure among the study sites and confirmed the genetic isolation and high clonality of the innermost site. In all, 14 outlier loci were identified and annotated with several proteins including those relate to plant stress response, protein transport and regulators of plant-specific developmental events. Especially, five outlier loci showed maximum allele frequency at the innermost site, likely reflecting adaptation to the extreme temperature and salinity regimes, possibly due to the selection of more resistant genotypes and the progressive restriction of gene flow. Overall, this study helps us to disentangle the genetic basis of seagrass adaptation to local environmental conditions and may support future works on assisted evolution in seagrasses.

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.

Ten years of marine evolutionary biology-Challenges and achievements of a multidisciplinary research initiative

The Centre for Marine Evolutionary Biology (CeMEB) at the University of Gothenburg, Sweden, was established in 2008 through a 10-year research grant of 8.7 m€ to a team of senior researchers. Today, CeMEB members have contributed >500 scientific publications, 30 PhD theses and have organised 75 meetings and courses, including 18 three-day meetings and four conferences. What are the footprints of CeMEB, and how will the centre continue to play a national and international role as an important node of marine evolutionary research? In this perspective article, we first look back over the 10 years of CeMEB activities and briefly survey some of the many achievements of CeMEB. We furthermore compare the initial goals, as formulated in the grant application, with what has been achieved, and discuss challenges and milestones along the way. Finally, we bring forward some general lessons that can be learnt from a research funding of this type, and we also look ahead, discussing how CeMEB's achievements and lessons can be used as a springboard to the future of marine evolutionary biology.

Single-molecule real-time sequencing of the full-length transcriptome of Halophila beccarii

Ecologically, Halophila beccarii Asch. is considered as a colonizing or a pioneer seagrass species and a “tiny but mighty” seagrass species, since it may recover quickly from disturbance generally. The use of transcriptome technology can provide a better understanding of the physiological processes of seagrasses. To date, little is known about the genome and transcriptome information of H. beccarii. In this study, we used single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterize the transcriptome structure. A total of 11,773 of the 15,348 transcripts were successfully annotated in seven databases. In addition, 1573 long non-coding RNAs, 8402 simple sequence repeats and 2567 transcription factors were predicted in all the transcripts. A GO analysis showed that 5843 transcripts were divided into three categories, including biological process (BP), cellular component (CC) and molecular function (MF). In these three categories, metabolic process (1603 transcripts), protein-containing complex (515 transcripts) and binding (3233 transcripts) were the primary terms in BP, CC, and MF, respectively. The major types of transcription factors were involved in MYB-related and NF-YB families. To the best of our knowledge, this is the first report of the transcriptome of H. beccarii using SMRT sequencing technology.

A Pleistocene legacy structures variation in modern seagrass ecosystems

Sustaining biodiversity and ecosystems in the long term depends on their adjustment to a rapidly changing climate. By characterizing the structure of the marine plant eelgrass and associated communities at 50 sites across its broad range, we found that eelgrass growth form and biomass retain a legacy of Pleistocene range shifts and genetic bottlenecks that in turn affect the biomass of algae and invertebrates that fuel coastal food webs. The ecosystem-level effects of this ancient evolutionary legacy are comparable to or stronger than effects of current environmental forcing, suggesting that this economically important ecosystem may be unable to keep pace with rapid global change.

Distribution of Earth’s biomes is structured by the match between climate and plant traits, which in turn shape associated communities and ecosystem processes and services. However, that climate–trait match can be disrupted by historical events, with lasting ecosystem impacts. As Earth’s environment changes faster than at any time in human history, critical questions are whether and how organismal traits and ecosystems can adjust to altered conditions. We quantified the relative importance of current environmental forcing versus evolutionary history in shaping the growth form (stature and biomass) and associated community of eelgrass (Zostera marina), a widespread foundation plant of marine ecosystems along Northern Hemisphere coastlines, which experienced major shifts in distribution and genetic composition during the Pleistocene. We found that eelgrass stature and biomass retain a legacy of the Pleistocene colonization of the Atlantic from the ancestral Pacific range and of more recent within-basin bottlenecks and genetic differentiation. This evolutionary legacy in turn influences the biomass of associated algae and invertebrates that fuel coastal food webs, with effects comparable to or stronger than effects of current environmental forcing. Such historical lags in phenotypic acclimatization may constrain ecosystem adjustments to rapid anthropogenic climate change, thus altering predictions about the future functioning of ecosystems.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

Evaluating Seagrass Meadow Dynamics by Integrating Field-Based and Remote Sensing Techniques

Marine phanerogams are considered biological sentinels or indicators since any modification in seagrass meadow distribution and coverage signals negative changes in the marine environment. In recent decades, seagrass meadows have undergone global losses at accelerating rates, and almost one-third of their coverage has disappeared globally. This study focused on the dynamics of seagrass meadows in the northern Adriatic Sea, which is one of the most anthropogenically affected areas in the Mediterranean Sea. Seagrass distribution data and remote sensing products were utilized to identify the stable and dynamic parts of the seagrass ecosystem. Different seagrass species could not be distinguished with the Sentinel-2 (BOA) satellite image. However, results revealed a generally stable seagrass meadow (283.5 Ha) but, on the other hand, a stochastic behavior in seagrass meadow retraction (90.8 Ha) linked to local environmental processes associated with anthropogenic activities or climate change. If systemized, this proposed approach to monitoring seagrass meadow dynamics could be developed as a spatial decision support system for the entire Mediterranean basin. Such a tool could serve as a key element for decision makers in marine protected areas and would potentially support more effective conservation and management actions in these highly productive and important environments.

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.

Overexpression of seagrass DnaJ gene ZjDjB1 enhances the thermotolerance of transgenic arabidopsis thaliana

Seagrass meadows are one of the most important marine resources that grow along the coast. They provide habitat and a food source for animals. They also protect the coast, fix sediment and purify seawater. In the current period of global climate change, anomalies in coastal water temperatures are increasing. A sudden increase in water temperature owing to a heat wave can have a profound effect on seagrass. Zostera japonica is a type of intertidal seagrasses, which is exposed to the air at low tide. High temperatures in the summer often lead to a decline in seagrass meadows. DnaJ proteins, also known as J proteins, are a family of conserved chaperone proteins. They are designated as J proteins because they contain a highly conserved J domain. They function as chaperones of heat shock proteins in organisms. In this study, the role of DnaJ protein (ZjDjB1) of Z. japonica under heat stress was studied. ZjDjB1 was localized to the cytoplasm and nucleus. The overexpression of ZjDjB1 in Arabidopsis thaliana results in an increase in thermotolerance and a decrease in the accumulation of reactive oxygen species and also a reduction in membrane damage. ZjDjB1 may achieve this goal by maintaining a low activity of proteolytic enzymes.

Gene body DNA methylation in seagrasses: inter- and intraspecific differences and interaction with transcriptome plasticity under heat stress

The role of DNA methylation and its interaction with gene expression and transcriptome plasticity is poorly understood, and current insight comes mainly from studies in very few model plant species. Here, we study gene body DNA methylation (gbM) and gene expression patterns in ecotypes from contrasting thermal environments of two marine plants with contrasting life history strategies in order to explore the potential role epigenetic mechanisms could play in gene plasticity and responsiveness to heat stress. In silico transcriptome analysis of CpGO/E ratios suggested that the bulk of Posidonia oceanica and Cymodocea nodosa genes possess high levels of intragenic methylation. We also observed a correlation between gbM and gene expression flexibility: genes with low DNA methylation tend to show flexible gene expression and plasticity under changing conditions. Furthermore, the empirical determination of global DNA methylation (5-mC) showed patterns of intra and inter-specific divergence that suggests a link between methylation level and the plants' latitude of origin and life history. Although we cannot discern whether gbM regulates gene expression or vice versa, or if other molecular mechanisms play a role in facilitating transcriptome responsiveness, our findings point to the existence of a relationship between gene responsiveness and gbM patterns in marine plants.

Local adaptation contributes to gene expression divergence in maize

Gene expression links genotypes to phenotypes, so identifying genes whose expression is shaped by selection will be important for understanding the traits and processes underlying local adaptation. However, detecting local adaptation for gene expression will require distinguishing between divergence due to selection and divergence due to genetic drift. Here, we adapt a QST−FST framework to detect local adaptation for transcriptome-wide gene expression levels in a population of diverse maize genotypes. We compare the number and types of selected genes across a wide range of maize populations and tissues, as well as selection on cold-response genes, drought-response genes, and coexpression clusters. We identify a number of genes whose expression levels are consistent with local adaptation and show that genes involved in stress response show enrichment for selection. Due to its history of intense selective breeding and domestication, maize evolution has long been of interest to researchers, and our study provides insight into the genes and processes important for in local adaptation of maize.

Intraspecific genetic variation matters when predicting seagrass distribution under climate change

Seagrasses play a vital role in structuring coastal marine ecosystems, but their distributional range and genetic diversity have declined rapidly in recent decades. To improve conservation of seagrass species, it is important to predict how climate change may impact their ranges. Such predictions are typically made with correlative species distribution models (SDMs), which can estimate a species' potential distribution under present and future climatic scenarios given species' presence data and climatic predictor variables. However, these models are typically constructed with species-level data, and thus ignore intraspecific genetic variability, which can give rise to populations with adaptations to heterogeneous climatic conditions. Here, we explore the link between intraspecific adaptation and niche differentiation in Thalassia hemprichii, a seagrass broadly distributed in the tropical Indo-Pacific Ocean and a crucial provider of habitat for numerous marine species. By retrieving and re-analysing microsatellite data from previous studies, we delimited two distinct phylogeographical lineages within the nominal species and found an intermediate level of differentiation in their multidimensional environmental niches, suggesting the possibility for local adaptation. We then compared projections of the species' habitat suitability under climate change scenarios using species-level and lineage-level SDMs. In the Central Tropical Indo-Pacific region, models for both levels predicted considerable range contraction in the future, but the lineage-level models predicted more severe habitat loss. Importantly, the two modelling approaches predicted opposite patterns of habitat change in the Western Tropical Indo-Pacific region. Our results highlight the necessity of conserving distinct populations and genetic pools to avoid regional extinction due to climate change and have important implications for guiding future management of seagrasses.

Seagrasses in an era of ocean warming: a review

Seagrasses are valuable sources of food and habitat for marine life and are one of Earth's most efficient carbon sinks. However, they are facing a global decline due to ocean warming and eutrophication. In the last decade, with the advent of new technology and molecular advances, there has been a dramatic increase in the number of studies focusing on the effects of ocean warming on seagrasses. Here, we provide a comprehensive review of the future of seagrasses in an era of ocean warming. We have gathered information from published studies to identify potential commonalities in the effects of warming and the responses of seagrasses across four distinct levels: molecular, biochemical/physiological, morphological/population, and ecosystem/planetary. To date, we know that although warming strongly affects seagrasses at all four levels, seagrass responses diverge amongst species, populations, and over depths. Furthermore, warming alters seagrass distribution causing massive die-offs in some seagrass populations, whilst also causing tropicalization and migration of temperate species. In this review, we evaluate the combined effects of ocean warming with other environmental stressors and emphasize the need for multiple-stressor studies to provide a deeper understanding of seagrass resilience. We conclude by discussing the most significant knowledge gaps and future directions for seagrass research.

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.

The Seagrass Methylome Is Associated With Variation in Photosynthetic Performance Among Clonal Shoots

Evolutionary theory predicts that clonal organisms are more susceptible to extinction than sexually reproducing organisms, due to low genetic variation and slow rates of evolution. In agreement, conservation management considers genetic variation as the ultimate measure of a population's ability to survive over time. However, clonal plants are among the oldest living organisms on our planet. Here, we test the hypothesis that clonal seagrass meadows display epigenetic variation that complements genetic variation as a source of phenotypic variation. In a clonal meadow of the seagrass Zostera marina, we characterized DNA methylation among 42 shoots. We also sequenced the whole genome of 10 shoots to correlate methylation patterns with photosynthetic performance under exposure to and recovery from 27°C, while controlling for somatic mutations. Here, we show for the first time that clonal seagrass shoots display DNA methylation variation that is independent from underlying genetic variation, and associated with variation in photosynthetic performance under experimental conditions. It remains unknown to what degree this association could be influenced by epigenetic responses to transplantation-related stress, given that the methylomes showed a strong shift under acclimation to laboratory conditions. The lack of untreated control samples in the heat stress experiment did not allow us to distinguish methylome shifts induced by acclimation from such induced by heat stress. Notwithstanding, the co-variation in DNA methylation and photosynthetic performance may be linked via gene expression because methylation patterns varied in functionally relevant genes involved in photosynthesis, and in the repair and prevention of heat-induced protein damage. While genotypic diversity has been shown to enhance stress resilience in seagrass meadows, we suggest that epigenetic variation plays a similar role in meadows dominated by a single genotype. Consequently, conservation management of clonal plants should consider epigenetic variation as indicator of resilience and stability.

Response of foundation macrophytes to near-natural simulated marine heatwaves

Marine heatwaves have been observed worldwide and are expected to increase in both frequency and intensity due to climate change. Such events may cause ecosystem reconfigurations arising from species range contraction or redistribution, with ecological, economic and social implications. Macrophytes such as the brown seaweed Fucus vesiculosus and the seagrass Zostera marina are foundation species in many coastal ecosystems of the temperate northern hemisphere. Hence, their response to extreme events can potentially determine the fate of associated ecosystems. Macrophyte functioning is intimately linked to the maintenance of photosynthesis, growth and reproduction, and resistance against pathogens, epibionts and grazers. We investigated morphological, physiological, pathological and chemical defence responses of western Baltic Sea F. vesiculosus and Z. marina populations to simulated near-natural marine heatwaves. Along with (a) the control, which constituted no heatwave but natural stochastic temperature variability (0HW), two treatments were applied: (b) two late-spring heatwaves (June, July) followed by a summer heatwave (August; 3HW) and (c) a summer heatwave only (1HW). The 3HW treatment was applied to test whether preconditioning events can modulate the potential sensitivity to the summer heatwave. Despite the variety of responses measured in both species, only Z. marina growth was impaired by the accumulative heat stress imposed by the 3HW treatment. Photosynthetic rate, however, remained high after the last heatwave indicating potential for recovery. Only epibacterial abundance was significantly affected in F. vesiculosus. Hence both macrophytes, and in particular F. vesiculosus, seem to be fairly tolerant to short-term marine heatwaves at least at the intensities applied in this experiment (up to 5°C above mean temperature over a period of 9 days). This may partly be due to the fact that F. vesiculosus grows in a highly variable environment, and may have a high phenotypic plasticity.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

Surface chemical defence of the eelgrass Zostera marina against microbial foulers

Plants rely on both mechanical and chemical defence mechanisms to protect their surfaces against microorganisms. The recently completed genome of the eelgrass Zostera marina, a marine angiosperm with fundamental importance for coastal ecosystems, showed that its re-adaptation from land to the sea has led to the loss of essential genes (for chemical communication and defence) and structural features (stomata and thick cuticle) that are typical of terrestrial plants. This study was designed to understand the molecular nature of surface protection and fouling-control strategy of eelgrass against marine epiphytic yeasts. Different surface extraction methods and comparative metabolomics by tandem mass spectrometry (LC-MS/MS) were used for targeted and untargeted identification of the metabolite profiles of the leaf surface and the whole tissue extracts. Desorption electrospray ionization-imaging mass spectrometry (DESI-IMS) coupled with traditional bioassays revealed, for the first time, the unique spatial distribution of the eelgrass surface-associated phenolics and fatty acids, as well as their differential bioactivity against the growth and settlement of epiphytic yeasts. This study provides insights into the complex chemical defence system of the eelgrass leaf surface. It suggests that surface-associated metabolites modulate biotic interactions and provide chemical defence and structural protection to eelgrass in its marine environment.

Rapid adaptation to cold in the invasive cane toad Rhinella marina

Understanding rapid adaptation to novel environments is essential as we face increasing climatic change. Invasive species are an ideal system for studying adaptation as they are typically introduced to novel environments where they must adapt if they are to persist. We used the invasive cane toad, Rhinella marina, to investigate the contribution of plasticity and evolution to rapid adaptation in a novel environment. Rhinella marina is a neotropical toad that has invaded areas with climates outside of its native environmental niche. The goal of this research was to understand how cane toads persist in northern Florida, the coldest region of their combined natural and invasive range, and originally thought to be beyond their thermal breadth. We measured Critical thermal minima in cane toads from the original, warm introduction location (Miami), and their northern range edge (Tampa) to determine whether northern toads were more cold-tolerant, and to examine the contribution of adaptive plasticity and evolution to any changes in tolerance. Our results show that following acclimation to cold temperatures, southern toads are less tolerant of cold than northern toads. This persistent population difference implies selection for cold-tolerance in northern populations. Differences in individual responses indicate that plasticity is also involved in this response. Our findings have implications for conservation because predatory cane toad invasions threaten local faunas, especially native amphibians. Characterizing specific adaptive mechanisms that allow R. marina to expand its range will identify evolutionary processes that shape a highly successful invasive species.

Differentiation in fitness-related traits in response to elevated temperatures between leading and trailing edge populations of marine macrophytes

The nature of species distribution boundaries is a key subject in ecology and evolution. Edge populations are potentially more exposed to climate-related environmental pressures. Despite research efforts, little is known about variability in fitness-related traits in leading (i.e., colder, high latitude) versus trailing (i.e., warmer, low latitude) edge populations. We tested whether the resilience, i.e. the resistance and recovery, of key traits differs between a distributional cold (Greenland) and warm (Portugal) range edge population of two foundation marine macrophytes, the intertidal macroalga Fucus vesiculosus and the subtidal seagrass Zostera marina. The resistance and recovery of edge populations to elevated seawater temperatures was compared under common experimental conditions using photosynthetic efficiency and expression of heat shock proteins (HSP). Cold and warm edge populations differed in their response, but this was species specific. The warm edge population of F. vesiculosus showed higher thermal resistance and recovery whereas the cold leading edge was less tolerant. The opposite was observed in Z. marina, with reduced recovery at the warm edge, while the cold edge was not markedly affected by warming. Our results confirm that differentiation of thermal stress responses can occur between leading and trailing edges, but such responses depend on local population traits and are thus not predictable just based on thermal pressures.

Adaptive responses along a depth and a latitudinal gradient in the endemic seagrass Posidonia oceanica

Seagrass meadows provide important ecosystem services and are critical for the survival of the associated invertebrate community. However, they are threatened worldwide by human-driven environmental change. Understanding the seagrasses' potential for adaptation is critical to assess not only their ability to persist under future global change scenarios, but also to assess the persistence of the associated communities. Here we screened a wild population of Posidonia oceanica, an endemic long-lived seagrass in the Mediterranean Sea, for genes that may be target of environmental selection, using an outlier and a genome-wide transcriptome analysis. We identified loci where polymorphism or differential expression was associated with either a latitudinal or a bathymetric gradient, as well as with both gradients in an effort to identify loci associated with temperature and light. We found the candidate genes underlying growth and immunity to be divergent between populations adapted to different latitudes and/or depths, providing evidence for local adaptation. Furthermore, we found evidence of reduced gene flow among populations including adjacent populations. Reduced gene flow, combined with low sexual recombination, small effective population size, and long generation time of P. oceanica raises concerns for the long-term persistence of this species, especially in the face of rapid environmental change driven by human activities.

Unique maternal and environmental effects on the body morphology of the Least Killifish, Heterandria formosa

An important step in diagnosing local adaptation is the demonstration that phenotypic variation among populations is at least in part genetically based. To do this, many methods experimentally minimize the environmental effect on the phenotype to elucidate the genetic effect. Minimizing the environmental effect often includes reducing possible environmental maternal effects. However, maternal effects can be an important factor in patterns of local adaptation as well as adaptive plasticity. Here, we report the results of an experiment with males from two populations of the poeciliid fish, Heterandria formosa, designed to examine the relative influence of environmental maternal effects and environmental effects experienced during growth and development on body morphology, and, in addition, whether the balance among those effects is unique to each population. We used a factorial design that varied thermal environment and water chemistry experienced by mothers and thermal environment and water chemistry experienced by offspring. We found substantial differences between the two populations in their maternal and offspring norms of reaction of male body morphology to differences in thermal environment and water chemistry. We also found that the balance between maternal effects and postparturition environmental effects differed from one thermal regime to another and among traits. These results indicate that environmental maternal effects can be decidedly population-specific and, as a result, might either contribute to the appearance of or blur evidence for local adaptation. These results also suggest that local adaptation might also occur through the evolution of maternal norms of reaction to important, and varying, environmental factors.

A new mechanistic understanding of light-limitation in the seagrass Zostera muelleri

In this study we investigated the effect of light-limitation (∼20 μmol photons m^-2 s^-1) on the southern hemisphere seagrass, Zostera muelleri. RNA sequencing, chlorophyll fluorometry and HPLC techniques were used to investigate how the leaf-specific transcriptome drives changes in photosynthesis and photo-pigments in Z. muelleri over 6 days. 1593 (7.51%) genes were differentially expressed on day 2 and 1481 (6.98%) genes were differentially expressed on day 6 of the experiment. Differential gene expression correlated with significant decreases in rETR_Max, I_k, an increase in Yi (initial photosynthetic quantum yield of photosystem II), and significant changes in pigment composition. Regulation of carbohydrate metabolism was observed along with evidence that abscisic acid may serve a role in the low-light response of this seagrass. This study provides a novel understanding of how Z. muelleri responds to light-limitation in the marine water column and provides potential molecular markers for future conservation monitoring efforts.

SeagrassDB: An open-source transcriptomics landscape for phylogenetically profiled seagrasses and aquatic plants

Seagrasses and aquatic plants are important clades of higher plants, significant for carbon sequestration and marine ecological restoration. They are valuable in the sense that they allow us to understand how plants have developed traits to adapt to high salinity and photosynthetically challenged environments. Here, we present a large-scale phylogenetically profiled transcriptomics repository covering seagrasses and aquatic plants. SeagrassDB encompasses a total of 1,052,262 unigenes with a minimum and maximum contig length of 8,831 bp and 16,705 bp respectively. SeagrassDB provides access to 34,455 transcription factors, 470,568 PFAM domains, 382,528 prosite models and 482,121 InterPro domains across 9 species. SeagrassDB allows for the comparative gene mining using BLAST-based approaches and subsequent unigenes sequence retrieval with associated features such as expression (FPKM values), gene ontologies, functional assignments, family level classification, Interpro domains, KEGG orthology (KO), transcription factors and prosite information. SeagrassDB is available to the scientific community for exploring the functional genic landscape of seagrass and aquatic plants at: http://115.146.91.129/index.php.

Antioxidant response to heat stress in seagrasses. A gene expression study

Seawater warming associated to the ongoing climate change threatens functioning and survival of keystone coastal benthic species such as seagrasses. Under elevated temperatures, the production of reactive oxygen species (ROS) is increased and plants must activate their antioxidant defense mechanisms to protect themselves from oxidative damage. Here we explore from a molecular perspective the ability of Mediterranean seagrasses to activate heat stress response mechanisms, with particular focus on antioxidants. The level of expression of targeted genes was analyzed in shallow and deep plants of the species Posidonia oceanica and in shallow plants of Cymodocea nodosa along an acute heat exposure of several days and after recovery. The overall gene expression response of P. oceanica was more intense and complete than in C. nodosa and reflected a higher oxidative stress level during the experimental heat exposure. The strong activation of genes with chaperone activity (heat shock proteins and a luminal binding protein) just in P. oceanica plants, suggested the higher sensitivity of the species to increased temperatures. In spite of the interspecific differences, genes from the superoxide dismutase (SOD) family seem to play a pivotal role in the thermal stress response of Mediterranean seagrasses as previously reported for other marine plant species. Shallow and deep P. oceanica ecotypes showed a different timing of response to heat. Shallow plants early responded to heat and after a few days relaxed their response which suggests a successful early metabolic adjustment. The response of deep plants was delayed and their recovery incomplete evidencing a lower resilience to heat in respect to shallow ecotypes. Moreover, shallow ecotypes showed some degree of pre-adaptation to heat as most analyzed genes showed higher constitutive expression levels than in deep ecotypes. The recurrent exposure of shallow plants to elevated summer temperatures has likely endowed them with a higher basal level of antioxidant defense and a faster responsiveness to warming than deep plants. Our findings match with previous physiological studies and supported the idea that warming will differently impact Mediterranean seagrass meadows depending on the species as well as on the depth (i.e. thermal regimen) at which the meadow grows. The increase in the incidence of summer heat waves could therefore produce a significant change in the distribution and composition of Mediterranean seagrass meadows with considerable consequences for the functioning of the whole ecosystem and for the socio-economic services that these ecosystems offer to the riverine populations.

Molecular Mechanisms behind the Physiological Resistance to Intense Transient Warming in an Iconic Marine Plant

The endemic Mediterranean seagrass Posidonia oceanica is highly threatened by the increased frequency and intensity of heatwaves. Meadows of the species offer a unique opportunity to unravel mechanisms marine plants activate to cope transient warming, since their wide depth distribution impose divergent heat-tolerance. Understanding these mechanisms is imperative for their conservation. Shallow and deep genotypes within the same population were exposed to a simulated heatwave in mesocosms, to analyze their transcriptomic and photo-physiological responses during and after the exposure. Shallow plants, living in a more unstable thermal environment, optimized phenotype variation in response to warming. These plants showed a pre-adaptation of genes in anticipation of stress. Shallow plants also showed a stronger activation of heat-responsive genes and the exclusive activation of genes involved in epigenetic mechanisms and in molecular mechanisms that are behind their higher photosynthetic stability and respiratory acclimation. Deep plants experienced higher heat-induced damage and activated metabolic processes for obtaining extra energy from sugars and amino acids, likely to support the higher protein turnover induced by heat. In this study we identify transcriptomic mechanisms that may facilitate persistence of seagrasses to anomalous warming events and we discovered that P. oceanica plants from above and below the mean depth of the summer thermocline have differential resilience to heat.