Finding genes and pathways that underlie coral adaptation.

Mass coral bleaching is one of the clearest threats of climate change to the persistence of marine biodiversity. Despite the negative impacts of bleaching on coral health and survival, some corals may be able to rapidly adapt to warming ocean temperatures. Thus, a significant focus in coral research is identifying the genes and pathways underlying coral heat adaptation. Here, we review state-of-the-art methods that may enable the discovery of heat-adaptive loci in corals and identify four main knowledge gaps. To fill these gaps, we describe an experimental approach combining seascape genomics with CRISPR/Cas9 gene editing to discover and validate heat-adaptive loci. Finally, we discuss how information on adaptive genotypes could be used in coral reef conservation and management strategies.

Vertical distribution and seasonal dynamics of planktonic cyanobacteria communities in a water column of deep mesotrophic Lake Geneva.

Background: Temperate subalpine lakes recovering from eutrophication in central Europe are experiencing harmful blooms due to the proliferation of Planktothrix rubescens, a potentially toxic cyanobacteria. To optimize the management of cyanobacteria blooms there is the need to better comprehend the combination of factors influencing the diversity and dominance of cyanobacteria and their impact on the lake's ecology. The goal of this study was to characterize the diversity and seasonal dynamics of cyanobacteria communities found in a water column of Lake Geneva, as well as the associated changes on bacterioplankton abundance and composition. Methods: We used 16S rRNA amplicon high throughput sequencing on more than 200 water samples collected from surface to 100 meters deep monthly over 18 months. Bacterioplankton abundance was determined by quantitative PCR and PICRUSt predictions were used to explore the functional pathways present in the community and to calculate functional diversity indices. Results: The obtained results confirmed that the most dominant cyanobacteria in Lake Geneva during autumn and winter was Planktothrix (corresponding to P. rubescens). Our data also showed an unexpectedly high relative abundance of picocyanobacterial genus Cyanobium, particularly during summertime. Multidimensional scaling of Bray Curtis dissimilarity revealed that the dominance of P. rubescens was coincident with a shift in the bacterioplankton community composition and a significant decline in bacterioplankton abundance, as well as a temporary reduction in the taxonomic and PICRUSt2 predicted functional diversity. Conclusion: Overall, this study expands our fundamental understanding of the seasonal dynamics of cyanobacteria communities along a vertical column in Lake Geneva and the ecology of P. rubescens, ultimately contributing to improve our preparedness against the potential occurrence of toxic blooms in the largest lake of western Europe.

Seascape genomics reveals candidate molecular targets of heat stress adaptation in three coral species.

Anomalous heat waves are causing a major decline of hard corals around the world and threatening the persistence of coral reefs. There are, however, reefs that have been exposed to recurrent thermal stress over the years and whose corals appear to have been tolerant against heat. One of the mechanisms that could explain this phenomenon is local adaptation, but the underlying molecular mechanisms are poorly known. In this work, we applied a seascape genomics approach to study heat stress adaptation in three coral species of New Caledonia (southwestern Pacific) and to uncover the molecular actors potentially involved. We used remote sensing data to characterize the environmental trends across the reef system, and sampled corals living at the most contrasted sites. These samples underwent next generation sequencing to reveal single nucleotide polymorphisms (SNPs), frequencies of which were associated with heat stress gradients. As these SNPs might underpin an adaptive role, we characterized the functional roles of the genes located in their genomic region. In each of the studied species, we found heat stress-associated SNPs located in proximity of genes involved in pathways well known to contribute to the cellular responses against heat, such as protein folding, oxidative stress homeostasis, inflammatory and apoptotic pathways, and DNA damage-repair. In some cases, the same candidate molecular targets of heat stress adaptation recurred among species. Together, these results underline the relevance and the power of the seascape genomics approach for the discovery of adaptive traits that could allow corals to persist across wider thermal ranges.

Coral cover surveys corroborate predictions on reef adaptive potential to thermal stress.

As anomalous heat waves are causing the widespread decline of coral reefs worldwide, there is an urgent need to identify coral populations tolerant to thermal stress. Heat stress adaptive potential is the degree of tolerance expected from evolutionary processes and, for a given reef, depends on the arrival of propagules from reefs exposed to recurrent thermal stress. For this reason, assessing spatial patterns of thermal adaptation and reef connectivity is of paramount importance to inform conservation strategies. In this work, we applied a seascape genomics framework to characterize the spatial patterns of thermal adaptation and connectivity for coral reefs of New Caledonia (Southern Pacific). In this approach, remote sensing of seascape conditions was combined with genomic data from three coral species. For every reef of the region, we computed a probability of heat stress adaptation, and two indices forecasting inbound and outbound connectivity. We then compared our indicators to field survey data, and observed that decrease of coral cover after heat stress was lower at reefs predicted with high probability of adaptation and inbound connectivity. Last, we discussed how these indicators can be used to inform local conservation strategies and preserve the adaptive potential of New Caledonian reefs.

Seascape genomics as a new tool to empower coral reef conservation strategies: An example on north-western Pacific Acropora digitifera.

Coral reefs are suffering a major decline due to the environmental constraints imposed by climate change. Over the last 20 years, three major coral bleaching events occurred in concomitance with anomalous heatwaves, provoking a severe loss of coral cover worldwide. The conservation strategies for preserving reefs, as they are implemented now, cannot cope with global climatic shifts. Consequently, researchers are advocating for preservation networks to be set-up to reinforce coral adaptive potential. However, the main obstacle to this implementation is that studies on coral adaption are usually hard to generalize at the scale of a reef system. Here, we study the relationships between genotype frequencies and environmental characteristics of the sea (seascape genomics), in combination with connectivity analysis, to investigate the adaptive potential of a flagship coral species of the Ryukyu Archipelago (Japan). By associating genotype frequencies with descriptors of historical environmental conditions, we discovered six genomic regions hosting polymorphisms that might promote resistance against heat stress. Remarkably, annotations of genes in these regions were consistent with molecular roles associated with heat responses. Furthermore, we combined information on genetic and spatial distances between reefs to predict connectivity at a regional scale. The combination of these results portrayed the adaptive potential of this population: we were able to identify reefs carrying potential heat stress adapted genotypes and to understand how they disperse to neighbouring reefs. This information was summarized by objective, quantifiable and mappable indices covering the whole region, which can be extremely useful for future prioritization of reefs in conservation planning. This framework is transferable to any coral species on any reef system and therefore represents a valuable tool for empowering preservation efforts dedicated to the protection of coral reefs in warming oceans.

Correction to: Sex-specific changes in gene expression in response to estrogen pollution around the onset of sex differentiation in grayling (Salmonidae).

An amendment to this paper has been published and can be accessed via the original article.

Sampling strategy optimization to increase statistical power in landscape genomics: A simulation-based approach.

An increasing number of studies are using landscape genomics to investigate local adaptation in wild and domestic populations. Implementation of this approach requires the sampling phase to consider the complexity of environmental settings and the burden of logistical constraints. These important aspects are often underestimated in the literature dedicated to sampling strategies. In this study, we computed simulated genomic data sets to run against actual environmental data in order to trial landscape genomics experiments under distinct sampling strategies. These strategies differed by design approach (to enhance environmental and/or geographical representativeness at study sites), number of sampling locations and sample sizes. We then evaluated how these elements affected statistical performances (power and false discoveries) under two antithetical demographic scenarios. Our results highlight the importance of selecting an appropriate sample size, which should be modified based on the demographic characteristics of the studied population. For species with limited dispersal, sample sizes above 200 units are generally sufficient to detect most adaptive signals, while in random mating populations this threshold should be increased to 400 units. Furthermore, we describe a design approach that maximizes both environmental and geographical representativeness of sampling sites and show how it systematically outperforms random or regular sampling schemes. Finally, we show that although having more sampling locations (between 40 and 50 sites) increase statistical power and reduce false discovery rate, similar results can be achieved with a moderate number of sites (20 sites). Overall, this study provides valuable guidelines for optimizing sampling strategies for landscape genomics experiments.

Sex-specific changes in gene expression in response to estrogen pollution around the onset of sex differentiation in grayling (Salmonidae).

The synthetic 17α-ethinylestradiol (EE2) is a common estrogenic pollutant that has been suspected to affect the demography of river-dwelling salmonids. One possibility is that exposure to EE2 tips the balance during initial steps of sex differentiation, so that male genotypes show female-specific gene expression and gonad formation. Here we study EE2 effects on gene expression around the onset of sex differentiation in a population of European grayling (Thymallus thymallus) that suffers from sex ratio distortions. We exposed singly-raised embryos to one dose of 1 ng/L EE2, studied gene expression 10 days before hatching, at the day of hatching, and around the end of the yolk-sac stage, and related it to genetic sex (sdY genotype). We found that exposure to EE2 affects expression of a large number of genes, especially around hatching. These effects were strongly sex-dependent. We then raised fish for several months after hatching and found no evidence of sex reversal in the EE2-exposed fish. We conclude that ecologically relevant (i.e. low) levels of EE2 pollution do not cause sex reversal by simply tipping the balance at early stages of sex differentiation, but that they interfere with sex-specific gene expression.

Rapid identification and interpretation of gene-environment associations using the new R.SamBada landscape genomics pipeline.

samßada is a genome-environment association software, designed to search for signatures of local adaptation. However, pre- and postprocessing of data can be labour-intensive, preventing wider uptake of the method. We have now developed R.SamBada, an r-package providing a pipeline for landscape genomic analysis based on samßada, spanning from the retrieval of environmental conditions at sampling locations to gene annotation using the Ensembl genome browser. As a result, R.SamBada standardizes the landscape genomics pipeline and eases the search for candidate genes of local adaptation, enhancing reproducibility of landscape genomic studies. The efficiency and power of the pipeline is illustrated using two examples: sheep populations from Morocco with no evident population structure and Lidia cattle from Spain displaying population substructuring. In both cases, R.SamBada enabled rapid identification and interpretation of candidate genes, which are further discussed in the light of local adaptation. The package is available in the r CRAN package repository and on GitHub (github.com/SolangeD/R.SamBada).

Sex differentiation in grayling (Salmonidae) goes through an all-male stage and is delayed in genetic males who instead grow faster.

Fish populations can be threatened by distorted sex ratios that arise during sex differentiation. Here we describe sex differentiation in a wild grayling (Thymallus thymallus) population that suffers from distorted sex ratios. We verified that sex determination is linked to the sex determining locus (sdY) of salmonids. This allowed us to study sex-specific gene expression and gonadal development. Sex-specific gene expression could be observed during embryogenesis and was strong around hatching. About half of the fish showed immature testes around eleven weeks after fertilization. This phenotype was mostly replaced by the "testis-to-ovary" or "ovaries" phenotypes during development. The gonads of the remaining fish stayed undifferentiated until six months after fertilization. Genetic sexing revealed that fish with undifferentiated gonads were all males, who grew larger than the genetic females during the observational period. Only 12% of the genetic males showed testicular tissue six months after fertilization. We conclude that sex differentiation starts before hatching, goes through an all-male stage for both sexes (which represents a rare case of "undifferentiated" gonochoristic species that usually go through an all-female stage), and is delayed in males. During these juvenile stages males grow faster than females instead of developing their gonads.

The Genome of the Toluene-Degrading Pseudomonas veronii Strain 1YdBTEX2 and Its Differential Gene Expression in Contaminated Sand.

The natural restoration of soils polluted by aromatic hydrocarbons such as benzene, toluene, ethylbenzene and m- and p-xylene (BTEX) may be accelerated by inoculation of specific biodegraders (bioaugmentation). Bioaugmentation mainly involves introducing bacteria that deploy their metabolic properties and adaptation potential to survive and propagate in the contaminated environment by degrading the pollutant. In order to better understand the adaptive response of cells during a transition to contaminated material, we analyzed here the genome and short-term (1 h) changes in genome-wide gene expression of the BTEX-degrading bacterium Pseudomonas veronii 1YdBTEX2 in non-sterile soil and liquid medium, both in presence or absence of toluene. We obtained a gapless genome sequence of P. veronii 1YdBTEX2 covering three individual replicons with a total size of 8 Mb, two of which are largely unrelated to current known bacterial replicons. One-hour exposure to toluene, both in soil and liquid, triggered massive transcription (up to 208-fold induction) of multiple gene clusters, such as toluene degradation pathway(s), chemotaxis and toluene efflux pumps. This clearly underlines their key role in the adaptive response to toluene. In comparison to liquid medium, cells in soil drastically changed expression of genes involved in membrane functioning (e.g., lipid composition, lipid metabolism, cell fatty acid synthesis), osmotic stress response (e.g., polyamine or trehalose synthesis, uptake of potassium) and putrescine metabolism, highlighting the immediate response mechanisms of P. veronii 1YdBTEX2 for successful establishment in polluted soil.