rbz Revista Brasileira de Zootecnia R. Bras. Zootec. 1516-3598 1806-9290 Sociedade Brasileira de Zootecnia 00605 10.37496/rbz5220210149 Breeding and genetics Expression analysis of m6A-related genes in various tissues of Meishan pigs at different developmental stages 0000-0001-8386-2738 Cao Yanan Conceptualization. Formal analysis. Methodology. Visualization. Writing – original draft 1 0000-0003-3539-1526 Zhang Shuoshuo Formal analysis. Methodology. Visualization 1 0000-0002-7795-3016 Wang Guangzheng Formal analysis. Methodology 1 0000-0001-8507-3492 Zhang Shuai Writing – original draft 1 0000-0001-5486-2253 Bao Wenbin Writing – review & editing 1 2 0000-0002-0954-6312 Wu Shenglong Conceptualization. Funding acquisition. Supervision. Writing – review & editing 1 2 * Yangzhou University College of Animal Science and Technology Key Laboratory for Animal Genetics, Breeding, Reproduction and Molecular Design Yangzhou China Yangzhou University, College of Animal Science and Technology, Key Laboratory for Animal Genetics, Breeding, Reproduction and Molecular Design, Yangzhou, China. Yangzhou University Joint International Research Laboratory of Agriculture & Agri-Product Safety Yangzhou China Yangzhou University, Joint International Research Laboratory of Agriculture & Agri-Product Safety, Yangzhou, China. Corresponding author: pigbreeding@163.com

Conflict of Interest

The authors declare no conflict of interest.

 10 04 2023 2023 52 e20210149 20 08 2021 12 09 2022 This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. ABSTRACT

To characterize the N6-methyladenosine (m6A)-related gene expression profiles in various tissues of Meishan pigs at different stages, m6A modification-related genes (METTL3, METTL14, METTL16, WTAP, RBM15, and FTO) were detected from newborn to physical maturity of Meishan pigs at eight important developmental stages (1, 7, 14, 21, 28, 35, 134, and 158 days old). The expression of m6A-related genes was tissue-specific. Furthermore, the level of METTL3 messenger RNA (mRNA) was higher on day 35 than in other stages in most tissues, and the expression of METTL14 increased after day 35, and FTO exhibited a peak on day 14 in muscle, intestine, lymph nodes, thymus, and kidney. This study provided a reference for an in-depth study of the expression patterns of m6A modification-related genes in Meishan pigs.

 FTO METTL3 tissue expression profile weaned piglets National Natural Science Foundation of China 31772560 Key Research and Development Project BE2019344 Key Research and Development Project BE2019341 1. Introduction

The Chinese Meishan pig breed is well known as one of the most prolific breeds in the world. Meishan pig has an excellent high reproductive performance, including more litters, strong lactation ability, higher conception rate, and early maturity (Ellis et al., 1995; Désautés et al., 1999; Groenen et al., 2012; Zheng et al., 2020). Besides the great reproductive performance, the Meishan pig breed also has the advantage of tolerating rough feeds and environments (Christenson et al., 1993; Haley and Lee, 1993; Zheng et al., 2020). Hence, it is necessary to discover the reason why the Meishan pig breed possesses the superior phenotypes.

Several studies have focused on the identification of genetic diversity and population structure to unravel functional genes and found that the harbored gene insulin like growth factor 1 receptor (IGF1R) may contribute to the higher fertility of Meishan pigs. Nuclear factor-kappaB (NF-κB) as the key gene of NF-κB signaling, positively regulates the hyaluronan biosynthesis, which may explain the wrinkled skin and face of Meishan pigs (Zhao et al., 2018; Zhou et al., 2021). Previous research studying the protein expression patterns in endometrial tissue from Meishan and Duroc sows revealed that there were 114 differentially expressed proteins identified at day 49 and 98 differentially expressed proteins were identified at day 72 during pregnancy, suggesting that the differentially expressed proteins may be a major factor influencing the differences in embryonic loss between Meishan and Duroc sows during mid-gestation (Wang et al., 2019a). However, given the high complexity of reproduction performance and related phenotypes, the genetic basis of superior characteristics in Meishan pigs remains largely unknown.

In eukaryotes, N6-methyladenosine (m6A) is the most abundant form of methylation modification and is critical for various physiological and pathological processes. The m6A methylation is a dynamic process, which is added by methyltransferases (writers) or removed by demethylases (erasers) and exerts its function by m6A binding proteins (readers) (Wang et al., 2020a). Methyltransferases are a complex that contains methyltransferase-like protein 3 (METTL3) and METTL14. Other methyltransferases also have been identified, such as wilms tumor 1 associated protein (WTAP; Ping et al., 2014), Virilizer (KIAA1429; Jiang et al., 2021a), and RNA binding motif protein 15 (RBM15; Zhang et al., 2015). Conversely, two demethylases, fat mass and obesity-associated protein (FTO; Jia et al., 2011) and AlkB homolog 5 (ALKBH5; Zheng et al., 2013) have the function to remove the m6A modification. Additionally, the specific readers mediate the function of m6A modification on target mRNA. Members of the YTH family have been identified to specifically bind to m6A-containing precursor RNA. The members of YTH family contain YTHDF1, YTHDF2, and YTHDF3, which are located mainly in the cytoplasm (Zaccara and Jaffrey, 2020).

Research has shown that m6A modification plays critical roles in numerous biological functions such as the differentiation of stem cells (Zhang et al., 2016; Chen et al., 2019), maturation of egg cells (Kasowitz et al., 2018), and development of brain and metabolism of adipose (Wang et al., 2020b). METTL3 plays multiple roles in biological process, tumorigenesis, and cell proliferation (Li et al., 2017b; Liu et al., 2019). Additionally, METTL3 also has the function to modulate virus replication (Hao et al., 2019). METTL14 was identified as the other writer complex component and provides a platform with METTL3 for RNA recognition (Wang et al., 2017). Knocking down METTL14 leads to decreased m6A levels in human cell lines (Schwartz et al., 2014). Besides, knockdown of METTL14 inhibits the differentiation and promotes the proliferation of C2C12 myoblast cells (Zhang et al., 2020). FTO, a demethylase, except for being associated with obesity (Loos and Yeo, 2014), is a crucial component of m6A modification. FTO plays a critical role in occurrence, progression and treatment of various cancers, and even acts as a cancer oncogene in acute myeloid leukemia (Chen and Du, 2019). Previous study found that FTO level is increased in human melanoma and enhances melanoma tumorigenesis in mice, and that FTO regulates melanoma tumorigenicity and response to anti-PD-1 blockade (Yang et al., 2019). Another study also reported that m6A modification regulates the expression of the growth arrest and DNA damage-inducible 45B (GADD45B) gene by activation of P38/MAPK signaling pathway, in which is involved in myogenic differentiation (Deng et al., 2021). Recent study compared the differences in m6A methylation pattern between fat and lean broilers and found that the high m6A methylated genes (fat birds vs. lean birds) are primarily involved in fatty acid biosynthesis and fatty acid metabolism, while the low m6A methylated genes are mainly participated in processes associated with development (Cheng et al., 2021). Additionally, the m6A methylation regulates the embryo viability and germline development in mouse. The nuclear m6A reader YTHDC1is required for spermatogonial development in males and for oocyte growth and maturation in females (Kasowitz et al., 2018). However, there are few studies about the expression pattern of m6A modification-related genes in various tissues of Meishan pigs at different developmental stages.

This study focused on the specific expression of m6A modification-related genes in different tissues at the different stages of Meishan pigs and aimed to explore the relationship between m6A modification and the development of Meishan pigs at different stages.

 2. Material and Methods 2.1. Animals and collection of samples

All animal procedures were approved by the Institutional Animal Care and Use Committee (IACUC) of the Animal Experiments Ethics Committee, China (SYXK (Su) 2016-0019). All experimental procedures were carried out in accordance with the Regulations for the Administration of Affairs Concerning Experimental Animals. The Meishan pigs were obtained from the Meishan Pig National Breeding Farm (Taicang City, Jiangsu Province). Pigs were chosen based on the similarity of fur color, body shape, and weight from five litters at the various developmental stages (1, 7, 14, 21, 28, 35, 134, and 158 days old). Forty pigs (five pigs per group) were used; the weight information of pigs was listed in Table 1. All pigs were raised in the same condition, fed commercial fodder, and had free access to water.

Weight information of pigs
Number of pigs in group Sex Age (days) Weight (kg)
5 2 females and 3 males 1 1.04±0.21
5 3 females and 2 males 7 1.77±0.43
5 1 female and 4 males 14 3.23±0.41
5 3 females and 2 males 21 4.36±0.73
5 2 females and 3 males 28 5.80±0.49
5 3 females and 2 males 35 7.94±0.51
5 4 females and 1male 134 35.44+3.50
5 3 females and 2 males 158 68.10±4.05

Tissue samples (muscle, heart, duodenum, jejunum, ileum, stomach, liver, spleen, thymus, lymph node, kidney, and lung tissues) at the same position from different pigs were rapidly excised, the residual blood was washed with normal saline, and the tissue was dried with filter paper to remove surface liquid. The samples were then placed in 1.5 mL nuclease-free Eppendorf tubes, and immediately placed into liquid nitrogen, and stored at −80 ℃ until RNA extraction.

 2.2. Primers design

The primers (Table 2) were designed using software Primer Premier 5.0 referenced to the sequence from Genebank data.

Primer sequence of genes for qPCR
Gene Sequence of the primer Length of products (bp)
METTL3 F: 5’- CCCTATGGGACCCTGACAGA- 3’ R: 5’- TGACACCAACCAAGCAGTGT- 3’ 250
METTL14 F: 5’- GGGAGAGTGTGTTTACGCAAG- 3’ R: 5’- TGAAGTCCCCGTCTGTGCTA- 3’ 184
METTL16 F: 5’- ATTAAGAAAGAGGTTCTGGAGAGGC- 3’ R: 5’- AGGGCATGCCTTCAATCATCA- 3’ 126
WTAP F: 5’- TCCATTCGTCTTTCCTCTCCG- 3’ R: 5’- GCCTCACTCAGTCGAACCTTT- 3’ 130
RBM15 F: 5’- CATTGTCCGTGCAAAACTGGT- 3’ R: 5’- ACTTAAAACACCGGCATTGGC- 3’ 211
FTO F: 5’- GATCTCAATGCCACCCACCA- 3’ R: 5’- CCACTCAAACTCGACCTCGT- 3’ 237
GAPDH F: 5’- ACATCATCCCTGCTTCTACTGG- 3’ R: 5’- CTCGGACGCCTGCTTCAC- 3’ 187

mRNA level of METTL3 in heart and muscle (A); ileum, jejunum, duodenum, stomach, and liver (B); lymph nodes, thymus, and spleen (C); kidney and lung (D).

*P<0.05, **P<0.01.

2.3. qRT-PCR

Total RNA was extracted by RNAiso Plus (Takara, Japan) according to the manufacturer’s instructions. Briefly, tissues were homogenized by adding 1 mL of RNAiso Plus, then 200 µL of chloroform was mixed until the solution became milky, and was kept at room temperature for 5 min, then centrifuged at 12000 g for 15 mins at 4 ℃. The upper layer was transferred to a new centrifuge tube and equal volume of isopropanol was added, and the mixture was kept at room temperature for 10 min and centrifuged at 12000 g for 20 mins at 4 ℃. The supernatant was carefully removed and washed with 75% ethanol twice, the supernatant was then discarded. Finally, the precipitated RNA was dried and dissolved with appropriate amount DPEC-treated water. Reverse transcription of the RNA was performed by using reverse transcriptase (R312-01, Vazyme, China). Two microliter cDNA templates were added to the reaction system (Q111-02, Vazyme, China). The qRT-PCR were programmed as denaturation at 95 °C for 5 min, followed by 40 cycles of denaturation at 95 °C for 10 s, annealing/elongation at 65 °C for 30 s. From three independent experiments, different gene expressions were calculated by using the comparative Ct method and normalized against GAPDH level.

 2.4. Statistical analysis

The results are expressed as means ± standard error of mean (SEM) from at least three independent experiments. Statistical analysis was performed by SPSS version 16.0. The difference among different groups was analyzed by Student’s t-test and One-way Analysis of Variance (ANOVA), with *P<0.05 and **P<0.01.

 3. Results 3.1. <italic>METTL3</italic> mRNA levels in different tissues of Meishan pigs at different developmental stages

METTL3 is the core element of m6A modification. The levels of METTL3 mRNA in different tissues of Meishan pigs at different developmental stages were detected by qRT-PCR. The organs were classified according to their functions. METTL3 displayed a similar expression pattern in different tissues, and the results showed that the expression of METTL3 peaked at 35 days of age in almost all organs except the heart (Figure 1). In the digestive system, the level of METTL3 in the liver was higher on day 35 and increased after day 134 in ileum, jejunum, stomach, and liver (Figure 1B). As the largest immune organ, it was worth mentioning that the METTL3 gene was highly expressed in the spleen at the weaning stage (Figure 1C). Regardless of the developmental stage, the level of METTL3 mRNA is much lower in heart and muscle tissue compared with other tissues.

Expression of <italic>METTL3</italic> in different tissues at different stages by using the comparative Ct method.

mRNA level of METTL14 in heart and muscle (A); ileum, jejunum, duodenum, stomach, and liver (B); lymph nodes, thymus, and spleen (C); kidney and lung (D).

*P<0.05, **P<0.01.

 3.2. <italic>METTL14</italic> mRNA levels in different tissues of Meishan pigs at different developmental stages

METTL14 can form a stable heterodimer with METTL3, and both co-localize in the nucleus and can maintain the stability of each other’s protein. The mRNA level of METTL14 was also detected and displayed a different expression pattern from METTL3, and the expression of METTL14 in the heart at different stages remained low (Figure 2A). The changing trends of each part were consistent in intestinal tissues, and the expression of METTL14 in liver increased after day 35 (Figure 2B). Besides, a peak expression appeared in lymph nodes at seven days of age, and the mRNA level increased in lymphatic organs from 35 days old (Figure 2C). Additionally, the expression of METTL14 in the kidney and lung was up-regulated after day 134 (Figure 2D).

Expression of <italic>METTL14</italic> in different tissues at different stages by using the comparative Ct method.

mRNA level of METTL16 in heart and muscle (A); ileum, jejunum, duodenum, stomach, and liver (B); lymph nodes, thymus, and spleen (C); kidney and lung (D).

*P<0.05, **P<0.01.

 3.3. <italic>METTL16</italic> mRNA levels in different tissues of Meishan pigs at different development stages

The expression of METTL16 was higher in muscle at 35 days of age (Figure 3A). In addition, the expression of METTL16 was lower before 134 days old and became high at 158 days old in sections of small intestine and stomach (Figure 3B). The expression of METTL16 was higher on day 28 in the thymus than in other developmental stages (Figure 3C). However, the expression of METTL16 in the kidney presented a trend of first increasing and then decreasing, while in the lung, it fluctuated (Figure 3D).

Expression of <italic>METTL16</italic> in different tissues at different stages by using the comparative Ct method.

mRNA level of WTAP in heart and muscle (A); ileum, jejunum, duodenum, stomach, and liver (B); lymph nodes, thymus, and spleen (C); kidney and lung (D).

*P<0.05, **P<0.01.

 3.4. <italic>WTAP</italic> mRNA levels in different tissues of Meishan pigs at different development stages

WTAP is a conserved nuclear protein as the partner of Wilms’ tumor 1, and removal of WTAP has been approved to be embryonically lethal. WTAP is a member of methylases and involves in multiple biological functions. We found that the expression of WTAP escalated in most organs after day 134 (Figure 4A-4D). In addition, the expression in the immune organs, such as spleen, lymph nodes, and thymus, exhibited a small peak on day 35, indicating that the change in methylation level of the immune system may be caused by weaning of piglets (Figure 4C). Moreover, the expression of WTAP showed a peak in kidney on day 14 (Figure 4D).

Expression of <italic>WTAP</italic> in different tissues at different stages by using the comparative Ct method.

mRNA level of RBM15 in heart and muscle (A); ileum, jejunum, duodenum, stomach, and liver (B); lymph nodes, thymus, and spleen (C); kidney and lung (D).

*P<0.05, **P<0.01.

 3.5. <italic>RBM15</italic> mRNA levels in different tissues of Meishan pigs at different developmental stages

RBM15 is an RNA-binding protein that is part of the WTAP-METTL3 m6A methyltransferase complex. RBM15 and its paralogue RBM15B interact with WTAP to recruit the complex to target mRNA and are essential for the long non-coding RNA X-inactive specific transcript (XIST)-mediated gene silencing (Patil et al., 2016). In this study, the results showed that the expression of RBM15 is higher on day 14 in the heart, and then gradually decreased, whereas the RBM15 in muscle maintained a lower expression before day 35 and increased after day 134 (Figure 5A). The expression of RBM15 was relatively higher on day 14 in liver, duodenum, and stomach (Figure 5B). Besides, RBM15 mRNA level showed a short rise on day 28 in the thymus and spleen (Figure 5C). The RBM15 in the kidney and lung displayed a similar expression pattern at different developmental stages (Figure 5D). These results indicated that RBM15 showed a tissue-specific expression.

Expression of <italic>RBM15</italic> in different tissues at different stages by using the comparative Ct method.

mRNA level of FTO in heart and muscle (A); ileum, jejunum, duodenum, stomach, and liver (B); lymph nodes, thymus, and spleen (C); kidney and lung (D).

*P<0.05, **P<0.01.

 3.6. <italic>FTO</italic> mRNA levels in different tissues of Meishan pigs at different developmental stages

The FTO mRNA level was low after day 21 in muscle and heart (Figure 6A) and showed a higher expression than other stages on day 14 in various parts of the intestine (Figure 6B), lymph nodes, and thymus (Figure 6C). Additionally, we also found that the expression was up-regulated after weaning on day 35 in the lymph nodes, spleen, and kidney (Figure 6C and 6D). Intriguingly, the FTO mRNA level in the kidney and lung appeared in a converse situation (Figure 6D).

Expression of <italic>FTO</italic> in different tissues at different stages by using the comparative Ct method.

4. Discussion

Meishan pig breed is one of the most famous Chinese domestic pig breeds, which is renowned for the high prolificacy, high reproductive performance, great meat physical properties, and tolerance for rough feed (Zhao et al., 2018). To reveal the mechanisms of Meishan pigs with excellent economic traits, researchers have used the technique of molecular genetics to explore why Meishan pig breed has good productive performance. The m6A modification has been found to be involved in regulating a variety of physiological activities (Cao et al., 2016), but compared with the studies about humans and mice, there are few studies underscoring the influences of m6A modification on the productive performance of Meishan pigs.

In this study, we performed the expression changes of m6A modification-related genes by qRT-PCR in various tissues of Meishan pigs at different developmental stages. The results showed obvious tissue specificity and preference for the expression of m6A-related genes of Meishan pigs. Importantly, we found that the expression of METTL3 showed a similar trend in different organs at different developmental stages, which was higher on day 35. Moreover, METTL14 and WTAP mRNA levels were increased after day 134. As one of the demethylases, the expression of FTO was up-regulated in the heart, intestine, lymph nodes, thymus, kidney and lung on day 14.

The m6A methylation plays important roles in the development of organisms, which can modulate stem cell specification and cell development and cell fate (Zhang et al., 2017). m6A is installed by m6A methyltransferases, removed by m6A demethylases and recognized by reader proteins. m6A methylation regulates RNA metabolism including translation, splicing, export, degradation, and microRNA processing (He et al., 2019). Additionally, the m6A methylation regulates the growth and development of ovaries and testes. Jiang et al. (2021b) found that FTO-knockdown granulosa cells show faster aging-related phenotypes and FTO retards FOS-dependent ovarian aging. Study also found that the METTL3-mediated mRNA N6-methylademosine plays important roles in oocyte and follicle development of mice (Mu et al., 2021). Moreover, the m6A modification is essential for oogenesis, and researchers proved that there was significant enrichment of m6A methylation-related genes in several signaling pathways associated with the development of ovary. Besides, study also has confirmed that METTL3, METTL14, YTHDF1, YTHDF2, YTHDF3, and KIAA1429 were involved in ovulation (Chen et al., 2022). METTL3/METTL14 mediated mRNA N6-methyladenosine participates in murine spermatogenesis (Lin et al., 2017). Sun et al. (2020) also investigated the dynamic status of m6A during the development of ovary and testis, and found that the m6A level increases with age in both females and males. These studies all indicated that the m6A modification plays a vital role in the growth and development of animals.

In this study, we focused on the dynamic m6A-related genes expression patterns of Meishan pigs at different developmental stages. The expression levels of m6A-related genes were compared based on functions and cell composition of tissues. Interestingly, the different tissues exhibited distinct expression patterns of m6A-related genes, and we found that the expression of m6A-related genes in muscle and heart was low than in other tissues, maybe because these two organs are made up with muscle cells.

Moreover, previous studies indicated that the METTL3 has the capability to promote the activation of dendritic cells and regulates immune response of the host (Wang et al., 2019b). Herein, we found that the expression patterns of METTL3 in different tissues were similar, and the expression of METTL3 and METTL14 was higher on day 35, which indicated that environmental changes, including weaning stress, may lead to the changes. Besides, the expression of WTAP also showed similar trend in different tissues; the expression was lower on day 134 and then upregulated on day 158, and it was reported that WTAP may participate in fat deposition after day 134 (Heng et al., 2019). Combined with these results, we suggested that WTAP play an important role in finishing pigs.

FTO is a member of Fe (II)-and oxoglutarate-dependent AlkB dioxygenase family and is closely related to obesity and intellectual disability (Li et al., 2017a). Previous study found that loss of FTO leads to the brain size and body weight change (Li et al., 2017a). Moreover, FTO could regulate immune response. Studies showed that tumors use FTO as an epitranscriptomic regulator to escape immune surveillance (Su et al., 2020), and FTO promotes M1 and M2 macrophage activation (Gu et al., 2020). However, in this study, as the demethylase, the expression of FTO showed different trends in spleen, lymph nodes, and thymus, which was highly expressed in the lymph nodes at early stage and then decreased from day 21. Compared with that in lymph node, the FTO showed diverse expression patterns in spleen and thymus, which was lowly expressed in early stage and highly expressed in late stage; the composition of immune cells in different immune tissues may contribute to this phenomenon. In digestive system, FTO mRNA level was highly expressed on day 14 and then remained at relatively low level, which may be related to the role of FTO in regulating fat and weight. Additionally, although the expression of FTO is higher in patients with lung cancer and lower in patients with nephritis (Liu et al., 2018; Zhao et al., 2021), the diverse expression of FTO in lung and kidney may due to the fact that FTO play different roles the in kidney and lung. Apart from Meishan pigs, the m6A modification takes part in the growth and development of other pig breeds. He et al. (2017) profiled transcriptome-wide m6A in porcine liver at three developmental stages of Rongchang pigs: newborn (0 day), suckling (21 days), and adult (two years) and found that one third transcribed genes are modified by m6A. Wang et al. (2022) reported the m6A-methylation patterns of lncRNA via MeRIP-seq and found that m6A methylation regulates the muscle-fiber-type conversion of Duroc pigs. Taken together, the dynamic regulation of m6A modification may contribute to the excellent production performance of Meishan pigs.

In this study, three important points were associated with the changes of m6A-related genes, including 7 to 14 days old, 35 days old, and 134 days old. Generally, the newborn piglets are susceptible to kinds of pathogens, and the piglets are usually weaned on day 35. In those three different developmental stages, the expression of m6A-related genes was usually altered. The immune protection of piglets mainly derives from two aspects. The immunity obtained from breast milk is passive immunity (Myles and Datta, 2021), and the immune protection formed by the development of the autoimmune system is active immunity (Baxter, 2014).

In this study, we noticed a phenomenon that the methyltransferase METTL3 was lowly expressed in the thymus, spleen, and lymph nodes before 35 days of age, but highly expressed at 35 days of age. We suspected that weaning may be responsible for this phenomenon. The first two weeks after weaning is the key turning point of the early developmental stage of piglets. Moreover, weaning usually causes the interruption of maternal antibodies, while self-digestion and immune functions are not yet complete in piglets. Therefore, in this period, the metabolism of main immune organs like the thymus and spleen is very rapid. These organs are in the peak period of activity and are greatly affected by maternal factors. During the development of thymus and spleen, the high expression of methyltransferases is involved in the development of immunity to against weaning stress.

Furthermore, newborn piglets usually obtain immunoprotection from breast milk, and the immune efficacy reaches peak at seven days old. Besides, the autoimmune system of piglets works from 28 to 35 days old, which may explain the higher expression of METTL3 on day 35. Apart from affecting the immune function, weaning stress affects intake, production performance, and development of piglets. In the intestinal tissues, the expression of METTL3 also was upregulated on day 35. Hence, the dynamic regulation of m6A modification during weaning stress stage may contribute to the function of intestine and immunity of piglets.

 5. Conclusions

The expression of METTL3 is higher in different tissues on day 35, which is usually the weaning date. The expression of METTL14 increases after day 35. Compared with the expression of METTL3 on day 35, the expression of FTO exhibits a lower level, indicating that organism regulates m6A methylation in response to weaning stress. Collectively, our study provided a reference for the expression pattern of m6A methylation-related genes in different tissues in Meishan pigs.

 Acknowledgments

This research was funded by the National Natural Science Foundation of China (31772560), Key Research and Development Project (Modern Agriculture) of Jiangsu Province (BE2019344, BE2019341), and the Priority Academic Program Development of Jiangsu Higher Education Institutions.

 References Baxter, D. 2014. Active and passive immunization for cancer. Human Vaccines & Immunotherapeutics 10:2123-2129. https://doi.org/10.4161/hv.29604 Baxter D 2014 Active and passive immunization for cancer Human Vaccines & Immunotherapeutics 10 2123 2129 https://doi.org/10.4161/hv.29604 Cao, G.; Li, H.-B.; Yin, Z. and Flavell, R. A. 2016. Recent advances in dynamic m6A RNA modification. Open Biology 6:160003. https://doi.org/10.1098/rsob.160003 Cao G Li H.-B Yin Z Flavell R. A. 2016 Recent advances in dynamic m6A RNA modification Open Biology 6 160003 https://doi.org/10.1098/rsob.160003 Chen, J. and Du, B. 2019. Novel positioning from obesity to cancer: FTO, an m6A RNA demethylase, regulates tumour progression. Journal of Cancer Research and Clinical Oncology 145:19-29. https://doi.org/10.1007/s00432-018-2796-0 Chen J Du B 2019 Novel positioning from obesity to cancer: FTO, an m6A RNA demethylase, regulates tumour progression Journal of Cancer Research and Clinical Oncology 145 19 29 https://doi.org/10.1007/s00432-018-2796-0 Chen, J.; Fang, Y.; Xu, Y. and Sun, H. 2022. Role of m6A modification in female infertility and reproductive system diseases. International Journal of Biological Sciences 18:3592-3604. https://doi.org/10.7150/ijbs.69771 Chen J. Fang Y. Xu Y. Sun H. 2022 Role of m6A modification in female infertility and reproductive system diseases International Journal of Biological Sciences 18 3592 3604 https://doi.org/10.7150/ijbs.69771 Chen, J.; Wang, C.; Fei, W.; Fang, X. and Hu, X. 2019. Epitranscriptomic m6A modification in the stem cell field and its effects on cell death and survival. American Journal of Cancer Research 9:752-764. Chen J. Wang C. Fei W. Fang X. Hu X. 2019 Epitranscriptomic m6A modification in the stem cell field and its effects on cell death and survival American Journal of Cancer Research 9 752 764 Cheng, B.; Leng, L.; Li, Z.; Wang, W.; Jing, Y.; Li, Y.; Wang, N.; Li, H. and Wang, S. 2021. Profiling of RNA N6-methyladenosine methylation reveals the critical role of m6A in chicken adipose deposition. Frontiers in Cell and Developmental Biology 9:590468. https://doi.org/10.3389/fcell.2021.590468 Cheng B. Leng L. Li Z. Wang W. Jing Y. Li Y. Wang N. Li H. Wang S. 2021 Profiling of RNA N6-methyladenosine methylation reveals the critical role of m6A in chicken adipose deposition Frontiers in Cell and Developmental Biology 9 590468 https://doi.org/10.3389/fcell.2021.590468 Christenson, R. K.; Vallet, J. L.; Leymaster, K. A. and Young, L. D. 1993. Uterine function in Meishan pigs. Journal of Reproduction and Fertility Supplement 48:279-289. Christenson R. K Vallet J. L Leymaster K. A Young L. D 1993 Uterine function in Meishan pigs Journal of Reproduction and Fertility Supplement 48 279 289 Deng, K.; Fan, Y.; Liang, Y.; Cai, Y.; Zhang, G.; Deng, M.; Wang, Z.; Lu, J.; Shi, J.; Wang, F. and Zhang, Y. 2021. FTO-mediated demethylation of GADD45B promotes myogenesis through the activation of p38 MAPK pathway. Molecular Therapy: Nucleic Acids 26:34-48. https://doi.org/10.1016/j.omtn.2021.06.013 Deng K. Fan Y. Liang Y. Cai Y. Zhang G. Deng M. Wang Z. Lu J. Shi J. Wang F. Zhang Y. 2021 FTO-mediated demethylation of GADD45B promotes myogenesis through the activation of p38 MAPK pathway Molecular Therapy: Nucleic Acids 26 34 48 https://doi.org/10.1016/j.omtn.2021.06.013 Désautés, C.; Sarrieau, A.; Caritez, J.-C. and Mormède, P. 1999. Behavior and pituitary-adrenal function in Large White and Meishan pigs. Domestic Animal Endocrinology 16:193-205. Désautés C Sarrieau A Caritez J.-C Mormède P 1999 Behavior and pituitary-adrenal function in Large White and Meishan pigs Domestic Animal Endocrinology 16 193 205 Ellis, M.; Lympany, C.; Haley, C. S.; Brown, I. and Warkup, C. C. 1995. The eating quality of pork from Meishan and Large White pigs and their reciprocal crosses. Animal Science 60:125-131. https://doi.org/10.1017/S1357729800008225 Ellis M. Lympany C. Haley C. S. Brown I. Warkup C. C. 1995 The eating quality of pork from Meishan and Large White pigs and their reciprocal crosses Animal Science 60 125 131 https://doi.org/10.1017/S1357729800008225 Groenen, M. A. M.; Archibald, A. L.; Uenishi, H.; Tuggle, C. K.; Takeuchi, Y.; Rothschild, M. F.; Rogel-Gaillard, C.; Park, C.; Milan, D.; Megens, H.-J.; Li, S. T.; Larkin, D. M.; Kim, H.; Frantz, L. A. F.; Caccamo, M.; Ahn, H.; Aken, B. L.; Anselmo, A.; Anthon, C.; Auvil, L.; Badaoui, B.; Beattie, C. W.; Bendixen, C.; Berman, D.; Blecha, F.; Blomberg, J.; Bolund, L.; Bosse, M.; Botti, S.; Bujie, Z.; Bystrom, M.; Capitanu, B.; Carvalho-Silva, D.; Chardon, P.; Chen, C.; Cheng, R.; Choi, S. H.; Chow, W.; Clark, R. C.; Clee, C.; Crooijmans, R. P. M. A.; Dawson, H. D.; Dehais, P.; De Sapio, F.; Dibbits, B.; Drou, N.; Du, Z. Q.; Eversole, K.; Fadista, J.; Fairley, S.; Faraut, T.; Faulkner, G. J.; Fowler, K. E.; Fredholm, M.; Fritz, E.; Gilbert, J. G. R.; Giuffra, E.; Gorodkin, J.; Griffin, D. K.; Harrow, J. L.; Hayward, A.; Howe, K.; Hu, Z. L.; Humphray, S. J.; Hunt, T.; Hornshoj, H.; Jeon, J. T.; Jern, P.; Jones, M.; Jurka, J.; Kanamori, H.; Kapetanovic, R.; Kim, J.; Kim, J. H.; Kim, K. W.; Kim, T. H.; Larson, G.; Lee, K.; Lee, K. T.; Leggett, R.; Lewin, H. A.; Li, Y. R.; Liu, W. S.; Loveland, J. E.; Lu, Y.; Lunney, J. K.; Ma, J.; Madsen, O.; Mann, K.; Matthews, L.; McLaren, S.; Morozumi, T.; Murtaugh, M. P.; Narayan, J.; Nguyen, D. T.; Ni, P. X.; Oh, S. J.; Onteru, S.; Panitz, F.; Park, E. W.; Park, H. S.; Pascal, G.; Paudel, Y.; Perez-Enciso, M.; Ramirez-Gonzalez, R.; Reecy, J. M.; Rodriguez-Zas, S.; Rohrer, G. A.; Rund, L.; Sang, Y. M.; Schachtschneider, K.; Schraiber, J. G.; Schwartz, J.; Scobie, L.; Scott, C.; Searle, S.; Servin, B.; Southey, B. R.; Sperber, G.; Stadler, P.; Sweedler, J. V.; Tafer, H.; Thomsen, B.; Wali, R.; Wang, J.; Wang, J.; White, S.; Xu, X.; Yerle, M.; Zhang, G. J.; Zhang, J. G.; Zhang, J.; Zhao, S. H.; Rogers, J.; Churcher, C. and Schook, L. B. 2012. Analyses of pig genomes provide insight into porcine demography and evolution. Nature 491:393-398. https://doi.org/10.1038/nature11622 Groenen M. A. M Archibald A. L Uenishi H Tuggle C. K Takeuchi Y Rothschild M. F Rogel-Gaillard C Park C Milan D Megens H.-J Li S. T Larkin D. M Kim H Frantz L. A. F Caccamo M Ahn H Aken B. L Anselmo A Anthon C Auvil L Badaoui B Beattie C. W Bendixen C Berman D Blecha F Blomberg J Bolund L Bosse M Botti S Bujie Z Bystrom M Capitanu B Carvalho-Silva D Chardon P Chen C Cheng R Choi S. H Chow W Clark R. C Clee C. Crooijmans R. P. M. A. Dawson H. D Dehais P. De Sapio F Dibbits B Drou N Du Z. Q Eversole K Fadista J Fairley S Faraut T Faulkner G. J Fowler K. E Fredholm M Fritz E Gilbert J. G. R Giuffra E Gorodkin J Griffin D. K Harrow J. L Hayward A Howe K Hu Z. L Humphray S. J Hunt T Hornshoj H Jeon J. T Jern P Jones M Jurka J Kanamori H Kapetanovic R Kim J Kim J. H Kim K. W Kim T. H Larson G Lee K Lee K. T Leggett R Lewin H. A Li Y. R Liu W. S Loveland J. E Lu Y Lunney J. K Ma J Madsen O Mann K Matthews L McLaren S Morozumi T Murtaugh M. P Narayan J Nguyen D. T Ni P. X Oh S. J Onteru S Panitz F Park E. W Park H. S Pascal G Paudel Y Perez-Enciso M Ramirez-Gonzalez R Reecy J. M Rodriguez-Zas S Rohrer G. A Rund L Sang Y. M Schachtschneider K Schraiber J. G Schwartz J Scobie L Scott C Searle S Servin B Southey B. R Sperber G Stadler P Sweedler J. V Tafer H Thomsen B Wali R Wang J Wang J White S Xu X Yerle M Zhang G. J Zhang J. G Zhang J Zhao S. H Rogers J Churcher C Schook L. B 2012 Analyses of pig genomes provide insight into porcine demography and evolution Nature 491 393 398 https://doi.org/10.1038/nature11622 Gu, X.; Zhang, Y.; Li, D.; Cai, H.; Cai, L. and Xu, Q. 2020. N6-methyladenosine demethylase FTO promotes M1 and M2 macrophage activation. Cellular Signalling 69:109553. https://doi.org/10.1016/j.cellsig.2020.109553 Gu X. Zhang Y. Li D. Cai H. Cai L. Xu Q. 2020 N6-methyladenosine demethylase FTO promotes M1 and M2 macrophage activation Cellular Signalling 69 109553 https://doi.org/10.1016/j.cellsig.2020.109553 Haley, C. S. and Lee, G. J. 1993. Genetic basis of prolificacy in Meishan pigs. Journal of Reproduction and Fertility Supplement 48:247-259. Haley C. S Lee G. J 1993 Genetic basis of prolificacy in Meishan pigs Journal of Reproduction and Fertility Supplement 48 247 259 Hao, H.; Hao, S.; Chen, H.; Chen, Z.; Zhang, Y.; Wang, J.; Wang, H.; Zhang, B.; Qiu, J.; Deng, F. and Guan, W. 2019. N6-methyladenosine modification and METTL3 modulate enterovirus 71 replication. Nucleic Acids Research 47:362-374. https://doi.org/10.1093/nar/gky1007 Hao H. Hao S. Chen H. Chen Z. Zhang Y. Wang J. Wang H. Zhang B. Qiu J. Deng F. Guan W. 2019 N6-methyladenosine modification and METTL3 modulate enterovirus 71 replication Nucleic Acids Research 47 362 374 https://doi.org/10.1093/nar/gky1007 He, L.; Li, H.; Wu, A.; Peng, Y.; Shu, G. and Yin, G. 2019. Functions of N6-methyladenosine and its role in cancer. Molecular Cancer 18:176. https://doi.org/10.1186/s12943-019-1109-9 He L. Li H. Wu A. Peng Y. Shu G. Yin G. 2019 Functions of N6-methyladenosine and its role in cancer Molecular Cancer 18 176 https://doi.org/10.1186/s12943-019-1109-9 He, S.; Wang, H.; Liu, R.; He, M.; Che, T.; Jin, L.; Deng, L.; Tian, S.; Li, Y.; Lu, H.; Li, X.; Jiang, Z.; Li, D. and Li, M. 2017. mRNA N6-methyladenosine methylation of postnatal liver development in pig. PloS One 12:e0173421. https://doi.org/10.1371/journal.pone.0173421 He S. Wang H. Liu R. He M. Che T. Jin L. Deng L. Tian S. Li Y. Lu H. Li X. Jiang Z. Li D. Li M. 2017 mRNA N6-methyladenosine methylation of postnatal liver development in pig PloS One 12 e0173421 https://doi.org/10.1371/journal.pone.0173421 Heng, J.; Tian, M.; Zhang, W.; Chen, F.; Guan, W. and Zhang, S. 2019. Maternal heat stress regulates the early fat deposition partly through modification of m6A RNA methylation in neonatal piglets. Cell Stress and Chaperones 24:635-645. https://doi.org/10.1007/s12192-019-01002-1 Heng J. Tian M. Zhang W. Chen F. Guan W. Zhang S. 2019 Maternal heat stress regulates the early fat deposition partly through modification of m6A RNA methylation in neonatal piglets Cell Stress and Chaperones 24 635 645 https://doi.org/10.1007/s12192-019-01002-1 Jia, G.; Fu, Y.; Zhao, X.; Dai, Q.; Zheng, G.; Yang, Y.; Yi, C.; Lindahl, T.; Pan, T.; Yang, Y.-G. and He, C. 2011. N6-methyladenosine in nuclear RNA is a major substrate of the obesity-associated FTO. Nature Chemical Biology 7:885-887. https://doi.org/10.1038/nchembio.687 Jia G Fu Y Zhao X Dai Q Zheng G Yang Y Yi C Lindahl T Pan T Yang Y.-G He C 2011 N6-methyladenosine in nuclear RNA is a major substrate of the obesity-associated FTO Nature Chemical Biology 7 885 887 https://doi.org/10.1038/nchembio.687 Jiang, X.; Liu, B.; Nie, Z.; Duan, L.; Xiong, Q.; Jin, Z.; Yang, C. and Chen, Y. 2021a. The role of m6A modification in the biological functions and diseases. Signal Transduction and Targeted Therapy 6:74. https://doi.org/10.1038/s41392-020-00450-x Jiang X. Liu B. Nie Z. Duan L. Xiong Q. Jin Z. Yang C. Chen Y. 2021a The role of m6A modification in the biological functions and diseases Signal Transduction and Targeted Therapy 6 74 https://doi.org/10.1038/s41392-020-00450-x Jiang, Z.-X.; Wang, Y.-N.; Li, Z.-Y.; Dai, Z.-H.; He, Y.; Chu, K.; Gu, J.-Y.; Ji, Y.-X.; Sun, N.-X.; Yang, F. and Li, W. 2021b. The m6A mRNA demethylase FTO in granulosa cells retards FOS-dependent ovarian aging. Cell Death & Disease 12:744. https://doi.org/10.1038/s41419-021-04016-9 Jiang Z.-X Wang Y.-N Li Z.-Y Dai Z.-H He Y Chu K Gu J.-Y Ji Y.-X Sun N.-X Yang F Li W 2021b The m6A mRNA demethylase FTO in granulosa cells retards FOS-dependent ovarian aging Cell Death & Disease 12 744 https://doi.org/10.1038/s41419-021-04016-9 Kasowitz, S. D.; Ma, J.; Anderson, S. J.; Leu, N. A.; Xu, Y.; Gregory, B. D.; Schultz, R. M. and Wang, P. J. 2018. Nuclear m6A reader YTHDC1 regulates alternative polyadenylation and splicing during mouse oocyte development. PLoS Genetics 14:e1007412. Kasowitz S. D. Ma J. Anderson S. J. Leu N. A. Xu Y. Gregory B. D. Schultz R. M. Wang P. J. 2018 Nuclear m6A reader YTHDC1 regulates alternative polyadenylation and splicing during mouse oocyte development PLoS Genetics 14 e1007412 Li, L.; Zang, L.; Zhang, F.; Chen, J.; Shen, H.; Shu, L.; Liang, F.; Feng, C.; Chen, D.; Tao, H.; Xu, T.; Li, Z.; Kang, Y.; Wu, H.; Tang, H.; Zhang, P.; Jin, P.; Shu, Q. and Li, X. 2017a. Fat mass and obesity-associated (FTO) protein regulates adult neurogenesis. Human Molecular Genetics 26:2398-2411. https://doi.org/10.1093/hmg/ddx128 Li L. Zang L. Zhang F. Chen J. Shen H. Shu L. Liang F. Feng C. Chen D. Tao H. Xu T. Li Z. Kang Y. Wu H. Tang H. Zhang P. Jin P. Shu Q. Li X. 2017a Fat mass and obesity-associated (FTO) protein regulates adult neurogenesis Human Molecular Genetics 26 2398 2411 https://doi.org/10.1093/hmg/ddx128 Li, X.; Tang, J.; Huang, W.; Wang, F.; Li, P.; Qin, C.; Qin, Z.; Zou, Q.; Wei, J.; Hua, L.; Yang, H. and Wang, Z. 2017b. The m6A methyltransferase METTL3: Acting as a tumor suppressor in renal cell carcinoma. Oncotarget 8:96103-96116. https://doi.org/10.18632/oncotarget.21726 Li X. Tang J. Huang W. Wang F. Li P. Qin C. Qin Z. Zou Q. Wei J. Hua L. Yang H. Wang Z. 2017b The m6A methyltransferase METTL3: Acting as a tumor suppressor in renal cell carcinoma Oncotarget 8 96103 96116 https://doi.org/10.18632/oncotarget.21726 Lin, Z.; Hsu, P. J.; Xing, X.; Fang, J.; Lu, Z.; Zou, Q.; Zhang, K.-J.; Zhang, X.; Zhou, Y.; Zhang, T.; Zhang, Y.; Song, W.; Jia, G.; Yang, X.; He, C. and Tong, M. H. 2017. Mettl3-/Mettl14-mediated mRNA N6-methyladenosine modulates murine spermatogenesis. Cell Research 27:1216-1230. https://doi.org/10.1038/cr.2017.117 Lin Z Hsu P. J Xing X Fang J Lu Z Zou Q Zhang K.-J Zhang X Zhou Y Zhang T Zhang Y Song W Jia G Yang X He C Tong M. H. 2017 Mettl3-/Mettl14-mediated mRNA N6-methyladenosine modulates murine spermatogenesis Cell Research 27 1216 1230 https://doi.org/10.1038/cr.2017.117 Liu, J.; Ren, D.; Du, Z.; Wang, H.; Zhang, H. and Jin, Y. 2018. m6A demethylase FTO facilitates tumor progression in lung squamous cell carcinoma by regulating MZF1 expression. Biochemical and Biophysical Research Communications 502:456-464. https://doi.org/10.1016/j.bbrc.2018.05.175 Liu J. Ren D. Du Z. Wang H. Zhang H. Jin Y. 2018 m6A demethylase FTO facilitates tumor progression in lung squamous cell carcinoma by regulating MZF1 expression Biochemical and Biophysical Research Communications 502 456 464 https://doi.org/10.1016/j.bbrc.2018.05.175 Liu, Q.; Zhao, Y.; Wu, R.; Jiang, Q.; Cai, M.; Bi, Z.; Liu, Y.; Yao, Y.; Feng, J.; Wang, Y. and Wang, X. 2019. ZFP217 regulates adipogenesis by controlling mitotic clonal expansion in a METTL3-m6A dependent manner. RNA Biology 16:1785-1793. https://doi.org/10.1080/15476286.2019.1658508 Liu Q. Zhao Y. Wu R. Jiang Q. Cai M. Bi Z. Liu Y. Yao Y. Feng J. Wang Y. Wang X. 2019 ZFP217 regulates adipogenesis by controlling mitotic clonal expansion in a METTL3-m6A dependent manner RNA Biology 16 1785 1793 https://doi.org/10.1080/15476286.2019.1658508 Loos, R. J. F. and Yeo, G. S. H. 2014. The bigger picture of FTO—the first GWAS-identified obesity gene. Nature Reviews Endocrinology 10:51-61. https://doi.org/10.1038/nrendo.2013.227 Loos R. J. F Yeo G. S. H. 2014 The bigger picture of FTO—the first GWAS-identified obesity gene Nature Reviews Endocrinology 10 51 61 https://doi.org/10.1038/nrendo.2013.227 Mu, H.; Zhang, T.; Yang, Y.; Zhang, D.; Gao, J.; Li, J.; Yue, L.; Gao, D.; Shi, B.; Han, Y.; Zhong, L.; Chen, X.; Wang, Z. B.; Lin, Z.; Tong, M. H.; Sun, Q. Y.; Yang, Y. G. and Han, J. 2021. METTL3-mediated mRNA N6-methyladenosine is required for oocyte and follicle development in mice. Cell Death & Disease 12:989. https://doi.org/10.1038/s41419-021-04272-9 Mu H. Zhang T. Yang Y. Zhang D. Gao J. Li J. Yue L. Gao D. Shi B. Han Y. Zhong L. Chen X. Wang Z. B. Lin Z. Tong M. H. Sun Q. Y. Yang Y. G. Han J. 2021 METTL3-mediated mRNA N6-methyladenosine is required for oocyte and follicle development in mice Cell Death & Disease 12 989 https://doi.org/10.1038/s41419-021-04272-9 Myles, I. A. and Datta, S. K. 2021. Frontline science: Breast milk confers passive cellular immunity via CD8-dependent mechanisms. Journal of Leukocyte Biology 109:709-715. https://doi.org/10.1002/JLB.3HI0820-406RR Myles I. A Datta S. K 2021 Frontline science: Breast milk confers passive cellular immunity via CD8-dependent mechanisms Journal of Leukocyte Biology 109 709 715 https://doi.org/10.1002/JLB.3HI0820-406RR Patil, D. P.; Chen, C.-K.; Pickering, B. F.; Chow, A.; Jackson, C.; Guttman, M. and Jaffrey, S. R. 2016. m6A RNA methylation promotes XIST-mediated transcriptional repression. Nature 537:369-373. https://doi.org/10.1038/nature19342 Patil D. P Chen C.-K Pickering B. F Chow A Jackson C Guttman M Jaffrey S. R. 2016 m6A RNA methylation promotes XIST-mediated transcriptional repression Nature 537 369 373 https://doi.org/10.1038/nature19342 Ping, X.-L.; Sun, B.-F.; Wang, L.; Xiao, W.; Yang, X.; Wang, W.-J.; Adhikari, S.; Shi, Y.; Lv, Y.; Chen, Y.-S.; Zhao, X.; Li, A.; Yang, Y.; Dahal, U.; Lou, X. M.; Liu, X.; Huang, J.; Yuan, W. P.; Zhu, X. F.; Cheng, T.; Zhao, Y. L.; Wang, X. Q.; Danielsen, J. M. R.; Liu, F. and Yang, Y. G. 2014. Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase. Cell Research 24:177-189. https://doi.org/10.1038/cr.2014.3 Ping X.-L Sun B.-F Wang L Xiao W Yang X Wang W.-J Adhikari S Shi Y Lv Y Chen Y.-S Zhao X Li A Yang Y Dahal U Lou X. M Liu X Huang J Yuan W. P Zhu X. F Cheng T Zhao Y. L Wang X. Q Danielsen J. M. R Liu F Yang Y. G. 2014 Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase Cell Research 24 177 189 https://doi.org/10.1038/cr.2014.3 Schwartz, S.; Mumbach, M. R.; Jovanovic, M.; Wang, T.; Maciag, K.; Bushkin, G. G.; Mertins, P.; Ter-Ovanesyan, D.; Habib, N.; Cacchiarelli, D.; Sanjana, N. E.; Freinkman, E.; Pacold, M. E.; Satija, R.; Mikkelsen, T. S.; Hacohen, N.; Zhang, F.; Carr, S. A.; Lander, E. S. and Regev, A. 2014. Perturbation of m6A writers reveals two distinct classes of mRNA methylation at internal and 5′ sites. Cell Reports 8:284-296. https://doi.org/10.1016/j.celrep.2014.05.048 Schwartz S. Mumbach M. R. Jovanovic M. Wang T. Maciag K. Bushkin G. G. Mertins P. Ter-Ovanesyan D. Habib N. Cacchiarelli D. Sanjana N. E. Freinkman E. Pacold M. E. Satija R. Mikkelsen T. S. Hacohen N. Zhang F. Carr S. A. Lander E. S. Regev A. 2014 Perturbation of m6A writers reveals two distinct classes of mRNA methylation at internal and 5′ sites Cell Reports 8 284 296 https://doi.org/10.1016/j.celrep.2014.05.048 Su, R.; Dong, L.; Li, Y.; Gao, M.; Han, L.; Wunderlich, M.; Deng, X.; Li, H.; Huang, Y.; Gao, L.; Li, C.; Zhao, Z.; Robinson, S.; Tan, B.; Qing, Y.; Qin, X.; Prince, E.; Xie, J.; Qin, H.; Li, W.; Shen, C.; Sun, J.; Kulkarni, P.; Weng, H.; Huang, H.; Chen, Z.; Zhang, B.; Wu, X.; Olsen, M. J.; Müschen, M.; Marcucci, G.; Salgia, R.; Li, L.; Fathi, A. T.; Li, Z.; Mulloy, J. C.; Wei, M.; Horne, D. and Chen, J. 2020. Targeting FTO suppresses cancer stem cell maintenance and immune evasion. Cancer Cell 38:79-96. https://doi.org/10.1016/j.ccell.2020.04.017 Su R. Dong L. Li Y. Gao M. Han L. Wunderlich M. Deng X. Li H. Huang Y. Gao L. Li C. Zhao Z. Robinson S. Tan B. Qing Y. Qin X. Prince E. Xie J. Qin H. Li W. Shen C. Sun J. Kulkarni P. Weng H. Huang H. Chen Z. Zhang B. Wu X. Olsen M. J. Müschen M. Marcucci G. Salgia R. Li L. Fathi A. T. Li Z. Mulloy J. C. Wei M. Horne D. Chen J. 2020 Targeting FTO suppresses cancer stem cell maintenance and immune evasion Cancer Cell 38 79 96 https://doi.org/10.1016/j.ccell.2020.04.017 Sun, X.; Zhang, J.; Jia, Y.; Shen, W. and Cao, H. 2020. Characterization of m6A in mouse ovary and testis. Clinical and Translational Medicine 10:e141. https://doi.org/10.1002/ctm2.141 Sun X. Zhang J. Jia Y. Shen W. Cao H. 2020 Characterization of m6A in mouse ovary and testis Clinical and Translational Medicine 10 e141 https://doi.org/10.1002/ctm2.141 Wang, K.; Yang, K.; Xu, Q.; Liu, Y.; Li, W.; Bai, Y.; Wang, J.; Ding, C.; Liu, X.; Tang, Q.; Luo, Y.; Zheng, J.; Wu, K. and Fang, M. 2019a. Protein expression profiles in Meishan and Duroc sows during mid-gestation reveal differences affecting uterine capacity, endometrial receptivity, and the maternal–fetal interface. BMC Genomics 20:991. https://doi.org/10.1186/s12864-019-6353-2 Wang K. Yang K. Xu Q. Liu Y. Li W. Bai Y. Wang J. Ding C. Liu X. Tang Q. Luo Y. Zheng J. Wu K. Fang M. 2019a Protein expression profiles in Meishan and Duroc sows during mid-gestation reveal differences affecting uterine capacity, endometrial receptivity, and the maternal–fetal interface BMC Genomics 20 991 https://doi.org/10.1186/s12864-019-6353-2 Wang, H.; Hu, X.; Huang, M.; Liu, J.; Gu, Y.; Ma, L.; Zhou, Q. and Cao, X. 2019b. Mettl3-mediated mRNA m6a methylation promotes dendritic cell activation. Nature Communications 10:1898. https://doi.org/10.1038/s41467-019-09903-6 Wang H. Hu X. Huang M. Liu J. Gu Y. Ma L. Zhou Q. Cao X. 2019b Mettl3-mediated mRNA m6a methylation promotes dendritic cell activation Nature Communications 10 1898 https://doi.org/10.1038/s41467-019-09903-6 Wang, S.; Tan, B.; Xiao, L.; Zhao, X.; Zeng, J.; Hong, L.; Yang, J.; Cai, G.; Zheng, E.; Wu, Z. and Gu, T. 2022. Comprehensive analysis of long noncoding RNA modified by m6A methylation in oxidative and glycolytic skeletal muscles. International Journal of Molecular Sciences 23:4600. https://doi.org/10.3390/ijms23094600 Wang S. Tan B. Xiao L. Zhao X. Zeng J. Hong L. Yang J. Cai G. Zheng E. Wu Z. Gu T. 2022 Comprehensive analysis of long noncoding RNA modified by m6A methylation in oxidative and glycolytic skeletal muscles International Journal of Molecular Sciences 23 4600 https://doi.org/10.3390/ijms23094600 Wang, T.; Kong, S.; Tao, M. and Ju, S. 2020a. The potential role of rna N6-methyladenosine in cancer progression. Molecular Cancer 19:88. https://doi.org/10.1186/s12943-020-01204-7 Wang T. Kong S. Tao M. Ju S. 2020a The potential role of rna N6-methyladenosine in cancer progression Molecular Cancer 19 88 https://doi.org/10.1186/s12943-020-01204-7 Wang, X.; Huang, J.; Zou, T. and Yin, P. 2017. Human m6A writers: Two subunits, 2 roles. RNA Biology 14:300-304. https://doi.org/10.1080/15476286.2017.1282025 Wang X. Huang J. Zou T. Yin P 2017 Human m6A writers: Two subunits, 2 roles RNA Biology 14 300 304 https://doi.org/10.1080/15476286.2017.1282025 Wang, Y.; Gao, M.; Zhu, F.; Li, X.; Yang, Y.; Yan, Q.; Jia, L.; Xie, L. and Chen, Z. 2020b. METTL3 is essential for postnatal development of brown adipose tissue and energy expenditure in mice. Nature Communications 11:1648. https://doi.org/10.1038/s41467-020-15488-2 Wang Y. Gao M. Zhu F. Li X. Yang Y. Yan Q. Jia L. Xie L. Chen Z. 2020b METTL3 is essential for postnatal development of brown adipose tissue and energy expenditure in mice Nature Communications 11 1648 https://doi.org/10.1038/s41467-020-15488-2 Yang, S.; Wei, J.; Cui, Y.-H.; Park, G.; Shah, P.; Deng, Y.; Aplin, A. E.; Lu, Z.; Hwang, S.; He, C. and He, Y.-Y. 2019. m6A mRNA demethylase FTO regulates melanoma tumorigenicity and response to anti-PD-1 blockade. Nature Communications 10:2782. https://doi.org/10.1038/s41467-019-10669-0 Yang S Wei J Cui Y.-H Park G Shah P Deng Y Aplin A. E Lu Z Hwang S He C He Y.-Y. 2019 m6A mRNA demethylase FTO regulates melanoma tumorigenicity and response to anti-PD-1 blockade Nature Communications 10 2782 https://doi.org/10.1038/s41467-019-10669-0 Zaccara, S. and Jaffrey, S. R. 2020. A unified model for the function of YTHDF proteins in regulating m6A-modified mRNA. Cell 181:1582-1595.e18. https://doi.org/10.1016/j.cell.2020.05.012 Zaccara S Jaffrey S. R. 2020 A unified model for the function of YTHDF proteins in regulating m6A-modified mRNA Cell 181 1582 1595.e18 https://doi.org/10.1016/j.cell.2020.05.012 Zhang, C.; Chen, Y.; Sun, B.; Wang, L.; Yang, Y.; Ma, D.; Lv, J.; Heng, J.; Ding, Y.; Xue, Y.; Lu, X.; Xiao, W.; Yang, Y. G. and Liu, F. 2017. m6A modulates haematopoietic stem and progenitor cell specification. Nature 549:273-276. https://doi.org/10.1038/nature23883 Zhang C. Chen Y. Sun B. Wang L. Yang Y. Ma D. Lv J. Heng J. Ding Y. Xue Y. Lu X. Xiao W. Yang Y. G. Liu F. 2017 m6A modulates haematopoietic stem and progenitor cell specification Nature 549 273 276 https://doi.org/10.1038/nature23883 Zhang, C.; Samanta, D.; Lu, H.; Bullen, J. W.; Zhang, H.; Chen, I.; He, X. and Semenza, G. L. 2016. Hypoxia induces the breast cancer stem cell phenotype by HIF-dependent and ALKBH5-mediated m6A-demethylation of NANOG mRNA. Proceedings of the National Academy of Sciences 113:E2047-E2056. https://doi.org/10.1073/pnas.1602883113 Zhang C. Samanta D. Lu H. Bullen J. W. Zhang H. Chen I. He X. Semenza G. L 2016 Hypoxia induces the breast cancer stem cell phenotype by HIF-dependent and ALKBH5-mediated m6A-demethylation of NANOG mRNA Proceedings of the National Academy of Sciences 113 E2047 E2056 https://doi.org/10.1073/pnas.1602883113 Zhang, L.; Tran, N. T.; Su, H.; Wang, R.; Lu, Y.; Tang, H.; Aoyagi, S.; Guo, A.; Khodadadi-Jamayran, A.; Zhou, D.; Qian, K.; Hricik, T.; Côté, J.; Han, X.; Zhou, W.; Laha, S.; Abdel-Wahab, O.; Levine, R. L.; Raffel, G.; Liu, Y.; Chen, D.; Li, H.; Townes, T.; Wang, H.; Deng, H.; Zheng, Y. G.; Leslie, C.; Luo, M. and Zhao, X. 2015. Cross-talk between PRMT1-mediated methylation and ubiquitylation on RBM15 controls RNA splicing. Elife 4:e07938. https://doi.org/10.7554/eLife.07938 Zhang L. Tran N. T. Su H. Wang R. Lu Y. Tang H. Aoyagi S. Guo A. Khodadadi-Jamayran A. Zhou D. Qian K. Hricik T. Côté J. Han X. Zhou W. Laha S. Abdel-Wahab O. Levine R. L. Raffel G. Liu Y. Chen D. Li H. Townes T. Wang H. Deng H. Zheng Y. G. Leslie C. Luo M. Zhao X. 2015 Cross-talk between PRMT1-mediated methylation and ubiquitylation on RBM15 controls RNA splicing Elife 4 e07938 https://doi.org/10.7554/eLife.07938 Zhang, X.; Yao, Y.; Han, J.; Yang, Y.; Chen, Y.; Tang, Z. and Gao, F. 2020. Longitudinal epitranscriptome profiling reveals the crucial role of N6-methyladenosine methylation in porcine prenatal skeletal muscle development. Journal of Genetics and Genomics 47:466-476. https://doi.org/10.1016/j.jgg.2020.07.003 Zhang X. Yao Y. Han J. Yang Y. Chen Y. Tang Z. Gao F. 2020 Longitudinal epitranscriptome profiling reveals the crucial role of N6-methyladenosine methylation in porcine prenatal skeletal muscle development Journal of Genetics and Genomics 47 466 476 https://doi.org/10.1016/j.jgg.2020.07.003 Zhao, H.; Pan, S.; Duan, J.; Liu, F.; Li, G.; Liu, D. and Liu, Z. 2021. Integrative analysis of m6A regulator-mediated RNA methylation modification patterns and immune characteristics in lupus nephritis. Frontiers in Cell and Developmental Biology 9:724837. https://doi.org/10.3389/fcell.2021.724837 Zhao H. Pan S. Duan J. Liu F. Li G. Liu D. Liu Z. 2021 Integrative analysis of m6A regulator-mediated RNA methylation modification patterns and immune characteristics in lupus nephritis Frontiers in Cell and Developmental Biology 9 724837 https://doi.org/10.3389/fcell.2021.724837 Zhao, P.; Yu, Y.; Feng, W.; Du, H.; Yu, J.; Kang, H.; Zheng, X.; Wang, Z.; Liu, G. E.; Ernst, C. W.; Ran, X.; Wang, J. and Liu, J.-F. 2018. Evidence of evolutionary history and selective sweeps in the genome of Meishan pig reveals its genetic and phenotypic characterization. GigaScience 7:giy058. https://doi.org/10.1093/gigascience/giy058 Zhao P Yu Y Feng W Du H Yu J Kang H Zheng X Wang Z Liu G. E Ernst C. W Ran X Wang J Liu J.-F. 2018 Evidence of evolutionary history and selective sweeps in the genome of Meishan pig reveals its genetic and phenotypic characterization GigaScience 7 giy058 https://doi.org/10.1093/gigascience/giy058 Zheng, G.; Dahl, J. A.; Niu, Y.; Fedorcsak, P.; Huang, C.-M.; Li, C. J.; Vågbø, C. B.; Shi, Y.; Wang, W.-L.; Song, S.-H.; Lu, Z.; Bosmans, R. P. G.; Dai, Q.; Hao, Y. J.; Yang, X.; Zhao, W. M.; Tong, W. M.; Wang, X. J.; Bogdan, F.; Furu, K.; Fu, Y.; Jia, G.; Zhao, X.; Liu, J.; Krokan, H. E.; Klungland, A.; Yang, Y. G. and He, C. 2013. ALKBH5 is a mammalian RNA demethylase that impacts rna metabolism and mouse fertility. Molecular Cell 49:18-29. https://doi.org/10.1016/j.molcel.2012.10.015 Zheng G Dahl J. A Niu Y Fedorcsak P Huang C.-M Li C. J Vågbø C. B Shi Y Wang W.-L Song S.-H Lu Z Bosmans R. P. G Dai Q Hao Y. J Yang X Zhao W. M Tong W. M Wang X. J Bogdan F Furu K Fu Y Jia G Zhao X Liu J Krokan H. E Klungland A Yang Y. G He C 2013 ALKBH5 is a mammalian RNA demethylase that impacts rna metabolism and mouse fertility Molecular Cell 49 18 29 https://doi.org/10.1016/j.molcel.2012.10.015 Zheng, X.; Zhao, P.; Yang, K.; Ning, C.; Wang, H.; Zhou, L. and Liu, J. 2020. CNV analysis of Meishan pig by next-generation sequencing and effects of AHR gene CNV on pig reproductive traits. Journal of Animal Science and Biotechnology 11:42. https://doi.org/10.1186/s40104-020-00442-5 Zheng X. Zhao P. Yang K. Ning C. Wang H. Zhou L. Liu J. 2020 CNV analysis of Meishan pig by next-generation sequencing and effects of AHR gene CNV on pig reproductive traits Journal of Animal Science and Biotechnology 11 42 https://doi.org/10.1186/s40104-020-00442-5 Zhou, R.; Li, S. T.; Yao, W. Y.; Xie, C. D.; Chen, Z.; Zeng, Z. J.; Wang, D.; Xu, K.; Shen, Z. J.; Mu, Y.; Bao, W.; Jiang, W.; Li, R.; Liang, Q. and Li, K. 2021. The Meishan pig genome reveals structural variation‐mediated gene expression and phenotypic divergence underlying Asian pig domestication. Molecular Ecology Resources 21:2077-2092. https://doi.org/10.1111/1755-0998.13396 Zhou R. Li S. T. Yao W. Y. Xie C. D. Chen Z. Zeng Z. J. Wang D. Xu K. Shen Z. J. Mu Y. Bao W. Jiang W. Li R. Liang Q. Li K. 2021 The Meishan pig genome reveals structural variation‐mediated gene expression and phenotypic divergence underlying Asian pig domestication Molecular Ecology Resources 21 2077 2092 https://doi.org/10.1111/1755-0998.13396