An organoid-derived bronchioalveolar model for SARS-CoV-2 infection of human alveolar type II-like cells.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), which may result in acute respiratory distress syndrome (ARDS), multiorgan failure, and death. The alveolar epithelium is a major target of the virus, but representative models to study virus host interactions in more detail are currently lacking. Here, we describe a human 2D air-liquid interface culture system which was characterized by confocal and electron microscopy and single-cell mRNA expression analysis. In this model, alveolar cells, but also basal cells and rare neuroendocrine cells, are grown from 3D self-renewing fetal lung bud tip organoids. These cultures were readily infected by SARS-CoV-2 with mainly surfactant protein C-positive alveolar type II-like cells being targeted. Consequently, significant viral titers were detected and mRNA expression analysis revealed induction of type I/III interferon response program. Treatment of these cultures with a low dose of interferon lambda 1 reduced viral replication. Hence, these cultures represent an experimental model for SARS-CoV-2 infection and can be applied for drug screens.

Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS‐n class="Species">CoV‐2) has spread globally within several months after an initial outbreak in Wuhan, China, in December 2019 (Zhu et al, 2020). The World Health Organization (WHO) declared SARS‐CoV‐2 a pandemic on March 11, 2020. As of November 1, 2020, > 46,000,000 patients were confirmed including > 1,200,000 deaths (WHO, 2020). SARS‐CoV‐2 belongs to the Sarbecovirus subgenus (genus Betacoronavirus, family Coronaviridae), together with SARS‐CoV (Andersen et al, 2020; Coronaviridae Study Group of the International Committee on Taxonomy of V, 2020; Lu et al, 2020). SARS‐CoV‐2 causes coronavirus disease 2019 (COVID‐19), an influenza‐like illness associated with a broad spectrum of clinical respiratory syndromes, ranging from mild upper airway symptoms to a life‐threatening viral pneumonia (Chen et al, 2020; Wang et al, 2020; Wolfel et al, 2020; Xu et al, 2020). This pneumonia fulfills the radiographic and histological criteria for acute respiratory distress syndrome (ARDS). Classically, ARDS is caused by damage to the alveoli, more specifically to the alveolar type I cells that are critical in the oxygen transport from the alveoli to the blood, triggered through an inflammatory response to pathogens or toxins (Ware & Matthay, 2000; Thompson et al, 2017; Ackermann et al, 2020; Barton et al, 2020; Menter et al, 2020; Raptis et al, 2020).

The early events that lead to respiratory failure as a result of n class="Disease">coronavirus infection are initiated following virus replication (tenOever, 2016). After virus entry, which is determined by receptor availability, viral replication generates pathogen‐associated molecular patterns, specifically RNA structures, that can be detected by pattern recognition receptors (PRRs). PRR signaling subsequently leads to the transcription of antiviral genes; the response triggered varies between pulmonary cell types (preprint: Ravindra et al, 2020). Therefore, disease outcome can be greatly influenced by the viral target cell, which is determined largely by the presence of CoV entry receptor. In the lungs, angiotensin‐converting enzyme 2 (ACE2), the SARS‐CoV‐2 receptor (Hoffmann et al, 2020; Walls et al, 2020; Wrapp et al, 2020; Zhou et al, 2020), is expressed mainly in ciliated cells and alveolar type II cells (Jia et al, 2005; Hikmet et al, 2020; Qi et al, 2020). Ciliated cells, alveolar type II cells, but also alveolar type I cells, have been identified as SARS‐CoV‐2 target cells in animal models and in COVID‐19 patients (Hou et al, 2020; Rockx et al, 2020; Xu et al, 2020). In previous work, we confirmed ciliated cells as a viral target cell using organoid‐derived bronchial airway cultures (Lamers et al, 2020). In contrast to ciliated cells, alveolar cells are notoriously difficult to culture, thus limiting our understanding of COVID‐19, but also of other respiratory virus infections. Currently, standardized methods use immortal cell lines or primary alveolar cells to study disease (van den Bogaard et al, 2009). Immortal cell lines, however, do not fully recapitulate the complexity of the alveolar space and epithelium. Primary alveolar cells partially capture this complexity but remain incapable of undergoing passaging and quickly lose their in vivo phenotype which restricts the model by the availability of donor material (Logan & Desai, 2015; Zacharias et al, 2018). Importantly, primary adult alveolar cultures were recently shown to be poorly permissive to SARS‐CoV‐2 infection, emphasizing the urgent need to develop a susceptible SARS‐CoV‐2 alveolar infection model (Hou et al, 2020; Hui et al, 2020). Some induced pluripotent stem cell‐derived models have been able to show alveolar differentiation (Gotoh et al, 2014; Dye et al, 2015; Jacob et al, 2017; Yamamoto et al, 2017; de Carvalho et al, 2019; van Riet et al, 2020), yet differentiation from primary lung cells has remained challenging.

Here, we describe a bronchioalveolar‐like 2D air–liquid interface cell system by confocal and electron microscopy and mRNA expressioene">n aene">nalysis. Iene">n this model, alveolar cells, but also basal cells, aene">nd rare neuroeene">ndocrine cells, are growene">n from 3D self‐reene">newing luene">ng bud tip progeene">nitor orgaene">noids. Iene">n this system, SARS‐n class="Species">CoV‐2 replication competence, tropism, and induced host responses were determined and compared to 2D differentiated small airway epithelium.

Results

Replication competence and tropism of SARS‐CoV‐2 in 2D differentiated small airway epithelium

Human broene">nchial airway epithelial (HAE) cell cultures are aene">n established primary epithelial luene">ng cell model to study respin class="Species">ratory virus infections. In these cultures, primary basal cells can be differentiated at air–liquid interface into mature airway cell types. It has already been shown that 3D self‐renewing airway organoids can also be used as a source of basal cells in this culture system (Zhou et al, 2018; Sachs et al, 2019; Lamers et al, 2020). We first established this system using small airway basal cells (Appendix Fig S1A–C) and demonstrated that SARS‐CoV‐2 readily infects human small airway organoid‐derived epithelium cultured in 2D at air–liquid interface (ALI) (Fig 1). SARS‐CoV‐2 grew to relatively high titers on these cells, as shown by titration of VeroE6 cells (Fig 1A and B) and viral RNA quantification (Appendix Fig S1D). Shedding of virus occurred predominantly apically (Appendix Fig S1E and F). As we have previously shown for large airway (bronchial) cultures (Lamers et al, 2020), ciliated cells were extensively targeted by SARS‐CoV‐2 (Fig 1C). In addition, we noted rare infection of CC10+ club cells (Fig 1D), but no infection of MUC5AC+ goblet cells (Fig 1E).

Figure 1

SARS‐CoV‐2 infects human organoid‐derived 2D small airway cultures

Infectious virus titers can be observed by virus titrations on VeroE6 cells of apical washes at 2, 24, 48, and 72 h after infection at MOI 0.01 (A) or 0.1 (B) with SARS‐CoV‐2 (red). The dotted line indicates the lower limit of detection. Error bars represent SEM. N = 4. H p.i. = hours post‐infection.

Immunofluorescent staining of SARS‐CoV‐2 infected differentiated small airway cultures. Nucleoprotein (NP) stains viral capsid (red), which colocalized with the ciliated cell marker AcTUB (green; (C, D, E)) and club cell marker CC10 (blue; (D)). Phalloidin was included to stain actin (blue; (C)). Goblet cells are identified by MUC5AC (blue; (E)). Scale bars indicate 50 μm.

Data information: Nuclei are stained with Hoechst or TOPRO and shown in white (C‐E). Scale bars indicate 50 μm.

Infectious virus titers can be observed by virus titratioene">ns oene">n VeroE6 cells of apical washes at 2, 24, 48, aene">nd 72 h after n class="Disease">infection at MOI 0.01 (A) or 0.1 (B) with SARS‐CoV‐2 (red). The dotted line indicates the lower limit of detection. Error bars represent SEM. N = 4. H p.i. = hours post‐infection.

Immunofluorescent staining of SARS‐CoV‐2 n class="Disease">infected differentiated small airway cultures. Nucleoprotein (NP) stains viral capsid (red), which colocalized with the ciliated cell marker AcTUB (green; (C, D, E)) and club cell marker CC10 (blue; (D)). Phalloidin was included to stain actin (blue; (C)). Goblet cells are identified by MUC5AC (blue; (E)). Scale bars indicate 50 μm.

Establishment of a 2D differentiated bronchioalveolar‐like model

Another 3D lung organoid system was developed by Nikolic et al (2017), but this system has not yet beeene">n applied in virology. Iene">n this system, culture coene">nditioene">ns were established to support loene">ng‐term self‐reene">newal of multipoteene">nt n class="Gene">SOX2+SOX9+ lung bud tip progenitor cells which in vivo differentiate into both airway and alveolar cells. We grew these lung bud tip organoids (LBT) from canalicular stage human fetal lungs 16–17 pcw (post‐conception weeks). In expansion medium, which activates EGF, FGF, and WNT signaling, and inhibits BMP and TGFβ, the vast majority of cells were SOX2+SOX9+ (Fig 2A), but rare ATII‐L were also detected in a subpopulation of organoids using the HTII‐280 antibody which exclusively stains ATII cells in the human lung (Gonzalez et al, 2010) (Fig 2B). Organoid lines were maintained for > 14 passages without apparent change in morphology, or SOX2 and SOX9 expression. Next, we established a model that contains alveolar‐like cells, which could be accessed from the apical side for infection experiments. We plated SOX2+SOX9+ progenitor cells in Transwell inserts in differentiation medium. Based on recent literature, we reasoned that addition of human canalicular stage mesenchyme would provide cues for the survival of ATII‐L cells and differentiation toward an alveolar fate (Barkauskas et al, 2013; Nikolic et al, 2017; Nikolic & Rawlins, 2017; Leeman et al, 2019). After reaching confluency, cells were differentiated for at least 14 days at ALI (Fig 2C). The addition of canalicular stage mesenchymal cells in the bottom compartment of the Transwell increased the frequency of HTII‐280+ cells (Appendix Fig S2A and B). After 14 days of differentiation at ALI, cultures consisted of both multilayered and single‐layered squamous epithelium. The areas with a single layer of mostly thin epithelium contained cells expressing the ATI markers HOPX (Fig 2D) and HTI‐56 (Fig 2E), as well as the mature ATII marker HTII‐280 (Fig 2F) and LPCAT1 (Fig 2G) (Dobbs et al, 1999; Gonzalez et al, 2010; Gonzales et al, 2015). SFTPC+ cells were present in the single‐layered epithelium and the top layer of the multilayered areas (Fig 2H). The multilayered areas contained TP63+ basal cells as a bottom layer (Fig 2I). Canalicular stage mesenchyme from different donors resulted in similar expression of alveolar markers (Fig EV1, EV2, EV3, EV4, EV5), whereas some variation was observed in basal cell proliferation (Fig 1EVI and J). Using mass spectrometry on apical washes, we also detected secreted surfactant proteins SFTPA1, SFTPA2, SFTPB, and SFTPD (Appendix Fig S3). mRNA expression analysis of 3D LBT and confluent 2D cultures before (2D non‐ALI) and after (2D ALI) exposure to air indicated that lung bud tip progenitor markers SOX9 and SOX2 were decreased during differentiation (Fig EV2, EV3, EV4, EV5A and B), while alveolar cell markers SFTPA2, SFTPB, and HOPX gradually increased during differentiation (Fig EV2, EV3, EV4, EV5). A general increase in ATI‐ and ATII‐related genes from 3D LBT to 2D ALI cells was observed with intermediate levels at 2D non‐ALI (Fig EV2F, Tables EV1 and EV2). The composition of cells in 2D ALI cultures was different from small airway ALI cultures, as indicated in differential expression analysis (Fig 3A, Table EV3). While small airway ALI cultures were mainly expressing ciliated cell markers like DNAH10 (Fig 3A), FOXJ1, and SNTN (Fig 3B), alveolar‐like ALI cultures were enriched in ATI&II makers like HOPX and SFTPA1 (Fig 3). Moreover, 2D ALI cultures were enriched in markers for pulmonary neuroendocrine, tuft, and basal cells, while small airway cultures were enriched in goblet and ciliated cell markers (Fig 3B). All together, these data indicate that the 2D ALI cultures contain both alveolar‐like and bronchiolar‐like cells. These cultures could therefore be characterized as bronchioalveolar‐like by confocal microscopy and mRNA expression analysis. The presence of pulmonary neuroendocrine cells (PNEC) was confirmed by confocal imaging (Appendix Fig S6B). The club and ATII cell marker SCGB3A2 was also detected in both single‐ and multilayered areas (Appendix Fig S6C).

Figure 2

Self‐renewing human fetal lung bud tip progenitor organoids differentiate at air–liquid interface to alveolar type I‐ and type II‐like cells

Immunofluorescent staining of fetal lung bud tip progenitor organoids grown in expansion medium co‐expressing stem cell markers Sox2 (green) and Sox9 (red). Phalloidin (white) was used to stain actin.

Immunofluorescent staining of rare HTII‐280 + type II pneumocytes (green) in a subpopulation of fetal lung bud tip progenitor organoids grown in expansion medium.

Schematic of the 2D air–liquid interface bronchioalveolar‐like model.

Differentiated lung bud tip organoids at air–liquid interface in co‐culture with donor‐specific human fetal lung fibroblasts. After 14 days of differentiation at air–liquid interface, cells express alveolar type I (HOPX, green, (D); HTI‐56, green, (E)), type II cells (HTII‐280, green, (F); LPCAT1, green, (G); SFTPC (SPC), green, (H)), and basal cell (TP63 (P63), green, (I)) markers in areas containing one cell layer. Dotted lines indicate the barrier between multilayered and single‐layered epithelium.

Data information: Nuclei are stained with Hoechst (blue in (A, B) or white in (D–I)). Scale bars indicate 50 μm.

Figure EV1

A comparison of different mesenchymal donors on 2D bronchioalveolar‐like differentiation

Immunofluorescent staining of HTII‐280 (A, B), LPCAT1 (C, D), SPC (E, F), HOPX (G, H), and TP63 (I, J) after 2 weeks of differentiation at air–liquid interface in the presence of mesenchyme in the basal compartment of the Transwell. Nuclei are stained with Hoechst (blue). Scale bars represent 50 μm.

Figure EV2

Transcriptomic analysis of differentiation toward bronchioalveolar‐like cells from lung bud tip organoids

Line graphs indicating normalized read count of 3D lung bud tip organoids (3D LBT), airway organoids on transwell system with apical medium (2D non‐ALI), and 2D bronchioalveolar‐like air–liquid cultures (2D ALI). n = 2 and error bars represent stdev. (A, B) Graphs depicting fetal bud tip progenitor markers SOX9 and SOX2. (C, D) Graphs depicting ATII markers SFTPA2 and SFTPB. (E) Graphs depicting ATI marker HOPX.

Heatmaps depicting top 50 enriched genes in 2D bronchioalveolar‐like air–liquid cultures (2D ALI) compared to 3D lung bud tip organoids (3D LBT). N = 2 of 3D LBT, airway organoids on transwell system with apical medium (2D non‐ALI), and 2D ALI. Colored bars represent Z‐scores of log2‐transformed values.

Figure EV3

Comparison of bronchioalveolar culture to published datasets of AT2 cells

Heatmaps depicting marker genes in 2D small airway air–liquid cultures (small airway ALI), bronchioalveolar‐like air–liquid cultures (alveolar ALI), and published dataset of purified AT2 cells (purified AT2). Colored bars represent Z‐scores of log2‐transformed values.

Barplot indicating percentage of annotated cells within a cluster. Color represents annotated cell identity. Combined current dataset and of GSM2855485 are presented. Indicated annotation is based on separate clustering. Published dataset GSM2855485 is annotated as “Purified AT2”.

Violin plots visualizing expression levels of marker genes for AT2 cells. AT2‐like cells comprise cells annotated as AT2‐like cells in current dataset (clusters 0 and 2). Purified AT2 cells are from published dataset GSM2855485.

Data information: Error bars represent SEM. N = 4 for donor 1 and N = 3 for donor 2. H p.i. = hours post‐infection.

Figure EV4

Apical viral RNA and basal shedding in differentiated 2D bronchioalveolar‐like cultures

Viral SARS‐CoV‐2 RNA titers of apical washes at 2, 24, 48, and 72 h after infection at MOI 0.1 for donor 1 (A) and donor 2 (B).

Infectious virus titers as observed by virus titrations on VeroE6 cells (C) and qRT–PCR (D) in basal compartments. The dotted line indicates the lower limit of detection.

Data information: Error bars represent SEM. N = 4 for donor 1 and N = 3 for donor 2. H p.i. = hours post‐infection.

Figure EV5

Transmission electron microscopy analysis of 2D bronchioalveolar‐like (A–F) and 2D small airway (G–K) cultures after 72 h of SARS‐CoV‐2 infection

Pulmonary alveolar type II‐like (ATII‐L) cells were found in the 2D bronchioalveolar‐like culture (A, B). The early infection of the ATII‐L cell (C) is recognizable by the double‐membranes vesicles (Dmv) and virus particles (circles). (D) Next to the lamellar bodies (Lb), the ATII‐L also contained a lysosomal compartment with internal virus particles (circle). Extracellular virus particles were found in the intercellular space between the ATI‐L and ATII‐L cells (E), albeit at low concentration. Dashed circles indicate viral particles outside of the cell. The remnants of a dead detached cell (F) mainly showed apical‐shed virus particles (in areas marked by dotted lines), whereas the basal side was relatively free of virus particles.

Ciliated cells (CC; (H)) and goblet cells (GC; (K)) are visible in the cross section through the 2D small airway culture after 8 weeks of differentiation. SARS‐CoV‐2 replication (I) and virus shedding (J) were found to be associated particularly with ciliated cells, whereas no infections were found in goblet cells (K). Dashed circles indicate viral particles inside and outside the cell.

Data information: Scale bars represent 5 μm (A, G), 2.5 μm (B, H, K), 1 μm (F), 500 nm (D, I, J), and 250 nm (C, E).

Figure 3

Transcriptomic analysis of bronchioalveolar‐like and small airway air–liquid interface cultures

Heatmaps depicting top 30 up‐ and downregulated genes in 2D bronchioalveolar‐like air–liquid cultures (bronchioalveolar) compared to small airway air–liquid interface cultures (small airway). Colored bars represent Z‐scores of log2‐transformed values.

Heatmaps depicting pulmonary cell type marker genes in 2D bronchioalveolar‐like air–liquid cultures (bronchioalveolar) and small airway air–liquid interface cultures (small airway). Color‐coded bars indicate which cell type is marked by the presented gene. Colored bars represent Z‐scores of log2‐transformed values.

Immunofluorescent staining of fetal lung bud tip progenitor organoids grown in expansion medium co‐expressing stem cell markers Sox2 (greeene">n) aene">nd n class="Gene">Sox9 (red). Phalloidin (white) was used to stain actin.

Differentiated lung bud tip organoids at air–liquid interface in co‐culture with donor‐specific human fetal luene">ng fibroblasts. After 14 days of differeene">ntiatioene">n at air–liquid interface, cells express alveolar type I (n class="Gene">HOPX, green, (D); HTI‐56, green, (E)), type II cells (HTII‐280, green, (F); LPCAT1, green, (G); SFTPC (SPC), green, (H)), and basal cell (TP63 (P63), green, (I)) markers in areas containing one cell layer. Dotted lines indicate the barrier between multilayered and single‐layered epithelium.

Immunofluorescent staining of HTII‐280 (A, B), LPCAT1 (C, D), n class="Gene">SPC (E, F), HOPX (G, H), and TP63 (I, J) after 2 weeks of differentiation at air–liquid interface in the presence of mesenchyme in the basal compartment of the Transwell. Nuclei are stained with Hoechst (blue). Scale bars represent 50 μm.

Line graphs indicating normalized read count of 3D lung bud tip organoids (3D LBT), airway organoids on transwell system with apical medium (2D non‐ALI), and 2D bronchioalveolar‐like air–liquid cultures (2D ALI). n = 2 and error bars represent stdev. (A, B) Graphs depicting fetal bud tip progenitor markers SOX9 aene">nd n class="Gene">SOX2. (C, D) Graphs depicting ATII markers SFTPA2 and SFTPB. (E) Graphs depicting ATI marker HOPX.

Heatmaps depicting top 50 enriched genes in 2D bronchioalveolar‐like air–liquid cultures (2D ALI) compared to 3D lung bud tip organoids (3D LBT). N = 2 of 3D LBT, airway orgaene">noids oene">n traene">nswell system with apical medium (2D noene">n‐ALI), aene">nd 2D ALI. Colored bars represeene">nt Z‐scores of log2‐traene">nsformed values.

Barplot indicating percentage of annotated cells within a cluster. Color represeene">nts aene">nnotated cell ideene">ntity. Combined curreene">nt dataset aene">nd of n class="Chemical">GSM2855485 are presented. Indicated annotation is based on separate clustering. Published dataset GSM2855485 is annotated as “Purified AT2”.

Violin plots visualizing expression levels of marker genes for AT2 cells. AT2‐like cells comprise cells annotated as AT2‐like cells in current dataset (clusters 0 aene">nd 2). Purified AT2 cells are from published dataset n class="Chemical">GSM2855485.

Data information: Error bars represent SEM. N = 4 for doene">nor 1 aene">nd n class="Gene">N = 3 for donor 2. H p.i. = hours post‐infection.

Infectious virus titers as observed by virus titratioene">ns oene">n VeroE6 cells (C) aene">nd qRT–PCR (D) in basal compartmeene">nts. The dotted line indicates the lower limit of detectioene">n.

Data information: Error bars represent SEM. N = 4 for doene">nor 1 aene">nd n class="Gene">N = 3 for donor 2. H p.i. = hours post‐infection.

Pulmonary alveolar type II‐like (ATII‐L) cells were found in the 2D bronchioalveolar‐like culture (A, B). The early infection of the ATII‐L cell (C) is recogene">nizable by the double‐membraene">nes vesicles (Dmv) aene">nd virus particles (circles). (D) n class="Gene">Next to the lamellar bodies (Lb), the ATII‐L also contained a lysosomal compartment with internal virus particles (circle). Extracellular virus particles were found in the intercellular space between the ATI‐L and ATII‐L cells (E), albeit at low concentration. Dashed circles indicate viral particles outside of the cell. The remnants of a dead detached cell (F) mainly showed apical‐shed virus particles (in areas marked by dotted lines), whereas the basal side was relatively free of virus particles.

Ciliated cells (CC; (H)) and goblet cells (GC; (K)) are visible in the cross section through the 2D small airway culture after 8 weeks of differentiation. SARS‐CoV‐2 replicatioene">n (I) aene">nd virus shedding (J) were fouene">nd to be associated particularly with ciliated cells, whereas no n class="Disease">infections were found in goblet cells (K). Dashed circles indicate viral particles inside and outside the cell.

To further investigate the diversity and frequency of the cell types in the bronchioalveolar system, we performed single‐cell mRNA sequeene">ncing (Fig 4). This revealed a large alveolar‐like n class="Species">cluster (94.3% of cells), a smaller basal cell cluster (4.9% of cells), and a rare population of PNECs (0.9% of cells) (Fig 4A). The alveolar cluster contained ATII‐like cells (46.0% of cells, clusters 0 and 2) (Fig 4B–D; SFTPC, SFTPA1, LPCAT1; Appendix Fig S4A), proliferating ATII‐like cells (6.4% of alveolar cells, cluster 4) (Fig 4E), and alveolar cells that showed an increased expression of several ATI markers (6.5% of alveolar‐like cells, cluster 3) (Fig 4F–H; AGER, CAV1, HOPX; Appendix Fig S4B). The basal cell cluster (cluster 5) was characterized by KRT5 and TP63 expression (Fig 4I and J; Appendix Fig S4C). A rare PNEC population (cluster 7) expressed CHGA (Fig 4K; Appendix Fig S4D). To investigate whether the cultured cells were representative of adult alveolar cells, we compared our bulk RNA sequencing and single‐cell RNA sequencing dataset to previously published bulk (Fig EV3A) and single‐cell sequencing datasets (Figs 4L and M, and EV3B and C) of freshly isolated adult HTII‐280+ ATII cells. The freshly isolated adult ATII cells did not form a separate cluster in the t‐SNE map (Fig 4L and M), and overlapped with cultured alveolar, proliferating alveolar and ATII‐like cells (Figs 4M and EV3B). The expression of a set of ATII marker genes was also compared between the cultured ATII‐like cluster and the freshly purified ATII cells (Fig EV3C). The SARS‐CoV‐2 entry receptor ACE2 was detected at low levels in the alveolar cluster (Fig 4N, Appendix Fig S5A), and ACE2 expression in ATII‐like cells was confirmed by immunofluorescent imaging (Fig 4P). The activating protease TMPRSS2 was detected in all clusters with a higher expression in ATII‐like cells (Fig 4O). In addition, there was a slight trend toward a correlation between ACE2 and TMPRSS2 expression (Appendix Fig S5A–D). At the protein level, both ACE2 and TMPRSS2 were readily detected and often coexpressed (Fig 4Q). Other TMPRSS‐family members were detected in lower levels (Appendix Fig S5E–H). These data show that the bronchioalveolar model consists mainly of ATII‐like cells and potentially is susceptible to SARS‐CoV‐2 infection.

Figure 4

Single‐cell analysis of bronchioalveolar cultures

t‐SNE maps showing clustering based on differentially expression analysis and highlighting cell type identification of bronchioalveolar cultures based on cell makers. Each color represents a cluster or different cell type. Number indicates number of cells in cluster identified as cell type presented.

t‐SNE maps highlighting the expression of marker genes for each of the cell types. Transcript counts are given in a Log2 scale.

t‐SNE map showing clustering based on differentially expression analysis including dataset from (A) and published dataset of GSM2855485. Each color represents a single cluster.

t‐SNE map showing clustering of (L) with annotation of cell types. Cell type annotation is transferred from (A), and published dataset GSM2855485 is annotated as “Purified AT2”. Number indicates number of cells in cluster identified as cell type presented.

t‐SNE maps highlighting the expression of SARS‐CoV‐2 entry receptors ACE2 (N) and TMPRSS2 (O). Transcript counts are given in a Log2 scale.

Differentiated bronchioalveolar cultures express SARS‐CoV‐2 entry receptor ACE2 (green) on ATII‐like cells marked by HTII‐280 (red). Nuclei are stained with Hoechst and shown in white. Scale bar indicates 50 μm.

Bronchioalveolar cultures co‐express both entry receptors ACE2 (green) and TMPRSS2 (red). Nuclei are stained with Hoechst and shown in white. Scale bar indicates 50 μm.

t‐SNE maps showing n class="Species">clustering based on differentially expression analysis and highlighting cell type identification of bronchioalveolar cultures based on cell makers. Each color represents a cluster or different cell type. Number indicates number of cells in cluster identified as cell type presented.

t‐SNE map showing n class="Species">clustering of (L) with annotation of cell types. Cell type annotation is transferred from (A), and published dataset GSM2855485 is annotated as “Purified AT2”. Number indicates number of cells in cluster identified as cell type presented.

Differentiated bronchioalveolar cultures express SARS‐CoV‐2 eene">ntry receptor n class="Gene">ACE2 (green) on ATII‐like cells marked by HTII‐280 (red). Nuclei are stained with Hoechst and shown in white. Scale bar indicates 50 μm.

Bronchioalveolar cultures co‐express both entry receptors ACE2 (greeene">n) aene">nd n class="Gene">TMPRSS2 (red). Nuclei are stained with Hoechst and shown in white. Scale bar indicates 50 μm.

Replication competence and tropism of SARS‐CoV‐2 in 2D differentiated bronchioalveolar‐like cultures

In order to also assess the replication competence and tropism of SARS‐CoV‐2 in broene">nchioalveolar‐like cells, we n class="Disease">infected these cultures with SARS‐CoV‐2. SARS‐CoV‐2 readily infected the cells at a low multiplicity of infection (MOI) of 0.01 and reached slightly higher titers at a MOI of 0.1, as shown by titration of VeroE6 cells (Fig 5A and B) and viral RNA quantification (Fig EV4, EV5A and B). Shedding of virus occurred predominantly apically (Fig EV4, EV5C and D). In these cultures, we observed infection of alveolar cells (Fig 5C–E) as evidenced by staining against HOPX, HTII‐280, and SFTPC. Most infected cells appeared to be SFTPC+ (Fig 5E). Using a marker for dividing cells (Ki67), we showed that infected cells were not dividing (Appendix Fig S6C). We did not detect infected PNEC in the top compartment, nor infected mesenchymal cells in the basal compartment.

Figure 5

SARS‐CoV‐2 infects human organoid‐derived 2D bronchioalveolar‐like cultures

Infectious virus titers can be observed by virus titrations on VeroE6 cells of apical washes at 2, 24, 48, and 72 h after infection of bronchioalveolar‐like cells from donor 1 (A) and donor 2 (B) at MOI 0.01 (red) or 0.1 (black) with SARS‐CoV‐2. The dotted line indicates the lower limit of detection. Error bars represent SEM. N = 4 for all except donor 2 MOI 0.1, which is N = 3. H p.i. = hours post‐infection.

Immunofluorescent staining of SARS‐CoV‐2‐infected differentiated bronchioalveolar‐like cultures. Nucleoprotein (NP) stains viral capsid (red), which colocalized with HOPX (C), HTII‐280 (D), and SFTPC (E). Scale bars indicate 50 μm. Nuclei are stained with Hoechst and shown in white (C‐E).

Infectious virus titers can be observed by virus titratioene">ns oene">n VeroE6 cells of apical washes at 2, 24, 48, aene">nd 72 h after n class="Disease">infection of bronchioalveolar‐like cells from donor 1 (A) and donor 2 (B) at MOI 0.01 (red) or 0.1 (black) with SARS‐CoV‐2. The dotted line indicates the lower limit of detection. Error bars represent SEM. N = 4 for all except donor 2 MOI 0.1, which is N = 3. H p.i. = hours post‐infection.

Immunofluorescent staining of SARS‐CoV‐2‐n class="Disease">infected differentiated bronchioalveolar‐like cultures. Nucleoprotein (NP) stains viral capsid (red), which colocalized with HOPX (C), HTII‐280 (D), and SFTPC (E). Scale bars indicate 50 μm. Nuclei are stained with Hoechst and shown in white (C‐E).

Transmission electron microscopy analysis of SARS‐CoV‐2 infected bronchioalveolar‐like cells and small airway epithelial cells

Next, we performed traene">nsmissioene">n electroene">n microscopy (TEM) to coene">nfirm SARS‐n class="Species">CoV‐2 infection inside the 2D bronchioalveolar‐like (Fig EV4, EV5) and 2D small airway (Fig EV5G–K) cultures at high magnification. Upon inspection of the 2D bronchioalveolar‐like culture, ATII‐L were readily classified (Fig EV5A and B) based on the presence of lamellar bodies (Fig EV5D). SARS‐CoV‐2 infection was found in ATII‐L cells (Fig EV5C and D) and showed early infection onset with few double membrane vesicles and both electron dense (Fig EV5D) and lucent (Fig EV5C) compartments harboring virus particles. Interestingly, next to the ATII‐L cell, a fragmented dead cell is visible (Fig EV5F) with many extracellular virus particles still attached to the apical side and only few particles present on the basal side. This suggests that virus shedding mainly occurs on the air‐faced side inside lung tissue. Nonetheless, also in between the ATII‐L cells and an adjacent cell, some virus particles were detected, although it is not clear whether these sites are also the place of virion secretion or whether virions can diffuse between the cell layers. The latter is implausible due to the lack of detectable infectious virus in the basal compartment of the culture. Similarly, also the 2D small airway cultures (Fig EV5G–K) were examined with TEM and contained a confluent layer consisting mainly of ciliated cells (Fig EV5H–J) along with club cells and goblet cells (Fig EV5K). In accordance with confocal microscopy findings, at 72 h post‐infection, morphological evidence of SARS‐CoV‐2 infection, replication, and assembly was found inside ciliated cells (Fig EV5I). Secreted virions were again found at the apical side next to the still intact cilia (Fig EV5J), whereas no viruses were found at the basal side. Using confocal microscopy, we also detected rare infection of club cells (Fig 1D), but these rare events could not be visualized by TEM thus far.

Host responses induced by SARS‐CoV‐2 in bronchioalveolar‐like cells and small airway epithelial cells

To investigate the cellular response to infection, we stained broene">nchioalveolar‐like aene">nd small airway cultures for phosphorylated n class="Gene">STAT1 (pSTAT1) at 72 h post‐infection. STAT1 is phosphorylated in response to interferon (IFN) receptor signaling, which promotes translocation of STAT1 to the nucleus where it stimulates the expression of interferon‐stimulated genes (ISGs). We observed STAT1 phosphorylation and nuclear translocation in bystander cells of SARS‐CoV‐2 infected cultures in both bronchioalveolar‐like (Fig 6A) and small airway (Fig 6C) cultures, but not in infected cells, indicating potentially active blocking of the JAK/STAT pathway in infected cells. Active blocking of STAT1 nuclear import has been reported for SARS‐CoV. SARS‐CoV ORF6 antagonizes STAT1 function by sequestering nuclear import factors on the rough endoplasmic reticulum/Golgi membrane (Frieman et al, 2007; Kopecky‐Bromberg et al, 2007). Whether SARS‐CoV‐2 ORF6 has similar functions remains to be investigated, but recently the SARS‐CoV‐2 N protein was shown to suppress the translocation of STAT1 and STAT2 (Mu et al, 2020).

Figure 6

Host response analysis of SARS‐CoV‐2 infection in organoid‐derived bronchioalveolar‐like and small airway cells

pSTAT1 translocation in mock and SARS‐CoV‐2‐infected bronchioalveolar‐like cells visualized in green. Scale bars indicate 50 μm. Nucleoprotein (NP) stains viral capsid (red). Nuclei are stained with Hoechst (white). Scale bar indicates 50 µm.

Heatmaps depicting top 20 enriched genes based on adjusted P‐value in bronchioalveolar‐like cells after SARS‐CoV‐2 infection. Colored bars represent Z‐scores of log2‐transformed values.

pSTAT1 translocation in mock and SARS‐CoV‐2‐infected small airway cells visualized in green. Nucleoprotein (NP) stains viral capsid (red). Nuclei are stained with Hoechst (white). Scale bar indicates 50 μm.

Heatmaps depicting top 25 enriched genes based on adjusted p‐value and log2 fold changes in small airway cells after SARS‐CoV‐2 infection. Colored bars represent Z‐scores of log2‐transformed values.

GO term enrichment analysis for biological processes of the 20 most significantly up‐regulated genes upon SARS‐CoV‐2 infected in bronchioalveolar‐like cells. P‐values were determined using PANTHER Overrepresentation Test.

GO term enrichment analysis for biological processes of the 25 most significantly up‐regulated genes upon SARS‐CoV‐2 infected in small airway cells. P‐values were determined using PANTHER Overrepresentation Test.

Heatmaps depicting top 20 enriched genes based on adjusted P‐value in bronchioalveolar‐like cells after SARS‐CoV‐2 n class="Disease">infection. Colored bars represent Z‐scores of log2‐transformed values.

Heatmaps depicting top 25 enriched genes based on adjusted p‐value and log2 fold changes in small airway cells after SARS‐CoV‐2 n class="Disease">infection. Colored bars represent Z‐scores of log2‐transformed values.

GO term enrichment analysis for biological processes of the 20 most significantly up‐regulated genes upon SARS‐CoV‐2 n class="Disease">infected in bronchioalveolar‐like cells. P‐values were determined using PANTHER Overrepresentation Test.

GO term enrichment analysis for biological processes of the 25 most significantly up‐regulated genes upon SARS‐CoV‐2 n class="Disease">infected in small airway cells. P‐values were determined using PANTHER Overrepresentation Test.

Next, we performed mRn class="Gene">NA sequencing analysis to determine gene expression changes induced by SARS‐CoV‐2‐infection of bronchioalveolar‐like cultures. As suggested by pSTAT1 translocation, infection with SARS‐CoV‐2 elicited a broad signature of interferon‐stimulated genes (ISGs) attributed to type I/III interferon responses (Fig 6B and D), which was confirmed by Gene Ontology analysis (Fig 6E and F) in both bronchioalveolar‐like (Fig 6B and E, Table EV4) and small airway cultures (Fig 6D and F, Table EV5). This induction was similar to earlier reported SARS‐CoV‐2 infection in human bronchial epithelial cells (Blanco‐Melo et al, 2020) and indicated induction of strong chemotactic and inflammatory responses. Generally, the responses to SARS‐CoV‐2 appeared similar between the two pulmonary systems.

To test the suitability of the bronchioalveolar model for SARS‐CoV‐2 drug screeene">ns, we aene">nalyzed the effect of IFn class="Gene">N‐λ1 on viral replication. For this experiment, cultures were treated with IFN‐λ1 (50 ng/ml) 2 h prior to infection at an MOI of 0.1 or 24 h after infection. In the untreated control, SARS‐CoV‐2 replicated to high titers, approaching 106 TCID50 or 108 viral copies per ml (Fig 7A and B). Pre‐treatment with IFN‐λ1 abrogated viral replication entirely and treatment at 24 h post‐infection reduced infectious virus and RNA copies by ~ 5 logs and ~ 3 logs, respectively. These differences were also apparent by confocal imaging at 72 h post‐infection (Fig 7C). This bronchioalveolar culture system can therefore be applied to test COVID‐19 therapeutics.

Figure 7

Low‐dose interferon lambda 1 treatment dramatically reduces SARS‐CoV‐2 replication in bronchioalveolar‐like cultures

SARS‐CoV‐2 replication shown as infectious virus titers (A) and viral RNA copies (B) in apical washes at 2, 24, 48, and 72 h after infection at MOI 0.1. Cultures were either pretreated with 50 ng/ml IFN lambda 1 in the basal compartment 2 h prior to the infection (green), at 24 h p.i. or left untreated. The dotted line indicates the lower limit of detection. Error bars represent SEM. N = 2. H p.i. = hours post‐infection.

Whole well immunofluorescent staining of representative cultures in (A) and (B) using an anti‐NP antibody (red). Nuclei are stained with Hoechst (blue). Scale bar indicates 1 mm.

SARS‐CoV‐2 replicatioene">n showene">n as infectious virus titers (A) aene">nd viral Rn class="Gene">NA copies (B) in apical washes at 2, 24, 48, and 72 h after infection at MOI 0.1. Cultures were either pretreated with 50 ng/ml IFN lambda 1 in the basal compartment 2 h prior to the infection (green), at 24 h p.i. or left untreated. The dotted line indicates the lower limit of detection. Error bars represent SEM. N = 2. H p.i. = hours post‐infection.

Whole well immunofluorescent staining of representative cultures in (A) and (B) using an anti‐NP aene">ntibody (red). n class="Gene">Nuclei are stained with Hoechst (blue). Scale bar indicates 1 mm.

Discussion

In COVID‐19 n class="Species">patients, alveolar injury can trigger a cascade of events that leads to ARDS, limiting the transport of oxygen into the blood. A plethora of models are currently being developed to study this. In vivo models can be very useful, but are often extremely difficult and expensive to establish and may not represent conditions in humans. Most in vitro models utilize cell lines or primary cells that are often limited by the availability of donor materials and show donor–donor variation. In addition, primary human alveolar cultures are poorly susceptible to SARS‐CoV‐2 infection (Hou et al, 2020; Hui et al, 2020). Self‐renewing organoid models that contain stem cells which can differentiate into relevant cell types offer an elegant solution to develop human models that can be standardized. While a robust self‐renewing organoid model exists for the airway epithelium, the alveolar epithelium has proven very difficult to establish. As the alveolar epithelium is thought to play a major role in the disease progression of COVID‐19, the lack of an in vitro model greatly limits our understanding of this disease, but also of other respiratory virus infections.

This study shows that SARS‐CoV‐2 efficieene">ntly replicates in a n class="Species">human bronchioalveolar‐like model, targeting ATII‐L cells. This study is in accordance with clinical findings that SARS‐CoV‐2 infects alveolar cells in COVID‐19 patients. However, the variable incidence and severity of lower lung disease, and recent findings that alveolar cells have low levels of ACE2 expression in health, indicates that alveolar cells are unlikely to be the first cells infected through microaerosol inhalation of virus particles (Hou et al, 2020). Ciliated cells in the airways seem to express ACE2 more constitutively, with increasing expression toward the proximal part of the airway (Hikmet et al, 2020; Hou et al, 2020). Using an organoid‐derived small airway model, we show that ciliated cells are indeed the main target cell in the small airways, as was observed recently for primary small airway cells (Hou et al, 2020). Therefore, ciliated cells could be infected prior to the alveoli, with inflammatory or antiviral signals possibly leading to increased ACE2 expression as this gene was recently shown to be induced by interferon (Ziegler et al, 2020). Treatment of the bronchioalveolar cultures with a low dose of interferon lambda 1 was highly effective in reducing viral replication and dissemination, underscoring the high interferon sensitivity of this virus and indicating that interferons are a viable treatment option. Organoid‐derived alveolar‐like and airway models are ideal models to study these aspects of COVID‐19.

This model has several limitations. In human fetal developmeene">nt, alveolar cells are derived from bipoteene">nt progeene">nitor cells, whereas in adults, ATI cells are derived from a rare self‐reene">newing populatioene">n of ATII cells (Barkauskas et al, 2013; n class="Gene">Nikolic et al, 2017). In adult mice, an additional cell type, the bronchioalveolar stem cell (BASC) (Kim et al, 2005), differentiates into both bronchiolar and alveolar cells, but in adult humans, such bipotent progenitors have not been identified yet. Fetal and adult alveolar cells may therefore be phenotypically distinct and may respond differentially to infections. In addition, mesenchymal cells improved alveolar development in 2D and future work should focus on identifying which mesenchymal cues drive this. On the other hand, the presence of mesenchymal cells could offer an interesting opportunity to study differential expression of genes involved in lung fibrosis as well in these cells.

This model also has several advantages. The bronchioalveolar‐like model presented here shows a robust increase in infectious virus titers over time of ~ 5 logs, whereas other recently developed adult‐derived 3D alveolar organoid models show a more limited increase in viral titers (~ 1–2 logs) (Katsura et al, 2020; Salahudeen et al, 2020; Youk et al, 2020). This allows for a larger window to observe antiviral effects of potential drugs. In addition, in our model, pulmonary epithelial cells are infected from aene">nd shed virus apically, while in most 3D orgaene">noid models, the apical side of the cells is facing the inside of the orgaene">noid. Iene">n the broene">nchioalveolar‐like model preseene">nted here, the apical side of the cells is exposed to the air, which makes it is easily accessible, physiologically relevaene">nt aene">nd highly suitable for (high‐throughput) virology studies. Similar results may be obtained with the receene">ntly developed iPSC‐derived 2D ATII‐like cultures (Huaene">ng et al, 2020). Iene">n additioene">n, the modular Traene">nswell setup allows the additioene">n other cell types. For example, eene">ndothelial cells would be of interest, as the luene">ngs of deceased n class="Species">COVID‐19 patients were recently found to show distinctive vascular features, including endothelialitis, widespread thrombosis with microangiopathy, and angiogenesis (Ackermann et al, 2020). In addition, immune cells could be added, which may be relevant as severe COVID‐19 symptoms and hospitalization often occur ~ 8–15 days after disease onset (Lauer et al, 2020) indicating that acute lung injury may be the result of an overreaction of the immune system, rather than a direct effect of virus replication. As this model is derived from self‐renewing organoids, it is also amenable to genetic manipulation using CRISPR/Cas systems, which allows the study of specific host factors. Intrinsically, the model is also suitable for pulmonary cell and developmental biology studies, especially for studies on the interaction between epithelial and mesenchymal cells without the need for physical interaction.

In conclusion, bronchioalveolar‐like cells derived from self‐renewing human fetal n class="Chemical">lung organoids give rise to a 2D human bronchioalveolar‐like model with air‐exposed apical side. In these cultures, alveolar type II‐like cells are permissive to SARS‐CoV‐2 infection, which induces a type I/III interferon host response in vitro. SARS‐CoV‐2 replication is abrogated by low‐dose interferon lambda 1 treatment, showing that this model system can be used for COVID‐19 drug screens.

Materials and Methods

Reagents and Tools table

Methods and Protocols

Viruses and cells

Vero E6 cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemeene">nted with 10% fetal calf serum (FCS), n class="Chemical">HEPES, sodium bicarbonate, penicillin (100 IU/ml), and streptomycin (100 IU/ml) at 37°C in a humidified CO2 incubator. SARS‐CoV‐2 (isolate BetaCoV/Munich/BavPat1/2020; European Virus Archive Global #026V‐03883; kindly provided by Dr. C. Drosten) was propagated on VeroE6 (ATCC® CRL 1586TM) cells in Opti‐ MEM I (1×) + GlutaMAX (Gibco), supplemented with penicillin (100 IU/ml) and streptomycin (100 IU/ml) at 37°C in a humidified CO2 incubator. The SARS‐CoV‐2 isolate was obtained from a clinical case in Germany, diagnosed after returning from China. Stocks were produced by infecting VeroE6 cells at a multiplicity of infection (MOI) of 0.01 and incubating the cells for 72 h. The culture supernatant was cleared by centrifugation and stored in aliquots at −80°C. Stock titers were determined by preparing 10‐fold serial dilutions in Opti‐MEM I (1×) + GlutaMAX. Aliquots of each dilution were added to monolayers of 2 × 104 Vero E6 cells in the same medium in a 96‐well plate. Twenty‐four replicates were performed per virus stock. Plates were incubated at 37°C for 5 days and then examined for cytopathic effect. The TCID50 was calculated according to the method of Spearman & Kärber. All work with infectious SARS‐CoV and SARS‐CoV‐2 was performed in a Class II Biosafety Cabinet under BSL‐3 conditions at Erasmus Medical Center.

SARS‐CoV‐2 infection

2D differentiated cultures were washed twice with 200 µl Advanced DMEM/F12 (Gibco) with n class="Chemical">Glutamax (Gibco), HEPES (Lonza), penicillin (100 IU/ml), and streptomycin (100 IU/ml) (Lonza) (AdDF+++) before inoculation from the apical side at a MOI of 0.01 or 0.1 (as indicated) in 200 μl AdDF+++ per well. The number of cells in each culture was calculated by trypsinizing and counting the cells from 1 well. Next, cultures were incubated at 37°C 5% CO2 for 2 h before washing three times in 200 μl AdDF+++. At the indicated timepoints, virus was collected from the cells by adding 200 μl AdDF+++ apically, incubating 10 min at 37°C 5% CO2, and storing the supernatant at −80°C. Where indicated, cultures were treated basally with 50 ng/ml IFN‐L1 (Peprotech) at the indicated time pre‐ or post‐infection. Prior to determining the viral titer, all samples were centrifuged at 2,000 g for 5 min. Infectious virus titers were determined using the Spearman & Kärber TCID50 method on VeroE6 cells. All work was performed in a Class II Biosafety Cabinet under BSL‐3 conditions at Erasmus Medical Center.

Fixed immunofluorescence microscopy

Organoids and Transwell inserts were fixed in 4% paraformaldehyde for 20 min, permeabilized in 70% n class="Chemical">ethanol or 0.1% Triton X‐100 in PBS, and blocked for 60 min in 10% normal goat serum or 3% bovine serum albumin (BSA) in PBS (blocking buffer). For organoids, 0.1% Triton X‐100 was added to the blocking buffer to increase antibody penetration. Cells were incubated with primary antibodies overnight at 4°C in blocking buffer, washed twice with PBS, incubated with corresponding secondary antibodies Alexa 488‐ and Alexa 594‐conjugated secondary antibodies (1:200; Invitrogen) in blocking buffer for 2 h at room temperature, washed two times with PBS, incubated with indicated additional stains (TO‐PRO3, phalloidin‐633 (Santa Cruz, sc‐363796), Hoechst, DAPI), washed twice with PBS, and mounted in Prolong Antifade (Invitrogen) mounting medium.

Viral nucleocapsid was stained with mouse‐aene">nti‐SARS‐n class="Species">CoV NP (40143‐MM05, 1:1,000, Sino Biological) or rabbit‐anti‐SARS‐CoV NP (40143‐T62, 1:1,000, Sino biological), alveolar type 2 cells with anti‐HTII‐280 (TB‐27AHT2‐280, 1:150, Terrace Biotech), anti‐SFTPC (FL‐197, 1:200, Santa Cruz), anti‐LPCAT1 (HPA012501, 1:200, Sigma), alveolar type 1 cells with anti‐HTI‐56 (TB‐29 AHTI‐56, 1:30, Terrace Biotech) and anti‐HOPX (sc‐30216, 1:200, Santa Cruz Biotechnology), goblet cells with anti‐MUC5AC (MA5‐12178, 1:100, Invitrogen), club cells with anti‐CC10 (sc‐390313 AF594, 1:100, Santa Cruz Biotechnology) or anti‐SCGB3A2 (AF3545, 1:100, R&D systems), ciliated cells with anti‐AcTUB (sc‐23950 AF488, 1:100, Santa Cruz Biotechnology), STAT1 translocation anti‐pSTAT1 (MA5‐15071, 1:100, Thermo Fisher), proliferating cells with anti‐Ki67 (MIB‐1, GAG2661‐2, 1:100, DAKO), tight junctions with anti‐ZO‐I (ZO1‐1A12, 1:100, Invitrogen), multipotent lung bud tip stem cells with anti‐SOX2 (sc‐365823 AF488, 1:100, Santa Cruz Biotechnology) and anti‐SOX9 (sc‐166505 AF594, 1:100, Santa Cruz Biotechnology), basal cells with anti‐TP63 (ab735, 1:500, Abcam), and pulmonary neuroendocrine cells with anti‐SYP (sc‐17750, 1:100, Santa Cruz Biotechnology). ACE2 was stained with anti‐ACE2 (HPA000288, 1:200, Atlas antibodies), and TMPRSS2 was stained with anti‐TMPRSS2 (sc‐515727, 1:200, Santa Cruz). Samples were imaged on a LSM700 confocal microscope using ZEN software (Zeiss). Representative images were acquired and shown as Z‐projections, single slices or XZ cross sections.

Procurement of human material and informed consent

Adult lung tissue was obtained from residual, tumor‐free, material obtained at luene">ng resectioene">n surgery for n class="Disease">lung cancer. The Medical Ethical Committee of the Erasmus MC Rotterdam granted permission for this study (METC 2012‐512). Human fetal lung tissue was obtained from legally terminated second trimester pregnancies (15–20 weeks) by the HIS‐Mouse Facility of Academic Medical Center (AMC; Amsterdam, The Netherlands), after written informed consent of the mother for the tissue’s use in research and with approval of the Medical Ethical Review Board of the AMC (MEC: 03/038). Study procedures were performed according to the Declaration of Helsinki, and in compliance with relevant Dutch laws and institutional guidelines. The tissues obtained were anonymized and non‐traceable to the donor. On request by the researchers, only gender and gestational age is provided.

Isolation, culture, and differentiation of human small airway stem cells

Isolation of human small airway stem cells was performed using a protocol adapted from n class="Disease">Sachs et al (2019). One adult small airway donor was used in this study. Distal regions of adult lungs (lung parenchyma) were cut into ~ 2 mm sections, washed in AdDF+++, and incubated with dispase (Corning; #354235) mixed with airway organoid (AO) medium (Sachs et al, 2019) at a 1:1 ratio for 1 h at 37°C. The digested tissue suspension was sequentially sheared using 10‐ and 5‐ml plastic and flamed glass Pasteur pipettes. After shearing, the suspension was strained over a 100‐μm filter. Two percent FCS was added to the strained suspension before centrifugation at 400 g. The pellet was resuspended in 10 ml AdDF+++ and centrifuged again at 400 g. In case of a visible red pellet, erythrocytes were lysed in 2 ml red blood cell lysis buffer (Roche‐11814389001) for 5 min at room temperature before the addition of 10 ml AdDF+++ and centrifugation at 400 g. Lung cell pellets were then resuspended in 200 μl of growth factor reduced Matrigel (Corning; #356231) and plated in ~ 30 μl droplets in a 48‐well tissue culture plate. Plates were placed at 37°C with 5% CO2 for 10 min to solidify the Matrigel droplets upon which 250 µl of media was added to each well. Plates were incubated under standard tissue culture conditions (37°C, 5% CO2).

Airway organoids were passaged by mechanical disruption at 1:2 to 1:4 ratios after 10–14 days of culture. Medium was aspin class="Species">rated and ice‐cold AdDF+++ added to each well. The Matrigel in each well was disrupted by pipetting and sucked into a P1000 pipette tip and transferred into a 15‐ml tube. Ice‐cold AdDF+++ was added up to 10 ml after which the sample was centrifuged at 100 g for 5 min. The medium was aspirated and replaced by 1 ml TrypLE express (Gibco). Next, organoids were incubated at 37°C for 5 min after which they were triturated using a flame‐polished glass Pasteur pipette. Cold AdDF+++ containing 2% FCS was added up to 10 ml, and the sample was centrifuged at 400 g for 3 min. Medium was aspirated, washed once in cold AdDF+++, and the pellet incubated on ice for 2 min before being resuspended in Matrigel and plated in ~ 30 μl droplets in a 48‐well tissue culture. Plates were placed at 37°C with 5% CO2 for 20 min to solidify the Matrigel droplets, upon which 250 μl AO medium was added per well. Medium was changed every 4 days.

To obtain differentiated organoid‐derived cultures, organoids were dissociated into single cells using TrypLE express. Cells were seeded on Transwell membranes (Corning) coated with rat tail collageene">n type I (Fisher Scieene">ntific). Single cells were seeded in AO growth medium:complete base medium (CBM; Stemcell Pene">neumacult‐ALI; #05001) at a 1:1 n class="Species">ratio. After 2–4 days, confluent monolayers were cultured at air–liquid interface in CBM. Medium was changed every 5 days. After 8 weeks of differentiation, cultures were used for infection experiments.

Isolation, culture, and differentiation of human fetal lung bud tips

Isolation of human fetal luene">ng epithelial bud tips was performed using a protocol adapted from Miller et al (2018). Fetal luene">ng doene">nors were 16–19 weeks post‐coene">nceptioene">n (pcw). Tissues were shipped in Hiberene">nate A (Gibco; #A1247501) oene">n ice before processing oene">n the day the specimeene">n was obtained. Distal regioene">ns of fetal luene">ngs were cut into ~ 2 mm sectioene">ns aene">nd incubated with dispase oene">n ice for 30–60 min. Tissue pieces were theene">n incubated in 100% FBS oene">n ice for 15 min, after which they were traene">nsferred to a 10% FBS solutioene">n in n class="Chemical">AdDF+++. Lung bud tips were separated from mesenchymal cells through repeated pipetting and washing in AdDF+++ until mesenchymal cells were no longer visible. From each donor, donor‐specific mesenchymal cells were isolated by plating single cells present in the supernatant from the first wash. Mesenchymal cells were maintained in AdDF+++ supplemented with 10% FCS. After washing, bud tips were then resuspended in 200 μl of growth factor reduced Matrigel and plated in ~ 30 μl droplets in a 48‐well tissue culture plate. Plates were placed at 37°C with 5% CO2 for 20 min to solidify the Matrigel droplets upon which 250 µl of fetal lung (FL) medium (Nikolic et al, 2017) was added to each well. Plates were incubated under standard tissue culture conditions (37°C, 5% CO2).

Fetal lung bud tip organoids were passaged by mechanical disruption at 1:2 to 1:4 ratios after 10–14 days of culture. Medium was aspin class="Species">rated and ice‐cold AdDF+++ added to each well. The Matrigel in each well was disrupted by pipetting and sucked into a P1000 pipette tip and transferred into a 15‐ml tube. Ice‐cold AdDF+++ was added up to 10 ml after which the sample was centrifuged at 100 g for 5 min. The medium was aspirated and replaced by 1 ml AdDF+++. Next, organoids were triturated using a flame‐polished glass Pasteur pipette. Cold AdDF+++ was added up to 10 ml, and the sample was centrifuged at 400 g for 3 min. Medium was aspirated and the pellet incubated on ice for 2 min before being resuspended in Matrigel and plated in ~ 30 μl droplets in a 48‐well tissue culture. Plates were placed at 37°C with 5% CO2 for 20 min to solidify the Matrigel droplets, upon which 250 μl FL medium was added per well. Medium was changed every 4 days.

To obtain differentiated organoid‐derived cultures, 105 donor‐specific fetal lung mesenchymal cells were seeded in 12‐well plates in AdDF+++ coene">ntaining 10% FCS. The next day, orgaene">noids were dissociated into small clumps using TrypLE express. Cells were seeded oene">n Traene">nswell membraene">nes (Corene">ning‐3460). Clumps were seeded in FL:BEpiCM (broene">nchial epithelial cell growth medium; Scieene">ncell; Cat. n class="Gene">No. 3211) at a 1:1 ratio. BEpiCM was made by mixing basal medium (Sciencell; Cat. No. 3211‐b) and DMEM at a 1:1 ratio and adding bronchial epithelial cell growth supplement (BEpiCGS, Sciencell; Cat. No. 3262), penicillin (100 IU/ml), streptomycin (100 IU/ml), and 0.001 μM retinoic acid (Sigma). ROCK inhibitor was added for the first 5 days to prevent apoptosis (Y‐27632; Sigma; Y0503; 10 µM). After 4–7 days, confluent monolayers were cultured at air–liquid interface in BEpiCM. Eighty percent of the medium was changed every 4–5 days. After 14 days of differentiation, cultures were used for infection experiments.

Determination of virus titers using qRT–PCR

Supernatant and organoid samples were thawed and centrifuged at 2,000 g for 5 min. Sixty μl supernatant was lysed in 90 μl MagnaPure LC Lysis buffer (Roche) at room temperature for 10 min. Rn class="Gene">NA was extracted by incubating samples with 50 μl Agencourt AMPure XP beads (Beckman Coulter) for 15 min at RT, washing beads twice with 70% ethanol on a DynaMag‐96 magnet (Invitrogen), and eluting in 30 μl MagnaPure LC elution buffer (Roche). Viral titers (Viral RNA copies per ml) were determined by qRT–PCR using primers targeting the E gene (35) and compared to a counted RNA standard curve.

Mass spectrometry

Proteins were digested with trypsin using the SP3 method (Hughes et al, 2019) with minor modifications. Briefly, proteins secreted by alveolar‐like cultures were collected from the cells by adding 200 μl AdDF+++ apically, incubating 10 min at 37°C 5% n class="Chemical">CO2, and storing the supernatant at −80°C in Laemmli buffer (Bio‐Rad). Proteins in 30 μl Laemmli buffer (Bio‐Rad) were reduced with 5 mM dithiothreitol (Sigma, D8255) for 30 min at 50°C and alkylated with 10 mM iodoacetamide (Sigma, I6125) for 10 min in the dark. Ten µg Sera‐Mag SpeedBeads (Fisher Scientific) in 20 µl milli‐Q and ethanol were added to a final concentration of 50%, and the solution was mixed for 10 min at RT. Beads were immobilized using a magnetic rack and washed three times with 100 µl 80% ethanol. One 1 µg trypsin and 0.5 µg Lys‐C in 100 µl 50 mM Tris/HCl pH 8.3 were added to the beads, and the sample was incubated overnight at 37°C. The tryptic digest was acidified with TFA and desalted using a StageTip. Peptides were eluted with 100 µl 40% Acetonitrile and 0.1% formic acid and dried in a Speedvac. Before analysis by LC‐MS, peptides were dissolved in 20 µl 2% acetonitrile and 0.1% formic acid.

Naene">noflow liquid chromatography–taene">ndem mass spectrometry (nLC‐MS/MS) was performed oene">n aene">n EASY‐nLC 1200 coupled to aene">n Orbitrap Eclipse Tribrid mass spectrometer (Thermo Fisher Scieene">ntific) open class="Species">rating in positive mode. Peptide mixtures were trapped on a 2 cm × 100 μm Pepmap C18 column (Thermo Fisher 164564) and then separated on an in‐house packed 50 cm × 75 μm capillary column with 1.9 μm Reprosil‐Pur C18 beads (Dr. Maisch) at a flowrate of 250 nL/min, using a linear gradient of 0–32% acetonitrile (in 0.1% formic acid) during 90 min. The eluate was directly sprayed into the electrospray ionization (ESI) source of the mass spectrometer. Spectra were acquired in continuum mode; fragmentation of the peptides was performed in data‐dependent mode by HCD. Mass spectrometry data were analyzed with Proteome Discoverer 2.4 with Mascot version 2.6.2 as the search engine. The target FDR for both PSMs and peptides was set to 0.01. MS/MS spectra were searched against a combined database composed of all UniProt entries taxonomy Homo Sapiens (release October 2019) and all SARS‐CoV‐2 protein sequences (https://covid‐19.uniprot.org/) concatenated with the reversed versions of all sequences. For precursor ion quantification based on the Minora algorithm in Proteome Discoverer, only unique and razor peptides were used. Normalization was done on total peptide amounts for all samples.

Bulk RNA sequencing

Library preparatioene">n was performed at Single Cell Disn class="Species">coveries (Utrecht, The Netherlands), using an adapted version of the CEL‐seq protocol. In brief:

Total RNA was extracted using the staene">ndard n class="Chemical">TRIzol (Invitrogen) protocol and used for library preparation and sequencing. mRNA was processed as described previously, following an adapted version of the single‐cell mRNA seq protocol of CEL‐Seq (Simmini et al, 2014; Hashimshony et al, 2016). In brief, samples were barcoded with CEL‐seq primers during a reverse transcription and pooled after second strand synthesis. The resulting cDNA was amplified with an overnight in vitro transcription reaction. From this amplified RNA, sequencing libraries were prepared with Illumina TruSeq small RNA primers. Paired‐end sequencing was performed on the Illumina NextSeq 500 platform using barcoded 1 × 75 nt read setup. Read 1 was used to identify the Illumina library index and CEL‐Seq sample barcode. Read 2 was aligned to the CRCh38 human RefSeq transcriptome, with the addition of the SARS‐CoV‐2 (Ref‐SKU: 026V‐03883) genome, similar to earlier study (Lamers et al, 2020), using BWA using standard settings (Li & Durbin, 2010). Reads that mapped equally well to multiple locations were discarded. Mapping and generation of count tables was done using the MapAndGo script. Samples were normalized using RPM normalization.

Differential gene expression analysis was performed using the DESeq2 package (Love et al, 2014). SARS‐mapping reads were removed before analyzing the different datasets.

For bronchioalveolar‐like samples, two donors were sequenced twice as technical replicates which were combined in the DESeq2 methods when comparing to the remaining datasets. We compared the bulk mRNA sequeene">ncing of the broene">nchioalveolar culture system to freshly purified adult ATII cells. Previously published profiles of purified n class="Species">human ATII cultures were downloaded from NCBI GEO (Accession Nos. GSM3818009 and GSM3818010). Normalized counts were used in generating heatmaps.

The top genes which were most strongly induced in response to SARS‐CoV‐2 (Padj < 0.05, raene">nked by fold chaene">nge) were subjected to fuene">nctioene">nal eene">nrichmeene">nt aene">nalysis for a biological process using ToppFuene">n oene">n the ToppGeene">ne Suite (https://toppgeene">ne.cchmc.org/eene">nrichmeene">nt.jsp) as described before (Kaimal et al, 2010). The five biological processes with highest eene">nrichmeene">nt (after FDR correctioene">n aene">nd a P‐value cutoff of 0.05) for each virus are displayed with the correspoene">nding GO term, corrected P‐value, aene">nd FDR.

10× Single‐cell sequencing—library preparation

Single‐cell mRNA sequeene">ncing was performed at Single Cell Disn class="Species">coveries according to standard 10× Genomics 3’ V3.1 chemistry protocol. In short, single‐cell suspensions were methanol fixed in 80% methanol and frozen at −80°C. Prior to loading the cells on the 10× Chromium controller, cells were rehydrated in rehydration buffer and counted to assess cell concentration. Cells were loaded, and the resulting sequencing libraries were prepared following standard 10× Genomics protocol. The DNA libraries were paired‐end sequenced on an Illumina NovaSeq S4, with a 2 × 150 bp Illumina kit.

10× Single‐cell sequencing—Mapping of sequencing data and clustering and differential gene expression analysis

BCL files resulting from sequencing were transformed to FASTQ files with 10× Genomics Cell Ranger mkfastq. FASTQ files were mapped with Cell Ranger count. During sequencing, Read 1 was assigned 150 base pairs, but only the first 28 bp were used for identification of the Illumina library barcode, cell barcode, and UMI. R2 was used to map to a custom made reference transcriptome of Hg38 with additioene">nal sequeene">nces for SARS‐n class="Species">CoV‐2 added. Filtering of empty barcodes was done Cell Ranger. In total, approximately 5,015 cells were detected. Unsupervised clustering and differential gene expression analysis was done with the Seurat1 R toolkit (Butler et al, 2018). We compared the ATII‐L cells in the bronchioalveolar culture system to purified ATII cells. Previously published single‐cell transcriptomic profiles of purified ATII cultures were downloaded from NCBI GEO (Accession No GSM2855485). Raw counts were normalized, integrated with the abovementioned single‐cell dataset, and further analyzed using the Seurat R toolkit.

Transmission electron microscopy

2D cultures were chemically fixed for 3 h at room temperature aene">nd prepared for imaging within the traene">nsmissioene">n electroene">n microscope as described in (Lamers et al, 2020). All TEM data were collected with aene">n Eagle detector oene">n a FEI Tecnai T12 microscope open class="Species">rated at 120 kV.

Statistical analysis

For single‐cell sequencing data, clustering aene">nd differeene">ntial geene">ne expressioene">n aene">nalysis was performed with Seun class="Species">rat version 3 (Butler et al, 2018), following standard parameters.

Author contributions

Conceptualization, MML, JvdV, JB, HC, BLH; Methodology, MML, JVDV, KK, SR, JAAD, MJM, PJP, HC, BLH; Investigation, MML, JvdV, KK, SR, TIB, AZM, DS, KB, CDK, NG, RBGR, HQD, JAAD; Resources, GMGMV, MPGK, HC; Writing‐Original Draft, MML, JVDV, BLH; Writing‐Review & Editing, PJP, HC, BLH; Fuene">nding acquisitioene">n, PJP, HC, BLH.

Conflict of interest

H.C. is inventor on patents held by the Royal Netherlaene">nds Academy of Arts aene">nd Scieene">nces that n class="Species">cover organoid technology. H.C.’s full disclosure is given at https://www.uu.nl/staff/JCClevers/.

Appendix

Click here for additional data file.

Expanded View Figures PDF

Click here for additional data file.

Table EV1

Click here for additional data file.

Table EV2

Click here for additional data file.

Table EV3

Click here for additional data file.

Table EV4

Click here for additional data file.

Table EV5

Click here for additional data file.

Review Process File

Click here for additional data file.

Reagent/Resource	Reference or Source	Identifier or Catalog Number
Experimental Models
SARS‐CoV‐2 BavPat1	C. Drosten	European Virus Archive Global #026V‐03883
Small airway organoids	This study
Lung bud tip organoids Donor 1 (#19162)	This study
Lung bud tip organoids Donor 2 (#2)	This study
VeroE6 cells	ATCC	CRL 1586TM
Recombinant DNA
N/A
Antibodies
Mouse‐anti‐SARS‐CoV NP (monoclonal)	Sino Biological	40143‐MM05
Rabbit‐anti‐SARS‐CoV NP (polyclonal)	Sino Biological	40143‐T62
Mouse‐anti‐HTII‐280 (monoclonal)	Terrace Biotech	TB‐27AHT2‐280
Rabbit‐anti‐SFTPC (polyclonal)	Santa Cruz	FL‐197
Rabbit‐anti‐HOPX (polyclonal)	Santa Cruz	sc‐30216
Mouse‐anti‐HTI‐56 (monoclonal)	Terrace Biotech	TB‐29 AHTI‐56
Mouse‐anti‐MUC5AC (monoclonal)	Invitrogen	MA5‐12178
Anti‐LPCAT1 (polyclonal)	Sigma	HPA012501
Mouse‐anti‐CC10 (monoclonal)	Santa Cruz	sc‐390313 AF594
Goat‐anti‐SCGB3A2 (polyclonal)	R&D systems	AF3545
Mouse‐anti‐AcTUB (monoclonal)	Santa Cruz	sc‐23950 AF488
Rabbit‐anti‐pSTAT1 (monoclonal)	Thermo Fisher	MA5‐15071
Mous‐anti‐Ki67 (monoclonal)	DAKO	GAG2661‐2
Mouse‐anti‐ZO‐I (monoclonal)	Invitrogen	ZO1‐1A12
Mouse‐anti‐SOX2 (monoclonal)	Santa Cruz	sc‐365823 AF488
Mouse‐anti‐SOX9 (monoclonal)	Santa Cruz	sc‐166505 AF594
Mouse‐anti‐TP63 (monoclonal)	Abcam	ab735
Mouse‐anti‐SYP	Santa Cruz	sc‐17750
Rabbit‐anti‐ACE2	Atlas antibodies	HPA000288
Mouse‐anti‐TMPRSS2	Santa Cruz	sc‐515727
Goat‐anti‐mouse IgM 594	Invitrogen	A21044
Goat‐anti‐mouse IgG 488	Invitrogen	A11029
Goat‐anti‐mouse IgG2a 488	Invitrogen	A21131
Goat‐anti‐rabbit IgG 488	Invitrogen	A32731
Rabbit‐anti‐goat IgG 594	Invitrogen	A11080
Rabbit‐anti‐goat IgG 488	Invitrogen	A11078
Goat‐anti‐mouse IgG1 594	Invitrogen	A21125
Goat‐anti‐rabbit IgG 594	Invitrogen	A11012
Goat‐anti‐mouse IgM 488	Invitrogen	A21042
Oligonucleotides and other sequence‐based reagents
SARS‐CoV‐2 E gene taqman	Corman et al (2020)
CEL‐seq2‐based SORT‐seq primers	Hashimshony et al (2016), van den Brink et al (2017)
Illumina Truseq small RNA primers	Illumina
10x genomics primers and reagents V3.1	10× Genomics
Chemicals, enzymes and other reagents
Advanced DMEM/F12	Gibco	12634010
Glutamax	Gibco	35050061
AO medium	Sachs et al (2019)
FL medium	Nikolic et al (2017)
Opti‐MEM I (1×) + GlutaMAX	Gibco	51985‐042
HEPES	Lonza	BE17‐737E
Penicillin‐Streptomycin Mixture (Pen/Strep)	Lonza	DE17‐602E
Hoechst	Thermo Fisher	H1399
DAPI	Sigma	D9542
TO‐PRO‐3	Thermo Fisher	T3605
Phalloidin	Santa Cruz	ec‐363796
Prolong Diamond Antifade	Invitrogen	P36961
red blood cell lysis buffer	Roche	11814389001
Matrigel (GFR)	Corning	#356231
Dispase	Corning	#354235
TrypLE express	Gibco	12604‐013
Stemcell Pneumacult‐ALI	Stemcell	#05001
Hibernate A	Gibco	#A1247501
BEpiCM	Sciencell	#3211
Basal medium	Sciencell	#3211‐b
BEpiCGS	Sciencell	#3262
Retinoic acid	Sigma	R2625
ROCK inhibitor	Sigma	Y0503
MagnaPure LC Lysis buffer	Roche	05323738001
Agencourt AMPure XP beads	Beckman Coulter	A63882
MagnaPure LC elution buffer	Roche	05323738001
TaqMan™ Fast Virus 1‐Step Master Mix	Applied Biosystems	4444436
Laemmli buffer	Bio Rad	#161‐0747
Dithiothreitol	Sigma	D8255
iodoacetamide	Sigma	I6125
Sera‐Mag SpeedBeads	Fisher Scientific	09‐981‐123
TRIzol	Thermo Fisher	15596026
Pepmap C18 column	Thermo Fisher	164564
IFN‐L1	Peprotech	300‐02L
Software
ZEN software	ZEISS
Prism 8	Graphpad
Cell Ranger 4.0.0	10× Genomics
Rstudio v 1.1.463	Rstudio
Seurat v 3	Satija lab, NYU Genome Centre
DynaMag‐96	Invitrogen	12331D
LSM700	Zeiss
Transwell inserts	Corning	3260 & 3270
Orbitrap Eclipse Tribrid mass spectrometer	Thermo Fisher
Illumina Nextseq500	Illumina
Tecnai T12 microscope	FEI
10× Genomics chromium controller	10× Genomics