A revised airway epithelial hierarchy includes CFTR-expressing ionocytes

The airways of the lung are the primary sites of disease in asthma and cystic fibrosis. Here we study the cellular composition and hierarchy of the mouse tracheal epithelium by single-cell RNA-sequencing (scRNA-seq) and in vivo lineage tracing. We identify a rare cell type, the Foxi1+ pulmonary ionocyte; functional variations in club cells based on their location; a distinct cell type in high turnover squamous epithelial structures that we term ‘hillocks’; and disease-relevant subsets of tuft and goblet cells. We developed ‘pulse-seq’, combining scRNA-seq and lineage tracing, to show that tuft, neuroendocrine and ionocyte cells are continually and directly replenished by basal progenitor cells. Ionocytes are the major source of transcripts of the cystic fibrosis transmembrane conductance regulator in both mouse (Cftr) and human (CFTR). Knockout of Foxi1 in mouse ionocytes causes loss of Cftr expression and disrupts airway fluid and mucus physiology, phenotypes that are characteristic of cystic fibrosis. By associating cell-type-specific expression programs with key disease genes, we establish a new cellular narrative for airways disease.Single-cell RNA sequencing analysis identifies cell types and lineages in airway epithelium, including the pulmonary ionocyte, a new cell type predominantly expressing the cystic fibrosis gene CFTR.

[1]  Timothy J. Mitchison,et al.  A chemical method for fast and sensitive detection of DNA synthesis in vivo , 2008, Proceedings of the National Academy of Sciences.

[2]  W. Wong,et al.  Annexin‐1‐deficient mice exhibit spontaneous airway hyperresponsiveness and exacerbated allergen‐specific antibody responses in a mouse model of asthma , 2011, Clinical and experimental allergy : journal of the British Society for Allergy and Clinical Immunology.

[3]  P. Tsao,et al.  Epithelial Notch signaling regulates lung alveolar morphogenesis and airway epithelial integrity , 2016, Proceedings of the National Academy of Sciences.

[4]  Martin Rosvall,et al.  Maps of random walks on complex networks reveal community structure , 2007, Proceedings of the National Academy of Sciences.

[5]  Freddy Radtke,et al.  Multiple roles of Notch signaling in the regulation of epidermal development. , 2008, Developmental cell.

[6]  R. Dixon,et al.  Requirement of a 5-lipoxygenase-activating protein for leukotriene synthesis , 1990, Nature.

[7]  Matthew D. Young,et al.  Gene ontology analysis for RNA-seq: accounting for selection bias , 2010, Genome Biology.

[8]  K. Gunsalus,et al.  Network modeling links breast cancer susceptibility and centrosome dysfunction. , 2007, Nature genetics.

[9]  Daniel J. Jackson,et al.  Cadherin-related family member 3, a childhood asthma susceptibility gene product, mediates rhinovirus C binding and replication , 2015, Proceedings of the National Academy of Sciences.

[10]  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.

[11]  J. Gardecki,et al.  An autoregulatory mechanism governing mucociliary transport is sensitive to mucus load. , 2014, American journal of respiratory cell and molecular biology.

[12]  Hui Liu,et al.  AnimalTFDB: a comprehensive animal transcription factor database , 2011, Nucleic Acids Res..

[13]  Berthold Göttgens,et al.  Clonal Dynamics Reveal Two Distinct Populations of Basal Cells in Slow-Turnover Airway Epithelium , 2015, Cell reports.

[14]  Sylvie Breton,et al.  The Forkhead Transcription Factor Foxi1 Is a Master Regulator of Vacuolar H+-ATPase Proton Pump Subunits in the Inner Ear, Kidney and Epididymis , 2009, PloS one.

[15]  E. Morrisey,et al.  Distinct Mesenchymal Lineages and Niches Promote Epithelial Self-Renewal and Myofibrogenesis in the Lung , 2017, Cell.

[16]  Fabian J Theis,et al.  The Human Cell Atlas , 2017, bioRxiv.

[17]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[18]  Pak Chung Sham,et al.  GWASdb v2: an update database for human genetic variants identified by genome-wide association studies , 2015, Nucleic Acids Res..

[19]  Gergely Katona,et al.  Gram-positive bacteria are held at a distance in the colon mucus by the lectin-like protein ZG16 , 2016, Proceedings of the National Academy of Sciences.

[20]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[21]  Yarden Katz,et al.  A single-cell survey of the small intestinal epithelium , 2017, Nature.

[22]  Sean C. Bendall,et al.  Data-Driven Phenotypic Dissection of AML Reveals Progenitor-like Cells that Correlate with Prognosis , 2015, Cell.

[23]  Shunsuke Kimura,et al.  Uptake through glycoprotein 2 of FimH+ bacteria by M cells initiates mucosal immune response , 2009, Nature.

[24]  M. Krasnow,et al.  Single-cell Wnt signaling niches maintain stemness of alveolar type 2 cells , 2018, Science.

[25]  Kirk M Druey,et al.  Suppression of immunoglobulin E–mediated allergic responses by regulator of G protein signaling 13 , 2008, Nature Immunology.

[26]  Richard C. Boucher,et al.  Regulation of Murine Airway Surface Liquid Volume by CFTR and Ca2+-activated Cl− Conductances , 2002, The Journal of general physiology.

[27]  E. Hoffman,et al.  Impaired mucus detachment disrupts mucociliary transport in a piglet model of cystic fibrosis , 2014, Science.

[28]  Theo Stijnen,et al.  Random effects meta‐analysis of event outcome in the framework of the generalized linear mixed model with applications in sparse data , 2010, Statistics in medicine.

[29]  Douglas W. Ball,et al.  Notch Signaling Induces Rapid Degradation of Achaete-Scute Homolog 1 , 2002, Molecular and Cellular Biology.

[30]  Vladimir Vinarsky,et al.  Ciliated cells of pseudostratified airway epithelium do not become mucous cells after ovalbumin challenge. , 2013, American journal of respiratory cell and molecular biology.

[31]  N. Neff,et al.  Reconstructing lineage hierarchies of the distal lung epithelium using single cell RNA-seq , 2014, Nature.

[32]  Hans-Peter Kriegel,et al.  A Density-Based Algorithm for Discovering Clusters in Large Spatial Databases with Noise , 1996, KDD.

[33]  Martin Donnelley,et al.  Airway disease phenotypes in animal models of cystic fibrosis , 2018, Respiratory Research.

[34]  James M. Wilson,et al.  Submucosal glands are the predominant site of CFTR expression in the human bronchus , 1992, Nature Genetics.

[35]  Colin N. Dewey,et al.  RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome , 2011, BMC Bioinformatics.

[36]  Fabian J. Theis,et al.  Diffusion maps for high-dimensional single-cell analysis of differentiation data , 2015, Bioinform..

[37]  Scott H. Randell,et al.  Basal cells as stem cells of the mouse trachea and human airway epithelium , 2009, Proceedings of the National Academy of Sciences.

[38]  Alexandro E. Trevino,et al.  Genome-scale transcriptional activation by an engineered CRISPR-Cas9 complex , 2014, Nature.

[39]  Fabian J Theis,et al.  Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells , 2015, Nature Biotechnology.

[40]  David K Meyerholz,et al.  Disease phenotype of a ferret CFTR-knockout model of cystic fibrosis. , 2010, The Journal of clinical investigation.

[41]  Andrew E. Jaffe,et al.  Bioinformatics Applications Note Gene Expression the Sva Package for Removing Batch Effects and Other Unwanted Variation in High-throughput Experiments , 2022 .

[42]  L. Ostrowski,et al.  Expression of CFTR from a ciliated cell-specific promoter is ineffective at correcting nasal potential difference in CF mice , 2007, Gene Therapy.

[43]  C. Nurse,et al.  Epithelial mitochondria‐rich cells and associated innervation in adult and developing zebrafish , 2006, The Journal of comparative neurology.

[44]  Aleksandra A. Kolodziejczyk,et al.  Accounting for technical noise in single-cell RNA-seq experiments , 2013, Nature Methods.

[45]  YanZiying,et al.  Optimization of Recombinant Adeno-Associated Virus-Mediated Expression for Large Transgenes, Using a Synthetic Promoter and Tandem Array Enhancers , 2015 .

[46]  Amy Y. Chen,et al.  Notch2 is required for inflammatory cytokine-driven goblet cell metaplasia in the lung. , 2015, Cell reports.

[47]  F. Finkelman,et al.  Distinct roles for IL-13 and IL-4 via IL-13 receptor α1 and the type II IL-4 receptor in asthma pathogenesis , 2008, Proceedings of the National Academy of Sciences.

[48]  Ajit Varki,et al.  Acidic pH increases airway surface liquid viscosity in cystic fibrosis. , 2016, The Journal of clinical investigation.

[49]  Ann B. Lee,et al.  Geometric diffusions as a tool for harmonic analysis and structure definition of data: multiscale methods. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[50]  Allon M. Klein,et al.  A single cell atlas of the tracheal epithelium reveals the CFTR-rich pulmonary ionocyte , 2018, Nature.

[51]  Linbo Liu,et al.  A functional anatomic defect of the cystic fibrosis airway. , 2014, American journal of respiratory and critical care medicine.

[52]  Kai Long,et al.  Cochlin produced by follicular dendritic cells promotes antibacterial innate immunity. , 2013, Immunity.

[53]  Andrew Gelman,et al.  Data Analysis Using Regression and Multilevel/Hierarchical Models , 2006 .

[54]  K. Kawakami,et al.  Mechanism of development of ionocytes rich in vacuolar-type H(+)-ATPase in the skin of zebrafish larvae. , 2009, Developmental biology.

[55]  Ian K Quigley,et al.  Specification of ion transport cells in the Xenopus larval skin , 2011, Development.

[56]  R. Shivdasani,et al.  Transcription factor foxq1 controls mucin gene expression and granule content in mouse stomach surface mucous cells. , 2008, Gastroenterology.

[57]  Ivana V. Yang,et al.  Muc5b is required for airway defence , 2013, Nature.

[58]  Evan Z. Macosko,et al.  Comprehensive Classification of Retinal Bipolar Neurons by Single-Cell Transcriptomics , 2016, Cell.

[59]  Kevin W Eliceiri,et al.  NIH Image to ImageJ: 25 years of image analysis , 2012, Nature Methods.

[60]  Fan Wang,et al.  The role of Scgb1a1+ Clara cells in the long-term maintenance and repair of lung airway, but not alveolar, epithelium. , 2009, Cell stem cell.

[61]  K. Chu,et al.  Method for Quantitative Study of Airway Functional Microanatomy Using Micro-Optical Coherence Tomography , 2013, PloS one.

[62]  Geoffrey E. Hinton,et al.  Visualizing Data using t-SNE , 2008 .

[63]  Cheng Li,et al.  Adjusting batch effects in microarray expression data using empirical Bayes methods. , 2007, Biostatistics.

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

[65]  Guillermo J Tearney,et al.  Development of an airway mucus defect in the cystic fibrosis rat. , 2018, JCI insight.

[66]  Laurens van der Maaten,et al.  Accelerating t-SNE using tree-based algorithms , 2014, J. Mach. Learn. Res..

[67]  Hongmei Mou,et al.  Dual SMAD Signaling Inhibition Enables Long-Term Expansion of Diverse Epithelial Basal Cells. , 2016, Cell stem cell.

[68]  M. Frick,et al.  A new role for P2X4 receptors as modulators of lung surfactant secretion , 2013, Front. Cell. Neurosci..

[69]  Borja Saez,et al.  Parent stem cells can serve as niches for their own daughter cells , 2015, Nature.

[70]  Scott M. Palmer,et al.  LungMAP: The Molecular Atlas of Lung Development Program , 2017, American journal of physiology. Lung cellular and molecular physiology.

[71]  A. Buja,et al.  Remarks on Parallel Analysis. , 1992, Multivariate behavioral research.

[72]  J. Engelhardt,et al.  Species-specific differences in mouse and human airway epithelial biology of recombinant adeno-associated virus transduction. , 2006, American journal of respiratory cell and molecular biology.

[73]  Linbo Liu,et al.  Ion Channels and Transporters in Lung Function and Disease Combination therapy with cystic fibrosis transmembrane conductance regulator modulators augment the airway functional microanatomy , 2016 .

[74]  Åsa K. Björklund,et al.  Full-length RNA-seq from single cells using Smart-seq2 , 2014, Nature Protocols.

[75]  David K. Meyerholz,et al.  Airway acidification initiates host defense abnormalities in cystic fibrosis mice , 2016, Science.