The aim of the present study was to explore the function and interaction of differentially expressed genes (DEGs) in pulmonary embolism (PE). The gene expression profile GSE13535, was downloaded from the Gene Expression Omnibus database. The DEGs 2 and 18 h post‑PE initiation were identified using the affy package in R software. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of the DEGs were analyzed using Database for Annotation Visualization and Integrated Discovery (DAVID) online analytical tools. In addition, protein‑protein interaction (PPI) networks of the DEGs were constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins. The PPI network at 18 h was modularized using Clusterone, and a functional enrichment analysis of the DEGs in the top three modules was performed with DAVID. Overall, 80 and 346 DEGs were identified 2 and 18 h after PE initiation, respectively. The KEGG pathways, including chemokine signaling and toll‑like receptor signaling, were shown to be significantly enriched. The five highest degree nodes in the PPI networks at 2 or 18 h were screened. The module analysis of the PPI network at 18 h revealed 11 hub nodes. A Gene Ontology terms analysis demonstrated that the DEGs in the top three modules were associated with the inflammatory, defense and immune responses. The results of the present study suggest that the DEGs identified, including chemokine‑related genes TFPI2 and TNF, may be potential target genes for the treatment of PE. The chemokine signaling pathway, inflammatory response and immune response were explored, and it may be suggested that these pathways have important roles in PE.

Related Articles