The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange. Three-dimensional in vitro human distal lung culture systems would strongly facilitate the investigation of pathologies such as interstitial lung disease, cancer and coronavirus disease 2019 (COVID-19) pneumonia caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here we describe the development of a long-term feeder-free, chemically defined culture system for distal lung progenitors as organoids derived from single adult human alveolar epithelial type II (AT2) or KRT5+ basal cells. AT2 organoids were able to differentiate into AT1 cells, and basal cell organoids developed lumens lined with differentiated club and ciliated cells. Single-cell analysis of KRT5+ cells in basal organoids revealed a distinct population of ITGA6+ITGB4+ mitotic cells, whose offspring further segregated into a TNFRSF12Ahi subfraction that comprised about ten per cent of KRT5+ basal cells. This subpopulation formed clusters within terminal bronchioles and exhibited enriched clonogenic organoid growth activity. We created distal lung organoids with apical-out polarity to present ACE2 on the exposed external surface, facilitating infection of AT2 and basal cultures with SARS-CoV-2 and identifying club cells as a target population. This long-term, feeder-free culture of human distal lung organoids, coupled with single-cell analysis, identifies functional heterogeneity among basal cells and establishes a facile in vitro organoid model of human distal lung infections, including COVID-19-associated pneumonia.
The distal lung performs essential respiratory functions that can be compromised by inflammatory, neoplastic or infectious disorders such as COVID-19 pneumonia. The study of these conditions would be facilitated by robust in vitro models based on human cells. The identities of the stem cells that regenerate distal lung epithelium in vivo over the human lifespan have been inferred from mouse studies, despite differences between these species1. In humans, basal stem cells span the entire airway tree, whereas in mouse, club cells2 and/or secretoglobin 1A1 (SCGB1A1)-expressing non-club cells3 renew the distal bronchioles during ageing. In the gas exchange region, mouse alveolar epithelial type II (AT2) cells generate AT1 and AT2 cells during homeostasis4,5, and additional progenitors are recruited in response to severe injury3,6,7,8,9. The presence and/or roles of facultative progenitors in human lung are unknown. Human AT2 cells can be differentiated into AT1 cells, but these cultures are short-lived and do not demonstrate long-term self-renewal or enable expansion for disease modelling4,10,11; furthermore, their requirements for feeder cells impede the definition of specific niche components12,13. We have established long-term organoid culture of distal human lung, including AT2 and basal stem cells, and used this system to validate functional progenitor cells and model SARS-CoV-2 infection.
Clonogenic alveolar and basal organoids
We empirically defined medium conditions to support the clonal expansion of distal human lung progenitors from 136 individuals in collagen/laminin extracellular matrix (ECM) (Fig. 1a, Supplementary Table 1). Together, EGF and the BMP antagonist NOGGIN supported growth of disaggregated distal lung cells or purified epithelial fractions thereof (Extended Data Fig. 1a–c). Single cells (Extended Data Fig. 1d–g) expanded into either SFTPC+HTII-280+ AT2 cystic organoids (Fig. 1b–e) or KRT5+ solid organoids expressing the basal cell marker KRT5 (Fig. 1b, f–h).
Single cell RNA sequencing (scRNA-seq) of distal lung organoids confirmed the presence of populations of SFTPC+ AT2 cells, KRT5+ basal cells and SCGB1A1+ club cells (Fig. 1i, j, Extended Data Fig. 2). Cells that co-expressed KRT5 and SCGB1A1 were found in both basal and club cell populations, suggesting a transitional state (Fig. 1j, k, Extended Data Fig. 2) that was supported by the results of SPADE14 and Monocle15 trajectory analysis (Extended Data Fig. 3, Supplementary Table 2).
Characterization of human AT2 cell organoids
We generated pure AT2 organoids from mixed cultures via uptake of lysosomal dyes into lamellar bodies (Extended Data Fig. 4a, b, Supplementary Fig. 1). EPCAM+LysoTracker+ AT2 cells expanded clonally for up to 180 days (Fig. 1l, m, Extended Data Fig. 4c, d). Chemically defined EGF/NOGGIN medium was sufficient for baseline clonal proliferation of AT2 organoids, which was attenuated by blocking global endogenous WNT biosynthesis (Extended Data Fig. 4e), consistent with the requirement for autocrine WNT signalling in mouse AT2 cells16. Growth was enhanced by adding fibroblast conditioned medium containing serum and WNT agonists (Extended Data Fig. 4f). Transmission electron microscopy revealed microvilli and lamellar bodies characteristic of mature AT2 cells (Fig. 1n, Extended Data Fig. 4g). AT2 organoids showed upregulation of AT1 markers when cultured on glass with serum (Fig. 1o, p).
Single-cell RNA-seq of mixed distal lung organoids or purified alveolar organoids revealed uniformly high expression of canonical AT2 cell markers in alveolar populations (Fig. 1i–k, q, Extended Data Fig. 2, Supplementary Data 1). AT2 cell subsets were not readily observed, and cell cycle-related mRNAs did not localize to a specific AT2 subpopulation (Fig. 1q, Supplementary Data 1, Supplementary Table 3).
Differentiation of human basal organoids
Basal cell organoids in mixed distal lung culture grew more rapidly than alveolar organoids and initially formed solid KRT5+ spheroids (Fig. 1a, b, f–h). However, by one month of culture, about 50% of organoids had developed single or occasionally multiple lumens (Fig. 2a, b, Extended Data Fig. 5a–c) with club (SCGB1A1+KRT5–) and ciliated (AcTUB+KRT5–) cells lining the interior surface (Fig. 2a, b). Basal cultures purified by density sedimentation exhibited serial clonal outgrowth, dependence on EGF and NOGGIN, and cavitation lined by luminal SCGB1A1+ club and AcTUB+ ciliated cells (Fig. 2c, d, Extended Data Fig. 5d–h, Supplementary Videos 1, 2). Similar differentiation occurred when organoids were transferred into 2D air–liquid interface (ALI) cultures (Fig. 2e).
Distinct subtypes of airway basal cells
Single-cell RNA-seq clustering of organoid KRT5+ basal cells from multiple individuals reproducibly identified two populations: Basal 1 and Basal 2 (Fig. 2f, g, Extended Data Fig. 6a, Supplementary Data 1). Basal 1 included an actively cycling subpopulation (Basal 1.2) that was enriched in proliferation markers (PCNA, CDK1) and gene cell enrichment analysis (GSEA) cell-cycle-related processes (Fig. 2f, g, Extended Data Fig. 6a–c, Supplementary Table 3). Basal 1, but not Basal 2, expressed canonical lung basal cell mRNAs such as TP6317, integrin-α6 (ITGA6) and integrin-β4 (ITGB4, which encodes a binding partner of ITGA618 and marks mouse lineage-negative epithelial progenitors (LNEPs))7 (Fig. 2h). Basal 2 was enriched in expression of genes related to vesicular transport, endoplasmic reticulum processes and squamous markers (Fig. 2g, Supplementary Table 3).
Characterization of TNFRSF12A+ progenitors
We examined the differentially expressed Basal 1 gene TNFRSF12A (also known as Fn14, TWEAKR), which encodes a membrane receptor (Fig. 2g, h, Supplementary Table 4) because of its potential utility for fluorescence-activated cell sorting (FACS) isolation and its homology to the intestinal stem cell marker TNFRSF1919. Unbiased pseudotime analysis revealed a continuous single-cell trajectory connecting KRT5+ Basal 1 cells to SCGB1A1+ club cells, with TNFRSF12A mRNA being strongly associated with the proliferation marker gene MKI67 (Fig. 2i, Extended Data Fig. 3e). When Basal 1 cells co-expressing EPCAM, ITGA6 and ITGB4 were divided into three fractions (expressing low, medium and high levels of TNFRSF12A mRNA), a proliferative gene module20 was significantly enriched in the highest (TNFRSF12Ahi) versus lowest fraction (TNFRSF12Alo) (Extended Data Fig. 6d–f). To determine whether this enrichment reflected intrinsic proliferative potential, we fractionated total distal lung organoids by FACS into EPCAM+ITGA6+ITGB4+ Basal 1 cells and then into TNFRSF12Ahi and TNFRSF12Aneg subsets (Fig. 2j). When cultured, the TNFRSF12Ahi subset showed 4–12 times greater clonogenic organoid-forming capacity than the TNFRSF12Aneg subset (Fig. 2k, l).
We examined lineage relationships between Basal 1 and Basal 2 by fractionating density-sedimented KRT5+ basal organoids by FACS into EPCAM+ITGA6+ITGB4+TNFRSF12Ahi (Basal 1) and EPCAM+ITGA6−ITGB4−TNFRSF12Aneg (Basal 2) populations (Extended Data Fig. 7a). Clonogenic organoid formation was strongly enriched in Basal 1 versus Basal 2 from three separate individuals (Extended Data Fig. 7b–d). The basal 2-enriched genes SPRR1B and TMSB4X (Fig. 2g, Supplementary Table 3) were transiently induced in Basal 1 cell organoids (Extended Data Fig. 7e, f), suggesting that Basal 2 cells might differentiate from Basal 1.
The NOTCH target gene HES1 was one of the most differentially expressed loci in Basal 1 organoids, in which gene networks included NOTCH1, NOTCH2 and JAG1 (Fig. 2g, Supplementary Table 3). Inhibition of NOTCH significantly increased basal organoid proliferation from TNFRSF12AhiEPCAM+ITGA6+ITGB4+ cells (Extended Data Fig. 7g, h), suggesting that NOTCH signalling restrains growth in these cells. Conversely, NOTCH agonism did not affect proliferation, but did induce expression of SCGB1A1, similar to upper airway cells21,22 (Extended Data Fig. 7i).
Characterization of TNFRSF12A+ cells in distal lung
Immunostaining of human distal lung tissue showed TNFRSF12A+ basal cells that were enriched intermittently at the tips or bases of bronchiolar furrows (Fig. 3a, Extended Data Fig. 8a, b); the latter are recognized as a goblet cell niche23. We found TNFRSF12A in diverse lung stromal and epithelial cells, but it clearly marked a minor population of KRT5+ and p63+ basal cells (Fig. 3a, Extended Data Fig. 8a–c). The TNFRSF12A+ subset of KRT5+ basal cells had a higher mitotic index than total KRT5+ cells in vivo (Fig. 3b, c), consistent with the results of organoid scRNA-seq (Fig. 2i). FACS analysis of human distal lung cells confirmed that TNFRSF12A was expressed in 10.9% of basal cells (Extended Data Fig. 8d, top).
Upon prospective culture directly from human lung without an organoid intermediate, FACS-isolated EPCAM+ITGA6+ITGB4+TNFRSF12Ahi cells (TNFRSF12Ahi Basal 1 cells; Extended Data Fig. 8d, bottom) showed a 15-fold increase in KRT5+ organoid formation compared with EPCAM+ITGA6+ITGB4+TNFRSF12Aneg cells (TNFRSF12neg Basal 1 cells; Fig. 3d, e). TNFRSF12Ahi basal organoids also differentiated into SCGB1A1+ club and AcTUB+ ciliated cells in prolonged culture (Fig. 3f) or when grown as 2D ALI monolayers (Fig. 3g).
SARS-CoV-2 and H1N1 infection of organoids
Influenza virus H1N1 avidly infected distal lung organoids, which also expressed influenza receptors (Extended Data Fig. 9a–d), similar to proximal airway organoids24,25. Infection of organoids with influenza H1N1 PR8 was inhibited by nucleoside analogues, consistent with previous studies26 (Extended Data Fig. 9e), and screening of diverse antiviral compound classes in 48-well format revealed differential activity (Extended Data Fig. 9f), suggesting that this system could be used for scalable therapeutics discovery.
In COVID-19 pneumonia, severe SARS-CoV-2 infection of the distal lung induces alveolar damage and respiratory failure27. Organoid scRNA-seq at time points before ciliated differentiation identified expression of the SARS-CoV-2 receptor ACE2 and processing protease TMPRSS2 mRNAs predominantly in club and AT2 cells (Extended Data Fig. 10a), consistent with expression of ACE2 in KRT5− differentiated interior luminal cells (Fig. 2a). To facilitate access of SARS-CoV-2 to ACE2-expressing luminal cells, we adapted an apical-out suspension culture polarization method28 to distal lung organoids (Extended Data Fig. 10b). Within 48 h in suspension, organoids reorganized into apical-out epithelial spheroids with microvilli, apical junctions, and motile cilia facing the organoid exterior. Differentiation of outwardly oriented ciliated cells accelerated over five days and progressed over weeks (Fig. 4a, Extended Data Fig. 10c–f, Supplementary Video 3). Apical-out basal organoids also showed an increase in outwardly facing club cells with apical secretory granules (Fig. 4a, Extended Data Fig. 10g), and apical-out AT2 organoids showed differentiation into AT1 cells (Extended Data Fig. 10h–j). Crucially, in apical-out organoids, ACE2 localized to apical cell membranes on the external organoid surface (Fig. 4b, Extended Data Fig. 10k).
SARS-CoV-2 infected apical-out mixed distal lung organoids, with induction of unspliced SARS-CoV-2 genomic RNA reaching levels similar to those of the abundantly expressed U3 small nucleolar RNA (Fig. 4c, left). In addition, infected organoids showed replication-specific SARS-CoV-2 spliced subgenomic RNA (sgRNA; Fig. 4c, right) and production of infectious virions with VeroE6 cell plaque formation (35 PFU ml−1 from organoid lysates and 65 PFU ml−1 from organoid supernatants). In basal organoids infected with SARS-CoV-2, double-stranded RNA (dsRNA) appeared by 48 h after infection (Fig. 4d) and SARS-CoV-2 nucleocapsid protein (NP) by 96 h (Fig. 4e). Approximately 10% of AT2 organoids displayed prominent SARS-CoV-2 NP expression in SFTPC+ cells; the remaining organoids were devoid of infection (Fig. 4f). Similarly, SARS-CoV-2 infected about 10% of basal organoids. In 2,621 total distal airway basal cells representing cultures from four individuals (Supplementary Table 5), SARS-CoV-2 infection was not detected in KRT5+ basal or AcTUB+ ciliated cells (odds ratio 0, P < 0.05), in contrast to infection of upper airway ciliated cells by SARS-CoV-2 in 2D ALI culture29,30. However, SARS-CoV-2 NP and dsRNA immunofluorescence were primarily present in SCGB1A1+ club cells (Fig. 4g, h) which were strongly associated with and accounted for 79% of NP- or dsRNA-positive cells (odds ratio 19.33, P < 0.0001); 21% of infected cells lacked SCGB1A1 (Fig. 4g, h, Supplementary Table 5). Overall, these studies indicate that AT2 cells were directly infected by SARS-CoV-2, and suggest that club cells are a distal lung target population.
We have applied long-term human distal lung organoid cultures to progenitor discovery and modelling of infectious disease. Our findings extend upon earlier short-term and feeder-dependent adult lung culture methods4,10,11 and present an alternative to techniques involving the differentiation of inducible pluripotent stem cells31,32,33,34. Organoids contained two related subtypes of KRT5+ human distal lung basal cells—Basal 1 and Basal 2. Notably, the TNFRSF12A-expressing Basal 1 fraction possessed enriched clonogenic progenitor activity, establishing functional precedent for a proliferation-enriched basal cell subtype. Although TNFRSF12A is by no means exclusively present in basal cells, within the basal layer it often localizes to a postulated niche in airway furrow bases and tips23, extending recent notions that the lung epithelium shows spatial specialization35. It is possible that TNFRSF12A or analogous markers could distinguish basal cell progenitor subsets in other tissues. Our organoids also enable facile exploration of SARS-CoV-2 distal lung infection, which is relevant to COVID-19-associated pneumonia27, and implicate SCGB1A1+ club cells as a target whose infection could compromise protective lung glycosaminoglycans and precipitate a vicious infection cycle. We did not observe infection of ciliated cells, in contrast to 2D ALI lung studies29,30; this could require alternative culture conditions. SCGB1A1-negative populations were also infected and are under further investigation; for example, bronchial transient secretory cells express ACE2 and TMPRSS236.
Overall, single-cell analysis of organoid cultures, as described here, may represent a general strategy for identifying and functionally validating candidate stem cells in slowly proliferating tissues. The culture of progenitors for all adult human distal lung epithelial lineages, including alveoli, should substantially enable the modelling of diseases such as neoplastic and interstitial lung conditions12 and allow tissue engineering and precision medicine applications. Finally, this organoid system should facilitate diverse investigations of pulmonary pathogens, including the SARS-CoV-2 distal lung infection that is associated with fulminant respiratory failure.
Additional experimental details are in the Supplementary Methods. No statistical methods were used to predetermine sample size. The experiments were not randomized and the investigators were not blinded to allocation during experiments and outcome assessment.
Human tissue procurement and processing
All experiments performed in this work were approved by the Stanford University School of Medicine Institutional Review Board and performed under protocol no. 28908. Standard informed consent for research was obtained in writing from all patients who contributed to this study before tissue procurement and all experiments followed relevant guidelines and regulations. Peripheral lung tissue within 1 cm of the visceral pleura was obtained from discarded surgical tissue from lobectomies. For patients with suspected lung cancer, cases with clinical T4 (American Joint Cancer Committee 6th edition) disease (for example, features such as bronchial invasion or parenchymal satellite nodule/metastases) were deferred. Normal tissue was removed from the lung margin most anatomically distal to palpably well-defined lesions, or from uninvolved lobes in the case of pneumonectomies. Samples with tumours containing ill-defined margins were deferred. Tissue was either processed fresh or stored at 4 °C overnight and processed the following morning.
Mixed distal lung organoid culture
To isolate distal airway cells, lung parenchyma from 1 cm from the visceral pleura was mechanically dissociated with Castro scissors, washed and incubated with 5 units per ml porcine elastase (Worthington), 100 Kunitz units per ml DNase I (Worthington), and Normocin (InvivoGen), and resuspended in two tissue volumes of lung organoid medium, comprising advanced DMEM/F12 (Invitrogen) supplemented with 10 mM nicotinamide, N-acetyl cysteine, 1× B27 supplement minus vitamin A, recombinant human NOGGIN (100 ng ml−1, R&D Systems), recombinant human EGF (50 ng ml−1, R&D Systems), and TGF-β inhibitor A83-01 (100 nM, Tocris). This lung organoid medium was used for all experiments except for those shown in Extended Data Fig. 4f. The tissue was then agitated for 1 h at 37 °C and the resultant cell suspension was filtered through 100–40-μm cell strainers and subjected to ammonium chloride red blood cell lysis. The cell pellet was then washed and resuspended in 10 volumes of reduced growth factor Basement Membrane Extract II (Trevigen). Cells in matrix were then plated in 24-well plates in 50-μl droplets, and warm medium was added after the droplets had solidified for 10 min at room temperature Medium was changed every 3–4 days and organoids were passaged every 3–4 weeks by dissociation with TrypLE. Passaging was based on ECM durability and integrity and estimated organoid confluency, judged by estimated organoid volume to volume of the ECM droplet. Distal lung organoids could be passaged for about 6 months with basal organoids initially exhibiting 6–7 doublings every 2 weeks. Alveolar organoids expanded more slowly with an initial rate of 3–4 doublings every 2 weeks but predominated over basal organoids after several months. Calculated from initial cell division rates, the upper limits of basal and alveolar expansion were 219 (524,288-fold) and 216 (65,536-fold), respectively. To rule out contamination by malignant cells, long-term cultures were systematically evaluated for the presence of dysplasia or carcinoma by a board-certified pathologist. In addition, five long-term organoid cultures (2–6 months) underwent targeted next generation sequencing to exclude pathogenic nucleotide variants (see below). Full details are provided in the Supplementary Methods.
Tandem MACS stromal depletion and EPCAM purification of distal lung cells
Distal lung was dissociated as above, and all incubation steps were carried out on ice. We incubated 107 cells with Fc Block (Biolegend 422301) and diluted them 1:100 in FACS buffer (2 mM EDTA and 0.2% fetal calf serum in 1× PBS, pH 7.4), for 10 min. The cells were then mixed with APC-conjugated anti-CD45 antibodies at 1 μg ml−1 in FACS buffer for 30 min, washed, and subjected to two rounds of depletion with magnetic beads according to the manufacturer’s protocol (Miltenyi: anti-human fibroblast 130-050-601, anti-CD31 130-091-935, anti-APC 130-090-855, LS column 130-042-401). Unlabelled cells were then centrifuged at 300g and labelled with a cocktail of 1 μg ml−1 of PerCP-Cy5.5 anti-EPCAM antibody and Zombie Aqua viability stain (Biolegend 423101) diluted 1:400 from stock concentration in FACS buffer.
Organoid cryopreservation and recovery
For cryopreservation and recovery, ECM droplets were dissociated by pipetting in three volumes of PBS with 5 mM EDTA and then incubated on ice for 1 h. Cells were pelleted at 300g for 5 min and resuspended in freezing medium (fetal calf serum (Gibco), 10% v/v DMSO), placed into cryovials and then into Mr. Frosty (Thermo Fisher) containers and stored at –80 °C overnight, followed by transfer to liquid nitrogen vapour phase for long-term storage. Organoids were recovered by quick thaw in a 37 °C water bath followed by washing in organoid medium and plating in ECM with organoid medium plus 10 μM ROCK inhibitor Y-27632 (Tocris).
Screening exogenous growth factors in organoid culture
Distal airway cells were isolated and plated as above with the following exceptions: advanced DMEM/F12 was used instead of organoid medium during elastase digestion of lung tissue, and cells were serially diluted and filtered through a 40-μm cell strainer and counted with a haemocytometer. One thousand viable epithelial cells (by Trypan blue exclusion, size, and morphology) per μl ECM were plated per 5 μl Matrigel droplet per well. Base medium consisted of organoid medium lacking A83-01, EGF, NOGGIN, WNT3A or RSPO1. EGF (final 50 ng ml−1, R&D), NOGGIN (final 100 ng ml−1, R&D), WNT3A (final 100 ng ml−1, R&D), RSPO1 (final 500 ng ml−1, Peprotech) or the PORCUPINE inhibitor C59 (final 1 μM, Biogems) were added singly or in combination to base medium. Images were obtained ten days after primary plating with an inverted light microscope at 5× magnification. Each condition was plated in quadruplicate and organoid formation was quantified using the analyse particle (threshold, 4902 pixels) plugin in ImageJ (see Supplementary Methods).
Single-cell RNA-seq of unfractionated organoid cultures
Lung organoid cultures from separate individuals were dissociated 4 weeks after primary plating and subjected to droplet-based scRNA-seq with the 10x Genomics Single Cell 3′ platform with a 5-nucleotide UMI, according to the manufacturer’s protocol. Cell capture, library preparation, and sequencing were performed as previously described37. For scRNA-seq analysis in Fig. 1k, a modified Kruskal–Wallis rank sum test was performed to determine the significance of differential marker gene expression for AT2 (SFTPC), basal (KRT5), and club (SCGB1A1) cells, with all P < 0.001. Principal component analysis, t-SNE, unsupervised graph-based clustering, statistical testing, and the pseudotime trajectory for all scRNA-seq analyses are described in Supplementary Methods and Supplementary Data 1.
Single-cell RNA-seq of purified AT2 organoid cultures
LysoTracker+ AT2 cells38 from unfractionated organoids were purified by FACS and cultured for two months with one passage. These were dissociated and subjected to droplet-based scRNA-seq with the 10x Genomics Chromium Single Cell 3′ platform v2 according to the manufacturer’s protocol. The library was sequenced using paired-end sequencing (26 bp read 1 and 98 bp read 2) with a single sample index (8 bp) on an Illumina NextSeq 500. Data preprocessing and principal component analysis were carried out with CellRanger v1.2. Subsequent analysis is described in Supplementary Methods and Supplementary Data 1.
Organoid cultures were fixed in ECM with 2.5% glutaraldehyde in 0.1 M cacodylate buffer (pH 7.4), dehydrated, embedded in epoxy resin and visualized with a JEOL (model JEM1400) transmission-electron microscope with a LaB6 emitter at 120 kV.
Histology and immunocytochemistry
Organoids were fixed with 2% paraformaldehyde at 4 °C overnight, paraffin embedded and sectioned (10–20 μm) as previously described37. Sections were deparaffinized and stained with H&E for histological analysis. Antibodies used for immunocytochemistry staining following standard staining protocols39 are listed in Supplementary Methods and images were acquired on a Leica-SP8 confocal microscope.
RNA fluorescent in situ hybridization
Whole-mount organoid confocal immunofluorescence microscopy
Intact, uninfected organoids were fixed in 2% paraformaldehyde in 100 mM phosphate buffer (pH 7.4) (4% paraformaldehyde for infected organoids) for 1 h at room temperature, washed with PBS with 100 mM glycine, permeabilized in 0.5% Triton X-100 in PBS for 1 h, then incubated in staining buffer (4% BSA, 0.05% Tween-20 in PBS pH 7.4, 10% goat/donkey serum) for an additional hour, followed by incubation with primary antibody for 24 h at room temperature in staining buffer. Whole mounts were then washed with PBS-T and incubated with fluorescent secondary antibodies, phalloidin and DAPI, for 4 h at room temperature in staining buffer. Following additional washes, whole mounts were submerged in mounting medium (Vectashield, Vector Laboratories) and mounted on chambered coverslips for imaging in four channels using Zeiss LSM 700 or 900 confocal microscopes. 3D rendering of confocal image stacks was performed using Volocity Image Analysis software (Quorum Technologies Inc., Guelph, Ontario). For Fig. 4h, which required five colours, cilia were distinguished by staining with two fluorescent secondary antibodies and merging the colocalized voxels into a pseudocoloured channel using Volocity software. Lectin staining (FITC-Sambuca Nigrin, Vector Labs FL-1301; Biotin-Maackia Amurensis, Vector Labs FL-1301) was carried out according to the manufacturer’s protocol after fixation of organoids with 0.1% paraformaldehyde in PBS for 1 h at room temperature followed by blocking with avidin/biotin (Vector Labs SP-2001). Biotin-Maackia Amurensis lectin was labelled with streptavidin–PE conjugate (Thermo Fisher SA10041) and after washing lectin staining was imaged in a Keyence BZ-X700.
Next generation sequencing of organoid cultures
Ten organoid cultures were sequenced using a commercial targeted resequencing assay with end-to-end coverage of 131 cancer genes and companion software (TOMA COMPASS Tumour Mutational Profiling System, Foster City, CA) to determine the presence of oncogenic mutations in long-term organoid cultures. Libraries were sequenced on an Illumina NextSeq 500. Nonsynonymous variants are listed in Supplementary Table 6. Variant call files are provided in Supplementary Data 2.
Density sedimentation of basal cells
Organoid cultures within 2–3 weeks of primary plating were dissociated with 1 U ml−1 neutral protease (Worthington, Cat LS02100) and 100 KU DNase I in lung organoid medium. Basal organoids were then collected by gravity sedimentation and the supernatant was either aspirated or collected for downstream use. Basal organoids were then further fractionated on a custom Ficoll-Paque gradient (4 vol Ficoll-Paque to 1 vol PBS) and centrifuged at 300g for 10 min at room temperature. Supernatant was aspirated and the organoid pellet was resuspended in 10 ml PBS in a 15-ml conical tube, collected by gravity sedimentation, and plated into ECM as described above.
FACS isolation and culture of AT2 cells
Organoids were dissociated with TrypLE followed by neutralization with 10% volume fetal calf serum, subjected to DNase at 100 KU ml−1, washed with lung organoid medium and then incubated with 100 cell pellet volumes of lung organoid medium with 10 nM LysoTracker Red DND-99 (Thermo Fisher L7528) at 37 °C for 30 min. Cells were then washed and resuspended in FACS buffer as described above, incubated with Fc block, and then incubated on ice with labelling cocktail consisting of 1 μg ml−1 of PerCP-Cy5.5 anti-EPCAM antibody and Zombie Aqua viability stain (Biolegend 423101) diluted 1:400 from stock concentration in FACS buffer. EPCAMhi and LysoTrackerhi cells were sorted into lung organoid medium with 10 μM Y-27632 (Tocris 1254) and cultured in ECM and lung organoid medium with Y-27632 for 24 h, followed by lung organoid medium without Y-27632. Pure AT2 organoid growth was enhanced by the addition of 1:1 vol:vol serum-containing L-cell conditioned medium (L-WRN CM) containing WNT3A, R-SPONDIN3 and NOGGIN and supplemented with recombinant EGF41. The full gating strategy is provided in Supplementary Fig. 1. All FACS antibodies were purchased from Biolegend. Qualitatively identical results could also be obtained with anti-HTII-280 (AT2 marker, Terrace Biotech) FACS purification in lieu of LysoTracker. AT2 organoid cells were transdifferentiated into AT1 cells by TrypLE dissociation from ECM and seeding onto chambered glass coverslips followed by culturing with advanced DMEM/F12 and 5% fetal calf serum42.
Colour mixing studies with lentivirally transduced GFP and mCherry
FACS EPCAM+ stroma-depleted organoids at D14 were infected with lentivirus at an estimated MOI of 0.9 as described previously43 with third-generation lentiviral vectors (PGK-GFP T2A Puro, SBI cat. no. CD550A-1; mCherry modified from pLentiCRISPRv1 (Addgene no. 49545) to incorporate an EF-1a-mCherry P2A Puro cassette, a gift from Paul Rack). Ninety-six hours after infection, organoids were treated with puromycin at a concentration of 600 ng ml−1 for 48 h to select for transduced cells. Two weeks after selection, GFP-expressing organoids or mCherry-expressing organoids were dissociated into single cells and mixed in a 1:1 ratio and scored as monochromatic or mixed after 28 days of each passage. The same approach was employed for purified AT2 and basal cultures after respective isolation strategies from an initial FACS-purified, EPCAM+, stroma-depleted organoid starter culture.
Flow cytometry analysis of resident basal cells from adult human lung
Adult human lung tissue was procured and dissociated as above but cells were labelled with Zombie Aqua live:dead stain as above, washed with FACS buffer, and then fixed in 2% PFA in PBS overnight at 4 °C. Cells were then stained using the whole-mount procedure as described above with the omission of PBS-glycine washing. Fixed and permeabilized cells were then incubated with 1:400 dilution of Alexa Fluor 647 conjugated mouse anti-human cytokeratin 5 antibody (Abcam) for 24 h at 4 °C in permeabilization buffer. Cells were then washed with FACS buffer and labelled with PE conjugated mouse anti-human TNFRSF12A antibody (clone ITEM-4, Biolegend) for 30 min on ice, followed by washing and analysis on a BD Aria Fusion instrument. The full gating strategy and qPCR validation of the ITEM-4 antibody are detailed in Supplementary Fig. 1.
FACS isolation of TNFRSF12Ahi and TNFRSF12Aneg basal cells
Single-cell suspensions from either fresh human distal lung or primary organoid culture at approximately 4 weeks of culture were dissociated as above, treated with Fc Block (BioLegend), and incubated in FACS buffer with Zombie Aqua 1:400, 1 μg ml−1 PerCP-Cy5.5 anti-human EpCAM (CD326), 1 μg ml−1 APC anti-human ITGA6 (CD49f), 2 μg ml−1 FITC anti-human ITGB4 (CD104), and 1 μg ml−1 PE anti-human TNFRSF12A (CD266). Thirty minutes after labelling, the cells were washed twice with FACS buffer and sorted for EPCAMhi, ITGA6/ITGB4hi, TNFRSF12Ahi and TNFRSF12Aneg. The full gating strategy is provided in Supplementary Fig. 1. More than 5,000 cells were sorted into Eppendorf tubes with lung organoid medium and 10 μM ROCK inhibitor Y-27632. All FACS antibodies were purchased from Biolegend.
Culture of TNFRSF12Ahi and TNFRSF12Aneg basal cells
Cells were seeded in ECM and submerged in lung organoid medium with 10 μM ROCK inhibitor Y-27632. The seeding density for cells FACS-isolated from organoid culture was 1,000 cells per well at a density of 100 cells per μl ECM. The seeding density for cells FACS-isolated from fresh human distal lung was 3,000 cells per well at a density of 300 cells per μl ECM. After 24 h, the medium was changed to remove ROCK inhibitor, thereafter it was changed every 72 h. Organoid formation was manually quantified 14 days after plating by two independent observers.
H1N1 organoid influenza assay
Unfractionated cultures containing AT2, basal, and club cell types at 2–3 weeks were infected in triplicate with H1N1 strain PR8 modified to express GFP upon viral replication44 after 24 h of pretreatment with antiviral compounds. ECM was dispersed by addition of 5 mM EDTA in PBS, followed by washing and inoculation with PR8-H1N1–GFP reporter virus at an estimated MOI of 1 in medium containing either vehicle or antivirals. After 12 h (one influenza infection cycle), intact organoid GFP expression was either visualized by fluorescence microscopy with a Keyence BZ-X700 automated microscope, or dissociated to single cells and fixed with 0.1% PFA in PBS followed by FACS quantification of GFP+ cells (gating strategy is provided in Supplementary Fig. 1). Antiviral dose response curves were generated using four-parameter nonlinear regression curve fitting with GraphPad Prism 7 (GraphPad Software, San Diego, CA). H1N1 tropism was assessed in a manner similar to above with the exception that Ficoll-sedimented basal cell fraction versus non-basal fractions were dissociated into single cells, counted, and infected with an estimated MOI of 1 in organoid medium for 1 h at 37 °C, followed by washing and reseeding into ECM, culturing for 16 h, and then being subjected to dissociation and FACS as above.
Suspension culture to generate apical-out polarity in lung organoids
Lung organoids grown embedded in 50 μl ECM droplets were transferred to suspension culture as described28 with modifications. In brief, ECM-embedded organoids were dislodged gently by pipetting using sterile LoBind tips (Eppendorf 22493008) and placed in 15 ml LoBind conical tubes (Eppendorf 30122216) containing ice-cold 5 mM EDTA in PBS. Five millilitres of EDTA solution was used per ECM droplet (3 ECM droplets per 15-ml conical tube) rotating for 1 h at 4 °C on a rotating platform. Organoids were centrifuged at 200g for 3 min at 4 °C and the supernatant was removed. The pellet was resuspended in growth medium in ultra-low attachment six-well tissue culture plates (Corning Costar 3471). Suspended organoids were incubated at 37 °C with 5% CO2 for different times (range 0–30 days) to characterize apical-out polarity, ciliogenesis and differentiation, and to prepare apical-out organoids for infection experiments with SARS-CoV-2.
SARS-CoV2 infection of human distal lung organoids
VeroE6 cells were obtained from ATCC as mycoplasma-free stocks and maintained in supplemented DMEM with 10% FBS. SARS-CoV-2 (USA-WA1/2020) was passaged in VeroE6 cells in DMEM with 2% FBS. Titres were determined by plaque assay on VeroE6 cells using Avicel (FMC Biopolymer) and crystal violet (Sigma), the viral genome sequence was verified, and all infections were performed with passage 3 virus. Organoids were counted and passaged into suspension medium for 6–8 days and then resuspended in virus medium or an equal volume of mock medium at an MOI of 1 relative to total organoid cells in the sample, and then incubated at 37 °C under 5% CO2 for 2 h. Organoids were then plated in suspension in lung organoid medium (apical-out organoids). At the indicated time points, apical-out organoids were washed with lung organoid medium and PBS and resuspended in TRIzol LS (Thermo Fisher), freshly-made 4% PFA in PBS, or 250 μl lung organoid medium. Cells resuspended in lung organoid medium were lysed by freezing at −80 °C. Culture supernatants were preserved in TRIzol LS or added directly to plaque assay monolayers. All SARS-CoV-2 work was performed in a class II biosafety cabinet under BSL3 conditions at Stanford University.
qPCR analysis of SARS-CoV-2 RNA
RNA from SARS-CoV-2-infected organoids was extracted by adding 750 μl TRIzol (Thermo Fisher Scientific), incubating at 55 °C for 5 min and then adding 150 μl chloroform. After mixing each sample by vortexing for 7 s, the samples were incubated at 25 °C for 5 min and then centrifuged at 12,000 rpm for 15 min at 4 °C. The aqueous layer was carefully removed from each sample, mixed with two volumes of 100% ethanol and purified using an RNA Clean & Concentrator-25 kit (Zymo Research) as per the manufacturer’s instructions. All RNA samples were treated with DNase (Turbo DNA-free kit, Thermo Fisher Scientific). The Brilliant II SYBR Green QRT–PCR 1-Step Master Mix (VWR) was used to convert RNA into cDNA and to amplify specific RNA regions on the CFX96 Touch real-time PCR detection system (Bio-Rad). The reverse transcription reaction was performed for 30 min at 50 °C, 10 min at 95 °C, followed by two-step qPCR with 95 °C for 10 s and 55 °C for 30 s, for a total of 40 cycles. Two primer sets were used to amplify either non-spliced SARS-CoV-2 genomic RNA (gRNA) spanning nucleotide positions 14221–14306, or spliced SARS-CoV-2 sgRNA30. Primer sequences are in Supplementary Table 7.
TNFRSF12A immunostaining of intact distal lung
Optimal staining of human distal lung tissue was achieved from specimens fixed within 30 min of primary surgical resections in 4% paraformaldehyde in PBS. Specimens were incubated in fixative overnight at 4 °C, transferred to 30% sucrose, and embedded into OCT. Frozen sections were cut at 10 μm, subjected to citrate-based antigen retrieval (Vector Labs) at 70 °C for 30 min, and then blocked for 1 h with 10% goat serum in IF wash buffer as described above. Mouse anti-TNFRSF12A (clone ITEM-4, Biolegend) was used for Fig. 3a and polyclonal rabbit anti-TNFRSF12A (ThermoFisher PA5-20275) was used for Fig. 3b, f and Extended Data Fig. 845.
Live-imaging and confocal microscopy of immobilized apical-out lung organoids
Live organoids were held between two coverslips in a viewing chamber (Lab-Tek II two-chambered coverglass) and filmed using a Nikon TE2000E microscope using differential interference contrast (DIC) microscopy with a 63× objective. Samples were kept at 37 °C with 5% CO2 during imaging. Digital videos were collected by a Hamamatsu high-resolution ORCA-285 digital camera and rendered using OpenLab 5.5.2 software (Improvision). After recording, samples were fixed and stained without removal from the chambers and transferred to the confocal microscope for immunofluorescence microscopy.
Statistics and reproducibility
Unless stated otherwise, all data are representative of at least two independent experiments with each independent experiment carried out using an organoid culture derived from one individual. Box plot bounds span first through third quartiles, horizontal lines represent median values, and whiskers represent data range minima or maxima or, in the case of outliers, 1.5 times the interquartile range with outliers represented by data points. t-tests were two-tailed and P values are denoted as *P < 0.05, **P ≤ 0.01, and ***P ≤ 0.001. Full details are provided in Supplementary Methods.
Further information on research design is available in the Nature Research Reporting Summary linked to this paper.
Scripts to perform analyses of scRNA-seq data are provided with this paper. Custom code is available on GitHub (https://github.com/ameen-salahudeen/lung_organoid).
Hogan, B. & Tata, P. R. Cellular organization and biology of the respiratory system. Nat. Cell Biol. https://doi.org/10.1038/s41556-019-0357-7 (2019).
Rawlins, E. L. et al. The role of Scgb1a1+ Clara cells in the long-term maintenance and repair of lung airway, but not alveolar, epithelium. Cell Stem Cell 4, 525–534 (2009).
Kathiriya, J. J., Brumwell, A. N., Jackson, J. R., Tang, X. & Chapman, H. A. Distinct airway epithelial stem cells hide among club cells but mobilize to promote alveolar regeneration. Cell Stem Cell 26, 346–358.e4 (2020).
Barkauskas, C. E. et al. Type 2 alveolar cells are stem cells in adult lung. J. Clin. Invest. 123, 3025–3036 (2013).
Desai, T. J., Brownfield, D. G. & Krasnow, M. A. Alveolar progenitor and stem cells in lung development, renewal and cancer. Nature 507, 190–194 (2014).
Liu, Q. et al. Lung regeneration by multipotent stem cells residing at the bronchioalveolar-duct junction. Nat. Genet. 51, 728–738 (2019).
Vaughan, A. E. et al. Lineage-negative progenitors mobilize to regenerate lung epithelium after major injury. Nature 517, 621–625 (2015).
Zuo, W. et al. p63+Krt5+ distal airway stem cells are essential for lung regeneration. Nature 517, 616–620 (2015).
Juul, N. H., Stockman, C. A. & Desai, T. J. Niche cells and signals that regulate lung alveolar stem cells in vivo. Cold Spring Harb. Perspect. Biol. https://doi.org/10.1101/cshperspect.a035717 (2020).
Sucre, J. M. S. et al. Successful establishment of primary type II alveolar epithelium with 3D organotypic coculture. Am. J. Respir. Cell Mol. Biol. 59, 158–166 (2018).
Zacharias, W. J. et al. Regeneration of the lung alveolus by an evolutionarily conserved epithelial progenitor. Nature 555, 251–255 (2018).
Nikolić, M. Z. & Rawlins, E. L. Lung organoids and their use to study cell-cell interaction. Curr. Pathobiol. Rep. 5, 223–231 (2017).
Evans, K. V. & Lee, J. H. Alveolar wars: the rise of in vitro models to understand human lung alveolar maintenance, regeneration, and disease. Stem Cells Transl. Med. 9, 867–881 (2020).
Anchang, B. et al. Visualization and cellular hierarchy inference of single-cell data using SPADE. Nat. Protoc. 11, 1264–1279 (2016).
Cao, J. et al. The single-cell transcriptional landscape of mammalian organogenesis. Nature 566, 496–502 (2019).
Nabhan, A. N., Brownfield, D. G., Harbury, P. B., Krasnow, M. A. & Desai, T. J. Single-cell Wnt signaling niches maintain stemness of alveolar type 2 cells. Science 359, 1118–1123 (2018).
Rock, J. R. et al. Basal cells as stem cells of the mouse trachea and human airway epithelium. Proc. Natl Acad. Sci. USA 106, 12771–12775 (2009).
Kajiji, S., Tamura, R. N. & Quaranta, V. A novel integrin (alpha E beta 4) from human epithelial cells suggests a fourth family of integrin adhesion receptors. EMBO J. 8, 673–680 (1989).
Stange, D. E. et al. Differentiated Troy+ chief cells act as reserve stem cells to generate all lineages of the stomach epithelium. Cell 155, 357–368 (2013).
Whitfield, M. L., George, L. K., Grant, G. D. & Perou, C. M. Common markers of proliferation. Nat. Rev. Cancer 6, 99–106 (2006).
Rock, J. R. et al. Notch-dependent differentiation of adult airway basal stem cells. Cell Stem Cell 8, 639–648 (2011).
Pardo-Saganta, A. et al. Injury induces direct lineage segregation of functionally distinct airway basal stem/progenitor cell subpopulations. Cell Stem Cell 16, 184–197 (2015).
Quinton, P. M. Both ways at once: keeping small airways clean. Physiology (Bethesda) 32, 380–390 (2017).
Zhou, J. et al. Differentiated human airway organoids to assess infectivity of emerging influenza virus. Proc. Natl Acad. Sci. USA 115, 6822–6827 (2018).
Imai, M. & Kawaoka, Y. The role of receptor binding specificity in interspecies transmission of influenza viruses. Curr. Opin. Virol. 2, 160–167 (2012).
Kumaki, Y., Day, C. W., Smee, D. F., Morrey, J. D. & Barnard, D. L. In vitro and in vivo efficacy of fluorodeoxycytidine analogs against highly pathogenic avian influenza H5N1, seasonal, and pandemic H1N1 virus infections. Antiviral Res. 92, 329–340 (2011).
Zhu, N. et al. A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 382, 727–733 (2020).
Co, J. Y. et al. Controlling epithelial polarity: a human enteroid model for host-pathogen interactions. Cell Rep. 26, 2509–2520.e4 (2019).
Lamers, M. M. et al. SARS-CoV-2 productively infects human gut enterocytes. Science 369, 50–54 (2020).
Hou, Y. J. et al. SARS-CoV-2 reverse genetics reveals a variable infection gradient in the respiratory tract. Cell 182, 429–446.e14 (2020).
Dye, B. R. et al. In vitro generation of human pluripotent stem cell derived lung organoids. eLife 4, e05098 (2015).
Chen, Y. W. et al. A three-dimensional model of human lung development and disease from pluripotent stem cells. Nat. Cell Biol. 19, 542–549 (2017).
Jacob, A. et al. Differentiation of human pluripotent stem cells into functional lung alveolar epithelial cells. Cell Stem Cell 21, 472–488.e10 (2017).
Yamamoto, Y., Korogi, Y., Hirai, T. & Gotoh, S. A method of generating alveolar organoids using human pluripotent stem cells. Methods Cell Biol. 159, 115–141 (2020).
Montoro, D. T. et al. A revised airway epithelial hierarchy includes CFTR-expressing ionocytes. Nature 560, 319–324 (2018).
Lukassen, S. et al. SARS-CoV-2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells. EMBO J. 39, e105114 (2020).
Yan, K. S. et al. Non-equivalence of Wnt and R-spondin ligands during Lgr5+ intestinal stem-cell self-renewal. Nature 545, 238–242 (2017).
Van der Velden, J. L., Bertoncello, I. & McQualter, J. L. LysoTracker is a marker of differentiated alveolar type II cells. Respir. Res. 14, 123 (2013).
Chang, J. et al. Gpr124 is essential for blood-brain barrier integrity in central nervous system disease. Nat. Med. 23, 450–460 (2017).
Nagendran, M., Riordan, D. P., Harbury, P. B. & Desai, T. J. Automated cell-type classification in intact tissues by single-cell molecular profiling. eLife 7, e30510 (2018).
Neal, J. T. et al. Organoid modeling of the tumor immune microenvironment. Cell 175, 1972–1988.e16 (2018).
Dobbs, L. G., Williams, M. C. & Brandt, A. E. Changes in biochemical characteristics and pattern of lectin binding of alveolar type II cells with time in culture. Biochim. Biophys. Acta 846, 155–166 (1985).
Van Lidth de Jeude, J. F., Vermeulen, J. L., Montenegro-Miranda, P. S., Van den Brink, G. R. & Heijmans, J. A protocol for lentiviral transduction and downstream analysis of intestinal organoids. J. Vis. Exp. 98, e52531 (2015).
Manicassamy, B. et al. Analysis of in vivo dynamics of influenza virus infection in mice using a GFP reporter virus. Proc. Natl Acad. Sci. USA 107, 11531–11536 (2010).
Karaca, G. et al. TWEAK/Fn14 signaling is required for liver regeneration after partial hepatectomy in mice. PLoS ONE 9, e83987 (2014).
We thank members of the Kuo and Desai laboratories for discussions; the Stanford Tissue Bank, J. Shrager, M. Berry and W. Trope for tissue acquisition; S. Plevritis for trajectory analysis; and the Stanford Stem Cell FACS Facility, P. Chu, A. McCormick, D. Mendoza, F. de la Vega and J. Zengel for technical expertise. SARS-Related Coronavirus 2, Isolate USA-WA1/2020, NR-52281 was deposited by the CDC and obtained through BEI Resources, NIAID, NIH. Fellowships supporting authors are as follows: A.A.S.: A.P. Giannini, ECOG-ACRIN, P. Carbone, Stanford Cancer Institute; S.S.C.: Stanford Medical Scientist Training Program grant T32 GM007365-44; J.Z.: Stanford Graduate Fellowship.; S.M.d.l.O.: CIRM Bridges. Funding support is as follows: A.R.: NIH grant T32 AI007502-23; R.A.F.: Damon Runyon Cancer Research Foundation (DRG-2286-17); V.v.U.: Netherlands Organization for Scientific Research Rubicon grant (452181214); C.S. and J.Z.: NSF DMS 1712800 and the Stanford Discovery Innovation Fund; K.C.G. and M.M.D.: Howard Hughes Medical Institute; C.A.B.: Burroughs Wellcome Fund Investigators in the Pathogenesis of Infectious Disease Grant 1016687. This work was also supported by CIRM DISC2-09637 to C.J.K. and T.J.D.: Bill and Melinda Gates Foundation OPP1113682 to C.J.K., M.R.A., and C.A.B.; Novo Nordisk Foundation Challenge Grant to M.R.A and M.M.-C.; Mathers Foundation Covid Fund to K.C.G.; and NIH grants K08DE027730 to A.A.S., U19AI057229 to M.M.D., R56AI111460 to J.S.G., 5R01HL14254902 to T.J.D., DK11572802 to C.J.K. and K.C.G., R01AI157155 to R.S.B., and U19AI116484, U01DK085527, U01CA217851, U01CA176299 and U01DE025188 to C.J.K.. C.A.B is the Tashia and John Morgridge Faculty Scholar and Chan Zuckerberg Biohub Investigator. T.J.D. is the Woods Family Endowed Faculty Scholar in Pediatric Translational Medicine. C.J.K. is the Maureen Lyles D’Ambrogio Professor of Medicine.
C.J.K., A.A.S., S.S.C., C.A.B., A.R., M.R.A., M.M.-C., S.M.d.l.O. and T.U. are listed as inventors on provisional patent 63/053,079 describing the methods in this paper. C.J.K. is a founder of Surrozen Inc. All other authors declare no competing interests.
Peer review information Nature thanks the anonymous reviewers for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
a, Schematic of culture initiation from human distal lung. b, Brightfield microscopy evaluation of required exogenous growth factors and automated organoid quantitation after day 10 of chemically defined organoid culture with specified recombinant growth factors, N = Noggin, E = EGF, W = WNT3A, R = RSPO1, n = 4 per condition, data are mean ± s.e.m., * = P < 0.05, two-tailed student’s t-test, scale bar = 500 μm. c, Top, purification schema to isolate epithelial cells from distal human lung involving negative MACS bead depletion of CD45+ haematopoietic cells, endothelial cells and fibroblasts, followed by positive FACS selection for EPCAM+ epithelium. Bottom left, representative FACS demonstrating >99.9% EPCAM+ purity (orange) upon re-analysis versus unstained controls (grey). Bottom right, proliferation of EPCAM+ cells purified from distal lung cultures after day 10 of organoid culture with specified growth factors, N = Noggin, E = EGF, W = WNT3A, R = RSPO1, n = 3 per condition, data are mean ± s.e.m., ** = P < 0.01, two-tailed student’s t-test. d, Time lapse transmission confocal images of solid and cystic organoids originating from single dissociated human distal lung cells, scale bar = 100 μm. e-g, Clonality mixing studies. e, Schema of mixing studies of lentivirus-GFP- and lentivirus-mCherry-expressing cells to determine clonality. f, Representative live fluorescent imaging of resultant green and red organoids from (e), scale bar = 500 μm. g, Quantitation of red, green, or chimaeric, distal lung organoid cultures from two individuals (1, 2) after initial and serial passaging (P1 = passage 1).
Extended Data Fig. 2 scRNA-seq of human distal lung organoids cultured from three individuals reveals reproducible basal, club, and AT2 populations.
a–c, Unsupervised clustering of total cell populations demonstrates consistency in top differentially expressed genes corresponding to basal (KRT5/6), club (SCGB1A1), and AT2 (SFTPC) cells. The epithelial fraction from these cultures ranged from approximately 90–99% of all cells with the remainder being either fibroblasts (VIM+) or mononuclear cells (HLA-DR+, likely alveolar macrophages). d–f, t-SNE visualization and violin plots for marker genes corresponding to each population. Note, a unique population enriched for SPRR genes, which have been described to mark squamous metaplasia, were exclusively found in the organoid culture of Lung 3, derived from an individual who was an active smoker.
a, SPADE plot of pooled cells where each point represents cell states that are more related on the same or adjacent branches of a minimum spanning tree. Note: AT2 cells exist on a branch distal to basal and club cells, suggesting no lineage hierarchy between AT2, basal, and club cells. b, SPADE plots of pooled scRNA-seq samples after excluding AT2, VIM+ and HLA-DR+ cells support lineage relationships between basal (blue) and club (red) populations by club cell branches emanating from basal cells. c, SPADE plots of Basal 1, Basal 2, and club populations. d, Left, gene expression of SCGB1A1 shows higher expression in club versus basal cell lineages (left) as compared to KRT5 (middle). Right, median gene expression of TNFRSF12A, showing a high (orange outline) and a low (blue outline) within basal cell branches and inferring a potential lineage relationship. e, Monocle 3 pseudotime trajectory analysis of single cell transcriptomes of Basal 1 connected to club cells depicted with UMAP (cell number, n = 3,721), colored by (left) cell cluster and (right) pseudotime gradient.
a, Left, confocal images of a live AT2 organoid at 67 days of culture labelled with Hoescht nuclear stain and LysoTracker Red DND-99. Top right, isolation of purified AT2 organoids. Representative FACS plots showing LysoTracker Red AT2 purification from unfractionated organoid cultures. Bottom right, immunostaining of cytospin of LysoTracker-sorted AT2 cells show high purity (100/100 cells SPC+ SCGB1A1- KRT5); scale bar = 50 μm). b, Schema of FACS isolation of AT2 cells from human mixed distal lung organoids as EPCAM+LysoTracker+ AT2 cells followed by long-term clonogenic organoid culture. c, Representative image of clonal mixing studies from stroma-depleted, EPCAM+LysoTracker+, lentivirally-marked AT2 cells demonstrating presence of completely mCherry+ or GFP+ but not chimaeric organoids carried out as in Extended Data Fig. 1e-g, passage 1 after lentiviral infection, scale bar = 200 μm. d, Quantitation of red, green, or chimaeric AT2 organoid cultures as in (c) from two individuals (1, 2) after initial and serial passaging (P1 = passage 1). e, AT2 organoid proliferation with differing combinations of recombinant niche factors and PORCUPINE inhibitor C59 (1 mM), NOGGIN (N), EGF (E), WNT3A (W), R-SPONDIN1 (R). n = 3 per condition, data are mean ± s.e.m., * = P < 0.05, *** = P < 0.001, two-tailed student’s t-test. f, Brightfield microscopy comparing pure AT2 organoid growth enhancement in chemically defined lung organoid media (EN) versus serum-containing L-cell conditioned media containing WNT3A, NOGGIN, and R-SPONDIN3 (L-WRN CM) supplemented with recombinant EGF, one experiment. Scale bar = 200 μm. g, Transmission electron microscopy image of representative AT2 organoid at 28 days of culture. Note apical microvilli (black arrows) and lamellar bodies (red arrows); scale bar = 10 μm.
a–c, Basal organoids in mixed culture progressively form internal lumens, which is not associated with apoptosis. a, KRT5 IF, day 26 culture, scale bar = 200 μm. b. Lumen quantitation, d12 versus d26 culture, single determination. c, Absence of apoptosis in d26 basal cell organoid internal lumen, cleaved caspase IF, from Fig. 2b, scale bar = 20 μm. d, e, Isolation of purified basal cell organoids via differential sedimentation in Ficoll. d, Schema and enrichment to >90% KRT5+ cells as measured by intracellular KRT5 FACS of sedimented basal organoid cells; scale bar = 100 μm. e, Serial time lapse microscopy of sedimented basal organoids reveals spontaneous cavitation within two weeks post passage or within four weeks of culture initiation; scale bar = 25 μm. f, g, Clonal mixing studies from stroma-depleted, Ficoll-purified and lentivirally-marked basal organoid cells demonstrating fully mCherry+ or GFP+ but not chimaeric organoids as in Extended Data Fig. 1f, passage 1 after lentiviral infection, scale bar = 200 μm. f, Representative clonal mixing image study. g, Quantitation. h, Growth factor evaluation for basal organoids after d14 sedimentation, enzymatic dissociation and clonogenic culture. Growth was not affected by the PORCUPINE inhibitor C59 (1 μM). n = 4 per condition, data are mean ± s.e.m., *** = P < 0.001, two-tailed student’s t-test.
Extended Data Fig. 6 scRNA-seq identifies an active basal cell subpopulation across three individual patient organoid cultures.
a, High resolution clustering analysis identifies a reproducible active basal cell subpopulation with significantly higher expression of mRNAs for TNFRSF12A, the NOTCH pathway marker HES1, and the proliferation marker MKI67. Modified Kruskal–Wallis Rank Sum Test two tail p-values: TNFRSF12A 4.15 × 10−8; HES1 2.4 × 10−10; MKI67 3.4 × 10−3. b, Fine resolution clustering of KRT5+ populations identifies two Basal 1 sub-clusters, Basal 1.1 and 1.2. c, Gene Ontology PANTHER overrepresentation of differentially expressed genes enriched in Basal 1.2 versus 1.1 show the majority of Basal 1.2 processes involve cell cycle (asterisks). Complete analysis is provided in Supplementary Table 3. d, Violin plot of scRNA-seq analysis for Fig. 2f–h depicting KRT5 mRNA expression among triply EPCAM+ITGA6+ITGB4+ mRNA-expressing single cells (purple, that is, tandem mRNA expression of all three genes) versus the remainder of cells (grey), P < 0.001 two tailed Kruskal–Wallis Rank Sum Test. e, t-SNE visualization of TNFRSF12A and ITGA6 expression from d among cells with EPCAM+ITGA6+ITGB4+ gene expression and subdivision by high (top quartile, orange), medium (pink) and low (bottom quartile, navy blue) mRNA expression. f, Proliferation-associated gene expression is progressively enriched for scRNA-seq cell fractions of in EPCAM+ITGA6+ITGB4+ cells that are stratified for low, medium, or high expression of TNFRSF12A mRNA as in e. Data in f represent cell population fractions from a single experiment. *** = P < 0.001, two tailed Chi-square test.
Extended Data Fig. 7 Evaluation of Basal 1 lineage relationship to Basal 2 and the influence of NOTCH signalling on Basal 1 renewal and differentiation.
a, Isolation of Basal 1 and Basal 2 via differential sedimentation of KRT5+ cells followed by FACS sorting of EPCAM+ITGA6+ITGB4+TNFRSF12A+ (Basal 1) versus EPCAM-ITGA6-ITGB4-TNFRSF12A- (Basal 2). b, Intracellular FACS measurement of KRT5 protein expression in Basal 1 and 2 fractions from a. c, Representative brightfield of day 14 organoid cultures from a, b. d, Quantitation of 3 independent experiments from a–c, box plot represents first quartile, median, third quartile, and whiskers represent minimum and maximum. *** P < 0.001, two tailed student’s t-test. e, qPCR measurement of two differentially upregulated Basal 2 genes from the three scRNA-seq biological replicates (Extended Data Fig. 6, SPRR1B, TMSB4X) after prolonged culture of FACS isolated Basal 1 cells. Data are relative mean ± s.e.m. of cultures from three independent experiments, ** = P < 0.01, two tailed student’s t-test. f, RNA FISH demonstrating TMSB4X and SPRR1B cellular transcripts within organoids originating from Basal 1 cells (arrows), scale bar = 25 μm. g, KRT5 immunostaining and SFTPC and SCGB1A1 RNA FISH of FACS isolated TNFRSF12Ahi Basal 1 cells under vehicle, NOTCH agonism (JAG1 peptide), or NOTCH antagonism with the Delta-like ligand mutant 4 (DLL4E12 ; E12) or the gamma secretase inhibitor DBZ; scale bar = 50 μm. h, Fluorescent quantitation of resazurin dye reduction to estimate relative cellular proliferation in g, data are normalized to vehicle (V) and represent mean ± s.e.m. from five independent experiments, * P < 0.05, two tailed student’s t-test. i, Quantitation of SCGB1A1 and SFTPC gene expression by RNA FISH in the context of NOTCH agonism or antagonism from three independent experiments, ** P < 0.01, *** P < 0.001, two tailed student’s t-test. SFTPC mRNA upregulation was not accompanied by lamellar body or SFTPC protein production (data not shown).
a, b, Immunostaining of KRT5 and TNFRSF12A in human distal airways from two individuals, scale bar = 100 μm. c, Immunostaining of KRT5, TNFRSF12A, and p63 in human distal airway, scale bar = 100 μm. d, FACS analysis from freshly fixed human distal lung with anti-KRT5 (intracellular) and monoclonal anti-TNFRSF12A (cell surface) (top), or sequential FACS isolation from freshly dissociated human distal lung of EPCAM+ITGA6+ITGB4+ cells followed by fractionation into TNFRSF12Ahi or TNFRSF12Aneg subsets (bottom), pre-gated on live singlets and used for culture experiments in Fig. 3d–g.
a, b, Distal lung organoid modelling of H1N1 influenza infection. a, Merged transmission and GFP confocal images of purified basal (left) and purified AT2 organoids (right) 12 h after infection with PR8-GFP H1N1 influenza virus, quantified by FACS for % GFP+ cells. Bar plot represents the mean percentage of infected cells from three technical replicates, P = 0.57, Chi-square test. Scale bars = 50 μm. b, Viral genome quantitation over time of mixed distal lung organoid culture supernatants subjected to initial infection of wild-type H1N1 at an estimated multiplicity of infection (MOI) of 0.01, qRT–PCR, data represent the mean of three independent experiments ± s.e.m. c, d, Lectin staining with M. amurensis (a2-3) and S. nigra (a2-6) lectins or no lectin negative controls to characterize sialic acid residues which serve as surface molecules for influenza virus host cell entry. AT2 organoids (c) and basal organoids (d). Scale bar = 25 μm. e, Dose response curves for two different classes of antiviral drugs on influenza infectivity and replication. The nucleoside analogue FdC demonstrated dose dependent activity with IC50 of 340 nM as compared to neuraminidase inhibitor zanamivir, which only impairs viral shedding but not infectivity and replication. n = 3 per condition, data represents mean ± s.e.m. f, Fluorescence micrograph of multiwell screening of selected various antiviral agents after H1N1 PR8-GFP organoid infection in 48 well format. FdC = nucleoside analogue 2'-deoxy-2'-fluorocytidine. Cpd = compound #.
Extended Data Fig. 10 Apical-out polarization and multi-lineage differentiation of distal lung organoids upon suspension culture.
a, scRNA-seq plots of ACE2 and TMPRSS2 gene expression in ECM-embedded mixed distal lung organoids as in Fig. 1a–h. b, Diagram of ECM removal and suspension culture leading to apical-out polarity of lung organoids. c, Representative confocal microscopy showing reorganization of microfilaments (phalloidin) and acetylated microtubules (AcTUB) upon ECM removal. Scale bar = 10 μm. d–f, Polarization and accelerated ciliary differentiation of apical-out basal organoids. d, Confocal 3D sections (top panels) and surface reconstructions (bottom panels) of apical-out lung organoids at different days after ECM removal. At day 0 (d0) microfilament (green, phalloidin) and microtubule (red, acetylated tubulin) organization is not polarized while junctional strands (ZO-1, white) are polarized. By day 2 in suspension (d2) ZO-1 (white) forms junctional rings in the apical periphery of each cell facing the external side of the organoids and the actin cytoskeleton forms microvilli (green) facing outward (apical-out polarity). Also, at d2 some cells initiate microtubule polarization. By day 5 (d5) many more cells have motile cilia facing outwards. Mature motile cilia can be observed for several weeks, example at day 14 (d14). e, 3D confocal reconstruction of an organoid embedded in ECM consisting mostly of basal stem cells (KRT5+, white). f, As apical-out polarity is established in suspension culture and ciliogenesis begins, KRT5+ basal cells are found underneath the polarized epithelium. g, SCGB1A1+ Club cells with apical-out polarity are present on the exterior surface. In all panels nuclei are stained blue with DAPI, and actin microfilament organization visualized with phalloidin (green). Scale bars = 10 μm. h–j, Prolonged suspension culture of AT2 organoids (day 10 post-suspension) induces apical-out polarization and AT1 differentiation. h, Optical sections through alveolar-derived organoids after 10 days in suspension culture show decreased abundance of AT2 cells while individual cuboidal cells begin to express the AT1 marker HTI-56 (red), a transmembrane protein specific to the apical membrane of alveolar type 1 pneumocytes (AT1). Scale bars = 10 μm. i, j, Side views of alveolar organoids after 10 days of suspension culture reveal thin AT1 cells with phalloidin-reactive apical junctional complexes facing outwards (apical-out) (i) and expression of HTI-56 on the apical membrane (j). Scale bars = 10 μm. k, Representative confocal microscopy immunofluorescence of apical-out human basal organoid after 10 days in suspension expressing SARS-CoV-2 receptor ACE2 (green), cilia (AcTUB, red) and DAPI (blue). Scale bar = 20 μm.
Gating strategies. a, EPCAM+Lysotracker+ gating strategy to purify AT2 cells. Purity was confirmed with 100/100 sorted cells staining positive for SFTPC protein (corresponding to Extended Data Fig. 4a). b, EPCAM+ITGA6+ITGB4+TNFRSF12Ahi/neg gating strategy. Purity was confirmed with 100/100 pooled sorted cells staining positive for KRT5 protein (corresponding to Fig. 2j-k). c, Gating strategy to analyze KRT5+ cells and corresponding TNFRSF12A+ populations within cells freshly dissociated from distal human lungs (corresponding to Extended Data Fig. 8d). d, RT-PCR correlation of TNFRSF12A and other mRNAs from cells sorted by FACS into TNFRSF12A neg, med, and hi fractions (Fig. 3d). e, Gating strategy to estimate H1N1 PR8 GFP infectivity in human lung organoids (corresponding to Extended Data Fig. 9a).
Clinical demographics of 136 lung tissue donors used in this study. Each row corresponds to a unique individuals from whom normal lung was obtained and columns correspond to features such as age, sex, and lung anatomic location.
Genes used for SPADE analysis. Genes were utilized in SPADE based on graph-based clustering annotation of cells from scRNA-seq analysis of three individuals. Full SPADE details are provided in Supplementary Methods.
Gene Set Enrichment Analysis (GSEA) of top differentially expressed genes from scRNA-seq populations. Populations correspond to Basal 1, Basal 2, Basal 1.1, Basal 1.2, and Cluster 1 of purified AT2 scRNA-seq data sets. Top differentially expressed genes for each population are listed in each tab. Bold = genes used for GSEA analysis (Panther overrepresentation test, with Bonferroni correction, p values are two tailed).
List of Basal 1 membrane protein genes. Genes were identified by ontology analysis of top differentially expressed Basal 1 genes with GO term 0031224, (intrinsic component of membrane).
Analysis of SARS-CoV-2-infection in lung organoid cultures. Positively infected cells were quantified by dsRNA and SARS-CoV-2 NP staining among KRT5+ basal cells, SCGB1A1+ club cells, and AcTUB+ ciliated cells using confocal microscopy. Statistical analysis (chi-square or Fisher’s exact test) were performed as indicated.
Next Generation Sequencing of five distal lung organoid cultures. Summary and annotations of nonsynonymous variants with allele frequency > 0.10 for 130 cancer genes using the TOMA Tumor Profiling System (TOMA, Foster City, CA).
Tables of reagents. Individual tables of organoid media components, immunostaining and fluorescence microscopy reagents, RT-PCR reagents, RT-PCR primers, and RNA hybridization probes used in this study.
Guide to lung organoid scRNA-seq analysis. scRNA-seq analysis for lungs 1, 2, 3 and purified AT2 organoids provided in .rmd and .html format.
Lung organoid culture sequencing. Next Generation Sequencing data for each of five distal lung organoid cultures (TOMA Tumor Profiling System) are provided in Variant Call Format.
Differentiated basal organoids have functional cilia, Part 1. Transmission confocal microscopy video of beating cilia (25X magnification).
Differentiated basal organoids have functional cilia, Part 2. Bright-field microscopy video of beating cilia (10x magnification).
Apical-out differentiated basal organoids have functional cilia. DIC and 3D confocal microscopy video of motile cilia in apical-out basal organoids at suspension culture day 5 and 14.
About this article
Cite this article
Salahudeen, A.A., Choi, S.S., Rustagi, A. et al. Progenitor identification and SARS-CoV-2 infection in human distal lung organoids. Nature 588, 670–675 (2020). https://doi.org/10.1038/s41586-020-3014-1
Veterinary Research (2021)
Cell Regeneration (2021)
Experimental & Molecular Medicine (2021)
Nature Cell Biology (2021)