Whole transcriptome sequencing identifies key lncRNAs,circRNAs, and mRNAs for exploring the pathogenesis and therapeutic target of mouse pneumoconiosis

Pneumoconiosis is a kind of lung dysfunction caused by the inhalation of mineral dust. However, the potential molecular mechanism of pneumoconiosis have not been fully elucidated. In this study, the silica-treated pneumoconiosis mice model was constructed and the transcriptome sequencing data including lncRNA, circRNA, and mRNA were obtained. Firstly, differentially expressed lncRNA, circRNA, and mRNA (DElncRNA, DEcircRNA, DEGs) between control and pneumoconiosis/silicosis samples were screened, the target miRNAs (co-pre-miRNAs) were obtained by intersecting the miRNAs predicted by DElncRNA and DEcircRNA, respectively, and the target mRNAs (co-mRNA) were obtained by intersecting the mRNAs predicted by target miRNA and DEGs. Then, the lncRNA/circRNA-miRNA-mRNA networks were constructed by Cytoscape. Next, the key mRNAs were obtained by protein-protein interaction (PPI) analysis, and the key lncRNAs/circRNAs were selected by correlation analysis. Moreover, the expression of the key lncRNAs, circRNAs and mRNAs on chromosome were studied by the “circlize” package. Furthermore, the TFs-miRNA-mRNA network was constructed and the function of DEGs were explored by Ingenuity Pathway Analysis (IPA). To demonstrate the feasibility and value of the constructed ceRNA networks, we validated key genes and mmu-miR-682 pathway. Finally, We used the Drug-Gene Interaction database to predict potential drugs that could interfere with key genes,which may help to find promising treatment. There were 427 DElncRNAs, 107 DEcircRNAs and 1,597 DEGs between silicosis and control groups. Totals of 77 co-pre-miRNAs and 96 co-mRNA were screened, and the lncRNA/circRNA-miRNA-mRNA networks were constructed with 27 lncRNA/25 circRNAs, 74 miRNAs and 96 mRNAs. Then, 6 key mRNAs including Igf1, Klf4, Ptgs2, Epas1, Gnao1, and Il1a were obtained by PPI, and all of these key mRNAs and 10 key lncRNAs and 8 circRNAs were significantly different between the pneumoconiosis and normal groups, in which 10 lncRNAs and 9 circRNA that have not been previously studied in pneumoconiosis/silicosis can be used as new potential therapeutic targets. Moreover, the TFs-miRNA-mRNA network were constructed with 11 TFs, 1 key miRNA (mmu-miR-682) and 3 key mRNAs (Igf1, Epas1, Ptgs2). And the validation of key genes revealing by RNA-seq through experimental approaches shows the the predictive power of this study. Finally, IPA results indicated that 41 pathways were activated and 2 pathways were suppressed in pneumoconiosis/silicosis groups, and Pathogen Induced Cytokine Storm Signaling Pathway was the most significant pathway affected by pneumoconiosis/silicosis. In addition, 93 drugs were screened out by Drug-Gene Interaction database. Among them, Hydroxychloroquine was a kind of drug which associated with Il1a and Ptgs2, may be a promising treatment. This study constructed the lncRNA/circRNA-miRNA-mRNA and TFs-miRNA-mRNA networks, which could deepen the potential molecular regulatory mechanism of pneumoconiosis/silicosis.