Human airway cells prevent SARS-CoV-2 multibasic cleavage site cell culture adaptation

  1. Mart M Lamers
  2. Anna Z Mykytyn
  3. Tim I Breugem
  4. Yiquan Wang
  5. Douglas C Wu
  6. Samra Riesebosch
  7. Petra B van den Doel
  8. Debby Schipper
  9. Theo Bestebroer
  10. Nicholas C Wu
  11. Bart L Haagmans  Is a corresponding author
  1. Viroscience Department, Erasmus Medical Center, Netherlands
  2. Department of Biochemistry, University of Illinois at Urbana-Champaign, United States
  3. Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, United States
  4. Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, United States
8 figures, 1 table and 1 additional file

Figures

Figure 1 with 1 supplement
SARS-CoV-2 rapidly acquires multibasic cleavage site mutations when propagated on VeroE6 cells.

(A–C) Deep-sequencing analysis of VeroE6 passage 2 (A), passage 3 (B), and passage 4 (C) virus stocks. In each graph, the amino acid sequence logo of the multibasic cleavage site is shown. (D–F) Sanger sequencing chromatograms of VeroE6 passage 2 (D), passage 3 (E), and passage 4 (F) viruses. Multibasic cleavage site mutations identified by deep-sequencing are indicated with arrows. Translated sequences are indicated below Sanger reads. (G) Plaque size analysis of VeroE6 passage 2–4 virus stocks on VeroE6 cells. Red arrow heads indicate small plaques. Scale bar indicates 1 cm.

Figure 1—figure supplement 1
Deep-sequencing analysis of VeroE6 passage 1 virus multibasic cleavage site and full genome deep-sequencing analysis of passage 1–4 viruses.

(A) Deep-sequencing analysis of the VeroE6 passage 1 virus stock. In each graph, the amino acid sequence logo of the multibasic cleavage site is shown. (B) Full genome deep-sequencing analysis of VeroE6 passage 1, 2, 3, and passage 4 viruses.

Mutations in the multibasic cleavage site and the adjacent serine residue (S686) abrogate S1/S2 cleavage.

(A) Analysis of S1/S2 cleavage by S1 immunoblot of SARS-CoV-2 S (WT), multibasic cleavage site (MBCS) mutant and S686G mutant pseudoviruses. (B) Quantification of S1 cleavage from four independent pseudovirus productions. (C) Analysis of S1/S2 cleavage by multiplex S1 (red) and S2 (green) immunoblot of SARS-CoV-2 S (WT) and S686G mutant pseudoviruses. S0 indicates uncleaved spike; S1 indicates the S1 domain of cleaved spike; VSV-N indicates VSV nucleoprotein (production control). Numbers indicate the molecular weight (kDa) of bands of the protein standard. (D) Quantification of S2 cleavage from four independent pseudovirus productions. Error bars indicate SD. EV = empty vector. WT = wild type. kDa = kilo dalton.

The SARS-CoV-2 multibasic cleavage site and the adjacent serine residue (S686) enhance infectivity and serine protease mediated entry on Calu-3 and VeroE6-TMPRSS2 cells.

(A–B) SARS-CoV-2 (WT), multibasic cleavage site (MBCS) mutant and S686G pseudovirus infectious titers on (A) VeroE6 and (B) Calu-3 cells. (C) Fold change in SARS-CoV-2, MBCS mutant and S686G pseudovirus infectious titers on Calu-3 cells over infectious titers on VeroE6 cells. (D) SARS-CoV-2, MBCS mutant and S686G pseudovirus infectious titers on VeroE6-TMPRSS2 cells. (E) Fold change in SARS-CoV-2, MBCS mutant and S686G pseudovirus infectious titers on VeroE6-TMPRSS2 cells over infectious titers on VeroE6 cells. One-way ANOVA was performed for statistical analysis comparing all groups with WT. (F–I) SARS-CoV-2, MBCS mutant and S686G pseudovirus entry into (F and G) VeroE6 cells or (H and I) VeroE6-TMPRSS2 cells pre-treated with a concentration range of either (F and H) camostat mesylate or (G and I) E64D. Two-way ANOVA, followed by a bonferroni post hoc test was performed for statistical analysis comparing all groups to WT. WT pseudovirus entry into VeroE6 cells treated with 10 µM E64D was significantly different from del-RRAR, R682A, R685A and S686G pseudovirus entry. * indicates statistical significance (p<0.05) compared to WT (A–E). * indicates statistical significance (p<0.05) compared to WT at the highest inhibitor concentration (F–I). Experiments were performed in triplicate. Representative experiments from at least two independent experiments are shown. Error bars indicate SD. WT = wild type.

Multibasic cleavage site mutations and the adjacent serine residue (S686) impair spike protein fusogenicity.

(A–C) Fusogenicity of wild type SARS-CoV-2 spike and spike mutants was assessed after 18 hr by measuring the sum of all GFP+ pixels per well in a GFP-complementation fusion assay on VeroE6-GFP1-10 (A), VeroE6-TMPRSS2-GFP1-10 (B), and Calu-3-GFP1-10 (C) cells. The experiment was performed in triplicate. A representative experiment from two independent experiments is shown. Statistical analysis was performed by one-way ANOVA. * indicates a significant difference compared to WT (p<0.05). Error bars indicate SD. EV = empty vector. WT = wild type.

Figure 5 with 1 supplement
SARS-CoV-2 propagation in Calu-3 cells efficiently prevents SARS-CoV-2 cell culture adaptation.

(A) Deep-sequencing analysis of Calu-3 passage 2 virus from a VeroE6 passage 1. (B) Deep-sequencing analysis of Calu-3 passage 3 virus from the Calu-3 passage 2 in A. (C) Deep-sequencing analysis of Calu-3 passage 3 virus grown from a VeroE6 passage 2 stock (Figure 1A). Deep-sequencing analysis of Calu-3 passage 5 virus from a Calu-3 passage 3 stock in C. In each graph, the amino acid sequence logo of the multibasic cleavage site is shown.

Figure 5—figure supplement 1
Multibasic cleavage site deep-sequencing analysis of passage 4 Calu-3 viruses from an adapted VeroE6 P3 stock and full genome deep-sequencing analysis of Calu-3-propagated viruses.

(A) Deep-sequencing analysis of Calu-3 passage 4 virus from a VeroE6 passage 3 stock (from Figure 1B). (B) Deep-sequencing analysis of Calu-3 passage 4 virus from a VeroE6 passage 3 stock produced in the presence of 10 μM E64D. In each graph, the amino acid sequence logo of the multibasic cleavage site is shown. (C) Full genome deep-sequencing analysis of Calu-3-propagated viruses.

Figure 6 with 1 supplement
Serine protease expression prevents MBCS mutations.

(A–B) Deep-sequencing analysis of VeroE6 passage 4 virus from a VeroE6 passage 3 (A is a redisplay of Figure 1C) mock-treated or treated with 10 μM camostat. (C–D) Deep-sequencing analysis of VeroE6-TMPRSS2 passage 4 virus from a VeroE6 passage 3 mock-treated or treated with 10 μM camostat. In each graph the amino acid sequence logo of the multibasic cleavage site is shown.

Figure 6—figure supplement 1
Multibasic cleavage site and full genome deep-sequencing analysis of passage 4 VeroE6 and VeroE6-TMPRSS2 viruses.

(A–B) Deep-sequencing analysis of VeroE6 passage 4 virus from a VeroE6 passage 3, trypsin-treated (A) or treated with 10% FBS (B). (C–D) Deep-sequencing analysis of VeroE6-TMPRSS2 passage 4 virus from a VeroE6 passage 3, trypsin-treated (C) or treated with 10% FBS (D). In each graph, the amino acid sequence logo of the multibasic cleavage site is shown. (E) Full genome deep-sequencing analysis of VeroE6 and VeroE6-TMPRSS2-propagated viruses. In E VeroE6 P4 (mock) is a redisplay of VeroE6 P4 in Figure 1—figure supplement 1B.

A 2D air-liquid interface human airway organoid model for SARS-CoV-2 propagation.

(A) Human airway organoids were dissociated and plated onto 12 mm transwell inserts. After an 8–12 week differentiation period at air-liquid interface cultures contained ciliated, non-ciliated and basal cells as shown on a hematoxylin-eosin stain. (B) Air-exposed cells, but not basal cells, expressed the priming protease TMPRSS2 as shown by immunohistochemistry. (C) Immunofluorescent staining indicated that in these cultures, ciliated cells (acetylated tubulin+ or AcTUB+ cells) were infected by SARS-CoV-2. (D and E) At 5 days post-infection, whole-well confocal imaging indicated the infection was widespread (D) and cytopathic effects, including cilia damage (D and E) and syncytial cells (E) were visible. Scale bars indicate 20 µm in A, B, C; 2 mm in D; and 100 µm in E.

Figure 8 with 1 supplement
2D air-liquid interface human airway organoids produce high titer stocks without multibasic cleavage site mutations.

(A–B) Deep-sequencing analysis (A) and Sanger chromatogram (B) of Organoid passage 3 virus from a VeroE6 passage 2 stock (Figure 1A). The amino acid sequence logo of the multibasic cleavage site is shown. The translated sequence is indicated below the Sanger read. Arrows indicate where cell culture adaptations to VeroE6 cells occur. (C) Plaque size analysis of VeroE6 passage 2 and Organoid passage 3 virus (the VeroE6 data is a redisplay of Figure 1G). Red arrow heads indicate large plaques. Scale bar indicates 1 cm. (D) Full genome deep-sequencing analysis of VeroE6 passage 2 and organoid passage three stocks. In D VeroE6 P2 is a redisplay of VeroE6 P2 in Figure 1—figure supplement 1B. (E) Immunoblot analysis of VeroE6 passage 2 and 3, Calu-3 passage 3 and Organoid passage 3 stocks. S0 indicates uncleaved spike; S1 indicates the S1 domain of cleaved spike; NP indicates nucleoprotein. Numbers indicate the molecular weight (kDa) of bands of the protein standard. (F) Quantification of cleavage from three immunoblots. Error bars indicate SD.

Figure 8—figure supplement 1
Schematic workflow for the production of SARS-CoV-2 stocks on 2D air-liquid interface differentiated airway organoids.

Step 1. 3D self-renewing airway organoids are grown from human lung tissue. Next, these are dissociated to single cells and differentiated at air-liquid interface for 4–12 weeks. Step 2. Differentiated cultures are infected at a multiplicity of infection of 0.05 and washed daily for 5 days. The washes from day 2 to day 5 are collected and stored at 4°C. Step 3. Virus collections are cleared by centrifugation and filtered to remove debris larger than 0.45 μm. Next, the medium is exchanged three times using Amicon columns to remove cytokines and debris smaller than 100 kDa. Purified virus preparations are then stored at −80°C in aliquots. Step 4. Stocks can be characterized using plaque assays, Sanger sequencing and deep-sequencing. Created with BioRender.com.

Tables

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional information
AntibodyRabbit-anti-SARS-CoV NP (polyclonal)Sino BiologicalCat# 40143-T62IF (1:1000)
AntibodyMouse anti-TMPRSS2 (monoclonal)Santa CruzCat# sc-515727IHC (1:200)
AntibodyGoat-anti-mouseDakoCat# P0260IF (1:400)
AntibodyGoat anti-rabbit IgG (H+L) Alexa Fluor Plus 594InvitrogenCat# A32740IF (1:400)
AntibodyGoat anti-mouse IgG (H+L) Alexa Fluor 488InvitrogenCat# A11029IF (1:2000)
AntibodyMouse-anti-AcTub IgG2A Alexa Fluor 488 (monoclonal)Santa Cruz BiotechnologyCat# sc-23950 AF488IF (1:100)
AntibodyMouse anti-nucleocapsidSinobiologicalCat# 40143-MM05IF (1:1000)
AntibodyRabbit anti-SARS-CoV S1 (polyclonal)SinobiologicalCat# 40150-T62WB (1:1000)
AntibodyMouse-anti-SARS-CoV-2 S2 (monoclonal)GenetexCat# GTX632604WB (1:1000)
AntibodyMouse-anti-VSV-N (monoclonal)Absolute AntibodyCat# Ab01403-2.0WB (1:1000)
Biological sample (Homo sapiens)Airway organoidsMykytyn et al., 2021
Cell line (Cercopithecus aethiops)VeroE6ATCCCRL 1586TM
Cell line (Cercopithecus aethiops)VeroE6 TMPRSS2Mykytyn et al., 2021
Cell line (Cercopithecus aethiops)VeroE6 GFP1-10Mykytyn et al., 2021
Cell line (Cercopithecus aethiops)VeroE6 GFP1-10 TMPRSS2Mykytyn et al., 2021
Cell line (Homo sapiens)Calu-3ATCCHTB 55
Cell line (Homo sapiens)Calu-3 GFP1-10Mykytyn et al., 2021
Chemical compound, drugE64DMedChemExpressCat# HY-100229
Chemical compound, drugCamostat mesylateSigmaCat# SML0057
Chemical compound, drugPolyethylenimine linearPolysciencesCat# 23966
Chemical compound, drugHygromycin BInvitrogenCat# 10843555001
Chemical compound, drugGeneticinInvitrogenCat# 10131035
Chemical compound, drugAvicelFMC biopolymers-
Chemical compound, drugLaemmliBioRadCat# 1610747
Commercial assay or kitSuperScript
IV Reverse Transcriptase
InvitrogenCat# 18090200
Commercial assay or kitPfu Ultra II Fusion HS DNA PolymeraseAgilent TechnologiesCat# 600674
Commercial assay or kitQiaquick PCR Purification KitQIAGENCat# 28104
Commercial assay or kitBigDye
Terminator v3.1 Cycle Sequencing Kit
Applied BiosystemsCat# 4337456
Commercial assay or kitProtoScript II Reverse TranscriptaseNew England BioLabsCat# NEB M0368X
Commercial assay or kitKAPA HyperPlusRocheCat# 7962428001
OtherAmicon Ultra-15 Centrifugal Filter Unit with Ultracel-100 membraneMilliporeCat#
UFC910024
OtherOpti-MEM I (1X) + GlutaMAXGibcoCat# 51985–042
OtherAdvanced DMEM/F12Thermo Fisher scientificCat# 12634–010
OtherAO mediumSachs et al., 2019N/A
OtherPneumacult ALI mediumStemcellCat # 05001
OtherTryplEThermo Fisher scientificCat# 12605010
OtherBasement membrane extractR and D SystemsCat# 3533-005-02
OtherTranswell insertsCorningCat# 3460
OtherCollagen Type I, High concentration Rat tailCorningCat# 354249
Other0.45 μm low protein binding filterMilliporeCat# SLHV033RS
OtherHoechstThermo FisherCat# H1399
OtherAmpure XP BeadsBeckman CoulterCat# A63882
OtherIllumina sequencer V3 MiSeq flowcellIllumina
OtherABI PRISM 3100 Genetic AnalyzerApplied Biosystems
OtherOdyssey CLxLicor
OtherAmersham Typhoon Biomolecular ImageGE Healthcare
OtherAmersham Imager 600GE Healthcare
OtherLSM700 confocal microscopeZeiss
OtherCarl ZEISS Vert.A1Zeiss
Software, algorithmZENZeiss
Software, algorithmImageQuant TL 8.2GE Healthcare
Software, algorithmStudio Lite Ver 5.2Licor
Software, algorithmGraphPad PRISM 8, 9GraphPad
Software, algorithmIllustratorAdobe inc
Strain, strain background (SARS-CoV-2)SARS-CoV-2 BavPat-1Dr. Christian DrostenEuropean Virus Archive Global #026 V-03883

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Mart M Lamers
  2. Anna Z Mykytyn
  3. Tim I Breugem
  4. Yiquan Wang
  5. Douglas C Wu
  6. Samra Riesebosch
  7. Petra B van den Doel
  8. Debby Schipper
  9. Theo Bestebroer
  10. Nicholas C Wu
  11. Bart L Haagmans
(2021)
Human airway cells prevent SARS-CoV-2 multibasic cleavage site cell culture adaptation
eLife 10:e66815.
https://doi.org/10.7554/eLife.66815