Fig. 4

RNA-seq analysis identifies unique transcriptional signatures in lungs of BLM-induced fibrotic mice. (A) Principal components analysis (PCA) of RNA-seq data from control and BLM-induced fibrotic mice. (B) Spearman correlation heatmap across all samples. (C) Volcano plot, visualizing the statistical difference of the differentially expression regions. (D-E) Pathway enrichment analysis of differentially expression regions between control and BLM-induced fibrotic mice. (F) Pathway activity inference of differentially expression regions between control and BLM-induced fibrotic mice. The observed activation of TNFa (G) and TGFb (H) is due to the fact that majority of its target genes with positive weights have positive t-values (1st quadrant), and the majority of the ones with negative weights have negative t-values (3d quadrant). (I) The dotplot illustrates the combined score and corresponding odds ratio for various signaling pathways, with a color gradient indicating the adjusted p-value. Yellow dots represent higher scores and odds ratios, while black dots indicate lower scores and ratios. (J) This GSEA plot presents a gene set enrichment analysis for the term “EPITHELIAL_MESENCHYMAL_TRANSITION,” focusing on the enrichment score and ranked metric across different ranks. The enrichment score shows a significant peak at the beginning, indicating high enrichment. (K) The barplot illustrates the activity levels of transcript factor in response to BLM, with each column representing a different transcript factor and each row indicating the activity level. Higher activity level is represented by red, while lower activity level is depicted in blue. (L) The TF-Gene co-regulatory network for Spi1, Spic, Jun and Mafb