NF-κB drives epithelial-mesenchymal mechanisms of lung fibrosis in a translational lung cell model

In the progression phase of idiopathic pulmonary fibrosis (IPF), the normal alveolar structure of the lung is lost and replaced by remodeled fibrotic tissue and by bronchiolized cystic airspaces. Although these are characteristic features of IPF, knowledge of specific interactions between these pathological processes is limited. Here, the interaction of lung epithelial and lung mesenchymal cells was investigated in a coculture model of human primary airway epithelial cells (EC) and lung fibroblasts (FB). Single-cell RNA sequencing revealed that the starting EC population was heterogenous and enriched for cells with a basal cell signature. Furthermore, fractions of the initial EC and FB populations adopted distinct pro-fibrotic cell differentiation states upon cocultivation, resembling specific cell populations that were previously identified in lungs of patients with IPF. Transcriptomic analysis revealed active NF-κB signaling early in the cocultured EC and FB, and the identified NF-κB expression signatures were found in “HAS1 High FB” and “PLIN2+ FB” populations from IPF patient lungs. Pharmacological blockade of NF-κB signaling attenuated specific phenotypic changes of EC and prevented FB-mediated interleukin-6, interleukin-8, and CXC chemokine ligand 6 cytokine secretion, as well as collagen α-1(I) chain and α–smooth muscle actin accumulation. Thus, we identified NF-κB as a potential mediator, linking epithelial pathobiology with fibrogenesis.

[1]  Malte D. Luecken,et al.  An integrated cell atlas of the human lung in health and disease , 2022, bioRxiv.

[2]  M. Efremova,et al.  CellPhoneDB: inferring cell–cell communication from combined expression of multi-subunit ligand–receptor complexes , 2020, Nature Protocols.

[3]  Mark D. Robinson,et al.  pipeComp, a general framework for the evaluation of computational pipelines, reveals performant single cell RNA-seq preprocessing tools , 2020, Genome Biology.

[4]  Manuel Stritt,et al.  Orbit Image Analysis: An open-source whole slide image analysis tool , 2019, bioRxiv.

[5]  Stefano Monti,et al.  hypeR: An R Package for Geneset Enrichment Workflows , 2019, bioRxiv.

[6]  R. Satija,et al.  Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression , 2019, Genome Biology.

[7]  A. Butte,et al.  Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage , 2018, Nature Immunology.

[8]  R. Welford,et al.  Novel high–throughput myofibroblast assays identify agonists with therapeutic potential in pulmonary fibrosis that act via EP2 and EP4 receptors , 2018, PloS one.

[9]  Luke Zappia,et al.  Clustering trees: a visualization for evaluating clusterings at multiple resolutions , 2018, bioRxiv.

[10]  Hannah A. Pliner,et al.  Reversed graph embedding resolves complex single-cell trajectories , 2017, Nature Methods.

[11]  Russell B. Fletcher,et al.  Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics , 2017, bioRxiv.

[12]  Alexander Lex,et al.  UpSetR: an R package for the visualization of intersecting sets and their properties , 2017, bioRxiv.

[13]  Mariella G. Filbin,et al.  Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma , 2016, Nature.

[14]  Andrew D. Rouillard,et al.  Enrichr: a comprehensive gene set enrichment analysis web server 2016 update , 2016, Nucleic Acids Res..

[15]  P. Linsley,et al.  MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data , 2015, Genome Biology.

[16]  Xiao-Jun Ma,et al.  Quantitative ultrasensitive bright-field RNA in situ hybridization with RNAscope. , 2014, Methods in molecular biology.

[17]  Edward Y. Chen,et al.  Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool , 2013, BMC Bioinformatics.

[18]  Davis J. McCarthy,et al.  Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation , 2012, Nucleic acids research.

[19]  J. Flanagan,et al.  RNAscope: a novel in situ RNA analysis platform for formalin-fixed, paraffin-embedded tissues. , 2012, The Journal of molecular diagnostics : JMD.

[20]  Helga Thorvaldsdóttir,et al.  Molecular signatures database (MSigDB) 3.0 , 2011, Bioinform..

[21]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[22]  Lokesh Kumar,et al.  Mfuzz: A software package for soft clustering of microarray data , 2007, Bioinformation.

[23]  Matthias E. Futschik,et al.  Noise-robust Soft Clustering of Gene Expression Time-course Data , 2005, J. Bioinform. Comput. Biol..

[24]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[25]  F. Speleman,et al.  Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes , 2002, Genome Biology.