Diagnostic and predictive values of m5C‑associated genes in idiopathic pulmonary fibrosis

In patients with idiopathic pulmonary fibrosis (IPF), the role of 5‑methylcytosine (m5C)‑associated genes in the pathogenesis and development of the disease remains unclear. The present study aimed to identify reliable diagnostic markers based on the expression of m5C‑associated genes for the early detection of IPF. Count data were obtained by screening the IPF genome‑wide assay in the Gene Expression Omnibus database, followed by a comparison of m5C gene expression in patients with IPF and controls. The GSE150910 and GSE173355 datasets yielded a total of 23 differentially expressed m5C‑associated genes, which were then investigated for their functions. A diagnostic model was built using eight m5C genes and validated with training sets and the GSE124685 dataset. IPF subtypes were identified based on expression of m5C‑related genes as well as clinical and immunological characteristics. Furthermore, a pulmonary fibrosis model was established in mice by administering bleomycin into the trachea. Lungs were harvested and analyzed using quantitative PCR to determine the expression of m5C‑related genes. The Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that these genes were significantly enriched in 'base excision repair'. Immunoassay results revealed that 13 immune cell markers (naive, memory and B cell plasma, T cell CD4 naive, T cell CD4 memory resting, T cell follicular helper, T cell regulatory Tregs, NK cell resting, Monocyte, Macrophage M0, Mast cell activated, Eosinophil, and Neutrophil) were significantly associated with IPF. Patients with IPF had lower levels of resting memory CD4+ T cells, which were positively associated with Tet methylcytosine dioxygenase2 (TET2) and Thymine‑DNA glycosylase (TDG) but negatively correlated with NOP2/Sun RNA methyltransferase5 (NSUN5) expression. All samples were classified into based on the levels of the eight diagnostic m5C genes. Samples with high m5C scores are subtype 1, and those with low m5C scores are subtype 2. In subtype 2, male patients had lower levels of CD27 and CD70 but higher levels of CD274, CD86, Cytotoxic T‑lymphocyte‑associated protein4 and Hepatitis A virus cellular receptor2 (HAVCR2). When compared with normal mouse lung tissue samples, expression levels of NOP2/Sun RNA methyltransferase6 (NSUN6), Ubiquitin‑like with PHD and RING Finger Domains1, TDG and TET2 in lung fibrosis tissue samples were significantly higher, while expression levels of NSUN5, NTH‑like DNA glycosylase1, DNA (cytosine‑5‑)‑methyltransferase3 β and Methyl‑CpG binding domain protein 3) were lower. It is possible that m5C‑associated genes play an important role in the diagnosis and typing of IPF. These genes may facilitate investigation of the pathophysiology of IPF and identification of potential treatment targets.