Genome-wide methylation analysis identified sexually dimorphic methylated regions in hybrid tilapia

Development of DNA Insertion-specific Markers Based on the Intergenic Region of Oplegnathus punctatus Cdkn1/srsf3 for Sex Identification

Spotted knifejaw (Oplegnathus punctatus) is a marine economic fish with high food and ecological value, and its growth process has obvious male and female sexual dimorphism, with males growing significantly faster than females. However, the current sex identification technology is not yet mature, which will limit the growth rate of O. punctatus aquaculture and the efficiency of separate sex breeding, so the development of efficient sex molecular markers is imperative. This study identified a 926 bp DNA insertion fragment in the cdkn1/srsf3 intergenic region of O. punctatus males through whole-genome scanning, comparative genomics, and structural variant analysis. A pair of primers was designed based on the insertion information of the Y chromosome intergenic region in male individuals. Agarose gel electrophoresis revealed the amplification of two DNA fragments, 1118 bp and 192 bp, in male O. punctatus individuals. The 926 bp fragment was identified as the insertion in the intergenic region of cdkn1/srsf3 in males, while only a single 192 bp DNA fragment was amplified in females. The biological sex of the individuals identified in this manner was consistent with their known phenotypic sex. In this study, we developed a method to detect DNA insertion variants in the intergenic region of O. punctatus. Additionally, we introduced a new DNA marker for the rapid identification of the sex of O. punctatus, which enhances detection efficiency. The text has important reference significance and application value in sex identification, all-male breeding, and lineage selection. It provides new insights into the regulation of variation in the intergenic region of cdkn1/srsf3 genes and the study of RNA shearing.

Integrative analysis of miRNA-mRNA expression in the brain during high temperature-induced masculinization of female Nile tilapia (Oreochromis niloticus)

Temperature is one of the most important non-genetic sex differentiation factors for fish. The technique of high temperature-induced sex reversal is commonly used in Nile tilapia (Oreochromis niloticus) culture, although the molecular regulatory mechanisms involved in this process remain unclear. The brain is an essential organ for the regulation of neural signals involved in germ cell differentiation and gonad development. To investigate the regulatory roles of miRNAs-mRNAs in the conversion of female to male Nile tilapia gender under high-temperature stress, we compared RNA-Seq data from brain tissues between a control group (28 °C) and a high temperature-treated group (36 °C). The result showed that a total of 123,432,984 miRNA valid reads, 288,202,524 mRNA clean reads, 1128 miRNAs, and 32,918 mRNAs were obtained. Among them, there were 222 significant differentially expressed miRNAs (DE miRNAs) and 810 differentially expressed mRNAs (DE mRNAs) between the two groups. Eight DE miRNAs and eight DE mRNAs were randomly selected, and their expression patterns were validated by qRT-PCR. The miRNA-mRNA co-expression network demonstrated that 40 DE miRNAs targeted 136 protein-coding genes. Functional enrichment analysis demonstrated that these genes were involved in several gonadal differentiation pathways, including the oocyte meiosis signaling pathway, progesterone-mediated oocyte maturation signaling pathway, cell cycle signaling pathway and GnRH signaling pathway. Then, an interaction network was constructed for 8 miRNAs (mir-137-5p, let-7d, mir-1388-5p, mir-124-4-5p, mir-1306, mir-99, mir-130b and mir-21) and 10 mRNAs (smc1al, itpr2, mapk1, ints8, cpeb1b, bub1, fbxo5, mmp14b, cdk1 and hrasb) involved in the oocyte meiosis signaling pathway. These findings provide novel information about the mechanisms underlying miRNA-mediated sex reversal in female Nile tilapia.

Exploring the Effects of Rearing Densities on Epigenetic Modifications in the Zebrafish Gonads

Rearing density directly impacts fish welfare, which, in turn, affects productivity in aquaculture. Previous studies have indicated that high-density rearing during sexual development in fish can induce stress, resulting in a tendency towards male-biased sex ratios in the populations. In recent years, research has defined the relevance of the interactions between the environment and epigenetics playing a key role in the final phenotype. However, the underlying epigenetic mechanisms of individuals exposed to confinement remain elucidated. By using zebrafish (Danio rerio), the DNA methylation promotor region and the gene expression patterns of six genes, namely dnmt1, cyp19a1a, dmrt1, cyp11c1, hsd17b1, and hsd11b2, involved in the DNA maintenance methylation, reproduction, and stress were assessed. Zebrafish larvae were subjected to two high-density conditions (9 and 66 fish/L) during two periods of overlapping sex differentiation of this species (7 to 18 and 18 to 45 days post-fertilization, dpf). Results showed a significant masculinization in the populations of fish subjected to high densities from 18 to 45 dpf. In adulthood, the dnmt1 gene was differentially hypomethylated in ovaries and its expression was significantly downregulated in the testes of fish exposed to high-density. Further, the cyp19a1a gene showed downregulation of gene expression in the ovaries of fish subjected to elevated density, as previously observed in other studies. We proposed dnmt1 as a potential testicular epimarker and the expression of ovarian cyp19a1a as a potential biomarker for predicting stress originated from high densities during the early stages of development. These findings highlight the importance of rearing densities by long-lasting effects in adulthood conveying cautions for stocking protocols in fish hatcheries.

Estrodiol-17β and aromatase inhibitor treatment induced alternations of genome-wide DNA methylation pattern in Takifugu rubripes gonads

Estradiol-17β (E2) and aromatase inhibitor (AI) exposure can change the phenotypic sex of fish gonads. To investigated whether alterations in DNA methylation is involved in this process, the level of genome-wide DNA methylation in Takifugu rubripes gonads was quantitatively analyzed during the E2-induced feminization and AI-induced masculinization processes in this study. The methylation levels of the total cytosine (C) in control-XX(C-XX), control-XY (C-XY), E2-treated-XY (E-XY) and AI-treated-XX (AI-XX) were 9.11%, 9.19%, 8.63% and 9.23%, respectively. In the C-XX vs C-XY comparison, 4,196 differentially methylated regions (DMRs) overlapped with the gene body of 2,497 genes and 608 DMRs overlapped with the promoter of 575 genes. In the E-XY vs C-XY comparison, 6,539 DMRs overlapped with the gene body of 3,416 genes and 856 DMRs overlapped with the promoter of 776 genes. In the AI-XX vs C-XX comparison, 2,843 DMRs overlapped with the gene body of 1,831 genes and 461 DMRs overlapped with the promoter of 421 genes. Gonadal genomic methylation mainly occurred at CG sites and the genes that overlapped with DMRs on CG context were most enriched in the signaling pathways related to gonad differentiation, such as the Wnt, TGF-β, MAPK, CAM and GnRH pathways. The DNA methylation levels of steroid synthesis genes and estrogen receptor genes promoter or gene body were negative correlated with their expression. After bisulfite sequencing verification, the DNA methylation level of the amhr2 promoter in XY was increased after E2 treatment, which consistent with the data from the genome-wide DNA methylation sequencing. In C-XY group, the expression of amhr2 was significantly higher than that in E-XY (p < 0.05). Additionally, dnmt1, which is responsible for methylation maintenance, expressed at similar level in four groups (p > 0.05). dnmt3, tet2, and setd1b, which were responsible for methylation modification, expressed at significantly higher levels in E-XY compared to the C-XY (p < 0.05). Dnmt3 and tet2 were expressed at significantly higher levels in AI-XX than that in C-XX (p < 0.05). These results indicated that E2 and AI treatment lead to the aberrant genome-wide DNA methylation level and expression level of dnmt3, tet2, and setd1b in T. rubripes gonad.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Comparative genome-wide methylation analysis reveals epigenetic regulation of muscle development in grass carp (Ctenopharyngodon idellus) fed with whole faba bean

Crisp grass carp (CGC), the most representative improved varieties of grass carp (GC), features higher muscle hardness after feeding faba bean (Vicia faba L.) for 90-120 days. DNA methylation, a most widely studied epigenetic modification, plays an essential role in muscle development. Previous studies have identified numerous differentially expressed genes (DEGs) between CGC and GC. However, it remains unknown if the expression levels of these DEGs are influenced by DNA methylation. In the present study, we performed a comprehensive analysis of DNA methylation profiles between CGC and GC, and identified important candidate genes related to muscle development coupled with the transcriptome sequencing data. A total of 9,318 differentially methylated genes (DMGs) corresponding to 155,760 differentially methylated regions (DMRs) were identified between the two groups under the CG context in promoter regions. Combined with the transcriptome sequencing data, 14 key genes related to muscle development were identified, eight of which (gsk3b, wnt8a, wnt11, axin2, stat1, stat2, jak2, hsp90) were involved in muscle fiber hyperplasia, six of which (tgf-β1, col1a1, col1a2, col1a3, col4a1, col18a1) were associated with collagen synthesis in crisp grass carp. The difference of methylation levels in the key genes might lead to the expression difference, further resulting in the increase of muscle hardness in crisp grass carp. Overall, this study can help further understand how faba bean modulates muscle development by the epigenetic modifications, providing novel insights into the texture quality improvement in other aquaculture fish species by nutritional programming.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Genome-wide DNA methylation patterns harbour signatures of hatchling sex and past incubation temperature in a species with environmental sex determination

Conservation of thermally sensitive species depends on monitoring organismal and population-level responses to environmental change in real time. Epigenetic processes are increasingly recognized as key integrators of environmental conditions into developmentally plastic responses, and attendant epigenomic data sets hold potential for revealing cryptic phenotypes relevant to conservation efforts. Here, we demonstrate the utility of genome-wide DNA methylation (DNAm) patterns in the face of climate change for a group of especially vulnerable species, those with temperature-dependent sex determination (TSD). Due to their reliance on thermal cues during development to determine sexual fate, contemporary shifts in temperature are predicted to skew offspring sex ratios and ultimately destabilize sensitive populations. Using reduced-representation bisulphite sequencing, we profiled the DNA methylome in blood cells of hatchling American alligators (Alligator mississippiensis), a TSD species lacking reliable markers of sexual dimorphism in early life stages. We identified 120 sex-associated differentially methylated cytosines (DMCs; FDR < 0.1) in hatchlings incubated under a range of temperatures, as well as 707 unique temperature-associated DMCs. We further developed DNAm-based models capable of predicting hatchling sex with 100% accuracy (in 20 training samples and four test samples) and past incubation temperature with a mean absolute error of 1.2°C (in four test samples) based on the methylation status of 20 and 24 loci, respectively. Though largely independent of epigenomic patterning occurring in the embryonic gonad during TSD, DNAm patterns in blood cells may serve as nonlethal markers of hatchling sex and past incubation conditions in conservation applications. These findings also raise intriguing questions regarding tissue-specific epigenomic patterning in the context of developmental plasticity.

Effects of salinity stress on methylation of the liver genome and complement gene in large yellow croaker (Larimichthys crocea)

Salinity is an important environmental factor that affects the yield and quality of large yellow croaker (Larimichthys crocea) during aquaculture. Here, whole-genome bisulfite sequencing (WGBS), RNA-seq, bisulfite sequencing PCR (BSP), quantitative real-time PCR (qPCR), and dual luciferase reporter gene detection technologies were used to analyze the DNA methylation characteristics and patterns of the liver genome, the expression and methylation levels of important immune genes in large yellow croaker in response to salinity stress. The results of WGBS showed that the cytosine methylation of CG type was dominant, CpGIsland and repeat regions were important regions where DNA methylation occurred, and the DNA methylation in upstream 2k (2000bp upstream of the promoter) and repeat regions had different changes in the liver tissue of large yellow croaker in the response to the 12‰, 24‰, 36‰ salinity stress of 4 w (weeks). In the combined analysis of WGBS and transcriptome, the complement and coagulation cascade pathways were significantly enriched, in which the complement-related genes C7, C3, C5, C4, C1R, MASP1, and CD59 were mainly changed in response to salinity stress. In the studied area of MASP1 gene promoter, the methylation levels of many CpG sites as well as total cytosine were strongly negatively correlated with mRNA expression level. Methylation function analysis of MASP1 promoter further proved that DNA methylation could inhibit the activity of MASP1 promoter, indicating that salinity may affect the expressions of complement-related genes by DNA methylation of gene promoter region.

Identification of Male-Specific Molecular Markers by Recombination of RhoGEF10 Gene in Spotted Knifejaw (Oplegnathus punctatus)

The spotted knifejaw (Oplegnathus punctatus) is a marine economic fish with high ecological value, food value, and fishing value, and its growth has obvious sex dimorphism. The rapid identification of its sex is beneficial to the development of sex determination and breeding. In this study, the method of comparative genomics and PCR amplification was used to further establish a rapid detection method for the recombinant RhoGEF10 gene in O. punctatus, which can quickly, accurately, and efficiently identify the sex of the O. punctatus to be tested. The homologous comparison results of male and female individuals showed that the DNA fragment length of the RhoGEF10 gene on the X1 chromosome was 326 bp, and the DNA fragment length on the Y chromosome was 879 bp. Therefore, it can be concluded that there is an insert fragment of 553 bp on the Y chromosome. PCR amplification results showed that the two DNA fragments of 879 bp and 326 bp were amplified in the Y chromosome and X1 chromosome of the male O. punctatus (X1X2Y), respectively, and the 879 bp fragment was a unique marker fragment of the recombinant RhoGEF10 gene; The female O. punctatus (X1X1X2X2) only a single DNA fragment of 326 bp was amplified. At the same time, the inserted fragment of the male individual resulted in partial inactivation of the RhoGEF10 protein, which in turn resulted in a slowing of peripheral nerve conduction velocity and thinning of the myelin sheath in male O. punctatus. The method shortens the time for accurate identification of the O. punctatus RhoGEF10 gene recombination and improves the detection efficiency. It is of great significance and application value in the research of nerve conduction and myelin development, male and female sex identification, the preparation of high male seedlings, and family selection based on the RhoGEF10 gene in the O. punctatus.

The Role of Transposable Elements in Sexual Development

Up to 50% of most mammalian genomes are made up of transposable elements (TEs) that have the potential to mobilize around the genome. Despite this prevalence, research on TEs is only beginning to gain traction within the field of neuroscience. While TEs have long been regarded as “junk” or parasitic DNA, it has become evident that they are adaptive DNA and RNA regulatory elements. In addition to their vital role in normal development, TEs can also interact with steroid receptors, which are key elements to sexual development. In this review, we provide an overview of the involvement of TEs in processes related to sexual development- from TE activity in the germline to TE accumulation in sex chromosomes. Moreover, we highlight sex differences in TE activity and their regulation of genes related to sexual development. Finally, we speculate on the epigenetic mechanisms that may govern TEs’ role in sexual development. In this context, we emphasize the need to further the understanding of sexual development through the lens of TEs including in a variety of organs at different developmental stages, their molecular networks, and evolution.

Integrated Analyses of DNA Methylation and Gene Expression of Rainbow Trout Muscle under Variable Ploidy and Muscle Atrophy Conditions

Rainbow trout, Oncorhynchus mykiss, is an important cool, freshwater aquaculture species used as a model for biological research. However, its genome reference has not been annotated for epigenetic markers affecting various biological processes, including muscle growth/atrophy. Increased energetic demands during gonadogenesis/reproduction provoke muscle atrophy in rainbow trout. We described DNA methylation and its associated gene expression in atrophying muscle by comparing gravid, diploid females to sterile, triploid females. Methyl Mini-seq and RNA-Seq were simultaneously used to characterize genome-wide DNA methylation and its association with gene expression in rainbow trout muscle. Genome-wide enrichment in the number of CpGs, accompanied by depleted methylation levels, was noticed around the gene transcription start site (TSS). Hypermethylation of CpG sites within ±1 kb on both sides of TSS (promoter and gene body) was weakly/moderately associated with reduced gene expression. Conversely, hypermethylation of the CpG sites in downstream regions of the gene body +2 to +10 kb was weakly associated with increased gene expression. Unlike mammalian genomes, rainbow trout gene promotors are poor in CpG islands, at <1% compared to 60%. No signs of genome-wide, differentially methylated (DM) CpGs were observed due to the polyploidy effect; only 1206 CpGs (0.03%) were differentially methylated, and these were primarily associated with muscle atrophy. Twenty-eight genes exhibited differential gene expression consistent with methylation levels of 31 DM CpGs. These 31 DM CpGs represent potential epigenetic markers of muscle atrophy in rainbow trout. The DM CpG-harboring genes are involved in apoptosis, epigenetic regulation, autophagy, collagen metabolism, cell membrane functions, and Homeobox proteins. Our study also identified genes explaining higher water content and modulated glycolysis previously shown as characteristic biochemical signs of rainbow trout muscle atrophy associated with sexual maturation. This study characterized DNA methylation in the rainbow trout genome and its correlation with gene expression. This work also identified novel epigenetic markers associated with muscle atrophy in fish/lower vertebrates.

Identification of an Epigenetically Marked Locus within the Sex Determination Region of Channel Catfish

Channel catfish has an XY sex determination system. However, the X and Y chromosomes harbor an identical gene content of 950 genes each. In this study, we conducted comparative analyses of methylome and transcriptome of genetic males and genetic females before gonadal differentiation to provide insights into the mechanisms of sex determination. Differentially methylated CpG sites (DMCs) were predominantly identified on the sex chromosome, most notably within the sex determination region (SDR), although the overall methylation profiles across the entire genome were similar between genetic males and females. The drastic differences in methylation were located within the SDR at nucleotide position 14.0–20.3 Mb of the sex chromosome, making this region an epigenetically marked locus within the sex determination region. Most of the differentially methylated CpG sites were hypermethylated in females and hypomethylated in males, suggesting potential involvement of methylation modification in sex determination in channel catfish. Along with the differential methylation in the SDR, a number of differentially expressed genes within the SDR were also identified between genetic males and females, making them potential candidate genes for sex determination and differentiation in channel catfish.

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.

Differences in DNA methylation between slow and fast muscle in Takifugu rubripes

Fish skeletal muscle is comprised of fast muscle (FM) and slow muscle (SM), which constitutes 60% of total the body mass. Fish skeletal muscle can affect fish swimming activity, which is important for aquaculture due to its growth-potentiating effects. DNA methylation can influence gene expression level. We previously identified multiple differentially expressed genes (DEGs) between FM and SM in Takifugu rubripes. However, it is unknown if the expression levels of these DEGs are influenced by DNA methylation. In the present study, we used DNA methylation sequencing to study the DNA methylation profiles of FM and SM in T. rubripes. SM had higher overall methylation levels than FM. A total of 8479 differentially methylated genes (DMGs) and 3407 DMGs containing differentially methylated regions (DMRs) in the promoter regions between FM and SM were identified. After enrichment analysis, we found functionally relevant DMGs between FM and SM, including Kapca, Plcd3a, Plcd1, Pi3k, Tsp4b and Pgfrb in the hedgehog signaling pathway and phosphatidylinositol (PI)-related pathways. Due to the different methylation levels of these genes between FM and SM, the expression levels of Kapca, Plcd3a, Plcd1, Pi3k, and Tsp4b were higher in FM and Pgfrb was higher in SM. There were differences in the hedgehog signaling pathway and PI-related pathways between FM and SM. In SM, the cytokine-cytokine receptor interaction promoted focal adhesion, while ECM-receptor interactions promoted focal adhesion in FM. These results provide information regarding the difference between FM and SM in T. rubripes.

Genome-wide DNA methylation profiles provide insight into epigenetic regulation of red and white muscle development in Chinese perch Siniperca chuatsi

Fish skeletal muscles are composed of spatially well-separated fiber types, namely, red and white muscles with different physiological functions and metabolism. To compare the DNA methylation profiles of the two types of muscle tissues and identify potential candidate genes for the muscle growth and development under epigenetic regulation, genome-wide DNA methylation of the red and white muscle in Chinese perch Siniperca chuatsi were comparatively analyzed using bisulfate sequencing methods. An average of 0.9 billion 150-bp paired-end reads were obtained, of which 86% were uniquely mapped to the genome. Methylation mostly occurred at CG sites at a ratio of 94.43% in the red muscle and 93.16% in the white muscle. The mean methylation levels at C-sites were 5.95% in red muscle and 5.83% in white muscle, whereas the mean methylation levels of CG, CHG, and CHH were 73.23%, 0.62%, and 0.67% in red muscle, and 71.01%, 0.62%, and 0.67% in white muscle, respectively. A total of 4192 differentially methylated genes (DMGs) were identified significantly enriched in cell signaling pathways related to skeletal muscle differentiation and growth. Various muscle-related genes, including myosin gene isoforms and regulatory factors, are differentially methylated in the promoter region between the red and white muscles. Further analysis of the transcriptional expression of these genes showed that the muscle regulatory factors (myf5, myog, pax3, pax7, and twitst2) and myosin genes (myh10, myh16, myo18a, myo7a, myo9a, and myl3) were differentially expressed between the two kinds of muscles, consistent with the DNA methylation analysis results. ELISA assays confirmed that the level of 5mC in red muscle was significantly higher than in white muscle (P < 0.05). The RT-qPCR assays revealed that the expression levels of the three DNA methylation transferase (dnmt) subtypes, dnmt1, dnmt3ab, and dnmt3bb1, were significantly higher in red muscle than in white muscle. The higher DNA methylation levels in the red muscle may result from higher DNA methylation transferase expression in the red muscles. Thus, this study might provide a theoretical foundation to better understand epigenetic regulation in the growth and development of red and white muscles in animals, at least in Chinese perch fish.

Epigenetic Distribution of Recombinant Plant Chromosome Fragments in a Human-Arabidopsis Hybrid Cell Line

Methylation systems have been conserved during the divergence of plants and animals, although they are regulated by different pathways and enzymes. However, studies on the interactions of the epigenomes among evolutionarily distant organisms are lacking. To address this, we studied the epigenetic modification and gene expression of plant chromosome fragments (~30 Mb) in a human-Arabidopsis hybrid cell line. The whole-genome bisulfite sequencing results demonstrated that recombinant Arabidopsis DNA could retain its plant CG methylation levels even without functional plant methyltransferases, indicating that plant DNA methylation states can be maintained even in a different genomic background. The differential methylation analysis showed that the Arabidopsis DNA was undermethylated in the centromeric region and repetitive elements. Several Arabidopsis genes were still expressed, whereas the expression patterns were not related to the gene function. We concluded that the plant DNA did not maintain the original plant epigenomic landscapes and was under the control of the human genome. This study showed how two diverging genomes can coexist and provided insights into epigenetic modifications and their impact on the regulation of gene expressions between plant and animal genomes.

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.

Holo-Omics: Integrated Host-Microbiota Multi-omics for Basic and Applied Biological Research

From ontogenesis to homeostasis, the phenotypes of complex organisms are shaped by the bidirectional interactions between the host organisms and their associated microbiota. Current technology can reveal many such interactions by combining multi-omic data from both hosts and microbes. However, exploring the full extent of these interactions requires careful consideration of study design for the efficient generation and optimal integration of data derived from (meta)genomics, (meta)transcriptomics, (meta)proteomics, and (meta)metabolomics. In this perspective, we introduce the holo-omic approach that incorporates multi-omic data from both host and microbiota domains to untangle the interplay between the two. We revisit the recent literature on biomolecular host-microbe interactions and discuss the implementation and current limitations of the holo-omic approach. We anticipate that the application of this approach can contribute to opening new research avenues and discoveries in biomedicine, biotechnology, agricultural and aquacultural sciences, nature conservation, as well as basic ecological and evolutionary research.

Early leucine programming on protein utilization and mTOR signaling by DNA methylation in zebrafish (Danio rerio)

Background

Early nutritional programming affects a series of metabolism, growth and development in mammals. Fish also exhibit the developmental plasticity by early nutritional programming. However, little is known about the effect of early amino acid programming on growth and metabolism.

Methods

In the present study, zebrafish (Danio rerio) was used as the experimental animal to study whether early leucine stimulation can programmatically affect the mechanistic target of rapamycin (mTOR) signaling pathway, growth and metabolism in the later life, and to undercover the mechanism of epigenetic regulation. Zebrafish larvas at 3 days post hatching (dph) were raised with 1.0% leucine from 3 to 13 dph during the critical developmental stage, then back to normal water for 70 days (83 dph).

Results

The growth performance and crude protein content of zebrafish in the early leucine programming group were increased, and consistent with the activation of the mTOR signaling pathway and the high expression of genes involved in the metabolism of amino acid and glycolipid. Furthermore, we compared the DNA methylation profiles between the control and leucine-stimulated zebrafish, and found that the methylation levels of CG-differentially methylated regions (DMGs) and CHH-DMGs of genes involved in mTOR signaling pathway were different between the two groups. With quantitative PCR analysis, the decreased methylation levels of CG type of Growth factor receptor-bound protein 10 (Grb10), eukaryotic translation initiation factor 4E (eIF4E) and mTOR genes of mTOR signaling pathway in the leucine programming group, might contribute to the enhanced gene expression.

Conclusions

The early leucine programming could improve the protein synthesis and growth, which might be attributed to the methylation of genes in mTOR pathway and the expression of genes involved in protein synthesis and glycolipid metabolism in zebrafish. These results could be beneficial for better understanding of the epigenetic regulatory mechanism of early nutritional programming.

Sex Chromosomes and Sex Phenotype Contribute to Biased DNA Methylation in Mouse Liver

Sex biases in the genome-wide distribution of DNA methylation and gene expression levels are some of the manifestations of sexual dimorphism in mammals. To advance our understanding of the mechanisms that contribute to sex biases in DNA methylation and gene expression, we conducted whole genome bisulfite sequencing (WGBS) as well as RNA-seq on liver samples from mice with different combinations of sex phenotype and sex-chromosome complement. We compared groups of animals with different sex phenotypes, but the same genetic sexes, and vice versa, same sex phenotypes, but different sex-chromosome complements. We also compared sex-biased DNA methylation in mouse and human livers. Our data show that sex phenotype, X-chromosome dosage, and the presence of Y chromosome shape the differences in DNA methylation between males and females. We also demonstrate that sex bias in autosomal methylation is associated with sex bias in gene expression, whereas X-chromosome dosage-dependent methylation differences are not, as expected for a dosage-compensation mechanism. Furthermore, we find partial conservation between the repertoires of mouse and human genes that are associated with sex-biased methylation, an indication that gene function is likely to be an important factor in this phenomenon.

Population epigenetic divergence exceeds genetic divergence in the Eastern oyster Crassostrea virginica in the Northern Gulf of Mexico

Populations may respond to environmental heterogeneity via evolutionary divergence or phenotypic plasticity. While evolutionary divergence occurs through DNA sequence differences among populations, plastic divergence among populations may be generated by changes in the epigenome. Here, we present the results of a genome-wide comparison of DNA methylation patterns and genetic structure among four populations of Eastern oyster (Crassostrea virginica) in the northern Gulf of Mexico. We used a combination of restriction site-associated DNA sequencing (RADseq) and reduced representation bisulfite sequencing (RRBS) to explore population structure, gene-wide averages of F ST, and DNA methylation differences between oysters inhabiting four estuaries with unique salinity profiles. This approach identified significant population structure despite a moderately low F ST (0.02) across the freshwater boundary of the Mississippi river, a finding that may reflect recent efforts to restore oyster stock populations. Divergence between populations in CpG methylation was greater than for divergence in F ST, likely reflecting environmental effects on DNA methylation patterns. Assessment of CpG methylation patterns across all populations identified that only 26% of methylated DNA was intergenic; and, only 17% of all differentially methylated regions (DMRs) were within these same regions. DMRs within gene bodies between sites were associated with genes known to be involved in DNA damage repair, ion transport, and reproductive timing. Finally, when assessing the correlation between genomic variation and DNA methylation between these populations, we observed population-specific DNA methylation profiles that were not directly associated with single nucleotide polymorphisms or broader gene-body mean F ST trends. Our results suggest that C. virginica may use DNA methylation to generate environmentally responsive plastic phenotypes and that there is more divergence in methylation than divergence in allele frequencies.

Immune genes, IL1β and Casp9, show sexual dimorphic methylation patterns in zebrafish gonads

There is crosstalk between the immune and reproductive systems in which sexual dimorphism is a common pattern in vertebrates. In recent years, epigenetics has emerged as a way to study the molecular mechanisms involved in gonadal development, those responsible for integrating environmental information that contribute to assigning a specific sexual phenotype (either an ovary or a testis). The knowledge of epigenetic mechanisms in certain molecular processes allows the development of epigenetic markers. In fish gonads, the existence of reproduction-immune system interactions is known, although the epigenetic mechanisms involved are far from clear. Here, we used the zebrafish (Danio rerio) as a model to study the DNA methylation patterns in gonads of two well-known innate immune genes: IL1β and Casp9. DNA methylation levels were studied by a candidate gene approach at single nucleotide resolution and gene expression analyses were also carried out. Results showed that there was clear sexual dimorphism in the DNA methylation levels of the two immune genes studied, being significantly higher in the testes when compared to the ovaries. In summary, and although further research is needed, this paper presents sexual dimorphic methylation patterns of two immune-related genes, thus sex-biased differences in methylation profiles should considered when analyzing immune responses in fish. Data showed here can help to develop epimarkers with forthcoming applications in livestock and fish farming production, for example, in immune fish diseases or sexual control programs as epigenetic molecular tools to predict environmental pressure in the gonads.

Changes of DNA Methylation Pattern in Metabolic Pathways Induced by High-Carbohydrate Diet Contribute to Hyperglycemia and Fat Deposition in Grass Carp (Ctenopharyngodon idellus)

Although studies have determined that epigenetics plays an essential role in regulating metabolism in mammals, research on nutrition-related DNA methylation remains to be lacking in teleosts. In the present study, we provided a hepatic whole-genome DNA methylation analysis in grass carp fed with moderate- or excessive-carbohydrate-level diet. Although a high-carbohydrate (HC) diet significantly changed the mRNA expression levels of metabolic genes, it did not affect the global genomic DNA methylation levels in grass carp liver. However, compared with the control group, 3,972 genes were hyper-methylated and 2,904 genes were hypo-methylated in the promoter region. Meanwhile, 10,711 genes were hyper-methylated and 6,764 genes were hypo-methylated in the gene body region in the HC group. These differentially methylated genes (DMGs) were enriched in multiple pathways, including carbohydrate metabolism, insulin pathway, lipid metabolism, and adipocytokine signaling pathway. In addition, the variations in DNA methylation significantly regulated the transcription levels of key genes of metabolism, which could affect the glucose concentrations and the lipid deposition of grass carp. Furthermore, we compared the DNA methylation alterations of genes in glucose metabolism and obesity pathways of grass carp with those of mammalian models in different nutritional states. The results showed that most of the DMGs in grass carp were also regulated by DNA methylation in mammals when the nutritional state changed. The findings revealed more differentially methylated regions and candidate genes for glucose metabolism and broken species boundaries.

Identification and Characterization of a Non-muscular Myostatin in the Nile Tilapia

The growth and differentiation factor Myostatin (MSTN, also known as GDF8) negatively regulates skeletal muscle development and growth in vertebrates. Most fish genomes contain two or more mstn genes, which are expressed in muscle and other tissues. Yet, in the genome of Nile tilapia (Oreochromis niloticus), which is one of the world's most important aquaculture fish species, only one mstn gene has previously been identified. Here, we identify a second mstn gene in Nile tilapia. We show that it clusters phylogenetically with other piscine mstn2 genes and that it shares chromosomal synteny with the human and zebrafish orthologs. We further show that mstn2 is not expressed in red or white muscles of Nile tilapia, but rather that its main site of expression is the brain. To determine which physiological functions are correlated with mstn expression, adult Nile tilapia were exposed to various environmental conditions and their effect on mstn1 and mstn2 expression in the brain and muscles was measured using real-time PCR. We found that the centrally- and muscle-expressed mstn genes differ in their responsiveness to diverse challenges, suggesting differential gene- and tissue-specific regulation of their expression. Metabolic and stress marker analyses showed that the altered mstn expression is not regulated by classical stress response. Taken together, our findings expand the understanding of the MSTN system in Nile tilapia and provide evolutionary insight into its function.

Differences in DNA Methylation Between Disease-Resistant and Disease-Susceptible Chinese Tongue Sole (Cynoglossus semilaevis) Families

DNA methylation, the most widely studied and most well-understood epigenetic modification, has been reported to play crucial roles in diverse processes. Although it has been found that DNA methylation can modulate the expression of immune-related genes in teleosts, a systemic analysis of epigenetic regulation on teleost immunity has rarely been performed. In this research, we employed whole-genome bisulfite sequencing to investigate the genome-wide DNA methylation profiles in select disease-resistant Cynoglossus semilaevis (DR-CS, family 14L006) and disease-susceptible C. semilaevis (DS-CS, family 14L104) against Vibrio harveyi infection. The results showed that following selective breeding, DR-CS had higher DNA methylation levels and different DNA methylation patterns, with 3,311 differentially methylated regions and 6,456 differentially methylated genes. Combining these data with the corresponding transcriptome data, we identified several immune-related genes that exhibited differential expression levels that were modulated by DNA methylation. Specifically, DNA methylation of tumor necrosis factor–like and lipopolysaccharide-binding protein-like was significantly correlated with their expression and significantly contributed to the disease resistance of the selected C. semilaevis family. In conclusion, we suggest that artificial selection for disease resistance in Chinese tongue sole causes changes in DNA methylation levels in important immune-related genes and that these epigenetic changes are potentially involved in multiple immune responses in Chinese tongue sole.

Changes in DNA methylation during epigenetic-associated sex reversal under low temperature in Takifugu rubripes

DNA methylation has frequently been implicated in sex determination and differentiation in teleost species. In order to detect the DNA methylation patterns established during sexual differentiation in tiger pufferfish T. rubripes, we performed comprehensive whole genome methylation sequencing and analyses of the gonads of male, female, and pseudo male. We obtained a total of 33.12, 32.44, and 31.60 Gb clean data for male, female, and pseudo male, with a sequencing depth of 66.44×, 60.47× and 54.86×, respectively. The methylation level of cytosine (C) residues in the genomic DNA from gonads was 11.016%, 10.428%, and 11.083% in male, female, and pseudo male, respectively. More than 65% of C methylation was at CpG sites, and less than 1% was at CHG and CHH sites. In each regulatory element, there were low methylation levels on both sides of the transcription start site, and higher methylation levels in exons, introns, and downstream of genes. The highest mCpG was on chromosome 8 and the lowest mCpG was on chromosome 5. Comparisons of whole-genome DNA methylation between pairs of samples revealed that there were 3,173 differentially methylated regions (DMRs) between female and male, and 3,037 DMRs between male and pseudo male, corresponding to 0.232% and 0.223% of the length of the genome, respectively. There were only 1,635 DMRs between female and pseudo male, representing 0.127% of the length of the genome. A number of differentially methylated genes (DMGs) related to sex determination and differentiation were selected, such as amhr2 and pfcyp19a. After Bisulfite Sequencing PCR (BSP) verification, amhr2 was exhibited low methylation level in normal males and pseudo male, and high methylation level in normal females but pfcyp19a showed low methylation level in normal females and high methylation level in normal males and pseudo males. These results provide information about the molecular epigenetic mechanisms of DNA methylation during low-temperature induced masculinization of tiger pufferfish, and increase our understanding of the mechanisms of sex determination and differentiation in this important aquaculture fish species.

Differences in the fast muscle methylome provide insight into sex-specific epigenetic regulation of growth in Nile tilapia during early stages of domestication

Growth is a complex trait whose variability within a population cannot be explained solely by genetic variation. Epigenetic regulation is often suggested as an important factor shaping the phenotype, but its association with growth can be highly context- and species-dependent. Nevertheless, the mechanisms involved in epigenetic regulation of growth in fish are poorly understood. We have used reduced representation bisulphite sequencing to determine the genome-wide CpG methylation patterns in male and female Nile tilapia of different sizes but at the same early stage of domestication. The average CpG methylation level in the reduced genome representation was 63% across groups but many sites displayed group-specific methylation patterns. The number of differentially methylated (DM) CpGs was much higher when the interaction between sex and weight was included rather than when these factors were considered separately. There were 1128 DM CpGs between large and small females and 970 DM CpGs between large and small males. We have found many growth-related genes associated with DM CpGs, namely map3k5 and akt3 in females and gadd45g and ppargc1a in males. Only 5% of CpG locations associated with growth were common to both sexes. In particular, the autophagy-related gene atg14 displayed a high association of methylation with growth exclusively in males. The sexually dimorphic association between atg14 methylation and growth may uncover novel metabolic mechanisms at play during mouth brooding in Nile tilapia females. Taken together, our data suggest that epigenetic regulation of growth in Nile tilapia involves different gene networks in males and females.

Transcriptome Analysis Identified Genes for Growth and Omega-3/-6 Ratio in Saline Tilapia

Growth and omega-3/-6 ratio are important traits in aquaculture. The mechanisms underlying quick growth and high omega-3/-6 ratio in fish are not fully understood. The consumption of the meat of tilapia suffers a bad reputation due to its low omega-3/-6 ratio. To facilitate the improvement of these traits and to understand more about the mechanisms underlying quick growth and high omega-3/-6 ratio, we conducted transcriptome analysis in the muscle and liver of fast- and slow-growing hybrid saline tilapia generated by crossing Mozambique tilapia and red tilapia. A transcriptome with an average length of 963 bp was generated by using 486.65 million clean 100 bp paired-end reads. A total of 42,699 annotated unique sequences with an average length of 3.4 kb were obtained. Differentially expressed genes (DEGs) in the muscle and liver were identified between fast- and slow-growing tilapia. Pathway analysis classified these genes into many pathways. Ten genes, including foxK1, sparc, smad3, usp38, crot, fadps, sqlea, cyp7b1, impa1, and gss, from the DEGs were located within QTL for growth and omega-3, which were previously detected content in tilapia, suggesting that these ten genes could be important candidate genes for growth and omega-3 fatty acid content. Analysis of SNPs in introns 1 and 2 of foxK1 revealed that the SNPs were significantly associated with growth and omega-3/-6 ratio. This study lays the groundwork for further investigation of the molecular mechanisms underlying the phenotypic variation of these two traits and provides SNPs for selecting these traits at fingerling stage.

The DNA Methylation Landscape of Stickleback Reveals Patterns of Sex Chromosome Evolution and Effects of Environmental Salinity

Epigenetic mechanisms such as DNA methylation are a key component of dosage compensation on sex chromosomes and have been proposed as an important source of phenotypic variation influencing plasticity and adaptive evolutionary processes, yet little is known about the role of DNA methylation in an ecological or evolutionary context in vertebrates. The threespine stickleback (Gasterosteus aculeatus) is an ecological and evolutionary model system that has been used to study mechanisms involved in the evolution of adaptive phenotypes in novel environments as well as the evolution heteromorphic sex chromosomes and dosage compensation in vertebrates. Using whole genome bisulfite sequencing, we compared genome-wide DNA methylation patterns between threespine stickleback males and females and between stickleback reared at different environmental salinities. Apparent hypermethylation of the younger evolutionary stratum of the stickleback X chromosome in females relative to males suggests a potential role of DNA methylation in the evolution of heteromorphic sex chromosomes. We also demonstrate that rearing salinity has genome-wide effects on DNA methylation levels, which has the potential to lead to the accumulation of epigenetic variation between natural populations in different environments.

Evolutionary conservation of DNA methylation in CpG sites within ultraconserved noncoding elements

Ultraconserved noncoding elements (UCNEs) constitute less than 1 Mb of vertebrate genomes and are impervious to accumulating mutations. About 4000 UCNEs exist in vertebrate genomes, each at least 200 nucleotides in length, sharing greater than 95% sequence identity between human and chicken. Despite extreme sequence conservation over 400 million years of vertebrate evolution, we show both ordered interspecies and within-species interindividual variation in DNA methylation in these regions. Here, we surveyed UCNEs with high CpG density in 56 species finding half to be intermediately methylated and the remaining near 0% or 100%. Intermediately methylated UCNEs displayed a greater range of methylation between mouse tissues. In a human population, most UCNEs showed greater variation than the LINE1 transposon, a frequently used epigenetic biomarker. Global methylation was found to be inversely correlated to hydroxymethylation across 60 vertebrates. Within UCNEs, DNA methylation is flexible, conserved between related species, and relaxed from the underlying sequence selection pressure, while remaining heritable through speciation.

Applied Hologenomics: Feasibility and Potential in Aquaculture

Aquaculture will play an essential role in feeding a growing human population, but several biological challenges impede sustainable growth of production. Emerging evidence across all areas of life has revealed the importance of the intimate biological interactions between animals and their associated gut microbiota. Based on challenges in aquaculture, we leverage current knowledge in molecular biology and host microbiota interactions to propose an applied holo-omic framework that integrates molecular data including genomes, transcriptomes, epigenomes, proteomes, and metabolomes for analyzing fish and their gut microbiota as interconnected and coregulated systems. With an eye towards aquaculture, we discuss the feasibility and potential of our holo-omic framework to improve growth, health, and sustainability in any area of food production, including livestock and agriculture.

Phylogeny, expression patterns and regulation of DNA Methyltransferases in early development of the flatfish, Solea senegalensis

Background

The identification of DNA methyltransferases (Dnmt) expression patterns during development and their regulation is important to understand the epigenetic mechanisms that modulate larval plasticity in marine fish. In this study, dnmt1 and dnmt3 paralogs were identified in the flatfish Solea senegalensis and expression patterns in early developmental stages and juveniles were determined. Additionally, the regulation of Dnmt transcription by a specific inhibitor (5-aza-2′-deoxycytidine) and temperature was evaluated.

Results

Five paralog genes of dnmt3, namely dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1 and dnmt3bb.2 and one gene for dnmt1 were identified. Phylogenetic analysis revealed that the dnmt gene family was highly conserved in teleosts and three fish-specific genes, dnmt3aa, dnmt3ba and dnmt3bb.2 have evolved. The spatio-temporal expression patterns of four dnmts (dnmt1, dnmt3aa, dnmt3ab and dnmt3bb.1) were different in early larval stages although all of them reduced expression with the age and were detected in neural organs and dnmt3aa appeared specific to somites. In juveniles, the four dnmt genes were expressed in brain and hematopoietic tissues such as kidney, spleen and gills. Treatment of sole embryos with 5-aza-2′-deoxycytidine down-regulated dntm1 and up-regulated dntm3aa. Moreover, in lecithotrophic larval stages, dnmt3aa and dnmt3ab were temperature sensitive and their expression was higher in larvae incubated at 16 °C relative to 20 °C.

Conclusion

Five dnmt3 and one dnmt1 paralog were identified in sole and their distinct developmental and tissue-specific expression patterns indicate that they may have different roles during development. The inhibitor 5-aza-2′-deoxycytidine modified the transcript abundance of dntm1 and dntm3aa in embryos, which suggests that a regulatory feedback mechanism exists for these genes. The impact of thermal regime on expression levels of dnmt3aa and dnmt3ab in lecithotrophic larval stages suggests that these paralogs might be involved in thermal programing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-017-0154-0) contains supplementary material, which is available to authorized users.