1. Introduction
Epigenetic mechanisms are cellular processes that do not change DNA sequences but regulate gene expression [1]. DNA methylation modification is the most common epigenetic phenomenon and has become one of the most studied hotspots in the field of epigenetic research in recent years. Whole genome methylation sequencing is a high-definition sequencing technology that detects the methylation status of cytosine bases. Studies have shown that DNA methylation has a very important role in aquatic animal growth, development, sex reversal, and environmental adaptation. In aquatic animals, a large number of genome-wide methylation studies have been carried out, including on E. sinensis [2], Daphnia magna [3], Pacific oysters (Crassostrea gigas) [4], Takifugu rubripes [5], and Monopterus albus [6].
Intersex is a biological phenomenon characterized by the simultaneous presence of male and female traits [7]. Intersex may occur for a variety of reasons, including environmental contamination [8], parasitism [9], and genetic abnormalities [10]. These factors may disrupt the development and function of the gonads, as well as alter the homeostasis of sex hormones. In aquatic animals, including amphibians [11], fish [12], and crustaceans [13,14,15], the phenomenon of intersex has been reported. In the half-smooth tongue sole (Cynoglossus semilaevis), temperature affects DNA methylation and regulates sex determination-related genes, leading to masculinization in fish with a female genotype [16]. Environmental stimuli have been found to increase the methylation levels of the aromatase gene in bluehead wrasse (Thalassoma bifasciatum), inhibiting the maintenance of female gonads and ultimately causing sex reversal [17]. Besides DNA methylation modifications, histone demethylation can also trigger sex reversal in temperature-dependent sex determination in red-eared slider turtles (Trachemys scripta elegans) [18]. Thus, epigenetic regulatory mechanisms play a crucial role in sexual plasticity.
E. sinensis is regarded as an important aquatic resource due to its exceptional economic and ecological value. The occurrence of intersex individuals within this species may lead to significant economic losses in the aquaculture industry. However, DNA methylation in the formation of intersex individuals is still poorly understood. We have conducted research on the morphology, histology, sex genes, and transcriptomics of intersex E. sinensis [19,20,21]. In this study, we employed EM-seq technology to systematically analyze the genome-wide DNA methylation profiles of intersex crabs and normal crabs (female and male), thereby expanding the bioinformatics resources of the DNA methylation landscape in E. sinensis. Considering the biological specificity of intersex, studying DNA methylation in intersex crabs could offer new insights into its formation mechanisms.
2. Results
2.1. Sequencing Data Quality Control and Reference Genome Alignment
To assess the DNA methylation profiles in intersex (cm) and normal crabs (fm and mm), we conducted EM-seq. After quality control and sequencing data processing, valid data for all samples were obtained. The Q20 (%) and Q30 (%) were similar across samples. The content of guanine + cytosine (GC) ranged from 21.93% to 22.31% and did not differ significantly for cm, fm, and mm samples. The raw datum for cm was 49.51 Gb, which was lower than that of fm (57.8 Gb) and mm (54.8 Gb). After eliminating the low-quality sequencing data, 38.65 Gb of valid data remained for cm, which was also lower than that of fm (46.28 Gb) and mm (44.02 Gb) (Table 1). After enzyme treatment, the enzyme conversion of cm was 98.35%, which was between fm and mm. After mapping these data to the reference genome, the alignment ratio of cm was 43.33%, which was lower than fm and mm. The average methylation rate of C sites in cm was 3.5%, which was lower than fm (4.5%), and higher than mm (2.6%) (Table 2). The above results indicated that our data were up to standard and can be used for subsequent studies.
2.2. Genome-Wide DNA Methylation Profiles
The cytosine methylation types were classified as CG, CHH, and CHG. H can be any of A, T or C. In intersex crabs, the average methylation rates of CG, CHH, and CHG were 48%, 42.8%, and 9.2%, respectively. In female crabs, the average methylation rates of CG, CHH, and CHG were 39.2%, 49.2%, and 11.6%, respectively. In male crabs, the average methylation rates of CG, CHH, and CHG were 59.5%, 33.6%, and 6.9%, respectively (Figure 1). The C-base methylation levels of different functional elements in the CG, CHG, and CHH sequences were analyzed. The functional elements (from 2 kb upstream of the transcription start site to 2 kb downstream of the transcription termination site) were categorized into promoter, exon, intron, and downstream regions (intergenic region) (Figure 2). The results showed that methylation in the intersex E. sinensis mainly happened in the CG sequence context, which was lower than that in males but higher than in females. However, intersex crabs were lower than that of females but higher than males for CHG and CHH. We found that the methylation levels in the promoter region were higher than those in other functional element regions.
We counted the number of DMRs in the cm vs. fm and cm vs. mm comparison groups. In the analysis of the cm vs. fm group, we identified 394,511 hyper-methylated DMRs and 3,393,915 hypo-methylated DMRs. In the analysis of the cm vs. mm group, we identified 2,553,753 hyper-methylated DMRs and 436,082 hypo-methylated DMRs (Table S1). It also showed that the methylation levels of intersex crabs were lower than those of female crabs and higher than those of male crabs.
2.3. Functional Enrichment of Promoter-Region DMGs
DMGs obtained from cm vs. fm and cm vs. mm groups were analyzed by GO (Figure S1) and KEGG enrichment. In both comparison groups, KEGG enrichment of the top 10 pathways in terms of number of genes was shown. Metabolism pathways, peroxisomes, RNA polymerases, and ribosome biogenesis in eukaryotes were common to the cm vs. fm and cm vs. mm groups. The MAPK signaling pathway and phagosomes were unique to the cm vs. fm group. Endocytosis and phosphatidylinositol signaling systems were unique to the cm vs. mm group (Figure 3).
2.4. PPI Networks of DMGs in the Promoter Region of Intersex and Normal Crabs
In cm vs. fm, 460 DMGs in the promoter region were identified and enriched in specific pathways (Table S2). In the cm vs. mm group, 378 DMGs in the promoter region were found and enriched to specific pathways (Table S3). A PPI network containing 149 nodes (each representing a gene) and 279 edges (interactions between nodes) was constructed by analyzing 460 DMGs in the cm vs. fm group (Figure 4a). By using MCODE to analyze 149 genes, four key clusters of gene networks were identified, which were mainly associated with heat shock proteins, metabolism (fatty acid synthesis/degradation, purine metabolism, glycolysis/gluconeogenesis, fructose/mannose metabolism, amino acid/nucleotide metabolism and pyruvate metabolism), and the ribosomes (Figure 4b). The 378 DMGs screened in the cm vs. mm group were analyzed to form a PPI network containing 110 nodes and 130 edges (Figure 4c). Three key gene network clusters were screened, which were mainly involved in metabolism (fatty acid elongation/degradation, purine metabolism, glycolysis/gluconeogenesis, fructose/mannose metabolism, galactose metabolism, and pentose phosphate metabolism) and the ribosomes (Figure 4d). Among the metabolic pathways enriched, we found fatty acid elongation/degradation, purine metabolism, glycolysis/gluconeogenesis, and fructose/mannose metabolism were common to the cm vs. fm and cm vs. mm groups. However, amino acid/nucleotide metabolism and pyruvate metabolism were specific to the cm vs. fm group while galactose metabolism and pentose phosphate metabolism were specific to the cm vs. mm group.
2.5. Regulation of Gene Expression by Differential Methylation of Promoter Region
To investigate the impact of DNA methylation on the expression levels of key genes, we selected six key genes from the identified networks for validation. In the fm vs. cm group, heat shock cognate protein 71 kDa (Hsc70), heat shock protein 90 (Hsp90), 60S ribosomal protein L18 (Rpl18), and long-chain acyl-CoA synthetase 1 (Acsl1) were selected. In the mm vs. cm group, 60S ribosomal protein L7 (Rpl7) and Yippee-interacting protein 2 (Yip2) were selected. EM-seq data indicated that the methylation levels of the promoter region of these key genes were consistently higher in normal crabs compared with intersex crabs (Figure 5). The methylation levels obtained from BSP and EM-seq displayed a similar trend (Figure 6). Additionally, we employed qRT-PCR to assess the expression profiles of the selected key DMGs (Figure 7). The results demonstrated that the expression levels of these key DMGs were significantly upregulated in intersex crabs.
3. Discussion
DNA methylation is a form of epigenetic regulation with important roles in gene expression and tissue development [22]. As is known, arthropods exhibit low global genomic cytosine methylation levels, with some variability existing between crustaceans and insects (0–1% of total cytosines) [23]. However, the methylation levels of swimming crab (Portunus trituberculatus) range from 1 to 1.5% [24]. In this study, the genome-wide methylation levels of E. sinensis ranging from 2.6% to 4.5% were even higher than those of the swimming crab. DNA methylation levels in vertebrates are 5–10% [25]. The genome-wide methylation levels of E. sinensis were lower than those of vertebrates. Studies have shown that differences in the ability of species to adapt to their environments and the levels of evolution of species are important causes of differences in DNA methylation levels between species [26].
To our knowledge, this study represents the first comprehensive comparison of genome-wide DNA methylation patterns between normal and intersex crabs. There are differences among the three forms of cytosine methylation, with mCG being the predominant methylation type, which is consistent with the previous studies in other species [27,28,29,30]. By counting the number of DMRs and genome-wide methylation levels, we found that females are more methylated than intersex, and intersex are more methylated than males. The differences in DNA methylation levels between samples suggest that the intersex crab is an intermediate state between female and male crabs, which also corresponds to the sex marker identification as female. This difference may be the result of environmental stress and plays an important role in the formation of intersex crabs [31].
Heat shock proteins (HSPs) are a class of functionally similar proteins that alter gene expression in response to environmental stresses [32,33]. This family of genes is involved in the regulation of sex determination in several temperature-dependent sex-determining species [34,35,36]. After heat treatment, Hsc70 overexpression in seminiferous cells of half-smooth tongue sole can affect the expression of the sex-related genes sox9a and cyp19a1a. Hsc70 acts as a key regulator linking external high-temperature signaling to in vivo sex differentiation [37]. Similarly, Hsc70 is a hub gene for pseudomale half-smooth tongue sole formation [38]. Hsc70 and Hsp90 can interact with steroid hormone receptors [39,40]. In our study, promoter regions of Hsc70 and Hsp90 show hypermethylation in female and hypomethylation in intersex crabs. Their expression levels is an opposite trend. These findings suggest that Hsc70 and Hsp90 genes may be regulated by DNA methylation. DNA methylation modifications of Hsc70 and Hsp90 may be involved in the changes in female crabs into intersex crabs by regulating the expression of these genes in response to environmental factors.
KEGG enrichment and PPI network analysis revealed that many DMGs were enriched in the ribosomal pathway between the normal and intersex crabs. A number of ribosomal protein genes as key candidate genes were identified. Studies have shown that ribosomal proteins can be used as candidate target genes for various environmental stresses (heavy metals, temperature) in aquaculture [41,42,43]. For instance, Rpl7 promotes the proliferation of myoblasts and inhibits their differentiation by affecting ribosomal translation [44]. During the embryonic development of Chironomus riparius, the expression levels of Rpl11 and Rpl13 significantly increase, indicating the crucial role of ribosomal proteins in development [45]. In our study, we identified seven differentially methylated ribosomal protein genes (RplP0, Rpl1, Rpl7, Rps16, Rpl18, Rpl18A, and Rps20) in intersex crabs compared to normal crabs. Rpl7 and Rpl18 are hypomethylated in intersex crabs compared to normal crabs, and their gene expression levels are higher in intersex crabs. These findings suggest that these genes may have important significance in the formation of intersex crabs. The exact mechanisms remain unclear and warrant further investigation. Ribosomal proteins play a crucial biological role in regulating ribosomal structure and protein synthesis to maintain cellular homeostasis, and processes such as replication, transcription, RNA synthesis, and DNA repair [46,47]. The upregulation of these ribosomal genes may be due to increased ribosome synthesis rates required to support the normal growth and development of intersex crabs [48].
Many DMGs in the PPI network were also enriched in metabolic pathways. Among the metabolic pathways enriched in the cm vs. fm and cm vs. mm groups, both common and distinct pathways are present, but all of these pathways are relevant to energy metabolism. This implies that the energy metabolism of intersex crabs is active and higher than that of normal crabs. Some studies have shown that crustaceans increase their energy metabolism in response to various stressors [49,50,51]. Among these DMGs, Acsl1 and Yip2 were verified, and promoter region of Acsl1 and Yip2 genes showed hypermethylation in normal crabs and hypomethylation in intersex crabs. Their expression levels showed an opposite trend. The Acsl1 gene encodes acyl-CoA synthetase, which functions in the conversion of free long-chain fatty acids into fatty acyl-CoA esters. It participates in the heat stress response by regulating fatty acid metabolism, which has been demonstrated in fish [52,53]. Yip2 is a homolog of mammalian acetyl-CoA acyltransferase 2 (ACAA2), which catalyzes mitochondrial fatty acid β-oxidation [54]. This suggests that the intersex crabs may be in a state of stress, and the above pathways are adapted to this state of stress by providing energy to the intersex crabs.
4. Materials and Methods
4.1. Ethics Statement
Animal handling and experimental procedures were performed according to the Guide for the Use of Experimental Animals of Ningbo University. The study protocol was approved by the Animal Ethics Committee of Ningbo University.
4.2. Sample Collection and Classification
Intersex individuals (cm) of E. sinensis used in this study were collected from Xinghua, Jiangsu province, China. All samples had part of the sexual characteristics of a female or male. The external sex characteristics of these samples included a penis, pleopods, gonopods, abdominal appendage, gonopores, and ostioles of unknown function in the abdomen. These intersex crabs were previously identified as genetically female by our laboratory’s sex marker [19]. Normal mature female (fm) and male (mm) crabs were selected as control. Three crabs were selected for each sample. The muscle tissues from the three crabs were pooled into one test tube and rapidly frozen in liquid nitrogen until further use.
4.3. Total Genomic DNA Extraction
Total genomic DNA (gDNA) was extracted from muscle tissue. Total DNA was extracted using QIAamp Fast DNA Tissue Kit (Qiagen, Dusseldorf, Germany) following the manufacturer’s procedure. The quantity of DNA was measured by spectrophotometer. When the A260/280 ratio is between 1.8 and 2.0, DNA is available.
4.4. Enzymatic Methyl Sequencing (EM-seq)
The gDNA was acoustically sheared to an average size of 200–280 bp (peak approximately 250 bp) using the Covaris instrument. Two internal controls were added to each sample: unmethylated lambda and methylated pUC19 CG sites. Tet methylcytosine dioxygenase 2 (TET2) and T4-BGT enzymes were used to protect 5 mC and 5 hmC from deamination. Subsequently, APOBEC3A converted cytosines (but not protected forms of cytosines) into uracil. The converted DNA was then prepared into single-stranded DNA libraries using the Accel-NGS DNA Library Kit (Swift Biosciences, Ann Arbor, MI, USA) according to the manufacturer’s protocol. Briefly, truncated adapter 1 was ligated to the 3′ end of the single-stranded DNA. A new DNA strand was generated through an extension step, followed by the addition of truncated adapter 2 to the 5′ end of the DNA. Finally, PCR amplification was performed using full-length adapter 1 and adapter 2 to increase the yield of DNA molecules. The final libraries were sequenced on the NovaSeq 6000 platform at LC Sciences (Houston, TX, USA), using paired-end 2 × 150 bp sequencing.
4.5. Identification of Differentially Methylated Regions (DMRs)
DMRs were analyzed separately for cm vs. fm and cm vs. mm. Firstly, adapter contamination, low-quality bases, and undetermined bases were removed from the reads using Cutadapt [55] and Perl scripts. Sequence quality was then validated with FastQC (v0.11.9). Reads that passed quality control were mapped to the reference genome using Bismark [56]. Post-alignment, duplicate reads were further removed using Samtools [57]. DNA methylation levels at each cytosine site were calculated as the C/(C + T) ratio (methylated reads divided by total reads) using the software of custom scripts and MethPipe (v3.4.3) [58]. DMRs were calculated using the R package MethylKit (v1.26.0) [59] (sliding window 1000 bp, overlap 500 bp, p-value < 0.05).
4.6. Functional Enrichment Analysis of Differentially Methylated Genes (DMGs) in Promoter Region
The promoter region of genes are methylated, resulting in changes in gene expression levels. We performed functional enrichment analysis of genes that were differentially methylated in the overall promoter region (2500 bp upstream of the transcription start site). All DMGs in the following text refer to genes that are differentially methylated in the overall promoter region. Gene Ontology (GO) enrichment analysis for DMGs was conducted using the GOseq R package [60], with a p-value threshold of <0.05. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of DMGs was assessed using KOBAS software (v3.0) [61].
4.7. PPI Analysis of DMGs in Promoter Region
We separately counted genes in the cm vs. fm and cm vs. mm comparison groups that had differential methylation levels in the promoter region and were enriched for specific GO and KEGG pathways. After obtaining the relevant genetic information, we conducted PPI analysis of DMGs in the promoter region obtained in comparison groups of cm vs. fm and cm vs. mm. The STRING database (v12.0) was employed for the construction of the PPI of proteins encoded by the intersection of DMGs. STRING is a database of known and predicted PPIs. The interactions include direct (physical) and indirect (functional) associations. The interactions were procured from genomic context, high-throughput experiments, co-expression, and previous knowledge. For a representative choice of database, either “Drosophila” or “Mus musculus” was selected, as no species-specific protein database is available for E. sinensis in STRING. Networks were therefore built to be representative of the phylum Arthropoda (with Drosophila showing most homology protein hits to E. sinensis) and Mammalia (using Mus musculus), respectively. Because the Drosophila database has the largest number of protein sequences and a closer evolutionary relationship, we finally chose Drosophila melanogaster as the reference organization. Parameters applied in STRING were “basic settings” and “medium confidence”. Cytoscape (v3.10.0) was applied to construct the DMG interaction network using the degree value as a criterion. A graph of theoretical clustering algorithm, molecular complex detection (MCODE) was utilized to analyze densely connected region. MCODE is part of the Cytoscape software. Further screening of the subnetwork containing hub DMGs was performed by the MCODE plugin of Cytoscape and the results were regarded as key functional DMGs. Module identification criteria included a degree cutoff of 2, node score cutoff of 0.2, k-core of 2, and maximum depth of 100. Significant modules were identified with MCODE score > 4.
4.8. Bisulfite Sequencing PCR (BSP)
Bisulfite sequencing PCR was used to validate the DNA methylation levels of selected candidate genes. Genomic DNA was treated with sodium bisulfite using the EZ DNA Methylation Kit (Zymo Research, Los Angeles, CA, USA). Following this, PCR amplification of the bisulfite-converted gDNA was carried out using Zymo Taq™ DNA Polymerase (Zymo Research). Details of the BSP primers are provided in Table 3. The PCR products were purified using a gel extraction kit (Sangon, Shanghai, China), then ligated and cloned into the pMD19-T vector (Sangon). Fifteen random clones from each sample were selected for DNA sequencing. Methylation analysis was performed using quantitative tools for bisulfite sequencing data (QUMA; v1.1.16).
4.9. Quantitative Real-Time PCR Analysis
Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was utilized to detect the expression profiles of key functional DMGs to further understand the possible molecular regulatory mechanisms. The qRT-PCR was performed using the 7500 Real-Time PCR system (Applied Biosystems, Foster City, CA, USA) with TB Green Premix Dimer Eraser (2×) (TaKaRa, Kyoto, Japan). The qRT-PCR was carried out in a total volume of 20 μL, containing 10 μL TB Green Premix Dimer Eraser (2×), 6.4 μL RNase-free water, 2 μL diluted cDNA, 0.5 μL forward/reverse primers (1 mM), and 0.4 μL ROX Reference Dye II. Each sample was run in triplicates. All primers used for the qRT-PCR are listed in Table 4. A standard curve was used to calculate the amplification efficiency and ensure the primer specificity of each primer pair. The relative expression levels of target genes were calculated by the 2−ΔΔCt method [62]. The 18S rRNA gene was used as an internal reference.
5. Conclusions
In summary, we have systematically characterized the whole-genome DNA methylation patterns in female, male, and intersex E. sinensis for the first time, and identified many DMRs and DMGs. Specifically, we identified 149 key DMGs in the promoter region between the cm and fm groups, and 110 key DMGs between the cm and mm groups. The key DMGs associated with the ribosomal pathway and metabolic pathway were enriched in both the cm vs. fm and cm vs. mm groups. Interestingly, only heat shock protein-associated DMGs were enriched in cm vs. fm. Our study provides information on the DNA methylation differences among females, males, and intersex crabs, which enriches the genome-wide DNA methylation library of E. sinensis.
Conceptualization, S.-J.F. and Z.-X.C.; methodology, S.-J.F., S.-C.S., M.-Q.N. and Y.-N.Y.; software, S.-J.F., S.-C.S. and M.-Q.N.; validation, S.-J.F.; formal analysis, S.-J.F.; investigation, S.-J.F.; resources, S.-J.F.; data curation, S.-J.F.; writing—original draft preparation, S.-J.F.; writing—review and editing, Z.-X.C.; visualization, S.-J.F.; supervision, Z.-X.C.; project administration, Z.-X.C.; funding acquisition, Z.-X.C. All authors have read and agreed to the published version of the manuscript.
This study strictly adhered to the “Guidelines for the Care and Use of Laboratory Animals” (No. 2011-2) published by the China Association for Laboratory Animal Science.
Not applicable.
The data used in this article have been submitted to the NCBI database. (SUB15104681).
The authors declare no conflicts of interest.
Footnotes
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.
Enlarge this image.Figure 1. Overall methylation levels for fm, cm, and mm. The proportions of mCG, mCHG, and mCHH in mC. The red, blue, and green colors represent methylated (m) CG/mC, mCHH/mC, and mCHG/mC, respectively. mC represents cytosine methylation. H can be any of A, T, or C. fm, female crabs. mm, male crabs. cm, intersex crabs.
Enlarge this image.Figure 2. DNA methylation levels of three methylated forms of CG, CHH, and CHG among genomic functional elements. (a) CG, (b) CHH, and (c) CHG. H can be any of A, T, or C. The horizontal coordinates represent the different functional elements on the genome: A, B, and C are promoters. D, F, and I are exons. E, G, and H represent introns. J represents downstream. The left vertical coordinates indicate their methylation levels in their respective backgrounds. TSS is the transcription start site. TES is the transcription end site. fm, female crabs. mm, male crabs. cm, intersex crabs.
Enlarge this image.Figure 3. Kyoto encyclopedia of genes and genomes (KEGG) analysis of differentially methylated genes (DMGs) in the promoter region of cm vs. fm and cm vs. mm. (a) Scatter plot of KEGG pathways in the top 10 of cm vs. fm. (b) KEGG pathways for the top 10 of cm vs. mm. The ordinate represents the enriched pathway and the abscissa represents the Rich factor of the corresponding pathway. The size of the spots represents the number of DMGs enriched in each pathway, and the color of the spots represents the corrected p-value of each pathway. The Rich factor represents the ratio of the number of DMGs mapped to a pathway to the total number of genes mapped to that pathway. The larger the enrichment factor, the higher the degree of enrichment. fm, female crabs. mm, male crabs. cm, intersex crabs.
Enlarge this image.Figure 4. Protein–protein interaction (PPI) network and KEGG enrichment information for the cm vs. fm and cm vs. mm comparison groups. (a) The cm vs. fm PPI network of DMGs in the promoter region. (b) The core network of cm vs. fm identification using the MCODE plugin and enrichment analysis via the KEGG. (c) The cm vs. mm PPI network of DMGs in the promoter region. (d) The core network of cm vs. mm identification using the MCODE plugin and enrichment analysis via the KEGG. fm, female crabs. mm, male crabs. cm, intersex crabs.
Enlarge this image.Figure 5. Methylation levels of DMGs in the promoter region. (a) Methylation levels of fm and cm key functional DMGs. (b) Methylation levels of mm and cm key functional DMGs. Hsc70, heat shock cognate protein 71 kDa. Hsp90, heat shock protein 90. Rpl18, 60S ribosomal protein L18. Acsl1, long-chain acyl-CoA synthetase 1. Rpl7, 60S ribosomal protein L7. Yip2, Yippee-interacting protein 2. fm, female crabs. mm, male crabs. cm, intersex crabs.
Enlarge this image.Figure 6. Validation of enzymatic methyl sequencing (EM-seq) data by bisulfite sequencing PCR (BSP). (a) BSP results of Hsc70 gene between fm and cm. (b) BSP results of Hsp90 gene between fm and cm. (c) BSP results of Rpl7 gene between mm and cm. (d) BSP results of Yip2 gene between mm and cm. Black and white circles indicate methylated and unmethylated CG, respectively. Cross indicates mismatch or gap in alignment. fm, female crabs. mm, male crabs. cm, intersex crabs.
Enlarge this image.Figure 7. The relative expression of mRNA for the candidate gene in the muscle tissue of E. sinensis was detected using the qRT-PCR. (a) qRT-PCR of fm and cm key function DMGs (b) qRT-PCR of mm and cm key function DMGs. 18S was used as an internal control. The values are represented as the mean ± SEM of three replicates. Hsc70, heat shock cognate protein 71 kDa. Hsp90, heat shock protein 90. Rpl18, 60S ribosomal protein L18. Acsl1, long-chain acyl-CoA synthetase 1. Rpl7, 60S ribosomal protein L7. Yip2, Yippee-interacting protein 2. fm, female crabs. mm, male crabs. cm, intersex crabs.
Sequencing data quality control results.
| Sample | Raw Data | Valid Data | Valid% | Q20% | Q30% | GC% | ||
|---|---|---|---|---|---|---|---|---|
| Read | Base | Read | Base | |||||
| fm | 385,311,926 | 57.80 G | 374,807,550 | 46.28 G | 97.27 | 97.19 | 93.44 | 22.31 |
| mm | 365,349,640 | 54.80 G | 351,897,242 | 44.02 G | 96.32 | 97.09 | 93.21 | 21.93 |
| cm | 330,057,373 | 49.51 G | 321,911,061 | 38.65 G | 97.54 | 97.24 | 93.60 | 22.13 |
Note: fm, female crabs; mm, male crabs; cm, intersex crabs.
Sequencing data mapping to the reference genome.
| Sample | Total Reads | Mapped Reads | Mapping Rate (%) | EM Conversion Rate (%) | Total mC (%) |
|---|---|---|---|---|---|
| fm | 374,807,550 | 171,443,461 | 45.74 | 97.13 | 4.5 |
| mm | 351,897,242 | 161,322,433 | 45.84 | 99.06 | 2.6 |
| cm | 321,911,061 | 139,530,777 | 43.33 | 98.35 | 3.5 |
Note: fm, female crabs; mm, male crabs; cm, intersex crabs. mC represents cytosine methylation.
Bisulfite Sequencing PCR primers.
| Gene | Nucleotide Sequence (From 5′ to 3′) |
|---|---|
| Hsc70 | inner-F: GTTATTTAGAGAGATAAAGATYGAGAG |
| inner-R: CATTCAAAAAAACAACAATTCT | |
| outer-F: TTTTGAGGTTGTTTTGTGTAGG | |
| outer-R: CCAAACACACACATTTAACAAAC | |
| Hsp90 | inner-F: TGTTTTTTAGTTTGTGGTTTGT |
| inner-R: AAAAATTTCCTCAAAATCAAAT | |
| outer-F: AGAYGTTTTTATTTTTATTTTAGAGGGG | |
| outer-R: CCCTAAAAAAACAATAAACATCCAAAAA | |
| Rpl7 | F: ATTGAAGGTTTGTAAATTGGAGT |
| R: TTATAAAAAAAAACACATCCCAA | |
| Yip2 | F: TAGTGTATATTTGAAGGTGTTGAA |
| R: AATAATAATAACCCAAACAAACAC |
Note: Hsc70, heat shock cognate protein 71 kDa. Hsp90, heat shock protein 90. Rpl7, 60S ribosomal protein L7. Yip2, Yippee-interacting protein 2.
Quantitative real-time PCR primers.
| Gene | Nucleotide Sequence (From 5′ to 3′) |
|---|---|
| Hsc70 | F: CGAGACCAAGTCGTTCTACCC |
| R: GCAGCACATTGAGACCAGAGA | |
| Hsp90 | F: TCGCAGTTCATTGGCTATCC |
| R: CCTCAATCTTGGGCTTCTCAT | |
| Rpl18 | F: GAACCCAAGTCGCAGGATG |
| R: CACCAGCCCAGAGAGTGAGAT | |
| Acsl1 | F: TCCGTAACCTTGTGGATGACTA |
| R: TTCTCCTTATCACACCCCTTTC | |
| Rpl7 | F: GCTCCGTTATTGTCCCTGC |
| R: GTAGCCCCAGGCAATGAAG | |
| Yip2 | F: CTGTCAGGAACATTCGCTTTG |
| R: CTTCATTCTGGGCGGTAGC |
Note: Hsc70, heat shock cognate protein 71 kDa. Hsp90, heat shock protein 90. Rpl18, 60S ribosomal protein L18. Acsl1, long-chain acyl-CoA synthetase 1. Rpl7, 60S ribosomal protein L7. Yip2, Yippee-interacting protein 2.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Allis, C.D.; Jenuwein, T. The Molecular Hallmarks of Epigenetic Control. Nat. Rev. Genet.; 2016; 17, pp. 487-500. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27346641]
2. Li, G.-L.; Xu, Y.-J.; Huang, X.-M.; Xiao, J.; Nong, S.; Li, C.-G. MeDIP-Seq Reveals the Features of Mitochondrial Genomic Methylation in Immature Testis of Chinese Mitten Crab Eriocheir sinensis. Mitochondrial DNA Part A; 2018; 29, pp. 335-339. [DOI: https://dx.doi.org/10.1080/24701394.2016.1278537] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28129733]
3. Hearn, J.; Pearson, M.; Blaxter, M.; Wilson, P.J.; Little, T.J. Genome-Wide Methylation Is Modified by Caloric Restriction in Daphnia magna. BMC Genom.; 2019; 20, 197. [DOI: https://dx.doi.org/10.1186/s12864-019-5578-4]
4. Wang, X.; Li, A.; Wang, W.; Que, H.; Zhang, G.; Li, L. DNA Methylation Mediates Differentiation in Thermal Responses of Pacific Oyster (Crassostrea gigas) Derived from Different Tidal Levels. Heredity; 2021; 126, pp. 10-22. [DOI: https://dx.doi.org/10.1038/s41437-020-0351-7]
5. Wang, H.; Cui, J.; Qiu, X.; Wang, X. Differences in DNA Methylation between Slow and Fast Muscle in Takifugu rubripes. Gene; 2021; 801, 145853. [DOI: https://dx.doi.org/10.1016/j.gene.2021.145853] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34274464]
6. 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. [DOI: https://dx.doi.org/10.1186/s13293-022-00484-6]
7. Ford, A.T. Intersexuality in Crustacea: An Environmental Issue?. Aquat. Toxicol.; 2012; 108, pp. 125-129. [DOI: https://dx.doi.org/10.1016/j.aquatox.2011.08.016]
8. Takahashi, T.; Araki, A.; Nomura, Y.; Koga, M.; Arizono, K. The Occurrence of Dual-Gender Imposex in Japanese Freshwater Crab. J. Health Sci.; 2000; 46, pp. 376-379. [DOI: https://dx.doi.org/10.1248/jhs.46.376]
9. Rodgers-Gray, T.P.; Smith, J.E.; Ashcroft, A.E.; Isaac, R.E.; Dunn, A.M. Mechanisms of Parasite-Induced Sex Reversal in Gammarus duebeni. Int. J. Parasitol.; 2004; 34, pp. 747-753. [DOI: https://dx.doi.org/10.1016/j.ijpara.2004.01.005]
10. Parnes, S.; Khalaila, I.; Hulata, G.; Sagi, A. Sex Determination in Crayfish: Are Intersex Cherax quadricarinatus (Decapoda, Parastacidae) Genetically Females?. Genet. Res.; 2003; 82, pp. 107-116. [DOI: https://dx.doi.org/10.1017/S0016672303006372]
11. Abdel-moneim, A.; Coulter, D.P.; Mahapatra, C.T.; Sepúlveda, M.S. Intersex in Fishes and Amphibians: Population Implications, Prevalence, Mechanisms and Molecular Biomarkers. J. Appl. Toxicol.; 2015; 35, pp. 1228-1240. [DOI: https://dx.doi.org/10.1002/jat.3204]
12. Kidd, K.A.; Blanchfield, P.J.; Mills, K.H.; Palace, V.P.; Evans, R.E.; Lazorchak, J.M.; Flick, R.W. Collapse of a Fish Population after Exposure to a Synthetic Estrogen. Proc. Natl. Acad. Sci. USA; 2007; 104, pp. 8897-8901. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17517636]
13. Inui, N.; Oguchi, K.; Shinji, J.; Okanishi, M.; Shimomura, M.; Miura, T. Parasitism-Induced Intersexuality in a Sexually Dimorphic Varunid Crab, Ptychognathus ishii (Decapoda: Varunidae). Zool. Sci.; 2021; 38, pp. 416-426.
14. Souissi, A.; Souissi, S.; Devreker, D.; Hwang, J.-S. Occurence of Intersexuality in a Laboratory Culture of the Copepod Eurytemora affinis from the Seine Estuary (France). Mar. Biol.; 2010; 157, pp. 851-861.
15. Short, S.; Yang, G.; Kille, P.; Ford, A.T. A Widespread and Distinctive Form of Amphipod Intersexuality Not Induced by Known Feminising Parasites. Sex. Dev.; 2012; 6, pp. 320-324. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23154238]
16. Shao, C.; Li, Q.; Chen, S.; Zhang, P.; Lian, J.; Hu, Q.; Sun, B.; Jin, L.; Liu, S.; Wang, Z. Epigenetic Modification and Inheritance in Sexual Reversal of Fish. Genome Res.; 2014; 24, pp. 604-615.
17. Todd, E.V.; Ortega-Recalde, O.; Liu, H.; Lamm, M.S.; Rutherford, K.M.; Cross, H.; Black, M.A.; Kardailsky, O.; Marshall Graves, J.A.; Hore, T.A. Stress, Novel Sex Genes, and Epigenetic Reprogramming Orchestrate Socially Controlled Sex Change. Sci. Adv.; 2019; 5, eaaw7006.
18. 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, pp. 645-648.
19. Lu, W.; Aamoum, O.; Feng, T.; Mo, N.; Han, R.; Shao, S.; Cui, Z. The Comparison of Morphological Characters, Ovarian Histology and Gene Expression Levels between Intersex and Normal Eriocheir sinensis (Decapoda, Varunidae). Crustaceana; 2024; 97, pp. 1303-1322.
20. Zhang, P.; Yang, Y.; Xu, Y.; Cui, Z. Analyses of the Dmrt family in a decapod crab, Eriocheir sinensis uncover new facets on the evolution of DM domain genes. Front. Physiol.; 2023; 14, 1201846.
21. Zhu, D.; Feng, T.; Mo, N.; Han, R.; Lu, W.; Shao, S.; Cui, Z. New insights for the regulatory feedback loop between type 1 crustacean female sex hormone (CFSH-1) and insulin-like androgenic gland hormone (IAG) in the Chinese mitten crab (Eriocheir sinensis). Front. Physiol.; 2022; 13, 1054773. [DOI: https://dx.doi.org/10.3389/fphys.2022.1054773]
22. Zhang, W.; Zhang, S.; Xu, Y.; Ma, Y.; Zhang, D.; Li, X.; Zhao, S. The DNA Methylation Status of Wnt and Tgfβ Signals Is a Key Factor on Functional Regulation of Skeletal Muscle Satellite Cell Development. Front. Genet.; 2019; 10, 220. [DOI: https://dx.doi.org/10.3389/fgene.2019.00220]
23. Cribiu, P.; Chaumot, A.; Geffard, O.; Ravanat, J.-L.; Bastide, T.; Delorme, N.; Quéau, H.; Caillat, S.; Devaux, A.; Bony, S. Natural Variability and Modulation by Environmental Stressors of Global Genomic Cytosine Methylation Levels in a Freshwater Crustacean, Gammarus fossarum. Aquat. Toxicol.; 2018; 205, pp. 11-18. [DOI: https://dx.doi.org/10.1016/j.aquatox.2018.09.015] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30300817]
24. Guo, J. Elucidating the Role of DNA Methylation in Low-Salinity Adaptation of the Swimming Crab Portunus trituberculatus. Aquaculture; 2025; 595, 741707.
25. Feng, S.; Cokus, S.J.; Zhang, X.; Chen, P.-Y.; Bostick, M.; Goll, M.G.; Hetzel, J.; Jain, J.; Strauss, S.H.; Halpern, M.E. et al. Conservation and Divergence of Methylation Patterning in Plants and Animals. Proc. Natl. Acad. Sci. USA; 2010; 107, pp. 8689-8694.
26. Head, J.A. Patterns of DNA Methylation in Animals: An Ecotoxicological Perspective. Integr. Comp. Biol.; 2014; 54, pp. 77-86. [DOI: https://dx.doi.org/10.1093/icb/icu025] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24785828]
27. Pan, Y.; Chen, L.; Cheng, J.; Zhu, X.; Wu, P.; Bao, L.; Chu, W.; He, S.; Liang, X.; Zhang, J. Genome-Wide DNA Methylation Profiles Provide Insight into Epigenetic Regulation of Red and White Muscle Development in Chinese Perch Siniperca chuatsi. Comp. Biochem. Physiol. Part B Biochem. Mol. Biol.; 2021; 256, 110647.
28. Liu, Q.; Liu, X.; Wang, G.; Wu, F.; Hou, Y.; Liu, H. Genome-Wide DNA Methylation Analysis of Astragalus and Danshen on the Intervention of Myofibroblast Activation in Idiopathic Pulmonary Fibrosis. BMC Pulm. Med.; 2023; 23, 325. [DOI: https://dx.doi.org/10.1186/s12890-023-02601-6]
29. Fan, Y.; Liang, Y.; Deng, K.; Zhang, Z.; Zhang, G.; Zhang, Y.; Wang, F. Analysis of DNA Methylation Profiles during Sheep Skeletal Muscle Development Using Whole-Genome Bisulfite Sequencing. BMC Genom.; 2020; 21, 327.
30. Begue, G.; Raue, U.; Jemiolo, B.; Trappe, S. DNA Methylation Assessment from Human Slow-and Fast-Twitch Skeletal Muscle Fibers. J. Appl. Physiol.; 2017; 122, pp. 952-967.
31. Vogt, G. Studying Phenotypic Variation and DNA Methylation across Development, Ecology and Evolution in the Clonal Marbled Crayfish: A Paradigm for Investigating Epigenotype-Phenotype Relationships in Macro-Invertebrates. Sci. Nat.; 2022; 109, 16.
32. Boopathy, L.R.A.; Jacob-Tomas, S.; Alecki, C.; Vera, M. Mechanisms Tailoring the Expression of Heat Shock Proteins to Proteostasis Challenges. J. Biol. Chem.; 2022; 298, 101796.
33. Jacob, P.; Hirt, H.; Bendahmane, A. The Heat-shock Protein/Chaperone Network and Multiple Stress Resistance. Plant Biotechnol. J.; 2017; 15, pp. 405-414. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27860233]
34. Kohno, S.; Katsu, Y.; Urushitani, H.; Ohta, Y.; Iguchi, T.; Guillette, L.J., Jr. Potential Contributions of Heat Shock Proteins to Temperature-Dependent Sex Determination in the American Alligator. Sex. Dev.; 2010; 4, pp. 73-87.
35. Uchimura, T.; Hara, S.; Yazawa, T.; Kamei, Y.; Kitano, T. Involvement of Heat Shock Proteins on the Transcriptional Regulation of Corticotropin-Releasing Hormone in Medaka. Front. Endocrinol.; 2019; 10, 529.
36. Wang, C.; Chen, X.; Dai, Y.; Zhang, Y.; Sun, Y.; Cui, X. Comparative Transcriptome Analysis of Heat-Induced Domesticated Zebrafish during Gonadal Differentiation. BMC Genom. Data; 2022; 23, 39.
37. Liu, Q.; Wang, Y.; Tan, L.; Ma, W.; Zhao, X.; Shao, C.; Wang, Q. The Role of the Heat Shock Cognate Protein 70 Genes in Sex Determination and Differentiation of Chinese Tongue Sole (Cynoglossus semilaevis). Int. J. Mol. Sci.; 2023; 24, 3761. [DOI: https://dx.doi.org/10.3390/ijms24043761]
38. Wang, Q.; Liu, K.; Feng, B.; Zhang, Z.; Wang, R.; Tang, L.; Li, W.; Li, Q.; Piferrer, F.; Shao, C. Transcriptome of Gonads From High Temperature Induced Sex Reversal During Sex Determination and Differentiation in Chinese Tongue Sole, Cynoglossus semilaevis. Front. Genet.; 2019; 10, 1128. [DOI: https://dx.doi.org/10.3389/fgene.2019.01128]
39. Knee, D.A.; Froesch, B.A.; Nuber, U.; Takayama, S.; Reed, J.C. Structure-Function Analysis of Bag1 Proteins: Effects on Androgen Receptor Transcriptional Activity. J. Biol. Chem.; 2001; 276, pp. 12718-12724.
40. Knorr, E.; Vilcinskas, A. Post-Embryonic Functions of HSP90 in Tribolium castaneum Include the Regulation of Compound Eye Development. Dev. Genes Evol.; 2011; 221, pp. 357-362.
41. Khan, M.M.; Ashiqul Islam, S.M.; Mustafa, A. Identification of Stress Related Molecular Biomarkers in Zebrafish Employing an In-Silico Approach to Access Toxicity Based Risks in Aquaculture. Poult. Fish. Wildl. Sci.; 2015; 3, 1000137.
42. Quinn, N.L.; McGowan, C.R.; Cooper, G.A.; Koop, B.F.; Davidson, W.S. Ribosomal Genes and Heat Shock Proteins as Putative Markers for Chronic, Sublethal Heat Stress in Arctic Charr: Applications for Aquaculture and Wild Fish. Physiol. Genom.; 2011; 43, pp. 1056-1064.
43. Vieira, J.C.S.; Braga, C.P.; de Oliveira, G.; Padilha, C. do C.F.; de Moraes, P.M.; Zara, L.F.; Leite, A. de L.; Buzalaf, M.A.R.; Padilha, P. de M. Mercury Exposure: Protein Biomarkers of Mercury Exposure in Jaraqui Fish from the Amazon Region. Biol. Trace Elem. Res.; 2018; 183, pp. 164-171.
44. Yan, J.; Yang, Y.; Fan, X.; Tang, Y.; Tang, Z. Sp1-Mediated circRNA circHipk2 Regulates Myogenesis by Targeting Ribosomal Protein Rpl7. Genes; 2021; 12, 696. [DOI: https://dx.doi.org/10.3390/genes12050696]
45. Martínez-Guitarte, J.L.; Planelló, R.; Morcillo, G. Characterization and Expression during Development and under Environmental Stress of the Genes Encoding Ribosomal Proteins L11 and L13 in Chironomus riparius. Comp. Biochem. Physiol. Part B Biochem. Mol. Biol.; 2007; 147, pp. 590-596. [DOI: https://dx.doi.org/10.1016/j.cbpb.2007.03.015]
46. de la Cruz, J.; Karbstein, K.; Woolford, J.L., Jr. Functions of Ribosomal Proteins in Assembly of Eukaryotic Ribosomes in Vivo. Annu. Rev. Biochem.; 2015; 84, pp. 93-129.
47. Zhou, X.; Liao, W.-J.; Liao, J.-M.; Liao, P.; Lu, H. Ribosomal Proteins: Functions beyond the Ribosome. J. Mol. Cell Biol.; 2015; 7, pp. 92-104. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25735597]
48. Hu, W.; Zhu, Q.-L.; Zheng, J.-L.; Wen, Z.-Y. Cadmium Induced Oxidative Stress, Endoplasmic Reticulum (ER) Stress and Apoptosis with Compensative Responses towards the up-Regulation of Ribosome, Protein Processing in the ER, and Protein Export Pathways in the Liver of Zebrafish. Aquat. Toxicol.; 2022; 242, 106023.
49. Wang, J.; Liu, Q.; Zhang, X.; Gao, G.; Niu, M.; Wang, H.; Chen, L.; Wang, C.; Mu, C.; Wang, F. Metabolic Response in the Gill of Portunus Trituberculatus Under Short-Term Low Salinity Stress Based on GC-MS Technique. Front. Mar. Sci.; 2022; 9, 881016.
50. Jiang, Y.; Liu, X.; Shang, Y.; Li, J.; Gao, B.; Ren, Y.; Meng, X. Physiological and Transcriptomic Analyses Provide Insights into Nitrite Stress Responses of the Swimming Crab Portunus trituberculatus. Mar. Biotechnol.; 2024; 26, pp. 1040-1052.
51. Zhang, X.; Yuan, J.; Zhang, X.; Yu, Y.; Li, F. Comparative Transcriptomic Analysis Unveils a Network of Energy Reallocation in Litopenaeus vannamei Responsive to Heat-Stress. Ecotoxicol. Environ. Saf.; 2022; 238, 113600. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35526454]
52. Zhao, T.; Ma, A.; Huang, Z.; Liu, Z.; Sun, Z.; Zhu, C.; Yang, J.; Li, Y.; Wang, Q.; Qiao, X. et al. Transcriptome Analysis Reveals That High Temperatures Alter Modes of Lipid Metabolism in Juvenile Turbot (Scophthalmus maximus) Liver. Comp. Biochem. Physiol. Part D Genom. Proteom.; 2021; 40, 100887.
53. Li, Q.-Q.; Zhang, J.; Wang, H.-Y.; Niu, S.-F.; Wu, R.-X.; Tang, B.-G.; Wang, Q.-H.; Liang, Z.-B.; Liang, Y.-S. Transcriptomic Response of the Liver Tissue in Trachinotus ovatus to Acute Heat Stress. Animals; 2023; 13, 2053. [DOI: https://dx.doi.org/10.3390/ani13132053] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37443851]
54. Sun, X.; Shen, J.; Perrimon, N.; Kong, X.; Wang, D. The Endoribonuclease Arlr Is Required to Maintain Lipid Homeostasis by Downregulating Lipolytic Genes during Aging. Nat. Commun.; 2023; 14, 6254.
55. Martin, M. Cutadapt Removes Adapter Sequences from High-Throughput Sequencing Reads. EMBnet. J.; 2011; 17, pp. 10-12.
56. Krueger, F.; Andrews, S.R. Bismark: A Flexible Aligner and Methylation Caller for Bisulfite-Seq Applications. Bioinformatics; 2011; 27, pp. 1571-1572.
57. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. 1000 Genome Project Data Processing Subgroup the Sequence Alignment/Map Format and SAMtools. Bioinformatics; 2009; 25, pp. 2078-2079. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp352]
58. Song, Q.; Decato, B.; Hong, E.E.; Zhou, M.; Fang, F.; Qu, J.; Garvin, T.; Kessler, M.; Zhou, J.; Smith, A.D. A Reference Methylome Database and Analysis Pipeline to Facilitate Integrative and Comparative Epigenomics. PLoS ONE; 2013; 8, e81148.
59. Akalin, A.; Kormaksson, M.; Li, S.; Garrett-Bakelman, F.E.; Figueroa, M.E.; Melnick, A.; Mason, C.E. methylKit: A Comprehensive R Package for the Analysis of Genome-Wide DNA Methylation Profiles. Genome Biol.; 2012; 13, R87.
60. 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.
61. 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, pp. 3787-3793. [DOI: https://dx.doi.org/10.1093/bioinformatics/bti430] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15817693]
62. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods; 2001; 25, pp. 402-408. [DOI: https://dx.doi.org/10.1006/meth.2001.1262] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11846609]
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.