Characterization of Genetic and Epigenetic Variation in Sperm and Red Blood Cells from Adult Hatchery and Natural-Origin Steelhead, Oncorhynchus mykiss

Single generation epigenetic change in captivity and reinforcement in subsequent generations in a delta smelt (Hypomesus transpacificus) conservation hatchery

A refugial population of the endangered delta smelt (Hypomesus transpacificus) has been maintained at the Fish Conservation and Culture Laboratory (FCCL) at UC Davis since 2008. Despite intense genetic management, fitness differences between wild and cultured fish have been observed at the FCCL. To investigate the molecular underpinnings of hatchery domestication, we used whole-genome bisulfite sequencing to quantify epigenetic differences between wild and hatchery-origin delta smelt. Differentially methylated regions (DMRs) were identified from 104 individuals by comparing the methylation patterns in different generations of hatchery fish (G1, G2, G3) with their wild parents (G0). We discovered a total of 132 significant DMRs (p < .05) between G0 and G1, 132 significant DMRs between G0 and G2, and 201 significant DMRs between G0 and G3. Our results demonstrate substantial differences in methylation patterns emerged between the wild and hatchery-reared fish in the early generations in the hatchery, with a higher proportion of hypermethylated DMRs in hatchery-reared fish. The rearing environment was found to be a stronger predictor of individual clustering based on methylation patterns than family, sex or generation. Our study indicates a reinforcement of the epigenetic status with successive generations in the hatchery environment, as evidenced by an increase in methylation in hypermethylated DMRs and a decrease in methylation in hypomethylated DMRs over time. Lastly, our results demonstrated heterogeneity in inherited methylation pattern in families across generations. These insights highlight the long-term consequences of hatchery practices on the epigenetic landscape, potentially impacting wild fish populations.

Strengths and limitations of reduced representation bisulfite sequencing (RRBS) in the perspective of DNA methylation analysis in fish: a case-study on rainbow trout spermatozoa

DNA methylation in CpG dinucleotides is an important epigenetic mark in fish spermatozoa since it has been shown that some sperm methylome features are transmitted to the offspring. Reduced representation bisulfite sequencing (RRBS) is one genome-scale methods developed to assess DNA methylation at CpG sites. It allows the sequencing of a reduced fraction of the genome expected to be enriched in CpGs. The aim of this study is to characterize the extent of the CpG sites that can be identified in the RRBS-reduced sequenced fraction of rainbow trout spermatozoa, in order to evaluate the potential of RRBS for sperm DNA methylation studies. We observed that RRBS did provide a reduced amount of genomic data, the sum of the CpGs analyzed on 12 males spanning 9% of the total genomic CpGs. CpGs were only slightly enriched in the RRBS data (×1.7 times the sequenced nucleotides), the possible causes being linked to trout genome structure and sequenced fragments size. All genomic functional features were represented in our CpG dataset, with a noticeable enrichment in exons but, strikingly, not in promoters. The number of CpGs shared between biological replicates was low, but this proportion reached workable values from six biological replicates (46% of the analyzed cytosines) on. The choices that are to be made regarding fragment size selection and the options during bioinformatic data processing are discussed. In all, RRBS is a relevant first-approach method to scan the CpG DNA methylation status of spermatozoa along rainbow trout genome, although in a very reduced pattern among biological replicates.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Epigenetic and Genetic Differentiation Between Coregonus Species Pairs

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation.

Cryopreservation effect on DNA methylation profile in rainbow trout spermatozoa

Spermatozoa are the cells that are most commonly used for cryopreservation of valuable genetic resources in aquaculture. It is known that fish spermatozoa transmit to the embryo not only their genetic but also their epigenetic profile, especially DNA methylation. Therefore, any alteration of the DNA methylation profile in spermatozoa induces the risk of transmitting epigenetic alterations to the offspring. The aim of this study was to assess the effect of cryopreservation on DNA methylation in rainbow trout spermatozoa. To trigger variable cellular response after freezing-thawing, spermatozoa from mature males were cryopreserved with dimethyl sulfoxide, methanol or glycerol as cryoprotectant. We observed that dimethyl sulfoxide was the best to preserve thawed spermatozoa functions. Methanol only slightly preserved all the cellular parameters, while glycerol failed to protect motility and fertilization ability. The consequences on DNA methylation were assessed using Reduced Representation Bisulfite Sequencing (RRBS). Sperm cryopreservation did not thoroughly impact DNA methylation, although 335-564 differentially methylated cytosines were characterized depending on the cryoprotectant. Very few of them were shared between cryoprotectants, and no correlation with the extent of cellular damage was found. Our study showed that DNA methylation was only slightly altered after sperm cryopreservation, and this may render further analysis of the risk for the progeny very challenging.

Sperm quality parameters, fertilizing potential, metabolites, and DNA methylation in cold-stored and cryopreserved milt from Atlantic salmon (Salmo salar L.)

Cold storage and freezing/thawing of milt may affect sperm functionality and the subsequent fertilization ability of milt. This study aimed to investigate sperm quality parameters and fertilization potential of Atlantic salmon milt, stored cold and subsequently cryopreserved, using different storage conditions. The objective was also to assess if analysis of milt metabolites and sperm DNA methylation signatures could be applicable to further elucidate sperm quality and fertilization following preservation. Milt samples were collected from eight mature Atlantic salmon males and stored for 4 days at 2°C and 8°C. Samples were taken on day one of storage at 2°C and on day four of storage at 2°C and 8°C. Storage for 4 days at 8°C is expected to be detrimental to sperm quality, and was included to create contrasts. Correspondingly, aliquots of cold-stored milt were prepared for cryopreservation, resulting in a total of six experimental conditions. Samples from all six experimental conditions were used in fertilization trials and analyzed for sperm viability, motility, ATP content, DNA fragmentation index, and High DNA stainability. In addition, milt samples from four of the males were analyzed for targeted metabolites and DNA methylation signatures by reduced representation bisulfite sequencing. The fertilization trials were performed using sperm:egg ratios of 75 × 103 and 500 × 103, respectively. Storage duration, temperature, and cryopreservation of cold-stored milt influenced several sperm quality parameters, metabolites, and DNA methylation signatures. The total motility, progressive motility, ATP, and velocity parameters were the sperm parameters with the strongest correlation to fertilization rates (p < 0.01). Several metabolites were correlated with fertility rates in both cold-stored and cryopreserved samples (p < 0.05). The fertilizing capacity of cold-stored milt was significantly reduced after 4 days of storage at 8°C, while corresponding cryopreserved milt showed reduced fertilization at both storage temperatures (2°C and 8°C) (p < 0.05). The results indicate that cryopreservation of milt stored for 1 day does not compromise either fertilization ability or DNA methylation signatures.

Captive rearing effects on the methylome of Atlantic salmon after oceanic migration: Sex-specificity and intergenerational stability

Captive rearing in salmon hatcheries can have considerable impacts on both fish phenotype and fitness within a single generation, even in the absence of genetic change. Evidence for hatchery-induced changes in DNA methylation is becoming abundant, though questions remain on the sex-specificity of these effects, their persistence until spawning and potential for transmission to future generations. Here we performed whole genome methylation sequencing of fin tissue for 16 hatchery and 16 wild Atlantic salmon (Salmo salar) returning to spawn in the Rimouski River, Québec, Canada. We identified two cohorts of hatchery-reared salmon through methylation analysis, one of which was epigenetically similar to wild fish, suggesting that supplementation efforts may be able to minimize the epigenetic effects of hatchery rearing. We found considerable sex-specific effects of hatchery rearing, with few genomic regions being affected in both males and females. We also analysed the methylome of 32 F1 offspring from four groups (pure wild, pure hatchery origin and reciprocal hybrids). We found that few epigenetic changes due to parental hatchery rearing persisted in the F1 offspring though the patterns of inheritance appear to be complex, involving nonadditive effects. Our results suggest that the epigenetic effects of hatchery rearing can be minimal in F0 . There may also be minimal epigenetic inheritance and rapid loss of epigenetic changes associated with hatchery rearing. However, due to sex-specificity and nonadditive patterns of inheritance, methylation changes due to captive rearing are rather complex and the field would benefit from further research on minimizing the epigenetic effects of captive rearing in conservation efforts.

Genome-wide DNA methylation predicts environmentally driven life history variation in a marine fish

Epigenetic modifications are thought to be one of the molecular mechanisms involved in plastic adaptive responses to environmental variation. However, studies reporting associations between genome-wide epigenetic changes and habitat-specific variations in life history traits (e.g., lifespan, reproduction) are still scarce, likely due to the recent application of methylome resequencing methods to non-model species. In this study, we examined associations between whole genome DNA methylation and environmentally driven life history variation in 2 lineages of a marine fish, the capelin (Mallotus villosus), from North America and Europe. In both lineages, capelin harbor 2 contrasting life history tactics (demersal vs. beach-spawning). Performing whole genome and methylome sequencing, we showed that life history tactics are associated with epigenetic changes in both lineages, though the effect was stronger in European capelin. Genetic differentiation between the capelin harboring different life history tactics was negligible, but we found genome-wide methylation changes in both lineages. We identified 9,125 European and 199 North American differentially methylated regions (DMRs) due to life history. Gene ontology (GO) enrichment analysis for both lineages revealed an excess of terms related to neural function. Our results suggest that environmental variation causes important epigenetic changes that are associated with contrasting life history tactics in lineages with divergent genetic backgrounds, with variable importance of genetic variation in driving epigenetic variation. Our study emphasizes the potential role of genome-wide epigenetic variation in adaptation to environmental variation.

Population Epigenetics: The Extent of DNA Methylation Variation in Wild Animal Populations

Population epigenetics explores the extent of epigenetic variation and its dynamics in natural populations encountering changing environmental conditions. In contrast to population genetics, the basic concepts of this field are still in their early stages, especially in animal populations. Epigenetic variation may play a crucial role in phenotypic plasticity and local adaptation as it can be affected by the environment, it is likely to have higher spontaneous mutation rate than nucleotide sequences do, and it may be inherited via non-mendelian processes. In this review, we aim to bring together natural animal population epigenetic studies to generate new insights into ecological epigenetics and its evolutionary implications. We first provide an overview of the extent of DNA methylation variation and its autonomy from genetic variation in wild animal population. Second, we discuss DNA methylation dynamics which create observed epigenetic population structures by including basic population genetics processes. Then, we highlight the relevance of DNA methylation variation as an evolutionary mechanism in the extended evolutionary synthesis. Finally, we suggest new research directions by highlighting gaps in the knowledge of the population epigenetics field. As for our results, DNA methylation diversity was found to reveal parameters that can be used to characterize natural animal populations. Some concepts of population genetics dynamics can be applied to explain the observed epigenetic structure in natural animal populations. The set of recent advancements in ecological epigenetics, especially in transgenerational epigenetic inheritance in wild animal population, might reshape the way ecologists generate predictive models of the capacity of organisms to adapt to changing environments.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus)

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10 000 years, offers an interesting ecological model to address such effects. Using bisulfite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic - limnetic axis), both at single CpG sites and in 1,000bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

Fish as Model Systems to Study Epigenetic Drivers in Human Self-Domestication and Neurodevelopmental Cognitive Disorders

Modern humans exhibit phenotypic traits and molecular events shared with other domesticates that are thought to be by-products of selection for reduced aggression. This is the human self-domestication hypothesis. As one of the first types of responses to a novel environment, epigenetic changes may have also facilitated early self-domestication in humans. Here, we argue that fish species, which have been recently domesticated, can provide model systems to study epigenetic drivers in human self-domestication. To test this, we used in silico approaches to compare genes with epigenetic changes in early domesticates of European sea bass with genes exhibiting methylation changes in anatomically modern humans (comparison 1), and neurodevelopmental cognitive disorders considered to exhibit abnormal self-domestication traits, i.e., schizophrenia, Williams syndrome, and autism spectrum disorders (comparison 2). Overlapping genes in comparison 1 were involved in processes like limb morphogenesis and phenotypes like abnormal jaw morphology and hypopigmentation. Overlapping genes in comparison 2 affected paralogue genes involved in processes such as neural crest differentiation and ectoderm differentiation. These findings pave the way for future studies using fish species as models to investigate epigenetic changes as drivers of human self-domestication and as triggers of cognitive disorders.

Reduced relative fitness in hatchery‐origin Pink Salmon in two streams in Prince William Sound, Alaska

Previous studies generally report that hatchery‐origin Pacific Salmon (Oncorhynchus spp.) have lower relative reproductive success (RRS) than their natural‐origin counterparts. We estimated the RRS of Pink Salmon (O. gorbuscha) in Prince William Sound (PWS), Alaska using incomplete pedigrees. In contrast to other RRS studies, Pink Salmon have a short freshwater life history, freshwater habitats in PWS are largely unaltered by development, and sampling was conducted without the aid of dams or weirs resulting in incomplete sampling of spawning individuals. Pink Salmon released from large‐scale hatchery programs in PWS have interacted with wild populations for more than 15 generations. Hatchery populations were established from PWS populations but have subsequently been managed as separate broodstocks. Gene flow is primarily directional, from hatchery strays to wild populations. We used genetic‐based parentage analysis to estimate the RRS of a single generation of stray hatchery‐origin Pink Salmon in two streams, and across the odd‐ and even‐year lineages. Despite incomplete sampling, we assigned 1745 offspring to at least one parent. Reproductive success (RS), measured as sampled adult offspring that returned to their natal stream, was significantly lower for hatchery‐ vs. natural‐origin parents in both lineages, with RRS ranging from 0.03 to 0.47 for females and 0.05 to 0.86 for males. Generalized linear modeling for the even‐year lineage indicated that RRS was lower for hatchery‐origin fish, ranging from 0.42 to 0.60, after accounting for sample date (run timing), sample location within the stream, and fish length. Our results strongly suggest that hatchery‐origin strays have lower fitness in the wild. The consequences of reduced RRS on wild productivity depend on whether the mechanisms underlying reduced RRS are environmentally driven, and likely ephemeral, or genetically driven, and likely persistent across generations.

Early fish domestication affects methylation of key genes involved in the rapid onset of the farmed phenotype

Animal domestication is a process of environmental modulation and artificial selection leading to permanent phenotypic modifications. Recent studies showed that phenotypic changes occur very early in domestication, i.e., within the first generation in captivity, which raises the hypothesis that epigenetic mechanisms may play a critical role on the early onset of the domestic phenotype. In this context, we applied reduced representation bisulphite sequencing to compare methylation profiles between wild Nile tilapia females and their offspring reared under farmed conditions. Approximately 700 differentially methylated CpG sites were found, many of them associated not only with genes involved in muscle growth, immunity, autophagy and diet response but also related to epigenetic mechanisms, such as RNA methylation and histone modifications. This bottom-up approach showed that the phenotypic traits often related to domestic animals (e.g., higher growth rate and different immune status) may be regulated epigenetically and prior to artificial selection on gene sequences. Moreover, it revealed the importance of diet in this process, as reflected by differential methylation patterns in genes critical to fat metabolism. Finally, our study highlighted that the TGF-β1 signalling pathway may regulate and be regulated by several differentially methylated CpG-associated genes. This could be an important and multifunctional component in promoting adaptation of fish to a domestic environment while modulating growth and immunity-related traits.

Generational effects of a chronic exposure to a low environmentally relevant concentration of glyphosate on rainbow trout, Oncorhynchus mykiss

In the past few decades, glyphosate became the most used herbicide substance worldwide. As a result, the substance is ubiquitous in surface waters. Concerns have been raised about its ecotoxicological impact, but little is known about its generational toxicity. In this study, we investigate the impact of an environmentally relevant concentration of glyphosate and its co-formulants on an F2 generation issued from exposed generations F0 and F1. Trans, inter and multigenerational toxicity of 1 μgL-1 of the active substance was evaluated on early stages of development and juvenile rainbow trout (Oncorhynchus mykiss) using different molecular, biochemical, immuno-hematologic, and biometric parameters, behavior analysis, and a viral challenge. Reproductive parameters of generation F1 were not affected. However, developmental toxicity in generation F2 due to glyphosate alone or co-formulated was observed with head size changes (e.g. head surface up to +10%), and metabolic disruptions (e.g. 35% reduction in cytochrome-c-oxidase). Moreover, larvae exposed transgenerationally to Viaglif and intergenerationally to glyphosate and Roundup presented a reduced response to light, potentially indicating altered escape behavior. Overall methylation was, however, not altered and further experiments using gene-specific DNA metylation analyses are required. After several months, biochemical parameters measured in juvenile fish were no longer impacted, only intergenerational exposure to glyphosate drastically increased the susceptibility of rainbow trout to hematopoietic necrosis virus. This result might be due to a lower antibody response in exposed fish. In conclusion, our results show that generational exposure to glyphosate induces developmental toxicity and increases viral susceptibility. Co-formulants present in glyphosate-based herbicides can modulate the toxicity of the active substance. Further investigations are required to study the specific mechanisms of transmission but our results suggest that both non-genetic mechanisms and exposure during germinal stage could be involved.

Environment-driven reprogramming of gamete DNA methylation occurs during maturation and is transmitted intergenerationally in Atlantic Salmon

An epigenetic basis for transgenerational plasticity in animals is widely theorized, but convincing empirical support is limited by taxa-specific differences in the presence and role of epigenetic mechanisms. In teleost fishes, DNA methylation generally does not undergo extensive reprogramming and has been linked with environmentally induced intergenerational effects, but solely in the context of early life environmental differences. Using whole-genome bisulfite sequencing, we demonstrate that differential methylation of sperm occurs in response to captivity during the maturation of Atlantic Salmon (Salmo salar), a species of major economic and conservation significance. We show that adult captive exposure further induces differential methylation in an F1 generation that is associated with fitness-related phenotypic differences. Some genes targeted with differential methylation were consistent with genes differential methylated in other salmonid fishes experiencing early-life hatchery rearing, as well as genes under selection in domesticated species. Our results support a mechanism of transgenerational plasticity mediated by intergenerational inheritance of DNA methylation acquired late in life for salmon. To our knowledge, this is the first-time environmental variation experienced later in life has been directly demonstrated to influence gamete DNA methylation in fish.

Epigenomic modifications induced by hatchery rearing persist in germ line cells of adult salmon after their oceanic migration

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human-induced genomic variation, human-induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild-born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine-scale analyses (local and chromosomal) allows the detection of parallel early-life hatchery-induced epimarks that differentiate wild from hatchery-reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild-born Coho salmon. Moreover, those early-life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early-life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

Rapid adaptation through genomic and epigenomic responses following translocations in an endangered salmonid

Identifying the molecular mechanisms facilitating adaptation to new environments is a key question in evolutionary biology, especially in the face of current rapid and human-induced changes. Translocations have become an important tool for species conservation, but the attendant small population sizes and new ecological pressures might affect phenotypic and genotypic variation and trajectories dramatically and in unknown ways. In Scotland, the European whitefish (Coregonus lavaretus) is native to only two lakes and vulnerable to extirpation. Six new refuge populations were established over the last 30 years as a conservation measure. In this study, we examined whether there is a predictable ecological and evolutionary response of these fishes to translocation. We found eco-morphological differences, as functional traits relating to body shape differed between source and refuge populations. Dual isotopic analyses suggested some ecological release, with the diets in refuge populations being more diverse than in source populations. Analyses of up to 9117 genome-mapped SNPs showed that refuge populations had reduced genetic diversity and elevated inbreeding and relatedness relative to source populations, though genomic differentiation was low (F ST = 0.002-0.030). We identified 14 genomic SNPs that showed shared signals of a selective response to translocations, including some located near or within genes involved in the immune system, nervous system and hepatic functions. Analysis of up to 120,897 epigenomic loci identified a component of consistent differential methylation between source and refuge populations. We found that epigenomic variation and genomic variation were associated with morphological variation, but we were not able to infer an effect of population age because the patterns were also linked with the methodology of the translocations. These results show that conservation-driven translocations affect evolutionary potential by impacting eco-morphological, genomic and epigenomic components of diversity, shedding light on acclimation and adaptation process in these contexts.

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene-targeted PCR-based assay for next-generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single-gene Mantel tests comparing methylation and pairwise F ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.

Population Genomics Training for the Next Generation of Conservation Geneticists: ConGen 2018 Workshop

The increasing availability and complexity of next-generation sequencing (NGS) data sets make ongoing training an essential component of conservation and population genetics research. A workshop entitled “ConGen 2018” was recently held to train researchers in conceptual and practical aspects of NGS data production and analysis for conservation and ecological applications. Sixteen instructors provided helpful lectures, discussions, and hands-on exercises regarding how to plan, produce, and analyze data for many important research questions. Lecture topics ranged from understanding probabilistic (e.g., Bayesian) genotype calling to the detection of local adaptation signatures from genomic, transcriptomic, and epigenomic data. We report on progress in addressing central questions of conservation genomics, advances in NGS data analysis, the potential for genomic tools to assess adaptive capacity, and strategies for training the next generation of conservation genomicists.

Genome-wide CpG density and DNA methylation analysis method (MeDIP, RRBS, and WGBS) comparisons

Genome-wide DNA methylation analysis is one of the most common epigenetic processes analysed for genome characterization and differential DNA methylation assessment. Previous genome-wide analysis has suggested an important variable in DNA methylation methods involves CpG density. The current study was designed to investigate the CpG density in a variety of different species genomes and correlate this to various DNA methylation analysis data sets. The majority of all genomes had >90% of the genome in the low density 1-3 CpG/100 bp category, while <10% of the genome was in the higher density >5 CpG/100 bp category. Similar observations with human, rat, bird, and fish genomes were observed. The methylated DNA immunoprecipitation (MeDIP) procedure uses the anti-5-methylcytosine antibody immunoprecipitation followed by next-generation sequencing (MeDIP-Seq). The MeDIP procedure is biased to lower CpG density of <5 CpG/100 bp, which corresponds to >95% of the genome. The reduced representation bisulphite (RRBS) protocol generally identifies DMRs in higher CpG density regions of ≥3 CpG/100 bp which corresponds to approximately 20% of the genome. The whole-genome bisulphite (WGBS) analyses resulted in higher CpG densities, often greater than 10 CpG/100bp. WGBS generally identifies ≥2 CpG/100bp, which corresponds to approximately 50% of the genome. Limitations and potential optimization approaches for each method are discussed. None of the procedures can provide complete genome-wide assessment of the genome, but MeDIP-Seq provides coverage of the highest percentage. Observations demonstrate that CpG density is a critical variable in DNA methylation analysis, and different molecular techniques focus on distinct genomic regions.

Family-effects in the epigenomic response of red blood cells to a challenge test in the European sea bass (Dicentrarchus labrax, L.)

Abstract

Background

In fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from ‘omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments.

Results

We retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs.

Conclusion

Our modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07420-9.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis

Unprecedented rate of increased CO₂ level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis

Unprecedented rate of increased CO₂ level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.

Differential DNA methylation in somatic and sperm cells of hatchery vs wild (natural-origin) steelhead trout populations

Environmental factors such as nutrition, stress, and toxicants can influence epigenetic programming and phenotypes of a wide variety of species from plants to humans. The current study was designed to investigate the impacts of hatchery spawning and rearing on steelhead trout (Oncorhynchus mykiss) vs the wild fish on a molecular level. Additionally, epigenetic differences between feeding practices that allow slow growth (2 years) and fast growth (1 year) hatchery trout were investigated. The sperm and red blood cells (RBC) from adult male slow growth/maturation hatchery steelhead, fast growth/maturation hatchery steelhead, and wild (natural-origin) steelhead were collected for DNA preparation to investigate potential alterations in differential DNA methylation regions (DMRs) and genetic mutations, involving copy number variations (CNVs). The sperm and RBC DNA both had a large number of DMRs when comparing the hatchery vs wild steelhead trout populations. The DMRs were cell type specific with negligible overlap. Slow growth/maturation compared to fast growth/maturation steelhead also had a larger number of DMRs in the RBC samples. A number of the DMRs had associated genes that were correlated to various biological processes and pathologies. Observations demonstrate a major epigenetic programming difference between the hatchery and wild natural-origin fish populations, but negligible genetic differences. Therefore, hatchery conditions and growth/maturation rate can alter the epigenetic developmental programming of the steelhead trout. Interestingly, epigenetic alterations in the sperm allow for potential epigenetic transgenerational inheritance of phenotypic variation to future generations. The impacts of hatchery exposures are not only important to consider on the fish exposed, but also on future generations and evolutionary trajectory of fish in the river populations.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

DNA Methylation Changes in the Sperm of Captive-Reared Fish: A Route to Epigenetic Introgression in Wild Populations

Interbreeding between hatchery-reared and wild fish, through deliberate stocking or escapes from fish farms, can result in rapid phenotypic and gene expression changes in hybrids, but the underlying mechanisms are unknown. We assessed if one generation of captive breeding was sufficient to generate inter- and/or transgenerational epigenetic modifications in Atlantic salmon. We found that the sperm of wild and captive-reared males differed in methylated regions consistent with early epigenetic signatures of domestication. Some of the epigenetic marks that differed between hatchery and wild males affected genes related to transcription, neural development, olfaction, and aggression, and were maintained in the offspring beyond developmental reprogramming. Our findings suggest that rearing in captivity may trigger epigenetic modifications in the sperm of hatchery fish that could explain the rapid phenotypic and genetic changes observed among hybrid fish. Epigenetic introgression via fish sperm represents a previously unappreciated mechanism that could compromise locally adapted fish populations.

Epimutations in Developmental Genes Underlie the Onset of Domestication in Farmed European Sea Bass

Domestication of wild animals induces a set of phenotypic characteristics collectively known as the domestication syndrome. However, how this syndrome emerges is still not clear. Recently, the neural crest cell deficit hypothesis proposed that it is generated by a mildly disrupted neural crest cell developmental program, but clear support is lacking due to the difficulties of distinguishing pure domestication effects from preexisting genetic differences between farmed and wild mammals and birds. Here, we use a farmed fish as model to investigate the role of persistent changes in DNA methylation (epimutations) in the process of domestication. We show that early domesticates of sea bass, with no genetic differences with wild counterparts, contain epimutations in tissues with different embryonic origins. About one fifth of epimutations that persist into adulthood are established by the time of gastrulation and affect genes involved in developmental processes that are expressed in embryonic structures, including the neural crest. Some of these genes are differentially expressed in sea bass with lower jaw malformations, a key feature of domestication syndrome. Interestingly, these epimutations significantly overlap with cytosine-to-thymine polymorphisms after 25 years of selective breeding. Furthermore, epimutated genes coincide with genes under positive selection in other domesticates. We argue that the initial stages of domestication include dynamic alterations in DNA methylation of developmental genes that affect the neural crest. Our results indicate a role for epimutations during the beginning of domestication that could be fixed as genetic variants and suggest a conserved molecular process to explain Darwin’s domestication syndrome across vertebrates.

Domestication may affect the maternal mRNA profile in unfertilized eggs, potentially impacting the embryonic development of Eurasian perch (Perca fluviatilis)

Domestication is an evolutionary process during which we expect populations to progressively adapt to an environment controlled by humans. It is accompanied by genetic and presumably epigenetic changes potentially leading to modifications in the transcriptomic profile in various tissues. Reproduction is a key function often affected by this process in numerous species, regardless of the mechanism. The maternal mRNA in fish eggs is crucial for the proper embryogenesis. Our working hypothesis is that modifications of maternal mRNAs may reflect potential genetic and/or epigenetic modifications occurring during domestication and could have consequences during embryogenesis. Consequently, we investigated the trancriptomic profile of unfertilized eggs from two populations of Eurasian perch. These two populations differed by their domestication histories (F1 vs. F7+-at least seven generations of reproduction in captivity) and were genetically differentiated (FST = 0.1055, p<0.05). A broad follow up of the oogenesis progression failed to show significant differences during oogenesis between populations. However, the F1 population spawned earlier with embryos presenting an overall higher survivorship than those from the F7+ population. The transcriptomic profile of unfertilized eggs showed 358 differentially expressed genes between populations. In conclusion, our data suggests that the domestication process may influence the regulation of the maternal transcripts in fish eggs, which could in turn explain differences of developmental success.

Anadromy, potamodromy and residency in brown trout Salmo trutta: the role of genes and the environment

Brown trout Salmo trutta is endemic to Europe, western Asia and north-western Africa; it is a prominent member of freshwater and coastal marine fish faunas. The species shows two resident (river-resident, lake-resident) and three main facultative migratory life histories (downstream-upstream within a river system, fluvial-adfluvial potamodromous; to and from a lake, lacustrine-adfluvial (inlet) or allacustrine (outlet) potamodromous; to and from the sea, anadromous). River-residency v. migration is a balance between enhanced feeding and thus growth advantages of migration to a particular habitat v. the costs of potentially greater mortality and energy expenditure. Fluvial-adfluvial migration usually has less feeding improvement, but less mortality risk, than lacustrine-adfluvial or allacustrine and anadromous, but the latter vary among catchments as to which is favoured. Indirect evidence suggests that around 50% of the variability in S. trutta migration v. residency, among individuals within a population, is due to genetic variance. This dichotomous decision can best be explained by the threshold-trait model of quantitative genetics. Thus, an individual's physiological condition (e.g., energy status) as regulated by environmental factors, genes and non-genetic parental effects, acts as the cue. The magnitude of this cue relative to a genetically predetermined individual threshold, governs whether it will migrate or sexually mature as a river-resident. This decision threshold occurs early in life and, if the choice is to migrate, a second threshold probably follows determining the age and timing of migration. Migration destination (mainstem river, lake, or sea) also appears to be genetically programmed. Decisions to migrate and ultimate destination result in a number of subsequent consequential changes such as parr-smolt transformation, sexual maturity and return migration. Strong associations with one or a few genes have been found for most aspects of the migratory syndrome and indirect evidence supports genetic involvement in all parts. Thus, migratory and resident life histories potentially evolve as a result of natural and anthropogenic environmental changes, which alter relative survival and reproduction. Knowledge of genetic determinants of the various components of migration in S. trutta lags substantially behind that of Oncorhynchus mykiss and other salmonines. Identification of genetic markers linked to migration components and especially to the migration-residency decision, is a prerequisite for facilitating detailed empirical studies. In order to predict effectively, through modelling, the effects of environmental changes, quantification of the relative fitness of different migratory traits and of their heritabilities, across a range of environmental conditions, is also urgently required in the face of the increasing pace of such changes.

Temporal Dynamics of DNA Methylation Patterns in Response to Rearing Juvenile Steelhead (Oncorhynchus mykiss) in a Hatchery versus Simulated Stream Environment

Genetic selection is often implicated as the underlying cause of heritable phenotypic differences between hatchery and wild populations of steelhead trout (Oncorhynchus mykiss) that also differ in lifetime fitness. Developmental plasticity, which can also affect fitness, may be mediated by epigenetic mechanisms such as DNA methylation. Our previous study identified significant differences in DNA methylation between adult hatchery- and natural-origin steelhead from the same population that could not be distinguished by DNA sequence variation. In the current study, we tested whether hatchery-rearing conditions can influence patterns of DNA methylation in steelhead with known genetic backgrounds, and assessed the stability of these changes over time. Eyed-embryos from 22 families of Methow River steelhead were split across traditional hatchery tanks or a simulated stream-rearing environment for 8 months, followed by a second year in a common hatchery tank environment. Family assignments were made using a genetic parentage analysis to account for relatedness among individuals. DNA methylation patterns were examined in the liver, a relatively homogeneous organ that regulates metabolic processes and somatic growth, of juveniles at two time points: after eight months of rearing in either a tank or stream environment and after a subsequent year of rearing in a common tank environment. Further, we analyzed DNA methylation in the sperm of mature 2-year-old males from the earlier described treatments to assess the potential of environmentally-induced changes to be passed to offspring. Hepatic DNA methylation changes in response to hatchery versus stream-rearing in yearling fish were substantial, but few persisted after a second year in the tank environment. However, the early rearing environment appeared to affect how fish responded to developmental and environmental signals during the second year since novel DNA methylation differences were identified in the livers of hatchery versus stream-reared fish after a year of common tank rearing. Furthermore, we found profound differences in DNA methylation due to age, irrespective of rearing treatment. This could be due to smoltification associated changes in liver physiology after the second year of rearing. Although few rearing-treatment effects were observed in the sperm methylome, strong family effects were observed. These data suggest limited potential for intergenerational changes, but highlight the importance of understanding the effects of kinship among studied individuals in order to properly analyze and interpret DNA methylation data in natural populations. Our work is the first to study family effects and temporal dynamics of DNA methylation patterns in response to hatchery-rearing.