Comprehensive epigenomic profiling of human alveolar epithelial differentiation identifies key epigenetic states and transcription factor co-regulatory networks for maintenance of distal lung identity

Single-nucleus chromatin accessibility identifies a critical role for TWIST1 in idiopathic pulmonary fibrosis myofibroblast activity

BACKGROUND

In idiopathic pulmonary fibrosis (IPF), myofibroblasts are key effectors of fibrosis and architectural distortion by excessive deposition of extracellular matrix and their acquired contractile capacity. Single-cell RNA-sequencing (scRNA-seq) has precisely defined the IPF myofibroblast transcriptome, but identifying critical transcription factor activity by this approach is imprecise.

METHODS

We performed single-nucleus assay for transposase-accessible chromatin sequencing on explanted lungs from patients with IPF (n=3) and donor controls (n=2) and integrated this with a larger scRNA-seq dataset (10 IPF, eight controls) to identify differentially accessible chromatin regions and enriched transcription factor motifs within lung cell populations. We performed RNA-sequencing on pulmonary fibroblasts of bleomycin-injured Twist1-overexpressing COL1A2 Cre-ER mice to examine alterations in fibrosis-relevant pathways following Twist1 overexpression in collagen-producing cells.

RESULTS

TWIST1, and other E-box transcription factor motifs, were significantly enriched in open chromatin of IPF myofibroblasts compared to both IPF nonmyogenic (log2 fold change (FC) 8.909, adjusted p-value 1.82×10-35) and control fibroblasts (log2FC 8.975, adjusted p-value 3.72×10-28). TWIST1 expression was selectively upregulated in IPF myofibroblasts (log2FC 3.136, adjusted p-value 1.41×10- 24), with two regions of TWIST1 having significantly increased accessibility in IPF myofibroblasts. Overexpression of Twist1 in COL1A2-expressing fibroblasts of bleomycin-injured mice resulted in increased collagen synthesis and upregulation of genes with enriched chromatin accessibility in IPF myofibroblasts.

CONCLUSIONS

Our studies utilising human multiomic single-cell analyses combined with in vivo murine disease models confirm a critical regulatory function for TWIST1 in IPF myofibroblast activity in the fibrotic lung. Understanding the global process of opening TWIST1 and other E-box transcription factor motifs that govern myofibroblast differentiation may identify new therapeutic interventions for fibrotic pulmonary diseases.

Motif and conserved module analysis in DNA (promoters, enhancers) and RNA (lncRNA, mRNA) using AlModules

Nucleic acid motifs consist of conserved and variable nucleotide regions. For functional action, several motifs are combined to modules. The tool AIModules allows identification of such motifs including combinations of them and conservation in several nucleic acid stretches. AIModules recognizes conserved motifs and combinations of motifs (modules) allowing a number of interesting biological applications such as analysis of promoter and transcription factor binding sites (TFBS), identification of conserved modules shared between several gene families, e.g. promoter regions, but also analysis of shared and conserved other DNA motifs such as enhancers and silencers, in mRNA (motifs or regulatory elements e.g. for polyadenylation) and lncRNAs. The tool AIModules presented here is an integrated solution for motif analysis, offered as a Web service as well as downloadable software. Several nucleotide sequences are queried for TFBSs using predefined matrices from the JASPAR DB or by using one's own matrices for diverse types of DNA or RNA motif discovery. Furthermore, AIModules can find TFBSs common to two or more sequences. Demanding high or low conservation, AIModules outperforms other solutions in speed and finds more modules (specific combinations of TFBS) than alternative available software. The application also searches RNA motifs such as polyadenylation site or RNA-protein binding motifs as well as DNA motifs such as enhancers as well as user-specified motif combinations ( https://bioinfo-wuerz.de/aimodules/ ; alternative entry pages: https://aimodules.heinzelab.de or https://www.biozentrum.uni-wuerzburg.de/bioinfo/computing/aimodules ). The application is free and open source whether used online, on-site, or locally.

Impaired energy metabolism and altered functional activity of alveolar type II epithelial cells following exposure of rats to nitrogen mustard

Nitrogen mustard (NM) is a cytotoxic vesicant known to cause acute lung injury which progresses to fibrosis. Alveolar Type II cells are primarily responsible for surfactant production; they also play a key role in lung repair following injury. Herein, we assessed the effects of NM on Type II cell activity. Male Wistar rats were administered NM (0.125 mg/kg) or PBS control intratracheally. Type II cells, lung tissue and BAL were collected 3 d later. NM exposure resulted in double strand DNA breaks in Type II cells, as assessed by expression of γH2AX; this was associated with decreased expression of the DNA repair protein, PARP1. Expression of HO-1 was upregulated and nitrotyrosine residues were noted in Type II cells after NM exposure indicating oxidative stress. NM also caused alterations in Type II cell energy metabolism; thus, both glycolysis and oxidative phosphorylation were reduced; there was also a shift from a reliance on oxidative phosphorylation to glycolysis for ATP production. This was associated with increased expression of pro-apoptotic proteins activated caspase-3 and -9, and decreases in survival proteins, β-catenin, Nur77, HMGB1 and SOCS2. Intracellular signaling molecules important in Type II cell activity including PI3K, Akt2, phospho-p38 MAPK and phospho-ERK were reduced after NM exposure. This was correlated with dysregulation of surfactant protein production and impaired pulmonary functioning. These data demonstrate that Type II cells are targets of NM-induced DNA damage and oxidative stress. Impaired functioning of these cells may contribute to pulmonary toxicity caused by mustards.

FOXO1 Couples KGF and PI-3K/AKT Signaling to NKX2.1-Regulated Differentiation of Alveolar Epithelial Cells

NKX2.1 is a master regulator of lung morphogenesis and cell specification; however, interactions of NKX2.1 with various transcription factors to regulate cell-specific gene expression and cell fate in the distal lung remain incompletely understood. FOXO1 is a key regulator of stem/progenitor cell maintenance/differentiation in several tissues but its role in the regulation of lung alveolar epithelial progenitor homeostasis has not been evaluated. We identified a novel role for FOXO1 in alveolar epithelial cell (AEC) differentiation that results in the removal of NKX2.1 from surfactant gene promoters and the subsequent loss of surfactant expression in alveolar epithelial type I-like (AT1-like) cells. We found that the FOXO1 forkhead domain potentiates a loss of surfactant gene expression through an interaction with the NKX2.1 homeodomain, disrupting NKX2.1 binding to the SFTPC promoter. In addition, blocking PI-3K/AKT signaling reduces phosphorylated FOXO-1 (p-FOXO1), allowing accumulated nuclear FOXO1 to interact with NKX2.1 in differentiating AEC. Inhibiting AEC differentiation in vitro with keratinocyte growth factor (KGF) maintained an AT2 cell phenotype through increased PI3K/AKT-mediated FOXO1 phosphorylation, resulting in higher levels of surfactant expression. Together these results indicate that FOXO1 plays a central role in AEC differentiation by directly binding NKX2.1 and suggests an essential role for FOXO1 in mediating AEC homeostasis.

A single-cell regulatory map of postnatal lung alveologenesis in humans and mice

Ex-utero regulation of the lungs' responses to breathing air and continued alveolar development shape adult respiratory health. Applying single-cell transposome hypersensitive site sequencing (scTHS-seq) to over 80,000 cells, we assembled the first regulatory atlas of postnatal human and mouse lung alveolar development. We defined regulatory modules and elucidated new mechanistic insights directing alveolar septation, including alveolar type 1 and myofibroblast cell signaling and differentiation, and a unique human matrix fibroblast population. Incorporating GWAS, we mapped lung function causal variants to myofibroblasts and identified a pathogenic regulatory unit linked to lineage marker FGF18, demonstrating the utility of chromatin accessibility data to uncover disease mechanism targets. Our regulatory map and analysis model provide valuable new resources to investigate age-dependent and species-specific control of critical developmental processes. Furthermore, these resources complement existing atlas efforts to advance our understanding of lung health and disease across the human lifespan.