Investigation of circular RNAs and related genes in pulmonary fibrosis based on bioinformatics analysis

Pulmonary fibrosis is a lethal inflammatory disease. In this study, we aimed to explore the potential‐related circular RNAs (circRNAs) and genes that are associated with pulmonary fibrosis. Pulmonary fibrosis rat models were constructed and the fibrosis deposition was detected using hematoxylin and eosin and Masson staining. The differentially expressed circRNAs were obtained through RNA sequencing. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were further performed to uncover the key function and pathways in pulmonary fibrosis. The interaction networks between circRNAs and their downstream micro RNAs (miRNAs) and genes were constructed by Cytoscape Software. The quantitative polymerase chain reaction was performed to validate the expression of 10 candidate circRNAs and five of them were performed ringwise sequencing in pulmonary fibrosis rats. We further selected five candidate circRNAs target miRNAs and messenger RNAs and validated by real‐time polymerase chain reaction. The pulmonary fibrosis models were successfully constructed according to the pathological examination. circRNAs were differentially expressed between the pulmonary fibrosis and normal pulmonary tissues. GO analysis verified that the differentially expressed circRNAs were significantly clustered in the cellular component, molecular function, and biological process. In the KEGG analysis, circRNAs were enriched in the following pathways: antigen processing and presentation, phagosome, PI3K‐AKt signaling pathway, HTLV‐I infection, and Herpes simplex infection. After validation in pulmonary fibrosis rat models, it was found that five of those circRNAs (chr9:113534327|113546234 [down], chr1:200648164|200672411 [down], chr5:150850432|150865550 [up], chr20:14319170|14326640 [down], and chr10:57634023|57634588 [down]) showed a relatively consistent trend with predictions. Validation of these circRNAs target miRNAs and genes showed that chr9:113534327|113546234, chr20:14319170|14326640, and chr10:57634023|57634588 were implicated in Notch1 activated transforming growth factor‐β (TGF‐β) signaling pathway. The study demonstrated that a series of circRNAs are differentially expressed in pulmonary fibrosis rats. These circRNAs, especially TGF‐β‐ and Notch1‐related circRNAs might play an important role in regulating pulmonary fibrogenesis.