Abstract
Rosacea is a common skin disease that predominantly affects individuals aged between 30 and 50 years. While the exact cause of the disease remains unclear, various factors have been shown to trigger or exacerbate its symptoms. N6-methyladenosine (m(6)A) modification is one of the most abundant epigenetic methylation modification in messenger RNA (mRNA) and non-coding RNA (ncRNA), plays a crucial role in RNA splicing, export, stability, and translation. In this study, we aimed to characterize m(6)A genes in rosacea, identify molecular subtypes based on m(6)A gene expression, characterize the immune features among subtypes, explore key molecules based on co-expression analysis, and identify potential targets and drugs. To achieve our objectives, we first compared the expression pattern and immune regulation of m(6)A genes between healthy and diseased groups. Then, we performed clustering to stratify disease samples into different subtypes and analyzed immune regulation and functional enrichment among the subtypes. Then, we conducted differential analysis between subtypes and applied weighted gene co-expression network analysis (WGCNA) in three subtypes. We found hub differential expression analysis (DEG) genes and their potential drug based on the WGCNA results and the drug-gene interaction database (DGIdb). Finally, in vivo and in vitro studies showed significant differences in m(6)A methyltransferase METTL3 levels in rosacea mice and control mice, as well as in the skin of rosacea patients and healthy people, while reducing METTL3 significantly inhibited the inflammatory response of human fibroblasts (HDFs) stimulated by LL37, suggesting that METTL3 may be associated with changes in overall m(6)A levels in rosacea. Taken together, our findings provide valuable insights into therapeutic targets and drug predictions for rosacea.