Analysis of mRNA m(6)A modification and mRNA expression profiles in middle ear cholesteatoma

Introduction: Middle ear cholesteatoma is characterized by the hyperproliferation of keratinocytes. In recent decades, N(6)-methyladenosine (m(6)A) modification has been shown to play an essential role in the pathogenesis of many proliferative diseases. However, neither the m(6)A modification profile nor its potential role in the pathogenesis of middle ear cholesteatoma has currently been investigated. Therefore, this study aimed to explore m(6)A modification patterns in middle ear cholesteatoma. Materials and methods: An m(6)A mRNA epitranscriptomic microarray analysis was performed to analyze m(6)A modification patterns in middle ear cholesteatoma tissue (n = 5) and normal post-auricular skin samples (n = 5). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to predict the potential biological functions and signaling pathways underlying the pathogenesis of middle ear cholesteatoma. Subsequently, m(6)A modification levels were verified by methylated RNA immunoprecipitation–qPCR (MeRIP–qPCR) in middle ear cholesteatoma tissue and normal skin samples, respectively. Results: A total of 6,865 distinctive m(6)A-modified mRNAs were identified, including 4,620 hypermethylated and 2,245 hypomethylated mRNAs, as well as 9,162 differentially expressed mRNAs, including 4,891 upregulated and 4,271 downregulated mRNAs, in the middle ear cholesteatoma group relative to the normal skin group. An association analysis between methylation and gene expression demonstrated that expression of 1,926 hypermethylated mRNAs was upregulated, while expression of 2,187 hypomethylated mRNAs and 38 hypermethylated mRNAs was downregulated. Moreover, GO analysis suggested that differentially methylated mRNAs might influence cellular processes and biological behaviors, such as cell differentiation, biosynthetic processes, regulation of molecular functions, and keratinization. KEGG pathway analysis demonstrated that the hypermethylated transcripts were involved in 26 pathways, including the Hippo signaling pathway, the p53 signaling pathway, and the inflammatory mediator regulation of transient receptor potential (TRP) channels, while the hypomethylated transcripts were involved in 13 pathways, including bacterial invasion of epithelial cells, steroid biosynthesis, and the Hippo signaling pathway. Conclusion: Our study presents m(6)A modification patterns in middle ear cholesteatoma, which may exert regulatory roles in middle ear cholesteatoma. The present study provides directions for mRNA m(6)A modification-based research on the epigenetic etiology and pathogenesis of middle ear cholesteatoma.