Abstract
The recent emergence of SARS-CoV-2 variants showing increased transmissibility and immune escape is a matter of global concern. Their origin remains unclear, but intra-host virus evolution during persistent infections could be a contributing factor. Here, we studied the long-term SARS-CoV-2 infection in an immunosuppressed organ transplant recipient. Frequent respiratory specimens were tested for variant viral genomes by RT-qPCR, next-generation sequencing (NGS), and virus isolation. Late in infection, several virus variants emerged which escaped neutralization by COVID-19 convalescent and vaccine-induced antisera and had acquired genome mutations similar to those found in variants of concern first identified in UK, South Africa, and Brazil. Importantly, infection of susceptible hACE2-transgenic mice with one of the patient’s escape variants elicited protective immunity against re-infection with either the parental virus, the escape variant or the South African variant of concern, demonstrating broad immune control. Upon lowering immunosuppressive treatment, the patient generated spike-specific neutralizing antibodies and resolved the infection. Our results indicate that immunocompromised patients are an alarming source of potentially harmful SARS-CoV-2 variants and open up new avenues for the updating of COVID-19 vaccines.
Introduction
Individuals infected with SARS-CoV-2 develop neutralizing spike-specific antibodies that persist for months and protect against reinfection1. Similarly, neutralizing antibodies generated after vaccination efficiently protect from COVID-192. However, the recent emergence of SARS-CoV-2 variants in the UK (B.1.1.7), South Africa (B.1.351) and Brazil (P.1)3-5 pose a global threat due to their increased transmissibility and resistance to neutralizing antibodies2. The origin of these virus variants remains unclear, but long-term-infected immunocompromised individuals are a likely source, allowing prolonged viral replication and unhindered adaption to the host6,7.
In Germany, the SARS-CoV-2 epidemic started with local outbreaks in February 2020. The city of Freiburg at the border to France and Switzerland was a hotspot due to multiple unrecognized infection events in March 2020. Therefore, immunocompromised patients were closely monitored, as these individuals were expected to have an increased risk of developing severe COVID-19 illness and to suffer from long-term persistent infection with prolonged viral shedding.
Here, we describe an immunosuppressed organ transplant recipient who acquired SARS-CoV-2 during the early phase of the COVID-19 pandemic. The patient had mild respiratory symptoms and was tested positive for SARS-CoV-2 for over 145 days. During this long period, virus variants with multiple amino acid substitutions and deletions in the spike protein evolved that showed increased resistance to neutralizing antibodies, suggesting partial immune escape. Interestingly, however, one of the late virus variants isolated from the patient elicited a broad protective immune response in experimentally infected mice, suggesting that convalescent individuals might become resistant against reinfection by emerging variants of concern.
Results
Clinical presentation of the organ transplant patient persistently infected with SARS-CoV-2
A male person in his 50’s was admitted in spring 2020 to the University Medical Center, Freiburg, for transplantation. The patient was treated with a cocktail of tacrolimus, mycophenolate and prednisone for 7 month (Fig. 1a and b).
Temporal overview of (a) hospitalization, (b) immunosuppressive treatment (daily dose in mg/day) and (c) antiviral therapy (daily dose in mg/day), Remdesivir was given 200 mg on the first day and 100 mg/day 2 to 10. (d) Diagnostic SARS-CoV-2 RT-qPCR cycle threshold (Ct) values of oropharyngeal swabs over time. Day 0 indicates the first positive RT-qPCR result, 12 days after kidney transplantation. The dotted line indicates the cutoff value (Ct ≥ 40) between positive and negative results. (e) Attempts of virus isolation from oropharyngeal swabs. (f) Detection of spike S1-subunit- and nucleoprotein (N) specific antibodies by ELISA. The dotted line indicates the anti-S1 ELISA cutoff at 1.1 arbitrary unit (AU).
The patient was tested positive for SARS-CoV-2 by reverse transcription quantitative polymerase chain reaction (RT-qPCR) which is indicated with day 0(Fig. 1d). The source of infection remained unknown but strict infection prevention measures were initiated. No respiratory symptoms at the time of diagnosis were observed but a CT scan showed mild ground-glass opacities and discrete bilateral pleural effusions on day 4. The patient remained SARS-CoV-2 positive in the following weeks and was therefore kept in isolation. Two month after transplantation the patient suffered from a urinary tract infection with E. coli, requiring antibiotic treatment. Furthermore, he was treated for 5-days (day 56 to 60) with Ivermectin (33 mg/day) (Fig. 1c), a broad-spectrum drug with anti-viral activity in cell culture against several viruses including SARS-CoV-28,9. While the bacterial urinary tract infection was controlled, the infection with SARS-CoV-2 was not.
RT-qPCR positive swab samples were used to successfully isolate virus on VeroE6 cells confirming shedding of infectious SARS-CoV-2 (Fig. 1e). Due to his critical condition, the patient stayed in the hospital until day 72 when he was discharged for home quarantine. However, he was re-hospitalized at day 106 to 126. Afterwards, the immunosuppressive regimen was modified by withdrawing mycophenolate mofetil treatment and by increasing the dose of prednisone (day 126) to allow for a better antiviral adaptive immune response. After re-admission at day 140 for control purposes, the patient was still RT-qPCR positive. As an attempt to control the infection, he was treated for 10 days with Remdesivir (200 mg on day 140, then 100 mg/daily, day 141 - 149), a nucleoside analog with anti-SARS-CoV-2 activity in vitro10 and in vivo7,11. Subsequently, negative RT-qPCR tests until day 189 and failed virus isolation attempts suggested that the infection had resolved (Fig. 1d and e) (Extended Table 1).
Nucleoprotein (N)-specific antibodies were detected already 12 days after the first positive qPCR result and afterwards surged rapidly (Fig. 1f). In contrast, IgG antibody levels specific for SARS-CoV-2 spike protein (S1) determined by ELISA constantly oscillated around the cutoff value between days 40 and 123. Only, when the patient was hospitalized at day 140, high levels of S1-specific IgG were detected and persisted at least until day 175 (Fig. 1f). Concomitantly with the increased spike-specific antibodies and the onset of Remdesivir treatment (day 140), RT-qPCR analysis showed steadily increasing Ct values, indicating diminishing virus replication (Fig. 1d and Extended Table 1). In summary, during the 25 weeks of infection with SARS-CoV-2 the patient had no severe respiratory or systemic symptoms and was finally able to clear the virus, likely due to the development of neutralizing antibodies and possibly due to the inclusion of antiviral Remdesivir treatment.
The table specifies the days of hospitalization post first detection of SARS-CoV-2 infection, SARS-CoV-2 specific RT-qPCR Ct values from throat swabs, anti-S-specific ELISA (cutoff value of 0.8 to 1.1 AU), anti-N-specific ELISA (cutoff value of 24 AU), attempts for virus isolation from throat swabs on VeroE6 cells, and titers of neutralizing antibodies against the prototype Muc-IMB-1 reference strain 20.
Genetic relationship of the patient’s SARS-CoV-2 variants with circulating strains
Full-length SARS-CoV-2 genome sequences were obtained from oropharyngeal swabs collected between day 0 and 140 and phylogenetic trees were constructed including sequences representative for the Freiburg area in spring 2020 (Fig. 2a) or a set of randomly selected GISAID sequences of isolates collected in Germany in spring 2020(Fig. 2b). The viral genomic sequences obtained from the patient clearly clustered to those of strains circulating in spring in Germany and in the Freiburg area (Fig. 2a and b). The phylogenetic analysis also demonstrated a close relationship with sequences obtained from two patients of the Medical Center Freiburg (Fig. 2a). However, no clear epidemiologic link was found between the immunosuppressed and the other two patients.
Phylogenetic analysis of the viral sequences obtained from patient swabs between day 0 to day 140, after the first positive RT-qPCR result. The sequences were aligned to a set of representative SARS-CoV-2 genome sequences from the Freiburg area (a) and from Germany (b) that were obtained in spring 2020 and have been deposited in the GISAID data bank (Extended data table 2 and 3). The circularized maximum-likelihood phylogenetic tree was constructed with IQ-Tree (GTR+F+I) and rooted on the Wuhan-Hu-1 reference sequence (NC_045512). The sequences obtained from the immunosuppressed patient are indicated as red dots. The scale represents nucleotide substitutions per site. (c) Schematic overview of the viral genome variations from patient swab samples (day 0-140) in comparison to the Wuhan-Hu-1 reference sequence. The heatmap summarizes the positions in the viral genome and the variant frequencies in the different samples (cut off values of 25% and 10% for the S gene analysis, respectively). The days of sampling are indicated at the right and the heatmap color intensity indicates variant frequencies. Stars denote non-synonymous mutations leading to amino acid substitutions in the spike protein (> 50 % of reads). (d) Schematic overview of the SARS-CoV-2 spike protein including the S1 and S2 cleavage products and functional domains such as the N-terminal domain (NTD), receptor-binding domain (RBD), receptor binding motif (RBM), S1/S2 proteolytic furin cleavage site, fusion peptide (FP), heptad repeat regions (HR1/HR2), transmembrane domain (TM) and C-terminal domain (CT). Selected non-synonymous changes in the spike (S) gene from panel c are indicated. (e) Summary of mutations found in the spike protein of the patient sequences obtained on d14, d105 and d140 (>50 % of reads) in comparison to circulating new variants of concern from the UK, B.1.1.73, South Africa, B.1.3514 and Brazil, P.15.
Sequence accession numbers of full-length SARS-CoV-2 genome sequences generated from swabs of patients at the University Medical Center, Freiburg, Germany between February and April and of the immunosuppressed patient.
We gratefully acknowledge the following authors from the originating laboratories responsible for obtaining the specimens, as well as the submitting laboratories where the genome data were generated and shared via GISAID, on which this research is based. All submitters of data may be contacted directly via www.gisaid.org. Authors are sorted alphabetically.
Analysis of the viral RNA that has been extracted from the patient samples revealed several nucleotide substitutions in ORF1ab, the spike gene, ORF3a, M and N genes in comparison to the Wuhan-Hu-1 reference genome (Fig. 2c). The continuous presence of nine common mutations in all sequences argues against a possible reinfection but is compatible with viral persistence. Within the first two weeks, no changes in the viral genome were observed, while from day 42 onward acquisition of several mutations occurred. Apart from low frequency mutations, some mutations accumulated over time indicative for the selection of distinct variants. The mutation 23403G resulting in the D614G substitution marks the genotype that now dominates worldwide12.
The most remarkable changes found in the S gene, which were confirmed by Sanger sequencing of cDNA clones, were in-frame deletions and non-synonymous substitutions in the N-terminal domain (NTD) and the receptor binding motif (RBM) (Fig. 2d), respectively. Interestingly, the two amino acid deletions in the NTD were associated with specific single amino acid substitutions in the RBM: del141-144 with F490L and del244-247 with E484G. Both deletions, which are located in two adjacent flexible loops of the NTD (Fig. 3), might affect the conformation of this subdomain and are targets of neutralizing antibodies13,14. The two amino acid substitutions (E484G and F490L) may lead to subtle conformational changes in the RBD, which is the main target of neutralizing antibodies and a known hotspot for mutations conferring escape from neutralizing antibodies15-17. Thus, the SARS-CoV-2 spike variants from the immunosuppressed patient share mutations with the escape variants of concern from the UK, South Africa and Brazil (Fig. 2e).
The spike structure (PDB accession number: 7BNM) with the most prominent mutations in the patient viral sequences is shown in surface presentation. The NTD is colored in blue, the RBD in red and the fusion peptide in gray. Close-ups of the single NTD and RBD regions defined by boxes are presented as ribbons. The location of the deletions in the NTD and amino acid substitutions in the RBD are indicated by black residues. Furthermore, the deletions in the NTD are displayed as amino acid alignments at the right.
Prolonged viral persistence resulted in impaired viral fitness
In the early phase of viral persistence (days 0 to 34) virus isolation was successful, indicating constant virus shedding (Fig. 1e). We repeatedly failed, however, to isolate virus thereafter when the Ct values increased above 25. On day 105, the Ct value dropped to 23 and virus isolation was again successful (Fig. 1e). The sequences of the two distinct virus isolates obtained at day 14 (d14) and day 105 (d105) were compared to the corresponding sequences obtained from swab samples of the same day. While the sequences of the d14 isolate and the d14 swab were identical, the sequence of the d105 isolate partially differed from the d105 swab sequence (Fig. 4a). The viral genome sequences of the d105 isolate and the swab samples contained the amino acid deletion 244-247 combined with the E484G mutation, while the amino acid deletion del141-144 and the F490L substitution were only found in the swab samples. Since the majority of the mutations in the d105 isolate was also detected in the swab samples, the d105 virus might represent an abundant genotype that was selected during persistence in the patient. However, we cannot exclude that some alterations were introduced during the process of virus isolation in cell culture including a frame shift deletion in ORF7b and a deletion of 21 nucleotides (del23601-23621) in about 50% of the reads that mapped to the spike gene and resulted in the replacement of the eight amino acids SPRRARSV by a single leucine upstream of the S1/S2 cleavage site (Fig. 4a).
(a) Schematic overview of the sequence variations in SARS-CoV-2 genomes detected in early (d14) or late (d105) swab samples and isolated viruses. The heatmap illustrates the positions and the frequency of major variations in the viral genome (cut off 15%). The days of isolation are indicated at the right. The heatmap colors represent the variant frequencies. In ORF7b, L14* indicates a frame shift mutation due to a deletion of two nucleotides. (b) Immunofluorescence analysis of SARS-CoV-2 infected cell cultures. VeroE6 cells were infected with one of the two virus isolates (d14 or d105) using a multiplicity of infection (moi) of 0.1 plaque forming units (pfu)/cell. At 8 h post infection, the cells were fixed and stained with SARS-CoV-2 N- and S-specific antibodies (red). In addition, F-actin (white) and nuclear DNA (DAPI, blue) were detected. (c) Western blot analysis of viral protein expression. Calu-3 cells were infected with a moi of 0.001 pfu/cell with either of the two virus isolates (d14 or d105). Cells were lysed 8h, 24h, 48h and 72h post-infection and analyzed using N- and S-specific antibodies. Detection of β-actin was used as loading control. (Fig.4 continued on the following page) (d) Growth of the two patient isolates in VeroE6 and Calu-3 cells. The cells were infected with either of the two patient isolates (d14 or d105) using a moi of 0.001 pfu/cell. At different time points post infection, cell culture supernatants were collected and viral titers were determined by plaque assay on VeroE6 cells. The log-transformed titers are shown as means ± SD of results from three independent experiments. Significance was determined via two-way ANOVA with a Sidak’s multiple comparison test, **p<0.0021, ***p<0.0002, ns=non significant. (E and F) Late SARS-CoV-2 isolate is attenuated in mice. Weight loss (e) and survival (f) of 8 to 12 weeks-old K18-hACE2 mice intranasally infected with 200 or 2000 pfu of d14 or d105 viruses. Signs of disease and body weight loss were monitored daily for 14 days. The weight loss is visualized as mean ± SEM. Significance was determined via two-way ANOVA with a Sidak’s multiple comparison test, *p<0.0332, **p<0.0021, ***p<0.0002.
The d14 and d105 isolates both showed accumulation of viral N- and spike proteins in infected cells by indirect immunofluorescence and Western blot analyses (Fig. 4b and c). In VeroE6 cells, both isolates grew to comparable titers, whereas replication of the d105 isolate was impaired in human lung adenocarcinoma Calu-3 cells (Fig. 4d). To confirm the attenuated phenotype of d105 in vivo, K18-hACE2 mice encoding the human SARS-CoV-2 receptor, angiotensin-converting enzyme type 2 (ACE2)18 were intranasally infected with two different doses of the d14 and d105 isolates (200 or 2000 plaque-forming units (pfu)). Animals infected with the d14 isolate showed a significant weight loss within 4 to 8 days. One out of eight animals infected with 200 pfu survived the infection, whereas all animals infected with 2000 pfu had to be euthanized due to severe disease symptoms or weight loss (Fig. 4e and f). In contrast, K18-hACE2 mice infected with the d105 isolate showed a transient weight loss but only two out of seven animals infected with 2,000 pfu got severely sick and had to be sacrificed (Fig. 4e and f), highlighting the attenuated phenotype of this isolate. Together these findings suggest that the mutations in the viral genome of the d105 isolate caused reduced viral fitness.
SARS-CoV-2 escape variants emerged during viral persistence
The amino acid deletions and substitutions in the spike proteins of the emerging viral variants could have been selected by the antiviral immune response of the host. To address this issue, serum samples of the patient were tested for neutralizing antibodies in a plaque reduction assay performed with either the d14 or the d105 virus isolate. Only sera collected from the patient after day 123 showed SARS-CoV-2 neutralizing activity (Fig. 5a), which coincided with the increase of S-specific IgG between day 123 and 140 (Fig. 1f). Intriguingly, while the d14 isolate was efficiently neutralized up to serum dilutions of 1:128, the d105 virus was poorly inhibited even at the lowest serum dilutions used (1:32) (Fig. 5a), suggesting that the substitutions in the d105 spike protein caused escape from neutralizing antibodies. The neutralization titers detected in sera from the immunosuppressed patient were generally lower than those detected in serum from an immunocompetent convalescent COVID-19 patient (Fig. 5a, positive control). Similarly, antisera of convalescent COVID-19 patients and of individuals previously vaccinated with the BNT162b2 mRNA vaccine (Pfizer/BioNTech) showed higher neutralizing activity against the d14 isolate as compared to the d105 isolate (Fig. 5b).
(a and b) Detection of neutralizing activity of immune sera against SARS-CoV-2 variants. 100 pfu of the d14 and d105 isolates were incubated for 60 min at room temperature with serial dilutions of the patient sera. Sera obtained from naïve (– ctrl) or convalescent individuals (+ ctrl) served as negative and positive controls. Virus neutralization was determined by plaque assay on VeroE6 cells. Virus titers are indicated as percentages (mean ±SD) of the titer of the untreated virus inoculum. The dotted lines indicate the cutoff value between positive (<50%) and negative (>50%) neutralization. (a) Sera from the immunocompromised patient. The times of blood withdrawal are indicated. (b) Convalescent sera from COVID-19 patients suffering from mild, moderate or severe disease or human post vaccination (BNT162b2 mRNA) sera. (c-e) Neutralization capacity of SARS-CoV-2 antisera using VSV*ΔG(FLuc) vector pseudotyped with the SARS-CoV-2 spike protein and coding for firefly luciferase. The pseudotyped viruses were incubated with serial dilutions of a COVID-19 convalescent serum prior to inoculation of VeroE6 cells. Pseudotyped virus infection was monitored 16 h post infection by measuring the firefly luciferase activity in the cell lysates. (c) Neutralization of VSV*ΔG(FLuc) pseudotyped with the early and late SARS-CoV-2 spike variants (d14 and d105) using serial dilutions of a COVID-19 convalescent serum. (d and e) Neutralization of VSV*ΔG(FLuc) pseudotyped with the d14 spike protein containing the individual or combined mutations characteristically found in the late d105 and d140 variants. Immune serum from a convalescent COVID-19 patient (d) or a vaccinated person (e) were analyzed. The neutralization was determined by calculating the NT50 via a non-linear regression (variable slope, four parameters). Shown are means ± SD (n=3). Statistics were calculated with a one-way ANOVA (Tukey’s multiple comparison test), ns = non-significant, *p<0.0332, **p<0.0021, ***p<0.0002.
To independently determine the neutralization titer 50 (NT50) of convalescent sera against the various spike variants found in the patient, we made use of a virus pseudotype system based on single-cycle vesicular stomatitis virus (VSV) vector encoding a firefly luciferase reporter, VSV*ΔG(FLuc). High NT50 values were observed with the d14 spike protein (Fig. 5c). In contrast, pseudotype virus bearing the d105 spike protein was neutralized with significantly lower efficacy as indicated by a 8.3-fold reduced NT50 value (Fig. 5c).
Next, we analyzed pseudotype virus harboring the d14 spike protein with single or combined mutations including del244-248 and E484G (for d105) and del141-144 and F490L (for d140) (Fig. 2d). Compared to pseudotype virus displaying the parental d14 spike protein, viruses pseudotyped with either the E484G or the F490L mutant spike proteins were equally well neutralized by a COVID-19 convalescent serum, whereas pseudotype virus bearing spike proteins with the amino acid deletions del141-144 or del244-247 were significantly less well neutralized (Fig. 5d). The combination of the amino acid deletions del141-144 or del244-247 with either E484G or F490L did not further reduce neutralization efficacy. A different pattern was observed when immune serum from a person who had been immunized with an mRNA-based SARS-CoV-2 vaccine was analyzed (Fig. 5e). Using this immune serum, virus pseudotyped with the E484G mutant spike was less well neutralized than virus bearing the parental spike protein. Furthermore, pseudotype virus displaying spike protein containing both the E484G substitution and the del141-144 deletion, was less efficiently neutralized than pseudotype viruses containing either E484G or del141-144, suggesting that the two mutations acted in a synergistic manner. In summary, the amino acid changes del141-144 and del244-247, both located in the NTD, and E484G in the RBD all affect crucial antigenic regions19 which were selected during viral persistence as they allow escape of SARS-CoV-2 from the humoral immune response.
SARS-CoV-2 escape variant d105 induces broad protective immunity in vivo
The reduced neutralization capacity of the patient’s sera against the d105 isolate raised the question whether the changes in the spike protein might have compromised the induction of an efficient antiviral immune response. To address this question, sera from K18-hACE2 mice that survived the infection with the d14 or the d105 isolates (Fig. 4e and f) were collected 21 days post infection or later. We also prepared sera from K18-hACE2 mice surviving infection with the Muc-IMB-1 isolate, lineage B.1, encoding an identical spike protein sequence like the d14 isolate20, in order to increase the number of convalescent sera of animals surviving infection with wild-type SARS-CoV-2. The sera of the animals that survived infection with the wild-type viruses had about 2-fold higher levels of SARS-CoV-2-specific IgG antibodies (assessed by an immunofluorescence-based assay), than the sera of mice that survived the milder d105 infection (Fig. 6a).
Sera were collected from convalescent K18-hACE2 mice at least 21 days post infection with Muc-IMB-1 (n=5), d14 (n=2) or d105 (n=13) virus isolates. (a) Anti-SARS-CoV-2 IgG titers of serially diluted sera (mean ±SD) were determined using virus-infected cells and indirect immunofluorescence analysis (IFA). (b and c) Neutralization of d14 and d105 virus isolates by convalescent mouse sera obtained after infection with wild-type SARS-CoV-2, d14 and Muc-IMB-1 (anti-wt sera, n=7) (b), or with variant d105 virus isolate (anti-d105 sera, n=13) (c). Neutralization capacity was determined by incubating 100 pfu of either virus isolate with serial dilutions of the mouse sera. The mixture was then applied to VeroE6 cells and infected cells were stained with N-specific antibodies. The highest dilution for each individual serum causing 50% neutralization was plotted as neutralizing titer. (d and e) Neutralization of B.1.1.7 and B.1.351 variants of concern by mouse convalescent sera was determined as described in panels b and c. 100 pfu of the variant viruses were incubated with dilutions of sera from mice infected with wild-type SARS-CoV-2 (wt) (d) or d105 isolate (e). Neutralization titers are meant as the highest dilution for each individual serum causing 50% reduction of infectivity. Each serum titer (b - e) is shown as mean out of two independent experiments. Significance was determined via a paired t test with *p<0.0332, **p<0.0021, ***p<0.0002. (f and g) Convalescent animals are protected against re-challenge infection. Weight loss (f) and survival (g) of convalescent K18-hACE2 mice (mean ± SEM), challenged one to four months after the prime infection. Animals primarily infected with d14 and Muc-IMB-1 viruses (pooled wt survivors, n=7), or with d105 virus (n=14) were intranasally challenged with 100,000 pfu of d14 or d105 viruses (2 to 7 mice per group, as indicated). As a control, naïve 8 weeks old K18-hACE2 mice were intranasally infected with 100,000 pfu of d14 or d105 isolate viruses (n=2 per group).
The capacity of sera obtained from the d14- and Muc-IMB-1-infected mice to neutralize the d105 isolate was some 4-fold lower than the neutralizing capacity against the d14 isolate (Fig. 6b), indicating partial immune escape of the d105 isolate. Intriguingly, the d105 mouse sera neutralized the d105 virus more efficiently than the d14 isolate (Fig. 6c). Additionally, we tested the neutralization capacity of the convalescent mouse sera against two recent German isolates of the UK (B.1.1.7) and the South Africa (B.1.351) variants of concern. Sera of wild-type-infected mice were more effective in neutralizing the B.1.1.7 virus variant than the B.1.351 variant (Fig. 6d). However, the opposite was observed using sera from d105-infected animals since the B.1.351 variant demonstrated a higher sensitivity to neutralization by the d105 sera than the B.1.1.7 variant (Fig. 6e). Finally, all convalescent animals, including those infected with Muc-IMB-1, were challenged with a lethal dose (100,000 pfu) of either the d14 or the d105 virus isolate one to four months after the first infection. In contrast to the naïve control animals, all convalescent mice survived the challenge infection without any signs of disease or weight loss (Fig. 6f and g), demonstrating that the animals were protected against challenge infection by both virus variants.
SARS-CoV-2 specific CD8+ T cells are not driving emergence of escape variants
Finally, we assessed whether the variations in the spike S1 domain also resulted in an escape from the CD8+ T cell response. We performed in silico prediction of CD8+ T cell epitopes within the SARS-CoV-2 S1 domain restricted by the HLA class I alleles of the immunosuppressed COVID-19 patient (HLA-A*02:01, HLA-A*03:01, HLA-B*51:01, HLA-B*56:01). Using ANN4.0 of the Immune Epitope Database website21 we identified one nonamer peptide within the NTD with good binding properties to HLA-A*02:01 (S133-141 FQFCNDPFL) and another nonamer peptide (S142-150 GVYYHKNNK) with binding to HLA-A*03:01, both overlapping with del141-144. In addition, we identified a decamer peptide (S240-249) with potential binding to HLA-A*02:01 that overlaps with del244-247. Further CD8+ T cell epitope peptides overlapping with the E484G or F490L substitutions could not be predicted. Subsequently, we tested whether the selected peptides represent SARS-CoV-2-specific CD8+ T cell epitopes by incubating PBMCs from the convalescent, immunosuppressed COVID-19 patient on day 233 or from immunocompetent, convalescent, HLA-A*02:01/HLA-A*03:01 positive COVID-19 patients with the peptides for 14 days. However, after stimulation with the selected peptides neither PBMCs from the immunosuppressed nor from the convalescent donors showed any IFNγ-positive CD8+ T cell response whereas a weak IFNγ-positive CD8+ T cell response targeting the non-overlapping epitope HLA-A*03:01/S378-386 was detectable in PBMCs from the immunosuppressed patient (Extended Data Fig. 1). This observation indicates that the predicted SARS-CoV-2 peptides that overlap with the mutated S1 regions do not represent CD8+ T cell epitopes. Hence, the mutations within the S1 regions of the d105 and d140 spike proteins were most probably not selected due to CD8+ T cell escape, but to escape from the humoral response.
Discussion
Circulating SARS-CoV-2 variants typically acquire only few mutations over time which accumulate at a relatively constant rate of about 1-2 mutations per month22. Accordingly, the predominant virus genotypes initially found in swab samples during the first weeks of the persistent infection of the immunocompromised patient were relatively stable and grouped into Nextstrain clade 20B, Pangolin lineage B.1.1, together with simultaneously circulating variants (Fig. 2a and b). However, from day 42 onward, synonymous and non-synonymous mutations accumulated in the viral genomes, including two in-frame amino acid deletions in the NTD (del141-144 and del244-247) as well as two single amino acid exchanges (E484G and F490L) in the RBM. Interestingly, the substitutions in the RBM were exclusively associated with one or the other of the deletions, namely F490L with del141-144 and E484G with del244-247, respectively. This striking interlinked coevolution may have been favored by the necessity to preserve the functionality or stability of the spike protein and to maintain or improve viral fitness. The loops N3 (130-150) and N5 (240-260) that are affected by the two NTD deletions are a preferred target of in-frame deletions of variable length and are therefore referred to as “recurrent deletion regions”14.
Of note, the transient and altering patterns of amino acid changes in the viral spike protein were most likely the result of an ongoing conflict between the persisting virus and the patient’s adaptive immune system. The rather weak spike-specific antibody response between days 40 to 123 presumably selected for spike escape variants, as demonstrated by their reduced sensitivity to inhibition by immune sera obtained from the patient at later time points, or from convalescent COVID-19 patients, as well as from COVID-19 vaccinated individuals. It is conceivable that prolonged viral replication in immunosuppressed patients can lead to the emergence of new immune-escape variants, such as the SARS-CoV-2 variants of concern like UK B.1.1.7, South Africa B.1.351 and Brazil P.1. They all share mutations in the same regions of the spike protein as the escape variants described in this study (Fig. 2e). Accumulation of amino acid substitutions or deletions in similar regions of the spike protein were reported before for isolates of persistently infected, immunosuppressed patients6,7 and also for isolates from patients treated with antibody cocktails and convalescent plasma7,15,23.
The reasons for the late but sudden rise of spike-specific antibodies in the patient serum between day 123 and 141 are not clear. We suspect that the pausing of mycophenolate mofetil from day 126 until day 175 favored broad, spike-specific antibody response that finally terminated the infection. To allow spike-specific antibody production, discontinuation of mycophenolate mofetil treatment of COVID-19-infected transplant recipients is advisable and in line with current clinical guidelines e.g. of the British Transplantation Society (https://bts.org.uk/information-resources/covid-19-information/, updated 22ndJanuary 2021). Adjusting immunosuppressive medications appears to be crucial for induction of an antiviral immune response and clearance of SARS-CoV-224.
Postinfection sera from mice that survived infection with wild-type virus demonstrated reduced neutralizing activity against the late d105 virus isolate, highlighting the antibody escape phenotype of this variant. Conversely, sera from mice previously infected with the d105 isolate more efficiently neutralized the d105 variant than the d14 virus. Hence, broadly neutralizing antibodies that were elicited by new immunogenic epitopes exposed on the mutated d105 spike protein may have controlled the escape variant in the persistently infected patient. To confirm such an extended neutralizing activity of the d105 immune sera, we used recent virus variants of concern and detected enhanced neutralization of the South Africa B.1.351 variant by the d105 antisera when compared to the UK B.1.1.7 variant. These findings match with recent analyses of convalescent plasma samples from patients that recovered from B.1.351 infections. Cele et al. showed efficient neutralization of an early 2020 isolate as well as of the late South African variant of concern by these convalescent antisera 25. Because these globally emerging viruses show a clear escape from vaccine induced humoral immunity2,26,27, our findings might be important for the redesign of future vaccines.
The SARS-CoV-2 d105 variant that was isolated from the patient throat swab at a late stage of his disease showed reduced viral fitness in Calu-3 cell cultures and in K18-hACE2 mice, but not in VeroE6 cells. The reason for attenuation and the causative mutations in the viral genome are presently not known, but the observed viral growth restriction might depend on host cell-intrinsic innate restriction factors. In addition, we detected a 21 nucleotide-long deletion in the multi basic cleavage site of the spike protein connecting the S1 and S2 subdomains. Additional mutations included a truncation of ORF7b, encoding an accessory protein with yet unknown biological function. Deletions in ORF7b have been described in an early Asian cluster of SARS-CoV-2 without a sign of attenuation28 as well as in a SARS-CoV/2003 clone29. As recently shown30, we assume that the loss of the multibasic furin cleavage site during isolation in VeroE6 cells contributed to the attenuation of the d105 isolate.
In summary, we detected SARS-CoV-2 variants in a persistently infected immunocompromised patient which partially escaped the humoral immune response. Such escape mutants could serve as initial seed for newly emerging variants with enhanced epidemic potential, especially if they overcome impaired viral fitness by further adaptation. Unexpectedly, the spike protein of this escape variant worked particularly well in immunization approaches and elicited broadly active neutralizing antibodies able to control SARS-CoV-2 variants of concern.
Author contributions
MH, GK, CNH, MP, DS, MS, SW and GZ designed the study and contributed to experimental design and data interpretation. DH, JK, HL, SR, YT and DW preformed patient recruitment, clinical management, evaluation of clinical data and sample collection. JF and LK performed bioinformatic analyses. JA, JB, JF, VF, RG, DH, JK, HL, DS, SU, SW and GZ performed experiments and analyzed and processed the data. JF, MH, GK, MP, MS and GZ wrote the manuscript. MH, GK, CNH, MP, MS and RT were involved in funding acquisition.
Competing interests
All authors declare to have no financial or other associations that might pose a potential or actual conflict of interest.
Methods
Case history
A patient in his 50’swas hospitalized at the University Medical Center, Freiburg, Germany for 7 monthHe underwent an organ transplantation and was treated with diverse regimens of immunosuppression, consisting of tacrolimus (4 mg/day), mycophenolat mofetil (2 mg/day) and prednisone (10 mg/day) (see figure 1). The patient received five days ivermectin treatment (33 mg/day) and later a ten days Remdesivir treatment (200 mg on day 1, 100 mg/day 2 to 10).
Ethical statements
The project has been approved by the University Medical Center, Freiburg, ethical committee. Written informed consent was obtained from all participants and the study was conducted according to federal guidelines, local ethics committee regulations (Albert-Ludwigs-Universität, Freiburg, Germany: No. F-2020-09-03-160428 and no. 322/20) and the Declaration of Helsinki (1975). All routine virological laboratory testing of patient specimens was performed in the Diagnostic Department of the Institute of Virology, University Medical Center, Freiburg (Local ethics committee no. 1001913).
Virus detection by qRT-PCR
SARS-CoV-2 RNA testing of respiratory tract samples was performed using the RealStar SARS-CoV-2 RT-PCR kit (Altona Diagnostics, Hamburg, Germany). RNA samples were extracted using QIAamp MinElute Virus Spin kit (Qiagen, Hilden, Germany). Tests were performed and interpreted according to the manufacturer’s instructions and semi-quantitative results reported in cycle threshold (Ct) values.
Serological testing
Convalescent sera of COVID-19 patients and sera from vaccinees after the second dose of the BNT162b2 mRNA vaccine (Pfizer/BioNTech) were obtained from the Hepatology-Gastroenterology-Biobank as part of the Freeze-Biobank Consortium at the University Medical Center Freiburg. Written informed consent was obtained from all blood donors prior to inclusion.
SARS-CoV-2 specific anti-spike protein (S1) IgG (Euroimmun, Medizinische Labordiagnostika AG, Lübeck, Germany) and anti-nucleoprotein (N) IgG ELISA (Mikrogen Diagnostik GmbH, Neuried, Germany) were performed according to manufacturer’s protocol. Results were evaluated semi-quantitatively as arbitrary units (AU) compared to the manufacturer’s calibrators.
To determine the SARS-CoV-2-specific antibodies in mouse sera, VeroE6 cells in 96-well plates were infected with the prototypic Muc-IMB-1virus isolate (kindly provided by Roman Woelfel, Bundeswehr Institute of Microbiology20). Fixed and permeabilized cells were incubated with dilutions of the post-infectious mouse sera and SARS-CoV-2-specific antibodies were detected by fluorescence-labeled secondary anti-mouse IgG antiserum. The serum dilution giving a clear fluorescence signal in the infected cells was interpreted as positive.
For the SARS-CoV-2 neutralization plaque reduction assay, serial serum dilutions were incubated with 100 plaque forming units (pfu) for 1 hour. The mixture was dispersed on VeroE6 cells in 12-well format and the cells were overlaid with 0.6% Oxoid-agar for 48 h at 37°C. The fixed cells were stained with Crystal violet. Number of plaques was compared with an untreated control without serum.
For the detection of neutralizing antibodies by indirect immunofluorescence, 400 pfu of SARS-CoV-2 were preincubated with serially diluted serum samples for 1h and the mixture was used to infect VeroE6 cells in 96-well plates. For each sample, one control without serum was included. Cells were fixed 20 h post infection and stained with anti-SARS-CoV nucleocapsid (N) rabbit antiserum (#200-401-A50, Rockland Immunochemicals). The plates were evaluated by fluorescence microscopy. The highest dilution of the serum that showed less than 50% of infected cells compared to a non-reactive control serum was classified as neutralization titer.
SARS-CoV-2 S1-specific T cell response
The S1 amino-acid sequence of SARS-CoV-2 (GenBank: MN908947.3) was analyzed for in silico prediction of peptide binding with ANN 4.0 on the Immune Epitope Database website21. The HLA-A*02:01-restricted 9-mer peptide, S133-141 FQFCNDPFL, and the HLA-A*03:01-restricted 9-mer peptide, S142-150 GVYYHKNNK, both overlapping with del141-144, the HLA-A*02:01-restricted 10-mer peptide, S244-247 TLLALHRSYL, overlapping with del244-247, as well as a HLA-A*03:01-restricted 9-mer peptide, S378-386KCYGVSPTK, of the S1 domain were selected and synthesized for further analysis. Subsequently, PBMCs (1−2 × 106) of the convalescent immunosuppressed COVID-19 patient and of four HLA-A*02:01/HLA-A*03:01 positive SARS-CoV-2 convalescent immunocompetent donors were stimulated with these peptides (5 µM) and anti-CD28 mAb (0.5 µg ml−1, BD Biosciences) and expanded for 14 days in complete RPMI culture medium containing rIL2 (20 IU ml−1, Miltenyi Biotec). IFNγ production was assessed after a 5h re-stimulation with the respective peptide 31.
Cell culture
Virus isolation, cell culture and mouse infection experiments with SARS-CoV-2 were performed under Biosafety Level 3 (BSL3) protocols at the Institute of Virology, Freiburg, approved by the Regierungspraesidium Tuebingen (No. 25-27/8973.10-18 and UNI.FRK.05.16-29). Filtered throat swab samples were inoculated on African green monkey kidney VeroE6 cells (ATCC CRL-1586) (2 x 106 cells) in 4 ml DMEM with 2% FCS and incubated at 37°C and 5% CO2 for 4 to 6 days until cytopathic effect was visible. The culture supernatant was cleared and stored at -80°C. Virus titers were determined by plaque assay on VeroE6 cells.
For neutralization assays we used established prototypic isolates: Muc-IMB-1, lineage B.120; UK variant B.1.1.7 (hCoV-19/Germany/NW-RKI-I-0026/2020; ID: EPI_ISL_751799) and South African variant B.1.351 (hCoV-19/Germany/NW-RKI-I-0029/2020; ID: EPI_ISL_803957). VeroE6 or human bronchial epithelium Calu-3 cells (ATCC-HTB-55) in 6-well plates, 1 x 106 cells, were infected with a moi of 0.001 for 1.5h. Cells were washed three times with PBS and overlaid with 2 mL DMEM with 2% FCS. The supernatants were taken at different time points after infection and titers were determined by plaque assay. For viral protein analysis, the cells were lysed with tissue protein extraction reagent (T-PER; Thermo Fisher Scientific), separated by 10% SDS-polyacrylamid gel electrophoresis and transferred on PVDF membranes. The membranes were stained with the SARS-CoV-2 N-specific, SARS-CoV-2 spike protein (RBD)-specific (200-401-A50 and 600-401-MS8, Rockland) or actin-specific (A5060, Sigma) rabbit antisera. Detection of the primary antibodies was performed with fluorescent-labeled (Li-COR) secondary antibodies.
Immunofluorescence analysis
VeroE6 cells seeded on glass coverslips were infected with SARS-CoV-2 isolates at a moi of 0.1. At 8 hours post infection, cells were fixed in 4% paraformaldehyde, permeabilized with 0.3% Triton X-100 and blocked in 10% fetal calf serum. SARS-CoV-2 N- and spike-specific primary antibodies and AF568-labeled goat-anti-rabbit (Invitrogen, #A11011, 1:400) secondary antibody as well as AF488-labeled Phalloidin (Hypermol, #8813-01, 1:250) were used for staining. The coverslips were embedded in Diamond Antifade Mountant with 4′,6-diamidino-2-phenylindole (DAPI) (ThermoFisher, #P36971). Fluorescence images were generated using a LSM800 confocal laser-scanning microscope (Zeiss) equipped with a 63X, 1.4 NA oil objective and Airyscan detector and processed with Zen blue software (Zeiss) and ImageJ/Fiji.
Virus pseudotype VSV*ΔG(FLuc) neutralization assay
cDNAs encoding the S protein were prepared from oropharyngeal swab samples of the COVID-19 patient obtained at days d14 and d105 and were cloned into the eukaryotic expression vector pCAGGS32. Single and double spike mutations were introduced into the pCAGGS-S(d14) construct. BHK-21 cells were transfected with the pCAGGS-S plasmids and later inoculated with 5 ffu/cell of VSV*ΔG(FLuc), coding for firefly luciferase, as described33. Cells were incubated in medium containing the monoclonal mAb I1 antibody (ATCC) directed against VSV-G. The supernatants containing the pseudotype viruses were harvested and stored at -70°C.
The pseudotyped virus neutralization test was performed as described recently34. Pseudotyped VSV*ΔG(FLuc) (200 ffu) were preincubated with twofold serial dilutions of convalescent sera in DMEM cell culture medium. The virus/serum mixture was transferred to VeroE6 cells in 96-well plates and incubated for 16 hours at 37°C. Thereafter, the cells were lysed and firefly luciferase activity was determined using ONE-GloTM substrate (Promega) and a GloMax® plate reader (Promega). The reciprocal antibody dilution causing 50% reduction of the luminescence signal was calculated and expressed as neutralization titer 50% (NT50).
Infection of K18-hACE2 transgenic mice
Transgenic (K18-hACE2)2Prlmn mice18 were purchased from The Jackson Laboratory and bred locally. Hemizygous 8-12-week-old animals of both sexes were used in accordance with the guidelines of the Federation for Laboratory Animal Science Associations and the National Animal Welfare Body. All experiments were in compliance with the German animal protection law and approved by the animal welfare committee of the Regierungspraesidium Freiburg (permit G-20/91).
Mice were anesthetized using isoflurane and infected intranasally (i.n.) with virus dilutions in 40 µl PBS containing 0.1 % BSA. Mice were monitored daily and euthanized if severe symptoms were observed or body weight loss exceeded 25 % of the initial weight. Serum samples were collected from the vena facialis. SARS-CoV-2 specific antibody titers were determined by indirect immunofluorescence as described above.
Whole genome sequencing
cDNA was produced from extracted RNA of oropharyngeal swab samples using random hexamer primers and Superscript III (ThermoFisher) followed by a PCR tiling the entire SARS-CoV-2 genome (ARTIC V3 primersets; https://github.com/artic-network/artic-ncov2019). This produced 400 bp long, overlapping amplicons that were subsequently used to prepare the sequencing library. Briefly, the amplicons were cleaned with AMPure magnetic beads (Beckman Coulter). Afterwards the QIAseq FX DNA Library Kit (Qiagen) was used to prepare indexed paired end libraries for Illumina sequencing. Normalized and pooled sequencing libraries were denatured with 0.2 N NaOH. This 9 pM library was sequenced on an Illumina MiSeq instrument using the 300-cycle MiSeq Reagent Kit v2.
For sequencing of virus stocks produced in cell culture, RNA was extracted with the Quick-RNA Viral Kit (Zymo Research) and paired end libraries without a previous PCR amplification were prepared using the TruSeq Stranded Total RNA Kit (Illumina). A total of 10 pM library was sequenced on the Illumina MiSeq instrument.
The de-multiplexed raw reads were subjected to a custom Galaxy pipeline, which is based on bioinformatics pipelines on usegalaxy.eu (doi: 10.5281/zenodo.3685264)35. The raw reads were pre-processed with fastp (https://www.biorxiv.org/content/10.1101/274100v2) and mapped to the SARS-CoV-2 Wuhan-Hu-1 reference genome (Genbank: NC_045512) using BWA-MEM (https://academic.oup.com/bioinformatics/article/25/14/1754/225615). For datasets, which were produced with the ARTIC v3 protocol, primer sequences were trimmed with ivar trim (https://andersen-lab.github.io/ivar/html/manualpage.html). Variants (SNPs and INDELs) were called with the ultrasensitive variant caller LoFreq (https://academic.oup.com/nar/article/40/22/11189/1152727), demanding a minimum base quality of 30 and a coverage of at least 5-fold. Afterwards, the called variants were filtered based on a minimum variant frequency of 10 % and on the support of strand bias. The effects of the mutations were automatically annotated in the vcf files with SnpEff (https://www.tandfonline.com/doi/full/10.4161/fly.19695). Finally, consensus sequences were constructed by bcftools (https://academic.oup.com/bioinformatics/article/25/16/2078/204688). Regions with low coverage or variant frequencies between 30 and 70 % were masked with Ns. The final consensus sequences have been deposited in the GISAID database (www.gisaid.org) (Extended data table 2).
Phylogenetic and variant analysis
All available sequences from Germany deposited in GISAID (http://gisaid.org/) between February and April 2020 were downloaded (as of 11th of February 2021) and 250 sequences randomly subsampled excluding sequences already deposited by the Virology in Freiburg (Extended data table 3). For the phylogenetic analysis, the sequences were first aligned with MAFFT (v7.45)36 and a tree was constructed with IQ-Tree (v2.1.2)37. The best fitting substitution model was automatically determined and the tree was calculated with 1000 bootstrap replicates. Branch support was approximated using the Shimodaira–Hasegawa [SH]-aLRT method (1000 replicates). The tree was rooted to the reference sequence NC_045512. The clades were classified with the webservers of Nextclade (clades.nextstrain.org) and Pangolin (pangolin.cog-uk.io). To visualize the phylogenetic tree a custom R script was written utilizing the ggtree (v2.2.4) (https://academic.oup.com/mbe/article-abstract/35/12/3041/5142656), treeio (v1.12.0) (https://academic.oup.com/mbe/article-abstract/37/2/599/5601621) and ggplot2 (v3.3.3) packages (https://link.springer.com/chapter/10.1007/978-3-319-24277-4_12). An in-house R script was also used to plot the variant frequencies that were detected by LoFreq as a heatmap (pheatmap package v1.0.12). Both scripts are publicly available (github.com/jonas-fuchs/SARS-CoV-2-analyses) and the variant frequency plot has been implemented as a galaxy tool (Variant Frequency Plot on usegalaxy.eu).
Visualization of the spike protein structure
The EM structure of the closed conformation of D614G SARS-CoV-2 spike protein was loaded from the protein data bank (10.2210/pdb7BNM/pdb) and visualized with UCSF ChimeraX version: 1.1 (2020-09-09).
Data availability
All necessary data and informations are given in the manuscript. The sequence data are submitted to GISAID data base and are publicly available (Extended data table 2). Further additional informations about the patient will not be shared due to protection of individuals’ privacy.
Code availability
The scripts are publicly available (github.com/jonas-fuchs/SARS-CoV-2-analyses) and the variant frequency plot has been implemented as a galaxy tool (Variant Frequency Plot on usegalaxy.eu).
HLA-A*02:01/HLA-A*03:01-restricted CD8+ T cell epitopes overlapping with the del141-144 and the del244-247 deletions as well as a control peptide 378-386 in the S1 domain were predicted in silico. Peptide-specific CD8+ T cell responses were assessed by intra-cellular IFN□ production following a 14-day in vitro peptide-specific CD8+ T cell expansion. CD8+ T cells targeting the in silico-predicted SARS-CoV-2 S1-specific, HLA-A*02:01/HLA-A*03:01-restricted CD8+ T cell epitopes were neither detectable in the convalescent, immunosuppressed COVID-19 patient (upper panel) nor in four additional HLA-A*02:01/HLA-A*03:01 positive convalescent, immunocompetent COVID-19 patients (donor 1 to 4). CD8+ T cell responses targeting the immunodominant A*03:01/S378-386 epitope are depicted as positive control, as well as mock-stimulated and PMA/Ionomycin-stimulated cells.
Acknowledgements
This work was supported by the Bundesministerium fuer Bildung und Forschung (BMBF) through the Deutsches Zentrum fuer Luft-und Raumfahrt, Germany, (DLR, grant number 01KI2077) to MP, RT and MS and by the Federal State of Baden-Wuerttemberg, Germany, MWK-Sonderfoerdermaßnahme COVID-19/AZ.:33-7533.-6-21/7/2 to MS and AZ33-7533-6-10/89/8 to CNH. The funders had no role in the study design, data analysis, data interpretation, and in the writing of this report. All authors had full access to the data in the study and accept responsibility to submit for publication.
Footnotes
↵8 These authors jointly supervised this work: Marcus Panning, Georg Kochs.
References
Subject Area
Reviews and Context
0
Comment
0
TRIP Peer Reviews
0
Community Reviews
0
Automated Services
23
Blogs/Media
Author Videos