Human liver organoid derived intra-hepatic bile duct cells support SARS-CoV-2 infection and replication and its comparison with SARS-CoV

Liver tissues

Pediatric liver tissues were obtained from non-tumour margin of hepatoblastoma (HB). Liver biopsies were obtained during operations with full informed consent from parents or patients, and the study was approved by Hong Kong West Cluster-Hong Kong University Cluster Research Ethics Committee/Institutional Review Board (UW 16-052).

Human pediatric liver organoid culture

Wedge liver biopsies (1-2 cm³) from the non-tumor margin of pediatric patients (n=4) with hepatoblastoma (HB)). Liver tissues were minced in cold wash medium (Advanced DMEM/F12; 1% GlutaMAX; 1% FBS; 1% Penicillin/Streptomycin (P/S)) and digested in digestion medium (5 ml, Multi Tissue Dissociation Kit 1; Miltenyi Biotec Inc. CA, USA) in a gentleMACS-C Tube (Miltenyi Biotec Inc. CA, USA) on gentleMACS™ Octo Dissociator (Miltenyi Biotec Inc. CA, USA) using program 37°C-Multi-A-01. After digestion and washing with cold wash medium (5 ml), the digested sample was passed through a 70 µm strainer and then through a 30 µm strainer before centrifugation (300 g; 10 min) to pellet the cells. Cell pellet was resuspended cells in 200 µl of MACS column buffer (PBS pH 7.2 + 0.5 % BSA), and incubated with CD326 (EPCAM) microbeads (1 µl microbeads/2×10⁶ cells; 130-090-500; Miltenyi Biotec Inc. CA, USA) in the dark at 4°C for 30 minutes. After incubation, column buffer (0.5 ml) was added and EpCAM positive cells were sorted on MS Column following the manufacturer’s protocol (Miltenyi Biotec Inc. CA, USA). The cells were pelleted (300 g; 10 min.) and resuspended in 200 µl of column buffer for cell counting. The cell pellet (approximately 1×10⁵ cells) was mixed with 50-60 µl of Matrigel (356231; Corning Biocoat) and seeded per well of a prewarmed (37°C) 4-well plate (176740; Nunc; Thermo Scientific). After Matrigel had solidified, organoid medium (500 µl) was added to each well. Organoid medium was based on Advanced DMEM/F12 (Invitrogen) supplemented with Penicillin/Streptomycin (Invitrogen), GlutaMax (Invitrogen), 25 mM HEPES (Invitrogen), 1% N2 (GIBCO), 1% B27 (GIBCO), 1.25 mM Acetylcysteine (Sigma), 10 nM gastrin (G9145; Sigma), 50 ng/ml EGF (PMG8043; PeproTech), 100 ng/ml FGF10 (100-26-25UG; PeproTech), 25 ng/ml HGF (100-39-10UG; PeproTech), 10 mM Nicotinamide (Sigma), 5 µM A83.01 (Tocris), 10 µM Forskolin (Tocris), 500 ng/ml R-Spondin 1 (7150-RS-025; R&D), 100 ng/ml Noggin (250-38-20UG; Peprotech), 100 ng/ml Wnt3a (1324-WN-010; R&D) and 250 ng/ml Amphotericin B (#15290018; GIBCO). For the first six days of culture, 10 µM ROCK inhibitor Y-27632 (Tocris) was added. Medium was changed once every two days.

Single cell RNA Sequencing and data analysis

Single cells from organoids were prepared. The organoid medium along with matrix gel containing organoids was transferred to a 15 ml tube, 1-2 ml cold Advanced DMEM/F12 was added and the mixture was incubated on ice for 10 min to dissolve matrix gel before centrifugation (300g; 5 min). The supernatant was aspirated until the organoid pellet and a layer of matrix gel remained. 1 ml 5X TrypLE Express (12604013; GIBCO) was added, mixed well and incubated at 37°C for 5 min. 1 ml FBS was added to the mixture, which was pipetted up and down for 40-50 times with a small circumference opening glass pipette (diameter 0.3-0.5 mm) to efficiently dissociate organoids into single cells. 5-10 ml cold Advanced DMEM/F12 medium was added, and the suspension was passed through a 30 µm strainer before centrifugation (300 g; 5 min) at 4°C. The supernatant was aspirated until only the pellet remained and the cells were counted following addition of 200 µl 10% FBS. The cells were then counted using Invitrogen Countess II FL Automated Cell Counter. Samples were prepared as outlined by the 10x Genomics Single Cell 3’ v2 Reagent Kit user guide. The samples with more than 1 x 10⁶ cells/ml along with >80% viability is then resuspended in 50-100ul, 10% FBS and are submitted for 10X Genomics single-cell sequencing at Centre for PanorOmic Sciences (CPOS)-HKU. Sample viability was double confirmed by Trypan Blue (Thermo Fisher) and using a haemocytometer (Thermo Fisher). Following counting, the appropriate volume, ∼33.8ul of each sample were used for cell encapsulation for a target capture of 4000 cells. Sequencing run were carried out using Illumina NovaSeq 6000 for those sequencing libraries passing the QC with good cDNA yield. Single cell transcriptomic data have been uploaded to NCBI Sequence Read Archive (SRA) and can be found with accession number: PRJNA609259.

The CellRanger (10X Genomics) analysis pipeline was used to generate digital gene expression matrix (UMI counts per gene per cell) from sequencing data by aligning to the human genome. The raw gene expression matrix was filtered, normalized and clustered using standard Seurat package procedures 16. The low quality cells were removed from the analysis using the following thresholds: cells with a very small library size or UMI counts per cell (nUMI >1500), genes detected per cell (nGene >1000), UMIs vs. genes detected (log10GenesPerUMI >0.8), mitochondrial counts ratio (mitoRatio <0.1). Cell-cycle phases were predicted using a function included in Seurat that scores each cell based on expression of canonical marker genes for S and G2/M phases, this was used to regress out the cell-cycle effects from the downstream analysis. Clusters were visualized using uniform manifold approximation and projection (UMAP) as implemented in Seurat. First 13 PCs were selected for UMAP based on the points where the principal components cumulatively contribute 90% of variation associated with entire data. The cell-type identities for each cluster in UMAP were determined using presence of known marker genes in each cluster.

Virus stock preparation

A SARS-CoV-2 (BetaCoV/Hong Kong/VM20001061/2020, SCoV2) was isolated from a COVID-19 patient in Hong Kong in 2020. For comparison, we used a SARS-CoV (strain HK39849, SCoV) isolated from a hospitalized SARS patient in Hong Kong in 2003. SARS-CoV-2 and SARS-CoV viruses were propagated in Vero-E6 cells and virus stock was titrated to determine tissue culture infection dose 50% (TCID₅₀) in Vero-E6 cells. The experiments were carried out in a Bio-safety level 3 (BSL-3) facility.

Coronavirus infection

The 3D liver organoids were sheared mechanically using syringe to expose the apical surface to the virus inoculum. Around 100-200 organoids were infected with each coronavirus at 5×10⁵ TCID₅₀/ml for 1 h at 37°C. The organoids were washed three times with culture medium, re-embedded in Matrigel at the same conditions with the same growth medium and incubated at 37°C with 5% CO₂. The viral titers in the culture supernatants were measured at 1, 24, 48, and 72 h after infection using the TCID₅₀ assay in Vero-E6 cells. Cell lysates were collected 72h after infection to assess mRNA expression of cytokines. Organoids were fixed 72 h after infection in paraformaldehyde for immunofluorescent staining.

Immunofluorescence

Organoids in matrigel (from 2 HB livers) were fixed in 4% paraformaldehyde (w/v) in PBS (phosphate-buffered saline, pH 7.2) for 48 h at 4°C, dehydrated in graded series of alcohol, and cleared in xylene before being embedded in paraffin. Sections (6 µm in thickness) were prepared and mounted onto TESPA-coated microscope glass. Sections were dewaxed in xylene, hydrated in a graded series of alcohol and finally in distilled water. Antigen retrieval was performed by incubation of slides in Citrate buffer (pH 6.0) at 95°C for 10 minutes. After blocking in PBS-T (PBS with 0.1% Triton) supplemented with 1% Bovine Serum Albumin for 1 h at room temperature, sections were incubated with antibody diluted in PBS-T/BSA for overnight at 4°C. After washing in PBS-T, sections were incubated with fluorescent tagged secondary antibodies in PBS-T/BSA at 37°C for 1 h. Details of the primary and secondary antibodies and their dilutions are shown in Table 1 (supplementary data). After PBS-T washings, sections were mounted in DAPI-containing anti-fade mounting fluid. Images were taken with Nikon Eclipse 80i microscope mounted with a SPOT RT3 microscope digital camera under fluorescence illumination. Photos were compiled using Adobe Photoshop CS6.

Quantitative RT-PCR

The RNA of infected cells was extracted at 72 h post infection using a MiniBEST universal RNA extraction kit (Takara Biotechnology). RNA was reverse-transcribed by using oligo-dT primers with RT-PCR kit (Takara). mRNA expression of target genes was performed using an ABI ViiA™7 real-time PCR system (Applied Biosystems). The gene expression profiles of cytokines and chemokines were quantified and normalized with β-actin as previously described [26].

Statistical analysis

Statistical analysis was done using GraphPad Prism software version 9. Experiments with the human organoids were performed independently in two different donors each with triplicate wells. Viral titers and area under the curve (AUC) derived from viral titers and mRNA expression were compared using one- or two-way ANOVA with Tukey’s multiple comparisons test. Mock infected organoids served as negative controls. Results shown in figures are the calculated mean and standard deviation of mean. Differences were considered significant at p < 0.05.

Expression of ACE2 and TMPRSS2 in human liver organoids

To investigate if human liver derived organoids can be used to establish an ex vivo SARS-CoV-2 infection model for intra-hepatic bile ducts, we first determined whether the organoid culture could preserve the cholangiocytes expressing ACE2 and TMPRSS2 ex vivo. We performed single-cell RNA sequencing (scRNA-seq) analysis of human liver organoids from 2 patients to interrogate the transcriptomic signatures of cells in human liver organoids. The organoids used in the study are derived from the hepatoblast (intrahepatic ducts) progenitor (CK19+) lineage, and thus recapitulate the intrahepatic cholangiocyte development. A total number of 5885 cells were analyzed and cell populations were visualized by uniform manifold approximation and projection (UMAP), partitioning the cells into three clusters (Figure 1). The common cholangiocyte markers epithelial cell adhesion molecule (EPCAM) and keratin 19 (CK19) were uniformly highly expressed in all the clusters, indicating the heterogeneity of cholangiocytes in these organoids was relatively low (Figure 1). Notably, we identified the SARS-CoV-2 receptor gene ACE2 expressed sparsely among all the cluster in unbiased preferences and was detectable in 7.32% cells of developing cholangiocytes, 4.78% cells in cholangiocytes and 17.48% cells in hepatoblasts (Figure 1B and 1C). Besides, TMPRSS2 expressed uniformly across all the clusters and accounted for 30.18% of developing cholangiocytes, 20.52% in cholangiocytes and 17.89% in hepatoblasts, it is worth mentioning that the ACE2+ cells were co-expressing TMPRSS2 (159 out of 233) (Figure 1B and 1C), making this cell population potentially highly vulnerable to SARS-CoV-2 infection. Interestingly, these same cells were also found to express another serine protease, TMPRSS4. Immunostaining further verified the presence of ACE2+, TMPRSS2+ cholangiocytes (CK19+) in human liver ductal organoids (Figure 1D). Taken together, our data demonstrate that human liver tissue derived organoid preserves the human-specific ACE2+/TMPRSS2+ population of intra-hepatic cholangiocytes.

• Download figure
• Open in new tab

Figure 1

(A) UMAP showing the assigned identity for each cluster (3 clusters identified) in 10X Genomics single-cell RNA sequencing analysis of human control liver organoids. (B & C) Feature plot and violin plot displaying single cell expression distributions for features (ACE2, TMPRSS2 and TMPRSS4) in each cluster. (E) Immunofluorescence staining of human liver organoids for ACE2, TMPRSS2 and TMPRSS4.

Intra-hepatic cholangiocytes are susceptible to SARS-CoV-2 infection

Next, we examined the susceptibility of human liver ductal organoids to SARS-CoV-2. We infected liver organoids with the SARS-CoV-2, and the SARS-CoV. The liver organoids from two individuals were inoculated with viruses for 1 h, re-embedded in Matrigel and cultured for 48 and 72 hrs. We performed immunostaining to identify the virus-positive cholangiocytes 48 h and 72 h post-infection (Figure 2 and data not shown). The expression of SARS-CoV nucleocapsid protein (Sco-Np) was detected in patchy areas of human liver organoids infected with SARS-CoV-2 and SARS-CoV, whereas no signal was found in uninfected control (Mock, Figure 2A).

• Download figure
• Open in new tab

Figure 2

Viral infection and replication kinetics of SARS-CoV-2 and SARS-CoV viruses in human liver organoids. (A) Immunofluorescence staining for SARS-CoV nucleocapsid protein (Sco-Np) of mock and viral infected organoids at 48 h post-infection. (B) Bar chart showing viral titres in culture supernatants of human liver organoids infected with SARS-CoV-2 and SARS-CoV viruses. The organoids were infected with 5×10⁵ TCID₅₀/mL at 37°C. Culture supernatants were harvested at the indicted times and virus titres were measured by TCID₅₀ assay. Bar-charts show the mean virus titre ± standard deviation (SD). The horizontal dotted line denotes the limit of detection in the TCID₅₀ assay. Statistical significances compared among viruses were calculated using two-way ANOVA with Tukey’s multiple comparisons test. ***p<0·001. (C) Bar-charts showing the area under the curve (AUC) derived from the viral titres from 24 h to 72 h. Data are presented as mean ± standard deviation (SD). Statistical significances compared among viruses were calculated using one-way ANOVA with Tukey’s multiple comparisons test. ***p<0·001.

Intra-hepatic cholangiocytes supports robust replication of SARS-CoV-2

Virus tropism and replication competence of SARS-CoV-2 with SARS-CoV were compared. Here, we showed that the viral titres of SARS-CoV-2 increased starting at 24h and with more than 2 log increase from 24h to 72h in bile duct organoids. SARS-CoV-2 replicated significantly higher than that of the SARS-CoV, with about 1 log increase in viral titer along at all-time points (Figure 2B). When comparing the area under the curve (AUC) derived from the viral titers from 24h to 72h, similar observations were found, in that SARS-CoV-2 had a significantly higher viral replication competence than SARS-CoV (Figure 2C). These data demonstrate that human bile duct organoids are susceptible to SARS-CoV-2 infection and support robust viral replication.