Identification of key genes and biological pathways in lung adenocarcinoma via bioinformatics analysis

Comprehensive Profiling and Therapeutic Insights into Differentially Expressed Genes in Hepatocellular Carcinoma

Background: Drug repurposing is a strategy that complements the conventional approach of developing new drugs. Hepatocellular carcinoma (HCC) is a highly prevalent type of liver cancer, necessitating an in-depth understanding of the underlying molecular alterations for improved treatment. Methods: We searched for a vast array of microarray experiments in addition to RNA-seq data. Through rigorous filtering processes, we have identified highly representative differentially expressed genes (DEGs) between tumor and non-tumor liver tissues and identified a distinct class of possible new candidate drugs. Results: Functional enrichment analysis revealed distinct biological processes associated with metal ions, including zinc, cadmium, and copper, potentially implicating chronic metal ion exposure in tumorigenesis. Conversely, up-regulated genes are associated with mitotic events and kinase activities, aligning with the relevance of kinases in HCC. To unravel the regulatory networks governing these DEGs, we employed topological analysis methods, identifying 25 hub genes and their regulatory transcription factors. In the pursuit of potential therapeutic options, we explored drug repurposing strategies based on computational approaches, analyzing their potential to reverse the expression patterns of key genes, including AURKA, CCNB1, CDK1, RRM2, and TOP2A. Potential therapeutic chemicals are alvocidib, AT-7519, kenpaullone, PHA-793887, JNJ-7706621, danusertibe, doxorubicin and analogues, mitoxantrone, podofilox, teniposide, and amonafide. Conclusion: This multi-omic study offers a comprehensive view of DEGs in HCC, shedding light on potential therapeutic targets and drug repurposing opportunities.

Identification of key genes and biological pathways in lung adenocarcinoma by integrated bioinformatics analysis

BACKGROUND

The objectives of this study were to identify hub genes and biological pathways involved in lung adenocarcinoma (LUAD) via bioinformatics analysis, and investigate potential therapeutic targets.

AIM

To determine reliable prognostic biomarkers for early diagnosis and treatment of LUAD.

METHODS

To identify potential therapeutic targets for LUAD, two microarray datasets derived from the Gene Expression Omnibus (GEO) database were analyzed, GSE3116959 and GSE118370. Differentially expressed genes (DEGs) in LUAD and normal tissues were identified using the GEO2R tool. The Hiplot database was then used to generate a volcanic map of the DEGs. Weighted gene co-expression network analysis was conducted to cluster the genes in GSE116959 and GSE118370 into different modules, and identify immune genes shared between them. A protein-protein interaction network was established using the Search Tool for the Retrieval of Interacting Genes database, then the CytoNCA and CytoHubba components of Cytoscape software were used to visualize the genes. Hub genes with high scores and co-expression were identified, and the Database for Annotation, Visualization and Integrated Discovery was used to perform enrichment analysis of these genes. The diagnostic and prognostic values of the hub genes were calculated using receiver operating characteristic curves and Kaplan-Meier survival analysis, and gene-set enrichment analysis was conducted. The University of Alabama at Birmingham Cancer data analysis portal was used to analyze relationships between the hub genes and normal specimens, as well as their expression during tumor progression. Lastly, validation of protein expression was conducted on the identified hub genes via the Human Protein Atlas database.

RESULTS

Three hub genes with high connectivity were identified; cellular retinoic acid binding protein 2 (CRABP2), matrix metallopeptidase 12 (MMP12), and DNA topoisomerase II alpha (TOP2A). High expression of these genes was associated with a poor LUAD prognosis, and the genes exhibited high diagnostic value.

CONCLUSION

Expression levels of CRABP2, MMP12, and TOP2A in LUAD were higher than those in normal lung tissue. This observation has diagnostic value, and is linked to poor LUAD prognosis. These genes may be biomarkers and therapeutic targets in LUAD, but further research is warranted to investigate their usefulness in these respects.

AURKA, TOP2A and MELK are the key genes identified by WGCNA for the pathogenesis of lung adenocarcinoma

The comprehensive analysis of single or multiple microarray datasets is currently available in Gene Expression Omnibus (GEO) databases, with several studies having identified genes strongly associated with the development of lung adenocarcinoma (LUAD). However, the mechanisms of LUAD development remain largely unknown and has not yet been systematically studied; thus, further studies are required in this field. In the present study, weighted gene co-expression network analysis (WGCNA) was used for the evaluation of key genes with potential high risk of LUAD, and to provide more reliable evidence concerning its pathogenesis. The GSE140797 dataset from the high-throughput GEO database was downloaded and was first analyzed using the Limma package in the R language in order to determine the differentially expressed genes. The dataset was then analyzed using the WGCNA package to analyze the co-expressed genes, and the modular genes with the highest correlation with the clinical phenotype were identified. Subsequently, the pathogenic genes shared in common between the result of the two analyses were imported into the STRING database for protein-protein interaction network analysis. The hub genes were screened out using Cytoscape, and then The Cancer Genome Atlas analysis, receiver operating characteristic analysis and survival analysis were subsequently performed. Finally, the key genes were evaluated using reverse transcription-quantitative PCR and western blot analysis. Bioinformatics analysis of the GSE140797 dataset revealed eight key genes: AURKA, BUB1, CCNB1, CDK1, MELK, NUSAP1, TOP2A and PBK. Finally, the AURKA, TOP2A and MELK genes were evaluated in samples from patients with lung cancer using WGCNA and RT-qPCR, western blot analysis experiments, providing basis for further research on the mechanisms of LUAD development and targeted therapy.

Integrative analysis of TP53 mutations in lung adenocarcinoma for immunotherapies and prognosis

BACKGROUND

The TP53 tumor suppressor gene is one of the most mutated genes in lung adenocarcinoma (LUAD) and plays a vital role in regulating the occurrence and progression of cancer. We aimed to elucidate the association between TP53 mutations, response to immunotherapies and the prognosis of LUAD.

METHODS

Genomic, transcriptomic, and clinical data of LUAD were downloaded from The Cancer Genome Atlas (TCGA) dataset. Gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, gene set enrichment analysis (GSEA). Gene set variation analysis (GSVA) were performed to determine the differences in biological pathways. A merged protein-protein interaction (PPI) network was constructed and analyzed. MSIpred was used to analyze the correlation between the expression of the TP53 gene, tumor mutation burden (TMB) and tumor microsatellite instability (MSI). CIBERSORT was used to calculate the abundance of immune cells. Univariate and multivariate Cox regression analyses were used to determine the prognostic value of TP53 mutations in LUAD.

RESULTS

TP53 was the most frequently mutated in LUAD, with a mutational frequency of 48%. GO and KEGG enrichment analysis, GSEA, and GSVA results showed a significant upregulation of several signaling pathways, including PI3K-AKT mTOR (P < 0.05), Notch (P < 0.05), E2F target (NES = 1.8, P < 0.05), and G2M checkpoint (NES = 1.7, P < 0.05). Moreover, we found a significant correlation between T cells, plasma cells, and TP53 mutations (R² < 0.01, P = 0.040). Univariate and multivariate Cox regression analyses revealed that the survival prognosis of LUAD patients was related to TP53 mutations (Hazard Ratio (HR) = 0.72 [95% CI, 0.53 to 0.98], P < 0.05), cancer status (P < 0.05), and treatment outcomes (P < 0.05). Lastly, the Cox regression models showed that TP53 exhibited good power in predicting three- and five-year survival rates.

CONCLUSIONS

TP53 may be an independent predictor of response to immunotherapy in LUAD, and patients with TP53 mutations have higher immunogenicity and immune cell infiltration.

Identification of potential biomarkers for papillary thyroid carcinoma by comprehensive bioinformatics analysis

To perform bioinformatics analysis on the papillary thyroid carcinoma (PTC) gene chip dataset to explore new biological markers for PTC. The gene expression profiles of GSE3467 and GSE6004 chip data were collected by GEO2R, and the differentially expressed genes (DEGs) were selected for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Protein-protein interaction (PPI) relationship analysis was achieved using STRING, and the hub genes were obtained using the Cytoscape software. GEPIA was used to validate the expressions of the hub genes in the normal and tumor tissues and to conduct survival analyses. Pertinent genetic pathology results were fetched using the HPA database. Finally, the key genes were clinically verified by reverse transcription-polymerase chain reaction. 97 genes were jointly up-regulated and 107 genes were jointly down-regulated in GSE3467 and GSE6004. GO function enrichment analysis revealed that the DEGs were involved in the regulation of calcium ion transport into cytosol, integrin binding, and cell adhesion molecule binding. KEGG pathway enrichment analysis indicated that the DEGs were chiefly associated with thyroid cancer and non-small cell lung cancer. According to the PPI network, 30 key target genes were identified. Only the expressions of ANK2, TLE1, and TCF4 matched between the normal and tumor tissues, and were associated with disease prognosis. When compared with the normal thyroid tissues, the protein and mRNA expressions of ANK2, TLE1, and TCF4 were down-regulated in PTC. Significant differences exist in overall gene expression between the thyroid tissues of patients with PTC and those of healthy people. Furthermore, the differential genes ANK2, TLE1, and TCF4 are expected to be reliable molecular markers for the mechanism study and diagnosis of PTC.

Network module function enrichment analysis of lung squamous cell carcinoma and lung adenocarcinoma

BACKGROUND

Lung squamous cell carcinoma (LUSC) and lung adenocarcinoma (LUAD) are the two major subtypes of non-small cell lung cancer that pose a serious threat to human health. However, both subtypes currently lack effective indicators for early diagnosis.

METHODS

To identify tumor-specific indicators and predict cancer-related signaling pathways, LUSC and LUAD gene weighted co-expression networks were constructed. Combined with clinical data, core genes in LUSC and LUAD modules were then screened using protein-protein interaction networks and their functions and pathways were analyzed. Finally, the effect of core genes on survival of LUSC and LUAD patients was evaluated.

RESULTS

We identified 12 network modules in LUSC and LUAD, respectively. LUSC modules "purple" and "green" and LUAD modules "brown" and "pink" are significantly associated with overall survival and clinical traits of tumor node metastasis, respectively. Eleven genes from LUSC and eight genes from LUAD were identified as candidate core genes, respectively. Survival analysis showed that high expression of SLIT3, ABI3BP, MYOCD, PGM5, TNXB, and DNAH9 are associated with decreased survival in LUSC patients. Furthermore, high expression of BUB1, BUB1B, TTK, and UBE2C are associated with lower patient survival.

CONCLUSIONS

We found biomarker genes and biological pathways for LUSC and LUAD. These network hub genes are associated with clinical characteristics and patient outcomes and they may play important roles in LUSC and LUAD.

Enhancer-associated regulatory network and gene signature based on transcriptome and methylation data to predict the survival of patients with lung adenocarcinoma

Accumulating evidence has proved that aberrant methylation of enhancers plays regulatory roles in gene expression for various cancers including lung adenocarcinoma (LUAD). In this study, the transcriptome and methylation data of The Cancer Genome Atlas (TCGA)-LUAD cohort were comprehensively analyzed with a five-step Enhancer Linking by Methylation/Expression Relationships (ELMER) process. Step 1: 131,371 distal (2 kb upstream from the transcription start site) probes were obtained. Step 2: 10,665 distal hypomethylated probes were identified in an unsupervised mode with the get.diff.meth function. Step 3: 699 probe-gene pairs with negative correlations were screened using the get.pair function in an unsupervised mode. Step 4: After mapping with probes, 768 motifs were obtained and 24 of them were enriched. Step 5: 127 transcription factors (TFs) with differential expressions and negative correlations with methylation levels were screened, which were corresponding to 21 motifs. After the ELMER process, a prognostic “TFs-motifs-genes” regulatory network was constructed. The Least absolute shrinkage and selection operator (LASSO) and Stepwise regression analyses were further applied to identify variables in the TCGA-LUAD cohort and an eight-gene signature was constructed for calculating the risk score. The risk score was verified in two independent validation cohorts. The area under curve values of receiver operating characteristic curves predicting 1-, 3-, and 5-years survival ranged from 0.633 to 0.764. With the increase of the risk scores, both the survival statuses and clinical traits showed a worse tendency. There were significant differences in the degrees of immune cell infiltration, TMB values, and TIDE scores between the high-risk and low-risk groups. Finally, a better-performing prognostic nomogram was integrated with the risk score and other clinical traits. In short, this multi-omics analysis demonstrated the application of ELMER in analyzing enhancer-associated regulatory network in LUAD, which provided promising strategies for epigenetic therapy and prognostic biomarkers.

ASPM, CDC20, DLGAP5, BUB1B, CDCA8, and NCAPG May Serve as Diagnostic and Prognostic Biomarkers in Endometrial Carcinoma

Uterine Corpus Endometrial Carcinoma (UCEC), the most common gynecologic malignancy in developed countries, remains to be a major public health problem. Further studies are surely needed to elucidate the tumorigenesis of UCEC. Herein, intersecting 203 differentially expressed genes (DEGs) were identified with the GSE17025, GSE63678, and The Cancer Genome Atlas-UCEC datasets. The Gene Ontology/Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis and protein-protein interaction (PPI) network were performed on those 203 DEGs. Intriguingly, 6 of the top 10 nodes in the PPI network were related to unfavorable prognosis, that is, ASPM, CDC20, DLGAP5, BUB1B, CDCA8, and NCAPG. The mRNA and protein expression levels of the 6 hub genes were elevated in UCEC tissues compared to normal tissues. Higher expression of the 6 hub genes was associated with poor prognostic clinicopathological characteristics. The receiver operating characteristic curve suggested the significant diagnostic ability of the 6 hub genes for UCEC. Then, underlying pathogeneses of UCEC including promoter methylation level, TP53 mutation status, genomic genetic variation, and immune cells infiltration were analyzed. The mRNA expression level of the 6 hub genes was also higher in cervical squamous cell carcinoma and endocervical adenocarcinoma, uterine carcinosarcoma, and ovarian serous cystadenocarcinoma tissues than in corresponding normal tissues. In conclusion, ASPM, CDC20, DLGAP5, BUB1B, CDCA8, and NCAPG may be considered diagnostic and prognostic biomarkers in UCEC.

Let-7c-5p Restrains Cell Growth and Induces Apoptosis of Lung Adenocarcinoma Cells via Targeting ESPL1

Lung adenocarcinoma (LUAD) is a predominant malignancy, and its high mortality prompts us to incessantly probe the relevant targeted treatment. This work intended to study the molecular mechanism of ESPL1 in LUAD. Bioinformatics analysis was performed for pan-cancer and prognosis analysis as well as target gene prediction. Expression of ESPL1 mRNA and let-7c-5p was determined via qRT-PCR, and western blot was employed to detect protein level of ESPL1. Dual-luciferase reporter gene method verified the interaction between ESPL1 and let-7c-5p. Thereafter, CCK-8, wound healing, Transwell, and flow cytometry assays were utilized to investigate proliferation, migration, and apoptosis of LUAD cells. The results revealed that ESPL1 was upregulated in LUAD, which was associated with poor prognosis. Overexpressed ESPL1 promoted LUAD cells to invade, proliferate, and migrate. Furthermore, ESPL1 was a target gene of let-7c-5p. Let-7c-5p was downregulated in LUAD cells, and played a suppressive role in LUAD malignant development, while reversed by ESPL1. Taken together, it was posited that let-7c-5p/ESPL1 may be underlying therapeutic targets of LUAD.

Multiple datasets to explore the tumor microenvironment of cutaneous squamous cell carcinoma

BACKGROUND

Cutaneous squamous cell carcinoma (cSCC) is one of the most frequent types of cutaneous cancer. The composition and heterogeneity of the tumor microenvironment significantly impact patient prognosis and the ability to practice precision therapy. However, no research has been conducted to examine the design of the tumor microenvironment and its interactions with cSCC.

MATERIAL AND METHODS

We retrieved the datasets GSE42677 and GSE45164 from the GEO public database, integrated them, and analyzed them using the SVA method. We then screened the core genes using the WGCNA network and LASSO regression and checked the model's stability using the ROC curve. Finally, we performed enrichment and correlation analyses on the core genes.

RESULTS

We identified four genes as core cSCC genes: DTYMK, CDCA8, PTTG1 and MAD2L1, and discovered that RORA, RORB and RORC were the primary regulators in the gene set. The GO semantic similarity analysis results indicated that CDCA8 and PTTG1 were the two most essential genes among the four core genes. The results of correlation analysis demonstrated that PTTG1 and HLA-DMA, CDCA8 and HLA-DQB2 were significantly correlated.

CONCLUSIONS

Examining the expression levels of four primary genes in cSCC aids in our understanding of the disease's pathophysiology. Additionally, the core genes were found to be highly related with immune regulatory genes, suggesting novel avenues for cSCC prevention and treatment.

Multimerin-1 and cancer: a review

Multimerin-1 (MMRN1) is a platelet protein with a role in haemostasis and coagulation. It is also present in endothelial cells (ECs) and the extracellular matrix (ECM), where it may be involved in cell adhesion, but its molecular functions and protein–protein interactions in these cellular locations have not been studied in detail yet. In recent years, MMRN1 has been identified as a differentially expressed gene (DEG) in various cancers and it has been proposed as a possible cancer biomarker. Some evidence suggest that MMRN1 expression is regulated by methylation, protein interactions, and non-coding RNAs (ncRNAs) in different cancers. This raises the questions if a functional role of MMRN1 is being targeted during cancer development, and if MMRN1’s differential expression pattern correlates with cancer progression. As a result, it is timely to review the current state of what is known about MMRN1 to help inform future research into MMRN1’s molecular mechanisms in cancer.

Identification of MAD2L1 as a Potential Biomarker in Hepatocellular Carcinoma via Comprehensive Bioinformatics Analysis

Background

Hepatocellular carcinoma (HCC) is widely acknowledged as a malignant tumor with rapid progression, high recurrence rate, and poor prognosis. At present, there is a paucity of reliable biomarkers at the clinical level to guide the management of HCC and improve patient outcomes. Our research is aimed at assessing the prognostic value of MAD2L1 in HCC.

Methods

Four datasets, GSE121248, GSE101685, GSE85598, and GSE62232, were selected from the GEO database to analyze differentially expressed genes (DEGs) between HCC and normal liver tissues. After functional analysis, we constructed a protein-protein interaction network (PPI) for DEGs and identified core genes in this network with high connectivity with other genes. We assessed the relationship between core genes and the pathogenesis and prognosis of HCC. Finally, we explored the gene regulatory signaling mechanisms involved in HCC pathogenesis.

Results

145 DEGs were screened from the intersection of the four GEO datasets. MAD2L1 was associated with most genes according to the PPI network and was selected as a candidate gene for further study. Survival analysis suggested that high MAD2L1 expression in HCC correlated with a worse prognosis. In addition, real-time quantitative PCR (RT-qPCR), western blot (WB), and immunohistochemistry (IHC) findings suggested that the expression of MAD2L1 was abnormally increased in HCC tissues and cells compared to paraneoplastic tissues and normal hepatocytes.

Conclusion

We found that high MAD2L1 expression in HCC was significantly associated with overall patient survival and clinical features. We also explored the potential biological properties of this gene.

Comprehensive Bioinformatics Analysis to Identify the Gene HMMR Associated With Lung Adenocarcinoma Prognosis and Its Mechanism of Action in Multiple Cancers

Background

Lung cancer is the third most frequently diagnosed cancer in the world, with lung adenocarcinoma (LUAD) as the most common pathological type. But studies on the predictive effect of a single gene on LUAD are limited. We aimed to discover new predictive markers for LUAD.

Methods

Differentially high-expressed genes at each stage were obtained from the TCGA and GTEx databases. The functions of these genes were investigated through GO enrichment and KEGG pathway analyses. Then, the key genes were selected by applying whole gene overall survival time. The expression of the key gene was studied in LUAD, and survival analysis was performed using Kaplan-Meier mapper, followed by univariate and multifactorial COX analysis. Finally, the gene expression and its prognostic significance in the pan-cancer were examined.

Results

A total of 10,106 DEGs were obtained from the two datasets. The top 266 differentially upregulated genes intersected with the top 1,497 overall survival-related genes, and 87 key genes were identified. High-expressed HMMR was associated with a poor prognosis of LUAD. Univariate and multifactorial Cox analysis showed that HMMR was an independent prognostic factor for LUAD patients. A high HMMR expression was strongly associated with the overall survival (OS) and disease-specific survival (DSS) in 11 cancer types and with poorer OS, DSS, and PFI in 10 cancer types.

Conclusion

HMMR may be an independent prognostic indicator and an important biomarker in diagnosing and predicting the survival of LUAD patients. Also, HMMR may be a key predictor of a variety of cancers.

NMU Is a Poor Prognostic Biomarker in Patients with Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is the most prevalent histologic type of lung cancer, associated with a high incidence rate and substantial mortality rate worldwide. Accumulating evidence shows that the aberrant expression of neuromedin U (NMU) contributes to the initiation and progression of cancer. Herein, we explored whether NMU could be adopted as a new diagnostic and therapeutic marker in LUAD. The UALCAN and GEPIA web resources were employed to assess data on the NMU expression in LUAD. The STRING web resource was used to develop the PPI (protein-protein interaction) network of NMU, whereas Cytoscape was applied for module analysis. The Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of NMU and the interacting proteins were examined using the WebGestalt tool. Survival analysis was performed with the Kaplan-Meier plotter tool. Results revealed that the NMU expression in LUAD was significantly higher than in the nonmalignant tissues. Moreover, higher NMU levels were dramatically related to shorter overall survival, first progression survival, and postprogression survival. The specific gene mutations G45V, R143T, and F152L of NMU occurred in LUAD samples and were associated with a worse prognosis in patients. KEGG and western blot analyses demonstrated an association of NMU with the cell cycle and the cAMP signaling cascade. Bioinformatic analysis and the in vitro experiments implicated NMU as a promising prognostic signature and treatment target for LUAD.

Prognostic and clinicopathological value of BUB1B expression in patients with lung adenocarcinoma: a meta-analysis

BACKGROUND

Abnormal BUB1B expression has been proven to be related to the poor prognosis of various tumors. This meta-analysis aimed to identify the prognostic role of BUB1B in patients with lung adenocarcinoma (LUAD).

RESEARCH DESIGN AND METHODS

Relevant studies from the PubMed, Embase, Web of Science, and Cochrane Library databases and two public databases that stored sequencing data were retrieved. The standardized mean difference (SMD) and 95% confidence intervals (CIs) for the association between the BUB1B expression level and clinical characteristics were calculated. Pooled hazard ratios (HRs) and 95% CIs were calculated to estimate the association between BUB1B expression and survival outcomes.

RESULTS

A total of 16 studies involving 2771 LUAD patients with BUB1B expression were included in this meta-analysis. Patients with older age showed low BUB1B expression. High BUB1B expression was associated with male sex, a smoking history, and an advanced TNM stage. High BUB1B expression was predictive of poor overall survival (OS) and progression-free survival (PFS). In addition, no publication bias was found.

CONCLUSIONS

This meta-analysis demonstrates that BUB1B is a significant biomarker for a poor prognosis and poor clinicopathological outcomes in patients with LUAD.

Integrated analysis of the molecular mechanisms in idiopathic pulmonary fibrosis

Rationale: Idiopathic pulmonary fibrosis (IPF) is one of the most aggressive forms of idiopathic interstitial pneumonia. Some miRNAs may be associated with IPF and may affect the occurrence and development of IPF in various pathways. Many miRNAs and genes that may be involved in the development of IPF have been discovered using chip and high throughput technologies. Methods: We analyzed one miRNA and four mRNA databases. We identified hub genes and pathways related to IPF using GO, KEGG enrichment analysis, gene set variation analysis (GSVA), PPI network construction, and hub gene analysis. A comprehensive analysis of differentially expressed miRNAs (DEMs), predicted miRNA target genes, and differentially expressed genes (DEGs) led to the creation of a miRNA-mRNA regulatory network in IPF. Results: We found 203 DEGs and 165 DEMs that were associated with IPF. The findings of enrichment analyses showed that these DEGs were mainly involved in antimicrobial humoral response, antimicrobial humoral immune response mediated by antimicrobial peptide, extracellular matrix organization, cell killing, and organ or tissue specific immune response. The VEGFA, CDH5, and WNT3A genes overlapped between hub genes and the miRNA-mRNA regulatory network. The miRNAs including miR-199b-5p, miR-140-5p, miR-199a-5p, miR-125A-5p, and miR-107 that we predicted would regulate the VEGFA, CDH5, and WNT3A genes, which were also associated with IPF or other fibrosis-related diseases. GSVA indicated that metabolic processes of UTP and IMP, immune response, regulation of Th2 cell cytokine production, and positive regulation of NK cell-mediated immunity are associated with the pathogenesis and treatment of IPF. These pathways also interact with VEGFA, CDH5, and WNT3A. Conclusion: These findings provide a new research direction for the diagnosis and treatment of IPF.