Lung fibrosis, including idiopathic pulmonary fibrosis (IPF), represents a spectrum of progressive interstitial lung diseases characterized by disrupted epithelial repair, fibroblast activation, and excess extracellular matrix accumulation. A central feature of fibrotic progression is the loss of alveolar type 2 epithelial (AT2) cell identity and the emergence of aberrant transitional states that fail to support normal regeneration. Based on our re-analysis of multiple publicly available scRNA-seq datasets from the IPF Cell Atlas, S100A2 expression is tightly associated with basal-like reprogramming of AT2 cells and is a top upregulated gene at both mRNA and protein levels in IPF lungs. Furthermore, scRNA-seq data from human lung organoid models reveal that AT2 cells co-cultured with fibroblasts acquire a basal-like phenotype and express high levels of S100A2. Functionally, overexpression of S100A2 in human iPSC-derived lung alveolar epithelial type 2 (iAT2) cells leads to loss of AT2 cellular identity, increased generation of reactive oxygen species, and activation of RAGE signaling pathway. We demonstrate that pharmacologic inhibition of RAGE using Azeliragon preserves AT2 cell populations, associated with reduced oxidative stress in iAT2 cells, and significantly attenuates collagen deposition in bleomycin-induced lung fibrosis models. Collectively, our results indicate that S100A2 is a key driver of epithelial dysfunction in lung fibrosis, promoting loss of AT2 identity, aberrant basal fate acquisition, and persistent epithelial injury. Targeting the unrecognized S100A2-RAGE pathway may offer a new therapeutic strategy to restore epithelial homeostasis in lung fibrosis.

Repetitive injury is hypothesized to lead to progressive tissue fibrosis and end-stage organ failure. Whether tissue-resident mesenchymal cell populations retain epigenetic memory of prior injuries that contribute to this pathological process is unknown. Here we used a genetic lineage labeling approach to mark the lung mesenchyme prior to injury, then performed multi-modal analyses on isolated lung mesenchyme during the initiation, progression and resolution of the fibrotic response. Our results demonstrate the remarkable epigenetic and transcriptional plasticity of the lung mesenchyme during fibrotic activation and de-activation. Despite this plasticity, we also find that the lung mesenchyme exhibits an enhanced fibrotic program upon re-injury. We identify RUNX1 as a critical driver of both fibrotic activation and fibrotic memory. Comparison of fresh isolated and cultured lung mesenchyme demonstrates that RUNX1 is spontaneously activated in standard culture conditions, previously masking these roles of RUNX1. Targeted knockdown of RUNX1 dampens fibrotic mesenchymal cell activation immediately after cell isolation, but with reduced efficacy after only days of culture, confirming its functional importance to both early activation and long-term memory. Collectively, our findings implicate RUNX1 in the initiation and memory of fibrotic mesenchymal cell activation that together prime enhanced mesenchymal cell responses upon repeated injury.

Abstract Pulmonary fibrosis is characterized by pathological accumulation of scar tissue in the lung parenchyma. Many of the processes that are implicated in fibrosis, including increased extracellular matrix synthesis, also occur following pneumonectomy (PNX), but PNX instead results in regenerative compensatory growth of the lung. As fibroblasts are the major cell type responsible for extracellular matrix production, we hypothesized that comparing fibroblast responses to PNX and bleomycin (BLM) would unveil key differences in the role they play during regenerative versus fibrotic lung responses. RNA‐sequencing was performed on flow‐sorted fibroblasts freshly isolated from mouse lungs 14 days after BLM, PNX, or sham controls. RNA‐sequencing analysis revealed highly similar biological processes to be involved in fibroblast responses to both BLM and PNX, including TGF‐β1 and TNF‐α. Interestingly, we observed smaller changes in gene expression after PNX than BLM at Day 14, suggesting that the fibroblast response to PNX may be muted by expression of transcripts that moderate pro‐fibrotic pathways. Itpkc, encoding inositol triphosphate kinase C, was a gene uniquely up‐regulated by PNX and not BLM. ITPKC overexpression in lung fibroblasts antagonized the pro‐fibrotic effect of TGF‐β1. RNA‐sequencing analysis has identified considerable overlap in transcriptional changes between fibroblasts following PNX and those overexpressing ITPKC.

Alveolar type I (ATI) cells cover >95% of the lung's distal surface and facilitate gas exchange through their exceptionally thin shape. ATI cells in vivo are replenished by alveolar type II cell division and differentiation, but a detailed understanding of ATI biology has been hampered by the challenges in direct isolation of these cells due to their fragility and incomplete understanding of the signaling interactions that promote differentiation of ATII to ATI cells. Here, we explored the signals that maintain ATII versus promote ATI fates in three-dimensional (3-D) organoid cultures and developed a human alveolar type I differentiation medium (hATIDM) suitable for generating ATI cells from either mixed distal human lung cells or purified ATII cells. This media adds bone morphogenetic protein 4 (BMP4) and removes epidermal growth factor (EGF), Wnt agonist CHIR99021, and transforming growth factor-beta (TGF-β) inhibitor SB431542 from previously developed alveolar organoid culture media. We demonstrate that BMP4 promotes expression of the ATI marker gene AGER and HOPX, whereas CHIR99021 and SB431542 maintain expression of the ATII marker gene SFTPC. The human ATI spheroids generated with hATIDM express multiple molecular and morphological features reminiscent of human ATI cells. Our results demonstrate that signaling interactions among BMP, TGF-β, and Wnt signaling pathways in alveolar spheroids and distal lung organoids including IPF-organoids coordinate human ATII to ATI differentiation.NEW & NOTEWORTHY Alveolar type I (ATI) epithelial cells perform essential roles in maintaining lung function but have been challenging to study. We explored the signals that promote ATI fate in 3-D organoid cultures generated from either mixed distal human lung cells or purified alveolar type II (ATII) cells. This work fills an important void in our experimental repertoire for studying alveolar epithelial cells and identifies signals that promote human ATII to ATI cell differentiation.

Fibrosis in the lung is thought to be driven by epithelial cell dysfunction and aberrant cell-cell interactions. Unveiling the molecular mechanisms of cellular plasticity and cell-cell interactions is imperative to elucidating lung regenerative capacity and aberrant repair in pulmonary fibrosis. By mining publicly available RNA-Seq data sets, we identified loss of CCAAT enhancer-binding protein alpha (CEBPA) as a candidate contributor to idiopathic pulmonary fibrosis (IPF). We used conditional KO mice, scRNA-Seq, lung organoids, small-molecule inhibition, and potentially novel gene manipulation methods to investigate the role of CEBPA in lung fibrosis and repair. Long-term (6 months or more) of Cebpa loss in AT2 cells caused spontaneous fibrosis and increased susceptibility to bleomycin-induced fibrosis. Cebpa knockout (KO) in these mice significantly decreased AT2 cell numbers in the lung and reduced expression of surfactant homeostasis genes, while increasing inflammatory cell recruitment as well as upregulating S100a8/a9 in AT2 cells. In vivo treatment with an S100A8/A9 inhibitor alleviated experimental lung fibrosis. Restoring CEBPA expression in lung organoids ex vivo and during experimental lung fibrosis in vivo rescued CEBPA deficiency-mediated phenotypes. Our study establishes a direct mechanistic link between CEBPA repression, impaired AT2 cell identity, disrupted tissue homeostasis, and lung fibrosis.

Background:In recent years, with the development of transcriptome sequencing, the molecular characteristics of tumors are gradually revealed. Because of the complexity of tumor transcriptome, there is a need to look for the molecular signatures which can be used to evaluate the tissue origin and cell stemness of tumors in order to promote the diagnosis and treatment of tumors.Methods:Tumor tissue-specific gene sets (TTSGs) consisting of 200 genes were selected using RNA expression data of 9,875 patients from 33 tumor types. t-distributed Stochastic Neighbor Embedding (t-SNE) was used for dimensionality reduction and visualization of TTSGs in each tumor type. To evaluate oncogenic dedifferentiation and loss of cell stemness, Euclidean distance from each sample to a human embryo single-cell RNA-seq dataset (GSE36552) of TTSGs was calculated as TTSGs index indicating dissimilarity of tumors and embryo. TTSGs index was evaluated for prognosis in each tumor type. Two published signature indexes, the mRNA signature index (mRNAsi) and CIBERSORT, were compared to assess the correlation between the TTSGs index with cell stemness and immune microenvironment. Finally, the difference of prognosis, immune microenvironment and radiotherapy outcomes were compared between patients with high and low TTSGs index.Results:In this study, all 33 tumor types in The Cancer Genome Atlas (TCGA) were embedded into isolated clusters by t-SNE and confirmed by k-nearest neighbors (kNN) algorithm. Clusters of squamous-cell carcinoma were adjacent to each other revealing similar histologic origin. Basal-like breast cancer was separated from luminal and HER-2-amplified subtypes and closed to squamous-cell carcinoma. TTSGs index was related to overall survival outcomes in cancers derived from liver, thyroid, brain, cervical and kidney. There was a positive correlation between mRNAsi and TTSGs index in pan-kidney and pan-neuronal cancers. Furthermore, cell fractions of M2 macrophages and total leukocytes increased in the group with higher TTSGs index. Patients with higher TTSGs index had longer overall survival time and less radiation therapy resistance compared to patients with lower TTSGs index.Conclusions:The signature of TTSGs is related to tumor expression features that distinguish tumors of different histologic origin using t-SNE. The signature also relates to prognosis of certain kinds of tumors.

The lungs have a remarkable capacity to repair. However, repetitive injury can lead to progressive fibrosis and end-stage organ failure. Whether tissue-resident mesenchymal cell populations retain epigenetic memory of prior injuries that contribute to this pathological process is unknown. Here we used a genetic lineage labeling approach to mark the lung mesenchyme prior to injury, then performed multi-modal analyses on isolated lung mesenchyme during the initiation, progression and resolution of the fibrotic response. Our results demonstrate the remarkable epigenetic and transcriptional plasticity of the lung mesenchyme during fibrogenic activation and de-activation. Despite this plasticity, we also find that the lung mesenchyme retains specific epigenetic traits (memory) of prior activation, resulting in amplified induction of a fibrogenic program upon re-injury. We identify Runx1 as a critical driver of both fibrogenic activation and epigenetic memory. Comparison of fresh isolated and cultured lung mesenchyme demonstrates that Runx1 is spontaneously activated in standard culture conditions, previously masking these roles of Runx1. Genetic and pharmacological targeting of Runx1 dampens fibrogenic mesenchymal cell activation in cell and tissue models, confirming its functional importance. Finally, publicly available scRNAseq data reveal selective expression of Runx1 in the fibrogenic cell subpopulations that emerge in mouse and human fibrotic lung tissue. Collectively, our findings implicate Runx1 in both the initiation and memory of fibrogenic mesenchymal cell activation that together prime amplified mesenchymal cell responses upon repeated injury.

Abstract Background Idiopathic pulmonary fibrosis (IPF) is characterized by pathological accumulation of scar tissue in the lung parenchyma. Many of the processes that are implicated in fibrosis, such as matrix deposition also occur following pneumonectomy (PNX). However, fibrosis does not occur. Since fibroblasts are the major cell type responsible for extracellular matrix production, we hypothesized that comparing their responses to PNX and bleomycin (BLM) would unveil key differences in the role they play during regulated and fibrotic lung growth. Methods RNA-sequencing was performed on flow cytometry-sorted fibroblasts from mouse lungs 14 days after bleomycin challenge, PNX, or sham. Pathway and transcription factor binding motif enrichment analysis were performed to characterize RNA-sequencing data. Normal human lung fibroblasts (NHLFs) were used as a cell culture model to validate targets. Results RNA-sequencing analysis revealed similar biological processes to be involved in both responses, including signaling by transforming growth factor-β (TGF-β1) and tumor necrosis factor-α. Transcription factor binding motif enrichment predicted erythroid transformation specific (ETS) superfamily members to play a key role in the response to BLM, whereas nuclear factor-κB (NF-κB) and activator protein-1 (AP-1) were predicted to orchestrate much of the response to PNX. Itpkc, encoding inositol triphosphate kinase C, was a gene uniquely up-regulated by PNX and a likely AP-1 target. ITPKC overexpression in NHLFs antagonized the fibrotic effect of TGF-β1. RNA-sequencing analysis of primary NHLFs overexpressing ITPKC further supported its role in repressing collagen production and predicted it as an activator of canonical NF-κB signaling. ITPKC overexpression showed considerable overlap with the innate immune signaling seen following PNX. Conclusion Taken together, our RNA-sequencing analysis suggests that during post-PNX lung growth, AP-1 activates ITPKC to promote canonical NF-κB signaling to prevent fibrogenesis. Future studies to understand this lack of ITPKC induction during fibrotic injury may identify attractive therapeutic targets.

Idiopathic pulmonary fibrosis (IPF) is an aging-associated disease characterized by myofibroblast accumulation and progressive lung scarring. To identify transcriptional gene programs driving persistent lung fibrosis in aging, we performed RNA-Seq on lung fibroblasts isolated from young and aged mice during the early resolution phase after bleomycin injury. We discovered that, relative to injured young fibroblasts, injured aged fibroblasts exhibited a profibrotic state characterized by elevated expression of genes implicated in inflammation, matrix remodeling, and cell survival. We identified the proviral integration site for Moloney murine leukemia virus 1 (PIM1) and its target nuclear factor of activated T cells-1 (NFATc1) as putative drivers of the sustained profibrotic gene signatures in injured aged fibroblasts. PIM1 and NFATc1 transcripts were enriched in a pathogenic fibroblast population recently discovered in IPF lungs, and their protein expression was abundant in fibroblastic foci. Overexpression of PIM1 in normal human lung fibroblasts potentiated their fibrogenic activation, and this effect was attenuated by NFATc1 inhibition. Pharmacological inhibition of PIM1 attenuated IPF fibroblast activation and sensitized them to apoptotic stimuli. Interruption of PIM1 signaling in IPF lung explants ex vivo inhibited prosurvival gene expression and collagen secretion, suggesting that targeting this pathway may represent a therapeutic strategy to block IPF progression.

Matrix stiffness is a central regulator of fibroblast function. However, the transcriptional mechanisms linking matrix stiffness to changes in fibroblast phenotype are incompletely understood. Here, we evaluated the effect of matrix stiffness on genome-wide chromatin accessibility in freshly isolated lung fibroblasts using ATAC-seq. We found higher matrix stiffness profoundly increased global chromatin accessibility relative to lower matrix stiffness, and these alterations were in close genomic proximity to known profibrotic gene programs. Motif analysis of these regulated genomic loci identified ZNF416 as a putative mediator of fibroblast stiffness responses. Genome occupancy analysis using ChIP-seq confirmed that ZNF416 occupies a broad range of genes implicated in fibroblast activation and tissue fibrosis, with relatively little overlap in genomic occupancy with other mechanoresponsive and profibrotic transcriptional regulators. Using loss- and gain-of-function studies, we demonstrated that ZNF416 plays a critical role in fibroblast proliferation, extracellular matrix synthesis, and contractile function. Together, these observations identify ZNF416 as novel mechano-activated transcriptional regulator of fibroblast biology.

Fibroblast activation is transient in successful wound repair but persistent in fibrotic pathologies. Understanding fibroblast deactivation during successful wound healing may provide new approaches to therapeutically reverse fibroblast activation. To characterize the gene programs that accompany fibroblast activation and reversal during lung fibrosis resolution, we used RNA sequencing analysis of flow sorted Col1 alpha 1-GFP-positive and CD45-, CD31-, and CD326-negative cells isolated from the lungs of young mice exposed to bleomycin. We compared fibroblasts isolated from control mice with those isolated at Days 14 and 30 after bleomycin exposure, representing the peak of extracellular matrix deposition and an early stage of fibrosis resolution, respectively. Bleomycin exposure dramatically altered fibroblast gene programs at Day 14. Principal component and differential gene expression analyses demonstrated the predominant reversal of these trends at Day 30. Upstream regulator and pathway analyses of reversing "resolution" genes identified novel candidate antifibrotic genes and pathways. Two genes from these analyses that were decreased in expression at Day 14 and reversed at Day 30, Aldh2 and Nr3c1, were selected for further analysis. Enhancement of endogenous expression of either gene by CRISPR activation in cultured human idiopathic pulmonary fibrosis fibroblasts was sufficient to reduce profibrotic gene expression, fibronectin deposition, and collagen gel compaction, consistent with roles for these genes in fibroblast deactivation. This combination of RNA sequencing analysis of freshly sorted fibroblasts and hypothesis testing in cultured idiopathic pulmonary fibrosis fibroblasts offers a path toward identification of novel regulators of lung fibroblast deactivation, with potential relevance to understanding fibrosis resolution and its failure in human disease.

The supplemental data for manuscript titled 'TBK1 regulates YAP/TAZ and fibrogenic fibroblast activation'.

Uterine fibroids (UFs) are benign myometrial neoplasms. The mechanical environment activates signaling through the Hippo pathway effectors Yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding domain (TAZ) in other fibrotic disorders. Here, we assess the differences in YAP/TAZ responsiveness to signals in UF compared with myometrium (Myo). Matched samples of UF and Myo were collected. Atomic force microscopy (AFM) was used to determine in situ stiffness. Cells were plated sparsely on hydrogels or at confluence. Ten nanomolars of estradiol (E2) and 100 nM progesterone (P4) were used. Immunostaining for YAP/TAZ and extracellular matrix (ECM) proteins was performed. Cells were incubated with control or YAP1 (YAP)/WWTR1 (TAZ) small interfering RNA (siRNA). Real time qPCR was completed for connective tissue growth factor (CTGF). Cells were treated with verteporfin (a YAP inhibitor) or Y27632 (a ROCK inhibitor), and ECM gene expression was analyzed. Paired t test and Wilcoxon sign-rank test were used. AFM-measured tissue stiffness and YAP/TAZ nuclear localization in situ and in confluent cells were higher in UF compared with Myo (p < 0.05). Decreasing substrate stiffness reduced YAP/TAZ nuclear localization for both Myo and UF (p = 0.05). Stimulating cells with E2 or P4 increased YAP/TAZ nuclear localization, but only in Myo (p = 0.01). UFs had increased FN, COLI, and COLIII deposition. Following siRNA targeting, CTGF was found to be statistically decreased. Verteporfin treatment reduced cell survival and reduced FN deposition. Treatment with Y27632 demonstrated better cell tolerance and a reduction in ECM deposition. The mechanosensitive pathway may be linked to YAP/TAZ function and involved in transducing fibroid growth.

Idiopathic pulmonary fibrosis (IPF) results in scarring of the lungs by excessive extracellular matrix (ECM) production. Resident fibroblasts are the major cell type involved in ECM deposition. The biochemical pathways that facilitate pathological fibroblast activation leading to aberrant ECM deposition are not fully understood. Tank binding protein kinase-1 (TBK1) is a kinase that regulates multiple signaling pathways and was recently identified as a candidate regulator of fibroblast activation in a large-scale small-interfering RNA (siRNA) screen. To determine the effect of TBK1 on fibroblast activation, TBK1 was inhibited pharmacologically (MRT-68601) and genetically (siRNA) in normal and IPF human lung fibroblasts. Reducing the activity or expression of TBK1 led to reduction in α-smooth muscle actin stress fiber levels by 40–60% and deposition of ECM components collagen I and fibronectin by 50% in TGF-β-stimulated normal and IPF fibroblasts. YAP and TAZ are homologous mechanoregulatory profibrotic transcription cofactors known to regulate fibroblast activation. TBK1 knockdown or inhibition decreased the total and nuclear protein levels of YAP/TAZ. Additionally, low cell-cell contact and increased ECM substrate stiffness augmented the phosphorylation and activation of TBK1, consistent with cues that regulate YAP/TAZ. The action of TBK1 toward YAP/TAZ activation was independent of LATS1/2 and canonical downstream TBK1 signaling mediator IRF3 but dependent on proteasomal machinery of the cell. This study identifies TBK1 as a fibrogenic activator of human pulmonary fibroblasts, suggesting TBK1 may be a novel therapeutic target in pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive disease thought to result from impaired lung repair following injury and is strongly associated with aging. While vascular alterations have been associated with IPF previously, the contribution of lung vasculature during injury resolution and fibrosis is not well understood. To compare the role of endothelial cells (ECs) in resolving and non‐resolving models of lung fibrosis, we applied bleomycin intratracheally to young and aged mice. We found that injury in aged mice elicited capillary rarefaction, while injury in young mice resulted in increased capillary density. ECs from the lungs of injured aged mice relative to young mice demonstrated elevated pro‐fibrotic and reduced vascular homeostasis gene expression. Among the latter, Nos3 (encoding the enzyme endothelial nitric oxide synthase, eNOS) was transiently upregulated in lung ECs from young but not aged mice following injury. Young mice deficient in eNOS recapitulated the non‐resolving lung fibrosis observed in aged animals following injury, suggesting that eNOS directly participates in lung fibrosis resolution. Activation of the NO receptor soluble guanylate cyclase in human lung fibroblasts reduced TGFβ‐induced pro‐fibrotic gene and protein expression. Additionally, loss of eNOS in human lung ECs reduced the suppression of TGFβ‐induced lung fibroblast activation in 2D and 3D co‐cultures. Altogether, our results demonstrate that persistent lung fibrosis in aged mice is accompanied by capillary rarefaction, loss of EC identity, and impaired eNOS expression. Targeting vascular function may thus be critical to promote lung repair and fibrosis resolution in aging and IPF.

Reciprocal epithelial-mesenchymal interactions are pivotal in lung development, homeostasis, injury, and repair. Organoids have been used to investigate such interactions, but with a major focus on epithelial responses to mesenchyme and less attention to epithelial effects on mesenchyme. In the present study, we used nascent organoids composed of human and mouse lung epithelial and mesenchymal cells to demonstrate that healthy lung epithelium dramatically represses transcriptional, contractile, and matrix synthetic functions of lung fibroblasts. Repression of fibroblast activation requires signaling via the bone morphogenetic protein (BMP) pathway. BMP signaling is diminished after epithelial injury in vitro and in vivo, and exogenous BMP4 restores fibroblast repression in injured organoids. In contrast, inhibition of BMP signaling in healthy organoids is sufficient to derepress fibroblast matrix synthetic function. Our results reveal potent repression of fibroblast activation by healthy lung epithelium and a novel mechanism by which epithelial loss or injury is intrinsically coupled to mesenchymal activation via loss of repressive BMP signaling.