- Research article
- Open Access
Family-effects in the epigenomic response of red blood cells to a challenge test in the European sea bass (Dicentrarchus labrax, L.)
BMC Genomics volume 22, Article number: 111 (2021)
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.
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.
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.
Because samples are easy to obtain, poorly invasive, and can be stored in large collections that may reflect variation in many parameters at both the individual and the population levels, blood is certainly the most commonly used tissue to check for and to monitor the response of cells, organs, or whole organism to environmental perturbations, to assess health status of organisms, and to diagnose metabolic impairments and dysfunctions in vertebrates. As a tissue subjected to systematic hormonal fluctuations by a centrally produced stress response, blood is especially used to monitor stress indicators at the molecular, cellular or physiological levels in teleost [1, 2]. Plasma cortisol (the main glucocorticoid hormone) as a primary physiological stress indicator and few metabolites such as glucose and lactate as secondary physiological indicators are certainly the most commonly assessed biomarkers of stress in fish . These plasma biomarkers combine interesting advantages for stress monitoring (e.g., cheap data generation, nonlethal). Nevertheless, because the response of fish to stressors requires the consideration of a complex regulatory network of non-linear actions that could not be fully integrated by few parameters, it has been proposed that new technologies should give rise to new biomarkers for fish biomonitoring, especially to improve welfare in the farmed environment . Indeed, the last decade has seen the emergence of a number of technologies for quantifying the molecular responses of fish to stressors at a genome-wide scale, including transcriptomics, proteomics, and epigenomics (e.g. [4,5,6,7,8,9,10,11]). Omics studies traditionally target key organs for stress monitoring such as the brain, the kidney, or the liver, but tissue sampling is generally lethal.
Because fish blood cells are nucleated and, apart from blood plasma in which cortisol, glucose, lactate and other metabolites are measured, also mobilized as part of the stress response in fish [12, 13], it is appealing to investigate if components of their genomic machinery may respond to environmental stressors and broaden the panel for poorly invasive stress monitoring. To data, the use of red blood cells (RBCs) in ‘omics fish studies has received little attention [14,15,16], and a single study specifically investigated RBC epigenome in steelhead (Onchorhynchus mykiss) .
After salmonids, the European sea bass (Dicentrarchus labrax) is certainly the most investigated marine fish species in Europe using molecular tools. It has been extensively studied over the last three decades, in both natural and farmed populations (reviewed in ). This includes the sequencing of its genome  and an increasing number of epigenetic studies [20,21,22,23,24,25,26,27]. In this economically important fish (approx. 200,000 t produced worldwide in 2018 ), epigenetic studies covered research areas important to fish farming including, e.g., sex determination [19, 24], the dynamics of epigenetic marks in sperm , the effects of temperature , or the epigenetic impacts of the onset of domestication . However, only one of these studies was carried out at the genome-wide scale , others focusing at modifications of epigenetic profiles for reduced gene sets. None of these studies explicitly targeted ‘stress’ (but see ), and stress monitoring in the European sea bass remains largely evaluated using blood plasma (or serum) parameters (e.g. [29,30,31,32]). Some authors proposed alternatives based on, e.g., gene expression, but, by traditionally targeting tissues such as liver, brain or kidneys, they are invasive and fish are sacrificed in most of the cases (e.g. ). How the RBC methylome analyzed in minimally invasive blood samples may capture components of the stress response is actually missing in sea bass.
In this study, we adapted the epiGenotyping By Sequencing (epiGBS) protocol originally proposed by Van Gurp et al.  to assess the genome-wide epigenomic variation in the RBCs of D. labrax submitted to periods of acute stress during a 3 month challenge test. EpiGBS targets variation in cytosine methylation – the covalent addition of a methyl group to cytosine nucleotides – that has long been accepted as an important epigenetic modification in many organisms [35, 36]. This modification integrates a second restriction enzyme and further multiplexing of individuals. Our aim was to explore the changes in the epigenomic landscape of sea bass RBCs in pre- and post-challenge fish to initiate and to motivate the use of differentially methylated cytosines (DMCs) as putative biomarkers of stress.
Twenty sea bass families were produced to initiate a 3 month test in 6 month-old individual sea bass. This challenge was seeded with 20 individuals of each family (N = 400), minimizing tank effects. During the full challenge, fish were regularly submitted to acute stress, then could recover (see Methods section for details). In order to evaluate if this challenge could induce genome-wide methylation changes in sea bass RBC, a total of seventy-four randomly caught individuals (37 pre- [T0] and 37 post-challenge [T4] out of the 400 fish) were considered in this study. All individuals were submitted to the challenge, no unstressed individuals were available (see Methods section). While developed on a family-based experimental design, we only compared methylation difference between pre- and post-challenge juvenile sea bass and did not compare families in this study. Indeed, random sampling induced uneven representation of families within and among samples, and only nineteen out of 20 sea bass families were represented by at least one individual among the 74 samples analyzed in this study. Fish number per family ranged from one (families A, D, N) to nine (family R) individuals. Except for the families with a single representative and family M with post-challenge fish only (four), both pre- and post-challenge individuals were present in the 15 remaining families. Also because of random sampling, four individuals from four distinct families were retained twice by chance (Fig. 1). They were thus analyzed for both pre- and post-challenge conditions. A total of 70 distinct fish has been analyzed in this study.
EpiGBS library construction and sequencing
We obtained 504,271,331 total sequencing reads of which 99.4% (501,108,033) were retrieved after trimming of our single library. After demultiplexing, read numbers per sample ranged from 2,284,915 to 16,314,759, with an average of 5,212,596 reads per sample (see Additional File 1). Demultiplexed samples were mapped against the D. labrax reference genome (~ 676 Mb) with a mean mapping efficiency of 74.5% (73.0% for unique best hits; Additional File 1). Sequencing reads mapped across all linkage groups (Additional File 2). The mean per base pair read depth was 250X.
Out of the 10,368,945 CG dinucleotides present in the MspI-SbfI reduced-representation of D. labrax genome we obtained, 47,983 CpG coordinates were extracted with a minimum of 30X read depth and presence in at least 20 individuals. They were filtered out using a 15% methylation difference threshold and a nominal cut-off value of q < 0.001. With these parameters, only a total of 57 cytosines in CpG context were defined as DMCs between pre- and post-challenge sea bass (Table 1). Methylation differences ranged up to 46.4% for hypermethylated cytosines, and down to − 27.5% for hypomethylated cytosines. Hyper-methylation was more frequently detected than hypomethylation (11 [19.30%] hypo- vs 46 [80.70%] hypermethylated DMCs) in post-challenge sea bass. DMCs were distributed on 17 out of 24 LG groups and in or close to 51 distinct genes. Further information is provided in Additional File 3 (e.g. gene annotations, CpG context).
Most identified DMCs were located within identified gene bodies (44 out of 57, 77.19%), one in the 3’UTR regions of the Solute Carrier family 22 Member 2 (SLC22A2) gene on LG17, and one in a repeated region (a non LTR Retrotransposon Element on LG20). In the remaining cases (n = 11), DMCs are intergenic and located in a window ranging from 0.9 kb to 51 kb to the closest gene (respectively: SASH1A on LG12 and TRMT11 on LG16; Table 1). Two pairs of overlapping, but inversely oriented genes share on their sense vs antisense strand an identical DMC: PLG and SLC22A2 on LG17, and SART3 and FICD on LG20 (Table 1).
When located on the same LG, DMCs were usually distant by at least 30 kb from each other. In only three instances, some DMCs were located close from each other (< 1.5 kb). These DMCs target the same gene (Table 1). This includes three hypo-methylated cytosines located ~ 1500 bp downstream of the predicted Kelch Repeat and BDB domain 13 (KBTBD13) gene with at most 88 bp between the cytosines. It also includes three hypomethylated cytosines (> 20%) on LG1A in the second intron of the forkhead box J3 (FOXJ3) gene. One other groups of two cytosines were found in the same exon (distant by 3 bp) of the BTR30 gene (Table 1). A single DMC was associated to a repeat region and two DMCs were found to refer to the same gene (homologous to the Gasterosteus aculeatus paralogue of COL4A5, a collagen gene of type IV mostly implicated in the protein network of the basement membrane) (Table 1). For this gene, one DMC is located in the first intron while the second is 16.5 kb upstream of the start codon.
Hierarchical clustering showed a strong family effect in methylation patterns (i.e. individuals within family clustered together; Fig. 1). The four individuals that were caught twice clustered together by pairs in all four cases. These individuals have the lowest levels of dissimilarity in hierarchical clustering, suggesting that family – not fully considered in our sampling scheme - may explain considerably more variation than treatment in their methylation profiles. Despite this strong family effect and clues of low impact of the challenge test on methylation, pre- and post-challenge groups can be distinguished based on their DMC profile in PCA. Mean loading scores of individuals were found significant among T0 and T4 for PC1 that explained 7.2% of total variance (Student t-test; p < 0.005, Fig. 2). No significant difference was found for loading scores along PC2 (3.0% of total variation; P = 0.404).
Database mining screening for specific protein interactions on the String server revealed few possible pairs of associations between DMC-related genes (n = 8). These associations involve ROBO3-CHN1, ROBO3-PRKCQ, ROBO3-LRRC3, DLG1-NCBP2, FURIN-PLG, PLG-MMRN2a, CELF-RRM2, and CRTC2-DENND4B (see Additional File 4), with some them possibly linked to stress. FURIN - a subtilisin-like protein proconvertase – and plasminogen (PLG) processed proBDNF to mature BDNF (Brain-Derived Neurotrophic Factor), one of the most important molecule in fear memory (see Discussion). ROBO3 and CHN1 have been shown to interact with poorly-understood implications of CHN1 in stress disorders . CELF2 (CUGBP Elav-like family member 2) and RRM2 (Ribonucleotide reductase M2 polypeptide) are both known to participate to messenger RNA (mRNA) metabolism . CELF2 acts to post-transcriptionally stabilize mRNAs by relocating them to stress granules in the cytosol. CELF2 interferes with RRM2 that modulates its splicing activity. As post-transcriptional activities are at the core of methylation studies, the detection of this association seems relevant to our study.
We showed that a modified epiGBS protocol originally proposed by Van Gurp et al.  was applicable to further analyze patterns of cytosine methylation in RBCs of D. labrax. This is the first use of epiGBS in fish and the second in an animal species (Canadian lynx ). Overall, RBC’s DNA methylation was shown to respond to the challenge test, but observed changes were found mainly explained by the genetic background of individuals resulting from family-based effects, and involved relatively few sites and DMC-related genes.
Mining the sea bass epigenome
The addition of a second restriction enzyme illustrates the flexibility of the epiGBS originally proposed by Van Gurp et al.  and more generally of reduced-representation bisulfite sequencing (RRBS) protocols for data acquisition and impact. The addition of a second restriction enzyme to a RRBS protocol in order to improve coverage and accuracy of CpG methylation profiling was however already shown , but hereby proposed in a context of improved multiplexing of samples.
The information provided in this study is based on the analysis of 47,983 distinct methylated sites distributed over all sea bass LGs. The mapping efficiency was high (74.5%) when compared to early values retrieved in human (~ 65%) , or in fish studies screening for genome-wide methylation (e.g. 55–60% in ; 40% in ). Other studies reported similar mapping efficiencies, but reported percentages of mapping for unique best hits that were generally lower. For example, in Kryptolebias marmoratus, Berbel-Filho et al.  reported a mean mapping efficiency of 74.2% but 61.1% unique best hits while, in this study, this latter percentage reached 73.0%. This reflects a more robust mapping of the DMCs we detected and significantly enlarge the breadth of the sites that can confidently exploit to retrieve functional information. Taking advantage of the epiGBS protocol that allow to process more samples , the number of individuals considered in this study is rather high (n = 70 distinct individuals), when most epigenomic studies in fish dealt with less than 30 individuals (range: n = 3 in ; n = 106 in  for a population study). In sea bass, Anastasiadi and Piferrer  previously reported a study that used 27 samples and as many libraries to be sequenced while our data were obtained from a unique library preparation. Our modified epiGBS protocol provides a considerable amount of information, certainly at a reasonable cost, to decipher methylation landscapes of sea bass or other species.
The operational and statistical thresholds used in the successive steps of this study are conservative, resulting in the discovery of a rather low total number of methylated sites, but certainly limiting the report of false positives. For example, a threshold of 30X and nominal cut-off value of 0.001 are quite conservative, when some studies might consider a threshold of 5X or 10X for a CpG to be analyzed and associated cut-off values of 0.05 or 0.01 (e.g. [26, 46, 50]). Relaxing thresholds would enable to retrieve more DMCs, but elevated thresholds should normally ensure that access to relevant information is reached. Thus, only 57 DMCs have been found in RBCs of pre- and post-challenge European sea bass. These DMCs were found mostly hypermethylated in post- compared to the pre-challenge individuals, and mostly located in gene bodies (i.e. the transcriptionally active portion of the genome) of fifty-one different genes. Differential methylation in gene bodies may regulate splicing and/or act as alternative promoters to reshape gene expression [51,52,53].
In addition to DMCs located in gene bodies, a dozen of DMCs were found in intergenic regions (21.0%). Intergenic cytosine methylation has been frequently described, including in response to stress , but its role remains poorly understood . While numbers of genic vs intergenic DMCs may greatly vary, a ratio of ~ 80% of DMCs located in gene bodies and ~ 20% located in other genomic regions has been reported in other fish studies (e.g. ).
An epigenomic perspective on stress biomarkers
Studies looking at the epigenomic landscape of RBCs in fish are scarce, and did not focus on response to a so-called stress challenge . Our study yielded mixed results regarding this issue. Negatively, this study did not considered controls (i.e. unstressed fish) and it is difficult to assess if cytosines that were shown to respond to the challenge test really reflect the impact of stress or other parameters. This notably includes growth and ontogenetic changes in first year juvenile sea bass, together with sexual differentiation. Sexual differentiation occurs between 150 and 250 days post-fertilization in sea bass (e.g. ) and differential methylation measured in gonads at few candidate genes has been reported over this period [20, 23]. As our challenge test covers this period, results might be partially influenced by sexual differentiation. This has to be investigated further. However, we are not aware of studies that showed that differential methylation recorded in gonads might translate to RBCs, and none of the candidate differentially methylated genes previously studied in sea bass gonads was detected in this study. Methylation variation accompanying ontogeny and/or aging is reported in fish [56, 57], and it has been shown to be modified with age in sea bass muscles . Unfortunately, methylation results also concerned candidate genes not detected as differentially methylated in this study. Relevant to our study, BMP3 , FURIN , NOL4B , Myf5 [61, 62], NPAS3 , and ROBO3  are engaged in the development of the anterior region and/or the craniofacial skeleton which is known to be modified during sea bass farming . Their roles were however studied in early development stages of mammals or zebrafish. As ‘epigenetic programming’ – apart of transgenerational inheritance - is mostly an early-life process that influence late-life effects (in fish, see, e.g., [5, 9, 38, 65]), we thus hypothesize that the epigenetic marks that could affect them would have been already present in 6 month-old sea bass that initiated the challenge test. Developmentally induced differential methylation acquired during the challenge test seems unlikely. Furthermore, some of them have clear relationships to stress exposure (e.g. FURIN, NPAS3; see below). However, in absence of dedicated study, we cannot totally rule out that methylation patterns observed in sea bass RBCs could also partly reflect the developmental or sexual regulation of a particular phenotype between pre- and post-challenge fish, rather than being directly related to the challenge.
Nevertheless - and more positively - some DMC-related genes detected in this study have been shown to be involved in the stress response in fish. Strikingly, 37 over 51 DMC related-genes were reported mainly from few neurotranscriptomics zebrafish studies that dealt either with reactive-proactive behavioural response to stress [37, 39] or with changes in social regulation that may promote stressful behaviour among congeners [40, 66] (Table 1). While not detected in RBCs but in brain tissues, correspondence across stress studies is interesting in this particular case. Indeed, while not investigated in fish so far, human stress studies have shown that blood cells responded to DNA methylation in the brain [67,68,69] (but see ), and specifically RBCs in birds . Furthermore, several of these DMC-related genes are involved in the maturation of proBDNF to mature BDNF or the regulation of its activities (ABLIM2, ADCY1b, CRTC2, FURIN, NPAS3, PLG, and possibly SLC22A2 and DLG1). BDNF consolidates both the within- and between-generation fear memory owing to epigenetic regulation [72, 73]. Its activity is strongly linked to glucocorticoid stress to imprint neurogenesis  and it acts as both a regulator and a target of stress hormone signaling . BDNF is one of the target genes in fish stress studies (zebrafish , sea bream , sea bass [33, 78,79,80]), but its methylation status could not be investigated in this study as no SbfI restriction site is present within or close to this gene. However, the above-mentioned DMC-related genes might be linked to a ‘BDNF network’. The adenyl (ate) cyclase (AC, ADCY1b gene) is a brain-specific signaling enzyme that synthesizes the cyclic AMP . This inducible signaling pathway participates to the synthesis of the active form of BDNF (proBDNF to mature BDNF) . ProBDNF is processed by furin and the plasminogen system , including processing steps that necessitate actions of actin-binding LIM kinases (ABLIM) . Stress imprinting at FURIN is likely and it has recently been shown that transgenerational epigenetic effects of furin activity were active in brain of mice . Furin has also been shown to modulate learning abilities and memory . ProBDNF cleavage by furin depends on brain AC and CREB (cAMP response element-binding protein) signaling [87, 88], but also plasminogen (PLG) . This activity is modulated by stress hormones (corticosteroids) [90, 91]. One interesting supplementary observation is that CREB signaling necessary to furin is associated to CRTC2 - a CREB co-activator. In mice, CRTC2 is known to act as a switch for BDNF and glucocorticoids to direct the expression of corticotropin-releasing hormone (CRH) in the hypothalamus . Additionally, in the brain, plasminogen encoded by PLG is converted to plasmin that cleaves BDNF in the extracellular synaptic domain [83, 93]. PLG has also been shown to regulate pro-opiomelanocortin (POMC) in the hypothalamic-pituitary axis, then the production of peptides hormones such as the adrenocorticotropic hormone (ACTH) .
Few other DMC-related genes should be mentioned. SLC22A2 - also known as OCT2 (organic cation transporter 2) – associated to the unique 3’UTR DMC found in this study is involved in numerous transmembrane transports , including at the blood-brain barrier . It was found involved in memory in mice [97, 98] or Drosophila . In this study, PLG and SLC22A2 are associated to the same DMC; the functional significance of this situation needs further investigation. NPAS3 (neuronal PAS domain containing protein 3) has a well-established action in memory [100, 101], and participate to a neural network that also includes BDNF . NPAS3 is also associated to the glial cell line-derived neurotrophic factor receptor-alpha2 gene (GFRA2) detected in this study and related to stress and anxiety . Finally, DLG1 (Disk-large homolog 1) plays a critical role in neural synapse formation, insulin secretion and glucose transport that are activated or modulated by stressors . Adrenergic modulation implying DLG1 was also found to correlate with emotional states and stress sensitivity in mice  and it indirectly participates to the regulation of BDNF as DLG1 activates the glutamate receptor 1 (GluR1 ;) that interacts with molecular processing of BDNF . A relationship of DLG1 with NCBP2 (nuclear cap binding protein 2) was detected in protein-protein interaction analyses. NCBP2 protect cellular RNA polymerase II transcripts from degradation and to guide them through the sequence of steps leading from transcription to translation . As the cellular response to environmental challenges requires immediate and precise regulation of transcriptional programs, differences in cytosine methylation among pre- and post-challenge sea bass close the NCBP2 gene could partly reflect the impact of the challenge.
While confounding factors may be present, results thus suggest that some DMCs reported in this study did not occur only by chance, are related to processes that regulate the hypothalamus-pituitary-interrenal (HPI) axis and hormones, but also to features that are expected to be developed or regulated during a challenge test (e.g. anxiety, fear memory, neurogenesis). While presumptive, these DMCs could effectively reflect the impact of the challenge test, and suggest that traditional blood plasma biomarkers could be potentially enriched by epigenetic marks to monitor welfare in cultured fish species. The link between brain and blood epigenomics remains however to be explored more deeply in fish and requires careful evaluation and validation to correct for tissue specificity, as requested in human . Recent results on chicken RBCs are encouraging .
Hereby, we focused on possible relationships among DMC-related genes detected in this study and expressed in stress-related studies in brain tissues of fish [37, 39, 40, 66]. It should be however important to note that some DMC-related genes could be related to other components of the stress response (e.g., immune response, glucose metabolism). For example, a role of PRKCQ (PKC-theta, a protein kinase C theta type) in the immune response is well-known in vertebrates (e.g. ). PRKCQ is also known to participate to glucose metabolism, including glucose homeostasis . One association with DLG1 is reported, also related to the immune response . The role of CRTC2 on glucose homeostasis when facing stress has also been repeatedly reported [112,113,114], notably in relation to glucocorticoid levels . We cannot expand further on this topic, but this suggests that DMC-related genes detected in this study may integrate several aspects of the stress response in fish.
A family effect, but the possible absence of individual response
As in other fish species [116,117,118], a family effect imprinting the methylome was detected in this study. In parallel, results showed that the epigenetic profiles of the four individuals that were analyzed in the pre- and post-challenge conditions clustered very closely from each other. While based on few observations of randomly sampled individuals, the challenge had little impact on sea bass cytosine methylation landscape in comparison to family effects. Nature and strength of family-based epigenomic variation are of considerable importance attention to engage future selection breeding improvements in cultured fish like sea bass, including issues about health and welfare . More generally, how transgenerational and within-generation stress-imprinting events may interact to shape both the plastic and the heritable component of the stress response in relation to environmental stimuli require in depth evaluation [120, 121]. To do so, far more complex and rigorous experimental designs that the one followed in the present study as to be adopted and temporal monitoring of the individual response of blood methylome to stress has to be promoted. In sea bass, such a research has been engaged for sex determination . Results showed that some epigenetic marks were more likely engaged in transgenerational inheritance, while others be related to within-generation differences acquired during early development .
Conclusions to this study are twofold. First, the European sea bass has become one of the most studied species in fish epigenetics [20,21,22,23,24,25,26,27], and for this species or for other cultured fish species, our modified version of the original epiGBS protocol seems to be a powerful and affordable method to screen a significant number of individuals with sufficient depth and coverage to reach meaningful conclusions. While this protocol should be compared to others (e.g. [44, 122, 123]), its use certainly deserves attention to design more integrated epigenomic-genomic studies [124, 125], as multi-omics investigations of stress, health and welfare [126, 127]. Second, as , we showed that RBCs are amenable to epigenomic investigations at a genome-wide scale in fish. RBC methylome revealed a family-based response to a challenge test. Efforts should be dedicated to validate some DMCs as stress biomarkers using improved experimental designs that would have, e.g., to set up baselines and to estimate context- and/or species-specific differences in order to enlarge the panel of diagnostic tools to monitor good practices in production setups .
Rearing, stress challenge, and blood sampling
Four hundred European sea bass were initially used in this study. Fish were produced in the hatchery facility of Nireus S.A. from breeders maintained in this company for scientific purposes. They resulted of a single crossing experiment (12 dams, 20 sires) that took place in January 2018. During the larval phase, 20 different families labelled from A to T were raised. Each family was reared separately in open circulation tanks at Nireus S.A. research facilities (Greece). To avoid environmental effects, fish were tagged at ~ 180 days post-hatch (July 10-13th 2018) and distributed in 20 tanks; each tank receiving 1 fish from each family (i.e., 20 fish per tank). Fish were fed twice a day, for 6 days a week, using a commercial diet (Blue Line 45:20 3.5 mm, Feedus S.A., Greece). Throughout the experimental period, the photoperiod was set at 12 L:12D, the water temperature and the salinity held constant (18.1 ± 0.2 °C and 28 ppt, respectively). Fish weights (mean ± SD) were 48.1 ± 12.8 g and 86.2 ± 24.1 g in pre- and post-challenge individuals, respectively.
Fish were submitted to one acute challenge test per month for three consecutive months, from July to October 2018. During this challenge, fish were exposed to high density stress by lowering water levels in the tank to 1/3 of the original volume, followed by chasing of the fish with a net for 5 min and a 30 min waiting period before sampling. These stressors are classical in sea bass studies regarding response to acute stress . This protocol took place in each tank then for each fish individual entering the experiment. It was repeated for 3 consecutive times at 20–21 day intervals (period long enough for fish to recover).
An initial blood sampling occurred 2 weeks prior to the implementation of the stress challenge (hereafter T0, July 19th, 2018, pre-stress/control group) then at the end of the challenge test (hereafter T4, October 5th, 2018). Blood samplings were performed in anesthetized fish. Specifically, fish were anesthetized in 2-phenoxyethanol (300 ppm; Merck; 807,291; USA). Recovery was performed in a separate tank with provision of air before fish returned to their holding tank. Blood from the caudal vessel was collected. Plasma and RBCs were separated by centrifugation at 2000 g for 10 min. Careful separation of the plasma and RBCs was performed using 200 μl pipettes. RBC extracts were heparinized (heparin sodium; Sigma-Aldrich), transferred in microtubes and conserved at − 20 °C in 1 ml of RNA later. DNA was extracted using the Macherey Nagel Nucleo Spin Tissue DNA kit and quantified using a Qubit fluorometer (Qubit dsDNA BR Assay Kit, Q32853, Invitrogen). Thirty-seven blood samples at T0 (pre-stress) and thirty-seven additional samples at T4 (post-challenge) have been randomly selected for the downstream epigenomic analysis.
Library preparation and sequencing
We followed the epiGBS protocol published by van Gurp et al. . As the method was developed for plants (with methylation occurring in CpG, CHH, and CHG context; H being any nucleotide but a cytosine) and used a single digestion approach, the protocol was modified to make it more suitable and straightforward for our vertebrate system. Particularly, a double digest instead of the single digest approach was implemented. We chose the restriction enzyme MspI, a standard choice in reduced representation bisulfite sequencing-like (RRBS) studies, as its recognition site targets CpG rich regions . A second enzyme (SfbI) with a recognition site length of 8 bp was used to reduce fragment numbers, and thus to increase read depth per fragment. The choice of this enzyme was guided by in silico digestion of the European sea bass genome  using simRAD . This genome is available at: https://www.ensembl.org/Dicentrarchus_labrax/Info/Index.
One single library was prepared for a set of 74 samples. For this library, 200 ng of DNA of each sample were digested in a 40 μl reaction, using 0.25 μl MspI (NEB, 20,000 U/ml R0106S), 0.25 μl SbfI-HiFi (NEB 20,000 U/ml R3642L) and 4 μl of 10X cutsmart buffer. The reaction was run overnight at 37 °C. Unique forward and reverse adapter combinations allow multiplexing samples in the library. We added forward and reverse adapters in unique combinations (1 μl of adapter, 2.5 μM), 0.5 μl T4 Ligase (NEB, 400,000 U/ml m0202L), 6 μl of 10X T4 ligase buffer and 11.5 μl water were added directly to the digested DNA. Sequences of adapters are provided as Additional File 1. Adapters were ligated for 3 h at 23 °C followed by 10 min of enzyme inactivation at 65 °C. After ligation, all samples were pooled and one third of the total volume was used in the following step. The mixture volume was reduced using a Qiaquick PCR purification kit (28,104 Qiagen). The resulting product was cleaned a second time to ensure the removal of small fragments and adapter remnants using CleanPCR paramagnetic beads (Proteigene, CPCR-0050) with a ratio of 0.8X [sample:beads].
Since oligonucleotides were not phosphorylated (see Additional File 5), a nick translation was performed to repair the nick of the DNA at the restriction site and to fully methylate the hemi-methylated adapters. To do so, 19.25 μl of the concentrated and purified ligation pool were used in a 25 μl reaction, including 0.75 μl DNA Polymerase I (E. coli, 10,000 U/ml M0209S), 2.5 μl of 10 mM 5-Methylcytosine dNTP mix (Zymo Research, D1030), and incubated at 15 °C for 1 h. The library was then treated with sodium bisulfite using the EZ DNA Methylation Gold Kit (D5005, Zymo Research) following manufacturer instructions to convert unmethylated cytosines to uracil, paying attention to the optimal DNA amount per reaction. After bisulfite conversion, 14 cycles of PCR (95 °C for 1 min, 91 °C for 10s, 65 °C for 15 s, 72 °C for 10s and a final elongation step at 72 °C for 5 min) were performed, followed by a paramagnetic bead clean up with a 0.8X ratio. The library quality (fragment size distribution, no adapters left, no primers left; fragment size range from 300 to 800 bp; see Additional File 6) was verified on an Agilent 5300 Fragment Analyzer (Santa Clara, USA). It was sequenced (paired-end, 150 bp) on one lane of a SP flow cell on an Illumina™ NovaSeq 6000 at the MGX sequencing facility in Montpellier, France.
Raw sequencing data (NCBI accession number: GSE153838) were trimmed of low quality reads and adapter residues using Trim Galore! (v.0.6.4; available at https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). The trimmed reads were then demultiplexed using the process_radtags command of Stacks . The option disable_rad_check was applied to avoid reads with bisulfite-modified restriction sites to be discarded. The demultiplexed reads were trimmed a second time on both 5′ and 3′ ends with the options --clip_R1 --clip_R2 and --three_prime_clip_R1 --three_prime_clip_R2 to remove introduced methylated cytosines during adapter ligation. Using the bismark_genome_preparation function of the program bismark , a bisulfite converted version of the European sea bass genome was prepared against which the demultiplexed and trimmed reads were mapped.
The R package MethylKit  was used to determine differential methylation between pre- and post-challenge fish. CpGs with less than 30 read depth and with coverage > 99.9% of the distribution of read counts were filtered out to account for PCR bias. Coverage was normalized across samples, and only CpGs present in at least 20 out of 37 samples per group (56%) were kept for further analysis. Differentially methylated CpG sites were determined between pre- and post-challenge individuals using logistic regression (calculateDiffMeth function). Pre-stress individuals were considered as the baseline. Cytosines were considered as differentially methylated when presenting at least 15% methylation difference (as in  for epigenomic variation in D. labrax), and a nominal q-value < 0.001 between pre- and post-challenge fish. The combination of thresholds, sample size, and previously mentioned depth coverage ensures detection of biologically meaningful differences. Reads presenting DMCs were extracted using Geneious (v.11.0; available at https://www.geneious.com/) and mapped along LGs of the European sea bass reference genome to produce a primary annotation of potential candidate genes. As the annotation of the European sea bass and the sea bream (Sparus aurata) reference genomes have been recently released on Ensembl (http://www.ensembl.org/index.html), gene names and models have been controlled, eventually including a TBLASTN search in case of discrepancy.
A principal component analysis (PCA) on the methylation profiles of individual samples at DMCs was performed to analyse the potential grouping structures within our data set. Differences in the distribution of mean individual loading scores for pre- and post-challenge fish were tested along each principal components axis using a Student t-test. Additionally, a hierarchical clustering was performed using Ward’s linkage method on Euclidean distances in order to explore whether sea bass families structure the data set.
We used String (https://string-db.org/) to investigate if DMC-related genes encoded for proteins are known to interact together. We used our gene list as input and the zebrafish genome as a reference for annotation. When one interaction was provided we specifically explored the literature for confirmation and relevant experimental evidence regarding stress.
Availability of data and materials
Individual Illumina raw reads (one round of trimming) and processed files have been deposited in NCBI’s Gene Expression Omnibus . Data set used in this study is accessible through Gene Expression Omnibus (GEO) Series accession number GSE153838 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE153838).
Actin-binding LIM kinase
Brain Derived Neurotrophic Factor
Bone morphogenetic protein 3
Bloodthirsty-related gene family member 30
Cyclic adenosine monophosphate
CUGBP Elav-like Family Member 2
Collagen type IV alpha 5 chain
CREB regulated transcription coactivator 2
cAMP response element-binding
Denn domain-containing protein 4b
Disk-large homolog 1
Differentially methylated cytosine
epiGenotyping By sequencing
FIC domain protein adenylyltransferase
Forkhead box J3
Kelch Repeat and BDB domain 13
Leucine rich repeat containing 3
Myogenic factor 5
National Center for Biotechnology Information
Nuclear Cap Binding Protein 2
Nucleolar protein 4B
Neuronal PAS domain protein 3
Principal component analysis
Protein kinase C theta
Red blood cell
Roundabout guidance receptor 3
Reduced representation bisulfite sequencing
Ribonucleotide reductase regulatory Subunit M2
Spliceosome associated factor 3
SAM and SH3 domain-containing 1a
Solute carrier family 22 member 2
TRNA methyltransferase 11 Homolog
Sopinka NM, Donaldson MR, O’Connor CM, Suski CD, Cooke SJ. Stress indicators in fish. In: Schreck CB, Tort L, Farrell AP, Brauner CJ, editors. Biology of Stress in Fish. Fish Physiology, Vol. 35. Cambridge: Academic Press; 2016. p. 405–62. https://doi.org/10.1016/B978-0-12-802728-8.00011-4.
Faught E, Aluru N, Vijayan MM. The molecular stress response. In: Schreck CB, Tort L, Farrell AP, Brauner CJ, editors. Biology of Stress in Fish. Fish Physiology, Vol. 35. Cambridge: Academic Press; 2016. p. 113–66. https://doi.org/10.1016/B978-0-12-802728-8.00004-7.
Prunet P, Overli Ø, Douxfils J, Bernardini G, Kestemont P, Baron D. Fish welfare and genomics. Fish Physiol Biochem. 2012;38:43–60. https://doi.org/10.1007/s10695-011-9522-z.
Aluru N, Vijayan MM. Stress transcriptomics in fish: a role for genomic cortisol signaling. Gen Comp Endocrinol. 2009;164(2–3):142–50. https://doi.org/10.1016/j.ygcen.2009.03.020.
Robinson NA, Johnsen H, Moghadam H, Andersen Ø, Tveiten H. Early developmental stress affects subsequent gene expression response to an acute stress in Atlantic salmon: an approach for creating robust fish for aquaculture? G3: genes, genomes. Genetics. 2019;9(5):1597–611. https://doi.org/10.1534/g3.119.400152.
Huang V, Butler AA, Lubin FD. Telencephalon transcriptome analysis of chronically stressed adult zebrafish. Sci Rep. 2019;9(1):1379. https://doi.org/10.1038/s41598-018-37761-7.
Raposo de Magalhães C, Schrama D, Farinha AP, Revets D, Kuehn A, Planchon S, et al. Protein changes as robust signatures of fish chronic stress: a proteomics approach to fish welfare research. BMC Genomics. 2020;21(1):309. https://doi.org/10.1186/s12864-020-6728-4.
Groh KJ, Suter MJ. Stressor-induced proteome alterations in zebrafish: a meta-analysis of response patterns. Aquat Toxicol. 2015;159:1–12. https://doi.org/10.1016/j.aquatox.2014.11.013.
Moghadam HK, Johnsen H, Robinson N, Andersen O, Jorgensen EH, Johnsen HK, et al. Impacts of early life stress on the methylome and transcriptome of Atlantic salmon. Sci Rep. 2017;7:5023. https://doi.org/10.1038/s41598-017-05222-2.
Todd EV, Ortega-Recalde O, Liu H, Lamm MS, Rutherford KM, Cross H, et al. Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change. Sci Adv. 2019;5(7):eaaw7006. https://doi.org/10.1126/sciadv.aaw7006.
Uren Webster TM, Rodriguez-Barreto D, Martin SAM, Van Oosterhout C, Orozco-terWengel P, Cable J, et al. Contrasting effects of acute and chronic stress on the transcriptome, epigenome, and immune response of Atlantic salmon. Epigenetics. 2018;13(12):1191–207. https://doi.org/10.1080/15592294.2018.1554520.
Tort L. Stress and immune modulation in fish. Dev Comp Immunol. 2011;35:1366–75. https://doi.org/10.1016/j.dci.2011.07.002.
Burgos-Aceves MA, Lionetti L, Faggio C. Multidisciplinary haematology as prognostic device in environmental and xenobiotic stress-induced response in fish. Sci Tot Environm. 2019;670:1170–83. https://doi.org/10.1016/j.scitotenv.2019.03.275.
Morera D, Roher N, Ribas L, Balasch JC, Donate C, Callol A, et al. RNA-Seq reveals an integrated immune response in nucleated erythrocytes. PLoS One. 2011;6:e26998. https://doi.org/10.1371/journal.pone.0026998.
Dahle MK, Wessel Ø, Timmerhaus G, Nyman IB, Jorgensen SM, Rimstad E, et al. Transcriptome analyses of Atlantic salmon (Salmo salar L.) erythrocytes infected with piscine orthoreovirus (PRV). Fish Shellfish Immunol. 2015;45(2):780–90. https://doi.org/10.1016/j.fsi.2015.05.049.
Puente-Marin S, Nombela I, Ciordia S, Mena MC, Chico V, Coll J, et al. In silico functional networks identified in fish nucleated red blood cells by means of transcriptomic and proteomic profiling. Genes. 2018;9(4):202. https://doi.org/10.3390/genes9040202.
Gavery MR, Nichols KM, Goetz G, Middleton MA, Swanson P. Characterization of genetic and epigenetic variation in sperm and red blood cells from adult hatchery and natural-origin steelhead, Oncorhynchus mykiss. G3: genes. Genomes, Genetics. 2018;8(11):3723–36. https://doi.org/10.1534/g3.118.200458.
Vandeputte M, Gagnaire PA, Allal F. The European sea bass: a key marine fish model in the wild and in aquaculture. Anim Genet. 2019;50(3):195–206. https://doi.org/10.1111/age.12779.
Tine M, Kuhl H, Gagnaire PA, Louro B, Desmarais E, Martins RST, et al. European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation. Nat Commun. 2014;5:5770. https://doi.org/10.1038/ncomms6770.
Navarro-Martín L, Viñas J, Ribas L, Diaz N, Gutierrez A, Di Croce L, et al. DNA methylation of the gonadal aromatase (cyp19a) promoter is involved in temperature-dependent sex ratio shifts in the European sea bass. PLoS Genet. 2011;7:e1002447. https://doi.org/10.1371/journal.pgen.1002447.
Geay F, Zambonino-Infante J, Reinhardt R, Kuhl H, Santigosa E, Cahu C, et al. Characteristics of fads2 gene expression and putative promoter in European sea bass (Dicentrarchus labrax): comparison with salmonid species and analysis of CpG methylation. Mar Genomics. 2012;5:7–13. https://doi.org/10.1016/j.margen.2011.08.003.
Terova G, Díaz N, Rimoldi S, Ceccotti C, Gliozheni E, Piferrer F. Effects of sodium butyrate treatment on histone modifications and the expression of genes related to epigenetic regulatory mechanisms and immune response in European sea bass (Dicentrarchus labrax) fed a plant-based diet. PLoS One. 2016;11(7):e0160332. https://doi.org/10.1371/journal.pone.0160332.
Anastasiadi D, Díaz N, Piferrer F. Small ocean temperature increases elicit stage-dependent changes in DNA methylation and gene expression in a fish, the European sea bass. Sci Rep. 2017;7(1):12401. https://doi.org/10.1038/s41598-017-10861-6.
Anastasiadi D, Vandeputte M, Sánchez-Baizán N, Allal F, Piferrer F. Dynamic epimarks in sex-related genes predict gonad phenotype in the European sea bass, a fish with mixed genetic and environmental sex determination. Epigenetics. 2018;13(9):988–1011. https://doi.org/10.1080/15592294.2018.1529504.
Barrachina F, Anastasiadi D, Jodar M, Castillo J, Estanyol JM, Piferrer F, et al. Identification of a complex population of chromatin-associated proteins in the European sea bass (Dicentrarchus labrax) sperm. Syst Biol Reprod Med. 2018;64(6):502–17. https://doi.org/10.1080/19396368.2018.1482383.
Anastasiadi D, Piferrer F. Epimutations in developmental genes underlie the onset of domestication in farmed European sea bass. Mol Biol Evol. 2019;36(10):2252–64. https://doi.org/10.1093/molbev/msz153.
Anastasiadi D, Piferrer F. A clockwork fish: age prediction using DNA methylation-based biomarkers in the European seabass. Mol Ecol Resour. 2020;20(2):387–97. https://doi.org/10.1111/1755-0998.13111.
FAO. FAO Yearbook of Fishery and Aquaculture Statistics. 2020. Available at: http://www.fao.org/fishery/statistics/yearbook/en (Accessed 27 Jan 2021).
Fanouraki E, Mylonas CC, Papandroulakis N, Pavlidis M. Species specificity in the magnitude and duration of the acute stress response in Mediterranean marine fish in culture. Gen Comp Endocrinol. 2011;173:313–22. https://doi.org/10.1016/j.ygcen.2011.06.004.
Vazzana M, Cammarata M, Cooper EL, Parrinello N. Confinement stress in sea bass (Dicentrarchus labrax) depresses peritoneal leukocyte cytotoxicity. Aquaculture. 2002;210:231–43. https://doi.org/10.1016/S0044-8486(01)00818-3.
Cœurdacier JL, Dutto G, Gasset E, Blancheton JP. Is total serum protein a good indicator for welfare in reared sea bass (Dicentrarchus labrax)? Aquat Living Resour. 2011;24:124–7. https://doi.org/10.1051/alr/2011130.
Samaras A, Dimitroglou A, Sarropoulou E, Papaharisis L, Kottaras L, Pavlidis M. Repeatability of cortisol stress response in the European sea bass (Dicentrarchus labrax) and transcription differences between individuals with divergent responses. Sci Rep. 2016;6:34858. https://doi.org/10.1038/srep34858.
Tognoli C, Rossi F, Di Cola F, Baj G, Tongiorgi E, Terova G, et al. Acute stress alters transcript expression pattern and reduces processing of proBDNF to mature BDNF in Dicentrarchus labrax. BMC Neurosci. 2010;11:4. https://doi.org/10.1186/1471-2202-11-4.
van Gurp TP, Wagemaker NCAM, Wouters B, Vergeer P, Ouborg JNJ, Verhoeven KJS. epiGBS: reference-free reduced representation bisulfite sequencing. Nat Meth. 2016;13:322–4. https://doi.org/10.1038/nmeth.3763.
Allis CD, Jenuwein T. The molecular hallmarks of epigenetic control. Nat Rev Genet. 2016;17(8):487–500. https://doi.org/10.1038/nrg.2016.59.
Schübeler D. Function and information content of DNA methylation. Nature. 2015;517:321–6. https://doi.org/10.1038/nature14192.
Rey S, Boltana S, Vargas R, Roher N, MacKenzie S. Combining animal personalities with transcriptomics resolves individual variation within a wild-type zebrafish population and identifies underpinning molecular differences in brain function. Mol Ecol. 2013;22:6100–15. https://doi.org/10.1111/mec.12556.
Liu J, Dias K, Plagnes-Juan E, Veron V, Panserat S, Marandel L. Long-term programming effect of embryonic hypoxia exposure and high-carbohydrate diet at first feeding on glucose metabolism in juvenile rainbow trout. J Exp Biol. 2017;220(Pt 20):3686–94. https://doi.org/10.1242/jeb.161406.
Wong RY, Lamm MS, Godwin J. Characterizing the neurotranscriptomic states in alternative stress coping styles. BMC Genomics. 2015;16(1):425. https://doi.org/10.1186/s12864-015-1626-x.
Greenwood AK. Peichel CL social regulation of gene expression in threespine sticklebacks. PLoS One. 2015;10(9):e0137726. https://doi.org/10.1371/journal.pone.0137726.
Walsh CA, Engle EC. Allelic diversity in human developmental neurogenetics: insights into biology and disease. Neuron. 2010;68(2):245–53. https://doi.org/10.1016/j.neuron.2010.09.042.
Vlasova-St Louis I, Dickson AM, Bohjanen PR, Wilusz CJ. 2013. CELFish ways to modulate mRNA decay. Biochim Biophys Acta. 2013;1829(6–7):695–707. https://doi.org/10.1016/j.bbagrm.2013.01.001.
Meröndun J, Murray DL, Shafer ABA. Genome-scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx. Mol Ecol. 2019;28(13):3186–96. https://doi.org/10.1111/mec.15131.
Wang J, Xia Y, Li L, Gong DS, Yao Y, Luo HJ, et al. Double restriction-enzyme digestion improves the coverage and accuracy of genome-wide CpG methylation profiling by reduced representation bisulfite sequencing. BMC Genomics. 2013;14:11. https://doi.org/10.1186/1471-2164-14-11.
Krueger F, Andrews SR. Bismark: A flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011;11:1571–2. https://doi.org/10.1093/bioinformatics/btr167.
Wan ZY, Xia JH, Lin G, Wang L, Lin VCL, Yue GH. Genome-wide methylation analysis identified sexually dimorphic methylated regions in hybrid tilapia. Sci Rep. 2016;6:35903. https://doi.org/10.1038/srep35903.
Berbel-Filho WM, Berry N, Rodríguez-Barreto D. Texeira SR, de Leaniz CG, Consuegra S. environmental enrichment induces intergenerational behavioural and epigenetic effects on fish. Mol Ecol. 2020;29(12):2888–299. https://doi.org/10.1111/mec.15481.
Wang J, Liu Y, Jiang S, Li WH, Gui L, Zhou T, et al. Transcriptomic and epigenomic alterations of Nile tilapia gonads sexually reversed by high temperature. Aquaculture. 2019;508:167–77. https://doi.org/10.1016/j.aquaculture.2019.04.073.
Heckwolf MJ, Meyer BS, Häsler R, Hoppner MP, Eizaguirre MC, Reusch TBH. Two different epigenetic information channels in wild three-spined sticklebacks are involved in salinity adaptation. Sci Adv. 2020;6(12):eaaz1138. https://doi.org/10.1126/sciadv.aaz1138.
Hu J, Pérez-Jvostov F, Blondel L, Barrett RDH. Genome-wide DNA methylation signatures of infection status in Trinidadian guppies (Poecilia reticulata). Mol Ecol. 2018;27(15):3087–102. https://doi.org/10.1111/mec.14771.
Ehrlich M, Lacey M. DNA methylation and differentiation: silencing, upregulation and modulation of gene expression. Epigenomics. 2013;5(5):553–68. https://doi.org/10.2217/epi.13.43.
Lev Maor G, Yearim A, Ast G. The alternative role of DNA methylation in splicing regulation. Trends Genet. 2015;31(5):274–80. https://doi.org/10.1016/j.tig.2015.03.002.
Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D’Souza C, Fouse SD, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010;466:253–7. https://doi.org/10.1038/nature09165.
Hing B, Braun P, Cordner ZA, Ewald ER, Moody L, McKane M, et al. 2018. Chronic social stress induces DNA methylation changes at an evolutionary conserved intergenic region in chromosome X. Epigenetics. 2018;13(6):627–41. https://doi.org/10.1080/15592294.2018.1486654.
Suzuki M, Oda M, Ramos MP, Pascual M, Lau K, Stasiek E, et al. Late-replicating heterochromatin is characterized by decreased cytosine methylation in the human genome. Genome Res. 2011;21(11):1833–40. https://doi.org/10.1101/gr.116509.110.
Gavery MR, Nichols KM, Berejikian BA, Tatara CP, Goetz GW, Dickey JT, et al. Temporal dynamics of DNA methylation patterns in response to rearing juvenile steelhead (Oncorhynchus mykiss) in a hatchery versus simulated stream environment. Genes. 2019;10:356. https://doi.org/10.3390/genes10050356.
Baumgart M, Groth M, Priebe S, Savino A, Testa G, Dix A, et al. RNA-seq of the aging brain in the short-lived fish N. furzeri - conserved pathways and novel genes associated with neurogenesis. Aging Cell. 2014;13(6):965–74. https://doi.org/10.1111/acel.12257.
Nie XG. Differential expression of Bmp2, Bmp4 and Bmp3 in embryonic development of mouse anterior and posterior palate. Chin Med J. 2005;118(20):1710–6. 16313756.
Walker MB, Miller CT, Coffin Talbot J, Stock DW, Kimmel CB. Zebrafish furin mutants reveal intricacies in regulating Endothelin1 signaling in craniofacial patterning. Dev Biol. 2006;295(1):194–205. https://doi.org/10.1016/j.ydbio.2006.03.028.
Kumari PK, Ali A, Singh SK, Chaurasia A, Raman R. Genetic heterogeneity in Van der Woude syndrome: identification of NOL4 and IRF6 haplotype from the noncoding region as candidates in two families. J Genet. 2018;97(1):275–85 PMID: 2966634.
Gustafsson MK, Pan H, Pinney DF, Liu YL, Lewandowski A, Epstein DJ, et al. 2002. Myf5 is a direct target of long-range Shh signaling and Gli regulation for muscle specification. Genes Dev. 2002;16(1):114–26. https://doi.org/10.1101/gad.940702.
Lin CY, Lee HC, Chen HC, Hsieh CC, Tsai HJ. Normal function of Myf5 during gastrulation is required for pharyngeal arch cartilage development in zebrafish embryos. Zebrafish. 2013;10(4):486–99. https://doi.org/10.1089/zeb.2013.09.
Friocourt F, Chédotal A. The Robo3 receptor, a key player in the development, evolution, and function of commissural systems. Dev Neurobiol. 2017;77(7):876–90. https://doi.org/10.1002/dneu.22478.
Kamm GB, Pisciottano F, Kliger R, Franchini LF. The developmental brain gene NPAS3 contains the largest number of accelerated regulatory sequences in the human genome. Mol Biol Evol. 2013;30(5):1088–102. https://doi.org/10.1093/molbev/mst023.
McGhee KE, Bell AM. Paternal care in a fish: epigenetics and fitness enhancing effects on offspring anxiety. Proc Biol Sci. 2014;281(1794):20141146. https://doi.org/10.1098/rspb.2014.1146.
Oliveira RF, Simões JM, Teles MC, Oliveira CR, Becker JD, Lopez JS. Assessment of fight outcome is needed to activate socially driven transcriptional changes in the zebrafish brain. Proc Natl Acad Sci U S A. 2016;113(5):E654–61. https://doi.org/10.1073/pnas.1514292113.
Auta J, Smith RC, Dong E, Tueting P, Sershen H, Boules S, et al. DNA-methylation gene network dysregulation in peripheral blood lymphocytes of schizophrenia patients. Schizophr Res. 2013;150(1):312–8. https://doi.org/10.1016/j.schres.2013.07.030.
Blaze J, Wang J, Ho L, Mendelev N, Haghighi F, Pasinetti GM. Polyphenolic compounds alter stress-induced patterns of global DNA methylation in brain and blood. Mol Nutr Food Res. 2018;62(8):e1700722. https://doi.org/10.1002/mnfr.20170072.
Vinkers CH, Geuze E, van Rooij SJH, Kennis M, Schür RR, Nispeling DM, et al. Successful treatment of post-traumatic stress disorder reverses DNA methylation marks. Mol Psychiatry. 2019. https://doi.org/10.1038/s41380-019-0549-3.
Arockiaraj AI, Liu D, Shaffer JR, Koleck TA, Crago EA, Weeks DE, Conley YP. Methylation data processing protocol and comparison of blood and cerebral spinal fluid following aneurysmal subarachnoid hemorrhage. Front Genet. 2020;11:671. https://doi.org/10.3389/fgene.2020.00671.
Pértille F, Ibelli AMG, Sharif ME, Poleti MD, Fröhlich AS, Rezaei S, et al. Putative epigenetic biomarkers of stress in red blood cells of chickens reared across different biomes. Front Genet. 2020;11:508809. https://doi.org/10.3389/fgene.2020.508809.
Lubin FD, Roth TL, Sweatt JD. Epigenetic regulation of BDNF gene transcription in the consolidation of fear memory. J Neurosci. 2008;28(42):10576–86. https://doi.org/10.1523/JNEUROSCI.1786-08.2008.
Coley EJL, Demaestri C, Ganguly P, Honeycutt JA, Peterzell S, Rose N, et al. Cross-generational transmission of early life stress effects on HPA regulators and Bdnf are mediated by sex, lineage, and upbringing. Front Behav Neurosci. 2019;13:101. https://doi.org/10.3389/fnbeh.2019.00101.
Numakawa T, Odaka H, Adachi N. Actions of brain-derived neurotrophic factor and glucocorticoid stress in neurogenesis. Int J Mol Sci. 2017;18(11):2312. https://doi.org/10.3390/ijms18112312.
Notaras M, van den Buuse M. Neurobiology of BDNF in fear memory, sensitivity to stress, and stress-related disorders. Mol Psychiatry. 2020;25:2251–74. https://doi.org/10.1038/s41380-019-0639-2.
Pavlidis M, Theodoridi A, Tsalafouta A. Neuroendocrine regulation of the stress response in adult zebrafish, Danio rerio. Prog Neuro-Psychopharmacol Biol Psychiatry. 2015;60:121–31. https://doi.org/10.1016/j.pnpbp.2015.02.014.
Vindas MA, Fokos S, Pavlidis M, Hoglund E, Dyonisopoulou S, Ebbesson LOE, et al. Early life stress induces long-term changes in limbic areas of a teleost fish: the role of catecholamine systems in stress coping. Sci Rep. 2018;8(1):5638. https://doi.org/10.1038/s41598-018-23950-x.
Sadoul B, Alfonso S, Bessa E, Bouchareb A, Blondeau-Bidet E, Clair P, et al. Enhanced brain expression of genes related to cell proliferation and neural differentiation is associated with cortisol receptor expression in fishes. Gen Comp Endocrinol. 2018;267:76–81. https://doi.org/10.1016/j.ygcen.2018.06.001.
Samaras A, Espirito Santo C, Papandroulakis N, Mitrizakis N, Pavlidis M, Hoglund E, et al. Allostatic load and stress physiology in European seabass (Dicentrarchus labrax L.) and gilthead seabream (Sparus aurata L.). Front Endocrinol. 2018;9:451. https://doi.org/10.3389/fendo.2018.00451.
Alfonso S, Sadoul B, Gesto M, Joassard L, Chatain B, Bégout ML. Coping styles in European sea bass: the link between boldness, stress response and neurogenesis. Physiol Behav. 2019;207:67–85. https://doi.org/10.1016/j.physbeh.2019.04.020.
Sassone-Corsi P. The cyclic AMP pathway. Cold Spring Harb Perspect Biol. 2012;4(12):a011148. https://doi.org/10.1101/cshperspect.a011148.
Lee EH, Seo SR. Neuroprotective roles of pituitary adenylate cyclase-activating polypeptide in neurodegenerative diseases. BMB Rep. 2014;47(7):369–75. https://doi.org/10.5483/bmbrep.2014.47.7.086.
Marosi K, Mattson MP. BDNF mediates adaptive brain and body responses to energetic challenges. Trends Endocrinol Metab. 2014;25(2):89–98. https://doi.org/10.1016/j.tem.2013.10.006.
Bennett MR, Lagopoulos J. Stress and trauma: BDNF control of dendritic-spine formation and regression. Prog Neurobiol. 2014;112:80–99. https://doi.org/10.1016/j.pneurobio.2013.10.005.
Buck JM, O’Neill HC, Stitzel JA. Developmental nicotine exposure elicits multigenerational disequilibria in proBDNF proteolysis and glucocorticoid signaling in the frontal cortices, striata, and hippocampi of adolescent mice. Biochem Pharmacol. 2019;168:438–51. https://doi.org/10.1016/j.bcp.2019.08.003.
Zhu B, Zhao L, Luo D, Xue DM, Tan T, Dong ZF, et al. Furin promotes dendritic morphogenesis and learning and memory in transgenic mice. Cell Mol Life Sci. 2018;75(13):2473–88. https://doi.org/10.1007/s00018-017-2742-3.
Marcinkiewicz M, Seidah NG, Chrétien M. Implications of the subtilisin/kexin-like precursor convertases in the development and function of nervous tissues. Acta Neurobiol Exp. 1996;56(1):287–98 PMID: 8787188.
Chen Y, Zhang J, Deng M. Furin mediates brain-derived neurotrophic factor upregulation in cultured rat astrocytes exposed to oxygen-glucose deprivation. J Neurosci Res. 2015;93(1):189–94. https://doi.org/10.1002/jnr.23455.
Yeh FC, Kao CF, Kuo PH. Explore the features of brain-derived neurotrophic factor in mood disorders. PLoS One. 2015;10(6):e0128605. https://doi.org/10.1371/journal.pone.0128605.
Pang PT, Teng HK, Zaitsev E, Woo NT, Sakata K, Zhen SH, et al. Cleavage of proBDNF by tPA/plasmin is essential for long-term hippocampal plasticity. Science. 2004;306:487–91. https://doi.org/10.1126/science.1100135.
McEwen BS. Central effects of stress hormones in health and disease: understanding the protective and damaging effects of stress and stress mediators. Eur J Pharmacol. 2008;583(2–3):174–85. https://doi.org/10.1016/j.ejphar.2007.11.071.
Jeanneteau FD, Lambert WM, Ismaili N, Bath KG, Lee FS, Garabedian MJ, et al. BDNF and glucocorticoids regulate corticotrophin-releasing hormone (CRH) homeostasis in the hypothalamus. Proc Natl Acad Sci U S A. 2012;109(4):1305–10. https://doi.org/10.1073/pnas.1114122109.
Gray K, Ellis V. Activation of pro-BDNF by the pericellular serine protease plasmin. FEBS Lett. 2008;582(6):907–10. https://doi.org/10.1016/j.febslet.2008.02.026.
Wang N, Zhang L, Miles L, Hoover-Plow J. Plasminogen regulates pro-opiomelanocortin processing. J Thromb Haemost. 2004;2(5):785–96. https://doi.org/10.1111/j.1538-7836.2004.00694.x.
Liu X. SLC family transporters. Adv Exp Med Biol. 2019;1141:101–202. https://doi.org/10.1007/978-981-13-7647-4_3.
Strazielle N, Ghersi-Egea JF. Efflux transporters in blood-brain interfaces of the developing brain. Front Neurosci. 2015;9:21. https://doi.org/10.3389/fnins.2015.00021.
Couroussé T, Bacq A, Belzung C, Guiard B, Balasse L, Louis F, et al. Brain organic cation transporter 2 controls response and vulnerability to stress and GSK3β signaling. Mol Psychiatry. 2015;20(7):889–900. https://doi.org/10.1038/mp.2014.86.
Couroussé T, Gautron S. Role of organic cation transporters (OCTs) in the brain. Pharmacol Ther. 2015;146:94–103. https://doi.org/10.1016/j.pharmthera.2014.09.008.
Gai Y, Liu Z, Cervantes-Sandoval I, Davis RL. Drosophila SLC22A transporter is a memory suppressor gene that influences cholinergic neurotransmission to the mushroom bodies. Neuron. 2016;90(3):581–95. https://doi.org/10.1016/j.neuron.2016.03.017.
Erbel-Sieler C, Dudley C, Zhou Y, Wu XL, Estill SJ, Han T, et al. Behavioral and regulatory abnormalities in mice deficient in the NPAS1 and NPAS3 transcription factors. Proc Natl Acad Sci U S A. 2004;101(37):13648–53. https://doi.org/10.1073/pnas.0405310101.
Luoma LM, Macintyre G, Tibbo PG, Wild TC, Colman I, Purdon SE. NPAS3 exonic SNP genotype is linked to working memory performance in healthy young adults. Psychiatry Res. 2018;265:263–4. https://doi.org/10.1016/j.psychres.2018.04.063.
Yang D, Zhang W, Padhiar A, Yue Y, Shi YH, Zheng TZ, et al. NPAS3 regulates transcription and expression of VGF: implications for neurogenesis and psychiatric disorders. Front Mol Neurosci. 2016;9:109. https://doi.org/10.3389/fnmol.2016.00109.
Lavedan C, Licamele L, Volpi S, Hamilton J, Heaton C, Mack K, et al. 2009. Association of the NPAS3 gene and five other loci with response to the antipsychotic iloperidone identified in a whole genome association study. Mol Psychiatry. 2009;14(8):804–19. https://doi.org/10.1038/mp.2008.56.
Walch L. Emerging role of the scaffolding protein Dlg1 in vesicle trafficking. Traffic. 2013;14:964–73. https://doi.org/10.1111/tra.12089.
Liu Y, Cui L, Schwarz MK, Fog Y, Schlüter OM. Adrenergic gate release for spike timing-dependent synaptic potentiation. Neuron 18. 2017;93(2):394–408. https://doi.org/10.1016/j.neuron.2016.12.039.
Jourdi H, Kabbaj M. Acute BDNF treatment upregulates GluR1-SAP97 and GluR2-GRIP1 interactions: implications for sustained AMPA receptor expression. PLoS One. 2013;8(2):e57124. https://doi.org/10.1371/journal.pone.0057124.
Singh MD, Jensen M, Lasser M, Huber E, Yusuff T, Pizzo L, et al. NCBP2 modulates neurodevelopmental defects of the 3q29 deletion in Drosophila and Xenopus laevis models. PLoS Genet. 2020;16(2):e1008590. https://doi.org/10.1371/journal.pgen.1008590.
Lauschke VM, Ivanov M, Ingelman-Sundberg M. Pitfalls and opportunities for epigenomic analyses focused on disease diagnosis, prognosis, and therapy. Trends Pharmacol Sci. 2017;38(9):765–70. https://doi.org/10.1016/j.tips.2017.05.007.
van Hamburg JP, Tas SW. Molecular mechanisms underpinning T helper 17 cell heterogeneity and functions in rheumatoid arthritis. J Autoimmun. 2018;87:69–81. https://doi.org/10.1016/j.jaut.2017.12.006.
Haasch D, Berg C, Clampit JE, Pederson T, Frost L, Kroeger P, et al. PKCtheta is a key player in the development of insulin resistance. Biochem Biophys Res Commun. 2006;343(2):361–8. https://doi.org/10.1016/j.bbrc.2006.02.177.
Zanin-Zhorov A, Kumari S, Hippen KL, Merkel SC, MacMillan ML, Blazar BR, Dustin ML. Human in vitro-induced regulatory T cells display Dlgh1dependent and PKC-θ restrained suppressive activity. Sci Rep. 2017;7(1):4258. https://doi.org/10.1038/s41598-017-04053-5 Erratum in: Sci Rep. 2020;10(1): 3317.
Erion DM, Kotas ME, McGlashon J, Yonemitsu S, Hsiao JJ, Nagai Y, et al. cAMP-responsive element-binding protein (CREB)-regulated transcription coactivator 2 (CRTC2) promotes glucagon clearance and hepatic amino acid catabolism to regulate glucose homeostasis. J Biol Chem. 2013;288(22):16167–76. https://doi.org/10.1074/jbc.M113.460246.
Anyamaneeratch K, Rojvirat P, Sukjoi W, Jitrapakdee S. Insights into transcriptional regulation of hepatic glucose production. Int Rev Cell Mol Biol. 2015;318:203–53. https://doi.org/10.1016/bs.ircmb.2015.05.004.
Hogan MF, Ravnskjaer K, Matsumura S, Huising MO, Hull RL, Kahn SE, et al. Hepatic insulin resistance following chronic activation of the CREB Coactivator CRTC2. J Biol Chem. 2015;290(43):25997–6006. https://doi.org/10.1074/jbc.M115.679266.
Hill MJ, Suzuki S, Segars JH, Kino T. CRTC2 is a Coactivator of GR and couples GR and CREB in the regulation of hepatic gluconeogenesis. Mol Endocrinol. 2016;30(1):104–17. https://doi.org/10.1210/me.2015-1237.
Metzger DCH, Schulte PM. The DNA methylation landscape of stickleback reveals patterns of sex chromosome evolution and effects of environmental salinity. Genome Biol Evol. 2018;10(3):775–85. https://doi.org/10.1093/gbe/evy034.
Sagonas K, Meyer BS, Kaufmann J, Lenz TL, Häsler R, Eizaguirre C. Experimental parasite infection causes genome-wide changes in DNA methylation. Mol Biol Evol. 2020;37(8):2287–99. https://doi.org/10.1093/molbev/msaa084.
Xiu Y, Shao C, Zhu Y, Li YZ, Gan T, Xu WT, et al. Differences in DNA methylation between disease-resistant and disease-susceptible Chinese tongue sole (Cynoglossus semilaevis) families. Front Genet. 2019;10:847. https://doi.org/10.3389/fgene.2019.00847.
Triantaphyllopoulos KA, Ikonomopoulos I, Bannister AJ. Epigenetics and inheritance of phenotype variation in livestock. Epigenetics Chromatin. 2016;9:31. https://doi.org/10.1186/s13072-016-0081-5.
Burggren WW. Dynamics of epigenetic phenomena: intergenerational and intragenerational phenotype 'washout'. J Exp Biol. 2015;218:80–7. https://doi.org/10.1242/jeb.107318.
Trollope AF, Mifsud KR, Saunderson EA, Reul JMHM. Molecular and epigenetic mechanisms underlying cognitive and adaptive responses to stress. Epigenomes. 2017;1(3):17. https://doi.org/10.3390/epigenomes1030017.
Suzuki M, Liao W, Wos F, Johnston AD, DeGrazia J, Ishii J, et al. Whole-genome bisulfite sequencing with improved accuracy and cost. Genome Res. 2018;28(9):1364–71. https://doi.org/10.1101/gr.232587.117.
Werner O, Prudencio ÁS, de la Cruz-Martinez E, Nieto-Lugilde M, Martinez-Gomez P, Ros RM. A cost-reduced variant of epi-genotyping by sequencing for studying DNA methylation in non-model organisms. Front Plant Sci. 2020;11:694. https://doi.org/10.3389/fpls.2020.00694.
Prudencio ÁS, Werner O, Martínez-García PJ, Dicenta F, Ros RM, Martinez-Gomez P. DNA methylation analysis of dormancy release in almond (Prunus dulcis) flower buds using epi-genotyping by sequencing. Int J Mol Sci. 2018;19(11):3542. https://doi.org/10.3390/ijms19113542.
Crotti M, Adams CE, Elmer KR. Population genomic SNPs from epigenetic RADs: gaining genetic and epigenetic data from a single established next-generation sequencing approach. Methods Ecol Evol. 2020;11:839–49. https://doi.org/10.1111/2041-210X.13395.
Suravajhala P, Kogelman LJ, Kadarmideen HN. Multi-omic data integration and analysis using systems genomics approaches: methods and applications in animal production, health and welfare. Genet Sel Evol. 2016;48(1):38. https://doi.org/10.1186/s12711-016-0217-x.
Schwartz TS. The promises and the challenges of integrating multi-omics and systems biology in comparative stress biology. Integr Comp Biol. 2020;60(1):89–97. https://doi.org/10.1093/icb/icaa026.
Gu H, Smith ZD, Bock C, Boyle P, Gnirke A, Meissner A. Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nat Protocols. 2011;6:468–81. https://doi.org/10.1038/nprot.2010.190.
Lepais O, Weir JT. SimRAD: an R package for simulation-based prediction of the number of loci expected in RADseq and similar genotyping by sequencing approaches. Mol Ecol Resour. 2014;14:1314–21. https://doi.org/10.1111/1755-0998.12273.
Catchen J, Hohenlohe P. Bassham S, Amoros A, Cresko WA. Stacks: an analysis tool set for population genomics. Mol Ecol 2013; 22(11): 3124–3140. doi: https://doi.org/10.1111/mec.12354.
Akalin A, Kormaksson M, Li S, Garrett-Bakelman FE, Figueroa ME, Melnick A, et al. MethylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012;13:R87. https://doi.org/10.1186/gb-2012-13-10-R87.
Edgar R, Domarachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–1.
We thank N. Wagemaker at the Radboud University (Nijmegen, The Netherlands) for his support, and Prof. C. Grunau at the University of Perpignan Via Domitia (France) for hosting M.V.K at various steps of this work. S. Ben Chehida helped during the last stages of manuscript preparation. We thank people from the GenSeq technical facility and from the Montpellier Bioinformatics Biodiversity (MBB) facility at ISEM for continuous support for early step of molecular analysis and computer resources availability, respectively. These facilities are supported by the Laboratoire d’Excellence (LabEx) CeMEB and by ANR “Investissements d’Avenir” program (ANR-10-LABX-04-01). Comments by two anonymous reviewers improved the final version of the manuscript.
This research was funded by the EU ERA-NET FP7 programme for the Cooperation in Fisheries, Aquaculture and Seafood Processing (COFASP) and the Agence Nationale pour la Recherche (ANR; grant number: ANR-16-0001-COFA). Montpellier GenomiX received financial support from the France Génomique National infrastructure, also funded as part of “Investissement d’Avenir” program managed by ANR (grant number: ANR-10-INBS-09). These funding sources had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All experiments were performed in accordance with relevant guidelines and regulations. Nireus S.A. research facilities are certified and have obtained the codes for the rearing and use of fish for scientific purposes (EL04-BIOexp-01). All procedures on fish used in this study were approved by the Departmental Animal Care Committee following the Three Rs principle, in accordance with Greek (PD 56/2013) and EU (Directive 63/2010) legislation on the care and use of experimental animals. All sample procedures were performed by FELASA certified researchers (http://www.felasa.eu/).
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Krick, M.V., Desmarais, E., Samaras, A. et al. Family-effects in the epigenomic response of red blood cells to a challenge test in the European sea bass (Dicentrarchus labrax, L.). BMC Genomics 22, 111 (2021). https://doi.org/10.1186/s12864-021-07420-9