Next Article in Journal
Uterine Aging and Reproduction: Dealing with a Puzzle Biologic Topic
Next Article in Special Issue
Temporal and Spatial Signatures of Scylla paramamosain Transcriptome Reveal Mechanistic Insights into Endogenous Ovarian Maturation under Risk of Starvation
Previous Article in Journal
Targeting TRP Channels for Pain, Itch and Neurogenic Inflammation
Previous Article in Special Issue
Role of Sox3 in Estradiol-Induced Sex Reversal in Pelodiscus sinensis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 1
10.3390/ijms25010321
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads

by
Mingxi Hou
1,
Qi Wang
1,
Ran Zhao
1,
Yiming Cao
1,
Jin Zhang
1,
Xiaoqing Sun
1,
Shuangting Yu
1,2,
Kaikuo Wang
3,
Yingjie Chen
3,
Yan Zhang
1 and
Jiongtang Li
1,*
1
Key Laboratory of Aquatic Genomics, Ministry of Agriculture and Rural Affairs, Beijing Key Laboratory of Fishery Biotechnology, Chinese Academy of Fishery Sciences, Beijing 100141, China
2
Chinese Academy of Agricultural Sciences, Beijing 100081, China
3
National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai 201306, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(1), 321; https://doi.org/10.3390/ijms25010321
Submission received: 27 November 2023 / Revised: 22 December 2023 / Accepted: 23 December 2023 / Published: 25 December 2023
(This article belongs to the Special Issue Molecular Advance on Reproduction and Fertility of Aquatic Animals)
Browse Figures
Review Reports Versions Notes

Abstract

:
Epigenetic modifications are critical in precisely regulating gene expression. The common carp (Cyprinus carpio) is an economically important fish species, and females exhibit faster growth rates than males. However, the studies related to epigenetic modifications in the common carp gonads are limited. In this study, we conducted the Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) and Bisulfite sequencing (BS-seq) to explore the roles of epigenetic modifications in the common carp gonads. We identified 84,207 more accessible regions and 77,922 less accessible regions in ovaries compared to testes, and some sex-biased genes showed differential chromatin accessibility in their promoter regions, such as sox9a and zp3. Motif enrichment analysis showed that transcription factors (TFs) associated with embryonic development and cell proliferation were heavily enriched in ovaries, and the TFs Foxl2 and SF1 were only identified in ovaries. We also analyzed the possible regulations between chromatin accessibility and gene expression. By BS-seq, we identified 2087 promoter differentially methylated genes (promoter-DMGs) and 5264 gene body differentially methylated genes (genebody-DMGs) in CG contexts. These genebody-DMGs were significantly enriched in the Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway, indicating that methylation in gene body regions could play an essential role in sex maintenance, just like methylation in promoter regions. Combined with transcriptomes, we revealed that the expression of dmrtb1-like, spag6, and fels was negatively correlated with their methylation levels in promoter regions. Our study on the epigenetic modifications of gonads contributes to elucidating the molecular mechanism of sex differentiation and sex maintenance in the common carp.

1. Introduction

Sex differentiation refers to the process by which undifferentiated gonads develop into ovaries or testes after sex determination [1]. Although the triggers for sex determination vary among species, ranging from genetic to environmental factors [2], the downstream genes and signaling pathways involved in sex differentiation are conserved, and their expression patterns show sexual dimorphism. For instance, dmrt1, gsdf, amh, and sox9 are highly expressed in testes, while foxl2 and cyp19a1a are predominantly expressed in ovaries [3,4,5]. In addition to ovary and testis development, sex differentiation is accompanied by a range of sexual dimorphisms, such as the growth rate, color, and size [6]. The mono-sex population is desirable in the aquaculture industry, as females and males exhibit differences in economic values [7].
Epigenetic modifications are momentous in sex determination and differentiation [8]. Chromatin accessibility refers to the degree of physical contact between nuclear macromolecules and DNA, of which altering is a kind of significant epigenetic modification [9]. Numerous factors, such as the conformation and composition of nucleosomes, histone architectural proteins, RNA polymerases, and transcription factors (TFs), can influence DNA accessibility [10,11]. Active regulatory elements and motifs tend to be depleted for nucleosomes and are accessible to transcription machinery and enzymes [12]. Therefore, insights from these open regions will advance our comprehension of transcription factor binding and the regulatory potential of genetic loci. The Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) is a sensitive method for obtaining the open chromatin regions [13]. Currently, the roles of chromatin accessibility in fish sex differentiation have only received limited attention.
DNA methylation is the most extensively studied and understood epigenetic modification in teleost [14]. Many previous studies have demonstrated that methylation in promoter regions is critical for programming sexual fate and maintaining sexual identity by negatively regulating sex-biased gene expression [15,16,17,18]. For instance, the methylation status of dmrt1 and foxl2 promoter regions governs sex differentiation in Japanese flounder (Paralichthys olivaceus) [19,20]. DNA methylation machinery displays dynamic shifts during sex reversal in bluehead wrasse (Thalassoma bifasciatum) [21].
The common carp (Cyprinus carpio) is a crucial fish cultivated worldwide for food consumption [22]. In 2021, the global output of common carp was about 4.41 million metric tons, accounting for approximately 5.2% of the freshwater fishery production [23]. The common carp conducts the XX/XY sex-determination system, with females growing faster than males [24,25]. Breeding the all-female carp population can improve efficiency and increase the yield for aquaculture industries. The study on reproduction is significant for improving productivity and sustainable production.
Previous studies have identified many differentially expressed genes between ovaries and testes in the common carp [26]. A detailed description of gonadal histological and sex-biased gene expression during early sex differentiation was provided [27]. The deletion of cyp17a1 results in XX common carps developing into males and being capable of normal spermatogenesis [28]. In addition, plenty of transcriptome sequencing has uncovered numerous sex-biased genes [23,29,30]. However, the function of epigenetic modifications in sex differentiation remains unknown in the common carp.
In the present study, we generated chromatin accessibility landscapes and DNA methylome profiles of the common carp gonads to understand the molecular mechanism of sex differentiation. We identified differential epigenetic modifications between ovaries and testes and explored the possible epigenetic-mediated regulation of gene expression. To the best of our knowledge, this is the first time the functions of epigenetic modifications in the common carp gonads have been investigated using ATAC-seq and BS-seq.

2. Results

2.1. Profiles of Chromatin Accessibility in Ovaries and Testes

We performed ATAC-seq to explore the chromatin accessibility of the common carp ovaries and testes. We generated 223.52 million raw reads. After quality filtering, 220.62 million clean reads were retained and mapped to the reference genome [31] (Table S1). The quality of ATAC-seq libraries can be assessed based on the distribution of reads near transcription start sites (TSS). As shown in Figure 1A, the reads were highly enriched in TSS, implying a low background in our profiles. Pearson’s correlation heatmap indicated a high level of consistency within groups but a low consistency across groups (Figure 1B).
Peak calling was conducted with an adjusted p ≤ 0.05. We finally identified 120,701 peaks in the ovaries and 193,469 in the testes. We grouped these peaks into four categories based on their genome regions: the promoter (3 kb upstream of TSS), exon, intron, and intergenic region. We found that 34.32% of peaks in the ovaries and 22.13% in the testes were enriched in promoter regions, much higher than the proportion of promoter regions in the whole genome. Moreover, 37.06% of peaks in the ovaries were enriched in gene bodies and 24.44% were enriched in intergenic regions, while 41.12% of peaks in testes were enriched in gene bodies and 32.71% were enriched in intergenic regions (Figure 1C).

2.2. Differential Peaks between Ovaries and Testes

We performed a comparative analysis of peaks between the ovaries and testis groups. Among all the peaks, we identified 67,753 peaks unique to the ovary group and 140,521 unique to the testis group (Figure 2A; Table S2). The reads enriched in peaks were normalized using the values of Reads Per Million mapped reads (RPM), and the RPM values of peaks in the ovaries and testes were present utilizing a scatter plot (Figure 2B). To obtain high-confidence differential peaks, we set a fold change of RPM greater than two as the threshold. Finally, we found that 84,207 peaks exhibited more enrichment in the ovaries and 77,922 peaks exhibited more enrichment in the testes. These genome regions in which differential peaks were located were regarded as differentially accessible regions (Figure 2C).

2.3. Functional Annotation of Chromatin Accessibility in Common Carp Gonads

The open chromatin regions around genes can predict the interaction sites between regulatory elements and genes [32]. All peaks were assigned to their nearest genes. The common carp is an allo-tetraploid species derived from recent whole genome duplication [33]. We found that 47.42% (123,881) of peaks were in the subgenome A, and 50.33% (131,482) were in the subgenome B. There were significantly enriched peaks in the promoter regions of sex-biased genes, such as sox9a and zp3. The details and annotations of all peaks are shown in Table S2.
Genes associated with differentially accessible regions were denoted as differentially accessible genes (DAGs). In comparing ovaries to testes, more accessible regions were associated with 34,037 genes, and less accessible regions were associated with 24,999 genes. The overlapping 17,605 genes were associated with increased and decreased chromatin accessibility regions simultaneously, implying that they could be co-regulated by several regulatory elements (Table S3). These DAGs were then subjected to enrichment analysis. The most significant enrichment GO terms were cellular, metabolic, and biosynthetic processes at the Biological Processes level, intracellular at the Cell Component level, and organic cyclic compound binding and heterocyclic compound at the Molecular Function level (Figure 3A). Although the KEGG analysis of all DAGs showed that no significant pathway was enriched, genes more accessible in testes showed significant enrichment in the TGF-beta signaling pathway, neuroactive ligand–receptor interaction, and calcium signaling pathway (Figure 3B).

2.4. Motif Enrichment Analysis

TFs regulate gene expression by recognizing and binding specific motifs. We analyzed TF binding motifs in the open chromatin regions to mine potential regulations, and the most enriched motifs of each sample were shown in a bubble diagram (Figure 4A). In comparing ovaries and testes, the top ten TFs enriched in ovary-biased peaks included CTCF (CCCTC-binding factor), ETS1 (E26 transformation-specific transcription factor 1), and GABPA (GA-binding protein A) (Figure 4B). These TFs were essential for embryonic development and cell proliferation. Meanwhile, the testis-biased peaks were highly enriched in TBP3 (TATA-binding protein 3), ZEB1 (Zinc finger E-Box binding homeobox 1), RXR (Retinoid X receptor), and other TFs (Figure 4C). All differential TF binding motifs between ovaries and testes are listed in Table S4 (p ≤ 0.01). The motifs of the well-known TFs FoxL2 (Forkhead box L2) (p = 1.00 × 10−178) and SF1 (Steroidogenic Factor 1) (p = 1.00 × 10−345), playing a pivotal role in regulating the synthesis of steroid hormones, were enriched in the ovaries (Table S4).

2.5. Relationship between Chromatin Accessibility and Gene Expression

We previously conducted RNA sequencing (RNA-seq) of ovary and testis samples [34]. These samples overlapped with those in the current study, allowing us to investigate the underlying relationships between chromatin accessibility and gene expression. We divided genes into five groups based on expression levels (Figure 5A). The Fragments Per Kilobase of exon per Million mapped fragments (FPKM) values obtained from RNA-seq were used as proxies of expression levels. Genes with higher FPKM values showed stronger ATAC signals near the TSS in both the ovary and testis groups (Figure 5A). RNA-seq showed that 11,404 genes had up-regulated expression and 15,523 genes had down-regulated expression in ovaries compared to testes (p ≤ 0.05, log2foldchange > 1). We constructed a quadrant diagram to show the possible relationship between chromatin accessibility and gene expression. Red indicates positive correlations, whereas blue indicates negative correlations (Figure 5B). In ovaries compared to testes, 4745 up-regulated genes and 4339 down-regulated genes were only associated with more accessible regions, while 1087 up-regulated genes and 2558 down-regulated genes were only associated with less accessible regions. Additionally, 4072 up-regulated and 5740 down-regulated genes were simultaneously associated with increased and decreased chromatin accessibility regions (Figure 5C). We listed overlaps between DAGs and differential expression genes (DEGs) in Table S5. We revealed that increased and decreased chromatin accessibility regions existed simultaneously near well-known genes participating in sex differentiation, such as foxl2a, amh, and sox9a.

2.6. Summary of Bisulfite Sequencing Data

To reveal DNA methylation modifications in common carp gonads, we performed BS-seq and generated 370.75 million raw reads. After quality filtering, 363.98 million high-quality clean reads were retained. The bisulfite conversion efficiency for all libraries was over 99.71%. More than 60.56% of clean reads were mapped uniquely to the reference genome, and these reads were submitted for the next analyses. The summary of BS-seq data is shown in Table S6. Of all genomic cytosines, the average percentages of methylated cytosines (mCs) were 3.97% and 4.07% in ovaries and testes, respectively. The methylation rates in CG contexts were 43.32% in ovaries and 44.54% in testes, while the methylation rates in CHG and CHH contexts (H represents T/A/C) were less than 0.2% in both ovaries and testes. Among all the mCs, more than 96.04% in ovaries and 96.26% in testes were mCGs (Table S7). Methylation occurred predominantly in cytosines of CG contexts.

2.7. DNA Methylation Profiles of Ovaries and Testes

We calculated the methylation levels of mCGs. The results showed that most mCGs displayed high methylation levels. Especially in testes, more than 90% of the mCGs demonstrated a methylation level of over 80%. The proportion of mCGs with methylation levels over 60% in the ovaries exceeded 85% (Figure 6A). We further investigated methylation levels in different genome elements. The 5′ UTR regions showed low methylation levels, whereas the repeats and 3′ UTR regions exhibited higher methylation levels overall in both ovaries and testes (Figure 6B). The methylation levels decreased with the proximity to the TSS in 2 kb upstream regions and increased dramatically after the TSS, and the gene body regions exhibited relatively high methylation levels (Figure 6C). We subsequently analyzed differentially methylated regions (DMRs) between ovaries and testes in CG contexts (Table S8). The number of DMRs in different genome regions was calculated. Plenty of DMRs were in gene bodies (exons and introns) and promoter regions, as shown in Figure 6D. A Circos plot was conducted to display the distribution and significance of DMRs. The DMRs were distributed relatively evenly throughout the genome (Figure 6E). There were 4392 DMRs in subgenome A and 5151 DMRs in subgenome B.
Those genes overlapping with DMRs in gene body regions or promoters were regarded as differentially methylated genes (DMGs). We finally identified 2087 promoter-DMGs and 5264 genebody-DMGs, respectively (Table S9). Among them, 1294 genes simultaneously exhibited differentially methylated genes in the promoters and gene bodies. The KEGG enrichment analysis revealed that promoter-DMGs were enriched in metabolic pathways, pentose and glucuronate interconversions, and terms related to cytochrome P450 (Figure 7A), and genebody-DMGs were enriched in the MAPK signaling pathway, Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway (Figure 7B). Genebody-DMGs seem to be enriched in the signaling pathways related to sex differentiation. We also performed GO enrichment analysis, as shown in Figure S1.
In addition, we identified nine promoter-DMGs and 95 genebody-DMGs in CHG contexts (Table S10) and 62 promoter-DMGs and 287 genebody-DMGs in CHH contexts (Table S11). Venn diagrams were made to show overlapping DMGs across different contexts (Figure S2).

2.8. Overall Relationship between DNA Methylation and Gene Expression

We classified the genes into four groups. Genes with FPKM < 1 were considered unexpressed, and the remaining genes were divided into low- (1 ≤ FPKM < 3.76 in ovaries and 1 ≤ FPKM < 2.99 in testes), medium- (3.76 ≤ FPKM ≤ 19.08 in ovaries and 2.99 ≤ FPKM ≤ 14.20 in testes), and high (FPKM > 19.08 in ovaries and FPKM > 14.20 in testes)-expression groups based on mRNA expression levels. We then compared DNA methylation levels in CG contexts among the different groups. In both ovaries and testes, all four groups showed decreased methylation levels from the 2 kb upstream of TSS to TSS, where promoters are located. The unexpressed group showed the relatively highest methylation levels among the four groups in the promoter regions. Notably, the levels of gene expression and DNA methylation showed a strictly negative correlation in the TSS upstream and TES downstream regions in the testis group (Figure 8A).
According to the methylation levels in promoter regions, we ranked promoter-methylated genes and then equally divided them into five groups. The methylation level gates of each group were 0.17, 0.42, 0.63, and 0.76 in the ovaries and 0.16, 0.47, 0.76, and 0.91 in the testes. From the first group to the fifth group, methylation levels increased sequentially. In ovaries, less than 40% of genes were unexpressed in the first and second groups, whereas over 60% were unexpressed in the fourth and fifth groups. Similar results were observed in the testis group (Figure 8B). These results are consistent with the generally acknowledged hypothesis that promoter methylation levels are inversely correlated with transcript abundance.

2.9. Combined Analysis of DMGs and DEGs

A scatter plot was made to display the association between DMRs and DEGs. Red dots represent negative correlations, and blue dots represent positive correlations (Figure 8C). Finally, we identified 441 (431 + 10) promoter-hypermethylated genes whose expression levels were concurrently down-regulated and 175 (170 + 5) promoter-hypomethylated genes whose expression levels were up-regulated (Figure 8D; Table S12). This result included many well-known genes, such as dmrtb1-like, spag6, and fels. Methylation in gene body regions showed more complex correlations in regulating gene expression. Overlaps between genebody-DMGs and DEGs were analyzed, and there were 1659 (837 + 130 + 74 + 618) genes whose methylation levels in gene body regions were negatively related to mRNA expression and 1671 (570 + 74 + 130 + 897) genes whose methylation levels in gene body regions were positively related to mRNA expression (Figure 8E; Table S13).

3. Discussion

The common carp exhibits evident sexual dimorphism in growth, which serves as a crucial catalyst for elucidating the mechanism of sex maintenance and consequently producing mono-sex populations. Despite the extensive efforts made by scholars to uncover the distinct transcripts between ovaries and testes, there is still much room for studying the differences in epigenetic modifications.

3.1. Chromatin Accessibility Profiles in Common Carp Gonads

Chromatin accessibility reflects the regulatory capacity of gene expression, which is involved in various biological processes. Chromatin accessibility varies considerably between ovaries and testes [35,36]. In the common carp, our ATAC-seq data analysis identified 120,701 peaks in ovaries and 193,469 peaks in testes, respectively. A total of 34.32% peaks in ovaries and 22.13% peaks in testes were enriched in promoter regions. The proportion was significantly higher than the proportion of promoter regions in the whole genome, revealing the presence of regulatory elements in these regions [37]. The distribution of peaks was consistent with many previous ATAC-seq data and indicated the reliability of our data [38,39].
Motif analysis showed that more TFs enriched in ovaries were associated with embryonic development and cell proliferation, which might be due to the high transcriptional activity of numerous maternal genes during oogenesis [40]. The most enriched TFs in testes were systemic and related to fundamental biological processes. A previous study has found that Foxl2 and SF1 can bind to the promoter region of cyp19a1a, which is essential for regulating estrogen synthesis [41]. We found that motifs associated with Foxl2 and SF1 were enriched only in the ovary group, indicating that target genes of Foxl2 and SF1 were more accessible in ovaries compared to testes. Many genes were linked with multiple peaks, indicating that several regulatory elements could synergistically regulate a single gene. Integrating the analysis of ATAC-seq and RNA-seq revealed that some genes are associated with both relatively increased and decreased peaks simultaneously and that open chromatin regions could either improve or repress gene expression, which contributed to forming complex regulatory networks. Just like some well-known genes related to sex differentiation, dmrt1 can stimulate the expression of sox9a and repress the expression of foxl2, whereas foxl2 can inhibit dmrt1 and sox9a [3,4].

3.2. Differential DNA Methylation between Ovaries and Testes

DNA methylation is one of the most well-understood epigenetic modifications and is crucial in various biological processes in teleost. We decoded the methylome of common carp gonads using bisulfite sequencing, presenting that over 96.04% of mCs were in CG contexts. The general methylation pattern in different genome elements was consistent with that of other fish [42,43]. For instance, methylation levels presented relatively lower degrees and decreased with proximity to TSS in the promoter regions, while the regions of repeats and 3′ UTR exhibited higher methylation levels overall.
We identified 2087 promoter-DMGs and 5264 genebody-DMGs between ovaries and testes in CG contexts. The promoter-DMGs were enriched in KEGG pathways related to cytochrome P450, which is essential for steroid hormone synthesis [44]. Several studies indicated that methylation in the gene body regions may correlate positively with gene expression or participate in alternative mRNA splicing, exhibiting a more complex mechanism in regulating gene expression [45,46]. The role of genebody methylation in sex differentiation is rarely reported. Our study showed that genebody-DMGs were enriched in the Wnt signaling pathway and GnRH signaling pathway. The Wnt signaling pathway is known to have a crucial role in ovarian differentiation and maintenance [47], and the GnRH signaling pathway is essential in modulating reproduction [48]. These findings suggest that DNA methylation, whether in promoter or genebody regions, might play a vital role in sex maintenance.
Previous studies have illustrated that methylation in the promoter regions could reversely modulate gene expression to govern sex differentiation [14]. Integrated DNA methylation and gene expression analysis revealed genes whose expressions were negatively correlated with the methylation levels in promoter regions. dmrt1b-like may be the mediated regulation of gene expression in this way and be involved in sex maintenance.

4. Materials and Methods

4.1. Sample Collection

The animals used in this experiment were reared at the Fangshan Experimental Station in Beijing, China. The common carp received MS222 (40 mg/L) anesthesia to alleviate their discomfort during sampling. The sexually mature common carp (18 months old) were dissected in October to obtain the gonad, which contains gametes in different development stages. The ovary and testis were subsequently frozen in liquid nitrogen for further use. The ovaries and testis surveyed in the current study had three replicates. All experiments were conducted in compliance with the principles of animal care and use for scientific purposes established by the Animal Care and Use Committee of the Chinese Academy of Fishery Sciences (protocol code ACUC-CAFS-20220615, approved on 15 June 2022).

4.2. ATAC Library Preparation and Sequencing

The library preparation and sequencing were conducted following previously published protocols [49]. The nucleus was extracted from digested gonad tissue and then incubated in the Tn5 transposase reaction mix. The transposase preferentially fragmented DNA in open chromatin regions, and adapter sequences were added simultaneously to the ends of the fragments. Then, PCR amplification was performed to complete the library construction. The library was sequenced on the Illumina Novaseq platform to generate 150 bp paired-end reads.

4.3. ATAC-Seq Reads Alignment

Adaptor sequences and low-quality bases were removed to prepare the data for downstream analysis. The clean reads were mapped to the reference genome using BWA (v0.7.12) with default parameters. The common carp reference genome and gene annotation files were downloaded from our previous study [31]. Reads aligned to mitochondrial DNA and PCR duplicates were discarded. High-quality reads (Mapping Quality ≥ 13) were used for further analysis.

4.4. Peak Calling and Motif Analysis

We used MACS2 (v2.1.0) to perform peak calling [50]. The dynamic Poisson distribution was employed to calculate the p-value of specific regions based on the uniquely mapped reads. Afterward, the Q-value (adjusted p-value) was derived from the p-value using the Benjamini–Hochberg method, and the region would be deemed as a peak when the Q-value is less than 0.05. Peaks were adjusted to the same size with summits centered, and then loci sequences were examined for motif discovery using HOMER software (v4.9.1) with default parameters [51].

4.5. Differential Peaks Analysis and Annotation

All alignment files were scaled to read coverage files normalized by RPM using deepTools (v3.0.2) to identify differential peaks between ovaries and testes. The differential peaks were those with a fold change of RPM greater than two (|log2foldchange| > 1). We retrieved gene annotation information surrounding peaks and analyzed the genomic distributions of peaks using ChIPseeker [52]. In ovaries compared to testes, the regions that enriched more/less peaks were more/less accessible regions.

4.6. Bisulfite Library Construction and Sequencing

Genomic DNA was extracted from gonads using the QIAamp Fast DNA Tissue Kit (Qiagen, Duesseldorf, Germany). The extracted high-quality DNA was mixed with positive control lambda DNA and then fragmented into 200–400 bp fragments using Covaris S220 (Covaris, Woburn, MA, USA). These fragments were treated with Bisulfite using an EZ DNA Methylation-Gold™ Kit (Zymo Research, Tustin, CA, USA). After bisulfite treatment, unmethylated cytosine was converted to uracil, while methylated cytosine remained unchanged. Afterward, we completed library construction using the Accel-NGS Methyl-Seq DNA Library Kit (Swift Biosciences, Ann Arbor, MI, USA). The library was then sequenced on the Novaseq platform, and 150 bp paired-end reads were generated.

4.7. Bisulfite Sequencing Reads Mapping

The raw data were processed to eliminate low-quality reads, adapters, and poly N sequences. The resulting clean reads were used for subsequent analysis. The reference genome and clean reads were first transformed into bisulfite-converted versions (C-to-T and G-to-A converted). The similarly converted versions of the reference genome and clean reads were aligned using Bismark software (v0.16.3) [53]. The unique and best alignment reads were then aligned back to the original reference genome to infer the methylation state of each cytosine position. The percentage of cytosines sequenced at cytosine reference positions in the lambda genome was calculated as the sodium bisulfite non-conversion rate.

4.8. Methylation Level Estimation and Differential Analysis

We conducted a binomial test for every cytosine site to identify methylated sites, and sites with Q-values (Benjamini–Hochberg method corrected p-values) less than 0.05 were thought to be methylated. The methylation level of each cytosine site was calculated as the number of reads methylated at the site divided by the total number of reads covering it, and the methylation level of a region was quantified as the average methylation level of all cytosines in this region. Differentially methylated regions between ovaries and testes were identified using DSS software (v2.12.0) [54,55].

4.9. Gene Functional Annotation

Enrichment analysis was conducted to aid in understanding gene functions in various biological processes using the ClusterProfiler R package (v3.9). This package included GO and KEGG pathway terms. Q-values less than 0.05 were deemed significant.

5. Conclusions

We generated chromatin accessibility landscapes and DNA methylomes of the common carp gonads and identified widespread alterations of epigenetic modifications between ovaries and testes. Combined with the previously published transcriptomes, we further explored the possible epigenetic-mediated regulation of gene expression. This study not only provides novel insights into the potential molecular mechanism of sex differentiation but also facilitates sex control breeding programs for the common carp.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25010321/s1.

Author Contributions

Conceptualization, Y.Z. and J.L.; Data curation, Q.W.; Formal analysis, M.H., K.W. and Y.C. (Yingjie Chen); Funding acquisition, M.H. and J.L.; Investigation, R.Z. and Y.C. (Yiming Cao); Methodology, M.H., R.Z., Y.C. (Yiming Cao), S.Y. and Y.C. (Yingjie Chen); Project administration, Y.Z. and J.L.; Resources, S.Y. and K.W.; Software, Q.W. and J.Z.; Supervision, J.L.; Validation, M.H., X.S. and Y.Z.; Visualization, Q.W. and J.Z.; Writing—original draft, M.H.; Writing—review and editing, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National High-Level Personnel of Special Support Program, Special Scientific Research Funds for Central Non-profit Institutes, Chinese Academy of Fishery Sciences (grant numbers 2023XT0102, 2023A003, and 2023TD25), Beijing Municipal Natural Science Foundation (grant number 6212033), Fisheries Innovation Team of Beijing Agriculture Innovation Consortium (grant number BAIC07-2023-03), and National Freshwater Genetic Resource Center (grant number FGRC: 18537).

Institutional Review Board Statement

The study was conducted following the principles of animal care and use for scientific purposes established by the Animal Care and Use Committee of the Chinese Academy of Fishery Sciences (protocol code ACUC-CAFS-20220615, approved on 15 June 2022).

Informed Consent Statement

Not applicable.

Data Availability Statement

The ATAC-seq and BS-seq data presented in the current study have been submitted to the National Center for Biotechnology Information (NCBI) BioProject database under the accession number PRJNA1042881.

Acknowledgments

We thank Yingjie Qi and Ning Ding for their help with the experimental materials. We also thank Zhengsheng Zhao for the fish culture.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ATAC-seq, Assay for Transposase Accessible Chromatin sequencing; BS-seq, bisulfite sequencing; TF, transcription factor; promoter-DMGs, promoter differentially methylated genes; genebody-DMGs, gene body differentially methylated genes; RPM, Reads Per Million mapped reads; DAGs, differentially accessible genes; RNA-seq, RNA sequencing; DEGs, differential expression genes; FPKM, Fragments Per Kilobase of exon per Million mapped fragments; DMRs, differentially methylated regions; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.

References

  1. Bachtrog, D.; Mank, J.E.; Peichel, C.L.; Kirkpatrick, M.; Otto, S.P.; Ashman, T.L.; Hahn, M.W.; Kitano, J.; Mayrose, I.; Ming, R.; et al. Sex determination: Why so many ways of doing it? PLoS Biol. 2014, 12, e1001899. [Google Scholar] [CrossRef] [PubMed]
  2. Holleley, C.E.; O‘Meally, D.; Sarre, S.D.; Marshall Graves, J.A.; Ezaz, T.; Matsubara, K.; Azad, B.; Zhang, X.; Georges, A. Sex reversal triggers the rapid transition from genetic to temperature-dependent sex. Nature 2015, 523, 79–82. [Google Scholar] [CrossRef] [PubMed]
  3. Li, M.; Sun, L.; Wang, D. Roles of estrogens in fish sexual plasticity and sex differentiation. Gen. Comp. Endocrinol. 2019, 277, 9–16. [Google Scholar] [CrossRef] [PubMed]
  4. Aharon, D.; Marlow, F.L. Sexual determination in zebrafish. Cell. Mol. Life Sci. 2021, 79, 8. [Google Scholar] [CrossRef] [PubMed]
  5. Rajendiran, P.; Jaafar, F.; Kar, S.; Sudhakumari, C.; Senthilkumaran, B.; Parhar, I.S. Sex determination and differentiation in teleost: Roles of genetics, environment, and brain. Biology 2021, 10, 973. [Google Scholar] [CrossRef] [PubMed]
  6. Martínez, P.; Viñas, A.M.; Sánchez, L.; Díaz, N.; Ribas, L.; Piferrer, F. Genetic architecture of sex determination in fish: Applications to sex ratio control in aquaculture. Front. Genet. 2014, 5, 340. [Google Scholar] [CrossRef] [PubMed]
  7. Budd, A.; Banh, Q.; Domingos, J.; Jerry, D. Sex control in fish: Approaches, challenges and opportunities for aquaculture. J. Mar. Sci. Eng. 2015, 3, 329–355. [Google Scholar] [CrossRef]
  8. Ge, C.; Ye, J.; Weber, C.; Sun, W.; Zhang, H.; Zhou, Y.; Cai, C.; Qian, G.; Capel, B. The histone demethylase KDM6B regulates temperature-dependent sex determination in a turtle species. Science 2018, 360, 645–648. [Google Scholar] [CrossRef]
  9. Klemm, S.L.; Shipony, Z.; Greenleaf, W.J. Chromatin accessibility and the regulatory epigenome. Nat. Rev. Genet. 2019, 20, 207–220. [Google Scholar] [CrossRef]
  10. Allis, C.D.; Jenuwein, T. The molecular hallmarks of epigenetic control. Nat. Rev. Genet. 2016, 17, 487–500. [Google Scholar] [CrossRef]
  11. Tafessu, A.; O’Hara, R.; Martire, S.; Dube, A.L.; Saha, P.; Gant, V.U.; Banaszynski, L.A. H3.3 contributes to chromatin accessibility and transcription factor binding at promoter-proximal regulatory elements in embryonic stem cells. Genome Biol. 2023, 24, 25. [Google Scholar] [CrossRef] [PubMed]
  12. Buenrostro, J.D.; Giresi, P.G.; Zaba, L.C.; Chang, H.Y.; Greenleaf, W.J. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat. Methods 2013, 10, 1213–1218. [Google Scholar] [CrossRef]
  13. Buenrostro, J.D.; Wu, B.; Chang, H.Y.; Greenleaf, W.J. ATAC-seq: A method for assaying chromatin accessibility genome-wide. Curr. Protoc. Mol. Biol. 2015, 109, 1–9. [Google Scholar] [CrossRef] [PubMed]
  14. Ortega-Recalde, O.; Goikoetxea, A.; Hore, T.A.; Todd, E.V.; Gemmell, N.J. The genetics and epigenetics of sex change in fish. Annu. Rev. Anim. Biosci. 2020, 8, 47–69. [Google Scholar] [CrossRef] [PubMed]
  15. Shao, C.; Li, Q.; Chen, S.; Zhang, P.; Lian, J.; Hu, Q.; Sun, B.; Jin, L.; Liu, S.; Wang, Z.; et al. Epigenetic modification and inheritance in sexual reversal of fish. Genome Res. 2014, 24, 604–615. [Google Scholar] [CrossRef] [PubMed]
  16. Zhou, H.; Zhuang, Z.; Sun, Y.; Chen, Q.; Zheng, X.; Liang, Y.; Mahboob, S.; Wang, Q. Changes in DNA methylation during epigenetic-associated sex reversal under low temperature in Takifugu rubripes. PLoS ONE 2019, 14, e0221641. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, Y.; Zhang, S.; Liu, Z.; Zhang, L.; Zhang, W. Epigenetic modifications during sex change repress gonadotropin stimulation of cyp19a1a in a teleost ricefield eel (Monopterus albus). Endocrinology 2013, 154, 2881–2890. [Google Scholar] [CrossRef] [PubMed]
  18. Wu, G.C.; Li, H.W.; Huang, C.H.; Lin, H.J.; Lin, C.J.; Chang, C.F. The testis is a primary factor that contributes to epigenetic modifications in the ovaries of the protandrous Black Porgy, Acanthopagrus schlegelii. Biol. Reprod. 2016, 94, 132. [Google Scholar] [CrossRef]
  19. Wen, A.Y.; You, F.; Sun, P.; Li, J.; Xu, D.D.; Wu, Z.H.; Ma, D.Y.; Zhang, P.J. CpG methylation of dmrt1 and cyp19a promoters in relation to their sexual dimorphic expression in the Japanese flounder Paralichthys olivaceus. J. Fish Biol. 2014, 84, 193–205. [Google Scholar] [CrossRef]
  20. Si, Y.; Ding, Y.; He, F.; Wen, H.; Li, J.; Zhao, J.; Huang, Z. DNA methylation level of cyp19a1a and Foxl2 gene related to their expression patterns and reproduction traits during ovary development stages of Japanese flounder (Paralichthys olivaceus). Gene 2016, 575, 321–330. [Google Scholar] [CrossRef]
  21. Todd, E.V.; Ortega-Recalde, O.; Liu, H.; Lamm, M.S.; Rutherford, K.M.; Cross, H.; Black, M.A.; Kardailsky, O.; Graves, J.A.M.; Hore, T.A.; et al. Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change. Sci. Adv. 2019, 5, eaaw7006. [Google Scholar] [CrossRef] [PubMed]
  22. Shabir, U.; Dar, J.S.; Bhat, A.H.; Nazir, R.; Ameen, F. The hidden world of fish fungal pathogens: Molecular identification and phylogenetic analysis in common carp, Cyprinus carpio. Arch. Microbiol. 2023, 205, 311. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, M.; Chen, L.; Zhou, Z.; Xiao, J.; Chen, B.; Huang, P.; Li, C.; Xue, Y.; Liu, R.; Bai, Y.; et al. Comparative transcriptome analysis of early sexual differentiation in the male and female gonads of common carp (Cyprinus carpio). Aquaculture 2023, 563, 738984. [Google Scholar] [CrossRef]
  24. Feng, X.; Yu, X.; Fu, B.; Wang, X.; Liu, H.; Pang, M.; Tong, J. A high-resolution genetic linkage map and QTL fine mapping for growth-related traits and sex in the Yangtze River common carp (Cyprinus carpio haematopterus). BMC Genom. 2018, 19, 230. [Google Scholar] [CrossRef] [PubMed]
  25. Kocour, M.; Linhart, O.; Gela, D.; Rodina, M. Growth performance of all-female and mixed-sex common carp Cyprinus Carpio L. populations in the central Europe climatic conditions. J. World Aquac. Soc. 2005, 36, 103–113. [Google Scholar] [CrossRef]
  26. Chen, L.; Xu, J.; Sun, X.; Xu, P. Research advances and future perspectives of genomics and genetic improvement in allotetraploid common carp. Rev. Aquac. 2021, 14, 957–978. [Google Scholar] [CrossRef]
  27. Jiang, M.; Jia, S.; Chen, J.; Chen, K.; Ma, W.; Wu, X.; Luo, H.; Li, Y.; Zhu, Z.; Hu, W. Timing of gonadal development and dimorphic expression of sex-related genes in gonads during early sex differentiation in the Yellow River carp. Aquaculture 2020, 518, 734825. [Google Scholar] [CrossRef]
  28. Chen, G.; Huang, J.; Jia, J.; Lou, Q.; Shi, C.; Yasheng, M.; Zhao, Y.; Yuan, Q.; Tang, K.; Liu, X.; et al. The food safety assessment of all-female common carp (Cyprinus carpio) (cyp17a1+/−; XX genotype) generated using genome editing technology. Food Chem. Toxicol. 2023, 181, 114103. [Google Scholar] [CrossRef]
  29. Chen, J.J.; Xia, X.H.; Wang, L.F.; Jia, Y.F.; Nan, P.; Li, L.; Chang, Z.J. Identification and comparison of gonadal transcripts of testis and ovary of adult common carp Cyprinus carpio using suppression subtractive hybridization. Theriogenology 2015, 83, 1416–1427. [Google Scholar] [CrossRef]
  30. Anitha, A.; Gupta, Y.R.; Deepa, S.; Ningappa, M.; Rajanna, K.B.; Senthilkumaran, B. Gonadal transcriptome analysis of the common carp, Cyprinus carpio: Identification of differentially expressed genes and SSRs. Gen. Comp. Endocrinol. 2019, 279, 67–77. [Google Scholar] [CrossRef]
  31. Li, J.T.; Wang, Q.; Huang Yang, M.D.; Li, Q.-S.; Cui, M.S.; Dong, Z.J.; Wang, H.W.; Yu, J.H.; Zhao, Y.J.; Yang, C.R.; et al. Parallel subgenome structure and divergent expression evolution of allo-tetraploid common carp and goldfish. Nat. Genet. 2021, 53, 1493–1503. [Google Scholar] [CrossRef] [PubMed]
  32. Yoshida, H.; Lareau, C.A.; Ramirez, R.N.; Rose, S.A.; Maier, B.; Wroblewska, A.; Desland, F.; Chudnovskiy, A.; Mortha, A.; Dominguez, C.; et al. The cis-regulatory atlas of the mouse immune system. Cell 2019, 176, 897–912.e20. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, P.; Xu, J.; Liu, G.; Chen, L.; Zhou, Z.; Peng, W.; Jiang, Y.; Zhao, Z.; Jia, Z.; Sun, Y.; et al. The allotetraploid origin and asymmetrical genome evolution of the common carp Cyprinus carpio. Nat. Commun. 2019, 10, 4625. [Google Scholar] [CrossRef] [PubMed]
  34. Hou, M.; Wang, Q.; Zhang, J.; Zhao, R.; Cao, Y.; Yu, S.; Wang, K.; Chen, Y.; Ma, Z.; Sun, X.; et al. Differential expression of miRNAs, lncRNAs, and circRNAs between ovaries and testes in common carp (Cyprinus carpio). Cells 2023, 12, 2631. [Google Scholar] [CrossRef] [PubMed]
  35. Wu, X.; Yang, Y.; Zhong, C.; Guo, Y.; Wei, T.; Li, S.; Lin, H.; Liu, X. Integration of ATAC-seq and RNA-seq unravels chromatin accessibility during sex reversal in Orange-Spotted Grouper (Epinephelus coioides). Int. J. Mol. Sci. 2020, 21, 2800. [Google Scholar] [CrossRef] [PubMed]
  36. Hu, Q.; Lian, Z.; Xia, X.; Tian, H.; Li, Z. Integrated chromatin accessibility and DNA methylation analysis to reveal the critical epigenetic modification and regulatory mechanism in gonadal differentiation of the sequentially hermaphroditic fish, Monopterus albus. Biol. Sex Differ. 2022, 13, 73. [Google Scholar] [CrossRef] [PubMed]
  37. He, H.H.; Meyer, C.A.; Shin, H.; Bailey, S.T.; Wei, G.; Wang, Q.; Zhang, Y.; Xu, K.; Ni, M.; Lupien, M.; et al. Nucleosome dynamics define transcriptional enhancers. Nat. Genet. 2010, 42, 343–347. [Google Scholar] [CrossRef] [PubMed]
  38. Liu, G.; Wang, W.; Hu, S.; Wang, X.; Zhang, Y. Inherited DNA methylation primes the establishment of accessible chromatin during genome activation. Genome Res. 2018, 28, 998–1007. [Google Scholar] [CrossRef]
  39. Li, Y.; Liu, Y.; Yang, H.; Zhang, T.; Naruse, K.; Tu, Q. Dynamic transcriptional and chromatin accessibility landscape of medaka embryogenesis. Genome Res. 2020, 30, 924–937. [Google Scholar] [CrossRef]
  40. Xu, P.; Meng, A. Maternal factors regulate embryonic development in zebrafish. Sci. Sin. Vitae 2014, 44, 954–964. [Google Scholar] [CrossRef]
  41. Navarro-Martín, L.; Viñas, J.; Ribas, L.; Díaz, N.; Gutiérrez, A.; Di Croce, L.; Piferrer, F. 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. [Google Scholar] [CrossRef] [PubMed]
  42. Murphy, P.J.; Cairns, B.R. Genome-wide DNA methylation profiling in zebrafish. In The Zebrafish-Genetics, Genomics, and Transcriptomics; Methods in Cell Biology; Academic Press: Salt Lake City, UT, USA, 2016; pp. 345–359. [Google Scholar] [CrossRef]
  43. He, L.; Liang, X.; Wang, Q.; Yang, C.; Li, Y.; Liao, L.; Zhu, Z.; Wang, Y. Genome-wide DNA methylation reveals potential epigenetic mechanism of age-dependent viral susceptibility in grass carp. Immun. Ageing 2022, 19, 28. [Google Scholar] [CrossRef] [PubMed]
  44. Tenugu, S.; Pranoty, A.; Mamta, S.-K.; Senthilkumaran, B. Development and organisation of gonadal steroidogenesis in bony fishes. Aquac. Fish. 2021, 6, 223–246. [Google Scholar] [CrossRef]
  45. Arechederra, M.; Daian, F.; Yim, A.; Bazai, S.K.; Richelme, S.; Dono, R.; Saurin, A.J.; Habermann, B.H.; Maina, F. Hypermethylation of gene body CpG islands predicts high dosage of functional oncogenes in liver cancer. Nat. Commun. 2018, 9, 3164. [Google Scholar] [CrossRef] [PubMed]
  46. Nanavaty, V.; Abrash, E.W.; Hong, C.; Park, S.; Fink, E.E.; Li, Z.; Sweet, T.J.; Bhasin, J.M.; Singuri, S.; Lee, B.H.; et al. DNA methylation regulates alternative polyadenylation via CTCF and the cohesin complex. Mol. Cell 2020, 78, 752–764. [Google Scholar] [CrossRef] [PubMed]
  47. Kossack, M.E.; High, S.K.; Hopton, R.E.; Yan, Y.-l.; Postlethwait, J.H.; Draper, B.W. Female sex development and reproductive duct formation depend on Wnt4a in zebrafish. Genetics 2019, 211, 219–233. [Google Scholar] [CrossRef] [PubMed]
  48. Zohar, Y.; Muñoz-Cueto, J.A.; Elizur, A.; Kah, O. Neuroendocrinology of reproduction in teleost fish. Gen. Comp. Endocrinol. 2010, 165, 438–455. [Google Scholar] [CrossRef]
  49. Corces, M.R.; Trevino, A.E.; Hamilton, E.G.; Greenside, P.G.; Sinnott-Armstrong, N.A.; Vesuna, S.; Satpathy, A.T.; Rubin, A.J.; Montine, K.S.; Wu, B.; et al. An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues. Nat. Methods 2017, 14, 959–962. [Google Scholar] [CrossRef]
  50. Zhang, Y.; Liu, T.; Meyer, C.A.; Eeckhoute, J.; Johnson, D.S.; Bernstein, B.E.; Nusbaum, C.; Myers, R.M.; Brown, M.; Li, W.; et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9, R137. [Google Scholar] [CrossRef]
  51. Heinz, S.; Benner, C.; Spann, N.; Bertolino, E.; Lin, Y.C.; Laslo, P.; Cheng, J.X.; Murre, C.; Singh, H.; Glass, C.K. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 2010, 38, 576–589. [Google Scholar] [CrossRef]
  52. Yu, G.; Wang, L.-G.; He, Q.Y. ChIPseeker: An R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics 2015, 31, 2382–2383. [Google Scholar] [CrossRef]
  53. Krueger, F.; Andrews, S.R. Bismark: A flexible aligner and methylation caller for bisulfite-seq applications. Bioinformatics 2011, 27, 1571–1572. [Google Scholar] [CrossRef]
  54. Wu, H.; Xu, T.; Feng, H.; Chen, L.; Li, B.; Yao, B.; Qin, Z.; Jin, P.; Conneely, K.N. Detection of differentially methylated regions from whole-genome bisulfite sequencing data without replicates. Nucleic Acids Res. 2015, 43, e141. [Google Scholar] [CrossRef]
  55. Park, Y.; Wu, H. Differential methylation analysis for BS-seq data under general experimental design. Bioinformatics 2016, 32, 1446–1453. [Google Scholar] [CrossRef]
Ijms 25 00321 g001
Figure 1. Overview of common carp gonadal chromatin accessibility profiles. (A) Distributions of reads around gene body regions in each sample. TSS, transcription start sites; TES, transcription end sites. (B) The heat map of correlations among different samples. (C) Proportions of peaks distributed in different genome elements.
Figure 1. Overview of common carp gonadal chromatin accessibility profiles. (A) Distributions of reads around gene body regions in each sample. TSS, transcription start sites; TES, transcription end sites. (B) The heat map of correlations among different samples. (C) Proportions of peaks distributed in different genome elements.
Ijms 25 00321 g001
Ijms 25 00321 g002
Figure 2. Identification of differential peaks. (A) The Venn diagram of peaks between ovaries and testes. (B) The Scatter plot of RPM for peaks in ovaries and testes. (C) Differential accessible regions between ovaries and testes.
Figure 2. Identification of differential peaks. (A) The Venn diagram of peaks between ovaries and testes. (B) The Scatter plot of RPM for peaks in ovaries and testes. (C) Differential accessible regions between ovaries and testes.
Ijms 25 00321 g002
Ijms 25 00321 g003
Figure 3. Functional annotation of differentially accessible genes. (A) GO enrichment analysis of DAGs between ovaries and testes. ‘*’ represents Q-values less than 0.05. (B) KEGG enrichment analysis of more accessible genes in the testes.
Figure 3. Functional annotation of differentially accessible genes. (A) GO enrichment analysis of DAGs between ovaries and testes. ‘*’ represents Q-values less than 0.05. (B) KEGG enrichment analysis of more accessible genes in the testes.
Ijms 25 00321 g003
Ijms 25 00321 g004
Figure 4. TF binding motifs enriched in open chromatin regions. (A) The bubble diagram of motifs enriched in each sample. (B) Top ten motifs enriched in ovary-biased peaks. (C) Top ten motifs enriched in testis-biased peaks.
Figure 4. TF binding motifs enriched in open chromatin regions. (A) The bubble diagram of motifs enriched in each sample. (B) Top ten motifs enriched in ovary-biased peaks. (C) Top ten motifs enriched in testis-biased peaks.
Ijms 25 00321 g004
Ijms 25 00321 g005
Figure 5. The relationship between chromatin accessibility and gene expression in common carp gonads. (A) The abundance of ATAC signals in gene bodies and flanking regions among gene groups with different expression levels. (B) The underlying relationship between chromatin accessibility and gene expression. Red dots indicate positive correlations; blue dots indicate negative correlations. (C) The number of overlapping genes between DAGs and differential expression genes (DEGs). Red indicates up-regulated genes; blue indicates down-regulated genes. ‘a’ represents genes only associated with more accessible regions, ‘b’ represents genes only associated with less accessible regions, and ‘c’ represents genes simultaneously associated with increased and decreased chromatin accessibility regions.
Figure 5. The relationship between chromatin accessibility and gene expression in common carp gonads. (A) The abundance of ATAC signals in gene bodies and flanking regions among gene groups with different expression levels. (B) The underlying relationship between chromatin accessibility and gene expression. Red dots indicate positive correlations; blue dots indicate negative correlations. (C) The number of overlapping genes between DAGs and differential expression genes (DEGs). Red indicates up-regulated genes; blue indicates down-regulated genes. ‘a’ represents genes only associated with more accessible regions, ‘b’ represents genes only associated with less accessible regions, and ‘c’ represents genes simultaneously associated with increased and decreased chromatin accessibility regions.
Ijms 25 00321 g005
Ijms 25 00321 g006
Figure 6. Methylation levels of mCGs in different genome regions and DMRs between the two groups. (A) Distribution of methylation levels in mCG contexts. (B) Methylation levels in different genome elements. (C) Methylation levels in gene bodies and flanking regions. (D) Number of DMRs in different genome elements. Red and blue represent hypermethylation and hypomethylation regions in ovaries compared to testes. (E) Distribution of DMRs in the genome. Red and blue dots represent hypermethylation and hypomethylation in ovaries compared to testes, respectively. CGI, CpG island; TSS, transcription start sites; TES, transcription end sites; TE, transposable elements.
Figure 6. Methylation levels of mCGs in different genome regions and DMRs between the two groups. (A) Distribution of methylation levels in mCG contexts. (B) Methylation levels in different genome elements. (C) Methylation levels in gene bodies and flanking regions. (D) Number of DMRs in different genome elements. Red and blue represent hypermethylation and hypomethylation regions in ovaries compared to testes. (E) Distribution of DMRs in the genome. Red and blue dots represent hypermethylation and hypomethylation in ovaries compared to testes, respectively. CGI, CpG island; TSS, transcription start sites; TES, transcription end sites; TE, transposable elements.
Ijms 25 00321 g006
Ijms 25 00321 g007
Figure 7. KEGG enrichment analyses of DMGs between ovaries and testes. (A) KEGG enrichment of promoter-DMGs. (B) KEGG enrichment of genebody-DMGs.
Figure 7. KEGG enrichment analyses of DMGs between ovaries and testes. (A) KEGG enrichment of promoter-DMGs. (B) KEGG enrichment of genebody-DMGs.
Ijms 25 00321 g007
Ijms 25 00321 g008
Figure 8. The relationship between DNA methylation and gene expression. (A) The methylation levels among gene groups with different expression levels. (B) The gene expression levels among groups with different methylation levels in promoter regions. Methylation levels increased sequentially from the first group to the fifth group. (C) Association analysis between DMRs and DEGs. The horizontal axis represents the methylation levels of DMRs, and the vertical axis represents gene expression levels (log2foldchange). Red dots represent negative correlations, and blue dots represent positive correlations. (D) The Venn diagram represents the overlaps of promoter-DMGs and DEGs. (E) The Venn diagram represents the overlaps of genebody-DMGs and DEGs.
Figure 8. The relationship between DNA methylation and gene expression. (A) The methylation levels among gene groups with different expression levels. (B) The gene expression levels among groups with different methylation levels in promoter regions. Methylation levels increased sequentially from the first group to the fifth group. (C) Association analysis between DMRs and DEGs. The horizontal axis represents the methylation levels of DMRs, and the vertical axis represents gene expression levels (log2foldchange). Red dots represent negative correlations, and blue dots represent positive correlations. (D) The Venn diagram represents the overlaps of promoter-DMGs and DEGs. (E) The Venn diagram represents the overlaps of genebody-DMGs and DEGs.
Ijms 25 00321 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hou, M.; Wang, Q.; Zhao, R.; Cao, Y.; Zhang, J.; Sun, X.; Yu, S.; Wang, K.; Chen, Y.; Zhang, Y.; et al. Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads. Int. J. Mol. Sci. 2024, 25, 321. https://doi.org/10.3390/ijms25010321

AMA Style

Hou M, Wang Q, Zhao R, Cao Y, Zhang J, Sun X, Yu S, Wang K, Chen Y, Zhang Y, et al. Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads. International Journal of Molecular Sciences. 2024; 25(1):321. https://doi.org/10.3390/ijms25010321

Chicago/Turabian Style

Hou, Mingxi, Qi Wang, Ran Zhao, Yiming Cao, Jin Zhang, Xiaoqing Sun, Shuangting Yu, Kaikuo Wang, Yingjie Chen, Yan Zhang, and et al. 2024. "Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads" International Journal of Molecular Sciences 25, no. 1: 321. https://doi.org/10.3390/ijms25010321

APA Style

Hou, M., Wang, Q., Zhao, R., Cao, Y., Zhang, J., Sun, X., Yu, S., Wang, K., Chen, Y., Zhang, Y., & Li, J. (2024). Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads. International Journal of Molecular Sciences, 25(1), 321. https://doi.org/10.3390/ijms25010321

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop