Idiopathic pulmonary fibrosis (IPF) is a chronic progressive lung disease with reduced quality of life and earlier mortality, but its pathogenesis and key genes are still unclear. In this investigation, bioinformatics was used to deeply analyze the pathogenesis of IPF and related key genes, so as to investigate the potential molecular pathogenesis of IPF and provide guidance for clinical treatment. Next generation sequencing dataset GSE213001 was obtained from Gene Expression Omnibus (GEO), the differentially expressed genes (DEGs) were identified between IPF and normal control group. The DEGs between IPF and normal control group were screened with the DESeq2 package of R language. The Gene Ontology (GO) and REACTOME pathway enrichment analyses of the DEGs were performed. Using the g:Profiler, the function and pathway enrichment analysis of DEGs were performed. Then, a protein protein interaction (PPI) network was constructed via the Integrated Interactions Database (IID) database. Cytoscape with Network Analyzer was used to identify the hub genes. miRNet and NetworkAnalyst database was used to construct the targeted microRNAs (miRNAs) and transcription factors (TFs) of the hub genes. Finally receiver operating characteristic (ROC) curve analysis was used to validates the hub genes. A total of 958 DEGs were screened out in this study, including 479 up regulated genes and 479 down regulated genes. Most of the DEGs were significantly enriched in response to stimulus, GPCR ligand binding, microtubule-based process and defective GALNT3 causes HFTC. In combination with the results of the PPI network, miRNA-hub gene regulatory network and TF-hub gene regulatory network, hub genes including LRRK2, BMI1, EBP, MNDA, KBTBD7, KRT15, OTX1, TEKT4, SPAG8 and EFHC2 were selected. Our findings will contribute to identification of potential biomarkers and novel strategies for the treatment of IPF, and provide a novel strategy for clinical therapy.