BackgroundsPulmonary fibrosis (PF) is a pathological state presenting at the progressive stage of heterogeneous interstitial lung disease (ILD). The current understanding of the molecular mechanisms involved is incomplete. This clinical toxicology study focused on the pulmonary fibrosis induced by paraquat (PQ), a widely-used herbicide. Using proteo-transcriptome analysis, we identified differentially expressed proteins (DEPs) derived from the initial development of fibrosis to the dissolved stage and provided further functional analysis. MethodsWe established a mouse model of progressive lung fibrosis via intratracheal instillation of paraquat. To acquire a comprehensive and unbiased understanding of the onset of pulmonary fibrosis, we performed time-series proteomics profiling (iTRAQ) and RNA sequencing (RNA-Seq) on lung samples from paraquat-treated mice and saline control. The biological functions and pathways involved were evaluated through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway analysis. Correlation tests were conducted on comparable groups 7 days and 28 days post-exposure. Differentially expressed proteins and genes following the same trend on the protein and mRNA levels were selected for validation. The functions of the selected molecules were identified in vitro. The protein level was overexpressed by transfecting gene-containing plasmid or suppressed by transfecting specific siRNA in A549 cells. The levels of endothlial-mesenchymal transition (EMT) markers, including E-cadherin, vimentin, FN1, and α-SMA, were determined via western blot to evaluate the fibrotic process. ResultsWe quantified 1358 DEPs on day 7 and 426 DEPs on day 28 post exposure (Fold change >1.2; Q value < 0.05). The top 5 pathways – drug metabolism-cytochrome P450, metabolism of xenobiotics by cytochrome P450, complement and coagulation cascades, chemical carcinogenesis, protein digestion and absorption — were involved on both day 7 and day 28. Several pathways, including tight junction, focal adhesion, platelet activation, and ECM-receptor interaction, were more enriched on day 28 than on day 7. Integrative analysis of the proteome and transcriptome revealed a moderate correlation of quantitative protein abundance ratios with RNA abundance ratios (Spearman R = 0.3950 and 0.2477 on days 7 and 28, respectively), indicating that post-transcriptional regulation plays an important role in lung injury and repair. Western blot identified that the protein expressions of FN1, S100A4, and RBM3 were significantly upregulated while that of CYP1A1, FMO3, and PGDH were significantly downregulated on day 7. All proteins generally recovered to baseline on day 28. qPCR showed the mRNA levels of Fn1, S100a4, Rbm3, Cyp1a1, Fmo3, and Hpgd changed following the same trend as the levels of their respective proteins. Further, in vitro experiments showed that RBM3 was upregulated while PGDH was downregulated in an EMT model established in human lung epithelial A549 cells. RBM3 overexpression and PGDH knockout could both induce EMT in A549 cells. RBM3 knockout or PGDH overexpression had no reverse effect on EMT in A549 cells. ConclusionsOur proteo-transcriptomic study determined the proteins responsible for fibrogenesis and uncovers their dynamic regulation from lung injury to repair, providing new insights for the development of biomarkers for diagnosis and treatment of fibrotic diseases.