- 1Yangtze River Fisheries Research Institute, Chinese Academy of Fisheries Science, Wuhan, China
- 2College of Fisheries and Life Science, Shanghai Ocean University, Shanghai, China
Sex dimorphism is a key feature of Chinese soft-shelled turtle (Pelodiscus sinensis). The males (M) have higher econosmic value than females (F) due to wider calipash and faster growth. Exogenous hormones like estradiol and methyltestosterone can induce sexual reversal to form new phenotypes (pseudo-female, PF; pseudo-male, PM) without changing the genotype. The possibility of inducing sexual reversal is particularly important in aquaculture breeding, but the underlying biological mechanisms remain unclear. Here we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in P. sinensis. Nanopore sequencing of the four gender types (M, F, PF, and PM) showed that the distribution of read length and gene expression was more similar between same-sex phenotypes than same-sex genotypes. Compared to turtles with an M phenotype, alternative splicing was more pronounced in F turtles, especially at alternative 3′ splice sites, alternative 5′ splice sites, and alternative first exons. Furthermore, the two RNA methylation modifications m5C and m6A were differentially distributed across gender phenotypes, with the M type having more modification sites in coding sequence regions, but fewer modification sites in 3′UTR regions. Quantitative analysis of enriched m6A RNAs revealed that the N6-methylated levels of Odf2, Pacs2, and Ak1 were significantly higher in M phenotype individuals, while the N6-methylated levels of Ube2o were reduced after sexual reversal from both M and F phenotypes. Taken together, these findings reveal an important role of epigenetics during sexual reversal in Chinese soft-shelled turtles.
Introduction
Chinese soft-shelled turtle, an economically valuable aquatic animal, exhibits obvious sexual dimorphism (Liang et al., 2018; Ma et al., 2021; Nagahama et al., 2021). Male individuals are characterized by wider calipash and faster growth rate, and are more popular in the market (Liang et al., 2019a; Zhou et al., 2021). The mechanism of sex determination in the Chinese soft-shelled turtle is genotypic sex determination (GSD), not temperature-dependent sex determination (TSD), and follows the ZZ/ZW chromosomal system where females have a pair of dissimilar ZW micro-sex chromosomes and males have two similar ZZ chromosomes (Kawai et al., 2007; Kawagoshi et al., 2009).
Exogenous hormones can induce sexual reversal in Chinese soft-shelled turtle during the period of sex differentiation (Schulz et al., 2007; Liang et al., 2019b). The balance between estrogens and androgens determines the direction of gonadal development in vertebrates (Ma et al., 2021). Exogenous hormones such as estradiol can induce the transformation from a male phenotype into a pseudo-female (PF) phenotype without changing the genetic gender. The transition from male to female is generally referred to as the male reversal process, while the transition from female to male is called the female reversal process. Male reversal of P. sinensis is of great significance to the aquaculture industry because a sexually mature PF can be crossed with a genetic male (ZZ) to obtain all-male offspring (Zhou et al., 2021). However, the biological underpinnings of sexual reversal in Pelodiscus sinensis are still little explored.
Transcriptome analysis is an effective tool for gaining insights into biological processes (Kukurba and Montgomery, 2015). Identification and quantification of various transcripts can improve our understanding of the gene functions and the complex mechanisms involved in sex control and sex differentiation, but it remains a challenging task (Li et al., 2020). Previous studies have typically used short-read (<200 nt) RNA sequencing (RNA-seq) followed by reassembly, making it difficult to sequence repetitive regions, resolve large structural variations and homologous gene segments, and differentiate isoforms.
An alternative approach is the long-range direct RNA-seq technology offered by Oxford Nanopore Technologies (ONT) which has no theoretical upper limit for read length and can sequence transcripts from end to end at a single-molecule level (Jenjaroenpun et al., 2018; Lubuta et al., 2019). ONT sequencing is a new generation of single-molecule real-time electrical signal sequencing technology based on nanopores, and its sequencing principles are the same on all platforms. The DNA/RNA double strand is led by a motor protein to bind to the nanopore embedded in a biological membrane. Under the action of the voltage difference between the two sides of the biological membrane, the DNA/RNA strand passes through the nanopore channel at a certain rate. Because of the differences in the chemical properties of different bases in the DNA/RNA chain, when a single base passes through the nanopore channel, it will cause specific changes to the electrical current that flows though the nanopore. These electrical disruptions are then decoded using base calling algorithms to determine the DNA/RNA sequence in real time (Lu et al., 2016; Karamitros and Magiorkinis, 2018; Dumschott et al., 2020).
This technology overcomes many of the shortcomings of traditional RNA-seq techniques as it allows the direct sequencing of RNA transcripts without the need for inverse transcription and amplification (Kono and Arakawa, 2019). The long sequencing reads provided by this method offer a promising tool for the detection of alternative splicing, fusion genes and novel isoforms as well as the accurate characterization of methylation modification sites (Roach et al., 2020). Compared to standard transcriptome analysis, direct RNA-seq is better suited to identify novel splicing isoforms (Li et al., 2020). Nanopore-based sequencing methods have been used to annotate transcriptome structures in a variety of organisms ranging from the relatively simple Saccharomyces cerevisiae to complex human cell lines (Byrne et al., 2017; Garalde et al., 2018; Bayega et al., 2021). Direct DNA sequencing can unlock a wealth of information, such as the nuances of alternative splicing, and thereby contribute to the understanding of molecular mechanisms (Garalde et al., 2018). Many novel alternative splicing events continue to be discovered in a variety of biological processes including sexual differentiation and gonadal development (Ule and Blencowe, 2019).
Epigenetics is the study of heritable changes in gene expression that are not the result of alterations in the nucleotide sequence of the gene (Ladd-Acosta and Fallin, 2015; Kaliman, 2019). Epigenetic changes are defined as heritable cellular information which are not genetic and can be transmitted through cell division (Wu et al., 2015). RNA methylation modification is an important branch of epigenetics (Wang and Wang, 2015). N6-methyladenosine (m6A) and C5-methylcytidine (m5C) are two common RNA post-transcriptional modifications in eukaryotes. RNA methylation plays an important role in regulating gene expression and in processes related to editing, stability and degradation of transcripts. In particular, dynamic methylation of adenine and cytosine is essential for cells to adapt to their environment and for the development of complex organisms from a single cell (Traube and Carell, 2017; Yang et al., 2017; Yang, Y. et al., 2017). For example, m6A can regulate the oocyte-to-embryo transition and the reprogramming of somatic cells into induced pluripotent stem cells (Yu et al., 2021). Wang et al. found that the presence of sex-specific RNA modifications may be involved in the regulation of gonadal development and gametogenesis in Paralichthys olivaceus (Wang et al., 2020). In the case of m5C, research has shown that m5C methyltransferases are closely associated with various human diseases including cancers (Liu et al., 2017), and m5C modifications are thought to play a role in intracellular and nucleocytoplasmic transport (Sanchez-Vasquez et al., 2018). Although these two RNA methylation modifications play a crucial role in the growth of animals, there are few reports on the epigenetic relationship between RNA methylation modifications and the sexual reversal of Chinese soft-shelled turtles.
In this study, we employed the Nanopore long-read direct RNA sequencing platform to investigate the transcriptomic complexity of P. sinensis genders. Moreover, we comprehensively analyzed alternative splicing events and different RNA methylation modifications, including m6A and m5C, which revealed an important role of epigenetics in Chinese soft-shelled turtle gonadal development.
Materials and Methods
P. sinensis Sample Collection
Female and male P. sinensis were cultured in ponds, which were located in Anhui Xijia Agricultural Development Co. Ltd. (Bengbu, China). The turtles are kept in natural ponds with commercial feeds (Jinjia, China) three times a day. The stocking density of juvenile turtles is about 4/m2. The transition to PF or PM was induced at the embryonic stage 15th day before sex differentiation by estradiol (E2) or methyltestosterone (MT). We injected 5 µL of 10 mg/mL E2 (E110147, Aladdin, shanghai, China) or MT (M830053, Macklin, shanghai, China) into the fertilized eggs whose eggshells were softened by hydrochloric acid (HCl) to induce the sex reversal of Chinese soft-shelled turtle (Zhou et al., 2021; Othman et al., 2022). The gonads of pseudo-female and pseudo-male individuals began to differentiate at the embryonic stage 15th day and mature after 2 years. After the turtles were 5 months old, the physiological gender and genetic gender were identified by the tail and sex specific markers (Liang et al., 2019b). The turtles were judged to be male if the tail exceeded the calipash edge, otherwise it was judged to be a female. The genetic gender was determined by sex-specific markers (4085-f/r, Supplementary Table S4) using PCR. Those individuals with inconsistent physiological gender and genetic gender had a sex reversal. After 2 years of cultivation of the turtles, the female, male, PF and PM individuals which every group contains three turtles were used in this study (Supplementary Figure S1). The gonad tissues were collected for further analysis after MS-222 (Sigma, St. Louis, MO, United States) treatment. The turtles were injected intraperitoneally with 600 mg/kg of buffered MS-222 (Keifer and Zheng, 2017). All samples were stored at −80°C.
RNA Extraction
Total RNA was isolated using the RNA prep pure kit (Tiangen, China), and DNA was removed with DNase I (Sigma, St. Louis, MO, United States) following the protocol Green and Sambrook, 2019. The RNA quality was assessed using NanoDrop One spectrophotometer (NanoDrop Technologies, Wilmington, DE, United States) and Qubit 3.0 Fluorometer (Life Technologies, Carlsbad, CA, United States).
Library Preparation and Nanopore Sequencing
Libraries of direct RNA sequencing were constructed following the ONT Direct RNA Sequencing protocol (Lu et al., 2016; Grünberger et al., 2022). Approximately 30 μg total RNA was enriched for mRNA using NEBNext Poly(A) mRNA Magnetic Isolation Module (E7490S).
The Nanopore RTA adapter (SQK-RNA002) was added to the Poly(A) mRNA using T4 DNA Ligase (NEB) at 25°C for 10 min, and then the first strand cDNA was synthesized using SuperScript III reverse transcriptase (1808004, Thermo Fisher Scientific, Waltham, MA, United States) to form RNA/DNA duplexes and relax the secondary structure of RNA. The RNA-cDNA hybrid products were purified using 1.8 × Agencourt RNA Clean XP beads and eluted using 23 µL Nuclease-free water (Promega, Madison, WI, United States).
Sequencing adapters were ligated to the mRNA using T4 DNA ligase (NEB) in a total volume of 40 μL, including 23 µL reverse-transcribed mRNA, 8 µL NEBNext Quick Ligation Reaction Buffer (NEB), 6 µL RNA Adapter (RMX; SQK-RNA002), and 3 µL T4 DNA ligase. Then, Agencourt RNA Clean XP beads were added and eluted using 41 µLElution Buffer, and 1 μL eluted sample was quantified by Qubit. Sequencing libraries (145 µL) were prepared by mixing 40 µL mRNA with sequencing adapters, 75 µL RRB (SQK-RNA002), and 30 µL Nuclease-free water. The libraries were loaded on ONT R9.4 flowcells and sequenced according to the standard protocol using the Nanopore PromethION sequencer (Oxford Nanopore Technologies, Oxford, United Kingdom) for 48 h. Three independent biological replicates were sequenced. The base calling was performed concurrent with sequencing based on the Albacore, which is integrated within the MinION software (MinKNOW, v1.10.23). Only “pass” reads were used for subsequent analyses.
Read Filtering
The reads were filtered with the Miniknow default cutoff of the adaptors, a minimum average q-score of 7 and the length <50 bp. And the data were written out in the FASTQ format using NanoFilt (version:2.6.0). Then the effective reads were further analyzed.
Mapping and Qualitative Analysis
The clean reads were mapped to the Chinese soft-shelled turtle genome (GCF_000230,535.1) using Minimap2. The parameters were set as “-ax splice -uf -k14”. The blasted transcripts were classified into “unmapped,” “Forward mapped,” and “Reverse mapped”.
Alternative Splicing Analysis
We utilized the SUPPA2 method (https://github.com/comprna/SUPPA) to analyze the alternative splicing of each individual. In total, seven types of alternative splicing events were detected, including skipping exon, mutually exclusive exons, alternative 5′ splice-site, alternative 3’ splice-site, retained intron, alternative first exon and alternative last exon. The analysis of alternative splicing events was carried out by SUPPA2.
Characterization of Poly(A) Tail Length
Poly(A) tail lengths of each read were calculated from raw signals using the program nanopolish with “polya” module (version 0.12.5) (Loman and Watson, 2015). Only poly(A) tail length estimates marked as “PASS” were considered for downstream analyses.
Identification of RNA Methylation Sites
The modifications sites of RNA sequences were identified with Tombo version 1.5 with alternative model. The models of “5mC” and “de novo” were implemented separately to detect possible modifications in each read. The score on each site indicated the fraction of a possible modification on a given site. For plotting the modification coverage along the gene body, the modification coverage was normalized for each isoform using a “w0” method with a bin size of 5 nt with EnrichedHeatmap (Gu et al., 2018). Only the isoforms with both 5′ and 3′UTRs longer than 50 nt were used in the calculation. The identification of m6A sites was carried out based on the de novo of Tombo (Lorenz et al., 2020).
Enrichment of m6A RNA From Total RNAs
The EpiQuik CUT&RUN m6A RNA Enrichment kit (EpiGentek, P-9018) was used to hatch the m6A RNAs from total RNAs which were extracted from gonad tissues from the different samples. In this reaction, RNA sequences in both ends of the target m6A-containing regions are cleaved/removed, and the target m6A-containing fragments are pulled down using a beads-bound m6A capture antibody. The enriched RNA is then released, purified, and eluted. The eluted RNAs were used for further quantitative analyses.
The subsequent quantification of cDNA products was performed as described previously (Zhou et al., 2021). The primers for target gene qRT-PCR are provided in Supplementary Table S4. Relative m6A mRNA expression levels were calculated by the 2−(∆∆Ct) method (Schmittgen, 2008).
Statistical Analysis
All the experimental data from at least three independent experiments were analyzed using GraphPad Prism 7.0 software (San Diego, CA, United States) and were presented as mean ± SD. Student’s t-test or One-Way ANOVA were performed to compare the pairwise differences between groups.
Results
Statistics of Direct RNA Sequencing Data of Different Samples
To assess the quality of the RNA sequencing data, we identified basic data indicators such as the read number, the assembly quality metric N50, and mean length of sequencing reads, and compared these between different physiological and genetic genders (Figure 1). Since the raw sequencing data may contain low-quality sequences and adapter sequences, we filtered the data to obtain clean reads and thereby ensure the reliability of the analysis results. After filtering, 62,672 clean reads were mapped to the Chinese soft-shelled turtle genome (version 1.1) using Minimap2 (v.2.17-r941). The results showed that forward mapped reads and reverse mapped reads occupied 99.19% of total clean reads (Figure 1A, Supplementary Table S1). A histogram displaying the distribution of read lengths across quality passing sequences is shown in Supplementary Figure S2. For each gender, we obtained >1.5 million reads with an average read length of 1,223, 844, 1,341, 979 for female (F), male (M), PF and PM, respectively. The N50 statistic of the four gender types (F, M, PF, and PM) were 1760, 1,197, 1811 and 1,339. The N50 and mean length of PM were close to those of M while far from those of F. On the contrary, the N50 and mean length of PF were close to those of F but far from those of M. These characteristics of basic data among four types of Chinese soft-shelled turtle showed that the changes brought about by sex reversal (Figure 1B, Supplementary Table S2). TPM box plot showed the median log2 (TPM) values for all transcripts in each group were different and the values of M and PM were lower than F and PF (Figure 1C). We conducted correlation analysis on the different samples, and found a high degree of correlation between most of them, meaning that they met the basic requirements for further analysis (Figure 1D). In general, the sequencing data had high depth and confidence, and the M and PM groups were different from the F and PF in read length and the overall expression level of transcripts.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876045/fcell-10-876045-HTML/image_m/fcell-10-876045-g001.jpg)
FIGURE 1. Statistics of direct RNA sequencing data of different samples. (A) The mapped diagram of sequencing data. (B) The histogram of the distribution of N50 and mean length. Data were expressed as mean ± SD. Different superscripts indicate significant difference (p < 0.05). (C) The TPM box plot of different samples. (D) The correlation analysis on the different samples. F1: Female-1, F2: Female-2, F3: Female-3, M1: Male-1, M2: Male-2, M3: Male-3, PF1: Pseudo-female-1, PF2: Pseudo-female-2, PF3: Pseudo-female-3, PM1: Pseudo-male-1, PM2: Pseudo-male-2, PM3: Pseudo-male-3.
Structure Analysis of P. sinensis Genders
The annotations of the selected reference genomes are often not accurate enough due to software and data limitations, making it necessary to optimize the structure of the original annotated transcripts. We used gffcompare to compare the annotated transcripts with the known transcripts of the genome, to find new transcripts, and to supplement transcript annotations. In total, we identified 5 new transcript types (O, J, X, I, and U) containing 16,077 transcripts (Figures 2A,B) Samples correlation matrix revealed that there was good agreement both between and within the sample groups (Supplementary Figure S3A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876045/fcell-10-876045-HTML/image_m/fcell-10-876045-g002.jpg)
FIGURE 2. Structure analysis of P. sinensis genders. (A) The structure of 5 new transcripts types. (B) The number of 5 new transcripts types. (C) The Wayne analysis of the number of long noncoding RNA transcripts identified by three software (CNC1, CPC2, and Pfam). (D) The distribution of the length of ploy(A) tails of different samples. (E) The number of alternative splice sites of females and pseudo-males. (F) The number of alternative splice sites of males and pseudo-females. Data were expressed as mean ± SD. *p < 0.05.
Long noncoding RNA is a type of RNA molecule whose transcript length exceeds 200 nt and does not encode a protein. We used three software (CNC1, CPC2, Pfam) to analyze the number of long noncoding RNA transcripts (Figure 2C). Differentially expressed long noncoding RNAs in M vs. PF and F vs. PM were enriched for many biological processes such as cellular process, single-organism process and biological regulation (Supplementary Figures S3B,C). Poly(A) tails are known regulators of mRNA stability and intracellular transport. Since the length of ploy(A) tails can affect the translation efficiency, the analysis of their length can uncover the dynamic processes of mRNA and protein translation. Based on the direct mRNA sequencing data, the global poly(A) profiles were analyzed, and length distributions were found to be different between the four gender populations but consistent within each group (Figure 2D). We next analyzed the alternative splicing types in each sample and counted the number of transcripts in each type. As shown in Figure 2E, the number of alternative splicing sites at positions such as alternative 3′ splice site (A3), alternative 5’ splice site (A5) and alternative first exons (AF) had significant differences between F and PM, which the number of alternative splicing sites was obviously reduced in the sex reversal process from F to PM. During the sex reversal process from M to PF, the number of alternative splicing sites at those positions was significantly increased (Figure 2F).
Differential Gene Expression Analysis Between Gender Types
To further evaluate the differences between the four P. sinensis genders, we conducted pairwise differential gene expression analyses. Approximately 15,000 differentially expressed genes were identified between M and F or M and PF, while 12,000 differential genes were found between female individuals and male phenotype individuals. The number of differentially expressed genes between the same gender (M and PM; F and PF) was small (Figure 3A). We divided these data into two groups (group I and group II) for analysis. Group I (i.e., F, M, PF) consisted of the differentially expressed genes related to the gender reversal process from male to pseudo-female, and group II (i.e., F, M, PM) consisted of the differential genes related to the F to PM transition. The numbers of differentially expressed genes in these two groups as determined by Wayne analysis are shown in Figures 3B,C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876045/fcell-10-876045-HTML/image_m/fcell-10-876045-g003.jpg)
FIGURE 3. Differential gene expression analysis between gender types. (A) The number of differential transcripts in pairs. (B) The Wayne analysis of differential genes among females, males and pseudo-females. (C) The Wayne analysis of differential genes among females, males and pseudo-males. (D) The cluster analysis of differential up-regulated genes of group I (i.e., F, M, and PF). (E) The cluster analysis of differential up-regulated genes of group II (i.e., F, M, and PM). (F) The cluster analysis of differential down-regulated genes of group I (i.e., F, M, and PF). (G) The cluster analysis of differential down-regulated genes of group II (i.e., F, M, and PM).
We performed cluster analysis on the differential genes in the two groups, and the up-regulated and down-regulated genes clusters identified for group I are shown in Figures 3D,F, while group II results are shown in 3E and 3G. All the differentially expressed genes were clustered into 10 categories, and among the up-regulated genes, those of the first category constituted the largest proportion, whereas the genes of the second category ranked first among the down-regulated differential genes in group I. GO enrichment analysis revealed that a large proportion of the up-regulated genes are involved in processes related to the endosome membrane, cellular calcium ion homeostasis and autophagosome (Supplementary Figure S4A). KEGG enrichment analysis showed that several signal pathways such as microRNA in cancer, protein processing in ER and influenza A were associated with the transition from M individuals to PF individuals (Supplementary Figure S4B). GO enrichment analysis of the second category among down-regulated differential genes showed processes related to motile cilium, axoneme and microtubule-based movement were significantly reduced during sex reversal (Supplementary Figure S4C). The down-regulated differential genes were mainly concentrated in the purine metabolism and Huntington disease signaling pathways (Supplementary Figure S4D). The differential genes of group II were likewise clustered into 10 categories, and cluster results indicated that the majority of up-regulated and down-regulated differential genes were in the first category (Figures 3E,G). GO and KEGG enrichment analysis results showed that the differential genes of the first category were mainly enriched in vesicle-mediated transport, motile cilium and multiple signaling pathways such as hippo signaling pathway and glycerophospholipid metabolism (Supplementary Figures S4E,F).
Methylation Analysis of P. sinensis Gender Types
The distribution and overall expression levels of transcripts with m6A modifications are shown in Figure 4A. The expression levels of the F phenotype were significantly higher than the expression level of M phenotype. Correlation analysis between polyA length and transcripts with m6A modification indicated a significant difference between female phenotype samples (F and PF) and male phenotype samples (M and PM; Figure 4B). We identified four different sequences of five bases each which are known m6A methylation sites, and compared the percentage of these sequences that harbored a m6A methylation. We found no significant difference between gender groups (Figure 4C). Then, we counted the total number of methylation modification sites including m5C and m6A, and found that the individuals induced by exogenous hormones (E2 and MT) had fewer methylation modification sites (Figure 4D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876045/fcell-10-876045-HTML/image_m/fcell-10-876045-g004.jpg)
FIGURE 4. Methylation analysis of P. sinensis gender types. (A) The distribution and overall expression levels of transcripts with m6A modifications of different samples. (B) The correlation analysis between poly(A) length and transcripts with m6A modification of different samples. (C) The percentage of different m6A methylation types of P. sinensis gender types. (D) The number of m5C and m6A modification sites of P. sinensis gender types. (E) The percentage of m5C methylation annotation in different positions. (F) The percentage of m6A methylation annotation in different positions. Data were expressed as mean ± SD. Different superscripts indicate significant difference (p < 0.05). ns, no significant difference.
Methylated gene regions (m5C and m6A) were annotated based on the reference genome of soft-shelled turtle. The results of anno-percent (Methylation annotation percentage) of m5C showed that the location of methylated modification sites had changed from CDS regions to 3′UTR regions during the process of sex reversal from male to pseudo-female. However, the location of many methylated modification sites (m5C) had changed from 3′UTR regions to CDS regions during female transition (Figure 4E). As with m5C modifications, the position of m6A modification sites also changed. Compared to female phenotype individuals, male phenotype individuals had more modification sites in the CDS regions and fewer modification sites in 3′UTR regions (Figure 4F). These results suggest that the PF is more similar to the F with respect to methylation modifications, while the PM is more similar to the M.
Analysis of m6A Modification Genes Across Gender Types
To further analyze the differentially expressed m6A modification genes among the four gender types, we conducted two sets of differential expression analyses: F vs. PM and PF vs. M (Figures 5A,B). In the first analysis, we identified 1,038 up-regulated and 1743 down-regulated m6A-modified genes, while in the second set we identified 1,279 up-regulated and 880 down-regulated m6A-modified genes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876045/fcell-10-876045-HTML/image_m/fcell-10-876045-g005.jpg)
FIGURE 5. Analysis of m6A modification genes across gender types. (A) The Volcano map of differential m6A methylation genes between PM and F. (B) The Volcano map of differential m6A methylation genes between PF and M. (C) The GO enrichment analysis of differential m6A methylation genes between PM and F. (D) The GO enrichment analysis of differential m6A methylation genes between PF and M. (E) The KEGG enrichment analysis of differential m6A methylation genes between PM and F. (F) The KEGG enrichment analysis of differential m6A methylation genes between PF and M.
GO enrichment analysis of the genes in the first set showed enrichment for the GO terms motile cilium, axoneme and mitochondrial outer membrane, indicating that these processes play an important role in the transition from F to PM (Figure 5C). In addition, KEGG enrichment analysis of the same set of genes identified several significant pathways such as phagosome, influenza A and protein processing in ER (Figure 5E). Likewise, GO and KEGG analyses of the differentially expressed genes in the PF vs. M set showed enrichment for many of the same GO terms, including motile cilium, axoneme and autophagosome (Figure 5D), as well as many of the same KEGG pathways, such as phagosome, influenza A and protein processing in ER (Figure 5F). These findings suggest that similar biological processes drive gender reversal in either direction.
The Chromosomal Location of Methylated Modification Sites
The genetic features of the 5 bp long m5C and m6A methylation sites were analyzed using the MEME Suite (Figures 6A,B). The results showed that the position of the second base in the m5C methylation site changed between male phenotype and female phenotype, while the m6A sequence features did not change (Figure 6B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876045/fcell-10-876045-HTML/image_m/fcell-10-876045-g006.jpg)
FIGURE 6. The chromosomal location of methylated modification sites. (A) The genetic features of the 5 bp long m5C methylation sites from male phenotype to female phenotype. (B) The genetic features of the 5 bp long m6A methylation sites from male phenotype to female phenotype. The abscissa represents the number of bases at the methylation site, the total height of each position is the sequence conservation of the bases, and the height of the base signal represents the relative frequency of the base at that position. (C) The chromosomal location of differential methylated gens both in m5C and m6A related sexual reversal process.
Given that the differential methylated genes seem to play a significant role in the transformation of male and female phenotypes, we analyzed and mapped some crucial genes to their chromosomal positions. Figure 6C shows the location of differential m5C and m6A modification genes that were identified in both sexual reversal pathways. These genes were mainly located on seven chromosomes (NW_005871042.1, NW_005870997.1, NW_005859001.1, NW_005855937.1, NW_005855804.1, NW_005854023.1, and NW_005853326.1), suggesting that these methylated chromosomes may be closely related to gender differentiation in both directions.
Verification of N6-Methylation Modification Sites of Sex-Related Genes
All the differentially expressed methylated genes located in these chromosomes are listed in Supplementary Table S3. We designed the primers corresponding to these differential genes and verified the expression of N6-methylated fragments. The methylated RNAs were captured by magnetic beads with m6A antibody and were reverse transcribed into cDNA for further analysis (Supplementary Table S4). The quantitative results showed that the N6-methylated levels of Odf2, Pacs2, Ak1 and Ube2o in the transition from M to PF were significantly reduced, while the N6-methylated levels of Odf2, Pacs2, Ak1 in the transition from F to PM increased and the N6-methylated level of Ube2o decreased (Figures 7A–F). The changes in RNA methylation of these genes can play a significant role in the sexual reversal process.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876045/fcell-10-876045-HTML/image_m/fcell-10-876045-g007.jpg)
FIGURE 7. Verification of N6-methylation modification sites of sex-related genes. (A–F) The mRNA expression level of Odf2, Pacs2-1, Pacs2-2, Pacs2-3, Ube2o, Ak1 related to actin among 4 types. Data were expressed as mean ± SD. Different superscripts indicate significant difference (p < 0.05). (G) The schematic of the changes of two sexual reversal processes.
A schematic summarizing the differences between the two sexual reversal ways is shown in Figure 7G. The schematic displays the most important changes in the sexual reversal process, including the number of different alternative splice sites, location of RNA methylations (m5C and m6A) and the N6-methylated levels of critical sex-related genes. Taken together, these results demonstrate that epigenetics plays an indispensable role in the sexual reversal process.
Discussion
Sex differentiation and differences in vertebrates have always attracted attention (Nagahama et al., 2021). Many animals, including many aquatic animals, show clear gender dimorphism (Dean and Mank, 2014), which has provided an important impetus for new approaches to aquatic breeding. In recent years, the new aquatic species of yellow catfish, Japanese flounder and tilapia were successfully bred based on gender dimorphism (Palaiokostas et al., 2015; Zhang et al., 2016). Chinese soft-shelled turtle is another economically important aquatic animal that displays a high degree of gender dimorphism (Zhang et al., 2018). Despite years of study, our understanding of the biological basis of sexual dimorphism of Chinese soft-shelled turtle remains incomplete. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that in addition to sex-changing fishes, gonochoristic vertebrates like mice maintain some degree of gonadal sexual plasticity into adulthood (Godwin, 2010; Li et al., 2019; Nagahama et al., 2021). Exogenous hormones such as estradiol and methyltestosterone can break the balance between estrogen and androgen, which can change the direction of gonadal differentiation (Balthazart et al., 2009; Carvalho et al., 2018; Toyota et al., 2020). In this study, we used exogenous hormones estradiol (E2) and methyltestosterone (MT) to stimulate the embryos of the Chinese soft-shelled turtle during the gonadal differentiation stage, and obtained PF and PM individuals.
Understanding sexual reversal requires insight into RNA transcription, processing, and modification (Parker et al., 2020). We employed nanopore direct RNA sequencing to examine the transcriptome of four gender types (F, M, PF, and PM) of Chinese soft-shelled turtle with a special focus on RNA methylation (m5C and m6A). Because of the long reads generated by nanopore RNA-seq (>1900 bp in this study), it can overcome some of the limitations of traditional next-generation RNA-seq, and thus reveal the complexity of RNA processing and modification in full-length single molecule reads.
Alternative splicing (AS) has emerged as a key event during gonad development in vertebrate (Domingos et al., 2018; Zhao et al., 2018). In this study, we showed that female turtles had more splicing events than male turtles, and that they were mainly concentrated in alternative 3′ splice sites (A3), alternative 5’ splice sites (A5) and alternative first exons (AF). Changes in alternative splicing events can cause evolutionary alterations in RNA binding proteins and related cis-elements (Ule and Blencowe, 2019). We found that the number of splicing events was significantly increased in female phenotype individuals, and that the related proteins and biological processes of the female pathway were consequently more complex. We speculated that these differential transcripts may have important regulatory effects on gonadal differentiation.
Epigenetics and genetics work together to affect the differentiation direction of the gonads in vertebrates (Gunes et al., 2016; Credendino et al., 2020; Ortega-Recalde et al., 2020). RNA methylation modification is an important type of epigenetics which has attracted much attention in recent years (Fu et al., 2019; Mongan et al., 2019). Methylation of the N (6) position of adenosine (m6A) is a posttranscriptional modification. The changes of RNA methylation sites may influence the progression of several diseases such as cancer (Meyer et al., 2012; Lin et al., 2019). During the process of sexual reversal, we found that the location of several methylation modification sites had changed. The methylation modification sites (m5C and m6A) on 3′UTRs were significantly increased in the transition from M individuals to PF individuals. By contrast, F individuals had strikingly fewer methylation sites on 3′UTRs compared to PM individuals. The changes of methylation sites revealed that the same genes may have different expression patterns in different gender types, even though their genotypes remain unchanged (Zhang et al., 2021). The transition of methylation position between CDS and 3′UTR proved that these changed genes are closely related to the differentiation of gonads.
Several studies have demonstrated an important role of epigenetics in the process of sexual reversal (Shao et al., 2014; Scaduto et al., 2017; Ortega-Recalde et al., 2020). For example, the X chromosome protects against bladder cancer in female mice via a KDM6A-dependent epigenetic mechanism (Kaneko and Li, 2018). Shao et al. found that epigenetic regulation plays multiple crucial roles in sexual reversal of tongue sole fish (Shao et al., 2014). However, the connection between RNA methylation and the sexual reversal process is rarely reported. In this study, we compared the RNA methylation levels of the m5C and m6A modifications between the four gender types of Chinese soft-shelled turtle, and identified ∼2,700 methylated genes that were differently expressed between F and PM and ∼2000 genes that were differentially expressed between M and PF. These genes were subsequently mapped to their chromosomal positions (Figure 6C, Supplementary Table S3). The quantitative results showed that the N6-methylated level of these differential methylated genes was significantly higher in male phenotypic individuals than in female phenotypic individuals (Figure 7). Taken together, our findings suggest that the direction of gender differentiation is determined by the degree of RNA methylation. Meanwhile, the sexual reversal induced by exogenous hormones such as estradiol and methyltestosterone will bring the changes of the location of RNA methylation. In summary, epigenetic regulation such as RNA methylation (m5C and m6A) and alternative splicing events play a key and irreplaceable role in the sexual reversal process of Chinese soft-shelled turtles.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: BioProject, accession number: PRJNA808173.
Ethics Statement
The animal study was reviewed and approved by the Animal Experimental Ethical Inspection of Laboratory Animal Centre of the Yangtze River Fisheries Research Institute, Chinese Academy of Fishery Sciences (Wuhan, China) (ID Number: 2022YFI-ZT-01).
Author Contributions
HL and TZ conceived the study. TZ conducted Chinese soft-shelled turtle cultivation and prepared the samples. GC performed the PCR and qPCR experiments. YW and MC performed RNA-seq data analyses. GZ provided the experimental consumables. HL and TZ prepared and revised the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by the National Freshwater Aquatic Germplasm Resource Center (FGRC18537) and Central Public-interest Scientific Institution Basal Research Fund, CAFS (No. 2020TD33).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.876045/full#supplementary-material
References
Balthazart, J., Cornil, C. A., Charlier, T. D., Taziaux, M., and Ball, G. F. (2009). Estradiol, a Key Endocrine Signal in the Sexual Differentiation and Activation of Reproductive Behavior in Quail. J. Exp. Zool. 311A (5), 323–345. doi:10.1002/jez.464
Bayega, A., Oikonomopoulos, S., Gregoriou, M.-E., Tsoumani, K. T., Giakountis, A., Wang, Y. C., et al. (2021). Nanopore Long-Read RNA-Seq and Absolute Quantification Delineate Transcription Dynamics in Early Embryo Development of an Insect Pest. Sci. Rep. 11 (1), 7878. doi:10.1038/s41598-021-86753-7
Byrne, A., Beaudin, A. E., Olsen, H. E., Jain, M., Cole, C., Palmer, T., et al. (2017). Nanopore Long-Read RNAseq Reveals Widespread Transcriptional Variation Among the Surface Receptors of Individual B Cells. Nat. Commun. 8, 16027. doi:10.1038/ncomms16027
Carvalho, L. B. d., Burusco, K. K., Jaime, C., Venâncio, T., Carvalho, A. F. S. d., Murgas, L. D. S., et al. (2018). Complexes between Methyltestosterone and β-cyclodextrin for Application in Aquaculture Production. Carbohydr. Polym. 179, 386–393. doi:10.1016/j.carbpol.2017.09.023
Credendino, S. C., Neumayer, C., and Cantone, I. (2020). Genetics and Epigenetics of Sex Bias: Insights from Human Cancer and Autoimmunity. Trends Genet. 36 (9), 650–663. doi:10.1016/j.tig.2020.06.016
Dean, R., and Mank, J. E. (2014). The Role of Sex Chromosomes in Sexual Dimorphism: Discordance between Molecular and Phenotypic Data. J. Evol. Biol. 27 (7), 1443–1453. doi:10.1111/jeb.12345
Domingos, J. A., Budd, A. M., Banh, Q. Q., Goldsbury, J. A., Zenger, K. R., and Jerry, D. R. (2018). Sex-specific Dmrt1 and Cyp19a1 Methylation and Alternative Splicing in Gonads of the Protandrous Hermaphrodite Barramundi. PLoS One 13 (9), e0204182. doi:10.1371/journal.pone.0204182
Dumschott, K., Schmidt, M. H.-W., Chawla, H. S., Snowdon, R., and Usadel, B. (2020). Oxford Nanopore Sequencing: New Opportunities for Plant Genomics? J. Exp. Bot. 71 (18), 5313–5322. doi:10.1093/jxb/eraa263
Fu, Y., Zorman, B., Sumazin, P., Sanna, P. P., and Repunte-Canonigo, V. (2019). Epitranscriptomics: Correlation of N6-Methyladenosine RNA Methylation and Pathway Dysregulation in the hippocampus of HIV Transgenic Rats. PLoS One 14 (1), e0203566. doi:10.1371/journal.pone.0203566
Garalde, D. R., Snell, E. A., Jachimowicz, D., Sipos, B., Lloyd, J. H., Bruce, M., et al. (2018). Highly Parallel Direct RNA Sequencing on an Array of Nanopores. Nat. Methods 15 (3), 201–206. doi:10.1038/nmeth.4577
Godwin, J. (2010). Neuroendocrinology of Sexual Plasticity in Teleost Fishes. Front. Neuroendocrinology 31 (2), 203–216. doi:10.1016/j.yfrne.2010.02.002
Green, M. R., and Sambrook, J. (2019). Removing DNA Contamination from RNA Samples by Treatment with RNase-free DNase I. Cold Spring Harb Protoc. 2019 (10), 101725. doi:10.1101/pdb.prot101725
Grünberger, F., Ferreira-Cerca, S., and Grohmann, D. (2022). Nanopore Sequencing of RNA and cDNA Molecules in Escherichia coli. RNA 28 (3), 400–417. doi:10.1261/rna.078937.121
Gu, Z., Eils, R., Schlesner, M., and Ishaque, N. (2018). EnrichedHeatmap: an R/Bioconductor Package for Comprehensive Visualization of Genomic Signal Associations. BMC Genomics 19 (1), 234. doi:10.1186/s12864-018-4625-x
Gunes, S. O., Metin Mahmutoglu, A., and Agarwal, A. (2016). Genetic and Epigenetic Effects in Sex Determination. Birth Defects Res. C Embryo Today 108 (4), 321–336. doi:10.1002/bdrc.21146
Jenjaroenpun, P., Wongsurawat, T., Pereira, R., Patumcharoenpol, P., Ussery, D. W., Nielsen, J., et al. (2018). Complete Genomic and Transcriptional Landscape Analysis Using Third-Generation Sequencing: a Case Study of Saccharomyces cerevisiae CEN.PK113-7D. Nucleic Acids Res. 46 (7), e38. doi:10.1093/nar/gky014
Kaliman, P. (2019). Epigenetics and Meditation. Curr. Opin. Psychol. 28, 76–80. doi:10.1016/j.copsyc.2018.11.010
Kaneko, S., and Li, X. (20182018). X Chromosome Protects against Bladder Cancer in Females via a KDM6A-dependent Epigenetic Mechanism. Sci. Adv. 4 (6), eaar5598. doi:10.1126/sciadv.aar5598
Karamitros, T., and Magiorkinis, G. (2018). Multiplexed Targeted Sequencing for Oxford Nanopore MinION: A Detailed Library Preparation Procedure. Methods Mol. Biol. 1712, 43–51. doi:10.1007/978-1-4939-7514-3_4
Kawagoshi, T., Uno, Y., Matsubara, K., Matsuda, Y., and Nishida, C. (2009). The ZW Micro-sex Chromosomes of the Chinese Soft-Shelled Turtle (Pelodiscus Sinensis, Trionychidae, Testudines) Have the Same Origin as Chicken Chromosome 15. Cytogenet. Genome Res. 125 (2), 125–131. doi:10.1159/000227837
Kawai, A., Nishida-Umehara, C., Ishijima, J., Tsuda, Y., Ota, H., and Matsuda, Y. (2007). Different Origins of Bird and Reptile Sex Chromosomes Inferred from Comparative Mapping of Chicken Z-Linked Genes. Cytogenet. Genome Res. 117 (1-4), 92–102. doi:10.1159/000103169
Keifer, J., and Zheng, Z. (2017). Cold Block of In Vitro Eyeblink Reflexes: Evidence Supporting the Use of Hypothermia as an Anesthetic in Pond Turtles. J. Exp. Biol. 220 (23), 4370–4373. doi:10.1242/jeb.168427
Kono, N., and Arakawa, K. (2019). Nanopore Sequencing: Review of Potential Applications in Functional Genomics. Develop. Growth Differ. 61 (5), 316–326. doi:10.1111/dgd.12608
Kukurba, K. R., and Montgomery, S. B. (2015). RNA Sequencing and Analysis. Cold Spring Harb Protoc. 2015 (11), pdb.top084970–969. doi:10.1101/pdb.top084970
Ladd-Acosta, C., and Fallin, M. D. (2016). The Role of Epigenetics in Genetic and Environmental Epidemiology. Epigenomics 8 (2), 271–283. doi:10.2217/epi.15.102
Li, M., Sun, L., and Wang, D. (2019). Roles of Estrogens in Fish Sexual Plasticity and Sex Differentiation. Gen. Comp. Endocrinol. 277, 9–16. doi:10.1016/j.ygcen.2018.11.015
Li, R., Ren, X., Ding, Q., Bi, Y., Xie, D., and Zhao, Z. (2020). Direct Full-Length RNA Sequencing Reveals Unexpected Transcriptome Complexity during Caenorhabditis elegans Development. Genome Res. 30 (2), 287–298. doi:10.1101/gr.251512.119
Liang, H., Meng, Y., Cao, L., Li, X., and Zou, G. (2019a). Effect of Exogenous Hormones on R-Spondin 1 (RSPO1) Gene Expression and Embryo Development in Pelodiscus Sinensis. Reprod. Fertil. Dev. 31 (9), 1425–1433. doi:10.1071/RD19045
Liang, H., Tong, M., Cao, L., Li, X., Li, Z., and Zou, G. (2018). Amino Acid and Fatty Acid Composition of Three Strains of Chinese Soft-Shelled Turtle (Pelodiscus Sinensis). Pjz 50 (3). doi:10.17582/journal.pjz/2018.50.3.1061.1069
Liang, H., Wang, L., Sha, H., and Zou, G. (2019b). Development and Validation of Sex-specific Markers in Pelodiscus Sinensis Using Restriction Site-Associated DNA Sequencing. Genes 10 (4), 302. doi:10.3390/genes10040302
Lin, X., Chai, G., Wu, Y., Li, J., Chen, F., Liu, J., et al. (2019). RNA m6A Methylation Regulates the Epithelial Mesenchymal Transition of Cancer Cells and Translation of Snail. Nat. Commun. 10 (1), 2065. doi:10.1038/s41467-019-09865-9
Liu, R.-J., Long, T., Li, J., Li, H., and Wang, E.-D. (2017). Structural Basis for Substrate Binding and Catalytic Mechanism of a Human RNA:m5C Methyltransferase NSun6. Nucleic Acids Res. 45 (11), 6684–6697. doi:10.1093/nar/gkx473
Loman, N. J., and Watson, M. (2015). Successful Test Launch for Nanopore Sequencing. Nat. Methods 12 (4), 303–304. doi:10.1038/nmeth.3327
Lorenz, D. A., Sathe, S., Einstein, J. M., and Yeo, G. W. (2020). Direct RNA Sequencing Enables m6A Detection in Endogenous Transcript Isoforms at Base-specific Resolution. RNA 26 (1), 19–28. doi:10.1261/rna.072785.119
Lu, H., Giordano, F., and Ning, Z. (2016). Oxford Nanopore minION Sequencing and Genome Assembly. Genomics, Proteomics & Bioinformatics 14 (5), 265–279. doi:10.1016/j.gpb.2016.05.004
Lubuta, P., Workman, M., Kerkhoven, E. J., and Workman, C. T. (2019). Investigating the Influence of Glycerol on the Utilization of Glucose in Yarrowia Lipolytica Using RNA-Seq-Based Transcriptomics. G3 (Bethesda) 9 (12), 4059–4071. doi:10.1534/g3.119.400469
Ma, X., Wang, G., Liu, H., Wu, L., Tian, X., Wang, L., et al. (2021). Characterization of Aromatase in Testes and Potential Functional Analysis with Letrozole Administration in Male Chinese Soft-Shelled Turtle (Pelodiscus Sinensis). Aquaculture 544, 737089. doi:10.1016/j.aquaculture.2021.737089
Meyer, K. D., Saletore, Y., Zumbo, P., Elemento, O., Mason, C. E., and Jaffrey, S. R. (2012). Comprehensive Analysis of mRNA Methylation Reveals Enrichment in 3′ UTRs and Near Stop Codons. Cell 149 (7), 1635–1646. doi:10.1016/j.cell.2012.05.003
Mongan, N. P., Emes, R. D., and Archer, N. (2019). Detection and Analysis of RNA Methylation. F1000Res 8, 559. doi:10.12688/f1000research.17956.1
Nagahama, Y., Chakraborty, T., Paul-Prasanth, B., Ohta, K., and Nakamura, M. (2021). Sex Determination, Gonadal Sex Differentiation, and Plasticity in Vertebrate Species. Physiol. Rev. 101 (3), 1237–1308. doi:10.1152/physrev.00044.2019
Ortega-Recalde, O., Goikoetxea, A., Hore, T. A., Todd, E. V., and Gemmell, N. J. (2020). The Genetics and Epigenetics of Sex Change in Fish. Annu. Rev. Anim. Biosci. 8, 47–69. doi:10.1146/annurev-animal-021419-083634
Othman, R., Ron, X. J., Yao, H., O'Bryant, P., Rapp, D., Pei, J. C., et al. (2022). The Effect of Methyltestosterone (MT) on Sex Differentiation and Growth in Juvenile Yellow Perch (Perca flavescens). Fish. Physiol. Biochem. 48, 161–171. doi:10.1007/s10695-021-01038-0
Palaiokostas, C., Bekaert, M., Khan, M. G., Taggart, J. B., Gharbi, K., McAndrew, B. J., et al. (2015). A Novel Sex-Determining QTL in Nile tilapia (Oreochromis niloticus). BMC Genomics 16, 171. doi:10.1186/s12864-015-1383-x
Parker, M. T., Knop, K., Sherwood, A. V., Schurch, N. J., Mackinnon, K., Gould, P. D., et al. (2020). Nanopore Direct RNA Sequencing Maps the Complexity of Arabidopsis mRNA Processing and m6A Modification. Elife 9. doi:10.7554/eLife.49658
Roach, N. P., Sadowski, N., Alessi, A. F., Timp, W., Taylor, J., and Kim, J. K. (2020). The Full-Length Transcriptome of C. elegans Using Direct RNA Sequencing. Genome Res. 30 (2), 299–312. doi:10.1101/gr.251314.119
Sánchez-Vásquez, E., Alata Jimenez, N., Vázquez, N. A., and Strobl-Mazzulla, P. H. (2018). Emerging Role of Dynamic RNA Modifications during Animal Development. Mech. Develop. 154, 24–32. doi:10.1016/j.mod.2018.04.002
Scaduto, C. M., Kabrawala, S., Thomson, G. J., Scheving, W., Ly, A., Anderson, M. Z., et al. (2017). Epigenetic Control of Pheromone MAPK Signaling Determines Sexual Fecundity in Candida Albicans. Proc. Natl. Acad. Sci. U.S.A. 114 (52), 13780–13785. doi:10.1073/pnas.1711141115
Schmittgen, T. D., and Livak, K. J. (2008). Analyzing Real-Time PCR Data by the Comparative C(T) Method. Nat. Protoc. 3 (1001), 1101–1108. doi:10.1038/nprot.2008.73
Schulz, R. W., Bogerd, J., Male, R., Ball, J., Fenske, M., Olsen, L. C., et al. (2007). Estrogen-induced Alterations in Amh and Dmrt1 Expression Signal for Disruption in Male Sexual Development in the Zebrafish. Environ. Sci. Technol. 41, 6305–6310. doi:10.1021/es070785+
Shao, C., Li, Q., Chen, S., Zhang, P., Lian, J., Hu, Q., et al. (2014). Epigenetic Modification and Inheritance in Sexual Reversal of Fish. Genome Res. 24 (4), 604–615. doi:10.1101/gr.162172.113
Toyota, K., Masuda, S., Sugita, S., Miyaoku, K., Yamagishi, G., Akashi, H., et al. (2020). Estrogen Receptor 1 (ESR1) Agonist Induces Ovarian Differentiation and Aberrant Müllerian Duct Development in the Chinese Soft-Shelled Turtle, Pelodiscus Sinensi. Zool Stud. 59, e54. doi:10.6620/ZS.2020.59-54
Traube, F. R., and Carell, T. (2017). The Chemistries and Consequences of DNA and RNA Methylation and Demethylation. RNA Biol. 14 (9), 1099–1107. doi:10.1080/15476286.2017.1318241
Ule, J., and Blencowe, B. J. (2019). Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution. Mol. Cel 76 (2), 329–345. doi:10.1016/j.molcel.2019.09.017
Wang, L., Wu, Z., Zou, C., Liang, S., Zou, Y., Liu, Y., et al. (2020). Sex-dependent RNA Editing and N6-Adenosine RNA Methylation Profiling in the Gonads of a Fish, the Olive Flounder (Paralichthys olivaceus). Front. Cel Dev. Biol. 8, 751. doi:10.3389/fcell.2020.00751
Wang, Y., and Wang, Z. (2015). Efficient Backsplicing Produces Translatable Circular mRNAs. RNA 21 (2), 172–179. doi:10.1261/rna.048272.114
Wu, Y., Sarkissyan, M., and Vadgama, J. V. (2015). Epigenetics in Breast and Prostate Cancer. Methods Mol. Biol. 1238, 425–466. doi:10.1007/978-1-4939-1804-1_23
Yang, X., Yang, Y., Sun, B.-F., Chen, Y.-S., Xu, J.-W., Lai, W.-Y., et al. (2017). 5-methylcytosine Promotes mRNA export - NSUN2 as the Methyltransferase and ALYREF as an m5C Reader. Cell Res 27 (5), 606–625. doi:10.1038/cr.2017.55
Yang, Y., Fan, X., Mao, M., Song, X., Wu, P., Zhang, Y., et al. (2017). Extensive Translation of Circular RNAs Driven by N6-Methyladenosine. Cel Res 27 (5), 626–641. doi:10.1038/cr.2017.31
Yu, T., Qi, X., Zhang, L., Ning, W., Gao, D., Xu, T., et al. (2021). Dynamic Reprogramming and Function of RNA N6-Methyladenosine Modification during Porcine Early Embryonic Development. Zygote 29 (6), 417–426. doi:10.1017/s0967199420000799
Zhang, J., Ma, W., He, Y., Wu, J., Dawar, F. U., Ren, F., et al. (2016). Sex Biased Expression of Ghrelin and GHSR Associated with Sexual Size Dimorphism in Yellow Catfish. Gene 578 (2), 169–176. doi:10.1016/j.gene.2015.12.017
Zhang, K., Zhuang, X., Dong, Z., Xu, K., Chen, X., Liu, F., et al. (2021). The Dynamics of N6-Methyladenine RNA Modification in Interactions between rice and Plant Viruses. Genome Biol. 22 (1), 189. doi:10.1186/s13059-021-02410-2
Zhang, Z., Song, R., Xing, X., Wang, L., and Niu, C. (2018). Division of Chinese Soft‐shelled Turtle Intestine with Molecular Markers Is Slightly Different from the Morphological and Histological Observation. Integr. Zoolog. 13 (1), 112–121. doi:10.1111/1749-4877.12261
Zhao, L., Wang, C., Lehman, M. L., He, M., An, J., Svingen, T., et al. (2018). Transcriptomic Analysis of mRNA Expression and Alternative Splicing during Mouse Sex Determination. Mol. Cell Endocrinol. 478, 84–96. doi:10.1016/j.mce.2018.07.010
Keywords: RNA methylation, nanopore sequencing, epigenetics, chinese soft-shelled turtle, sexual reversal process
Citation: Zhou T, Chen G, Chen M, Wang Y, Zou G and Liang H (2022) Direct Full-Length RNA Sequencing Reveals an Important Role of Epigenetics During Sexual Reversal in Chinese Soft-Shelled Turtle. Front. Cell Dev. Biol. 10:876045. doi: 10.3389/fcell.2022.876045
Received: 15 February 2022; Accepted: 10 March 2022;
Published: 25 March 2022.
Edited by:
Douglas Mark Ruden, Wayne State University, United StatesCopyright © 2022 Zhou, Chen, Chen, Wang, Zou and Liang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guiwei Zou, zougw@yfi.ac.cn; Hongwei Liang, lianghw@yfi.ac.cn