m6A mRNA methylation analysis provides novel insights into heat stress responses in the liver tissue of sheep

Adaptive integumentary features of beef cattle raised on afforested or non-shaded tropical pastures

We aimed to analyze the seasonal acclimatization process of Nelore and Canchim cattle raised on two production systems (non-shaded, NS, and integrated crop-livestock-forest, ICLF), based on the dynamics of the morphological and functional attributes of the hair coat and skin during winter and summer. The study was conducted in Brazil, in a low-altitude tropical climate region. A completely randomized 2 × 2 factorial design was adopted as follows: two production systems (NS and ICLF), two breeds (Nelore and Canchim) in a longitudinal structure, with measurements repeated over time through two stations (winter and summer). The experimental animals consisted of 32 Nelore (Bos indicus) and 32 Canchim (5/8 Bos taurus × 3/8 Bos indicus) bulls. The animals were equally distributed between two intensive rotational grazing systems. In both breeds, the hair coat was significantly thicker in winter but longer in summer, which increased epidermal protection. The Nelore bulls had shorter, wider, and thicker hairs, which are attributes that promote heat loss via conduction. The Canchim bulls showed significantly lower hair density and higher epithelium distance to sweat glands, which resulted in higher core temperature and respiratory rate. In turn, Nelore bulls had higher serum concentrations of triiodothyronine and lower serum concentrations of cortisol. However, Canchim bulls more frequently and intensely activated their thermoregulatory system and markedly adjusted their hair coat and hair features to reduce heat gain, especially in summer. Therefore, the anatomical plasticity and functional integumentary characteristics of Nelore and Canchim bulls reflect their acclimatization to tropical conditions.

Exploring candidate genes for heat tolerance in ovine through liver gene expression

Thermotolerance has become an essential factor in the prevention of the adverse effects of heat stress, but it varies among animals. Identifying genes related to heat adaptability traits is important for improving thermotolerance and for selecting more productive animals in hot environments. The primary objective of this research was to find candidate genes in the liver that play a crucial role in the heat stress response of Santa Ines sheep, which exhibit varying levels of heat tolerance. To achieve this goal, 80 sheep were selected based on their thermotolerance and placed in a climate chamber for 10 days, during which the average temperature was maintained at 36 °C from 10 a.m. to 4 p.m. and 28 °C from 4 p.m. to 10 a.m. A subset of 14 extreme animals, with seven thermotolerant and seven non-thermotolerant animals based on heat loss (rectal temperature), were selected for liver sampling. RNA sequencing and differential gene expression analysis were performed. Thermotolerant sheep showed higher expression of genes GPx3, RGS6, GPAT3, VLDLR, LOC101108817, and EVC. These genes were mainly related to the Hedgehog signaling pathway, glutathione metabolism, glycerolipid metabolism, and thyroid hormone synthesis. These enhanced pathways in thermotolerant animals could potentially mitigate the negative effects of heat stress, conferring greater heat resistance.

m6A mRNA Methylation Analysis Provides Novel Insights into Pigmentation in Sheep Skin

N6-methyladenosine (m⁶A) is the most universal post-transcriptional modification of mRNA which may play important roles in verious species. However, the potential roles of m⁶A in the pigmentation of skin are not completely understood. To explore the role of m6A modification in pigmentation of sheep skin, we used MeRIP-seq and RNA-seq to profile the skin transcriptome in black and white coat color (n=3). Our results showed that an average of 7701 m⁶A peaks were obtained for all samples and the average length was 305.89 bp. The GGACUU sequence was the most enrichment motif and shared in black skin and white skin. The m⁶A peaks were mainly enriched in the CDS, 3'UTR and 5'UTR, especially in CDS region near the stop codon of the transcript. 235 significantly differential peaks were found in black skin vs. white skin. The KEGG signaling pathways of downregulated and upregulated m⁶A peaks were mainly enriched in AGE-RAGE signaling pathway in diabetic complications, Viral carcinogenesis, Transcriptional misregulation in cancer, ABC transporters, Basal transcription factors and Thyroid hormone synthesis (P value <0.05). For RNA-seq, 71 differently expressed genes (DEGs) were scanned in black skin vs. white skin. DEGs were significantly enriched in tyrosine metabolism, melanogenesis, neuroactive ligand-receptor interaction pathway (P value <0.05). Combined m⁶A-seq and RNA-seq analysis showed that the hyper-up genes and hypo-up genes were both enriched in ErbB signaling pathway (P value <0.05). In conclusion, it provide a basis for further research into the functions of m⁶A methylation modifications in pigmentation.

Biological methylation: redefining the link between genotype and phenotype

The central dogma of molecular biology is responsible for the crucial flow of genetic information from DNA to protein through the transcription and translation process. Although the sequence of DNA is constant in all organs, the difference in protein and variation in the phenotype is mainly due to the quality and quantity of tissue-specific gene expression and methylation pattern. The term methylation has been defined and redefined by various scientists in the last fifty years. There is always huge excitement around this field because the inheritance of something is beyond its DNA sequence. Advanced gene methylation studies have redefined molecular genetics and these tools are considered de novo in alleviating challenges of animal disease and production. Recent emerging evidence has shown that the impact of DNA, RNA, and protein methylation is crucial for embryonic development, cell proliferation, cell differentiation, and phenotype production. Currently, many researchers are focusing their work on methylation to understand its significant role in expression, disease-resistant traits, productivity, and longevity. The main aim of the present review is to provide an overview of DNA, RNA, and protein methylation, current research output from different sources, methodologies, factors responsible for methylation of genes, and future prospects in animal genetics.

Expression profiling of N6-methyladenosine-modified mRNA in PC12 cells in response to unconjugated bilirubin

BACKGROUND

Abnormal methylation of N6-methyladenosine (m6A) is reportedly associated with central nervous system disorders. However, the role of m6A mRNA methylation in unconjugated bilirubin (UCB) neurotoxicity requires further research.

METHODS

Rat pheochromocytoma PC12 cells treated with UCB were used as in vitro models. After the PC12 cells were treated with UCB (0, 12, 18, and 24 µM) for 24 h, the total RNA m6A levels were measured using an m6A RNA methylation quantification kit. The expression of m6A demethylases and methyltransferases was detected through western blotting. We determined the m6A mRNA methylation profile in PC12 cells exposed to UCB (0 and 18 µM) for 24 h using methylated RNA immunoprecipitation sequencing (MeRIP-seq).

RESULTS

Compared with the control group, UCB (18 and 24 µM) treatment decreased the expression of the m6A demethylase ALKBH5 and increased the expression of the methyltransferases METTL3 and METTL14, which resulted in an increase in the total m6A levels in PC12 cells. Furthermore, 1533 m6A peaks were significantly elevated and 1331 peaks were reduced in the UCB (18 µM)-treated groups compared with those in the control group. Genes with differential m6A peaks were mainly enriched in protein processing in the endoplasmic reticulum, ubiquitin-mediated proteolysis, cell cycle, and endocytosis. Through combined analysis of the MeRIP-seq and RNA sequencing data, 129 genes with differentially methylated m6A peaks and differentially expressed mRNA levels were identified.

CONCLUSION

Our study suggests that the modulation of m6A methylation modifications plays a significant role in UCB neurotoxicity.

Transcriptome-wide high-throughput m6 A sequencing of differential m6 A methylation patterns in the decidual tissues from RSA patients

Recurrent spontaneous abortion (RSA) is characterized by two or more consecutive pregnancy losses in the first trimester of pregnancy, experienced by 5% of women during their reproductive age. As a complex pathological process, the etiology of RSA remains poorly understood. Recent studies have established that gene expression changes dramatically in human endometrial stromal cells (ESCs) during decidualization. N6-methyladenosine (m⁶ A) modification is the most prevalent epigenetic modification of mRNA in eukaryotic cells and it is closely related to the occurrence and development of many pathophysiological phenomena. In this study, we first confirmed that high levels of m⁶ A mRNA methylation in decidual tissues are associated with RSA. Then, we used m⁶ A-modified RNA immunoprecipitation sequence (m⁶ A-seq) and RNA sequence (RNA-seq) to identify the differentially expressed m⁶ A methylation in decidual tissues from RSA patients and identified the key genes involved in abnormal decidualization by bioinformatics analysis. Using m⁶ A-seq, we identified a total of 2169 genes with differentially expressed m⁶ A methylation, of which 735 m⁶ A hypermethylated genes and 1434 m⁶ A hypomethylated genes were identified. Further joint analysis of m⁶ A-seq and RNA-seq revealed that 133 genes were m⁶ A modified with mRNA expression. GO and KEGG analyses indicated that these unique genes were mainly enriched in environmental information processing pathways, including the cytokine-cytokine receptor interaction and PI3K-Akt signaling pathway. In summary, this study uncovered the transcriptome-wide m⁶ A modification pattern in decidual tissue of RSA, which provides a theoretical basis for further research into m⁶ A modification and new therapeutic strategies for RSA.

The alteration of N6-methyladenosine (m6A) modification at the transcriptome-wide level in response of heat stress in bovine mammary epithelial cells

Background

Heat stress has a substantial negative economic impact on the dairy industry. N6-methyladenosine (m6A) is the most common internal RNA modification in eukaryotes and plays a key role in regulating heat stress response in animals. In dairy cows, however, this modification remains largely unexplored. Therefore, we examined the effects of heat stress on the m6A modification and gene expression in bovine mammary epithelial cells to elucidate the mechanism of heat stress response. In this study, Mammary alveolar cells-large T antigen (MAC-T) cells were incubated at 37 °C (non-heat stress group, NH) and 40 °C (heat stress group, H) for 2 hours, respectively. HSP70, HSF1, BAX and CASP3 were up regulated in H group compared with those in the NH group.

Results

Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were conducted to identify m6A peaks and to produce gene expression data of MAC-T cells in the H and NH groups. In total, we identified 17,927 m6A peaks within 9355 genes in the H group, and 18,974 peaks within 9660 genes in the NH groups using MeRIP-seq. Compared with the NH group, 3005 significantly differentially enriched m6A peaks were identified, among which 1131 were up-regulated and 1874 were down-regulated. In addition, 1502 significantly differentially expressed genes were identified using RNA-seq, among which 796 were up-regulated and 706 were down-regulated in the H group compared to the NH group. Furthermore, 199 differentially expressed and synchronously differentially methylated genes were identified by conjoint analysis of the MeRIP-seq and RNA-seq data, which were subsequently divided into four groups: 47 hyper-up, 53 hyper-down, 59 hypo-up and 40 hypo-down genes. In addition, GO enrichment and KEGG analyses were used to analyzed the potential functions of the genes in each section.

Conclusion

The comparisons of m6A modification patterns and conjoint analyses of m6A modification and gene expression profiles suggest that m6A modification plays a critical role in the heat stress response by regulating gene expression.

Effects of High Heat Load Conditions on Blood Constituent Concentrations in Dorper, Katahdin, and St. Croix Sheep from Different Regions of the USA

Forty-six Dorper (DOR), 46 Katahdin (KAT), and 43 St. Croix (STC) female sheep (initial body weight of 58, 59, and 46 kg, respectively, SEM = 1.75; 3.3 ± 0.18 years of age, 2.6−3.7), derived from 45 commercial farms in four regions of the USA (Midwest, Northwest, Southeast, and central Texas), were used to evaluate responses in blood constituent concentrations to increasing heat load index (HLI) conditions. There were four sequential 2 weeks periods with target HLI during day/nighttime of 70/70 (thermoneutral zone conditions), 85/70, 90/77, and 95/81 in period 1, 2, 3 and 4, respectively. A 50% concentrate pelletized diet was fed at 53.3 g dry matter/kg body weight0.75. The analysis of most constituents was for samples collected on the last day of the second week of each period at 13:00 h; samples for cortisol, thyroxine, and heat shock protein were collected in week 2 and 8. Previously, it was noted that resilience to high HLI conditions was greatest for STC, lowest for DOR, and intermediate for KAT. There were few effects of region. Other than hemoglobin concentration, there were no interactions between breed and period. Blood oxygen concentration was greatest (p < 0.05) among breeds for STC (5.07, 5.20, and 5.53 mmol/L for DOR, KAT, and STC, respectively; SEM = 0.114) and differed among periods (4.92, 5.26, 5.36, and 5.52 mmol/L for period 1, 2, 3, and 4, respectively; SEM = 0.093). There were breed differences (i.e., main effects; p < 0.05) in glucose (50.0, 52.6, and 52.1 mg/dL; SEM = 0.76), urea nitrogen (17.2, 17.3, and 19.4 mg/dL; SEM = 0.33), creatinine (0.991, 0.862, and 0.802 mg/dL; SEM = 0.0151), total protein (6.50, 6.68, and 6.95 g/l; SEM = 0.017), triglycerides (28.4, 29.1, and 23.5 mg/dL; SEM = 0.87), and cortisol (6.30, 8.79, and 6.22 ng/mL for DOR, KAT, and STC, respectively; SEM = 0.596). Differences among periods (p < 0.05) were observed for lactate (27.9, 25.3, 27.8, and 24.0 mg/dL; SEM = 0.99), creatinine (0.839, 0.913, 0.871, and 0.917 mg/dL; SEM = 0.0128), total protein (6.94, 6.66, 6.60, and 6.65 g/l; SEM = 0.094), and cholesterol (60.2, 56.5, 58.3, and 57.6 mg/dL for period 1, 2, 3, and 4, respectively; SEM = 1.26). In addition, the concentration of cortisol (7.62 and 6.59 ng/mL; SEM = 0.404), thyroxine (5.83 and 5.00 µg/dL; SEM = 0.140), and heat shock protein (136 and 146 ng/mL for week 2 and 8, respectively; SEM = 4.0) differed between weeks (p < 0.05). In conclusion, the lack of interaction between breed and period with different HLI conditions suggests that levels of these blood constituents were not highly related to resilience to high HLI.

Identifying the heat resistant genes by multi-tissue transcriptome sequencing analysis in Turpan Black sheep

Heat stress not only affects the physical condition but also affects reproductive performance in sheep. A thorough understanding of the molecular and physiological mechanisms underlying heat stress would certainly improve livestock productivity and provide genetic evaluation ways for heat resistant breeds selection. In this study, 85 Turpan Black sheep, a breed exhibited excellent heat resistance from long-term artificial selection, and 85 heat sensitive Kazakh sheep in Turpan basin were tested for physiological and reproductive performance from July to August in summer. The results showed that the estrus rate was significantly higher in Turpan Black sheep (P < 0.05), while the heart rate and respiratory rate of Turpan Black sheep are significantly lower than that of Kazakh sheep (P < 0.05). Furthermore, to clarify genes participated in heat stress response, the pituitary, ovarian and hepatic tissues from three Turpan Black sheep and three Kazakh sheep were subjected to RNA-seq. The results indicated that 32, 49 and 69 genes were up-regulated, and 39, 60 and 145 genes were down-regulated in pituitary, ovarian and hepatic tissues in Turpan Black sheep compared with that of the Kazakh sheep, respectively. KEGG and gene set enrichment analysis showed that the differentially expressed genes were mainly involved in signal transduction pathways. In particular, the differentially expressed genes in hepar were enriched in the energy metabolism pathway, while the differentially expressed genes in ovarian tissue were enriched in the ovarium steroidogenesis pathway. In conclusion, our results implied that the pituitary-ovary axis might include hepar as downstream targeted organism in heat resistant regulation. Under heat stress, the signals released from pituitary would impact steroidogenesis in ovary, and further alter energy metabolism in hepar. As we know, this is the first comparative study to investigate the gene expression in multi-tissue in sheep under heat stress.

N6-Methyladenosine Methylation Analysis of Long Noncoding RNAs and mRNAs in IPEC-J2 Cells Treated With Clostridium perfringens beta2 Toxin

Background

The n6-methyladenosine (m6A) modification is present widely in mRNAs and long non-coding RNAs (lncRNAs), and is related to the occurrence and development of certain diseases. However, the role of m6A methylation in Clostridium perfringens type C infectious diarrhea remains unclear.

Methods

Here, we treated intestinal porcine jejunum epithelial cells (IPEC-J2 cells) with Clostridium perfringens beta2 (CPB2) toxin to construct an in vitro model of Clostridium perfringens type C (C. perfringens type C) infectious diarrhea, and then used methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) to identify the methylation profiles of mRNAs and lncRNAs in IPEC-J2 cells.

Results

We identified 6,413 peaks, representing 5,825 m6A-modified mRNAs and 433 modified lncRNAs, of which 4,356 m6A modified mRNAs and 221 m6A modified lncRNAs were significantly differential expressed between the control group and CPB2 group. The motif GGACU was enriched significantly in both the control group and the CPB2 group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analysis showed that the differentially methylated modified mRNAs were mainly enriched in Hippo signaling pathway and Wnt signaling pathway. In addition, the target genes of the differentially m6A modified lncRNAs were related to defense response to virus and immune response. For example, ENSSSCG00000042575, ENSSSCG00000048701 and ENSSSCG00000048785 might regulate the defense response to virus, immune and inflammatory response to resist the harmful effects of viruses on cells.

Conclusion

In summary, this study established the m6A transcription profile of mRNAs and lncRNAs in IPEC-J2 cells treated by CPB2 toxin. Further analysis showed that m6A-modified RNAs were related to defense against viruses and immune response after CPB2 toxin treatment of the cells. Threem6A-modified lncRNAs, ENSSSCG00000042575, ENSSSCG00000048785 and ENSSSCG00000048701, were most likely to play a key role in CPB2 toxin-treated IPEC-J2 cells. The results provide a theoretical basis for further research on the role of m6A modification in piglet diarrhea.

N 6-Methyladenosine Methylomic Landscape of Lung Tissues in Murine Acute Allergic Asthma

Allergic asthma is well known as a common respiratory disorder comprising an allergic inflammatory nature and excessive immune characteristic. N ⁶-methyladenosine (m6A) methylation is an RNA epigenetic modification that post-transcriptionally regulates gene expression and function by affecting the RNA fate. Currently, m6A methylation is gaining attention as a mechanism of immunoregulation. However, whether m6A methylation engages the pathological process of asthma remains uncertain. Here, we present the m6A methylomic landscape in the lung tissues of ovalbumin-induced acute asthma mice using MeRIP-seq and RNA-seq. We identified 353 hypermethylated m6A peaks within 329 messenger RNAs (mRNAs) and 150 hypomethylated m6A peaks within 143 mRNAs in the lung tissues of asthmatic mice. These differentially methylated mRNAs were found to be involved in several immune function-relevant signaling pathways. In addition, we predicted 25 RNA-binding proteins that recognize the differentially methylated peak sites by exploring public databases, and the roles of these proteins are mostly related to mRNA biogenesis and metabolism. To further investigate the expression levels of the differentially methylated genes, we performed combined analysis of the m6A methylome and transcriptome data and identified 127 hypermethylated mRNAs (107 high and 20 low expression) and 43 hypomethylated mRNAs with differential expressions (9 high and 34 low expression). Of these, there are a list of mRNAs involved in immune function and regulation. The present results highlight the essential role of m6A methylation in the pathogenesis of asthma.

Gamma-irradiation fluctuates the mRNA N6-methyladenosine (m6A) spectrum of bone marrow in hematopoietic injury

Humans benefit from nuclear technologies but consequently experience nuclear disasters or side effects of iatrogenic radiation. Hematopoietic system injury first arises upon radiation exposure. As an intricate new layer of genetic control, the posttranscriptional m⁶A modification of RNA has recently come under investigation and has been demonstrated to play pivotal roles in multiple physiological and pathological processes. However, how the m⁶A methylome functions in the hematopoietic system after irradiation remains ambiguous. Here, we uncovered the time-varying epitranscriptome-wide m⁶A methylome and transcriptome alterations in γ-ray-exposed mouse bone marrow. 4 Gy γ-irradiation rapidly (5 min and 2 h) and severely impaired the mouse hematopoietic system, including spleen and thymus weight, blood components, tissue inflammation and malondialdehyde (MDA) levels. The m⁶A content and expression of m⁶A related enzymes were altered. Gamma-irradiation triggered dynamic and reversible m⁶A modification profiles and altered mRNA expression, where both m⁶A fold-enrichment and mRNA expression most followed the (5 min_up/2 h_down) pattern. The CDS enrichment region preferentially upregulated m⁶A peaks at 5 min. Moreover, the main GO and KEGG pathways were closely related to metabolism and the classical radiation response. Finally, m⁶A modifications correlated with transcriptional regulation of genes in multiple aspects. Blocking the expression of m⁶A demethylases FTO and ALKBH5 mitigated radiation hematopoietic toxicity. Together, our findings present the comprehensive landscape of mRNA m⁶A methylation in the mouse hematopoietic system in response to γ-irradiation, shedding light on the significance of m⁶A modifications in mammalian radiobiology. Regulation of the epitranscriptome may be exploited as a strategy against radiation damage.

METTL3-Mediated lncRNA m6A Modification in the Osteogenic Differentiation of Human Adipose-Derived Stem Cells Induced by NEL-Like 1 Protein

OBJECTIVES

This study aimed to explore the regulatory mechanism of methyltransferase3 (METTL3) -mediated long non-coding RNA (lncRNA) N6-methyladenosine (m⁶A) modification in the osteogenic differentiation of human adipose-derived stem cells (hASCs) induced by NEL-like 1 protein (NELL-1).

MATERIALS AND METHODS

Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and high- throughput sequencing for RNA (RNA-seq) were performed on hASCs. Osteogenic ability was detected by alkaline phosphatase (ALP) staining, Alizarin Red S(ARS) staining, ALP quantification and Quantitative real-time polymerase chain reaction analysis (qRT-PCR). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis predicted the osteogenesis-related pathways enriched for the lncRNAs and identified the target lncRNAs. After overexpression and knockdown of METTL3, methylated RNA immunoprecipitation-qPCR (MeRIP-qPCR) and qRT-PCR were used to detect the levels of m⁶A modification and the expression of the target lncRNA, and the binding of both was confirmed by RNA binding protein immunoprecipitation (RIP) assay. The effects of lncRNA and METTL3 on phosphorylation of the key proteins of the pathway were detected by western blot analysis.

RESULTS

In vitro experiments showed that METTL3 can promote osteogenic differentiation and that its expression level is upregulated. KEGG pathway analysis predicted that lncRNAs with differentially upregulated methylated peaks were enriched mostly in the mitogen-activated protein kinase (MAPK) signaling pathway, in which Serine/threonine protein kinase 3 (STK3) was the predicted target gene of the lncRNA RP11-44 N12.5. The m⁶A modification and expression of RP11-44 N12.5 were both regulated by METTL3. Subsequently, lncRNA RP11-44 N12.5 and METTL3 were found to regulate the phosphorylation levels of three key proteins in the MAPK signaling pathway, ERK, JNK and p38.

CONCLUSIONS

This study shows, for the first time, that METTL3 can activate the MAPK signaling pathway by regulating the m⁶A modification and expression of a lncRNA, thereby enhancing the osteogenic differentiation of hASCs.

Comprehensive analysis of the transcriptome-wide m6A methylome in invasive malignant pleomorphic adenoma

Background

Invasive malignant pleomorphic adenoma (IMPA) is a highly invasive parotid gland tumor and lacks effective therapy. N6-Methyladenosine (m⁶A) is the most prevalent post-transcriptional modification of mRNAs in eukaryotes and plays an important role in the pathogenesis of multiple tumors. However, the significance of m⁶A-modified mRNAs in IMPA has not been elucidated to date. Hence, in this study, we attempted to profile the effect of IMPA in terms of m⁶A methylation in mRNA.

Methods

Methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) were utilized to acquire the first transcriptome-wide profiling of the m⁶A methylome map in IMPA followed by bioinformatics analysis.

Results

In this study, we obtained m⁶A methylation maps of IMPA samples and normal adjacent tissues through MeRIP-seq. In total, 25,490 m⁶A peaks associated with 13,735 genes were detected in the IMPA group, whereas 33,930 m⁶A peaks associated with 18,063 genes were detected in the control group. Peaks were primarily enriched within coding regions and near stop codons with AAACC and GGAC motifs. Moreover, functional enrichment analysis demonstrated that m⁶A-containing genes were significantly enriched in cancer and metabolism relevant pathways. Furthermore, we identified a relationship between the m⁶A methylome and the RNA transcriptome, indicating a mechanism by which m⁶A modulates gene expression.

Conclusions

Our study is the first to provide comprehensive and transcriptome-wide profiles to determine the potential roles played by m⁶A methylation in IMPA. These results may open new avenues for in-depth research elucidating the m⁶A topology of IMPA and the molecular mechanisms governing the formation and progression of IMPA.