Abstract
Free full text
SARS-CoV-2 hijacks neutralizing dimeric IgA for nasal infection and injury in Syrian hamsters1
ABSTRACT
Prevention of robust severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection in nasal turbinate (NT) requires in vivo evaluation of IgA neutralizing antibodies. Here, we report the efficacy of receptor binding domain (RBD)-specific monomeric B8-mIgA1 and B8-mIgA2, and dimeric B8-dIgA1, B8-dIgA2 and TH335-dIgA1 against intranasal SARS-CoV-2 challenge in Syrian hamsters. These antibodies exhibited comparable neutralization potency against authentic virus by competing with human angiotensin converting enzyme-2 (ACE2) receptor for RBD binding. While reducing viral loads in lungs significantly, prophylactic intranasal B8-dIgA unexpectedly led to high amount of infectious viruses and extended damage in NT compared to controls. Mechanistically, B8-dIgA failed to inhibit SARS-CoV-2 cell-to-cell transmission, but was hijacked by the virus through dendritic cell-mediated trans-infection of NT epithelia leading to robust nasal infection. Cryo-EM further revealed B8 as a class II antibody binding trimeric RBDs in 3-up or 2-up/1-down conformation. Neutralizing dIgA, therefore, may engage an unexpected mode of SARS-CoV-2 nasal infection and injury.
Introduction
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a member of the Betacoronavirus genus, is the causative agent of Coronavirus Disease 2019 (COVID-19) [1]. SARS-CoV-2 enters host cells through the binding of the receptor binding domain (RBD) of its surface trimeric spike (S) protein to the cellular angiotensin-converting enzyme-2 (ACE-2) receptor [2]. After the ACE2 binding, the S protein is cleaved into S1 and S2 subunits by host cellular proteases including the transmembrane protease serine 2 (TMPRSS2) to promote fusion of viral and cellular membranes for viral entry [3]. Apparently, these entry processes are similar to those of SARS-CoV, the causative agent of SARS, although these two coronaviruses share just 76% and 40% amino acid identity between their genomes and their RBD external subdomains, respectively [4]. SARS-CoV-2, however, has displayed remarkably higher transmissibility than SARS-CoV [5]. By the end of 2022, the rapid all-year-round transmission of SARS-CoV-2 has resulted in several waves of pandemic with over 667 million infections and 6.7 million deaths [6]. This is in great contrast to the SARS epidemic, which caused only 8096 cases and 774 deaths, and disappeared in 2003 [7]. Therefore, prevention of robust SARS-CoV-2 nasal infection remains a research priority.
The high transmissibility of SARS-CoV-2 is likely associated with multiple factors. First, unlike SARS-CoV-infected cases characterized by high fever and prominent respiratory symptoms, afebrile individuals with SARS-CoV-2 were often found upon diagnosis, allowing person-to-person transmissions by asymptomatic carriers including international travellers [8]. Second, while both coronaviruses employ the same receptor ACE2, highly conserved RBD residues or side chain properties of SARS-CoV-2 might account for increased ACE2 binding [9,10]. Third, the unique insertion of the PRRA sequence in SARS-CoV-2 S glycoprotein promotes higher virion infectivity and cell-to-cell transmission, leading to enhanced pathogenicity in vivo [11]. Fourth, by altering the S protein conformation, the D614G mutation increases the stability of S trimer to avoid premature S1 shedding, which results in a rapid dominance of this mutation globally [12]. The D614G mutation also induces higher infectious titres in nasal washes and the trachea of infected hamsters [13]. Subsequent S mutations in variants of concerns (VOCs) further improved ACE2 affinity for enhanced infectivity [14,15]. Multiple mechanisms, therefore, contribute to the high upper respiratory tract (URT) viral loads (robust nasal infection) characteristic of SARS-CoV-2 infection [16,17], which facilitates viral transmission in the human population. So far, the in vivo activities of various types of IgA remains understudied.
Since the outbreak of COVID-19, worldwide research efforts have led to the identification of many potent human neutralizing antibodies (HuNAbs) mainly in IgG form for preclinical and clinical developments [18]. Some studies also investigated IgA antibodies, which are known to play an important role in mucosal immunity, especially in their secretory form (SIgA) [19,20]. RBD-specific IgA antibodies were rapidly discovered in COVID-19 patients [21]. IgA and SIgA were even shown to dominate the early antibody response as compared to IgG and IgM in saliva and bronchoalveolar lavage fluids due to expansion of IgA plasmablasts with mucosal homing characteristics [22]. However, it was noted that COVID-19 patients with acute respiratory distress syndrome (ARDS) had higher SIgA in the airway mucosa for unknown reasons [23]. Moreover, a recent study reported that SARS-CoV-2 viral loads were closely associated with spike-specific IgA responses in the nasal samples of acute COVID-19 patients [24]. Since dIgA were shown to be about 15 times more potent than mIgA in vitro against the same target, dIgA were suggested to be particularly valuable for therapeutic application against SARS-CoV-2 [25]. It is, therefore, critical to investigate the potential of mIgA and dIgA in preventing SARS-CoV-2 nasal infection in vivo.
In this study, we used the technology of single B cell antibody gene cloning to generate a panel of SARS-CoV-2 RBD-specific monoclonal HuNAbs from the peripheral blood mononuclear cells (PBMCs) of one acute and three convalescent COVID-19 patients in Hong Kong. Since intramuscular or intranasal inoculation of several potent IgG HuNAbs cannot completely prevent SARS-CoV-2 infection in the nasal turbinate (NT) of Syrian hamsters [26], we sought to improve the efficacy of HuNAb by converting IgG to IgA. To achieve this goal, we engineered the potent B8-IgG1 into monomeric IgA1 (B8-mIgA1), monomeric IgA2 (B8-mIgA2), dimeric IgA1 (B8-dIgA1) and dimeric IgA2 (B8-dIgA2), as well as a human native dimeric IgA1 (TH335-dIgA1). We then determined their efficacies in the Syrian hamster model against the live intranasal ancestral SARS-CoV-2 challenge as we previously described [27].
Materials and methods
Human subjects
A total of 4 patients with COVID-19 including 3 convalescent cases and one acute case were recruited between February and May 2020. All patients were confirmed by reverse-transcription polymerase chain reaction (RT–PCR) as described previously [16]. Clinical and laboratory findings were entered into a predesigned database. Written informed consent was obtained from all patients. This study was approved by the Institutional Review Board of The University of Hong Kong/Hospital Authority Hong Kong West Cluster, the Hong Kong East Cluster Research Ethics Committee, and the Kowloon West Cluster Research Ethics Committee (UW 13-265, HKECREC-2018-068, KW/EX-20-038[144-26]).
Syrian hamsters
The animal experimental plan was approved by the Committee on the Use of Live Animals in Teaching and Research (CULATR 5359-20) of the University of Hong Kong (HKU). Male and female golden Syrian hamsters (Mesocricetus auratus) (aged 6–10 weeks) were purchased from the Chinese University of Hong Kong Laboratory Animal Service Centre through the HKU Laboratory Animal Unit (LAU). The animals were kept in Biosafety Level-2 housing and given access to standard pellet feed and water ad libitum following LAU’s standard operational procedures (SOPs). The viral challenge experiments were then conducted in our Biosafety Level-3 animal facility following SOPs strictly, with strict adherence to SOPs.
Cell lines
HEK293 T cells, HEK293T-hACE2 cells Vero-E6 cells, HK2 cells and Vero-E6-TMPRSS2 cells were maintained in DMEM containing 10% FBS, 2 mM L-glutamine, 100 U/ml penicillin and incubated at 37°C in a 5% CO2 setting [28]. Expi293FTM cells were cultured in Expi293TM Expression Medium (Thermo Fisher Scientific) at 37°C in an incubator with 80% relative humidity and a 5% CO2 setting on an orbital shaker platform at 125±5 rpm/min (New Brunswick innova™ 2100) according to the manufacturer’s instructions.
ELISA analysis of plasma and antibody binding to RBD and trimeric spike
The recombinant RBD and trimeric spike proteins derived from SARS-CoV-2 (Sino Biological) were diluted to final concentrations of 1 μg/ml, then coated onto 96-well plates (Corning 3690) and incubated at 4 °C overnight. Plates were washed with PBS-T (PBS containing 0.05% Tween-20) and blocked with blocking buffer (PBS containing 5% skim milk or 1% BSA) at 37 °C for 1 h. Serially diluted plasma samples or isolated monoclonal antibodies were added to the plates and incubated at 37 °C for 1 h. Wells were then incubated with a secondary goat anti-human IgG labelled with horseradish peroxidase (HRP) (Invitrogen) or with a rabbit polyclonal anti-human IgA alpha-chain labelled with HRP (Abcam) and TMB substrate (SIGMA). Optical density (OD) at 450 nm was measured by a spectrophotometer. Serially diluted plasma from healthy individuals or previously published monoclonal antibodies against HIV-1 (VRC01) were used as negative controls.
Isolation of RBD-specific IgG+single memory B cells by FACS
RBD-specific single B cells were sorted as previously described [29]. In brief, PBMCs from infected individuals were collected and incubated with an antibody cocktail and a His-tagged RBD protein for identification of RBD-specific B cells. The cocktail consisted of the Zombie viability dye (Biolegend), CD19-Percp-Cy5.5, CD3-Pacific Blue, CD14-Pacific Blue, CD56-Pacific Blue, IgM-Pacific Blue, IgD-Pacific Blue, IgG-PE, CD27-PE-Cy7 (BD Biosciences) and the recombinant RBD-His described above. Two consecutive staining steps were conducted: the first one used an antibody and RBD cocktail incubation of 30 min at 4 °C; the second staining involved staining with anti-His-APC and anti-His-FITC antibodies (Abcam) at 4 °C for 30 min to detect the His tag of the RBD. The stained cells were washed and resuspended in PBS containing 2% FBS before being strained through a 70-μm cell mesh filter (BD Biosciences). RBD-specific single B cells were gated as CD19+CD27+CD3-CD14-CD56-IgM-IgD-IgG+RBD+ and sorted into 96-well PCR plates containing 10 μL of RNAase-inhibiting RT–PCR catch buffer (1M Tris-HCl pH 8.0, RNase inhibitor, DEPC-treated water). Plates were then snap-frozen on dry ice and stored at −80 °C until the reverse transcription reaction.
Single B cell RT–PCR and antibody cloning
Single memory B cells isolated from PBMCs of infected patients were cloned as previously described [30]. Briefly, one-step RT–PCR was performed on sorted single memory B cell with a gene specific primer mix, followed by nested PCR amplifications and sequencing using the heavy chain and light chain specific primers. Cloning PCR was then performed using heavy chain and light chain specific primers containing specific restriction enzyme cutting sites (heavy chain, 5′-AgeI/3′-SalI; kappa chain, 5′-AgeI/3′-BsiWI). The PCR products were purified and cloned into the backbone of antibody expression vectors containing the constant regions of human Igγ1 or Igα1 and Igα2. The Igα1 and Igα2 vectors were purchased from InvivoGen (pfusess-hcha1 for IgA1 and pfusess-hcha2m1 for IgA2). The constructed plasmids containing paired heavy and light chain expression cassettes were co-transfected into 293 T cells (ATCC) grown in 6-well plates. Antigen-specific ELISA and pseudovirus-based neutralization assays were used to analyze the binding capacity to SARS-CoV-2 RBD and the neutralization capacity of transfected culture supernatants, respectively.
Genetic analysis of the BCR repertoire
Heavy chain and light chain germline assignment, framework region annotation, determination of somatic hypermutation (SHM) levels (in nucleotides) and CDR loop lengths (in amino acids) were performed with the aid of the IMGT/HighV-QUEST software tool suite (www.imgt.org/HighV-QUEST). Sequences were aligned using Clustal W in the BioEdit sequence analysis package (Version 7.2). Antibody clonotypes were defined as a set of sequences that share genetic V and J regions as well as an identical CDR3.
Antibody production and purification
The paired antibody VH/VL chains were cloned into Igγ, Igα1 or Igα2 and Ig
Size exclusion chromatography
The prepacked HiLoad 26/60 SuperdexTM 200pg (code No. 17-1071-01, Cytiva) column was installed onto the Amersham Biosciences AKTA FPLC system. After column equilibration with 2 column volumes (CV) of PBS, the concentrated IgA antibodies were applied onto the column using a 500-µl loop at a flow rate of 2 ml/min. Dimers of IgA1 or IgA2 were separated from monomers upon washing with 2 CV of PBS. The milli-absorbance unit at OD280 nm was recorded during the washing process. 2 ml-fractions were collected, pooled, concentrated and evaluated by western blot using mouse anti-IGJ monoclonal antibody [KT109] (Abcam) and rabbit anti-human IgA alpha chain antibody (Abcam).
Pseudovirus-based neutralization assay
The neutralizing activity of NAbs was determined using a pseudotype-based neutralization assay as we previously described [31]. Briefly, The pseudovirus was generated by co-transfection of 293 T cells with pVax-1-S-COVID19 and pNL4-3Luc_Env_Vpr, carrying the optimized spike (S) gene (QHR63250) and a human immunodeficiency virus type 1 backbone, respectively [31]. Viral supernatant was collected at 48 h post-transfection and frozen at −80 °C until use. The serially diluted monoclonal antibodies or sera were incubated with 200 TCID50 of pseudovirus at 37 °C for 1 hour. The antibody-virus mixtures were subsequently added to pre-seeded HEK 293T-ACE2 cells. 48 hours later, infected cells were lysed to measure luciferase activity using a commercial kit (Promega, Madison, WI). Half-maximal (IC50) or 90% (IC90) inhibitory concentrations of the evaluated antibody were determined by inhibitor vs. normalized response – 4 Variable slope using GraphPad Prism 6 or later (GraphPad Software Inc.).
Neutralization activity of monoclonal antibodies against authentic SARS-CoV-2
The SARS-CoV-2 focus reduction neutralization test (FRNT) was performed in a certified Biosafety level 3 laboratory. Neutralization assays against live SARS-CoV-2 were conducted using a clinical isolate (HKU-001a strain, GenBank accession no: MT230904.1) previously obtained from a nasopharyngeal swab from an infected patient [17]. Briefly, nasopharyngeal swab specimen was inoculated in serum-free MEM-cultured Vero-E6 cells. The viruses in a CPE positive sample were cloned by plaque purification. The plaque purified viral isolate was amplified by 1 additional passage in Vero-E6 cells to make working stocks of the virus. The tested antibodies were serially diluted, mixed with 50 μL of SARS-CoV-2 (1×103 focus forming unit/ml, FFU/ml) in 96-well plates, and incubated for 1 hour at 37°C. Mixtures were then transferred to 96-well plates pre-seeded with 1×104/well Vero E6 cells and incubated at 37°C for 24 hours. The culture medium was then removed and the plates were air-dried in a biosafety cabinet (BSC) for 20 mins. Cells were then fixed with a 4% paraformaldehyde solution for 30 min and air-dried in the BSC again. Cells were further permeabilized with 0.2% Triton X-100 and incubated with cross-reactive rabbit sera anti-SARS-CoV-2-N for 1 hour at RT before adding an Alexa Fluor 488 goat anti-rabbit IgG (H+L) cross-adsorbed secondary antibody (Life Technologies). The fluorescence density of SARS-CoV-2 infected cells were scanned using a Sapphire Biomolecular Imager (Azure Biosystems) and the neutralization effects were then quantified using Fiji software (NIH).
Antibody binding kinetics, and competition with the ACE2 receptor measured by Surface Plamon Resonance (SPR)
The binding kinetics and affinity of recombinant monoclonal antibodies for the SARS-CoV-2 spike protein (ACROBiosystems) were analysed by SPR (Biacore 8 K, GE Healthcare). Specifically, the spike protein was covalently immobilized to a CM5 sensor chip via amine groups in 10 mM sodium acetate buffer (pH 5.0) for a final RU around 500. SPR assays were run at a flow rate of 30 ml/min in HEPES buffer. For conventional kinetic/dose–response, serial dilutions of monoclonal antibodies were injected across the spike protein surface for 180s, followed by a 600s dissociation phase using a multi-cycle method. Remaining analytes were removed in the surface regeneration step with the injection of 10 mM glycine-HCl (pH 2.0) for 2×30s at a flow rate of 30 µl/min. Kinetic analysis of each reference subtracted injection series was performed using the Biacore Insight Evaluation Software (GE Healthcare). All sensorgram series were fit to a 1:1 (Langmuir) binding model of interaction. Before evaluating the competition between antibodies and the human ACE2 peptidase domain, both the saturating binding concentrations of antibodies and of the ACE2 protein (ACROBiosystems) for the immobilized SARS-CoV-2 spike protein were determined separately. In the competitive assay, antibodies at the saturating concentration were injected onto the chip with immobilized spike protein for 120s until binding steady-state was reached. ACE2 protein also used at the saturating concentration was then injected for 120s, followed by another 120s of injection of antibody to ensure a saturation of the binding reaction against the immobilized spike protein. The differences in response units between ACE2 injection alone and prior antibody incubation reflect the antibodies’ competitive ability against ACE2 binding to the spike protein.
Hamster experiments
In vivo evaluation of monoclonal antibody B8-IgG1, B8-mIgA1, B8-mIgA2, B8-dIgA1, B8-dIgA2 in the established golden Syrian hamster model of SARS-CoV-2 infection was performed as described previously, with slight modifications [27]. Approval was obtained from the University of Hong Kong (HKU) Committee on the Use of Live Animals in Teaching and Research. Briefly, 6–8-week-old male and female hamsters were housing with access to standard pellet feed and water ad libitum until live virus challenge in the BSL-3 animal facility. The hamsters were randomized from different litters into experimental groups. Experiments were performed in compliance with the relevant ethical regulations approved by the University of Hong Kong (CULATR 5359-20 and 5518-20). For prophylaxis studies, 24 hours before live virus challenge, three groups of hamsters were intraperitoneally (500 µl) or intranasally (100 µl through pipette) administered with one dose of test antibody in phosphate-buffered saline (PBS) at the indicated dose. At day 0, each hamster was intranasally inoculated with a challenge dose of 100 μL of Dulbecco’s Modified Eagle Medium containing 105 PFU of SARS-CoV-2 (HKU-001a strain, GenBank accession no: MT230904.1) under anaesthesia with intraperitoneal ketamine (200 mg/kg) and xylazine (10 mg/kg). For post-exposure treatment study, each hamster only received one 1.5 mg/kg dose of intraperitoneal B8-IgG1 at 24, 48, 72 hours (n=4 per group) after virus challenge. The hamsters were monitored twice daily for clinical signs of disease. Syrian hamsters typically clear virus within one week after SARS-CoV-2 infection. Accordingly, animals were sacrificed for analysis at day 4 after virus challenge with high viral loads [27]. Half the nasal turbinate, trachea, and lung tissues were used for viral load determination by quantitative SARS-CoV-2-specific RdRp/Hel RT-qPCR assay [16] and infectious virus titration by plaque assay [27]. For RT-qPCR, each sample was run in triplicate and the average cycle threshold (Ct) value was used for analysis of relative viral concentration.
Cryo-EM sample preparation and data acquisition
The purified SARS-CoV-2 S-B8 protein complexes were concentrated before being applied to the grids. Aliquots (4 μL) of the protein complex were placed on glow-discharged holey carbon grids (Quantifoil Au R1.2/1.3, 300 mesh). The grids were blotted and flash-frozen in liquid ethane cooled by liquid nitrogen with a Vitrobot apparatus (Mark IV, ThermoFisher Scientific). The grids sample quality was verified with an FEI Talos Arctica 200-kV electron microscope (Thermo Fisher Scientific). The verified grids with optimal ice thickness and particle density were transferred to a Titan Krios operating at 300 kV and equipped with a Cs corrector, a Gatan K3 Summit detector (Gatan Inc.) and a GIF Quantum energy filter (slit width 20 eV). Micrographs were recorded in the super-resolution mode with a calibrated pixel size of 0.54895 Å. Each movie has a total accumulated exposure of 50 e−/Å2 fractionated in 32 frames. The final image was binned 2-fold to a pixel size of 1.0979 Å. AutoEMation was used for the fully automated data collection. The defocus value of each image, which was set from −1.0 to −2.0 μm during data collection, was determined by Gctf. Data collection statistics are summarized in Table S11.
Cryo-EM data processing
The procedure for image processing of SARS-CoV-2 S-B8 complex is presented in Figure S2. In brief, Motion Correction (MotionCo2), CTF-estimation (GCTF) and non-templated particle picking (Gautomatch, http://www.mrc-lmb.cam.ac.uk/kzhang/) were automatically executed by the TsinghuaTitan.py program (developed by Dr. fang Yang). Sequential data processing was carried out on RELION 3.0 and RELION 3.1. Initially, 2,436,776 particles were auto-picked by Gautomatch or RELION 3.0 from 4213 micrographs. After several 2D classifications, 1,451,176 particles were selected and applied for 3D classification with one class. Two different states were obtained after further 3D classification: 3 RBD in up conformation bound with B8 Fab (3u), and 2 up RBDs and 1 down RBD with each bound to a B8 Fab (2u1d). 616,799 particles for the 2u1d state and 351,095 particles for the 3u state were subjected to 3D auto-refinement, yielding final resolutions at 3.21 and 3.06 Å, respectively. Further CTF refinement and Bayesian polishing improved the resolution to 2.65 Å (2u1d, C1 symmetry) and 2.67 Å (3u, C3 symmetry) with better map quality. To improve the RBD-B8 portion map density, focused local search classification was applied for each RBD-B8 portion with an adapted soft mask. The best classes for each RBD-Fab portion were selected and yielded a final resolution at 3.56 Å (RBD-Fab1, up), 3.34 Å (RBD-Fab2, up), 3.69 Å (RBD-Fab3, down), 3.87 Å (RBD-Fab3, up) from 479,305, 508,653, 656,429, and 136,482 particles, respectively. Further CTF refinement and Bayesian polishing improved the resolution of RBD-Fab2 to 3.11 Å with better map quality. RBD-Fab maps were fitted onto the whole structure map using Chimera, then combined using PHENIX combine_focused_maps. The reported resolutions were estimated with the gold-standard Fourier shell correlation (FSC) cutoff of 0.143 criterion. Data processing statistics are summarized in Table S11.
Model building and structure refinement
The spike model (PDB code: 6VSB) and the initial model of the B8 Fab generated by SWISS-Model were fitted into the EM density map, and further manually adjusted with Coot. Glocusides were built manually with carbohydrate tool in Coot. The atomic models were refinement using Phenix in real space with secondary structure and geometry restraints. The final structures were validated using Phenix.molprobity. UCSF Chimera, ChimeraX and PyMol were used for map segmentation and figure generation. Model refinement statistics are summarized in Table S11.
SARS-CoV-2 infection of reconstructed human nasal epithelia
MucilAirTM, corresponding to reconstructed human nasal epithelium cultures differentiated in vitro for at least 4 weeks, were purchased from Epithelix (Saint-Julien-en-Genevois, France). The cultures were generated from pooled nasal tissues obtained from 14 human adult donors. Cultures were maintained in air/liquid interface (ALI) conditions in transwells with 700 µL of MucilAirTM medium (Epithelix) in the basal compartment, and then kept at 37 °C under a 5% CO2 atmosphere. SARS-CoV-2 infection was performed as previously described [32]. Briefly, the apical side of ALI cultures was washed 20 min at 37 °C in MucilairTM medium to remove mucus. Cells were then incubated with 104 plaque-forming units (PFU) of the isolate BetaCoV/France/IDF00372/2020 (EVAg collection, Ref-SKU: 014V-03890; kindly provided by S. Van der Werf). The viral input was diluted in DMEM medium to a final volume 100 µL, and then left on the apical side for 4 h at 37 °C. Control wells were mock treated with DMEM medium (Gibco) for the same duration. Viral inputs were removed by washing twice with 200 µL of PBS (5 min at 37 °C) and once with 200 µL MucilairTM medium (20 min at 37 °C). The basal medium was replaced every 2–3 days. Apical supernatants were harvested every 2–3 days by adding 200 µL of MucilairTM medium on the apical side, with an incubation of 20 min at 37 °C prior to collection. For IgA treatment, cultures were washed once and then pretreated with antibodies added to the apical compartment for 1 h in 50 µL. Viral input was then directly added to reach a final volume of 100 µL. The antibodies were added again at day 2 d.p.i. in the apical compartment during an apical wash (20 min at 37 °C). To test the effect of dIgA treatment in the presence of mucus, dIgA were added directly to the apical compartment of MucilAirTM cultures without an initial wash. After IgA treatment for 1 h, the virus was added directly to the IgA/mucus mixture and left on the apical side for 4 h at 37°C. After viral inoculation, a single brief wash was made to remove the viral input while limiting mucus loss. The cultures were then maintained as in the no-mucus condition.
Viral RNA quantification in reconstructed human nasal epithelia
Apical supernatants were collected, stored at −80 °C until thawing and were then diluted 4-fold in PBS in a 96-well plate. Diluted supernatants were inactivated for 20 min at 80 °C. For SARS-CoV-2 RNA quantification, 1 µL of diluted supernatant was added to 4 µL of PCR reaction mix. PCR was carried out in 384-well plates using the Luna Universal Probe One-Step RT-qPCR Kit (New England Biolabs) with SARS-CoV-2 NP-specific primers (Forward 5'-TAA TCA GAC AAG GAA CTG ATT A-3’; Reverse 5'-CGA AGG TGT GAC TTC CAT G-3’) on a QuantStudio 6 Flex thermocycler (Applied Biosystems). A standard curve was established in parallel using purified SARS-CoV-2 viral RNA.
Histopathology and immunofluorescence (IF) staining
The lung and nasal turbinate tissues collected at necropsy were fixed in zinc formalin and then processed into paraffin-embedded tissue blocks. The tissue sections (4 μm) were stained with hematoxylin and eosin (H&E) for light microscopy examination as previously described with modifications [26]. For identification and localization of SARS-CoV-2 nucleocapsid protein (NP) in organ tissues, immunofluorescence staining was performed on deparaffinized and rehydrated tissue sections using a rabbit anti-SARS-CoV-2-NP protein antibody together with an AF488-conjugated anti-rabbit IgG (Jackson ImmunoResearch, PA, USA). Briefly, the tissue sections were first treated with antigen unmasking solution (Vector Laboratories) in a pressure cooker. After blocking with 0.1% Sudan black B for 15 min and 1% bovine serum albumin (BSA)/PBS at RT for 30 min, the primary rabbit anti-SARS-CoV-2-NP antibody (1:4000 dilution with 1% BSA/PBS) was incubated at 4°C overnight. This step was followed by incubation with a FITC-conjugated donkey anti-rabbit IgG (Jackson ImmunoResearch) for 30 min and the sections were then mounted in medium with 4’,6-diamidino-2-phenylindole (DAPI). For identification of DC-SIGN expression, we stained the NT slices with rabbit anti-DC-SIGN primary antibody (Abcam) and Alexa Fluor 488 goat anti-rabbit IgG (H+L) cross-adsorbed secondary antibody (Life Technologies) according to the manufacturer’s instructions. For identification of ACE2 expression, the goat anti-ACE2 primary antibody (R&D) and Alexa Fluor 568 donkey anti-goat IgG (H+L) secondary antibodies (Invitrogen) according to the manufacturer’s instructions. All tissue sections were examined, and the fluorescence images and whole section scanning were captured using 5×, 10× and 20× objectives with Carl Zeiss LSM 980. NP+ cells per field were quantified based on the mean fluorescence intensity (MFI) using the ZEN BLACK 3.0 and ImageJ (NIH).
Effects of B8-digA on SARS-CoV-2 infection in HK2 cells
HK2 cells were seeded into 24-well plates at the 40–50% confluency and cultured overnight. The B8-dIgA or control dIgA at the concentration of 1, 10, 100, 1000 ng/ml and then mixed with SARS-CoV-2 (1:10 TCID50) and incubated for 1 hour at room temperature. The antibody/virus mixture was then added to HK2 cells after the cell culture medium was removed and washed with PBS once and incubated for 1 hour at 37°C. The infectious medium was replaced with fresh medium containing respective concentration of antibody after washing 3 times with PBS. 24 h later, the infected cells were imaged under fluorescence microscope after staining with AF488-conjugated anti-SARS-CoV-2 NP antibody. Alternatively, the infected cells were lysed and blotted for SARS-CoV-2 NP protein to determine the extent of infection. Tubulin was blotted as the internal control.
B8-dIgA mediated enhancement via CD209
HEK293 T cells were seeded into 10-cm dish at 40% confluency and cultured overnight. The HEK293 T cells were transfected with human CD209 (Sino Biological) at 70–90% confluency. The expression of CD209 was measured by flow cytometry. The transfected HEK293T-CD209 cells were seeded into 96-well plates with 2.4×104 cells per well and cultured overnight. The HEK293T-CD209 cells were pre-treated with 10 ng/ml of B8-dIgA or control dIgA and incubated for 6 h prior SARS-CoV-2 infection (MOI=0.05). 24 h later, cells were then fixed with 4% paraformaldehyde solution for 30 min and air-dried in the BSC. Cells were further permeabilized with 0.2% Triton X-100 and incubated with cross-reactive rabbit sera anti-SARS-CoV-2-N for 1 hour at RT before adding Alexa Fluor 488 goat anti-rabbit IgG (H+L) cross-adsorbed secondary antibody (Life Technologies). The fluorescence density of SARS-CoV-2 infected cells was acquired using a Sapphire Biomolecular Imager (Azure Biosystems) and then the MFI of four randomly selected areas of each sample was quantified using Fiji software (NIH).
MoDC-mediated cell-to-cell transmission of SARS-CoV-2 pseudovirus
Monocytes were isolated from human healthy donor PBMCs by CD14+ microbeads and further differentiated into monocyte-derived dendritic cells (MoDCs) in RPMI1640 median containing 10% FBS, GM-CSF (40 ng/ml) and IL-4 (40 ng/ml) for 6 days. 25 ng/ml of hIgAs were treated with MoDCs for 1 hour prior to staining with Zombie aqua dye, BV421 anti-human CD1a, APC/Fire750 anti-human CD209 and PE anti-human IgA. Stained cells were acquired by ACEA NovoCyte Quanteon Flow Cytometry and analyzed with FlowJo software (v10.6) (BD Bioscience). For immunofluorescence assay, MoDCs were seeded on 8-well chamber slide and pre-treated with 25 ng/ml of hIgAs. At 1-hour post-treatment, MoDCs were washed twice by PBS prior to fixation by cold methanol for 20 minutes and incubated with blocking buffer (PBS with 5% FBS and 0.3% Triton X-100) for 1 hour at 4°C. MoDCs were further staining with rabbit anti-human CD209 antibody (Abcam, 1:500) overnight at 4°C. MoDCs were then stained with Goat AF488 anti-rabbit IgG, mouse PE anti-human IgA and DAPI for 1 hour and scanning using 20× objectives with LSM900 confocal microscope (ZEISS). For cell-to-cell transmission assay, 5×104 MoDCs were pre-treated with 25 ng/ml of different hIgAs for 1 hour and washed twice by PBS prior to SARS-CoV-2 pseudovirus infection. At 2-hour post-infection, MoDCs were washed twice by PBS and further co-cultured with 1.5×105 Vero-E6-TMPRSS2 cells. 48 hours later, infected cells were lysed to measure luciferase activity as described as above.
Effects of B8 antibodies on SARS-CoV-2 mediated cell–cell fusion
Vero-E6 TMPRSS2 cells were seeded into 48-well plates and cultured overnight. After treatment with B8 antibodies at the dose of 3000 ng/ml for 1 hour, HEK293 T cells transfected with SARS-CoV-2 spike-GFP were added into the treated Vero-E6 TMPRSS2 cells and co-cultured for 48 hours. The cell–cell fusion between Vero-E6 TMPRSS2 and HEK293T-Spike-GFP was then determined under a fluorescence microscope (Nikon ELIPSE) and the images of randomly selected region were captured using 4× and 10× objectives using the Nikon software.
Re-analysis of published nasal brushing single-cell data
The preprocessed scRNA-seq data from nasal brushing samples of 2 healthy controls and 4 COVID-19 patients were downloaded from Gene Expression Omnibus (GEO) database with accession numbers GSE171488 and GSE164547. Quality control metrics were consistent with the original article [33] and performed based on the R package Seurat (version 4.0.3) [34]. Harmony [35] was used to integrate the samples based on the top 4000 most variable genes obtained with the FindVariableFeatures function in Seurat. CD14+ (monocyte) cells were extracted for further analysis. The annotation of the cell type was performed by manually checking the marker genes of each cluster identified from the FindAllMarkers function in Seurat.
Quantification and statistical analysis
Statistical analysis was performed using PRISM 6.0 or later. Ordinary one-way ANOVA and multiple comparisons were used to compare group means and differences between multiple groups. Unpaired Student's t tests were used to compare group means between two groups only. A P-value <0.05 was considered significant. The number of independent experiments performed, the number of animals in each group, and the specific details of statistical tests are reported in the figure legends and the Methods section.
Results
Characterization of human monoclonal antibodies from COVID-19 patients
To isolate human monoclonal antibodies (MAbs), we obtained peripheral blood mononuclear cells (PBMCs) from one acute (P4) and three convalescent (P1-P3) COVID-19 patients in Hong Kong at a mean 73.5 (±25) days after symptoms onset (Table S1) between February and May 2020 when the ancestral SARS-CoV-2 was prevalent in Hong Kong. Enzyme-linked immunosorbent assay (ELISA) and pseudovirus neutralization assays revealed that all patient sera showed SARS-CoV-2 RBD- and spike-specific binding (Figure 1A and B) and neutralizing antibody (NAb) activities (Figure 1C). The mean NAb IC50 titer was 1:1753 with a range of 1:638-1:5701. Flow cytometry was then used to sort SARS-CoV-2-specific immunoglobulin positive (IgG+) memory B cells from individual PBMC samples using two fluorescent-conjugated RBD probes. The percentage of RBD-binding IgG+ memory B cells ranged from 0.19% to 0.52% (Figure S1A and S1B). We successfully cloned a total of 34 MAbs from these patients, including 3 from P1, 8 from P2, 17 from P3 and 6 from P4. We confirmed that 18 of these MAbs exhibited RBD-specific binding activities detected by ELISA (Figure S1C). No clear dominance of heavy (H) chain gene family was found among these 4 subjects by sequence analysis (Figure 1D, left panels). VLK1, however, was the most used variable gene family for the light (L) chain (Figure 1D, right panels). The average somatic hypermutation (SHM) rate ranged from 0% to 12.2% for the H chain and from 0.7% to 7.9% for the L chain (Figure 1E, left). The average complementarity-determining region 3 (CDR3) lengths ranged from 12.3 to 17.4 for the H chain and 8.4–9.4 for the L chain, respectively (Figure 1E, right). These results suggested overall comparable degrees of affinity maturation in these RBD-specific human MAbs obtained from individual memory B cells of four Hong Kong patients.
Specificity and potency of SARS-CoV-2-specific human neutralizing antibodies
To determine the antiviral activities of these 18 RBD-specific human MAbs, we performed binding and neutralization assays. Five of them, namely A6-IgG1, B4-IgG1, B7-IgG1, B8-IgG1 and C5-IgG1 displayed RBD- and spike-specific binding by ELISA (Figure 1F) and neutralizing activities against both pseudotyped and authentic viruses (Figure 1G and Figure S1D). Interestingly, the four most potent HuNAbs, B4-IgG1, B7-IgG1, B8-IgG1 and C5-IgG1, all came from patient P3 (Figure S1D). Sequence analysis revealed strong similarities between B7-, B8-, and C5-IgG1, which all contained an IGHV1-69 heavy chain gene and an IGKV3 kappa light chain gene. B4-IgG1 contained distinct IGHV3-66 and IGKV1-33 genes with CDR3 lengths of 12 amino acids (aa) and 9 aa, and somatic hypermutation (SHM) rate of 3.8% and 4.6%, respectively (Table S3). B7 and B8 were the most similar, as both contained IGHV1-69 and IGKV3-11 though B8 had a shorter CDR3 (14 aa vs 18 aa) and higher SHM (4.8% vs 0.0%) than those of B7 in the heavy chain. In the L chain, B7 and B8 shared a 9 aa CDR3, but B8 had a higher SHM rate than that of B7 (1.7% vs 0.7%) C5-IgG1 had a similar IGHV1-69 light chain gene with a 16 aa CDR3 and 2.4% SHM, but a different IGKV3-20 with a 9 aa CDR3 and 2.8% SHM. By ELISA, these four P3-derived HuNAbs bound to the SARS-CoV-2 RBD with half-maximal effective concentration (EC50) values ranging from 0.02 to 0.06 µg/ml, indicating stronger binding than that of the P4-derived A6-IgG1 (0.3 µg/ml) (Figure 1F, left and Table S4). P3-derived MAbs also exhibited stronger binding activities to the spike, with the EC50 values ranging from 0.018 to 0.06 µg/ml, compared to A6-IgG1 (17.94 µg/ml) (Figure 1F, right and Table S4). Neutralizing assays using ancestral SARS-CoV-2 pseudoviruses revealed that these four potent HuNAbs had IC50 values ranging from 0.0095 to 0.038 µg/ml, and IC90 values ranging from 0.046 to 0.136 µg/ml, respectively (Figure 1G, left and Table S5). Moreover, B8 proved to be the most potent HuNAb, capable of inhibiting authentic ancestral SARS-CoV-2 with an IC50 value of 0.013 µg/ml and an IC90 value of 0.032 µg/ml, respectively (Figure 1G, right and Table S5). We then determined if these HuNAbs could compete with ACE2 for RBD binding by surface plasmon resonance (SPR). We found that all of them strongly competed with ACE2 (Figure 1H). In line with the ELISA results, B8-IgG1 displayed the best KD value for RBD binding (169 pM) (Figure S2A, Table S6) and the strongest competition with ACE2 (Figure 1H), which explained its potent neutralizing activity. As B4-IgG1 displayed only partial competition for RBD binding with the other antibodies (Figure S2B), we further performed antibody synergy experiments using the pseudotype neuralization assay. No significant synergistic effects were found between any pairs of these four HuNAbs (Figure S2C). These results demonstrated that B4-, B7-, B8- and C5-IgG1 were all RBD-specific and competed with ACE2 for the same binding sites on the RBD.
B8-IgG1 pre-exposure prophylaxis and post-exposure treatment in Syrian hamsters
To determine the efficiency of B8-IgG1 in pre-exposure prophylaxis and post-exposure treatment against live intranasal SARS-CoV-2 infection, we administered B8-IgG1 intraperitoneally in golden Syrian hamsters, before or after viral challenge in our Biosafety Level-3 (BSL-3) animal laboratory (Figure 2A). In the pre-exposure prophylaxis group (G1, n=4), each hamster received a single intraperitoneal injection of 1.5 mg/kg B8-IgG1. In the post-exposure treatment groups, each hamster received a single intraperitoneal injection of 1.5 mg/kg B8-IgG1 at day 1 (G2, n=4), day 2 (G3, n=4) or day 3 (G4, n=4) postinfection (dpi), respectively (Figure 2A). The challenge dose was 105 plaque-forming units (PFU) of live ancestral SARS-CoV-2 (HKU-001a strain) [26,27]. Another group of hamsters (G0, n=4) received PBS injection as a no-treatment control. One G4 animal died accidentally during the procedure. Since Syrian hamsters recover quickly from SARS-CoV-2 infection, with resolution of clinical signs and clearance of virus shedding within one week after infection [13,27], we chose to sacrifice the animals at 4 dpi for HuNAb efficacy analysis, at a time when high viral loads and acute lung injury were consistently observed. At 4 dpi, NT and lung tissues were harvested to quantify infectious viruses by measuring PFUs, viral RNA loads by real-time reverse-transcription polymerase chain reaction (RT–PCR) and infected cells by immunofluorescence (IF) staining of viral nucleocapsid protein (NP)-positive cells as we described previously [26,27].
We found that infectious virus, measured by PFU, was readily detected in all tissue compartments of G0 hamsters but not in the lungs of 75% G1, 100% G2, 75% G3 and 0% G4 animals, nor in the NT of 50% G1, 75% G2, 50% G3 and 25% G4 animals (Figure 2B and and2E).2E). The decrease in PFU was of 2–3 orders of magnitude, suggesting efficient viral suppression in the lungs for the G1, G2, and G3 groups. A sensitive RT–PCR assay further demonstrated that viral RNA copy numbers were decreased in the lungs by 3 orders of magnitude in 50% of G1 hamsters (Figure 2C). In contrast, there was no significant viral RNA reduction in the NT of G1 animals (Figure 2F), suggesting lower efficacy of B8-IgG1 to prevent viral entry in the URT than in the lungs. There were slight but not significant viral RNA load decreases in both lungs and NT of G2 and some G3 animals. We then evaluated the number of infected cells or foci in these two tissue compartments by anti-NP antibody staining. A clear decrease of NP-positive cells or foci was observed in the lungs of G1 and G2 hamsters (Figure 2D). Abundant NP-positive cells with a diffuse distribution, however, were readily detected in the NT of all the challenged hamsters (Figure 2G). These results demonstrated that systemic B8-IgG1 injection was effective at reducing productive SARS-CoV-2 infection in the lungs when used for pre-exposure prophylaxis and early treatment especially within 48 hours post infection, but was insufficient to prevent viral infection in the NT.
To determine the correlates of B8-IgG1-mediated protection, we also measured the antibody concentration in serum, lung homogenate and NT homogenates at 0 and 4 dpi for all experimental animals. On average, 4,257 ng/ml and 2,101 ng/ml B8-IgG1 were found in animal sera at 0 and 4 dpi (Table S7). On 4 dpi, lung and NT homogenates contained 128 ng/ml and 20 ng/ml in G1, 238 ng/ml and 86 ng/ml in G2, 229 ng/ml and 93 ng/ml in G3, and 192 ng/ml and 46 ng/ml in G4 animals, respectively. These results demonstrated that most animals maintained higher peripheral B8-IgG1 antibody concentration, while a decreasing concentration gradient was observed in the lungs and NT during infection. The concentrations of B8-IgG1 measured in NT homogenates were close to the neutralization IC90 measured in vitro (32 ng/ml), explaining why infectious virus was undetectable in the PFU assay. These findings are in line with results obtained for other potent IgG HuNAbs administered systemically, as reported in our previous study [26]. The B8-IgG1 concentrations measured in the NT appeared insufficient to completely block infection in vivo, as indicated by the presence of NP-positive cells scattered throughout the NT.
Pre-exposure prophylaxis by monomeric B8-migA1 and B8-migA2 in Syrian hamsters
Since systemic administration of the RBD-specific neutralizing B8-IgG1 did not suppress SARS-CoV-2 nasal infection effectively, we sought to construct various types of IgA for mucosal intervention. For this purpose, we engineered B8-IgG1 into monomeric B8-mIgA1 and B8-mIgA2, and then into dimeric B8-dIgA1 and B8-dIgA2 by introducing the J chain. By in vitro characterization, purified B8-mIgA1, B8-mIgA2, B8-dIgA1 and B8-dIgA2 retained similar binding to RBD and spike by ELISA as compared to B8-IgG1 (Figure S3A-B, Table S8), and comparable antiviral activities based on neutralization IC50 and IC90 values (Figure S3C-D and Table S9). That said, B8-mIgA2 and B8-dIgA2 showed slightly more potent IC90 activities than B8-IgG1 in the pseudovirus neutralization assay using 293T-ACE2 cells as targets (Table S9). After introducing J chain, the proper dimer formation of B8-dIgA1 and B8-dIgA2 was confirmed by size exclusion chromatography analysis (Figure S3E-F). Furthermore, B8-mIgA1, B8-mIgA2, B8-dIgA1 and B8-dIgA2 also retained comparable competition with ACE2 for binding to spike by SPR analysis (Figure S3G-J). Therefore, the engineered IgA had the expected structural properties, and showed antiviral activities as potent as those of B8-IgG1 in vitro.
We then evaluated the monomeric B8-mIgA1 and B8-mIgA2 in the hamster model, using a higher 4.5 mg/kg dose via either intranasal or intraperitoneal injection, and used the same amount of intranasal B8-IgG1 inoculation as a control (Figure 3A–I). Interestingly, while changes in total RdRp and subgenomic sgNP viral RNA loads were not obvious (Figure 3B–C), B8-IgG1 and B8-mIgA1 (both i.n. and i.p.), but not B8-mIgA2, were able to significantly suppress infectious virus production (PFU) in the lungs of 75% infected hamsters by 2 orders of magnitude (Figure 3D). Sporadic infected cell foci were still detected in the lung sections by anti-NP staining (Figure 3E and Figure S4A), suggesting that protection conferred by B8-IgG1 and B8-mIgA1 was not complete. B8-mIgA2 was neither able to suppress viral RNA load nor PFU in the challenged hamsters, regardless of the route of antibody injection. On the other hand, like B8-IgG1, both B8-mIgA1 and B8-mIgA2 did not achieve significant viral suppression in the NT. After intranasal administration of either B8-mIgA1 or B8-mIgA2, some hamsters even showed a trend of slightly increased infectious virus production in NT, though this did not reach statistical significance (Figure 3F–H). Among these animals, the NP-positive cells were detected readily in the NT, as demonstrated by the whole section scanning, indicating the comparable distribution compared with the B8-IgG1 and no-treatment groups (Figure 3I). These results were consistent with many NP-positive cells observed in diffusely infected areas of NT by classic IF (Figure S4B). These results demonstrated that B8-mIgA1 was more potent than B8-mIgA2 at limiting SARS-CoV-2 infection in the lungs, but that both mIgA neither prevented nor significantly limited SARS-CoV-2 nasal infection.
B8-dIgA1 and B8-digA2 mediate robust SARS-CoV-2 nasal infection and injury in Syrian hamsters
Next, we tested the effects of dimeric B8-dIgA1 and B8-dIgA2 in Syrian hamsters. To improve protective efficacy of intranasal dIgA treatment, we included a 3-fold higher dosage group of 13.5 mg/kg besides the 4.5 mg/kg group (Figure 4A–J), and we shortened the interval between dIgA and virus inoculation from 24 to 12 hours (Figure 4A). Both RdRp and sgNP viral RNA loads dropped significantly in the lungs of hamsters that received the high dose of B8-dIgA2 compared to the no-treatment group (Figure 4B and C). Both B8-dIgA1 and B8-dIgA2 at the high dose also suppressed infectious viruses (PFU) in the lungs of 75% and 100% treated hamsters, respectively (Figure 4D). High dose B8-dIgA1 and B8-dIgA2 also decreased the number of NP-positive cells or foci in the lungs, with a more marked change for B8-dIgA2 (Figure 4I and Figure S5A). Unexpectedly, however, we observed significantly enhanced SARS-CoV-2 nasal infection and tissue damage in most infected hamsters included the low and high dose groups of both B8-dIgA1 and B8-dIgA2, in all the four assays used (Figure 4E–G). High dose administration of B8-dIgA1 or B8-dIgA2 resulted in increased PFU production in the NT by 37-fold and 81-fold, respectively, compared to the no-treatment group (Figure 4G). Since our model showed comparable NT PFU on day 2 and day 4 as described previously [27], this level of enhanced infection was unusual. It was also not observed with B8-IgG1 or monomeric B8-mIgA1 and B8-mIgA2 treatment, as described above. Moreover, the distribution of NP-positive cells in hamsters treated with dimeric B8-IgA was broader and reached deeper into NT tissue compared to the control group by the whole NT section scanning (Figure 4J and Figure S5B), which was associated with more severe and extensive epithelium desquamation and luminal cell debris (Figure 4H). The density of nasal NP+ cells was also significantly higher in B8-dIgA2-treated hamsters than in PBS-treated animals (Figure S5C-D). It is therefore conceivable that treatment with dimeric B8-IgA expanded the type and distribution of target cells in the nasal epithelium. Critically, we confirmed that control dimeric dIgA1 and dIgA2 did not enhance SARS-CoV-2 infection under the same experimental conditions (Figure S6).
To validate the role of B8-dIgA1 and B8-dIgA2, we then measured B8-IgA concentrations in the serum at day 0 and 4 dpi, and in the lung and NT homogenates at 4 dpi. B8-dIgA1 and B8-dIgA2 were primarily detected in lung homogenates at 4 dpi and were apparently undetectable in the serum and NT homogenates (Table S10). The enhanced viral replication in NT probably exhausted B8-dIgA1 and B8-dIgA2 locally, through antibody-virus complex formation and clearance [36]. To address this possibility, we treated separately five groups of naïve Syrian hamsters (n=4 per group) with each antibody at the 4.5 mg/kg dose. Twelve hours after the inoculation, antibody concentrations were readily detected in each tissues compartment (Table S11). The highest concentrations of dIgA antibodies were found in lung homogenates, followed by nasal washes, NT homogenates, and serum. These results suggested that B8-dIgA1 and B8-dIgA2 concentrations in both NT and lungs were still above their neutralization IC90 values at the time of viral challenge, indicating that the results obtained in our experiments could not be explained by the insufficient amounts of antibodies.
In addition to B8-dIgA1 and B8-dIgA2, we also cloned a native dimeric IgA1 (TH335-dIgA1) from one ancestral SARS-CoV-2-infected patient’s the IgA+ B cells to provide clinical relevance (Figure 4K). TH335-dIgA1 exhibited an IC50 value of 0.0019 µg/ml against authentic SARS-CoV-2 in vitro, which was about 7.8-fold stronger than B8-dIgA1 (Table S9). Like B8-dIgA1 and B8-dIgA2, TH335-dIgA1 suppressed viral infection in lungs potently (Figure 4L–O). No significant differences in viral RNA and infectious virus PFU were found between ctrl dIgA1 and TH335-dIgA1 (Figure 4P–R) but a higher number of infectious viruses was found in NT of one animal. Notably, the distribution of NP-positive cells in hamsters treated with TH335-dIgA1 was also broader and reached deeper into NT tissue compared to control hamsters by the whole NT section scanning (Figure 4S). Taken together, our results consistently demonstrated that pre-exposure dimeric B8-dIgA1, B8-dIgA2 and TH335-dIgA1 treatment did not suppress SARS-CoV-2 nasal infection and injury in Syrian hamsters.
B8-dIgA1- and B8-digA2-mediated SARS-CoV-2 infection via CD209
Since B8 antibodies share the same binding site to the RBD domain, we sought to investigate possible mechanisms of B8-dIgA1- and B8-dIgA2-mediated infection. First, we consistently found that 10 ng/ml B8-dIgA1 or B8-dIgA2 completely neutralized SARS-CoV-2 infection in human renal proximal tubule cells (HK-2), as measured in our previously reported immunofluorescence assay (Figure S7) [34]. We also tested B8-mIgA2 and B8-dIgA2 neutralizing capacity in the MucilAirTM model, consisting of a reconstructed human nasal epithelium, which contained goblet, ciliated, and basal cells (Figure 5A) [32]. Both B8-mIgA2 and B8-dIgA2 neutralized SARS-CoV-2 in a dose-dependent fashion, when compared to a dIgA2 control antibody. Similar experiments were carried out in the presence of the mucus naturally secreted by goblet cells to determine whether dIgA interaction with the mucus might alter their neutralization capacity. B8-dIgA2 showed the same neutralization capacity in the presence and absence of mucus. These results demonstrated that B8-dIgA1 and B8-dIgA2 did not enhance SARS-CoV-2 infection in either human HK-2 or primary airway epithelial cells, which primarily expressed human ACE2 as a viral receptor. We then turned our attention to ACE2-independent mechanisms that might be associated with dimeric IgA-mediated enhancement of SARS-CoV-2 infection. Considering that mucosal monocyte-derived dendritic cells (DC) could mediate SARS-CoV infection and dissemination in rhesus monkeys as early as 2 dpi as we previously reported [37], we sought to investigate the role of DC-expressed surface receptors. We focused on CD209 (DC-SIGN) because this lectin was previously shown to act as a cellular receptor for secretory IgA [38]. By IF staining, intranasal administration of B8-dIgA2 alone did not increase CD209 expression in the NT of treated hamsters (Figure 5B, left). Upon SARS-CoV-2 infection, however, we noted an increase in CD209-positive cells in olfactory epithelium devoid of ACE2 expression (Figure 5B, middle). Importantly, most CD209+ cells were positive for NP (Figure 5B, right), indicating that these CD209+ cells were likely permissive to SARS-CoV-2. We then determined whether B8-dIgA1 and B8-dIgA2 could enhance SARS-CoV-2 infection in 293T cells expressing human CD209 or CD299 but not ACE2. Using a low MOI of 0.05, we found that pre-incubation of B8-dIgA1 and B8-dIgA2 enhanced live SARS-CoV-2 infection significantly in cells expressing human CD209, as determined by increased viral NP production (Figure 5C). Interestingly, human CD299, a type II integral membrane protein that is 77% identical to CD209, did not show similar activities in the same experiment (Figure 5C). Control dIgA1 and dIgA2 did not show any enhancement in NP+ cell detection compared with virus only. Considering that CD209+ DCs promote HIV-1 trans-infection of CD4+ T cells via cell–cell contacts, we speculated that B8-dIgA1 and B8-dIgA2 might not be able to block the similar process for SARS-CoV-2. Indeed, by testing the B8 antibodies at concentrations 100-times higher than IC90 neutralization values (around 3000 ng/ml), none of B8-IgG1, B8-mIgA1, B8-mIgA2, B8-dIgA1 and B8-dIgA2 could block cell–cell fusion (Figure 5D). Taken together, our results demonstrated that B8-dIgA1- and B8-dIgA2-mediated SARS-CoV-2 infection likely involved viral capture and infection of mucosal CD209+ cells, followed by more robust trans-infection of ACE2+ epithelial cells through cell–cell transmission in NT.
SARS-CoV-2 hijacks B8-digA1 for human MoDC-mediated trans-infection of target cells
Studies have indicated that SARS-CoV-2 did not replicate in dendritic cells productively [39,40], but CD209 expressing on DC surface could capture and transfer the virus to permissive ACE2+ cells [41]. Since B8-dIgA showed the binding capacity to both CD209 and viruses, we sought to determine the effect of different B8-IgAs on human monocyte-derived dendritic cells (MoDCs)-mediated SARS-CoV-2 trans-infection of target cells. After incubating human MoDCs with B8-IgAs, we found that the subset of IgA+ MoDCs were all CD209 positive (Figure 6A). Moreover, B8-mIgA, B8-dIgA and Ctrl IgA could bind to MoDCs with similar mean frequencies (range 0.32%–1.01%) (Figure 6B). Confocal imaging then demonstrated the co-localization of these IgAs with CD209 on MoDCs (Figure 6C). These results showed that various types of IgA could interact with MoDCs via CD209. Next, we sought to determine MoDC-mediated trans-infection of Vero-E6 TMPRSS2 cells by B8-dIgA-bound pseudovirus. Human MoDCs were pretreated with each type of IgA and pulsed with pseudotyped SARS-CoV-2, respectively. After extensive washes, pulsed MoDCs were co-cultured with Vero-E6 TMPRSS2 cells (Figure 6D). As controls, MoDCs alone did not produce any positive signal with or without hIgAs treatment (Figure 6E). After co-culture with Vero-E6 TMPRSS2 cells, however, B8-dIgA1 promoted MoDC-mediated enhancement of trans-infection by a mean of 4.79-fold compared to Ctrl IgA1 (Figure 6F). B8-mIgA1, B8-mIgA2 or B8-dIgA2 did not show enhanced trans-infection (Figure 6F–G). These results indicated that B8-dIgA1 increased human MoDC-mediated SARS-CoV-2 trans-infection of target cells.
Cryo-EM analysis of the spike-B8 complex
To understand the potential mechanism of action of B8 HuNAb, we performed a cryo-EM single-particle analysis of B8 Fab bound to the SARS-CoV-2 spike ectodomain trimer (Figure S8). Two B8-spike complex structures were determined based on 351,095 and 616,799 particles collected, respectively (Table S12). One structure at 2.67 Å resolution contained the spike with all three RBDs adopting the “up” conformation (3u), where each “up” RBD was bound by one B8 Fab (Figure 7A and B). The other structure at 2.65 Å resolution contained one spike trimer with 2 RBDs in the “up” conformation and 1 RBD in the “down” conformation (2u1d), where each RBD was also bound by one B8 Fab despite the presence of two distinct RBD conformations (Figure 7C and D). After superimposing the “3u” and “2u1d” spikes, a ~53-degree rotation was observed between the “up” RBD (red color) in the 3u spike trimer and the “down” RBD (gray) in the 2u1d spike trimer (Figure 7E). The B8 Fab appeared to bind to the receptor-binding motif (RBM) of the RBD through its heavy chain for most of the interactions (Figure 7F and G). Therefore, the cellular receptor ACE2 would clash with the B8 Fab due to the overlap of their respective epitopes on the RBM (Figure 7F–H and Table S13). The elucidation of the epitope revealed that B8 could be grouped into the SARS-CoV-2 neutralizing antibody class II [42]. These structural findings were further supported by neutralization assays using a panel of pseudoviruses containing naturally occurring mutations. Indeed, the E484 K mutation from the South African SA9 strain, which is located within the B8-binding interface, caused a major loss of neutralizing potency for all the B8 isotypes tested: IgG1, mIgA1, mIgA2, dIgA1 and dIgA2 (Table S14). The comparable neutralization profiles of these NAbs against the full panel of viral variants also indicated that the conformation of key RBD-binding residues remained unchanged after engineering of the constant regions of these B8 isotypes.
Discussion
In this study, we investigated the preventive potential of a potent RBD-specific NAb B8 primarily in the forms of monomeric and dimeric IgA against live intranasal SARS-CoV-2 infection in the golden Syrian hamster model as compared with B8-IgG1 [26]. While these B8-IgA antibodies maintained neutralizing activities against SARS-CoV-2 like B8-IgG1, they displayed distinct in vivo effects, with clear differences in their capacity to modulate viral infection in the NT. Pretreatment by intranasal administration of monomeric B8-mIgA1 or B8-mIgA2 did not significantly reduce infectious virus production in the NT homogenates. On the contrary, the antibody isotype had a marked effect because intranasal administration of dimeric B8-dIgA1 or B8-dIgA2 paradoxically increased the number of infectious viruses in NT homogenates. TH335-dIgA1 also increased the number of infectious viruses in one animal together with consistent broader infection in NT. Such effect was not observed with several intranasal IgG HuNAbs previously tested by other groups or our team [26,40]. Mechanistically, instead of neutralization, virus-bound B8-dIgA1 and B8-dIgA2 used CD209 as an alternative receptor to infect non-ACE2 cells. CD209+ cells were increased and permissive to viral infection in the olfactory epithelium of Syrian hamsters upon SARS-CoV-2 infection, suggesting that this cell population could contribute to viral mucosal seeding. Indeed, we found that CD209 expressing cells could be infected in vitro by live SARS-CoV-2 at 0.05 MOI in the presence of B8-dIgA1 and B8-dIgA2. Since none of the B8-based MAbs could prevent SARS-CoV-2 cell-to-cell transmission, even at high concentrations in vitro, virus-laden mucosal CD209+ cells might trans-infect ACE2+ cells through cell-to-cell contacts in NT, resulting in broader infection and injury.
Cryo-EM analysis further indicated that B8 is a typical class II HuNAb that binds to the SRAS-CoV-2 spike RBD in either a 3u or a 2u1d mode. Our findings, therefore, reveal a previously unrecognized pathway for RBD-specific dimeric IgA-mediated SARS-CoV-2 nasal infection and injury in Syrian hamsters.
The role of dimeric IgA has been explored primarily for mucosal transmitted viruses. At the mucosal surface, the major IgA type is the secretory form, which is generated from dIgA by the acquisition of a secretory component upon endocytosis and secretion by epithelial cells. In the simian AIDS macaque model, neutralizing dIgA given directly into the rectal lumen can prevent viral acquisition in rhesus monkeys challenged via the mucosal route [43]. Although the administered dIgA did not contain the secretory component (SC), they might have associated with free SC, which is present in mucosal secretions such as human lung lavages [44]. Neutralizing dIgA1 and dIgA2 could be protective through several mechanisms, including direct virus neutralization, virion capture, or the inhibition of virion transcytosis across the epithelium [20]. In this macaque study, however, Watkins et al. demonstrated that the dimeric HGN194 dIgA2 protected only 1/6 animals in a rectal challenge model [43]. Recently, Taylor et al. found an increase in virion number and penetration depth in the transverse colon and mesenteric lymph nodes, after mucosal treatment with the HGN194 dIgA2 compared to a PBS control [45]. The authors suggested that virus-specific dIgA somehow mediated the delivery of virus immune complexes to the mesenteric lymph nodes for systemic infection. Here, we report that SARS-CoV-2 may subvert the action of potent neutralizing antibodies, as pretreatment with neutralizing B8-dIgA1 and B8-dIgA2 induced a robust nasal infection via a previously unrecognized mode of action.
SARS-CoV-2 engages CD209+ cells to evade ACE2-dependent neutralizing B8-dIgA1 and B8-dIgA2 for broader NT infection and injury. Previous studies have indicated various scenarios for antibody-dependent enhancement (ADE) in viral infections. The well-known dengue ADE has been associated with poorly neutralizing cross-reactive antibodies against a heterologous viral serotype, leading to increased infection of FcγR-expressing cells [46]. Recent findings suggested that an increase in afucosylated antibodies contribute to dengue ADE [47]. In contrast, vaccine-associated enhanced respiratory disease induced by respiratory syncytial virus has not been found to be antibody-dependent [48]. For SARS and MERS, ADE observed in vitro depended on binding of the antibody Fab to the virus and the binding of the Fc component to FcγR on target cells [49]. We previously found that spike-specific neutralizing IgG antibody abrogated wound-healing responses in SARS-CoV-infected Chinese macaques [28]. In the case of COVID-19, vaccination and passive immunization studies have not revealed ADE of disease severity [50]. Comprehensive studies, however, are necessary to define the clinical correlates of protective immunity against SARS-CoV-2, especially in the context of vaccine breakthrough infections. During the natural infection, one study indicated that the increase in afucosylated antibodies might contribute to COVID-19 severity [51]. To date, four classes of potent HuNAbs have been isolated from convalescent COVID-19 patients [42,49]. The molecular mechanism of neutralization for most potent HuNAbs was primarily through blocking the interaction between ACE2 and the spike RBD. Currently, systemic RBD-specific HuNAb treatment remains to be improved for therapeutic suppression of SARS-CoV-2 replication in the NT or URT, both in animal models and human trials [26,27,52]. One limitation is the insufficient amounts of HuNAbs distributed on the nasal mucosal surface for protection [26]. Other reasons might include alternative entry pathways engaged by SARS-CoV-2 to evade HuNAbs. To this end, Liu et al. reported recently that antibodies against the spike N-terminal domain (NTD) induced an open conformation of the RBD and thus enhanced the binding capacity of the spike to the ACE2 receptor, leading to increased viral infectivity [53]. Yeung et al. demonstrated that SARS-CoV-2 could engage soluble ACE2 (sACE2) and then bind alternate receptors for viral entry, through interaction between a spike/sACE2 complex with the angiotensin II AT1 receptor, or interaction between a spike/sACE2/vasopressin complex with the AVPR1B vasopressin receptor, respectively [34]. In this study, we found that, in the presence of neutralizing B8-dIgA1 or B8-dIgA2 antibodies, SARS-CoV-2 might use the cellular receptor CD209 for capture or attachment, which likely expanded the use of CD209+ cells, leading to NT infection and trans-infection. Interestingly, a preprint report suggests that cells expressing CD209 can be infected directly by SARS-CoV-2 through an interaction of the spike with the NTD instead of the RBD [54]. This mode of action, however, was unlikely to explain our findings, because no enhancement of SARS-CoV-2 nasal infection was found in presence of control dIgA1 and dIgA2. Our results rather suggested that the direct binding of virus-bound B8-dIgA1 or virus-bound B8-dIgA2 to CD209 might be a likely pathway, resulting in severe SARS-CoV-2 nasal infection and damage. In line with our results, a previous study demonstrated that dIgA itself can use CD209 as a cellular receptor [38]. During the viral entry process, since neither B8-dIgA1 nor B8-dIgA2 could prevent virus cell-to-cell transmission, infected mucosal CD209+ cells might enable a more efficient viral transmission to ciliated nasal epithelial cells, which show the highest expression of ACE2 and TMPRSS2 receptors [52]. In support of this notion, previous studies indicated that mucosal DCs can capture HIV-1 through binding of its envelope glycoproteins to CD209 and efficiently transfer the bound virions to CD4+ T cells, in a process called trans-enhancement or trans-infection [55]. The trans-infection markedly decreased the neutralization efficiency of potent NAbs directed at HIV [56]. Moreover, although monocyte-derived DCs (MoDCs) cannot support productive SARS-CoV-2 replication [40], we and others demonstrated that MoDCs could mediate efficient viral trans-infection of ACE2+ cells [41]. We further provided the experimental evidence that B8-dIgA1 could promote MoDC-mediated SARS-CoV-2 trans-infection of target cells. Since B8-dIgA2 did not show similar effects on MoDCs for trans-infection, the action of mucosal DCs in NT remained to be investigated. Nevertheless, our findings of dIgA-mediated nasal infection and injury may have significant implications for SARS-CoV-2 transmission, COVID-19 pathogenesis, and immune interventions.
Limitations of the study
As our in vivo findings were obtained in Syrian hamsters, it remains to be determined whether CD209+ DCs are abundantly recruited to the nasal mucosa in SARS-CoV-2-infected humans. It is, however, known that myeloid DCs are increased in the nasal epithelium upon infection [37,57]. Our preliminary analysis of the human nasal cytology data (under accession code EGAS00001004082) revealed the presence of increased CD209+ DCs in addition to abundant ACE2, TMPRSS2, and furin expression in the apical side of multiciliated cells of SARS-CoV-2-infected human subjects (Figure S9) [33]. It is also necessary to evaluate the effect of dimeric IgA in viral infection of CD209 knockout animals in the future. Another limitation is that we did not have an NTD-specific neutralizing dIgA for comparison. Besides the class II neutralizing dIgA such as B8-mIgA2 and B8-dIgA2, neutralizing dIgA belonging to other classes should also be investigated in future. We do not know whether other cellular receptor such as the polymeric immunoglobulin receptor (pIgR) plays a role in B8-dIgA-mediated SARS-CoV-2 nasal infection. The pIgR is known for transcytosis of soluble dIgAs and immune complexes from the basolateral to the apical epithelial surface. It remains uncertain whether intranasal B8-dIgA-mediated NT infection would lead to worse clinical outcomes, especially the Omicron and its subvariants have become the dominant strains globally with a faster and stronger viral replication compared to ancestral and Delta variant in the nasal and bronchial compartments [58]. It should be aware that current intramuscular vaccinations might not induce sufficient amounts of secretory dIgA at the nasal mucosal sites [59] and the extensive mutations in the spike protein of Omicron variants have conferred escape from ancestral virus-based vaccination and primary infection [60]. Finally, since our findings were primarily based on the testing of engineered neutralizing B8-IgA antibodies in hamster models, future studies remain necessary to determine the role of vaccine- or infection-induced dIgA in mediating breakthrough nasal infection in humans.
Supplementary Material
Author contributions
Conceptualization, Z.C.; HuNAb cloning, Z.B.; experimental design, Z.C., Z.B., R.Z., J.F.W.C.; hamster experiments, J.F.W.C., C.C.S.C., V.K.M.P., C.C.Y.C., K.K.H.C., and JC; cryoEM study, J.Z., J.G., Z.Y.W., X.W.; SPR experiments, Q.Z., S.S. and L.Z.; LL confocal imaging, L.L. and D.Z.; MucilAirTM experiment, R.R. and L.A.C.; HK-2 experiment, M.L.Y.; nasal cytology data analysis, M.Y. and R.S.; clinical specimens, K.K.W.T.; in vitro experiments, H.C., Z.D., K.K.A., H.H., H.O.M., J.C., C.L., J.Z.; manuscript preparation, Z.C., Z.B., R.Z., L.Z., K.Y.Y.; study supervision, Z.C., K.Y.Y. and L.Z.
Data availability statement
The data of this studies are available upon reasonable request and no custom computer code or algorithm used to generate results that are reported in the paper and central to its main claims.
References
Articles from Emerging Microbes & Infections are provided here courtesy of Taylor & Francis
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
Discover the attention surrounding your research
https://www.altmetric.com/details/152515573
Article citations
Ultrapotent class I neutralizing antibodies post Omicron breakthrough infection overcome broad SARS-CoV-2 escape variants.
EBioMedicine, 108:105354, 27 Sep 2024
Cited by: 0 articles | PMID: 39341153 | PMCID: PMC11470419
Breakthrough infections by SARS-CoV-2 variants boost cross-reactive hybrid immune responses in mRNA-vaccinated Golden Syrian hamsters.
PLoS Pathog, 20(1):e1011805, 10 Jan 2024
Cited by: 2 articles | PMID: 38198521 | PMCID: PMC10805310
Structural insights into broadly neutralizing antibodies elicited by hybrid immunity against SARS-CoV-2.
Emerg Microbes Infect, 12(1):2146538, 01 Dec 2023
Cited by: 4 articles | PMID: 36354024 | PMCID: PMC9817130
Persistent immune and clotting dysfunction detected in saliva and blood plasma after COVID-19.
Heliyon, 9(7):e17958, 04 Jul 2023
Cited by: 2 articles | PMID: 37483779 | PMCID: PMC10362241
A broadly neutralizing monoclonal antibody overcomes the mutational landscape of emerging SARS-CoV-2 variants of concern.
PLoS Pathog, 18(12):e1010994, 12 Dec 2022
Cited by: 7 articles | PMID: 36508467 | PMCID: PMC9779650
Go to all (7) article citations
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
European Genome-Phenome Archive
- (1 citation) European Genome-Phenome Archive - EGAS00001004082
GEO - Gene Expression Omnibus (2)
- (1 citation) GEO - GSE171488
- (1 citation) GEO - GSE164547
Nucleotide Sequences
- (2 citations) ENA - MT230904
Protein structures in PDBe (2)
-
(1 citation)
PDBe - 6VSBView structure
-
(1 citation)
PDBe - 6M0JView structure
Electron Microscopy Data Bank (EMDB) at PDBe (3)
Protocols & materials
Related Immune Epitope Information - Immune Epitope Database and Analysis Resource
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Funding
Funders who supported this work.
Consultancy Service for Enhancing Laboratory Surveillance of Emerging Infectious Diseases and Research Capability on Antimicrobial Resistance for Department of Health of the Hong Kong (1)
Grant ID: na
Friends of Hope Education Fund (1)
Grant ID: na
Hainan talent development project (1)
Grant ID: SRC200003
Health and Medical Research Fund (1)
Grant ID: COVID190123
High Level-Hospital Program, Health Commission of Guangdong Province, China (1)
Grant ID: na
Hong Kong Research Grants Council Collaborative Research Fund (1)
Grant ID: C7156-20G, C1134-20G and C5110-20G
Innovation and Technology Fund (1)
Grant ID: na
National Key Research and Development Project of China (1)
Grant ID: 2020YFC0860600
National Program on Key Research Project of China (1)
Grant ID: 2020YFA0707500, 2020YFA0707504 and 2021YFC0866100
Sanming Project of Medicine in Shenzhen, China (1)
Grant ID: SZSM201911014
Shenzhen Science and Technology Program (1)
Grant ID: JSGG20200225151410198
Theme-Based Research Scheme (1)
Grant ID: T11-706/18-N and T11/709/21-N
Urgence COVID-19 Fundraising Campaign of Institute Pasteur (1)
Grant ID: TROPICORO
Wellcome Trust (1)
Grant ID: P86433