Skip to main content
Download PDF

The m6A methylation regulates gonadal sex differentiation in chicken embryo

Journal of Animal Science and Biotechnology volume 13, Article number: 52 (2022) Cite this article

Abstract

Background

As a ubiquitous reversible epigenetic RNA modification, N6-methyladenosine (m6A) plays crucial regulatory roles in multiple biological pathways. However, its functional mechanisms in sex determination and differentiation during gonadal development of chicken embryos are not clear. Therefore, we established a transcriptome-wide m6A map in the female and male chicken left gonads of embryonic day 7 (E7) by methylated RNA immunoprecipitation sequencing (MeRIP-seq) to offer insight into the landscape of m6A methylation and investigate the post-transcriptional modification underlying gonadal differentiation.

Results

The chicken embryonic gonadal transcriptome was extensively methylated. We found 15,191 and 16,111 m6A peaks in the female and male left gonads, respectively, which were mainly enriched in the coding sequence (CDS) and stop codon. Among these m6A peaks, we identified that 1013 and 751 were hypermethylated in females and males, respectively. These differential peaks covered 281 and 327 genes, such as BMP2, SMAD2, SOX9 and CYP19A1, which were primarily associated with development, morphogenesis and sex differentiation by functional enrichment. Further analysis revealed that the m6A methylation level was positively correlated with gene expression abundance. Furthermore, we found that YTHDC2 could regulate the expression of sex-related genes, especially HEMGN and SOX9, in male mesonephros/gonad mingle cells, which was verified by in vitro experiments, suggesting a regulatory role of m6A methylation in chicken gonad differentiation.

Conclusions

This work provided a comprehensive m6A methylation profile of chicken embryonic gonads and revealed YTHDC2 as a key regulator responsible for sex differentiation. Our results contribute to a better understanding of epigenetic factors involved in chicken sex determination and differentiation and to promoting the future development of sex manipulation in poultry industry.

Background

Sex manipulation technologies are of great significance to the production performance and economic benefit of livestock [1]. During the production of laying hens, at least seven billion day-old male chicks are culled every year worldwide, which raises cost and animal welfare concerns [2,3,4]. Research on the mechanisms of sex determination and differentiation in chickens, especially during the embryonic stage, will contribute to achieving early sex manipulation of chickens and be beneficial for chicken breeding and resource protection strategies [5]. In addition, chicken embryonic gonadal differentiation is also an excellent model for studying key factors of vertebrate sex determination [6] and human sexual development disorders [7, 8]. Although many investigations have focused on chicken sex determination and differentiation [9,10,11,12,13,14,15,16,17,18,19,20], the mechanisms underlying gonadal differentiation are still elusive.

The sex determination and differentiation of chickens during the embryonic stage is how the gonad develops into the testis or ovary. The process is well-known to be affected by the regulation of genetic factors and hormone levels [10, 12, 18, 20]. Unlike mammals, homogametic males (ZZ) produce testes, and heterogametic females (ZW) develop ovaries in birds. Considering that no female sex-determining genes have been found on the W chromosome [21], multiple studies have recognized that the Z chromosome gene doublesex and mab-3–related transcription factor 1 (DMRT1) is the key factor for the formation of testes [12, 16, 20]. The key conserved Sertoli cell differentiation factor (SOX9) may be induced and activated by high levels of DMRT1 in males (ZZ) and initiate testicular differentiation during embryonic development [16, 20, 22]. Due to the low dosage of DMRT1, the FOXL2 signaling pathway in females is activated [16, 20]. Estrogen is a key factor in determining female sex, and the blockade of estrogen could directly cause female embryonic gonad masculinization [9,10,11, 18, 20]. Although several genes responsible for sex determination and differentiation were identified, multiple pieces of evidence indicate that epigenetic modification is involved in this process [23,24,25,26,27,28,29,30]. The noncoding RNA transcribed from the male hypermethylated region (MHM) has been suggested to be associated with chicken gonadal development and differentiation [24, 25, 28,29,30]. Temporarily feminized gonads maintain transcriptomic and epigenetic memories of genetic sex [23]. Therefore, in-depth exploration from epigenetic aspects would provide novel insights into sex differentiation.

N6-methyladenosine (m6A) is an important posttranscriptional mRNA modification discovered in the 1970s and mainly regulates mRNA metabolism, including RNA stability, translation, RNA splicing, and transport [31,32,33,34,35,36,37,38]. In mammals, multiple studies have shown that m6A modification plays essential roles during spermatogenesis [39,40,41,42,43]. Knockout of Mettl3, Mettl14 and Ythdc2 in mice inhibited gonadal development and resulted in infertility [39,40,41]. Evidence has illustrated that m6A methylation is also involved in oocyte meiosis in Xenopus laevis [44]. Ythdf2 and Ythdf3 deletion prevented female gonad formation in zebrafish [45]. In addition, the dynamic m6A modification is abundant during follicle development in pigs [46] or follicle selection in chickens [47]. Several studies have shown that m6A methylation involves numerous biological processes [39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54]. However, knowledge about the distribution and function of m6A in chicken sex determination and differentiation of gonads is limited.

To decipher the role of RNA modification in the gonadal differentiation of chicken embryos, we performed global MeRIP-seq in female and male left gonads on E7. The present work generated a high-resolution m6A methylation profile and explored the potential molecular mechanisms underlying sex differentiation. We found that many differential m6A methylation peaks existed between the female and male left gonads in chicken embryos, and most genes with differential m6A methylation were involved mainly sex differentiation and development. It should be noted that YTHDC2 can influence the expression of sex differentiation-related genes in males, which was confirmed by in vitro experiments. Our results provide distinctive insights into the epigenetic mechanism of chicken embryonic gonadal differentiation and will benefit further investigation of sex manipulation in poultry industry and research on human disorders of sexual development.

Methods and materials

Ethics statement

The experiments were approved by the Animal Welfare Committee of China Agricultural University and performed in accordance with the protocol outlined in the “Guide for Care and Use of Laboratory Animals” (China Agricultural University).

Embryo tissue collection and RNA isolation

Fertilized eggs were obtained from a pure line of White Leghorns raised in the Experimental Base of Poultry Genetic Resources and Breeding, College of Animal Science and Technology, China Agricultural University. The eggs were incubated in an automated egg incubator at 37.8 °C and 60% relative humidity with rotation every 2 h. Once the eggs reached E7, only the left gonads of the chicken embryo were collected, immediately put into the RNAlater™ Stabilization Solution (Invitrogen, Carlsbad, CA, USA), and stored at − 20 °C. The remaining embryo tissues were used to determine sex by direct PCR kit (TransGen Biotech, Beijing, China) with CHD1 primers (forward-5′-GTTACTGATTCGTCTACGAGA-3′, reverse-5′-TTGAAATGATCCAGTGCTTG-3′) [55]. Due to the RNA sample from a single individual gonad is not enough for MeRIP-seq, more than 1200 embryos were collected, and around 200 left gonads were mixed into a pool according to sex. Subsequently, six pools from three pools of female left gonads and three pools of males were prepared for RNA isolation. Total RNA was extracted with TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions. The RNA concentration and purity were measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The integrity of the RNA was determined using an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA).

MeRIP-seq

Poly(A) RNA was purified from 40 μg of total RNA using Dynabeads Oligo (dT) and then fragmented into 100-nucleotide-long fragments using Magnesium RNA Fragmentation Module (Illumina, Inc., CA, United States) at 86 °C for 7 min. Then, the fragmented mRNAs were incubated for 2 h at 4 °C with anti-m6A polyclonal antibody (Synaptic Systems, Goettingen, Germany) in IP buffer. Eluted m6A-containing fragments (IP) and untreated input control fragments were then concentrated to generate the final cDNA library. The libraries were qualified and absolutely quantified using an Agilent Bioanalyzer 2100 (Agilent Technologies, CA, USA). The prepared libraries were then sequenced on an Illumina NovaSeq 6000 (150 bp paired-end, PE150).

The analysis of sequencing data

The clean reads were aligned to the chicken genome sequences (GRCg6a) with HISAT2 [56]. SAM files were converted to the BAM format using Samtools [57], and PCR duplicates were removed using the Picard MarkDuplicates option to generate filtered BAM files. The m6A modification peaks were called using MACS2 [58] (−g 9.6e8 --nomodel --extsize 200 -q 0.05). Meanwhile, input data (high-throughput paired-end RNA sequencing, RNA-seq) were used as background. Referring to the proportion of the effective genome sizes of humans and mice, we adjusted the parameters of -g in chickens. The number of reads falling in the m6A peak in each sample was counted with bedtools [59]. Moreover, the putative peaks were annotated, and the motifs enriched in peak regions were analyzed using RNAmod [60]. The Htseq-count [61] tool was used to count the gene-level reads. DEGs were identified with the R package DESeq2 [62] between female and male gonads. The genes/peaks with fold change > 2 or < 0.5 and FDR < 0.05 were considered different.

Functional annotation

Using BioMart [63], we identified homologs of chicken DEGs in humans. Functional analysis of these homologs was performed using the Metascape online tool [64]. The Gene Ontology (GO) terms for biological process, cellular component, and molecular function categories, as well as Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, were enriched based on the Metascape online tool with default parameters.

Cell separation and culture

The tissues (mesonephros and gonads) were harvested and placed in PBS when the fertile chicken eggs reached E7. Samples were dispersed by incubation with trypsin 0.25% (GIBCO, Grand Island, NY, USA) at 37 °C for 10 min with constant shaking. The digestions were stopped by the addition of culture medium containing DMEM/F12 with 10% (v/v) fetal bovine serum and 1% penicillin/streptomycin (GIBCO, Grand Island, NY, USA), and the samples were passed through a cell strainer (40 μm) and collected into 50-mL tubes. Subsequently, the cells were washed with PBS and seeded in cell culture plates. Sex identification was carried out by PCR as described above. All the cells from individual males were mixed.

Small interfering RNA assays and qRT-PCR

The design and synthesis of YTHDC2-siRNA primer was performed by RiboBio (Guangzhou, China). For the YTHDC2-knockdown assay, the following sequence was used: 5′-CTCACAGATACCAAGTAT-3′. We applied YTHDC2-siRNA (RiboBio, Guangzhou, China) according to Fugene HD (Promega, Madison, WI, USA) transfection into male cells. After transfection for 48 h, the cells were collected for RNA extraction, and cDNA was synthesized using the First-Strand Synthesis kit (Takara, Japan) and PrimeScript™ RT reagent kit with gDNA Eraser (Takara, Japan). Gonadal cDNA was also synthesized as described above. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed with an ABI 7500 system (Applied Biosystems, Bedford, MA, USA) using the TB Green® Premix Ex Taq™ Kit (Takara, Japan) according to the manufacturer’s instructions. Primer sequences for β-actin, AMH, SOX9, FOXL2 and CYP19A1 have been published by our previous results [65]. The primer of YTHDC2, FTO, METTL3, ALKBH5, DMRT1 and HEMGN were designed by DNAMAN version 6.0 software. β-actin was used as the internal control, and the sequences of the gene-specific primers are listed in Table 1.

Table 1 Primers used in qRT-PCR
Full size table

Statistical analysis

In this study, statistical analyses were performed with SPSS 22.0 (SPSS, Chicago, IL, USA). The data are expressed as the mean ± SD (standard deviation) and were analyzed using a two-tailed Student’s t-test, and at least three replicates were conducted in multiple independent experiments. The differences were considered to be statistically significant at a P-value < 0.05.

Results

General features of m6A methylation in chicken embryonic gonads

Our previous RNA-seq data indicated that YTHDC2 showed differential expression levels between females and males in early chicken embryonic gonads (Additional file 1: Fig. S1a). Therefore, we hypothesized that m6A methylation may play significant roles in the process of sex determination and differentiation. To gain insight into this phenomenon, we collected the left gonads of both sexes at E7 in triplicate for MeRIP-seq. The mapping statistics of MeRIP-seq from all samples are displayed in Additional file 2: Table S1, suggesting that each sample showed high quality for the following analyses. We identified 15,191 and 16,111 m6A peaks in the female and male left gonads in chicken embryos, respectively. Overall, 12,515 m6A peaks overlapped between the two groups, along with 2676 and 3596 unique m6A peaks in female and male left gonads (Fig. 1a), indicating that m6A methylation in the male gonadal gene is relatively richer. Consistent with previous studies [66], gonadal m6A methylation peaks were generally found in the genomic features, in which exon regions accounted for the largest proportion (Additional file 1: Fig. S1b).

Fig. 1
figure 1

Overview of the m6A methylation profile in chicken embryonic gonads. a The overlap of m6A peaks between female and male gonads. b Distribution of m6A peaks along transcripts. Transcripts are grouped into 5′ UTR, CDS, and 3′ UTR methylation based on the identified m6A sites. c Analysis of gonadal m6A peaks in different genomic features. d Percentage of m6A methylated transcripts with different numbers of m6A peaks. e, f The top motifs enriched across m6A peaks identified from female and male gonads

Full size image

To further refine the topological pattern of m6A modification in chicken embryonic gonads, we investigated the distribution of m6A peaks. The results demonstrated that the peaks were markedly enriched in the CDS and stop codon (Fig. 1b, c). Based on the number of m6A peaks contained in each transcript, we further classified the methylated transcripts into five types and found that approximately 90% of transcripts contained one or two m6A peaks (Fig. 1d). Overall, there was high consistency in the distribution of m6A peaks between female and male chicken embryonic left gonads. Given that the motifs that bind to various transcription factors are the sites where RNA methylation and demethylation begin, we performed a comprehensive motif scan analysis on male and female m6A peaks, respectively. We found that GGACU was significantly enriched in the female and male chicken embryonic left gonads (Fig. 1e, f), which is in agreement with previous studies [67, 68].

Differential m6A methylation analysis

Three replicates of MeRIP-seq in the female and male groups showed high concordance using the Pearson correlation coefficient (Fig. 2a). A principal component analysis (PCA) plot based on the m6A methylation level displayed a clear separation between female and male left gonads (Fig. 2b), in which PC1 explained 87.5% of the phenotypic variance. To detect the differences in m6A methylation levels between female and male gonads, we assessed the differentially methylated m6A peaks (DMPs) marked by MeRIP signals. A total of 1764 DMPs were detected, including 1013 female- and 751 male-biased DMPs in chicken embryonic left gonads (Fig. 2c), which corresponded to 281 and 327 protein-coding genes (DMGs; Additional file 2: Table S2). The heatmap revealed that these sex-biased DMPs showed highly clear sex-specific patterns at E7 (Fig. 2c). The top 20 female- and male-biased DMPs are listed in Additional file 2: Table S3. The annotation results of DMPs on the chromosome distribution indicated that 56.3%, 1.7%, 39.0% and 3.1% female-biased peaks and 72.7%, 1.1%, 0.5% and 25.7% male-biased peaks were located on the autosome, unknown chromosome, W chromosome and Z chromosome (Additional file 1: Fig. S1c). More DMPs were located on Z chromosome in males than in females.

Fig. 2
figure 2

Analysis of differential m6A methylation in chicken embryonic gonads. a A heatmap of the sample correlation matrix of m6A showing high similarities between duplicates and dissimilarities between female and male left gonads. b PCA plot of MeRIP-seq of female and male left gonads. c Bandplots (top) and heatmaps (bottom) showing the quantification of MeRIP-seq data of female-biased and male-biased m6A peaks in gonads at E7. d Top 20 significantly enriched terms of female-biased (left) and male-biased (right) DMGs

Full size image

The distribution of DMGs on chromosomes was similar between female and male gonads (Additional file 1: Fig. S1d). Among these DMGs, we found that some genes were related to sex determination and differentiation, such as FOXL2, CYP19A1, AMH and SOX9. In addition, many development-related genes were also discovered, such as BMP2, BMPR2, FGF2 and SMAD2. To better understand the functional consequences of m6A methylation, we performed a functional enrichment analysis of sex-biased DMGs using Metascape software. The top 20 significantly enriched gene ontology (GO) terms were identified (Fig. 2d). The significant terms of the female-biased DMPs were related to ERK1 and ERK2 cascade, muscle tissue development, regulation of growth, receptor ligand activity, response to growth factor, sex differentiation, tissue morphogenesis, SMAD protein signal transduction and embryonic morphogenesis. Moreover, the terms driven by the male-biased DMPs were enriched mainly in cell junction organization, cellular component morphogenesis, urogenital system development, tissue morphogenesis, positive regulation of cell migration, MAPK cascade, smooth muscle cell migration, outflow tract morphogenesis, regulation of cell morphogenesis.

Correlation analysis of m6A methylation and DEGs

To clarify the functional consequences of the gene expression levels modified by m6A methylation, we profiled the global transcriptomic landscape of both gonads on the left using RNA-seq. The results of sequencing alignment of RNA-seq are listed in Additional file 2: Table S1. Similar to the m6A results, unsupervised hierarchical clustering of the top 1000 most variable genes revealed a distinct expression signature in female and male left gonads (Fig. 3a). RNA-seq analysis detected 902 differentially expressed genes (DEGs), including 457 female- and 445 male-biased DEGs (Fig. 3b). The majority of DEGs resided on the autosome (Fig. 3c). Subsequently, the global relationship between gene expression and m6A methylation was calculated. By assigning m6A methylation regions to the nearest genes, we found that most genes showed a positive correlation between the expression levels and the abundance of m6A peaks based on calculated fold changes, and only a few of them were negatively correlated (Fig. 3d). Moreover, we found 410 DEGs with differential m6A methylation levels. A total of 203 and 206 genes with m6A hypermethylation showed upregulated transcription levels in females and males, respectively, and one gene with m6A hypomethylation was downregulated in females (Fig. 3e). Interestingly, many sex-related genes showed differences in both m6A methylation and mRNA expression (Fig. 3f, Additional file 1: Fig. S1e-h, Fig. 4c). For instance, SOX9, a male gonadal developmental marker gene, showed significantly higher m6A methylation and mRNA levels in the male gonads (Figs. 3f, 4c), while CYP19A1, which is associated with female gonadal differentiation, exhibited stronger m6A methylation and mRNA expression levels in the female gonads (Additional file 1: Fig. S1h, Fig. 4c). Additionally, the genes of HEMGN, AMH and FOXL2 all have DMPs between female and male gonads (Additional file 1: Fig. S1e-g); however, we did not find DMPs in the DMRT1 gene body.

Fig. 3
figure 3

Integration analysis of MeRIP-seq and RNA-seq. a Unsupervised clustering analysis showing the expression profiles of the top 1000 most variable genes in the female and male gonads at E7. b Volcano plots of DEGs number between female and male left gonads. c The number and percentage of DEGs in each chromosomal allocation between female and male gonads. d Four-quadrant graph exhibiting the genes containing methylated m6A peaks in chicken embryonic gonads. e Venn diagrams show the shared and unique genes obtained from four groups: female DMGs, male DMGs, female DEGs and male DEGs. f The abundance of m6A peaks in the SOX9 gene of female and male gonads detected by MeRIP-seq

Full size image
Fig. 4
figure 4

Comprehensive analysis of gene expression. a qRT-PCR results of m6A methylation-related genes in chicken embryonic gonads. b The expression of YTHDC2 at E7 in chicken embryonic gonads. The y axis shows the rlog-normalized counts. c qRT-PCR results of sex differentiation-related genes in chicken embryonic gonads. d Changes in the relative expression of genes related to male gonadal markers in cells with YTHDC2 knockdown. * and ** indicate P < 0.05 and P < 0.01, respectively

Full size image

To explore the regulatory mechanism of m6A methylation in chicken gonadal differentiation, we performed a series of in vitro experiments for some differential m6A methylation-related genes and genes involved in gonadal differentiation. We found that the expression level of YTHDC2 was significantly higher in male gonads (Fig. 4a), which is consistent with the RNA-seq results (Fig. 4b). Some marker genes of sex differentiation were also differentially expressed (Fig. 4c). Moreover, we found that the stop codon of SOX9 and the CDS of HEMGN showed significant enrichment in IP samples (Fig. 3f; Additional file 1: Fig. S1e). Based on these results, we hypothesize that sex-related genes are affected by YTHDC2. To clarify the process, we used YTHDC2-siRNA to downregulate the expression of YTHDC2 in male mesonephros/gonad mingle cells. Interestingly, the expression of the sex-related genes SOX9 and HEMGN was downregulated in the YTHDC2-siRNA group (Fig. 4d). Based on these observations, we preliminarily concluded that YTHDC2 could affect sex differentiation-related gene expression by regulating epigenetic patterns in chicken embryonic gonads. Therefore, m6A methylation may play a role in sex determination and differentiation in chicken.

Discussion

The chicken embryo is a superior model in the study of cell biology and organ development [6, 69,70,71,72,73,74]. In the past few decades, the regulatory mechanisms for chicken sex determination and differentiation, such as transcription factors [12, 16, 20], DNA methylation [24, 25, 28,29,30] and hormones [9, 10, 18, 20, 23, 27], have been extensively studied. Recently, m6A methylation, a reversible apparent modification [31], has been reported to play an essential role in vertebrate physiology and reproductive processes like development, spermatogenesis, oogenesis and infertility [39,40,41,42, 44,45,46,47,48,49,50]. Notably, m6A methylation is essential in the sex determination and gonadal development process of many animals, including mice [39,40,41], fish [45, 75] and Drosophila [51,52,53]. However, the roles and molecular mechanisms of m6A modification in the sex differentiation of chicken embryos remain unclear. Here, we established a comprehensive m6A methylation profile in chicken embryonic gonads. Our results suggested that m6A methylation is significantly different between female and male gonads, and results in distinct expression patterns of many sex-related genes. Further, these marker genes involved in gonadal differentiation were confirmed to be regulated by m6A enzyme YTHDC2, which means that m6A methylation may be an important factor affecting chicken sex differentiation.

We found more m6A peaks in the male gonads than in the female gonads, implying that m6A methylation might be a vital impetus in sex determination and differentiation. The genomic distribution of m6A peaks in male and female embryonic gonads is globally similar. The percentage of m6A peaks located in the CDS was the largest, while the coverage density of m6A peaks was the highest at the stop codon in chicken embryonic gonads, which was consistent with research on mammals [67, 68, 76, 77]. Subsequently, we compared the DMPs between female and male gonads in chicken embryos. The quantitative advantage of DMPs in females is reflected mainly in the sex chromosomes, while the number of DMGs is equal between males and females, which may be caused by the poor annotation of the W chromosome [14]. Previous studies have also proven that m6A methylation plays an important role in the process of tissue development [78,79,80,81,82]. In this study, GO analysis showed that DMGs were enriched mainly in sex differentiation- and development-relevant items, including sex differentiation, SMAD protein signal transduction, urogenital system development and tissue morphogenesis, which suggests that DMGs were related to the dynamic process of gonadal development. A number of studies have reported the functions of these DMGs, like BMP2 [83], CYP19A1 [27], ESR1 [18], FOXL2 [6], JUN [84], SMAD2 [84, 85] and SOX9 [16, 20, 22]. Our findings revealed that the post-transcriptional modifications of these genes should be alternative factor driving sex difference, which means that m6A methylation may play a major role in the differentiation of gonads. This finding would contribute to understanding the mechanisms of sex determination and differentiation in chickens and providing a promising reference for further research of human gonad development related diseases.

The role and underlying mechanism of m6A modification in regulatory of gene expression is still uncharted territory. Collectively, m6A has been linked to reduce mRNA stability and promote mRNA degradation in various biological pathways [32, 86], suggesting that the m6A methylation level is negatively correlated with gene expression. In addition, several researches suggested a possible positive relationship between the extent of m6A methylation and the mRNA levels [87,88,89,90]. Integration analysis of RNA-seq and MeRIP-seq data in our work found that the expression level of most genes was positively correlated with the m6A methylation signal in the whole genome, which is also supported by many previous studies [87,88,89,90]. These findings mean that the roles of m6A methylation in transcriptional regulation need to be elucidated in future studies. The mRNA m6A modifications are well-known to be recognized and bound by m6A reader proteins, such as YTHDF1–3, YTHDC1 and YTHDC2 [36, 39, 45, 86, 91, 92]. In our study, RNA-seq analysis suggested that only YTHDC2 was differentially expressed between female and male gonads among m6A methylation-related genes, which was verified by qRT-PCR. Previous studies have shown that the loss of YTHDC2 downregulates the meiotic genes in mitotic spermatogonia [42]. The function of YTHDC2 can reduce the translation efficiency of target genes and reduce the mRNA abundance in the meiosis of germline cells [39]. Interestingly, our in vitro experiments of YTHDC2 knockdown in males suggested that YTHDC2 could downregulate HEMGN and SOX9 gene expression by reshaping their m6A patterns. HEMGN has been shown to be involved in sex determination in early chicken embryos [15]. SOX9 is a marker gene of testicular differentiation during chicken embryonic development [16, 20, 22]. Nearly half of the DEGs showed differential m6A methylation modifications in the gene body, which were involved in gonadal differentiation and development terms. The expression of sex differentiation marker genes was also verified by qRT-PCR. In chicken, the changes of these candidate regulators of gonadal sex differentiation induced sex reversal [9,10,11,12, 15,16,17]. Thus, m6A methylation may play roles during sex differentiation in chickens. However, further research at the single-cell level is required to validate this hypothesis and will be helpful to elucidate molecular basis underlying gonadal differentiation. Here, we provided a pioneering work demonstrating the function of YTHDC2 in chicken embryonic gonads. Although the roles of YTHDC2 have been studied in mammals [39, 42, 43, 54, 76, 93, 94], the precise transcriptional regulation molecular mechanism underlying multiple crucial roles of YTHDC2 remains to be determined in chicken. More evidence is needed to explore whether YTHDC2 can directly regulate the expression of sex-related genes.

Conclusions

This study first analyzed transcriptome-wide m6A methylation modification pattern in the gonads of chicken embryos and preliminarily explored roles of YTHDC2 in regulating key genes underlying sex differentiation. The m6A methylation profile and molecular basis underlying gonadal differentiation in chicken embryo provide a new avenue for studying RNA modification in vertebrate sex determination and differentiation and offer distinctive insights into the epigenetics mechanism of studies about human gonadal development and infertility.

Availability of data and materials

The sequencing datasets in this study are available at NCBI [PRJNA766306].

Abbreviations

CDS:

Coding sequence

DEGs:

Differentially expressed genes

DMGs:

Differential methylation related coding genes

DMPs:

Differential methylation m6A peaks

E7:

Embryonic Day 7

GO:

Gene ontology

m6A:

N6-methyladenosine

MeRIP-seq:

Methylated RNA immunoprecipitation sequencing

PCA:

Principal component analysis

RNA-seq:

High-throughput paired-end RNA sequencing

References

  1. Xie Y, Xu Z, Wu Z, Hong L. Sex manipulation technologies progress in livestock: a review. Front Vet Sci. 2020;7:481. https://doi.org/10.3389/fvets.2020.00481.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Gremmen B. Moral dilemmas of animal production systems. Anim Front. 2020;10(1):15–20. https://doi.org/10.1093/af/vfz051.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Krautwald-Junghanns ME, Cramer K, Fischer B, Förster A, Galli R, Kremer F, et al. Current approaches to avoid the culling of day-old male chicks in the layer industry, with special reference to spectroscopic methods. Poult Sci. 2018;97(3):749–57. https://doi.org/10.3382/ps/pex389.

    Article  CAS  PubMed  Google Scholar 

  4. de Haas EN, Oliemans E, van Gerwen M. The need for an alternative to culling day-old male layer chicks: a survey on awareness, alternatives, and the willingness to pay for alternatives in a selected population of dutch citizens. Front Vet Sci. 2021;8:662197. https://doi.org/10.3389/fvets.2021.662197.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Nakamura Y, Kagami H, Tagami T. Development, differentiation and manipulation of chicken germ cells. Develop Growth Differ. 2013;55(1):20–40. https://doi.org/10.1111/dgd.12026.

    Article  Google Scholar 

  6. Hirst CE, Major AT, Smith CA. Sex determination and gonadal sex differentiation in the chicken model. Int J Dev Biol. 2018;62(1–2-3):153–66. https://doi.org/10.1387/ijdb.170319cs.

    Article  CAS  PubMed  Google Scholar 

  7. Witchel SF. Disorders of sex development. Best Pract Res Clin Obstet Gynaecol. 2018;48:90–102. https://doi.org/10.1016/j.bpobgyn.2017.11.005.

    Article  PubMed  Google Scholar 

  8. Hiort O, Birnbaum W, Marshall L, Wünsch L, Werner R, Schröder T, et al. Management of disorders of sex development. Nat Rev Endocrinol. 2014;10(9):520–9. https://doi.org/10.1038/nrendo.2014.108.

    Article  PubMed  Google Scholar 

  9. Lambeth LS, Cummins DM, Doran TJ, Sinclair AH, Smith CA. Overexpression of aromatase alone is sufficient for ovarian development in genetically male chicken embryos. PLoS One. 2013;8(6):e68362. https://doi.org/10.1371/journal.pone.0068362.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Elbrecht A, Smith RG. Aromatase enzyme activity and sex determination in chickens. Science. 1992;255(5043):467–70. https://doi.org/10.1126/science.1734525.

    Article  CAS  PubMed  Google Scholar 

  11. Smith CA, Katz M, Sinclair AH. DMRT1 is upregulated in the gonads during female-to-male sex reversal in ZW chicken embryos. Biol Reprod. 2003;68(2):560–70. https://doi.org/10.1095/biolreprod.102.007294.

    Article  CAS  PubMed  Google Scholar 

  12. Smith CA, Roeszler KN, Ohnesorg T, Cummins DM, Farlie PG, Doran TJ, et al. The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. Nature. 2009;461(7261):267–71. https://doi.org/10.1038/nature08298.

    Article  CAS  PubMed  Google Scholar 

  13. Zhao D, McBride D, Nandi S, McQueen HA, McGrew MJ, Hocking PM, et al. Somatic sex identity is cell autonomous in the chicken. Nature. 2010;464(7286):237–42. https://doi.org/10.1038/nature08852.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Ayers KL, Davidson NM, Demiyah D, Roeszler KN, Grützner F, Sinclair AH, et al. RNA sequencing reveals sexually dimorphic gene expression before gonadal differentiation in chicken and allows comprehensive annotation of the W-chromosome. Genome Biol. 2013;14(3):R26. https://doi.org/10.1186/gb-2013-14-3-r26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Nakata T, Ishiguro M, Aduma N, Izumi H, Kuroiwa A. Chicken hemogen homolog is involved in the chicken-specific sex-determining mechanism. Proc Natl Acad Sci U S A. 2013;110(9):3417–22. https://doi.org/10.1073/pnas.1218714110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Lambeth LS, Raymond CS, Roeszler KN, Kuroiwa A, Nakata T, Zarkower D, et al. Over-expression of DMRT1 induces the male pathway in embryonic chicken gonads. Dev Biol. 2014;389(2):160–72. https://doi.org/10.1016/j.ydbio.2014.02.012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Lambeth LS, Morris K, Ayers KL, Wise TG, O'Neil T, Wilson S, et al. Overexpression of anti-Müllerian hormone disrupts gonadal sex differentiation, blocks sex hormone synthesis, and supports cell autonomous sex development in the chicken. Endocrinology. 2016;157(3):1258–75. https://doi.org/10.1210/en.2015-1571.

    Article  CAS  PubMed  Google Scholar 

  18. Guioli S, Zhao D, Nandi S, Clinton M, Lovell-Badge R. Oestrogen in the chick embryo can induce chromosomally male ZZ left gonad epithelial cells to form an ovarian cortex that can support oogenesis. Development. 2020;147(4):dev181693. https://doi.org/10.1242/dev.181693.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Estermann MA, Williams S, Hirst CE, Roly ZY, Serralbo O, Adhikari D, et al. Insights into gonadal sex differentiation provided by single-cell transcriptomics in the chicken embryo. Cell Rep. 2020;31(1):107491. https://doi.org/10.1016/j.celrep.2020.03.055.

    Article  CAS  PubMed  Google Scholar 

  20. Ioannidis J, Taylor G, Zhao D, Liu L, Idoko-Akoh A, Gong D, et al. Primary sex determination in birds depends on DMRT1 dosage, but gonadal sex does not determine adult secondary sex characteristics. Proc Natl Acad Sci U S A. 2021;118(10):e2020909118. https://doi.org/10.1073/pnas.2020909118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Hirst CE, Major AT, Ayers KL, Brown RJ, Mariette M, Sackton TB, et al. Sex reversal and comparative data undermine the W chromosome and support Z-linked DMRT1 as the regulator of gonadal sex differentiation in birds. Endocrinology. 2017;158(9):2970–87. https://doi.org/10.1210/en.2017-00316.

    Article  CAS  PubMed  Google Scholar 

  22. Jakob S, Lovell-Badge R. Sex determination and the control of Sox9 expression in mammals. FEBS J. 2011;278(7):1002–9. https://doi.org/10.1111/j.1742-4658.2011.08029.x.

    Article  CAS  PubMed  Google Scholar 

  23. Shioda K, Odajima J, Kobayashi M, Kobayashi M, Cordazzo B, Isselbacher KJ, et al. Transcriptomic and epigenetic preservation of genetic sex identity in estrogen-feminized male chicken embryonic gonads. Endocrinology. 2021;162(1):bqaa208. https://doi.org/10.1210/endocr/bqaa208.

    Article  PubMed  Google Scholar 

  24. Sun D, Maney DL, Layman TS, Chatterjee P, Yi SV. Regional epigenetic differentiation of the Z chromosome between sexes in a female heterogametic system. Genome Res. 2019;29(10):1673–84. https://doi.org/10.1101/gr.248641.119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Roeszler KN, Itman C, Sinclair AH, Smith CA. The long non-coding RNA, MHM, plays a role in chicken embryonic development, including gonadogenesis. Dev Biol. 2012;366(2):317–26. https://doi.org/10.1016/j.ydbio.2012.03.025.

    Article  CAS  PubMed  Google Scholar 

  26. Miao N, Wang X, Hou Y, Feng Y, Gong Y. Identification of male-biased microRNA-107 as a direct regulator for nuclear receptor subfamily 5 group a member 1 based on sexually dimorphic microRNA expression profiling from chicken embryonic gonads. Mol Cell Endocrinol. 2016;429:29–40. https://doi.org/10.1016/j.mce.2016.03.033.

    Article  CAS  PubMed  Google Scholar 

  27. Ellis HL, Shioda K, Rosenthal NF, Coser KR, Shioda T. Masculine epigenetic sex marks of the CYP19A1/aromatase promoter in genetically male chicken embryonic gonads are resistant to estrogen-induced phenotypic sex conversion. Biol Reprod. 2012;87(1):23, 1–12. https://doi.org/10.1095/biolreprod.112.099747.

  28. Yang X, Deng J, Zheng J, Xia L, Yang Z, Qu L, et al. A window of MHM demethylation correlates with key events in gonadal differentiation in the chicken. Sex Dev. 2016;10(3):152–8. https://doi.org/10.1159/000447659.

    Article  CAS  PubMed  Google Scholar 

  29. Teranishi M, Shimada Y, Hori T, Nakabayashi O, Kikuchi T, Macleod T, et al. Transcripts of the MHM region on the chicken Z chromosome accumulate as non-coding RNA in the nucleus of female cells adjacent to the DMRT1 locus. Chromosom Res. 2001;9(2):147–65. https://doi.org/10.1023/a:1009235120741.

    Article  CAS  Google Scholar 

  30. Caetano LC, Gennaro FG, Coelho K, Araújo FM, Vila RA, Araújo A, et al. Differential expression of the MHM region and of sex-determining-related genes during gonadal development in chicken embryos. Genet Mol Res. 2014;13(1):838–49. https://doi.org/10.4238/2014.February.13.2.

    Article  CAS  PubMed  Google Scholar 

  31. Desrosiers R, Friderici K, Rottman F. Identification of methylated nucleosides in messenger RNA from Novikoff hepatoma cells. Proc Natl Acad Sci U S A. 1974;71(10):3971–5. https://doi.org/10.1073/pnas.71.10.3971.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18(1):31–42. https://doi.org/10.1038/nrm.2016.132.

    Article  CAS  PubMed  Google Scholar 

  33. Wei CM, Moss B. Nucleotide sequences at the N6-methyladenosine sites of Hela cell messenger ribonucleic acid. Biochemistry. 1977;16(8):1672–6. https://doi.org/10.1021/bi00627a023.

    Article  CAS  PubMed  Google Scholar 

  34. Zhou J, Wan J, Gao X, Zhang X, Jaffrey SR, Qian SB. Dynamic m(6) a mRNA methylation directs translational control of heat shock response. Nature. 2015;526(7574):591–4. https://doi.org/10.1038/nature15377.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Roundtree IA, Evans ME, Pan T, He C. Dynamic RNA modifications in gene expression regulation. Cell. 2017;169(7):1187–200. https://doi.org/10.1016/j.cell.2017.05.045.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Xiao W, Adhikari S, Dahal U, Chen YS, Hao YJ, Sun BF, et al. Nuclear m(6) a reader YTHDC1 regulates mRNA splicing. Mol Cell. 2016;61(4):507–19. https://doi.org/10.1016/j.molcel.2016.01.012.

    Article  CAS  PubMed  Google Scholar 

  37. Hsu PJ, Shi H, Zhu AC, Lu Z, Miller N, Edens BM, et al. The RNA-binding protein FMRP facilitates the nuclear export of N (6)-methyladenosine-containing mRNAs. J Biol Chem. 2019;294(52):19889–95. https://doi.org/10.1074/jbc.AC119.010078.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Roundtree IA, Luo GZ, Zhang Z, Wang X, Zhou T, Cui Y, et al. YTHDC1 mediates nuclear export of N(6)-methyladenosine methylated mRNAs. Elife. 2017;6:e31311. https://doi.org/10.7554/eLife.31311.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Hsu PJ, Zhu Y, Ma H, Guo Y, Shi X, Liu Y, et al. Ythdc2 is an N(6)-methyladenosine binding protein that regulates mammalian spermatogenesis. Cell Res. 2017;27(9):1115–27. https://doi.org/10.1038/cr.2017.99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Xu K, Yang Y, Feng GH, Sun BF, Chen JQ, Li YF, et al. Mettl3-mediated m(6) a regulates spermatogonial differentiation and meiosis initiation. Cell Res. 2017;27(9):1100–14. https://doi.org/10.1038/cr.2017.100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Lin Z, Hsu PJ, Xing X, Fang J, Lu Z, Zou Q, et al. Mettl3−/Mettl14-mediated mRNA N(6)-methyladenosine modulates murine spermatogenesis. Cell Res. 2017;27(10):1216–30. https://doi.org/10.1038/cr.2017.117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Wojtas MN, Pandey RR, Mendel M, Homolka D, Sachidanandam R, Pillai RS. Regulation of m(6) a transcripts by the 3′→5′ RNA helicase ythdc2 is essential for a successful meiotic program in the mammalian germline. Mol Cell. 2017;68(2):374–87.e12. https://doi.org/10.1016/j.molcel.2017.09.021.

    Article  CAS  PubMed  Google Scholar 

  43. Bailey AS, Batista PJ, Gold RS, Chen YG, de Rooij DG, Chang HY, et al. The conserved RNA helicase YTHDC2 regulates the transition from proliferation to differentiation in the germline. Elife. 2017;6:e26116. https://doi.org/10.7554/eLife.26116.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Qi ST, Ma JY, Wang ZB, Guo L, Hou Y, Sun QY. N6-methyladenosine sequencing highlights the involvement of mRNA methylation in oocyte meiotic maturation and embryo development by regulating translation in Xenopus laevis. J Biol Chem. 2016;291(44):23020–6. https://doi.org/10.1074/jbc.M116.748889.

    Article  CAS  PubMed  Google Scholar 

  45. Kontur C, Jeong M, Cifuentes D, Giraldez AJ. Ythdf m(6) a readers function redundantly during zebrafish development. Cell Rep. 2020;33(13):108598. https://doi.org/10.1016/j.celrep.2020.108598.

    Article  CAS  PubMed  Google Scholar 

  46. Cao Z, Zhang D, Wang Y, Tong X, Avalos LFC, Khan IM, et al. Identification and functional annotation of m6a methylation modification in granulosa cells during antral follicle development in pigs. Anim Reprod Sci. 2020;219:106510. https://doi.org/10.1016/j.anireprosci.2020.106510.

    Article  CAS  PubMed  Google Scholar 

  47. Fan Y, Zhang C, Zhu G. Profiling of RNA N6-methyladenosine methylation during follicle selection in chicken ovary. Poult Sci. 2019;98(11):6117–24. https://doi.org/10.3382/ps/pez277.

    Article  CAS  PubMed  Google Scholar 

  48. Zheng G, Dahl JA, Niu Y, Fedorcsak P, Huang CM, Li CJ, et al. ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility. Mol Cell. 2013;49(1):18–29. https://doi.org/10.1016/j.molcel.2012.10.015.

    Article  CAS  PubMed  Google Scholar 

  49. Zhao T, Wang J, Wu Y, Han L, Chen J, Wei Y, et al. Increased m6a modification of RNA methylation related to the inhibition of demethylase FTO contributes to MEHP-induced Leydig cell injury. Environ Pollut. 2021;268(Pt A):115627. https://doi.org/10.1016/j.envpol.2020.115627.

    Article  CAS  PubMed  Google Scholar 

  50. Ding C, Zou Q, Ding J, Ling M, Wang W, Li H, et al. Increased N6-methyladenosine causes infertility is associated with FTO expression. J Cell Physiol. 2018;233(9):7055–66. https://doi.org/10.1002/jcp.26507.

    Article  CAS  PubMed  Google Scholar 

  51. Lence T, Akhtar J, Bayer M, Schmid K, Spindler L, Ho CH, et al. M(6) a modulates neuronal functions and sex determination in Drosophila. Nature. 2016;540(7632):242–7. https://doi.org/10.1038/nature20568.

    Article  CAS  PubMed  Google Scholar 

  52. Kan L, Grozhik AV, Vedanayagam J, Patil DP, Pang N, Lim KS, et al. The m(6) a pathway facilitates sex determination in Drosophila. Nat Commun. 2017;8(1):15737. https://doi.org/10.1038/ncomms15737.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Haussmann IU, Bodi Z, Sanchez-Moran E, Mongan NP, Archer N, Fray RG, et al. M(6) a potentiates Sxl alternative pre-mRNA splicing for robust Drosophila sex determination. Nature. 2016;540(7632):301–4. https://doi.org/10.1038/nature20577.

    Article  CAS  PubMed  Google Scholar 

  54. Jain D, Puno MR, Meydan C, Lailler N, Mason CE, Lima CD, et al. Ketu mutant mice uncover an essential meiotic function for the ancient RNA helicase YTHDC2. Elife. 2018;7:e30919. https://doi.org/10.7554/eLife.30919.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Fridolfsson A-K, Ellegren H. A simple and universal method for molecular sexing of non-ratite birds. J Avian Biol. 1999;30(1):116–21. https://doi.org/10.2307/3677252.

    Article  Google Scholar 

  56. Kim D, Langmead B, Salzberg SL. Hisat: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. https://doi.org/10.1038/nmeth.3317.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and samtools. Bioinformatics. 2009;25(16):2078–9. https://doi.org/10.1093/bioinformatics/btp352.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137. https://doi.org/10.1186/gb-2008-9-9-r137.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Quinlan AR. Bedtools: the swiss-army tool for genome feature analysis. Curr Protoc Bioinform. 2014;47:11.2.1–34. https://doi.org/10.1002/0471250953.bi1112s47.

    Article  Google Scholar 

  60. Liu Q, Gregory RI. RNAmod: an integrated system for the annotation of mRNA modifications. Nucleic Acids Res. 2019;47(W1):W548–55. https://doi.org/10.1093/nar/gkz479.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. https://doi.org/10.1093/bioinformatics/btu638.

    Article  CAS  PubMed  Google Scholar 

  62. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. https://doi.org/10.1186/gb-2010-11-10-r106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Haw RA, Croft D, Yung CK, Ndegwa N, D'Eustachio P, Hermjakob H, et al. The reactome biomart. Database (Oxford). 2011;2011:bar031. https://doi.org/10.1093/database/bar031.

    Article  CAS  Google Scholar 

  64. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523. https://doi.org/10.1038/s41467-019-09234-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Yang X, Zheng J, Xu G, Qu L, Chen S, Li J, et al. Exogenous cMHM regulates the expression of DMRT1 and ER alpha in avian testes. Mol Biol Rep. 2010;37(4):1841–7. https://doi.org/10.1007/s11033-009-9619-y.

    Article  CAS  PubMed  Google Scholar 

  66. Wang Y, Zheng Y, Guo D, Zhang X, Guo S, Hui T, et al. M6a methylation analysis of differentially expressed genes in skin tissues of coarse and fine type Liaoning cashmere goats. Front Genet. 2019;10:1318. https://doi.org/10.3389/fgene.2019.01318.

    Article  CAS  PubMed  Google Scholar 

  67. Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, et al. Topology of the human and mouse m6a RNA methylomes revealed by m6a-seq. Nature. 2012;485(7397):201–6. https://doi.org/10.1038/nature11112.

    Article  CAS  PubMed  Google Scholar 

  68. Meyer KD, Saletore Y, Zumbo P, Elemento O, Mason CE, Jaffrey SR. Comprehensive analysis of mRNA methylation reveals enrichment in 3′ UTRs and near stop codons. Cell. 2012;149(7):1635–46. https://doi.org/10.1016/j.cell.2012.05.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Kain KH, Miller JW, Jones-Paris CR, Thomason RT, Lewis JD, Bader DM, et al. The chick embryo as an expanding experimental model for cancer and cardiovascular research. Dev Dyn. 2014;243(2):216–28. https://doi.org/10.1002/dvdy.24093.

    Article  PubMed  Google Scholar 

  70. Vilches-Moure JG. Embryonic chicken (Gallus gallus domesticus) as a model of cardiac biology and development. Comp Med. 2019;69(3):184–203. https://doi.org/10.30802/aalas-cm-18-000061.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Vergara MN, Canto-Soler MV. Rediscovering the chick embryo as a model to study retinal development. Neural Dev. 2012;7(1):22. https://doi.org/10.1186/1749-8104-7-22.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Stern CD. The chick; a great model system becomes even greater. Dev Cell. 2005;8(1):9–17. https://doi.org/10.1016/j.devcel.2004.11.018.

    Article  CAS  PubMed  Google Scholar 

  73. Flentke GR, Smith SM. The avian embryo as a model for fetal alcohol spectrum disorder. Biochem Cell Biol. 2018;96(2):98–106. https://doi.org/10.1139/bcb-2017-0205.

    Article  CAS  PubMed  Google Scholar 

  74. Fauzia E, Barbhuyan TK, Shrivastava AK, Kumar M, Garg P, Khan MA, et al. Chick embryo: a preclinical model for understanding ischemia-reperfusion mechanism. Front Pharmacol. 2018;9:1034. https://doi.org/10.3389/fphar.2018.01034.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Wang L, Wu Z, Zou C, Liang S, Zou Y, Liu Y, et al. Sex-dependent RNA editing and N6-adenosine RNA methylation profiling in the gonads of a fish, the olive flounder (Paralichthys olivaceus). Front Cell Dev Biol. 2020;8:751. https://doi.org/10.3389/fcell.2020.00751.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Mao Y, Dong L, Liu XM, Guo J, Ma H, Shen B, et al. M(6) a in mRNA coding regions promotes translation via the RNA helicase-containing YTHDC2. Nat Commun. 2019;10(1):5332. https://doi.org/10.1038/s41467-019-13317-9.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Zhang H, Shi X, Huang T, Zhao X, Chen W, Gu N, et al. Dynamic landscape and evolution of m6a methylation in human. Nucleic Acids Res. 2020;48(11):6251–64. https://doi.org/10.1093/nar/gkaa347.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Cheng B, Leng L, Li Z, Wang W, Jing Y, Li Y, et al. Profiling of RNA N (6) -methyladenosine methylation reveals the critical role of m(6) a in chicken adipose deposition. Front Cell Dev Biol. 2021;9:590468. https://doi.org/10.3389/fcell.2021.590468.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Xi L, Carroll T, Matos I, Luo JD, Polak L, Pasolli HA, et al. M6a RNA methylation impacts fate choices during skin morphogenesis. Elife. 2020;9:e56980. https://doi.org/10.7554/eLife.56980.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Wang Y, Gao M, Zhu F, Li X, Yang Y, Yan Q, et al. METTL3 is essential for postnatal development of brown adipose tissue and energy expenditure in mice. Nat Commun. 2020;11(1):1648. https://doi.org/10.1038/s41467-020-15488-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Kmietczyk V, Riechert E, Kalinski L, Boileau E, Malovrh E, Malone B, et al. M(6)a-mRNA methylation regulates cardiac gene expression and cellular growth. Life Sci Alliance. 2019;2(2):e201800233. https://doi.org/10.26508/lsa.201800233.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Wang X, Wu R, Liu Y, Zhao Y, Bi Z, Yao Y, et al. M(6) a mRNA methylation controls autophagy and adipogenesis by targeting Atg5 and Atg7. Autophagy. 2020;16(7):1221–35. https://doi.org/10.1080/15548627.2019.1659617.

    Article  CAS  PubMed  Google Scholar 

  83. Kashimada K, Pelosi E, Chen H, Schlessinger D, Wilhelm D, Koopman P. FOXL2 and BMP2 act cooperatively to regulate follistatin gene expression during ovarian development. Endocrinology. 2011;152(1):272–80. https://doi.org/10.1210/en.2010-0636.

    Article  CAS  PubMed  Google Scholar 

  84. Zhang M, Xu P, Sun X, Zhang C, Shi X, Li J, et al. Jun promotes chicken female differentiation by inhibiting Smad2. Cytotechnology. 2021;73(1):101–13. https://doi.org/10.1007/s10616-020-00447-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Whyte J, Glover JD, Woodcock M, Brzeszczynska J, Taylor L, Sherman A, et al. FGF, insulin, and SMAD signaling cooperate for avian primordial germ cell self-renewal. Stem Cell Reports. 2015;5(6):1171–82. https://doi.org/10.1016/j.stemcr.2015.10.008.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Wang X, Lu Z, Gomez A, Hon GC, Yue Y, Han D, et al. N6-methyladenosine-dependent regulation of messenger RNA stability. Nature. 2014;505(7481):117–20. https://doi.org/10.1038/nature12730.

    Article  CAS  PubMed  Google Scholar 

  87. Yue B, Song C, Yang L, Cui R, Cheng X, Zhang Z, et al. METTL3-mediated N6-methyladenosine modification is critical for epithelial-mesenchymal transition and metastasis of gastric cancer. Mol Cancer. 2019;18(1):142. https://doi.org/10.1186/s12943-019-1065-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Wang YN, Jin HZ. Transcriptome-wide m(6) a methylation in skin lesions from patients with psoriasis vulgaris. Front Cell Dev Biol. 2020;8:591629. https://doi.org/10.3389/fcell.2020.591629.

    Article  PubMed  PubMed Central  Google Scholar 

  89. Liu L, Wang J, Sun G, Wu Q, Ma J, Zhang X, et al. M(6) a mRNA methylation regulates CTNNB1 to promote the proliferation of hepatoblastoma. Mol Cancer. 2019;18(1):188. https://doi.org/10.1186/s12943-019-1119-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Li B, Wang X, Li Z, Lu C, Zhang Q, Chang L, et al. Transcriptome-wide analysis of N6-methyladenosine uncovers its regulatory role in gene expression in the lepidopteran Bombyx mori. Insect Mol Biol. 2019;28(5):703–15. https://doi.org/10.1111/imb.12584.

    Article  CAS  PubMed  Google Scholar 

  91. Zhu T, Roundtree IA, Wang P, Wang X, Wang L, Sun C, et al. Crystal structure of the YTH domain of YTHDF2 reveals mechanism for recognition of N6-methyladenosine. Cell Res. 2014;24(12):1493–6. https://doi.org/10.1038/cr.2014.152.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Wang X, Zhao BS, Roundtree IA, Lu Z, Han D, Ma H, et al. N(6)-methyladenosine modulates messenger RNA translation efficiency. Cell. 2015;161(6):1388–99. https://doi.org/10.1016/j.cell.2015.05.014.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Abby E, Tourpin S, Ribeiro J, Daniel K, Messiaen S, Moison D, et al. Implementation of meiosis prophase I programme requires a conserved retinoid-independent stabilizer of meiotic transcripts. Nat Commun. 2016;7(1):10324. https://doi.org/10.1038/ncomms10324.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Soh YQS, Mikedis MM, Kojima M, Godfrey AK, de Rooij DG, Page DC. Meioc maintains an extended meiotic prophase I in mice. PLoS Genet. 2017;13(4):e1006704. https://doi.org/10.1371/journal.pgen.1006704.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We thank the researchers in our laboratory for their assistance in samples collection.

Funding

This work was funded in part by grants from China Agricultural Research System (CARS-40).

Author information

Authors and Affiliations

  1. National Engineering Laboratory for Animal Breeding and Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, China Agricultural University, Beijing, China

    Jianbo Li, Xiuan Zhang, Xiqiong Wang, Congjiao Sun, Jiangxia Zheng, Junying Li & Ning Yang

  2. Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China

    Guoqiang Yi

Authors
  1. Jianbo Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Xiuan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xiqiong Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Congjiao Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jiangxia Zheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Junying Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Guoqiang Yi
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Ning Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

NY and GY conceived and designed the project. JbL, XZ, XW, CS, JZ and JyL performed the experiments and collected samples. JbL, and GY analyzed the data and wrote the manuscript. NY revised the manuscript. All authors read and approved the final draft.

Corresponding authors

Correspondence to Guoqiang Yi or Ning Yang.

Ethics declarations

Ethics approval and consent to participate

All experiments on chicken embryos were performed in accordance with the protocol outlined in the “Guide for Care and Use of Laboratory Animals” (China Agricultural University).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Supplementary Information

Additional file 1: Fig. S1.

Characteristics of MeRIP-seq in chicken gonads. a Transcriptional change in YTHDC2 between the female and male gonads at four developmental stages: (embryonic Day 4.5, E4.5), E5.5, E7 and E10. b The distribution of sequencing reads on the genome. c-d Percentage of DMPs (c) and DMGs (d) in each chromosomal allocation between female and male gonads. e-h The abundance of m6A peaks in the HEMGN, AMH, FOXL2 and CYP19A1 genes of female and male gonads detected by MeRIP-seq.

Additional file 2: Table S1.

Summary of sequence data and read alignment statistics. Table S2. Summary of m6A peaks and annotation genes. Table S3. Summary of the female-biased and male-biased DMPs in chicken embryonic gonads.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, J., Zhang, X., Wang, X. et al. The m6A methylation regulates gonadal sex differentiation in chicken embryo. J Animal Sci Biotechnol 13, 52 (2022). https://doi.org/10.1186/s40104-022-00710-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40104-022-00710-6

Keywords

Journal of Animal Science and Biotechnology

ISSN: 2049-1891