Rationale:The respiratory tract constitutes an elaborated line of defense based on a unique cellular ecosystem. Single-cell profiling methods enable the investigation of cell population distributions and transcriptional changes along the airways. Methods:We have explored cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. 77,969 cells were collected by bronchoscopy at 35 distinct locations, from the nose to the 12 th division of the airway tree. Results:The resulting atlas is composed of a high percentage of epithelial cells (89.1%), but also immune (6.2%) and stromal (4.7%) cells with peculiar cellular proportions in different sites of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type specific gene expression was stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuroendocrine (PNEC) and brush cells, which are likely derived from a common population of precursor cells. We also report a population of KRT13 positive cells with a high percentage of dividing cells which are reminiscent of "hillock" cells previously described in mouse.Conclusions: Robust characterization of this unprecedented large single-cell cohort establishes an important resource for future investigations. The precise description of the continuum existing from nasal epithelium to successive divisions of lung airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed.

The upper airway epithelium, which is mainly composed of multiciliated, goblet, club and basal cells, ensures proper mucociliary function and can regenerate in response to assaults. In chronic airway diseases, defective repair leads to tissue remodeling. Delineating key drivers of differentiation dynamics can help understand how normal or pathological regeneration occurs. Using single-cell transcriptomics and lineage inference, we have unraveled trajectories from basal to luminal cells, providing novel markers for specific populations. We report that: (1) a precursor subgroup of multiciliated cells, which we have entitled deuterosomal cells, is defined by specific markers, such as DEUP1, FOXN4, YPEL1, HES6 and CDC20B; (2) goblet cells can be precursors of multiciliated cells, thus explaining the presence of hybrid cells that co-express markers of goblet and multiciliated cells; and (3) a repertoire of molecules involved in the regeneration process, such as keratins or components of the Notch, Wnt or BMP/TGFβ pathways, can be identified. Confirmation of our results on fresh human and pig airway samples, and on mouse tracheal cells, extend and confirm our conclusions regarding the molecular and cellular choreography at work during mucociliary epithelial differentiation.

248 words) Rationale: The respiratory tract constitutes an elaborated line of defense based on a unique cellular ecosystem. Single-cell profiling methods enable the investigation of cell population distributions and transcriptional changes along the airways. Methods:We have explored cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. 77,969 cells were collected by bronchoscopy at 35 distinct locations, from the nose to the 12 th division of the airway tree. Results:The resulting atlas is composed of a high percentage of epithelial cells (89.1%), but also immune (6.2%) and stromal (4.7%) cells with peculiar cellular proportions in different sites of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type specific gene expression was stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuroendocrine (PNEC) and brush cells, which are likely derived from a common population of precursor cells. We also report a population of KRT13 positive cells with a high percentage of dividing cells which are reminiscent of "hillock" cells previously described in mouse. Conclusions:Robust characterization of this unprecedented large single-cell cohort establishes an important resource for future investigations. The precise description of the continuum existing from nasal epithelium to successive divisions of lung airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed. Results Building a molecular cell atlas of the airways in healthy volunteers Data collectionCells were analyzed by droplet-based single-cell RNA sequencing (scRNA-seq), after isolation from 4 distinct locations using 2 sampling methods: (i) nasal biopsies (3 samples) and (ii) nasal brushings (4 samples), (iii) tracheal biopsies (carina, 1 st division, 9 samples), (iv) intermediate bronchial biopsies (5-6 th divisions, 10 samples), (v) distal brushings (9-12 th divisions, 9 samples) in 10 healthy volunteers ( Figure 1A, 1B, Figure E1A, Table E1). Optimized handling and dissociation protocols allowed the profiling of 77,969 single cells which were collected at 35 distinct positions of the airways, resulting in the detection of an average of 1,892 expressed genes per cell with 7,070 UMI per cell ( Figure E2A).Following batch correction and graph-based clustering, cell types were assigned to each cluster using well-established sets of marker genes ( Figure 1C, Figure E3). We identified 14 epithelial cell types, including 12 for the surface epithelium and 2 for submucosal glands, which collectively represented 89.1% of total cells ( Figure 1C-1E, Table E2; See also our interactive web tool https://www.genomique.eu/cellbrowser/HCA/?ds=HCA_airway_epithelium).Strom...

Single cell RNA sequencing approaches are instrumental in studies of cell-to-cell variability. 5΄ selective transcriptome profiling approaches allow simultaneous definition of the transcription start size and have advantages over 3΄ selective approaches which just provide internal sequences close to the 3΄ end. The only currently existing 5΄ selective approach requires costly and labor intensive fragmentation and cell barcoding after cDNA amplification. We developed an optimized 5΄ selective workflow where all the cell indexing is done prior to fragmentation. With our protocol, cell indexing can be performed in the Fluidigm C1 microfluidic device, resulting in a significant reduction of cost and labor. We also designed optimized unique molecular identifiers that show less sequence bias and vulnerability towards sequencing errors resulting in an improved accuracy of molecule counting. We provide comprehensive experimental workflows for Illumina and Ion Proton sequencers that allow single cell sequencing in a cost range comparable to qPCR assays.

The MIR449 genomic locus encompasses several regulators of multiciliated cell (MCC) formation (multiciliogenesis). The miR-449 homologs miR-34b/c represent additional regulators of multiciliogenesis that are transcribed from another locus. Here, we characterized the expression of BTG4, LAYN, and HOATZ, located in the MIR34B/C locus using single-cell RNA-seq and super-resolution microscopy from human, mouse, or pig multiciliogenesis models. BTG4, LAYN, and HOATZ transcripts were expressed in both precursors and mature MCCs. The Layilin/LAYN protein was absent from primary cilia, but it was expressed in apical membrane regions or throughout motile cilia. LAYN silencing altered apical actin cap formation and multiciliogenesis. HOATZ protein was detected in primary cilia or throughout motile cilia. Altogether, our data suggest that the MIR34B/C locus may gather potential actors of multiciliogenesis.