Abstract
The milder clinical manifestations of Omicron infection relative to pre-Omicron SARS-CoV-2 raises the possibility that extensive evolution results in reduced pathogenicity. To test this hypothesis, we quantified induction of cell fusion and cell death in SARS-CoV-2 evolved from ancestral virus during long-term infection. Both cell fusion and death were reduced in Omicron BA.1 infection relative to ancestral virus. Evolved virus was isolated at different times during a 6-month infection in an immunosuppressed individual with advanced HIV disease. The virus isolated 16 days post-reported symptom onset induced fusogenicity and cell death at levels similar to BA.1. However, fusogenicity was increased in virus isolated at 6 months post-symptoms to levels intermediate between BA.1 and ancestral SARS-CoV-2. Similarly, infected cell death showed a graded increase from earlier to later isolates. These results may indicate that, at least by the cellular measures used here, evolution in long-term infection does not necessarily attenuate the virus.
Introduction
The reduced incidence of severe disease reported with Omicron1 may result from the increasing immunity of the population because of vaccination and previous infections. Alternatively, the virus itself may have decreased its propensity to cause more severe disease2. If the ability of the Omicron virus itself to cause severe disease is attenuated independently of increased population immunity, the question arises of whether the virus is constrained to attenuate because of the evolutionary process through which variants emerge.
Mechanisms of variant formation may include reverse zoonosis3-12, the infection of an animal reservoir where the virus mutates to adapt to the new host species, then re-infection of a human host, or evolution in long-term infection in immunosuppressed individuals13-24. Evolution in long-term infection in immunosuppression is documented to occur in some people who are immunosuppressed because of advanced HIV disease13-15,24, defined as a CD4 T cell count < 200 cells/microliter in a person living with HIV.
Analysis of multiple long-term infections in immunosuppressed individuals has demonstrated recurrent mutations that are associated with escape from neutralizing antibodies22. However, mutations outside of spike are also common and may affect virus infection and cell-cell spread and may therefore impact pathogenicity22, with one possible outcome being that pathogenicity is reduced. In cell culture, viruses attenuate during long-term passaging, and such passaging is used to make live attenuated vaccines25.
SARS-CoV-2 infection can lead to disease in several ways26. One hallmark is the presence in the lung of syncytia, cells which underwent fusion through the interaction of the SARS-CoV-2 spike protein expressed on the infected cell surface with the human angiotensin converting enzyme 2 (ACE2) receptor on neighbouring cells27-30. Fusion results in multinucleated cells. This phenotype can be readily reproduced in cell culture infection with SARS-CoV-231-34.
The sequence of events which leads to the ability of the virus to enter cells by binding ACE2 on the plasma membrane is also required for infected cells to efficiently fuse to other cells. SARS-CoV-2 spike has two subunits, S1 and S2. The S1 subunit binds the ACE2 receptor, while S2 mediates membrane fusion35. Spike contains an S1/S2 cleavage site predominantly cleaved by the cellular furin protease28. This cleavage allows further cleavage at the S2’ site mediated by the cellular serine protease TMPRSS2, activating the S2 subunit for fusion35. A cathepsin-dependent alternate pathway for viral entry exists and enables TMPRSS2-independent viral infection28. The Omicron BA.1 subvariant does not have efficient S1/S2 cleavage and predominantly uses the alternative pathway to infect. Consequently, cell-cell fusion induced by BA.1 is lower than with ancestral SARS-CoV-2 and Delta variant infections27. Generally, pre-Omicron variants show enhanced fusogenicity compared to ancestral virus36-39, while Omicron subvariants BA.1 and BA.2 show decreased fusogenicity which correlates with decreased pathogenicity in hamster infections27,38,40,41. In the Omicron BA.4 and BA.5 subvariants, fusogenicity is higher relative to BA.1 and BA.238,40, and this is associated with higher pathogenicity in hamsters40.
SARS-CoV-2 infection also leads to death of infected cells42,43 which can be recapitulated in cellular assays44-46. Cell death may initiate an inflammatory response and lead to the trafficking of immune cells to the site of inflammation (reviewed in47,48). This causes increased lung fluid, cellularity, and later scarring, and in turn leads to less oxygen being able to diffuse into the blood, resulting in respiratory distress (reviewed in49). A reduction in cell death upon infection would therefore be predicted to lead to lower levels of inflammation and therefore lower disease severity.
We have previously reported on the evolution of SARS-CoV-2 from ancestral virus infection in an individual who was immunosuppressed because of advanced HIV disease13,14. Here we used timelapse microscopy and flow cytometry assays to determine changes in cell fusion and cell death induced by SARS-CoV-2 as it evolves over half a year of continuous infection. As a reference, we show cell fusion and cell death with Omicron BA.1 infection and ancestral virus with the D614G substitution.
We observed that virus from relatively early in the infection induced cell fusion and death at similar levels to Omicron BA.1. However, the virus which evolved over 6 months of infection showed fusogenicity intermediate between BA.1 and D614G and cell death induction more similar to D614G, indicating that, at least by these parameters, the virus did not attenuate during the course of evolution.
Results
We isolated live virus from an individual with advanced HIV disease13 (defined as CD4 T cell count of less than 200 cells per microliter) who was enrolled in our longitudinal cohort to investigate the immune response to SARS-CoV-2 infection. Participants in the study cohort were sampled as soon as practicable post-diagnosis, then within the first month post-diagnosis, and subsequently at 3-month intervals50. In the participant with advanced HIV, continuous SARS-CoV-2 infection was detected by qPCR from combined nasopharyngeal and oropharyngeal swabs for a period of 6 months from the date of diagnosis (Figure 1A), with virus isolation performed from day 6 post-diagnosis (isolate designated D6), the first available sample for isolation, up to day 190 post-diagnosis (designated D190). We have previously sequenced these viral isolates and tested them for escape from neutralizing antibodies elicited by SARS-CoV-2 infection13. D6 had low to moderate immune escape from plasma sampled from convalescent individuals previously infected with either ancestral SARS-CoV-2, the Beta, or the Delta variant. However, D190 showed more extensive immune escape from ancestral virus and Delta variant infection elicited neutralizing antibodies (Figure 1B, modified analysis with data from ref13). Phylogenetic analysis of the infection showed a pattern consistent with the evolution of one ancestral virus infection (Figure 1C). Viral isolates D6, day 20 (D20), day 34 (D34), day 106 (D106) and D190 showed mutations in spike and other viral genes relative to ancestral virus, with neutralization escape mutation E484K along with multiple other mutations already present in D6 and the neutralization escape mutations K417T and F490S51 being present in the D190 isolate (Figure 1D).
To test whether the virus has evolved other changes in addition to neutralizing antibody escape, we used the isolated viruses to infect the human H1299 lung cell line overexpressing the ACE2 receptor52. All results reported here are from live virus infections. This cell line has the endogenous histone H2AZ gene labelled with YFP by the insertion of the fluorophore sequence as an artificial exon into the first intron53, giving nuclear fluorescence. We used this nuclear fluorescence signal combined with automated image analysis to detect fused cells and cell number. We performed time-lapse microscopy with cells grown under controlled temperature and CO2 over 48 hours with images taken every 10 minutes. Uninfected cells grew until confluence with little evidence of cell death and cell-cell fusion (Video 1). In contrast, cells infected with an ancestral SARS-CoV-2 D614G isolate led to cell fusion, cell death and/or lack of cell division which became apparent about 12 hours post-infection (Video 2). These effects seemed less pronounced in cells infected with the Omicron BA.1 subvariant (Video 3) and in the D6 isolate from early after diagnosis (Video 4). However, the D190 isolate from 6 months later seemed to have increased fusogenicity and cytotoxic/cytostatic effects relative to BA.1 and D6 (Video 5).
We noted that cell nuclei become clustered together to form a contiguous region of fluorescence during cell fusion (Figure 2A). We used an automated image analysis pipeline (Figure 2 – figure supplement 1) to detect fused cells and cell number. Detection of fused cells was based on the observation that in the absence of fusion, individual cell nuclei are distinct even in confluent cell culture because they are separated by cellular cytoplasm. After fusion, cell nuclei are close together and form an area considerably larger than a single nucleus (Figure 2-figure supplement 1, Materials and methods).
Given that our image analysis pipeline was not designed to detect dead cells, we excluded the last 12 hours of the movies when extensive cell death occurred. Quantifying fusions over multiple independent experiments showed that uninfected cell culture had a low frequency of fusions which did not increase over time. Infection by the D614G virus showed an increasing fusion frequency, with about 40% of cell nuclei in fused cells by 36 hours post-infection. The frequency of fusions was lower with BA.1 infection throughout and reached less than half of that seen in the D614G infection. The pattern in D6 virus infections was similar to BA.1, while with D190 the frequency of fusions was intermediate between BA.1 and D614G infection (Figure 2B).
To quantify changes in cell number relative to the start of the movie, we used the number of nuclei as a measure of cell number Cell numbers in uninfected cultures increased until they were about 2.5-fold higher at 36 hours relative to the start of the movie. For the D614G infection, cell numbers stopped increasing about 20 hours post-infection and started to decrease. In BA.1 and D6 infection, the number of cells also stopped increasing after about 20 hours but did not decrease to the same extent as with ancestral virus. Infection with the evolved D190 virus from late in the infection showed a similar but less pronounced decline in cell numbers as D614G (Figure 2C). Excluding nuclei in fusions from the results gave a similar pattern (Figure 2-figure supplement 2).
We used a second assay to detect cell death at 24 hours post-infection. We used this relatively early timepoint because to avoid the effect of multiple infections per cell when infection is saturating54,55, which happens later (see for example Video 2). Infection was detected by staining for SARS-CoV-2 nucleocapsid, and the fraction of dead infected cells determined by co-staining with a death detection dye. The positive control was addition of ethanol (Figure 3A). We compared infection by ancestral virus to BA.1 and observed that, while the fraction of infected cells was similar and slightly higher in BA.1 relative to ancestral virus, the fraction of the infected cell population staining positive for the death detection dye was lower in BA.1 (Figure 3A).
We tested two independent isolates of ancestral D614G virus (D614G.1, D614G.2, see Materials and methods), Omicron subvariants BA.1 and BA.5, the D6 and D190 isolates, and isolates from study visits at day 20 (D20), 34 (D34), and 106 (D106) post-diagnosis. We then calculated the ratio of the fraction of dead infected cells to total infected cells at 24 hours post-infection. We observed that both isolates of ancestral virus caused cell death to a similar extent, with about 11-13% of infected cells being detected as dead. The frequency of death was about 3-fold lower for BA.1, at 4.2%, and was higher for BA.5, at 5.7%. The frequency of D190 induced cell death, at 9.6%, was significantly higher than both BA.1 and BA.5 and more like the ancestral virus isolates (Figure 3B). Among the isolates from the earlier timepoints, D6 infection led to 5.6% dead infected cells, a similar frequency of cell death induction as BA.5. The frequency of cell death increased to 6.8% for D20 and to 9.6% for D34 and D106. The cell death frequency was significantly lower in D6 and D20 relative to D190 infections (Figure 3C).
Because cell-cell fusion depends on cell surface spike expression, we investigated surface spike expression on the infected cells in infection foci (where each focus contains on the order of 10 cells, see27) formed by live virus on a confluent cell monolayer 18 hours post-infection. To ensure that we did not inadvertently permeabilize the cells as part of the staining procedure, thereby including intracellular spike, we also stained for SARS-CoV-2 nucleocapsid protein, not present on the cell surface (Figure 4A). We observed that in unpermeabilized cells, surface spike was readily detected. In contrast, nucleocapsid protein was not detected on the cell surface, although levels were high after permeabilization. We then measured mean intensity, which is independent of focus size, and normalized by the mean intensity of the ancestral SARS-CoV-2 isolate (D614G isolate 1) included in each of the independent experiments. We observed that the two ancestral virus isolates were similar in spike expression. BA.1 infection showed significantly reduced cell surface spike, but levels increased moderately in BA.5 infection. D190 infection had significantly lower spike than the ancestral virus isolates but higher surface spike relative to D6 infection (Figure 4B).
Discussion
Here we investigated the effect of evolutionary changes which happened in SARS-CoV-2 in prolonged infection in an individual immunosuppressed because of advanced HIV disease. We assayed a virus isolated from early in the course of infection at day 6 post-diagnosis, at day 190 post-diagnosis, as well as at the intermediate timepoints. The day 6 isolate already evolved a considerable number of mutations, or alternatively, the infecting virus was mutated. Between day 6 and day 190 post-diagnosis, the infection evolved immune escape from neutralizing antibodies, consistent with what is observed in long-term SARS-CoV-2 infections in individuals with immunosuppression22. However, the ability of the virus to induce cell fusion and cell death upon infection also changed through the course of evolution. While the early viral isolate induced lower cell fusion and death compared to ancestral virus, induction of fusion and cell death increased with evolution
We cannot determine whether the infecting virus evolved a reduced ability to cause cell fusion and death by day 6 post-diagnosis, or whether this was a property of the infecting virus. Symptom onset for the infection was reported to be 10 days pre-diagnosis14, which a seems a short time to evolve substantial differences relative to the ancestral virus, but we do not know for certain how long the individual was infected. In addition, SARS-CoV-2 has been documented to evolve relatively rapidly in some cases56. However, we have determined that the subsequently isolated viruses evolved an enhanced ability to fuse cells and cause cell death relative to the early isolate.
Consistent with the notion that highly mutated variants may not all have attenuated pathogenicity, Beta variant infections in South Africa did not show lower pathogenicity relative to ancestral virus infections50. This does not imply that there is no selective advantage in attenuation. For example, variants which do not lead to more severe disease may be more transmissible since the infected individual may continue daily activities and therefore transmit to more people.
The Omicron subvariants circulating since November 2021 have led to infections with a lower probability of developing severe disease1,2,57,58. Moreover, Omicron infection boosts the immune response against other variants as well as against other Omicron subvariants59, likely maintaining protective immunity through infection. The results presented here, with the limitations that we tested cellular parameters as measures of pathogenicity and viruses from only one case of long-term evolution, may indicate that SARS-CoV-2 evolution in long-term infection does not have to result in attenuation. It may indicate that a future variant could be more pathogenic than currently circulating Omicron strains.
Materials and methods
Informed consent and ethical statement
Swabs for the isolation of ancestral (D614G.1 and D614G.2), Beta, Delta, D6, D20, D34, D106, and D190 viruses and blood samples to test virus neutralization were obtained after written informed consent from adults with PCR-confirmed SARS-CoV-2 infection who were enrolled in a prospective cohort study at the Africa Health Research Institute approved by the Biomedical Research Ethics Committee at the University of KwaZulu–Natal (reference BREC/00001275/2020). The Omicron/BA.1 was isolated from a residual swab sample with SARS-CoV-2 isolation from the sample approved by the University of the Witwatersrand Human Research Ethics Committee (HREC) (ref. M210752). The sample to isolate Omicron/BA.5 was collected after written informed consent as part of the COVID-19 transmission and natural history in KwaZulu-Natal, South Africa: Epidemiological Investigation to Guide Prevention and Clinical Care in the Centre for the AIDS Programme of Research in South Africa (CAPRISA) study and approved by the Biomedical Research Ethics Committee at the University of KwaZulu–Natal (reference BREC/00001195/2020, BREC/00003106/2021).
Reagent availability statement
Isolates and raw image files are available upon reasonable request. Sequences of isolated SARS-CoV-2 used in this study have been deposited in GISAID with accession:
Whole-genome sequencing, genome assembly and phylogenetic analysis
RNA was extracted on an automated Chemagic 360 instrument, using the CMG-1049 kit (Perkin Elmer, Hamburg, Germany). The RNA was stored at −80?C prior to use. Libraries for whole genome sequencing were prepared using either the Oxford Nanopore Midnight protocol with Rapid Barcoding or the Illumina COVIDseq Assay. For the Illumina COVIDseq assay, the libraries were prepared according to the manufacturer’s protocol. Briefly, amplicons were tagmented, followed by indexing using the Nextera UD Indexes Set A. Sequencing libraries were pooled, normalized to 4 nM and denatured with 0.2 N sodium acetate. An 8 pM sample library was spiked with 1% PhiX (PhiX Control v3 adaptor-ligated library used as a control). We sequenced libraries on a 500-cycle v2 MiSeq Reagent Kit on the Illumina MiSeq instrument (Illumina). On the Illumina NextSeq 550 instrument, sequencing was performed using the Illumina COVIDSeq protocol (Illumina Inc, USA), an amplicon-based next-generation sequencing approach. The first strand synthesis was carried using random hexamers primers from Illumina and the synthesized cDNA underwent two separate multiplex PCR reactions. The pooled PCR amplified products were processed for tagmentation and adapter ligation using IDT for Illumina Nextera UD Indexes. Further enrichment and cleanup was performed as per protocols provided by the manufacturer (Illumina Inc). Pooled samples were quantified using Qubit 3.0 or 4.0 fluorometer (Invitrogen Inc.) using the Qubit dsDNA High Sensitivity assay according to manufacturer’s instructions. The fragment sizes were analyzed using TapeStation 4200 (Invitrogen). The pooled libraries were further normalized to 4nM concentration and 25 μL of each normalized pool containing unique index adapter sets were combined in a new tube. The final library pool was denatured and neutralized with 0.2N sodium hydroxide and 200 mM Tris-HCL (pH7), respectively. 1.5 pM sample library was spiked with 2% PhiX. Libraries were loaded onto a 300-cycle NextSeq 500/550 HighOutput Kit v2 and run on the Illumina NextSeq 550 instrument (Illumina, San Diego, CA, USA). For Oxford Nanopore sequencing, the Midnight primer kit was used as described by Freed and Silander55. cDNA synthesis was performed on the extracted RNA using LunaScript RT mastermix (New England BioLabs) followed by gene-specific multiplex PCR using the Midnight Primer pools which produce 1200bp amplicons which overlap to cover the 30-kb SARS-CoV-2 genome. Amplicons from each pool were pooled and used neat for barcoding with the Oxford Nanopore Rapid Barcoding kit as per the manufacturer’s protocol. Barcoded samples were pooled and bead-purified. After the bead clean-up, the library was loaded on a prepared R9.4.1 flow-cell. A GridION X5 or MinION sequencing run was initiated using MinKNOW software with the base-call setting switched off. We assembled paired-end and nanopore.fastq reads using Genome Detective 1.132 (https://www.genomedetective.com) which was updated for the accurate assembly and variant calling of tiled primer amplicon Illumina or Oxford Nanopore reads, and the Coronavirus Typing Tool56. For Illumina assembly, GATK HaploTypeCaller -- min-pruning 0 argument was added to increase mutation calling sensitivity near sequencing gaps. For Nanopore, low coverage regions with poor alignment quality (<85% variant homogeneity) near sequencing/amplicon ends were masked to be robust against primer drop-out experienced in the Spike gene, and the sensitivity for detecting short inserts using a region-local global alignment of reads, was increased. In addition, we also used the wf_artic (ARTIC SARS-CoV-2) pipeline as built using the nextflow workflow framework57. In some instances, mutations were confirmed visually with .bam files using Geneious software V2020.1.2 (Biomatters). The reference genome used throughout the assembly process was NC_045512.2 (numbering equivalent to MN908947.3). For lineage classification, we used the widespread dynamic lineage classification method from the ‘Phylogenetic Assignment of Named Global Outbreak Lineages’ (PANGOLIN) software suite (https://github.com/hCoV-2019/pangolin).
Cells
The H1299-E3 (H1299-ACE2, clone E3, H1299 originally from ATCC as CRL-5803) cell line was derived from H1299 as described in our previous work52,60. The H1299-E3 cells were propagated in growth medium consisting of complete Roswell Park Memorial Institute (RPMI) 1640 medium with 10% fetal bovine serum (Hyclone) containing 10mM of hydroxyethylpiperazine ethanesulfonic acid (HEPES), 1mM sodium pyruvate, 2mM L-glutamine and 0.1mM nonessential amino acids (all Sigma-Aldrich). Cells were passaged every second day. For virus isolation, Vero E6 cells (originally ATCC CRL-1586, obtained from Cellonex in South Africa) were propagated in complete growth medium consisting of Dulbecco’s Modified Eagle Medium (DMEM) with 10% fetal bovine serum (Hyclone) containing 10mM of HEPES, 1mM sodium pyruvate, 2mM L-glutamine and 0.1mM nonessential amino acids (all Sigma-Aldrich). Vero E6 cells were passaged every 3–4 days.
Virus isolation
All work with live virus was performed in Biosafety Level 3 containment using protocols for SARS-CoV-2 approved by the Africa Health Research Institute Biosafety Committee. ACE2-expressing H1299-E3 cells were seeded at 4.5 × 105 cells in a 6 well plate well and incubated for 18–20 h. After one Dulbecco’s phosphate-buffered saline (DPBS) wash, the sub-confluent cell monolayer was inoculated with 500 μL universal transport medium from swabs diluted 1:1 with growth medium filtered through a 0.45-μm filter. Cells were incubated for 1 h. Wells were then filled with 3 mL complete growth medium. After 4 days of infection (completion of passage 1 (P1)), cells were trypsinized (Sigma-Aldrich), centrifuged at 300 rcf for 3 min and resuspended in 4 mL growth medium. Then all infected cells were added to Vero E6 cells that had been seeded at 1.5 × 105 cells per mL, 20 mL total, 18–20 h earlier in a T75 flask for cell-to-cell infection. The coculture of ACE2-expressing H1299-E3 and Vero E6 cells was incubated for 1 h and the flask was filled with 20 mL of complete growth medium and incubated for 4 days. The viral supernatant from this culture (passage 2 (P2) stock) was used for experiments.
Time-lapse microscopy
6-well glass bottom plates (MatTek) were coated with 300 μL of 0.001% fibronectin (Sigma-Aldrich) in DPBS-/- (Gibco), incubated for 90 mins, then washed 3x with DPBS. H1299-E3 cells were then immediately plated at 60,000 cells per coated well. The next day the cells were infected at 1000 focus-forming units in 1 mL growth media per well. Cell–virus mixtures were incubated for 1 h at 37 °C, 5% CO2 then an additional 1 mL of growth media was added. Infections were imaged using a Metamorph controlled Nikon TiE motorized microscope (Nikon Corporation) in a Biosafety Level 3 Facility with a 20x, 0.75 NA phase objective. Images were captured using an 888 EMCCD camera (Andor). Temperature (37°C), humidity and CO2 (5%) were controlled using an environmental chamber (OKO Labs). Excitation source was 488 laser line and emission was detected through a Semrock Brightline quad band 440–40 /521-21/607-34/700-45 nm filter. For each well, 12 randomly chosen fields of view were captured every 10 minutes.
Image analysis
Timelapse microscopy images were analysed using custom MATLAB v.2019b (MathWorks) script and using the MATLAB Image Analysis Toolbox. For each frame in the movie, both the transmitted light and fluorescent images were used. A coarse segmentation was first performed using the transmitted light image of the cells. Images were flatfield corrected and contrast was enhanced by setting the top and bottom 1% of all pixel intensities to 1 and 0 respectively. The built-in function “rangefilt” was used to determine areas of high contrast (where high pixel intensities were immediately adjacent to low pixel intensities) in the image which corresponded to cell borders. Processed images were also median filtered and holes to filled within segmented objects. The image was then thresholded to remove background signal and obtain a mask of areas occupied by cells. The mask generated from the transmitted light image was then used to remove objects that were not in areas occupied by cells. In the fluorescence images corresponding to the transmitted light images. Fluorescence images were then processed using flatfield correction, contrast enhancement and median filtering. Fluorescent cell nuclei in the YFP channel were used to generate a binary mask for contiguous objects in each image after thresholding. Each object was categorized as multi-nucleate or uni-nucleate based on pixel area, where the threshold for a single nucleus was calculated as the mean area of objects/nuclei in the uninfected condition, at 12 hours post-movie start, + 3 standard deviations of the mean. The fraction of cells in fusions was calculated by dividing the total pixel area of objects above the single nucleus threshold by the total pixel area occupied by nuclei in the same frame. The number of nuclei was calculated by dividing the total pixel area occupied by nuclei by the mean pixel area of one nucleus in the uninfected condition at 12 hours post-movie start.
Detection of infected cell death
H1299-E3 cells were plated at 60,000 cells per well in 6-well plates (Corning) 1 day pre-infection. The next day the cells were infected at 1000 focus-forming units in 1 mL growth media per well. Cell–virus mixtures were incubated for 1 h at 37 °C, 5% CO2 then an additional 1 mL of growth media was added. 24 hours post-infection, cells were trypsinised (Sigma-Aldrich), collected and stained with Blue Live/Dead stain as per manufacturer instructions (L34961, ThermosScientific). The samples were then washed in 1 mL PBS-/- and resuspended in Cytofix/Cytoperm (BD Biosciences) for 20 min at 4°C in the dark. The samples were then stained with 0.5 μg/mL anti-SARS-CoV-2 nucleocapsid-PE (ab283244, Abcam) for 1 hour at 4°C in the dark. Cells were analysed on an Aria Fusion (BD). Data was analysed using FlowJo and Graphpad Prism 9.4.1 software.
Staining for cell-surface and total spike in a focus forming assay
H1299-E3 cells were plated in a 96-well plate (Corning) at 20,000 cells per well 1 day pre-infection. Virus stocks were used at 100 focus-forming units per microwell. Cells were infected with 100 μL of the virus for 1 h, then 100 μL of a 1X RPMI 1640 (Sigma-Aldrich, R6504), 1.5% carboxymethylcellulose (Sigma-Aldrich, C4888) overlay was added without removing the inoculum. Cells were fixed at 18 hours post-infection using 4% methanol-free formaldehyde (ThermoScientific) for 20 minutes. For staining of foci, a rabbit anti-SARS-CoV-2 spike monoclonal antibody (mAb BS-R2B12, GenScript A02058) at 0.5 ug/mL or a rabbit anti-SARS-CoV-2 nucleocapsid monoclonal (ab271180 Abcam) at 1 ug/mL were used as the primary detection antibody. Antibody was resuspended in either a permeabilization buffer containing 0.1% saponin (Sigma-Aldrich), 0.1% BSA (Sigma-Aldrich), and 0.05% Tween-20 (Sigma-Aldrich) in PBS+/+ or a non-permeabilization buffer containing 0.1% BSA and 0.05% Tween-20 in PBS+/+. Plates were incubated with primary antibody at room temperature for 2 hr with shaking, then washed with wash buffer containing 0.05% tween in PBS+/+. Secondary goat anti-rabbit horseradish peroxidase (Abcam ab205718) was added at 1 ug/mL in either permeabilization or non-permeabilization buffers as described above and incubated for 2 hours at room temperature with shaking. TrueBlue peroxidase substrate (SeraCare 5510-0030) was then added at 50uL per well and incubated for 15 minutes at room temperature. Plates were washed in distilled water and then dried for 2 hours and imaged in an ImmunoSpot Ultra-V S6-02-6140 Analyzer ELISPOT instrument with BioSpot Professional built-in image analysis (C.T.L). Data was analyzed using Graphpad Prism 9.4.1.
Statistics and fitting
Fitting was performed using MATLAB v.2019b. Neutralization data were fit to: Here Tx is the number of foci normalized to the number of foci in the absence of plasma on the same plate at dilution D and ID50 is the plasma dilution giving 50% neutralization. FRNT50 = 1/ID50. Values of FRNT50 <1 are set to 1 (undiluted), the lowest measurable value. We note that the most concentrated plasma dilution was 1:25 and therefore FRNT50 < 25 were extrapolated. Fold-change was calculated by dividing the FRNT50 obtained for the homologous virus (the virus which elicited the immunity, e.g. D614G) by the FRNT50 heterologous virus (e.g. D190) per participant, then calculating the geometric mean and 95% confidence intervals over all participant values. The 95% confidence intervals on the median in the time-lapse microscopy data was calculated by first ranking the values in ascending order, then finding the ranks of the lower and upper confidence interval by: Here LR is lower 95% interval index in the ranked vector, UR is the upper 95% interval index in the ranked vector. Values at indexes LR and UR were the lower 95% and upper 95% confidence intervals, respectively. Other statistical tests, measures of central tendency and confidence intervals were performed in GraphPad Prism version 9.4.1.
Data Availability
Isolates and raw image files are available upon reasonable request. Sequences of isolated SARS-CoV-2 used in this study have been deposited in GISAID.
Competing interest statement
AS received an honorarium for a talk given to Pfizer employees.
Acknowledgements
This study was supported by the Bill and Melinda Gates award INV-018944 (AS), National Institutes of Health award R01 AI138546 (AS), and South African Medical Research Council Award 6084COAP2020 (AS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.