Next Article in Journal
Evaluation of Gamma Aminobutyric Acid and Sodium Butyrate in Juvenile Red Seabream (Pagrus major) Diets Containing Graded Levels of Fish Meal and Soy Protein Concentrate
Previous Article in Journal
Updates and Current Challenges in Reproductive Microbiome: A Comparative Analysis between Cows and Women
Previous Article in Special Issue
Long Non-Coding RNAs Differentially Expressed in Swine Fetuses
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 14
Issue 13
10.3390/ani14131972
animals-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

The Whole Genome DNA Methylation Signatures of Hindlimb Muscles in Chinese Alligators during Hibernation and Active Periods

by
Jihui Zhang
1,2 and
Xiaobing Wu
2,*
1
School of Food Science and Biology Engineering, Wuhu Institute of Technology, Wuhu 241000, China
2
College of Life Sciences, Anhui Normal University, Wuhu 241000, China
*
Author to whom correspondence should be addressed.
Animals 2024, 14(13), 1972; https://doi.org/10.3390/ani14131972
Submission received: 11 May 2024 / Revised: 12 June 2024 / Accepted: 2 July 2024 / Published: 3 July 2024
(This article belongs to the Special Issue Epigenetics in Animal Development)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Hibernation is a physiological state for Chinese alligators to cope with cold weather. DNA methylation, an epigenetic factor discovered early and extensively studied, plays a crucial role in regulating gene expression, development, and stress responses. The role of DNA methylation in regulating gene expression related to hibernation in ectotherms remains unclear. Here, genome-wide DNA methylation maps of the hindlimb muscles of Chinese alligators were generated, and differential methylation genes and related metabolic pathways of the hindlimb muscles of Chinese alligators during hibernation and active periods were screened. These differential genes and metabolic pathways are closely related to lipid metabolism, energy metabolism, and amino acid metabolism. The findings suggest that DNA methylation may play a significant role in regulating the adaptive hibernation mechanisms of Chinese alligators. This study is helpful to further understand the molecular basis of hibernation adaptation in Chinese alligators.

Abstract

Many ectotherms hibernate to increase their chances of survival during harsh winter conditions. The role of DNA methylation in regulating gene expression related to hibernation in ectotherms remains unclear. Here, we employed whole-genome bisulfite sequencing (WGBS) technology to construct a comprehensive genome-wide DNA methylation landscape of the hindlimb muscles in the Chinese alligator during hibernation and active periods. The results indicated that methylation modifications were most abundant at CG sites, identifying 9447 differentially methylated regions (DMRs) and 2329 differentially methylated genes (DMGs). KEGG pathway enrichment analysis of the DMGs revealed significant enrichment in major pathways such as the neurotrophin signaling pathway, the MAPK signaling pathway, the GnRH signaling pathway, the biosynthesis of amino acids, and the regulation of the actin cytoskeleton, which are closely related to lipid metabolism, energy metabolism, and amino acid metabolism. Among these, 412 differentially methylated genes were located in promoter regions, including genes related to energy metabolism such as ATP5F1C, ATP5MD, PDK3, ANGPTL1, and ANGPTL2, and genes related to ubiquitin-proteasome degradation such as FBXO28, FBXO43, KLHL40, and PSMD5. These findings suggest that methylation in promoter regions may play a significant role in regulating the adaptive hibernation mechanisms in the Chinese alligator. This study contributes to a further understanding of the epigenetic mechanisms behind the hibernation of the Chinese alligator.

1. Introduction

Chinese alligators (Alligator sinensis), a species endemic to China, were once on the brink of extinction. They predominantly inhabit regions between 30.6° and 31.6° north latitude and 118° and 119.6° east longitude. As poikilotherms, Chinese alligators experience a decrease in core body temperature when ambient temperatures fall below 10 °C. Consequently, they cease feeding and enter a hibernation state, with oxygen consumption for respiration reduced to 20–25% of normal activity levels [1]. Prior researchers have demonstrated that the quality of hibernation directly impacts the hibernation survival rate and subsequent year’s reproductive rate of Chinese alligators [2,3]. Our team conducted a detailed comparative study on the blood of Chinese alligators during hibernation and the active phase and found that the red blood cells, white blood cells, cholesterol, and triglycerides in the blood during hibernation were significantly reduced, and the mean corpuscular volume and the mean corpuscular hemoglobin were significantly increased [4]. So far, the seasonal adaptation of the Chinese alligator and its intestinal flora have been previously studied by transcriptome and metagenomics analysis [1,2,5]. Building on this foundation, we hope to further investigate the regulatory role of methylation on differentially expressed genes in the muscles of Chinese alligators during their hibernation and active phases. Today, in addition to human hunting and habitat disruption, global climate change is becoming a critical threat to Chinese alligators, as indeed it is for many other hibernators [6]. Therefore, investigating the gene regulatory network underlying seasonal physiological changes is of paramount importance not only in elucidating how hibernating ectotherms adapt to harsh winter conditions and limited food resources in their natural habitats but also potentially contributing to their future conservation efforts.
DNA methylation, first discovered in 1975 [7,8], remains one of the most extensively investigated and well-studied epigenetic factors. DNA methylation is an ancient epigenetic modification in eukaryotic genomes that plays essential roles in various biological processes, including the regulation of gene expression, development, and stress responses [9,10,11,12,13]. Previous research evaluated global changes in DNA methylation in response to hibernation in the liver and skeletal muscle of thirteen-lined ground squirrels, along with changes in the expression of DNA methyltransferases (DNMT1/3B) and methyl binding domain proteins (MBDs) [14]. In the brown adipose tissue of hibernating thirteen-lined ground squirrels, global transcriptional suppression is achieved through epigenetic regulatory mechanisms [15]. However, previous studies largely focused on changes in overall DNA methylation levels or the methylation of specific genes. Epigenetic factors, including miRNA and DNA methylation, play crucial roles in regulating gene expression and serve as pivotal regulators of skeletal muscle development, energy metabolism, and various other biological processes [16,17,18,19]. Previously, we investigated the regulatory mechanism of miRNA in the hindlimb muscles of adult Chinese alligators with regards to their adaptability during hibernation [20]. In this study, we employed whole-genome bisulfite sequencing (WGBS) to generate single-base resolution DNA methylation maps of the hindlimb muscles in hibernating and active Chinese alligators, aiming to elucidate the methylation changes in the hindlimb muscles of Chinese alligators during hibernation and activity phases.
The vertebrate skeletal muscle is a substantial tissue in the body, primarily responsible for bodily movement, energy metabolism, protein storage, safeguarding internal organs, and fulfilling other vital functions [21,22]. The weight of the hindlimb muscles of quadrupedal tetrapod species accounts for a large proportion of body weight. Compared with the hindlimbs of other vertebrate groups, the hindlimbs of alligators have stronger compliant muscles, which enhances their mobility [23,24,25]. Alligators are unusual among tetrapods in their ability to use a wide range of hindlimb postures, even over a restricted range of speeds [26,27]. Research conducted on the skeletal muscles of hibernating animals demonstrates that their skeletal muscles exhibit well-maintained strength and anti-fatigue capabilities throughout the period of hibernation [28,29,30]. A similar phenomenon is observed in the skeletal muscles of hibernating crocodiles, which enables them to rapidly resume active behaviors such as predator avoidance, foraging, and mating when they awaken in the spring [31]. In addition, the muscles of amphibians also show adaptability to hibernation. For example, their skeletal muscles have been shown to degrade to produce glycogen to meet their hibernation needs [32]. Hibernation is characterized by prolonged periods of inactivity with concomitantly low nutrient intake, conditions that would typically result in muscle atrophy combined with a loss of oxidative fibers. Yet, hibernators consistently emerge from winter with very little atrophy. The investigation of the unique adaptability of skeletal muscles during reptile hibernation holds significant value for academic research, thus warranting further exploration. Although skeletal muscle adaptations have been extensively investigated in hibernating species such as Urocitellus richardsonii [28], Ursus arctos [33], and Spermophilus dauricus [34], the study of hindlimb muscle adaptations during hibernation in reptiles remains limited to a few species [2,35,36]. The findings from this study will provide a fundamental foundation for gaining further insights into the adaptive regulatory mechanisms of hindlimb muscles in Chinese alligators during hibernation.

2. Materials and Methods

2.1. Sample Collection

As the Chinese alligator is a rare protected species endemic to China, the availability of specimens for scientific research is extremely limited. We obtained four adult female Chinese alligators from the Chinese Alligator Nature Reserve in Xuancheng City, Anhui Province, China. These adult Chinese alligators exhibited similar weight and length characteristics (Supplementary Table S1A). Two female specimens were collected in January 2018 during deep hibernation, with an average body temperature of approximately 8.9 °C. The other two female specimens were collected in June 2018 during the active phase, with an average body temperature of ca. 23.7 °C. Euthanasia of the four adult Chinese alligators was conducted via intravenous injection of pentobarbital (350 mg/kg). Muscle tissue from the hindlimb (the quadriceps of the hindlimb) was dissected for analysis. The collected samples (Supplementary Table S1B) were promptly frozen in liquid nitrogen and stored at −80 °C until further use. All animal experiments were conducted with approval from the animal ethics committee of Anhui Normal University, under a license number of 2018012.

2.2. DNA Extraction, Library Construction, and Sequencing

Genomic DNA was extracted from hindlimb muscle tissue samples of Chinese alligators during both hibernation and active periods using a genomic DNA extraction kit (TIANamp, Beijing, China). Subsequently, the purity and concentration of the DNA samples were measured using 1% agarose gel electrophoresis and the Qubit 2.0 fluorometer (Life Technologies, Carlsbad, CA, USA). We utilized 1 µg of high-quality genomic DNA and supplemented it with 1 ng of unmethylated lambda DNA as a negative control. The genomic DNA underwent random shearing to yield fragments ranging from 200 to 400 bp using Covaris M220, followed by repair and adenylation of the fragmented DNA ends. Adapters containing methylated cytosines were then ligated to the A-tailed DNA fragments. Subsequently, the DNA fragments, extracted from gel according to the EZ DNA Methylation Gold Kit (Zymo Research, Irvine, CA, USA) instructions, underwent bisulfite treatment. Following treatment, methylated cytosines on the genomic DNA remained unaltered, while unmethylated cytosines were converted to uracils (subsequently transformed to thymines after PCR amplification). Finally, the DNA library was enriched through PCR amplification. The amplified library was purified by 0.85× AMPure XP. After purification, Qubit 2.0 was used for preliminary quantification, and the library was diluted to 1 ng/µL. Then the fragment range and effective concentration of the library were measured by Agilent 2100 (Agilent, Waldbroon, Germany) and qPCR to ensure the quality of the library, which can be used for subsequent experiments. After the paired-end library construction, the libraries were sequenced on an Illumina HiSeq 2500 platform (Illumina, San Diego, CA, USA).

2.3. Sequencing Data Analysis

The raw sequencing data were filtered using the fastqc software (v0.11.9, https://www.bioinformatics.babraham.ac.uk/projects/fastqc/, accessed on 20 January 2021) to remove low-quality junk data, retain clean data, and ensure its quality. We trimmed the sequencing adapters and low-quality data from the raw sequencing data to obtain clean data for subsequent analysis [37]. Then, the clean reads were aligned to the Chinese alligator genome [38] using Bismap software (V2.90) with the default parameters. Conversions of C to T and G to A were separately made in the reference genome and sequence reads before alignment, and the bisulfite-converted genome was indexed by Bowtie2 (v2.4.0). The bisulfite non-conversion rate is determined by the proportion of cytosines sequenced at the cytosine reference site. With coverage of ≥ 5× and a false discovery rate of < 0.05, we performed a binomial distribution test for each C site to detect methylated cytosines [39,40].

2.4. Identification of DMRs and DMGs

We employed the Metilene software (v0.2.8) for analyzing differentially methylated regions (DMRs). This software employs a binary segmentation algorithm in conjunction with dual statistical tests (MWU-test and 2D KS-test) to effectively detect DMRs between paired samples or two groups of samples. Subsequently, DMRs are refined through multiple testing corrections. In this study, CpG sites were utilized to identify differentially methylated regions. The criteria for selecting differentially methylated regions were as follows: an average methylation level difference > 0.1, a false discovery rate (FDR) < 0.05 for differential testing between the two regions, and the presence of at least five methylated CpG sites within each DMR, with the distance between CpG sites being less than 300 bp. Based on the DMR annotation results, genes overlapping with DMRs in the genebody and its upstream and downstream regions (upstream 2 K, genebody, and downstream 2 K) were screened and defined as differentially methylated genes (DMGs). We defined the upstream 2000 bp of genes as the promoter region and used the same screening method to identify differentially promoter-methylated genes (DPMGs).

2.5. Gene Ontology and KEGG Pathway Analysis of DMGs and DPMGs

The DMGs and DPMGs were screened and annotated by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. GO enrichment analysis was achieved using the GOseq R software (v3.19) package [41], in which GO items with a p-value < 0.05 were considered significantly enriched. KEGG can be used to analyze advanced functions and biological systems (such as cells, organisms, and ecosystems) at a molecular level (http://www.genome.jp/kegg/, accessed on 20 January 2021) [42]. KOBAS software (v3.0) was used to analyze the statistical enrichment of DMGs and DPMGs in KEGG pathways [43].

2.6. Bisulfite Sequencing PCR

Gene-specific DNA methylation was assessed by a next-generation sequencing-based BSP, according to a previously published method [44,45]. In brief, BSP primers were designed using the online MethPrimer software (v2.0) and listed in Supplementary Table S11. A total of 1 μg of genomic DNA was converted using the ZYMO EZ DNA Methylation-Gold Kit (Zymo Research, Irvine, CA, USA), and one twentieth of the elution products were used as templates for PCR amplification with 35 cycles using the KAPA 2G Robust HotStart PCR Kit (Kapa Biosystems, Wilmington, MA, USA). For each sample, BSP products of multiple genes were pooled equally, 5′-phosphorylated, 3′-dA-tailed, and ligated to the barcoded adapter using T4 DNA ligase (NEB). Barcoded libraries from all samples were sequenced on the Illumina platform.
For the bisulfite sequencing reads of each sample, adapters and low-quality reads were removed using the software Trimmomatic-0.36. After removing the adapter sequences and filtering out the low-quality reads, the clean sequencing reads were directly aligned to the target sequences using the software Bsmap (v2.73) with the default parameters, which combine genome hashing and bitwise masking to achieve fast and accurate bisulfite mapping. Methylation levels are defined as the fraction of read counts of ‘C’ in the total read counts of both ‘C’ and ‘T’ for each covered C site. On the basis of such a read fraction, methylated cytosine was shown using a binomial distribution as in the method described by Lister et al. [46], whereby a probability mass function is calculated for each methylation context (CpG, CHG, CHH). A two-tailed Fisher’s exact test was used to identify cytosines that are differentially methylated between two samples or groups. Only those CGs covered by at least 200 reads in at least one sample were considered for testing. p-value thresholds were selected such that the significance level was less than 0.001.

2.7. Supporting Data Information

The raw genome-wide methylation data have been submitted to the NCBI Sequence Read Archive under accession number PRJNA1088816.

3. Results

3.1. Summary of Methylome Sequencing

Through WGBS, we obtained the raw data of four hindlimb muscle samples from the hibernation and active periods, with each period’s samples comprising two female Chinese alligators. After quality control, more than 198 million reads were obtained from all samples, with base numbers greater than 4.61 × 1010 bp and Q20 greater than 97.93% (Supplementary Table S2), which was sufficient for subsequent analysis. Bisulfite conversion efficiency (BCE) reached 99%, suggesting the reliability of the methylome sequencing in this study.

3.2. DNA Methylation Patterns

We conducted coverage statistics for C sites, calculating the coverage of C sites in three sequence contexts: CG, CHH, and CHG (i.e., the number of reads supporting that context). The coverage of C sites at 5× was between 0.092909376 and 0.397458967, while at 10×, it ranged from 0.038056356 to 0.127042279. Please refer to Table 1 for specific statistical results. We observed an average methylation rate of 2.97% across all genomic C sites in the four samples. In vertebrates, DNA methylation predominantly occurs within three sequence contexts: CpG, CHG (where H represents A, C, or T), and CHH. The results of this study showed that the methylation ratio of CG sequence types in the hindlimb muscle of a Chinese alligator was 71.12–74.87%. The methylation ratio of CHH sequence type ranges from 0.09–0.84%, and that of CHG sequence type ranges from 0.09–0.83% (Supplementary Table S3). The proportion of methylated C sites (C sites that have at least one methylated read) to all methylated C sites in mCG, mCHG, and mCHH sequences was further analyzed. The results showed that the proportion of CG-methylated C sites was 76.98–96.85%, the proportion of mCHG was 0.75–5.27%, and the proportion of mCHH was 2.40–17.76% (Figure 1), indicating that the methylation of C sites in the mCG sequence context is the predominant type of genomic methylation in the hindlimb muscles of the Chinese alligator.

3.3. DNA Methylation Levels in Different Regions of the Gene

To further understand the methylation levels in different regions of the genome, this study targeted three sequence environments (mCG, mCHG, and mCHH) for C sites in the genebody and transcription start sites (TSS) upstream of 2 K. Average methylation levels at 2 K downstream of transcriptional termination sites (TES) were statistically analyzed. As shown in Figure 2, methylation levels of exons and introns were higher in CG sequences than in promoter regions upstream of transcription start sites (TSS) at 2 K, whereas in CHG and CHH sequences, exons and introns have lower methylation levels than the promoter region at 2 K upstream of the transcription start site (TSS). In the context of CHG and CHH sequences, the DNA methylation level in the 2 K region upstream of the transcription start site (TSS) was higher, and the DNA methylation level in the exon region was lower.

3.4. Differentially Methylated Regions and Differentially Methylated Gene Identification

In this study, CG sites were used to find differential methylation regions. A total of 9447 differentially methylated regions (DMRs) were detected in the hindlimb muscles of the Chinese alligator during the hibernation and active phase (HM vs. AM) (Supplementary Table S4), corresponding to 2329 genes (Supplementary Table S5). The average length of DMRs is 276, with 88.20% of DMRs being below 600 bp (Figure 3a). According to the average methylation level of the differential methylation region, the cluster heat map of DMRs was plotted (Figure 3b), and there were significant differences in methylation levels at different times. For the DMR anchor region, the DMR number distribution statistics were carried out, and most DMRs were located in the intron region, followed by the promoter region (Supplementary Table S6).

3.5. GO Annotation and KEGG Pathway Enrichment on the DMGs and DPMGs

In order to delve deeper into the potential links between differentially methylated genes (DMGs) and differentially promoter-methylated genes (DPMGs) and the hibernation adaptability of the Chinese alligator, we performed functional enrichment analysis using the clusterprofiler package. Our findings revealed a significant enrichment of DMGs across 79 GO terms (p < 0.05; Supplementary Table S7). From these, the most enriched terms included ATP binding, actin binding, calcium ion binding, actin cytoskeleton organization, intracellular signal transduction, olfactory receptor activity, signal transduction, presynaptic active zone, cell junction, regulation of GTPase activity, collagen trimer, extracellular matrix structural constituent, and G protein-coupled receptor activity (Figure 4a).
According to the KEGG analysis, the DMGs were significantly enriched in 36 pathways (p < 0.05; Supplementary Table S8). The results showed that pathway enrichment items mainly included the neurotrophin signaling pathway, adhesive patch, ErbB signaling pathway, Hippo signaling pathway, MAPK signaling pathway, oxytocin signaling pathway, GnRH signaling pathway, amino acid biosynthesis, actin cytoskeleton regulation, etc. (Figure 4b).
To further investigate the potential association of DPMGs with hibernation in the Chinese alligator, 426 DMRs corresponding to 412 genes were screened in the promoter regions of the genome. With regard to DPMGs, 106 GO terms were found to be significantly enriched (p < 0.05; Supplementary Table S9), including proton-transporting ATP synthase complex, fatty-acyl-CoA binding, glycerol-3-phosphate O-acyltransferase activity, hydroxyl-glutamine hydrolase activity, pyruvate dehydrogenase (acetyl-transferring) kinase activity, ATP synthesis-coupled proton transport, and other items closely related to energy metabolism. In addition, methyltransferase activity is also one of the significantly enriched GO items (Figure 5a). According to GO enrichment entries and gene functional annotations, genes related to energy metabolism were identified as ATP5F1C, ATP5MD, PDK3, ANGPTL1, and ANGPTL2. Genes related to ubiquitination and proteasome degradation were FBXO28, FBXO43, KLHL40, and PSMD5. These results suggest that promoter methylation may play an important role in the regulation of hibernation adaptability in Chinese alligators.
A total of 92 pathways were enriched (p < 0.05; Supplementary Table S10), in which the most significant pathways included graft-versus-host disease, allograft rejection, autoimmune thyroid disease, HTLV-I infection, influenza A, legionellosis, antigen processing and presentation, and phagosome (Figure 5b). These pathways are mostly related to immunity, which may be caused by the difference in disease resistance between hibernating and active Chinese alligators. In some major enrichment pathways, such as thyroid hormone synthesis, apoptosis, glutathione metabolism, and the Notch signaling pathway, there was no significant enrichment. But most of these pathways are related to immune activity, and these pathways may be related to the onset of disease and immune changes during hibernation in alligators.

3.6. Bisulfite Sequencing PCR (BSP) Verification Results

Six DMRs (three hypo-methylated DMRs and three hyper-methylated DMRs) were randomly selected from the HM group to verify the reliability of the WGBS data through BSP, namely, DMR_NW_005841875 (ACSL6), DMR_NW_005842223 (MICAL2), DMR_NW_005842074 (NDUFV2), DMR_NW_005842035 (PSMD1), DMR_NW_005841904 (RYR3), and DMR_NW_005841865 (STK4). Compared with the methylation levels observed in the AM group, the HM group exhibited hyper-methylation of ACSL6, MICAL2, and RYR3 (p < 0.05). Conversely, the HM group showed hypo-methylation of NDUFV2, PSMD1, and STK4 compared to the AM group (p < 0.05) (Figure 6). The methylation analysis of these six DMRs was consistent between BSP and WGBS, indicating that the WGBS data in this study were reliable.

4. Discussion

DNA methylation has been reported to play a significant role in regulating gene expression associated with mammalian hibernation [14,15,47]. In thirteen-lined ground squirrels, the global DNA methylation level was increased in brown adipose tissue [15] but decreased in skeletal muscle [14]. The CpG methylation level of the MEF2C promoter region was correlated with the downregulation of gene expression in the skeletal muscle of thirteen-lined ground squirrels [14]. In chipmunks, CpG methylation in the USF-binding site is crucial for liver-specific transcription of the hibernation-specific gene, HP-27 [47]. Although these studies provided a glimpse into DNA regulation in hibernation, they mainly focused on changes in overall DNA methylation levels or the methylation of specific genes.
After analyzing the transcriptome data of hindlimb muscle during hibernation and active periods in Chinese alligators, we observed that the muscle tissue exhibited the highest number of differentially expressed genes [20]. Furthermore, a significant proportion of these differentially expressed genes identified through GO enrichment analysis were found to be associated with the hibernation adaptability of Chinese alligators. Simultaneously, pathway analysis revealed significantly enriched pathways encompassing fatty acid metabolism, starch and sucrose metabolism, ketone body synthesis and degradation, the thyroid hormone signaling pathway, the pentose phosphate pathway, and fatty acid degradation in the muscles of Chinese alligators during hibernation [20]. These findings indicate a multitude of energy regulation processes occurring within their musculature. The differential gene expression in the muscle tissue of Chinese alligators during hibernation plays a critical role in coping with energy deficiency during this period. DNA methylation, as an important epigenetic factor regulating differential gene expression, is speculated to modulate the expression of key genes during physiological changes associated with hibernation. Therefore, in order to understand the association between methylation in the hindlimb muscles of Chinese alligators during hibernation and the physiological changes during this period, genome-wide DNA methylation maps with single-base resolution were generated.
The distribution of methylation sites in functional regions revealed that within the CG sequence context, the genebody region exhibited the highest level of DNA methylation, while the promoter region located at 2 K upstream of the transcriptional start site (TSS) displayed a relatively lower level of DNA methylation. Notably, the promoter proximal region near the TSS exhibited the lowest level of DNA methylation. These findings were consistent with observations in other species [48,49]. This observation suggests that different species exhibit similar DNA methylation patterns and possess certain conserved characteristics. A total of 2329 genes displayed differential methylation in the muscle comparison group during hibernation and active periods, with 412 of these genes exhibiting differential methylation in the promoter region, indicating that the majority of methylation differences occurred within the genebody region. This result is consistent with previous reports indicating widespread methylation within the genebody region [50,51]. The absence of a promoter region may be attributed to the hypomethylation status of CpG sites within the promoter region.
To explore the role of methylation modification in the muscle genome of the Chinese alligator hindlimb, pathway enrichment analysis was conducted on 2329 DMR-related genes. The main pathways of significant enrichment include the neurotrophin signaling pathway, adhesion plaque, ErbB signaling pathway, Hippo signaling pathway, MAPK signaling pathway, oxytocin signaling pathway, amino acid biosynthesis, and regulation of the actin cytoskeleton. The neurotrophin signaling pathway is a cellular signaling mechanism that plays a crucial role in the development, plasticity, and repair of the nervous system [52]. The ErbB signaling pathway regulates cell proliferation, migration, differentiation, apoptosis, and cell motility by mediating the PI3K/Akt signaling pathway, the JAK/STAT signaling pathway, and the MAPK signaling pathway [53]. The Hippo signaling pathway, composed of a set of conserved kinases, is a signaling pathway that inhibits cell growth [54]. The MAPK signaling pathway is primarily involved in the regulation of cell proliferation, growth, and differentiation in response to activation by various growth factors, cytokines, mitogens, and hormone receptors [55]. These pathways may be involved in protecting the brain, reducing metabolism, and controlling growth to conserve energy during hibernation in Chinese alligators.
It has been reported that the methylation modification of the gene promoter region is closely related to gene expression level and complex, and the correlation between them can be positive or negative [56,57,58]. Through GO enrichment entries and gene function annotation, some genes related to the differential methylation region in the promoter region were screened out that were related to energy metabolism and protein degradation. For example, PDK3, pyruvate dehydrogenase kinase 3, through phosphorylation, inhibits pyruvate dehydrogenase activity, thus regulating glucose metabolism and aerobic respiration, causing PDHA2 phosphorylation, reducing glucose utilization, and, at the same time, increasing fat metabolism in response to a long period of fasting, low nutrition level, and starvation state, and playing a role in maintaining glucose homeostasis and normal blood sugar levels [59,60,61]. ATP5F1C and ATP5MD are mitochondrial membrane ATP synthases, which catalyze ADP to produce ATP in the presence of a proton gradient on the cell membrane. ATP5MD is a small subunit of mitochondrial membrane ATP synthase that is required by the dimerization of ATP synthase and that regulates ATP synthesis in mitochondria [62,63]. Both FBXO28 and FBXO43 are members of the F-box protein family that recognize and bind some phosphorylated proteins and promote their ubiquitination and degradation [64]. PSMD5, the 26S subunit of the proteasome, acts as a chaperone protein during the assembly of the 26S proteasome and is directly involved in the ubiquitination proteasome degradation pathway [65,66]. As a key regulatory factor of skeletal muscle development and a gene related to the ubiquitination pathway of protein degradation, KLHL40 is differentially methylated in the promoter region [67], which may play a key role in the synthesis and degradation of hibernating Chinese alligator muscles. In conclusion, it is postulated that differential methylation of the promoter region plays a pivotal role in facilitating hibernation adaptability among Chinese alligators. Significantly enriched pathways in pathway analysis are mainly related to body immunity, which may be associated with the drastic changes in the immune defense capabilities of Chinese alligators during hibernation and active periods. This suggests that the decreased immunity of Chinese alligators during hibernation may exacerbate the survival pressure during this period.

5. Conclusions

This study employed whole-genome bisulfite sequencing (WGBS) technology to construct comprehensive DNA methylation maps for the hindlimb muscles of the Chinese alligator during both hibernation and active periods. By comparing the methylation patterns between these two states, 9447 differentially methylated regions (DMRs) and 2329 differentially methylated genes (DMGs) were identified, highlighting extensive epigenetic differences between hibernation and non-hibernation states that are likely related to the alligator’s physiological adjustments to cold environments. The KEGG enrichment analysis of DMGs revealed significant enrichment in several biological pathways, including energy metabolism, amino acid metabolism, and neurotrophin signaling pathways, suggesting that methylation may play a role in energy utilization, substance metabolism, and signal transduction during hibernation in Chinese alligators. This research provides valuable data and insights for understanding the epigenetic mechanisms behind the hibernation of the Chinese alligator.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani14131972/s1, Table S1A: Basic physiological information of the sample; Table S1B: Sample details and numbers; Table S2: Statistics and analysis of filtered data; Table S3: List of whole genome methylations; Table S4: Statistics of DMRs in the HM vs. AM comparison group; Table S5: Differentially methylated regions and differentially methylated genes; Table S6: DMR distribution among different genetic elements; Table S7: GO enrichment of DMR-related genes; Table S8: KEGG enrichment with DMR-related genes; Table S9: GO enrichment with DMR-related genes of the promoter region; Table S10: KEGG enrichment with DMR-related genes of the promoter region. Table S11: Bisulfite sequencing PCR primer information.

Author Contributions

Conceptualization, J.Z.; methodology, J.Z.; software, J.Z.; validation, J.Z.; formal analysis, J.Z.; investigation, J.Z.; resources, J.Z.; data curation, J.Z.; writing—original draft preparation, J.Z.; writing—review and editing, J.Z.; visualization, J.Z.; supervision, X.W.; project administration, X.W.; funding acquisition, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (No. 32170525), the Anhui University scientific research project (KJ2021A1323), and the Wuhu Vocational and Technical College outstanding youth top talent project (rc2022qnbj01).

Institutional Review Board Statement

The animal study protocol was approved by the Anhui Normal University Academic Ethics Committee (No. 2018012).

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and analyzed during the current study are available from the first author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lin, J.Q.; Huang, Y.Y.; Bian, M.Y.; Wan, Q.H.; Fang, S.G. A Unique Energy-Saving Strategy during Hibernation Revealed by Multi-Omics Analysis in the Chinese Alligator. iScience 2020, 23, 101202. [Google Scholar] [CrossRef] [PubMed]
  2. Sun, H.; Zuo, X.; Sun, L.; Yan, P.; Zhang, F.; Xue, H.; Li, E.; Zhou, Y.; Wu, R.; Wu, X. Insights into the seasonal adaptive mechanisms of Chinese alligators (Alligator sinensis) from transcriptomic analyses. Aust. J. Zool. 2018, 66, 93–102. [Google Scholar] [CrossRef]
  3. Zhang, R.; Hu, Y.; Wang, H.; Yan, P.; Zhou, Y.; Wu, R.; Wu, X. Molecular cloning, characterization, tissue distribution and mRNA expression changes during the hibernation and reproductive periods of estrogen receptor alpha (ESR1) in Chinese alligator, Alligator sinensis. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 2016, 200, 28–35. [Google Scholar] [CrossRef] [PubMed]
  4. Peng, F.; Chen, X.; Meng, T.; Li, E.; Zhou, Y.; Zhang, S. Hematology and serum biochemistry parameters of captive Chinese alligators (Alligator sinensis) during the active and hibernating periods. Tissue Cell 2018, 51, 8–13. [Google Scholar] [CrossRef] [PubMed]
  5. Tang, K.Y.; Wang, Z.W.; Wan, Q.H.; Fang, S.G. Metagenomics Reveals Seasonal Functional Adaptation of the Gut Microbiome to Host Feeding and Fasting in the Chinese Alligator. Front. Microbiol. 2019, 10, 2409. [Google Scholar] [CrossRef] [PubMed]
  6. Humphries, M.M.; Thomas, D.W.; Speakman, J.R. Climate-mediated energetic constraints on the distribution of hibernating mammals. Nature 2002, 418, 313–316. [Google Scholar] [CrossRef] [PubMed]
  7. Holliday, R.; Pugh, J.E. DNA modification mechanisms and gene activity during development. Science 1975, 187, 226–232. [Google Scholar] [CrossRef] [PubMed]
  8. Riggs, A.D. X inactivation, differentiation, and DNA methylation. Cytogenet. Cell Genet. 1975, 14, 9–25. [Google Scholar] [CrossRef]
  9. Breiling, A.; Lyko, F. Epigenetic regulatory functions of DNA modifications: 5-methylcytosine and beyond. Epigenetics Chromatin 2015, 8, 24. [Google Scholar] [CrossRef]
  10. Pelizzola, M.; Ecker, J.R. The DNA methylome. FEBS Lett. 2011, 585, 1994–2000. [Google Scholar] [CrossRef]
  11. Smith, Z.D.; Meissner, A. DNA methylation: Roles in mammalian development. Nat. Rev. Genet. 2013, 14, 204–220. [Google Scholar] [CrossRef]
  12. Su, Z.; Han, L.; Zhao, Z. Conservation and divergence of DNA methylation in eukaryotes: New insights from single base-resolution DNA methylomes. Epigenetics 2011, 6, 134–140. [Google Scholar] [CrossRef] [PubMed]
  13. Ilango, S.; Paital, B.; Jayachandran, P.; Padma, P.R.; Nirmaladevi, R. Epigenetic alterations in cancer. Front. Biosci. 2020, 25, 1058–1109. [Google Scholar] [CrossRef]
  14. Alvarado, S.; Mak, T.; Liu, S.; Storey, K.B.; Szyf, M. Dynamic changes in global and gene-specific DNA methylation during hibernation in adult thirteen-lined ground squirrels, Ictidomys tridecemlineatus. J. Exp. Biol. 2015, 218, 1787–1795. [Google Scholar] [CrossRef]
  15. Biggar, Y.; Storey, K.B. Global DNA modifications suppress transcription in brown adipose tissue during hibernation. Cryobiology 2014, 69, 333–338. [Google Scholar] [CrossRef]
  16. Guller, I.; Russell, A.P. MicroRNAs in skeletal muscle: Their role and regulation in development, disease and function. J. Physiol. 2010, 588, 4075–4087. [Google Scholar] [CrossRef]
  17. Svensson, K.; Handschin, C. MicroRNAs emerge as modulators of NAD+-dependent energy metabolism in skeletal muscle. Diabetes 2014, 63, 1451–1453. [Google Scholar] [CrossRef] [PubMed]
  18. Iswariya, G.T.; Paital, B.; Padma, P.R.; Nirmaladevi, R. microRNAs: Epigenetic players in cancer and aging. Front. Biosci. 2019, 11, 29–55. [Google Scholar] [CrossRef]
  19. Subaramaniyam, U.; Ramalingam, D.; Balan, R.; Paital, B.; Sar, P.; Ramalingam, N. Annonaceous acetogenins as promising DNA methylation inhibitors to prevent and treat leukemogenesis—An in silico approach. J. Biomol. Struct. Dyn. 2023, 1–14. [Google Scholar] [CrossRef]
  20. Zhang, J.; Cai, R.; Liang, J.; Izaz, A.; Shu, Y.; Pan, T.; Wu, X. Molecular mechanism of Chinese alligator (Alligator sinensis) adapting to hibernation. J. Exp. Zool. B Mol. Dev. Evol. 2021, 336, 32–49. [Google Scholar] [CrossRef]
  21. Egan, B.; Zierath, J.R. Exercise metabolism and the molecular regulation of skeletal muscle adaptation. Cell Metab. 2013, 17, 162–184. [Google Scholar] [CrossRef]
  22. Horak, M.; Novak, J.; Bienertova-Vasku, J. Muscle-specific microRNAs in skeletal muscle development. Dev. Biol. 2016, 410, 1–13. [Google Scholar] [CrossRef]
  23. Reilly, S.M.; Blob, R.W. Motor control of locomotor hindlimb posture in the American alligator (Alligator mississippiensis). J. Exp. Biol. 2003, 206, 4327–4340. [Google Scholar] [CrossRef]
  24. Reilly, S.M.; Willey, J.S.; Biknevicius, A.R.; Blob, R.W. Hindlimb function in the alligator: Integrating movements, motor patterns, ground reaction forces and bone strain of terrestrial locomotion. J. Exp. Biol. 2005, 208, 993–1009. [Google Scholar] [CrossRef]
  25. Iijima, M.; Munteanu, V.D.; Elsey, R.M.; Blob, R.W. Correction: Ontogenetic changes in limb posture, kinematics, forces and joint moments in American alligators (Alligator mississippiensis). J. Exp. Biol. 2022, 225, jeb242990. [Google Scholar] [CrossRef]
  26. Reilly, S.M. Sprawling locomotion in the lizard Sceloporus clarkii: Speed modulation of motor patterns in a walking trot. Brain Behav. Evol. 1998, 52, 126–138. [Google Scholar] [CrossRef]
  27. Blob, R.W.; Biewener, A.A. In vivo locomotor strain in the hindlimb bones of alligator mississippiensis and iguana iguana: Implications for the evolution of limb bone safety factor and non-sprawling limb posture. J. Exp. Biol. 1999, 202 Pt 9, 1023–1046. [Google Scholar] [CrossRef]
  28. Green, S.R.; Storey, K.B. Functional and post-translational characterization of pyruvate dehydrogenase demonstrates repression of activity in the liver but not skeletal muscle of the Richardson’s ground squirrel (Urocitellus richardsonii) during hibernation. J. Therm. Biol. 2021, 99, 102996. [Google Scholar] [CrossRef] [PubMed]
  29. Zhang, J.; Li, X.; Ismail, F.; Xu, S.; Wang, Z.; Peng, X.; Yang, C.; Chang, H.; Wang, H.; Gao, Y. Priority Strategy of Intracellular Ca(2+) Homeostasis in Skeletal Muscle Fibers During the Multiple Stresses of Hibernation. Cells 2019, 9, 42. [Google Scholar] [CrossRef] [PubMed]
  30. Andres-Mateos, E.; Mejias, R.; Soleimani, A.; Lin, B.M.; Burks, T.N.; Marx, R.; Lin, B.; Zellars, R.C.; Zhang, Y.; Huso, D.L.; et al. Impaired skeletal muscle regeneration in the absence of fibrosis during hibernation in 13-lined ground squirrels. PLoS ONE 2012, 7, e48884. [Google Scholar] [CrossRef] [PubMed]
  31. Cotton, C.J. Skeletal muscle mass and composition during mammalian hibernation. J. Exp. Biol. 2016, 219, 226–234. [Google Scholar] [CrossRef] [PubMed]
  32. Costanzo, J.P. Overwintering adaptations and extreme freeze tolerance in a subarctic population of the wood frog, Rana sylvatica. J. Comp. Physiol. B 2019, 189, 1–15. [Google Scholar] [CrossRef] [PubMed]
  33. Vella, C.A.; Nelson, O.L.; Jansen, H.T.; Robbins, C.T.; Jensen, A.E.; Constantinescu, S.; Abbott, M.J.; Turcotte, L.P. Regulation of metabolism during hibernation in brown bears (Ursus arctos): Involvement of cortisol, PGC-1alpha and AMPK in adipose tissue and skeletal muscle. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 2020, 240, 110591. [Google Scholar] [CrossRef] [PubMed]
  34. Yang, C.X.; He, Y.; Gao, Y.F.; Wang, H.P.; Goswami, N. Changes in calpains and calpastatin in the soleus muscle of Daurian ground squirrels during hibernation. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 2014, 176, 26–31. [Google Scholar] [CrossRef]
  35. Capraro, A.; O’Meally, D.; Waters, S.A.; Patel, H.R.; Georges, A.; Waters, P.D. Waking the sleeping dragon: Gene expression profiling reveals adaptive strategies of the hibernating reptile Pogona vitticeps. BMC Genom. 2019, 20, 460. [Google Scholar] [CrossRef]
  36. Capraro, A.; O’Meally, D.; Waters, S.A.; Patel, H.R.; Georges, A.; Waters, P.D. MicroRNA dynamics during hibernation of the Australian central bearded dragon (Pogona vitticeps). Sci. Rep. 2020, 10, 17854. [Google Scholar] [CrossRef]
  37. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  38. Wan, Q.H.; Pan, S.K.; Hu, L.; Zhu, Y.; Xu, P.W.; Xia, J.Q.; Chen, H.; He, G.Y.; He, J.; Ni, X.W.; et al. Genome analysis and signature discovery for diving and sensory properties of the endangered Chinese alligator. Cell Res. 2013, 23, 1091–1105. [Google Scholar] [CrossRef]
  39. Gifford, C.A.; Ziller, M.J.; Gu, H.; Trapnell, C.; Donaghey, J.; Tsankov, A.; Shalek, A.K.; Kelley, D.R.; Shishkin, A.A.; Issner, R.; et al. Transcriptional and epigenetic dynamics during specification of human embryonic stem cells. Cell 2013, 153, 1149–1163. [Google Scholar] [CrossRef]
  40. Habibi, E.; Brinkman, A.B.; Arand, J.; Kroeze, L.I.; Kerstens, H.H.; Matarese, F.; Lepikhov, K.; Gut, M.; Brun-Heath, I.; Hubner, N.C.; et al. Whole-genome bisulfite sequencing of two distinct interconvertible DNA methylomes of mouse embryonic stem cells. Cell Stem Cell 2013, 13, 360–369. [Google Scholar] [CrossRef]
  41. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef] [PubMed]
  42. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T.; et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36, D480–D484. [Google Scholar] [CrossRef] [PubMed]
  43. Mao, X.; Cai, T.; Olyarchuk, J.G.; Wei, L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef] [PubMed]
  44. Gao, F.; Liang, H.; Lu, H.; Wang, J.; Xia, M.; Yuan, Z.; Yao, Y.; Wang, T.; Tan, X.; Laurence, A.; et al. Global analysis of DNA methylation in hepatocellular carcinoma by a liquid hybridization capture-based bisulfite sequencing approach. Clin. Epigenetics 2015, 7, 86. [Google Scholar] [CrossRef] [PubMed]
  45. Pan, X.; Gong, D.; Nguyen, D.N.; Zhang, X.; Hu, Q.; Lu, H.; Fredholm, M.; Sangild, P.T.; Gao, F. Early microbial colonization affects DNA methylation of genes related to intestinal immunity and metabolism in preterm pigs. DNA Res. 2018, 25, 287–296. [Google Scholar] [CrossRef] [PubMed]
  46. Lister, R.; Pelizzola, M.; Dowen, R.H.; Hawkins, R.D.; Hon, G.; Tonti-Filippini, J.; Nery, J.R.; Lee, L.; Ye, Z.; Ngo, Q.M.; et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 2009, 462, 315–322. [Google Scholar] [CrossRef] [PubMed]
  47. Fujii, G.; Nakamura, Y.; Tsukamoto, D.; Ito, M.; Shiba, T.; Takamatsu, N. CpG methylation at the USF-binding site is important for the liver-specific transcription of the chipmunk HP-27 gene. Biochem. J. 2006, 395, 203–209. [Google Scholar] [CrossRef] [PubMed]
  48. Fu, Y.; Li, J.; Tang, Q.; Zou, C.; Shen, L.; Jin, L.; Li, C.; Fang, C.; Liu, R.; Li, M.; et al. Integrated analysis of methylome, transcriptome and miRNAome of three pig breeds. Epigenomics 2018, 10, 597–612. [Google Scholar] [CrossRef] [PubMed]
  49. Yang, Y.; Fan, X.; Yan, J.; Chen, M.; Zhu, M.; Tang, Y.; Liu, S.; Tang, Z. A comprehensive epigenome atlas reveals DNA methylation regulating skeletal muscle development. Nucleic Acids Res. 2021, 49, 1313–1329. [Google Scholar] [CrossRef]
  50. Jin, L.; Jiang, Z.; Xia, Y.D. Genome-wide DNA methylation changes in skeletal muscle between young and middle-aged pigs. BMC Genom. 2014, 15, 653. [Google Scholar] [CrossRef]
  51. Yuan, X.L.; Gao, N.; Xing, Y.; Zhang, H.B.; Zhang, A.L.; Liu, J.; He, J.L.; Xu, Y.; Lin, W.M.; Chen, Z.M.; et al. Profiling the genome-wide DNA methylation pattern of porcine ovaries using reduced representation bisulfite sequencing. Sci. Rep. 2016, 6, 22138. [Google Scholar] [CrossRef]
  52. Yano, H.; Chao, M.V. Neurotrophin receptor structure and interactions. Pharm. Acta Helv. 2000, 74, 253–260. [Google Scholar] [CrossRef]
  53. Yarden, Y.; Sliwkowski, M.X. Untangling the ErbB signalling network. Nat. Rev. Mol. Cell Biol. 2001, 2, 127–137. [Google Scholar] [CrossRef]
  54. Chan, S.W.; Lim, C.J.; Chen, L.; Chong, Y.F.; Huang, C.; Song, H.; Hong, W. The Hippo pathway in biological control and cancer development. J. Cell Physiol. 2011, 226, 928–939. [Google Scholar] [CrossRef]
  55. Ajibade, A.A.; Wang, H.Y.; Wang, R.F. Cell type-specific function of TAK1 in innate immune signaling. Trends Immunol. 2013, 34, 307–316. [Google Scholar] [CrossRef]
  56. Hu, Y.; Xu, H.; Li, Z.; Zheng, X.; Jia, X.; Nie, Q.; Zhang, X. Comparison of the genome-wide DNA methylation profiles between fast-growing and slow-growing broilers. PLoS ONE 2013, 8, e56411. [Google Scholar] [CrossRef]
  57. Mohn, F.; Weber, M.; Rebhan, M.; Roloff, T.C.; Richter, J.; Stadler, M.B.; Bibel, M.; Schubeler, D. Lineage-specific polycomb targets and de novo DNA methylation define restriction and potential of neuronal progenitors. Mol. Cell 2008, 30, 755–766. [Google Scholar] [CrossRef]
  58. Payer, B.; Lee, J.T. X chromosome dosage compensation: How mammals keep the balance. Annu. Rev. Genet. 2008, 42, 733–772. [Google Scholar] [CrossRef]
  59. Kato, M.; Chuang, J.L.; Tso, S.C.; Wynn, R.M.; Chuang, D.T. Crystal structure of pyruvate dehydrogenase kinase 3 bound to lipoyl domain 2 of human pyruvate dehydrogenase complex. EMBO J. 2005, 24, 1763–1774. [Google Scholar] [CrossRef] [PubMed]
  60. Korotchkina, L.G.; Patel, M.S. Site specificity of four pyruvate dehydrogenase kinase isoenzymes toward the three phosphorylation sites of human pyruvate dehydrogenase. J. Biol. Chem. 2001, 276, 37223–37229. [Google Scholar] [CrossRef] [PubMed]
  61. Lu, C.W.; Lin, S.C.; Chen, K.F.; Lai, Y.Y.; Tsai, S.J. Induction of pyruvate dehydrogenase kinase-3 by hypoxia-inducible factor-1 promotes metabolic switch and drug resistance. J. Biol. Chem. 2008, 283, 28106–28114. [Google Scholar] [CrossRef]
  62. Barca, E.; Ganetzky, R.D.; Potluri, P.; Juanola-Falgarona, M.; Gai, X.; Li, D.; Jalas, C.; Hirsch, Y.; Emmanuele, V.; Tadesse, S.; et al. USMG5 Ashkenazi Jewish founder mutation impairs mitochondrial complex V dimerization and ATP synthesis. Hum. Mol. Genet. 2018, 27, 3305–3312. [Google Scholar] [CrossRef]
  63. Ohsakaya, S.; Fujikawa, M.; Hisabori, T.; Yoshida, M. Knockdown of DAPIT (diabetes-associated protein in insulin-sensitive tissue) results in loss of ATP synthase in mitochondria. J. Biol. Chem. 2011, 286, 20292–20296. [Google Scholar] [CrossRef]
  64. Jin, J.; Cardozo, T.; Lovering, R.C.; Elledge, S.J.; Pagano, M.; Harper, J.W. Systematic analysis and nomenclature of mammalian F-box proteins. Genes. Dev. 2004, 18, 2573–2580. [Google Scholar] [CrossRef]
  65. Kaneko, T.; Hamazaki, J.; Iemura, S.; Sasaki, K.; Furuyama, K.; Natsume, T.; Tanaka, K.; Murata, S. Assembly pathway of the Mammalian proteasome base subcomplex is mediated by multiple specific chaperones. Cell 2009, 137, 914–925. [Google Scholar] [CrossRef]
  66. Roelofs, J.; Park, S.; Haas, W.; Tian, G.; McAllister, F.E.; Huo, Y.; Lee, B.H.; Zhang, F.; Shi, Y.; Gygi, S.P.; et al. Chaperone-mediated pathway of proteasome regulatory particle assembly. Nature 2009, 459, 861–865. [Google Scholar] [CrossRef]
  67. Ravenscroft, G.; Miyatake, S.; Lehtokari, V.L.; Todd, E.J.; Vornanen, P.; Yau, K.S.; Hayashi, Y.K.; Miyake, N.; Tsurusaki, Y.; Doi, H.; et al. Mutations in KLHL40 are a frequent cause of severe autosomal-recessive nemaline myopathy. Am. J. Hum. Genet. 2013, 93, 6–18. [Google Scholar] [CrossRef]
Animals 14 01972 g001
Figure 1. Proportional distribution of methylation C sites. (a) AM1: Samples of hindlimb muscles from active period individual 1; (b) AM2: Samples of hindlimb muscles from active period individual 2; (c) HM1: Samples of hindlimb muscles from hibernation period individual 1; (d) HM2: Samples of hindlimb muscles from hibernation period individual 2.
Figure 1. Proportional distribution of methylation C sites. (a) AM1: Samples of hindlimb muscles from active period individual 1; (b) AM2: Samples of hindlimb muscles from active period individual 2; (c) HM1: Samples of hindlimb muscles from hibernation period individual 1; (d) HM2: Samples of hindlimb muscles from hibernation period individual 2.
Animals 14 01972 g001
Animals 14 01972 g002
Figure 2. Distribution of sample methylation level in the genebody and upstream and downstream of 2 K. The horizontal axis represents different regions, while the vertical axis represents the level of methylation. Each region of each gene is divided into 20 bins, and then the average level of CpG sites in the corresponding bins of all regions is taken. Different colors represent different sequence contexts (CpG, CHG, and CHH). (a) Distribution of sample AM1 methylation level in the genebody and upstream and downstream of 2 K. (b) Distribution of sample AM2 methylation level in the genebody and upstream and downstream of 2 K. (c) Distribution of sample HM1 methylation level in the genebody and upstream and downstream of 2 K. (d) Distribution of sample HM2 methylation level in the genebody and upstream and downstream of 2 K.
Figure 2. Distribution of sample methylation level in the genebody and upstream and downstream of 2 K. The horizontal axis represents different regions, while the vertical axis represents the level of methylation. Each region of each gene is divided into 20 bins, and then the average level of CpG sites in the corresponding bins of all regions is taken. Different colors represent different sequence contexts (CpG, CHG, and CHH). (a) Distribution of sample AM1 methylation level in the genebody and upstream and downstream of 2 K. (b) Distribution of sample AM2 methylation level in the genebody and upstream and downstream of 2 K. (c) Distribution of sample HM1 methylation level in the genebody and upstream and downstream of 2 K. (d) Distribution of sample HM2 methylation level in the genebody and upstream and downstream of 2 K.
Animals 14 01972 g002
Animals 14 01972 g003
Figure 3. Differential methylation region analysis in the hibernating Chinese alligator. (a) Length distribution of DMRs; number of DMRs: 9447. (b) Heatmap of methylation levels in the hibernation and active hindlimb muscles comparison group; number of DMRs: 5999; HM includes two samples, HM1 and HM2; AM includes two samples, AM1 and AM2.
Figure 3. Differential methylation region analysis in the hibernating Chinese alligator. (a) Length distribution of DMRs; number of DMRs: 9447. (b) Heatmap of methylation levels in the hibernation and active hindlimb muscles comparison group; number of DMRs: 5999; HM includes two samples, HM1 and HM2; AM includes two samples, AM1 and AM2.
Animals 14 01972 g003
Animals 14 01972 g004
Figure 4. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment on the differentially methylated genes; number of genes: 2329. (a) TOP items by GO enrichment of DMR-related genes. (b) TOP items by KEGG enrichment with DMR-related genes.
Figure 4. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment on the differentially methylated genes; number of genes: 2329. (a) TOP items by GO enrichment of DMR-related genes. (b) TOP items by KEGG enrichment with DMR-related genes.
Animals 14 01972 g004
Animals 14 01972 g005
Figure 5. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment on the differentially promoter-methylated genes; number of genes: 412. (a) TOP items by GO enrichment with DMR-related genes of the promoter region. (b) KEGG enrichment of TOP entries with DMR-related genes in the promoter region.
Figure 5. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment on the differentially promoter-methylated genes; number of genes: 412. (a) TOP items by GO enrichment with DMR-related genes of the promoter region. (b) KEGG enrichment of TOP entries with DMR-related genes in the promoter region.
Animals 14 01972 g005
Animals 14 01972 g006
Figure 6. Verification of DMR-associated genes (DMGs) by bisulfite sequencing PCR (BSP). (a) Comparison of the degree of methylation of DMR_NW_005841875 (ACSL6) between BSP and whole-genome bisulfite sequencing (WGBS). For BSP, the ordinate represents the mean methylation rate of two samples in each group, while for WGBS, it represents the methylation level of the mean normalized DMR. (b) Comparison of the extent of methylation of DMR_NW_005842223 (MICAL2) between BSP and WGBS; (c) Comparison of the extent of methylation of DMR_NW_005841904 (RYR3) between BSP and WGBS; (d) Comparison of the extent of methylation of DMR_NW_005842074 (NDUFV2) between BSP and WGBS; (e) Comparison of the extent of methylation of DMR_NW_005842035 (PMSD1) between BSP and WGBS; (f) Comparison of the extent of methylation of DMR_NW_005841865 (STK4) between BSP and WGBS.
Figure 6. Verification of DMR-associated genes (DMGs) by bisulfite sequencing PCR (BSP). (a) Comparison of the degree of methylation of DMR_NW_005841875 (ACSL6) between BSP and whole-genome bisulfite sequencing (WGBS). For BSP, the ordinate represents the mean methylation rate of two samples in each group, while for WGBS, it represents the methylation level of the mean normalized DMR. (b) Comparison of the extent of methylation of DMR_NW_005842223 (MICAL2) between BSP and WGBS; (c) Comparison of the extent of methylation of DMR_NW_005841904 (RYR3) between BSP and WGBS; (d) Comparison of the extent of methylation of DMR_NW_005842074 (NDUFV2) between BSP and WGBS; (e) Comparison of the extent of methylation of DMR_NW_005842035 (PMSD1) between BSP and WGBS; (f) Comparison of the extent of methylation of DMR_NW_005841865 (STK4) between BSP and WGBS.
Animals 14 01972 g006
Table 1. Sequencing depth and coverage statistics of the C sites of different samples.
Table 1. Sequencing depth and coverage statistics of the C sites of different samples.
SampleContextTotal CytosineMean DepthEffect DepthCoverageCoverage (d ≥ 5)Coverage (d ≥ 10)
HM1CG469681044.158986756.0682552630.68536780.2935119540.091199125
CHG2303512874.227674685.6358259080.75014290.3510069510.112980272
CHH7007213714.438444995.7748822250.76857760.3683892890.121295292
HM2CG469698314.202580385.7008280760.73718770.3024494810.088304065
CHG2303576924.416765325.4937828460.8039570.3714694230.114224464
CHH7007380674.698717025.6962170240.82488380.3974589670.127042279
AM1CG469659112.517898035.1740293730.48664160.1621979180.045486097
CHG2303480122.701891854.8056406490.56223340.2091803380.061420725
CHH7007139763.032301925.0211832230.60390190.2370496720.072808576
AM2CG469573601.62910867.7369408480.21056240.0929093760.038056356
CHG2303196221.537019566.5489664630.23469650.1149334770.048410092
CHH7006391101.735148176.7677751480.25638380.1292970610.055583607
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

Zhang, J.; Wu, X. The Whole Genome DNA Methylation Signatures of Hindlimb Muscles in Chinese Alligators during Hibernation and Active Periods. Animals 2024, 14, 1972. https://doi.org/10.3390/ani14131972

AMA Style

Zhang J, Wu X. The Whole Genome DNA Methylation Signatures of Hindlimb Muscles in Chinese Alligators during Hibernation and Active Periods. Animals. 2024; 14(13):1972. https://doi.org/10.3390/ani14131972

Chicago/Turabian Style

Zhang, Jihui, and Xiaobing Wu. 2024. "The Whole Genome DNA Methylation Signatures of Hindlimb Muscles in Chinese Alligators during Hibernation and Active Periods" Animals 14, no. 13: 1972. https://doi.org/10.3390/ani14131972

APA Style

Zhang, J., & Wu, X. (2024). The Whole Genome DNA Methylation Signatures of Hindlimb Muscles in Chinese Alligators during Hibernation and Active Periods. Animals, 14(13), 1972. https://doi.org/10.3390/ani14131972

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