ABSTRACT
We studied a unique case of prolonged viral shedding in an immunocompromised patient that generated a series of SARS-CoV-2 immune escape mutations over a period of seven months. During the persisting SARS-CoV-2 infection seventeen non-synonymous mutations were observed, thirteen (13/17; 76.5%) of which occurred in the genomic region coding for spike. Fifteen (15/17; 88.2%) of these mutations have already been described in the context of variants of concern and include the prominent immune escape mutations S:E484K, S:D950N, S:P681H, S:N501Y, S:del(9), N:S235F and S:H655Y. Fifty percent of all mutations acquired by the investigated strain (11/22) are found in similar form in the Omicron variant of concern. The study shows the chronology of the evolution of intra-host mutations, which can be seen as the straight mutational response of the virus to specific antibodies and should therefore be given special attention in the rating of immune escape mutations of SARS-CoV-2.
INTRODUCTION
In December 2019 the Wuhan Municipial Health Commission (China) reported a cluster of cases of pneumonia of unknown etiology to the WHO China Country Office. By the beginning of 2020 it was confirmed that a novel coronavirus later named severe acute respiratory syndrome coronavirus (SARS-CoV-2), was the causative agent (1). SARS-CoV-2 spreads easily and effectively among human beings with a basic reproduction number (R0) of > 2 (2, 3). Following this rapid human-to-human transmission and intercontinental spread the WHO declared a global pandemic in March of 2020. The first cases in Austria were reported in Ischgl, Tyrol, as early as February 2020 – and East Tyrol was considered one of the first hotspot areas in Central Europe.
While mutations are common in RNA viruses and mostly will not make a significant difference, some mutations proved to provide SARS-CoV-2 with a selective advantage, such as increased transmissibility or increased escape from specific antibodies (4-8). Those variants with proven or suspected immune escape mutations were deemed variants of concern (VOC) or variants of interest (VOI), respectively, and require close monitoring (https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/). The spread of the first described variant of concern (Alpha variant, B.1.1.7, VOC) was confirmed early in Austria and increased from 0.7% in January to > 99% in April 2021 in the study area East Tyrol (https://www.ages.at/themen/krankheitserreger/coronavirus/sars-cov-2-varianten-in-oesterreich/). As of December 2020, the European Center for Disease prevention and control (ECDC) noticed a sharp decline of the Alpha variant in the European Union (EU; 14.5%) followed by a fast spread of the Delta (B.1.617.2) in spring, 2021, plus a few shares of the Gamma (P1 or B.1.1.28.1, 0.3%) and others virus variants (0.3%). By January 2022, these variants have been largely replaced by the Omicron variant (B.1.1.529) in the study area. The Omicron variant now accounts for > 95% of cases. Common in these successful variants are a few specific genomic changes that give them a decisive benefit in terms of transmission or replication, e.g. in the pharynx. The most prominent spike substitutions with immune escape effect can be found in several VOC or VOI. For example, N501Y in the Alpha, Beta and Omicron variant, a nucleotide exchange at position E484 to K in the Alpha- (in a subvariant), Beta and Gamma, to Q in a Delta subvariant or to A in both Omicron subvariants BA.1 and BA.2. Furthermore, P681H is found in the Alpha and Omicron variant as well as in the VOI B.1.1.238 and the spike deletions 144/145 described as recurrent deletion regions, since they multiply occurred (9, 10).
The origin of immune escape variants is still a matter of speculation. Several hypotheses take zoonotic origin, selective pressure during treatment with antiviral drugs, monoclonal antibodies or convalescent plasma into consideration and a few studies point to the significance of the exceptional intra-host environment of immunocompromised patients to explain the evolution of immune escape variants (11-13). Those cases showed especially long-lasting viable viral shedding of SARS-CoV-2 in immunosuppressed patients for a period of more than four months (11-13). Two of the patients treated with monoclonal and convalescent plasma showed unusually high numbers of nucleotide changes and deletion mutations (12, 13), among others the already described immune escape mutation S:del69/70 (13).
In November 2020 we became aware of a patient in her 60ies with lymphoma who showed persistingly high pharyngeal viral loads of SARS-CoV-2. Although the patient had detectable SARS-CoV-2 specific antibody responses, she was unable to clear the virus load. Therefore, close monitoring was performed by rtPCR of naso-pharyngeal swabs. This unusual case prompted us to perform a thorough sequential serological screening as well as investigation of virus mutant development. Viral titre and mutant development of SARS-CoV-2 was monitored by regular naso-pharyngeal swabs with subsequent rtPCR and Next-Generation Sequencing (NGS), which enabled us to monitor the chronology of the evolution of immune escape mutations. We discuss these findings in the light of potential new sources of intra-host escape mutations and with respect to adaptations of the vaccination strategy.
METHODSRESULTS
Clinical presentation of an immunocompromised individual persistently infected with SARS-CoV-2
In August 2015, the patient was diagnosed with stage IVa small cell lymphocytic lymphoma, complicated by a temporary reactivation of Epstein-Barr virus (EBV) with reactive splenomegaly and rapid nodal progression. From June 2016, she was given six cycles of Rituximab and Bendamustine, which led to remission. In October 2019, the patient suffered a relapse with washout and 90% bone marrow infiltration (B-CLL Binet B or RAI III), accompanied by pronounced B symptoms and antibody deficiency. Beginning in May 2020, another round of therapy with Rituximab and Bendamustine was administered. It was completed in November 2020 after six cycles. The leukocytes counted 4200/µL in the lower normal range, platelets 136,000/µL, the immunoglobulins were clearly reduced (IgG 249 mg/dL, IgA 3 mg/dL, IgM 12 mg/dL).
Four days after the last chemotherapy – mid November 2020 - the patient fell ill with fever, cough, headache and pain, but neither loss of taste nor smell. SARS-CoV-2 was detected in the throat swab by rtPCR. The patient was in quarantine for ten days; a final PCR control was not carried out. Due to persistent fatigue, recurrent fever episodes and persistent cough with little white expectorate, the patient was again admitted to the hospital in the middle of January, 2021 and rtPCR was again positive for SARS-CoV-2. At the same time there was a recurrence of EBV. The patient was enrolled for an inhalation therapy with N-chlorotaurine (3 times daily inhalation of 10ml of N-chlorotaurine for three minutes and 10 - 14 days) (Nagl et al., 2018; Cegolon et al., 2020) and received 15g intravenous immunoglobuline (IVIG, Intratect®, Biotest Pharma GmbH, Dreieich, Germany) on admission. Ten days later; no specific therapy with antibodies against SARS-CoV-2 having been carried out - the blood findings were unchanged (leukocytes 5500/µL, platelets 102,000/µL (decreased), IgG 301 mg/dL, IgA 3 mg/dL, IgM 25 mg/dL), but the leukocyte typing showed a decline by the B-CLL of 70% and a reduction of CD4+ helper T cells. The chest x-ray was normal. The patient suffered from impaired general condition, headache and sore throat and was not able to clear SARS-CoV-2 infection. The patient again received the IVIg therapy (IVIG, Intratect®, Biotest Pharma GmbH, Dreieich, Germany). No antiviral therapeutics were administered to the patient at any time of the infection. An increase of leukocytes was developed (21,600/µL).
Since the onset of the symptoms and the first rtPCR positive swab the patient was committed to home quarantine in accordance with Austrian law. When symptoms did not clear after a month, home quarantine was slightly lightened, yet, testing with naso-pharyngeal swabs and subsequent rtPCR was continued. Persistent viral shedding was determined on day 102, day 124, day 182 and day 205 since onset of symptoms and first positive rtPCR and the viral shedding from the upper respiratory tract has continued to date (Figure 1). The patient gave full written consent for the case to be attended and published.
Timeline of the course of disease in an immunocompromised female patient with almost continuous viral shedding throughout the study period of 207 days. Ct …Cycle threshold; IVIG therapy … intravenous immunoglobuline therapy; RTX … Rituximab and Bendamustine therapy.
Isolation
Isolation trials were performed from swab samples taken on day 73, 93, 99, 104, 109, 117, 127, 133 and 182 on VeroB4-cells and were successful on day 73, 93, 109 and 127. The isolation success correlated negatively with the Ct-values of the swab (k = -0.59).
Serology
At the same time points, specific antibodies were investigated using three different serological tools. An overview of the serological results is given in Table 1. On day 102 the IgG titres in general and neutralizing antibody titres in particular, were observed in a low degree with an IgG titre of 17.7 AU/mL in CLIA and a borderline titre of neutralizing antibodies of 1:4 in enzyme linked neutralisation assay (ELNA)
Humoral immune responses increased in the course of disease and yielded high values of 1320, > 2080 and 1750 AU/mL in chemiluminescent immunoassay (CLIA) and >1:32 in ELNA, on the days 124, 182 and 205, respectively. Although antibodies increased to positive activities at day 124 and were maintained, IgG avidity did not mature over time, as the RAI showed no significant increase and stayed in the low range <20%.
The analysis via immunoblot disclosed Spike 1 (S1) and the dedicated receptor-binding-domain (RBD) as the main epitopes of the IgG antibodies in the patient’s sera, whereas no IgG antibodies could be detected against the region Spike 2 (S2) or the nucleocapsid. No specific IgA antibodies were detected.
SARS-CoV-2 specific T-cell response
On day 193 no IFNg-producing SARS-CoV-2-specific immune cells could be detected in the ELISpot assay (SI = 0.86), although a significant positive reaction against pokeweed mitogen was demonstrated (mean of 213 SFU in the positive control versus mean of 1.4 SFU in the negative control and mean 1.2 SFU cells in the SARS-CoV-2-antigen stimulated wells).
Humoral immune response did not clear SARS-CoV-2 infection
The Ct values and numbers of PFU/mL were significantly lower after day 124. The high titre of IgG antibodies of 1320 AU/mL and a neutralizing antibody titre of 1:32 analysed by our in-house assay on day 124 was associated with a significant reduction of the viral load but could not clear the infection.
We therefore decided to undertake a detailed examination of the specific genetic background of the virus population present including potential intra-host mutational dynamics.
Mutational intra-host dynamics
Over the study period of 221 days, 14 haplotypes were sequenced out of naso-pharyngeal samples. The sequences were obtained on day 73, 93, 109, 129, 133, 136, 143, 158, 164, 171, 182, 192 and day 207 of the patient’s prolonged infection. The timeline of infection and a chronology of intra-host non-synonymous mutational events are given in Figure 4.
The calculation of the pairwise mutation distances did not show higher intra-host evolutionary rates in contrast to overall evolutionary rates of about (8-9) × 10−4 substitutions per year (Day et al., 2020; Dearlove et al., 2020). The pairwise distance between day 73 and day 171 was 4.4 × 10−4 in 98 days, implying a nucleotide substitution rate of 7.5 × 10−4.
All NGS sequences were shown to belong to the prevalent Pangolin lineage B.1.1. and the Nextstrain clade 20B.
We became aware of the prolonged viral shedding after about two months and started to regularly sequence the patient’s subsequent swabs as of day 73. The sequence derived from the swab of day 73 showed 18 mutations in comparison to the reference genome Wuhan (GenBank: MN908947.3, RefSeq: NC_045512.2), the sequence on day 171 had 26 substitutions which increased to 29 on day 207. A listing of all persistent and temporary non-synonymous mutations that the strain has accumulated intra-host and their concordance to variants of concern and variants of interest are given in Table 2. Overall, 22 non-synonymous mutations evolved over the study period of 221 days (7 months). Eleven (50%) of these non-synonymous mutations were persistent, whereas 11 (50%) occurred temporarily and were replaced by the wildtype or a different substitution. Seventeen of the 22 non-synonymous mutations evolved in the region coding for spike, eight of those were temporary. Seventeen of the 22 acquired non-synonymous mutations (77.3%) were issued as immune escape mutations by the WHO (https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/) (Figure 4). Among the persistent non-synonymous mutations in spike, as many as 88.2% are found in various VOIs or VOCs (Table 2). All those genetic changes occurred after the development of high antibody titres.
One region continuously showed diffuse mutational changes, with changing temporary adaptations of substitutions and deletions, which was ORF1b: position 709 – 716.
An overview of the intra-host mutational development of SARS-CoV-2 during the study period of 221 days is given in Table 3 and Figure 3.
Chronology of acquired mutations
In the underlying clinical case the substitutions emerged in the following chronological order: S:Y144-emerged immediately after the increase of the specific antibody titre at day 117 as a temporary mutation, followed by E484Q (day 129), which could not assert itself against E484K and was displaced at least 7 days later (day 136). Furthermore, we found the substitutions S:N354K (day 158, 164, 171 und 182), S:R346I (day 164) and ORF1a:T3284I (day 171), S:D950N (day 171) as well as the prominent S:P681H on day 182 (Fig. 2). Three of the six acquired substitutions (50%) have already been described as typical mutations acquired by diverse variants of concern.
Localization of the substitutions and deletions acquired by the strain EPI_ISL_2106199 of clade 20B in the course of the prolonged infection of an immunocompromised patient. All acquired and preserved genetic adaptations occurred in the regions ORF1a (n = 1), ORF1b (n = 1) and spike (n = 9).
Thirteen of the seventeen acquired substitutions (76.5%) occurred in the genomic region coding for spike, and one each in the regions coding for ORF1a:T3284, ORF3a:V255X (day 73), ORF8:Y73C* (day 73) and N:S235F* (day 136) (Fig. 2).
Other mutations appeared temporarily and were subsequently replaced by the wildtype variant. Five hitherto undescribed temporary mutations were observed on the days 73 (ORF3a:V255X), 117 (S:A831V), 117 - 123 (S:Y145X), day 129 (S:Y144H), as well as on day 129 (S:E484Q, S:Y144H) (Fig. 3).
The acquired and temporarily acquired mutations of strain EPI_ISL_2106199 in the course of a seven months lasting course of infection in an immunocompromised person in the region coding for spike. Overall, 17 persistent or temporary spike mutations were evolved, whereas 9 (52.9%) turned out to be temporary and were subsequently replaced by the wild-type variant. * … temporary mutations; S1 … spike 1; S2 … spike 2; hr … heptad repeat; RBD … receptor binding domain; The mutations marked in orange are also found in the Omicron variant (B.1.1.529) in similar or identical expression (10 out of 17), mutations marked in red are found in other variants of concern (3 of 17; 17.6%). Thirteen of the 17 mutations (76.5%) acquired in the course of the prolonged infectious phase are already described mutations in variants of concern.
Six temporary mutations have already been described previously, all of them prominent variations known in the context of variants of concern (https://covariants.org/shared-mutations; Center for Disease Control and Prevention; https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html): ORF8:Y73C and S:T716I evolved in the early stage of the infection and are described in the context of the Alpha variant,. S:N501Y, known from the Alpha, Beta, Gamma and Omicron variant, S:del(9) described in the Beta variant, N:S235F, known as a typical substitution of the Alpha variant and S:H655Y, described in the context of the Beta as well as the Omicron variants. All these temporarily recurring mutational events did not establish permanently but disappeared again or were dominated again by the wildtype variant. Figure 3 shows the acquired und temporarily acquired mutations of strain EPI_ISL_2106199 in the region coding for spike and demonstrates the high concordance of the acquired adaptations with described variants of concern, above all the Alpha and the Omicron variant (15/17; 88.2%). Figure 3 represents the adaptations in the region coding for spike. Further mutations in the genome were ORF1a:T3284I (day 171), ORF8:Y73C (day 73), ORF3a:V255X (day 73), N:S235F (day 136), ORF1b:L714- (day 158). Thirteen of the 17 mutations (76.5%) acquired in the course of the prolonged infectious phase are already described mutations in variants of concern. Ten of the 17 spike mutations occur in a similar or identical way in the Omicron variant (58.8%). The non-synonymous mutations S:del143, S:del144, S:N501Y, S:H655Y and S:P681H were developed in identical form in the Omicron variant. Further non-synonymous mutations occurred at the amino acid positions S:142, S:144, S:145, S:484 (twice) in strain EPI_ISL_2106199 as well as in the Omicron variant, which, however, led to different expressions (S:del142 instead of S:G142D, temporarily both S:Y144H versus S:144del and S:Y145X versus S:145del as well as S:E484Q and S:E484K instead of S:E484A). Overall, 76.5% of all mutations and 88.2% of the spike mutations acquired by strain EPI_ISL_21061 convergently evolved in other variants of concern, mainly in the Alpha and the Omicron variant.
Overall, SARS-CoV-2 developed eleven persistent mutations during the study period of 140 days as well as eleven temporary mutational events. The chronology of intra-host mutational events is displayed in Figure 4.
Chronology of the emergence of intra-host mutations in an immunocompromised patient with adequate humoral and lacking cellular immune response. The study period comprised 140 days of almost permanent viral shedding. High-quality next-generation sequences could be obtained at 14 time-points during the seven-month study period (starting on day 73, ending on day 207 with the last SARS-CoV-2 positive swab) and disclosed the chronological development of mutational events of SARS-CoV-2 as an answer to a unilateral immune response with strong antibody answer but lack of specific T-cells.
In the first swab sample, whole genome sequencing did not detect any spike mutations in the investigated strain compared to the reference genome. First spike variants appeared as E484K on day 133 as a heterozygotic mutation in 41.3% of the targeted reads. On day 136 the proportion of E484K increased to 76% and, after more than seven days (day 143) the new variant dominated with 100%, but decreased to 76.8% again on day 158. On day 171 the spike variant P681H was observed for the first time with a proportion of 24% and dominated within a couple of weeks reaching 100% on day 182.
Three of the six acquired substitutions (50%) are previously described substitutions of immune escape variants, namely: S:E484K, S:D950N and S:P681H.
The fluctuating occurrence of adaptive mutations
The emergence of adaptive mutations did not occur in a linear but fluctuating fashion. Frequently, new mutations arose at a certain time point to be later replaced by the wildtype variant, again. As shown in Figure 5, the mutation rate shows an oscillating course with peaks around day 125, increasing until day 182. Simultaneously, the viral load decreased continuously until the patient had several consecutive negative SARS-CoV-2 qPCR tests since day 232 and is therefore considered to be cured from COVID-19.
Chronology of the appearance of convergent intra-host mutations in a strain of the clade B.1.1., as they have been proven in identical form in the variants of concern. α … mutations are described for the Alpha variant B.1.1.7; β … Beta variant B.1.351; δ … Gamma variant B.1.1.28.1; γ … Delta variant B.1.617.2; o … Omicron variant B.1.1.529. The red line graph shows the mean frequency of all mutations at a given day. The occurrence of mutations did not occur straight-line, but in a fluctuating course, with frequent replacement by the wildtype variant.
Intra-host evolutionary history
The intra-host evolution of the strain EPI_ISL_2106199 from day 73 and the quasi-species arising from it in the course of the intra-host evolution form a distinct clade in the consensus tree and group together. The clade is embedded in a random composition of complete Austrian strains and variants of concern found in Europe and described and uploaded to GISAID platform in the same study period from January to May 2021 (Figure 6). Early Austrian sequences of the Omicron variant from December 2021 were included subsequently.
Outgroup-routed consensus tree based on 40 SARS-CoV-2 whole genome sequences with 29,806 nucleotide sites. The section highlighted in red shows the monophyletic clade of variants newly formed in an immunocompromised patient, embedded in prominent variants of concern and typical Austrian strains sequenced in the same investigation period, downloaded from GISAID and completed by early Austrian sequences of the Omicron variant in December 2021. Numbers at nodes indicate bootstrap support values (only values > 50 are shown).
DISCUSSION
In this unique case report we described the dynamics of intra-host mutational events in an immunocompromised patient during a seven-months period of prolonged viral shedding and proven infectivity. We considered the possible influence of a quantitatively strong but regarding binding capacities probably functionally ineffective humoral antibody response and a lacking cellular immune response on the site-directed mutagenesis of SARS-CoV-2.
Our sequencing approach resulted in high-confidence variant identification and robust genome-wide coverage and enabled the establishment of a chronology of immune escape mutations. In addition, previously undescribed site-directed base-exchanges, found in the regions ORF1a and b (n = 2;ORF1a:T3284I, ORF1b:L714-), ORF3a (ORF3a:V255X) and spike protein (n = 2; S:N354K, S:A831V), were described here. Four different well-established serological methods gave insights into the humoral immune response and demonstrated the inability to clear the SARS-CoV-2 infection despite positive antibody responses. This may be due to the relatively low neutralizing ability of the detected IgG which is also supported by the low avidity of the specific IgG and impaired avidity maturation over time. Administration of IVIG was not able to enhance the clearance of SARS-CoV-2. Cellular immunity was diminished in this patient and the lack of adapted T cell-mediated immune defence may have contributed to the inefficient clearance.
The substitution rate for SARS-CoV-2 was estimated as (8-9) × 10−4 nucleotides per site per year (32). This is comparable to previously reported substitution rates of SARS-CoV (8.0-23.8 * 10−4) (33) and MERS-CoV (11.2 * 10−4) (34, 35) and comparable to the reported substitution rates for Influenza A (4-5 * 10−3) and Influenza B (2 × 10−3) virus in the haemagglutinine gene (36). From this substitution rate it can be estimated that SARS-CoV-2 undergoes about one genetic change every other week (32). In comparison, the nucleotide substitution rate per site and per year for Ebola (EBOV Makona) is estimated to be ∼1.2 * 10−3 (37) and for HIV-1 (3.21-4.06) * 10−3 (38). Interestingly, the evolutionary rate stayed constant throughout the first months of infection but decreased slightly after the increase of specific antibodies on day 124. We did not find elevated intra-host substitution rates compared to the general rate reported for SARS-CoV-2 (32). This unaltered intra-host evolutionary rate compared to the global average evolutionary rate suggests that these mutations in an immunocompromised patient, driven by specific antibodies, do not lead to more frequent random genomic changes, but on the contrary to very specific targeted ones.
We assume that the presence of specific antibodies forced directional selection on retaining or regaining infectiousness and thereby strongly favoured directional mutations at particular sites, acting as immune escape mutations.
E484K, a substitution in the receptor binding domain (RBD) appeared early in the course of the infection and is described to impair neutralization resistance (39), potentially compromising vaccines effectiveness (4, 6, 40-47). E484K is a well-established distinction of the variants of concern (VOC) B.1.1.7 with E484K, P.1, P.2, B.1.315, B.1.525, B.1.526 as well as B.1.617.1 (Center for Disease Control and Prevention; https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html). At the same position E484, both subtypes of the Omicron variant have formed the alternative substitution alanine A (Center for Disease Control and Prevention; https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html). On day 158, the Omicron-specific mutation S:H655Y (48, 49) could also be detected as a temporary substitution. A further genomic change in the region coding for spike was identified on position S:N354K on day 158 and had never been described before. R346I was detected in the sequence of day 164. This mutation was previously described as a reaction of SARS-CoV-2 after monoclonal antibody treatment, seeming to maintain ACE2 binding activity (50) and has also developed in the VOI Mu, 21H, B.1.621 (https://covariants.org/shared-mutations; Center for Disease Control and Prevention; https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html). S:D950N arose around day 171 and is as adaptive mutation assigned to the Gamma and Delta variant (Center for Disease Control and Prevention; https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html).
The substitution P681H was observed for the first time in the sequence of day 171, whereby the amino acid histidine (H) first appeared as a polymorphism to become the dominant and finally fixed variant in the course of the next ten days. This transformation from proline (P) to histidine is relatively well studied and implicates a modification in the neighbouring furin cleavage site at the junction of the spike protein receptor-binding (S1) and fusion (S2) domains (51).
Another major transformation of the spike protein is the deletion of the amino acids 141 to 144. The deletions Y144/145-on the edge of the spike tip are modifications described in the VOC B.1.1.7 as the recurrent deletion region 2 (rdr2), occurring repeatedly in SARS-CoV-2 variants (9, 10).
In our case, the deletions were extended to three more deleted positions on S:141, 142 and 143. S:143del is another analogy to the Omicron variant B.1.1.592.
Nine permanent mutations were found in the spike-coding region. More precisely, four are located in the N-terminal domain (S:L141-, S:G142-, S:V143-, S:Y144-), three in the receptor-binding domain (S:R346I, S:N354K and S:E484K), both parts of S1 and three are positioned in the region encoding for S2, namely P681H in the immediate neighbourhood of the furin cleavage site and S:D905N near heptad repeat 1 as part of the fusion core region.
Of the eleven acquired adaptive mutations, only two were found outside the spike-coding regions, namely L714-in ORF1b (day 158) and T3284I in ORF1a (day 171). ORF1a and ORF1b are coding regions for non-structural proteins (nsp) (52). ORF1a:T3284I is located in the region encoding for nsp5. Nsp5 is regarded as the main protease, cleaves viral polyprotein and works closely with nsp12 and nsp13. Together, nsp5, 12 and 13 represent the replicase machinery (52-54). ORF1b:L714-is a deletion in the region coding for nsp13, the enzyme helicase, a main component of membrane-associated replication-transcription complexes (52, 55-57).
It is remarkable that nine of the eleven persistent mutations (81.8%) acquired in the course of the prolonged infection had previously been described in the context of immune escape and were assigned to diverse variants of concern. Our bioinformatic analyses revealed that 75% of the novel mutations in our investigated strain also occur in variants of concern, whereas the highest concordance was found between strain EPI_ISL_2106199 and the Omicron variant (50%). Furthermore, we found dynamic mutational events with fluctuations between the wildtype and the variational mutation. Nine of these temporary mutations (9 of 11; 81.8%) have also been described in the context of variants of concern (Center for Disease Control and Prevention; https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html). In the close proximity of the acquired deletion in ORF1b:L714-, which manifested homozygously, additional conspicuous polymorphic sequences in amino acid position ORF1b:708-716 were found. We measured frequent changes in substitutions and deletions. This certainly left us with the impression of mutational escape manoeuvres. Also, this hotspot in directional mutations encodes for nsp13, the helicase.
We thus suggest that the accumulated mutations are results of an increased selection pressure on spike, the key to entering the host cell. At the same time a second process takes place intra-host, which exerts increased pressure and enforces continual reconstructions in nsp13. The findings of these temporary mutations, which almost exclusively occurred in the spike region, also fit this pattern very well.
We managed to isolate SARS-CoV-2 from swabs at different time points, which is further evidence for the continuous viability of the virus over the study period, given the evolutionary dynamics of the different sequences. The isolation success correlated negatively with the Ct value, a fact that has already been observed in previous studies (22).
Treatment with Rituximab resulting in depletion of particularly memory and effector B cells by targeting CD20 is known to cause impaired antibody responses (58-60). As naïve B cell clones are less sensitive to Rituximab treatment due to their lower expression of CD20, a robust immune response can also be assumed for those patients. Immunosuppressive therapy as well as the lymphoma disease itself may have diminished the T cellular axis of immune defence against SARS-CoV-2, which targets infected cells, and loss of control by cytotoxic T cells may have caused the ongoing replication of SARS-CoV-2 in naso-pharyngeal epithelial cells (61, 62). Impaired T cell help may have contributed to the inefficient antibody maturation.
Meanwhile, there are more studies that shed light on the evolution of immune escape variants in immunocompromised patients and support the results of our study (12, 63-69).(12, 63-69) Nonetheless, our study not only shows the accumulation of an unusually high number of immune escape mutations in a single patient, which to a strikingly high degree evolved in parallel in various variants of concern. The chronology of mutation events during seven months of infection shows a rapid accumulation of non-synonymous mutations which in part were persistent, in part temporary or even repeatedly acquired and lost.
In summary, our case report documents the medical phenomenon of persisting SARS-CoV-2 infection in an immunocompromised patient with impaired humoral and cellular immune response. Potential interference of specific antibodies led to a significant reduction in the viral load, but at the same time generated sophisticated escape mechanisms while the cell-mediated immune defence for eradication of the infection is missing. With the aid of NGS, we witnessed the directed mutational changes of SARS-CoV-2, probably facilitated by insufficient humoral immune defence. This led to the formation of highly specific virus variants, highlighting the regions exposed to the highest intra-host selective pressure. Based on this observation one may hypothesize that immunocompromised patients pose a particular risk to establish a source of immune escape mutants of SARS-CoV-2. Our study also underlines the importance to protect these patients from SARS-CoV-2 infection by modified vaccination strategies as well as to reinforce vaccination efforts to increase herd immunity in general. Most importantly, the study points out the convergent evolution of specific mutations in SARS-CoV-2, both in VOCs, VOIs and intra-host in the strain we studied (EPI_ISL_2106199). Those specific, convergently evolving mutations reveal those neuralgic positions in the SARS-CoV-2 genome that on the one hand represent its highest fitness advantage, but on the other hand also uncovers its highest vulnerability and should be considered as the probably most important points of attack in future vaccine and therapeutics development.
METHODS
Immunological diagnostics
CLIA SARS-CoV-2 TrimericS IgG
Serological tests were performed using the LIAISON® SARS-CoV-2 TrimericS IgG (DiaSorin S.p.A., Saluggia, Italy) (LIAISON), an Immunoblot called ViraChip® assay (Viramed, Munich, Germany) and an in-house enzyme-linked neutralization assay (ELNA) (14) at day 102, 124, 182 and 205 after the first positive PCR.
The LIAISON® SARS-CoV-2 TrimericS IgG is a CLIA (Chemiluminescent Immunoassay) which detects IgG antibodies reactive with the spike protein (S1/S2 domain). The assay was performed on the LIAISON® XL Analyzer according to the manufacturer’s instructions and gives the arbitrary units per ml (AU/mL) according to the WHO International Standards for the Anti-SARS-CoV-2-immunoglobulin-binding activity (NIBSC 20-136).
Microarray immunoblots
The ViraChip® assay detects temporal antibody profiles of different immunoglobulin classes against S1, S2, and nucleocapsid (N) as well as against N of the four nonSARS human coronaviruses in a commercial, miniaturized 96 wells protein microarray. The ViraChip® assay is a useful tool to identify the epitope-specificity of IgG and IgA in serum samples. The quantitative antibody measurement was performed on a ViraChip® Scanner using ViraChip® Software.
Neutralization test
Neutralization ability of antibodies was determined performing an in-house enzyme-linked neutralization assay (ELNA) as described elsewhere (14).
Anti-IgG-SARS ELISA
Serum IgG antibodies against SARS-CoV-2 were determined by Serion agile SARS-CoV-2 ELISA with a sensitivity of 96.2% and a specificity of 100% according to manufacturer’s instructions (Virion/Serion, Wuerzburg, Germany). Antibody activities above 15U/mL were considered positive.
Anti-IgG-SARS-Avidity
Relative avidity index (RAI) was determined by a modification of the Serion agile SARS-CoV-2 IgG-SARS ELISA using 1M ammonium thiocyanate (NH4SCN) as a chaotropic agent as described previously (15-17). RAI values were considered as: RAI > 60% high avidity, 40% < RAI < 60% as moderate, and RAI < 40% as low avidity in reference to other viral infections (18).
SARS-CoV-2 specific T cell response
The ELISpot assay was performed using a commercially available precoated human SARS-CoV-2-specific IFN-γ ELISPOT kit according to the manufacturer’s protocol (AutoImmun Diagnostika, GmbH, Germany; Cat.no. ELSP 5500). Peripheral blood was collected into tubes coated with lithium-heparin (Vacuette, Greiner bio-one, Austria). PBMCs were separated from plasma and whole blood by gradient density (FicoLite® -H, Linaris, Germany). After washing with phosphate-buffered saline (PBS), depleting erythrocytes (RBD-Lyse Buffer Life Technologies, 1xRBC Lysis Buffer 200ml; Invitrogen eBioscience, USA REF: 00-4333-57) and washing again with PBS, cells were counted and resuspended in x-vivo medium (X-VIVO TM-10 Serum-free hematopoietic cell medium; BEBP02-055Q, Lonza, Switzerland).
Briefly, a total of 2 × 105 PBMCs were incubated in duplicate with x-vivo as a negative control, pokeweed mitogen (AutoImmun Diagnostika GmbH, Germany) as a positive control,15-20mer peptide pools for SARS-CoV-2 (AutoImmun Diagnostika GmbH, Germany) and PanCorona (AutoImmun Diagnostika GmbH, Germany) for the four nonSARS human coronaviruses 229E, HKU1, NL63 and OC43 as a control of possible cellular cross-reactive responses. After incubation at 37° C for 20 hours in a sterile and humidified atmosphere, plates were washed with washing buffer (AutoImmun Diagnostika GmbH, Germany) and stained with the kit-specific reagents according to the manufacturer’s protocol. Plates were then washed several times under running water and dried overnight. Spot forming units (SFU)/100.000 cells were counted using an automated AID ELISPOT reader system (AutoImmun Diagnostika GmbH, Germany).
The assessment criteria for the ELISpots were a minimum of 50 SFU in the positive control and a maximum of 10 SFU in the negative control according to the manufacturer’s definitions (19, 20). When those criteria were fulfilled, the stimulation index (SI) was calculated by dividing the mean SFU numbers in the antigen-specific wells with the mean SFU numbers of the negative control. The test was assessed negative with an SI < 2 according to previous determination of the cut-off by well-defined pre-pandemic PBMC samples and by PBMCs from SARS-CoV-2-naive individuals. The test was suggested to be poorly reactive with an SI between 2 and 7 and reactive with an SI ≥7 as defined by the manufacturer (19). According to standardized laboratory procedures, in each assay, a standard laboratory control sample of a high-reactive and a non-reactive PBMC sample, respectively, was run to determine inter-assay-variations. Only assays with less than two standard deviations of the high-reactive and the non-reactive PBMC control sample, respectively, were defined valid.
Sample collection
Nasopharyngeal swabs were taken in a standardized way in home quarantine in the context of primary care by a medical co-worker.
RNA extraction
Nucleic acids were isolated using the MagMAX™-96 Total RNA Isolation Kit (Thermo Fisher Scientific, Waltham, Massachusetts, USA; Cat. No. AM1830). Briefly, 200µl PBS were taken from patient swab sample and mixed with 265µL binding buffer, 5µl proteinase K (20mg/mL) and 5µL extraction control (Thermo Fisher Scientific, Waltham, Massachusetts, USA; Cat. No. AM1830) according to the KingFisher™ extraction protocol for 200µL sample volume (Thermo Fisher Scientific, Waltham, Massachusetts, USA). After incubation at room temperature for at least 15 minutes, samples were transferred from tubes into 96-well KingFisher deep well plates (Thermo Fisher Scientific, Waltham, Massachusetts, USA) containing 280µl isopropanol and 2µL Mag-Bind particles per well, using a KingFisher™ Flex purification system (Cat. No. 5400620).
rtPCR
qPCR extracts were tested for SARS-CoV-2 by qRT-PCR using the Bio-Rad CFX96 system (Bio-Rad, Germany) with a LightMix Modular Assay kit in accordance with the modified Charité guidelines (Corman et al., 2020). 10µL of extracted RNA were added into 15µL 4x Reliance One-Step Multiplex Supermix (Bio-Rad, Germany). Each 15µL mastermix contained 12.5µL buffer solution, 0.25µl enzyme mix, 1.75µL of nuclease-free water and 0.5µL primer probe wHCoV (E-Gene, as well as N-Gene and Rdrp-Gene for confirmation). Reactions were incubated at 55° C for 5 min and 95 C for 5 min in order to conduct reverse transcription of viral RNA, sample denaturation and enzyme activation. These steps were followed by PCR-amplification including 45 cycles at 95° C for 5 s, 60° C for 15 s and 72° C for 15 s. Cooling was implemented at 40° C for 30 s. Results were interpreted based on the Second Derivative Maximum (SDM) method. Positive results were confirmed by Rdrp and N-gene (21), samples with an initial Ct value lower than or equal to 37 were assigned to repeated testing including extraction. A Ct value higher than 40 was considered negative. Quantification of the viral load in the swabs was calculated via size standards of 1, 10, 100 and 1000 plaque-forming units (PFU)/mL. Standardization of viral stocks was carried out by virus titration. Isolation was performed on VeroB4 cells as described elsewhere (22).
Virus titration
Confluent VeroB4 cells were cultured in Medium199 including 5% FCS in T75 tissue culture flasks (Sarstedt, Germany) and transferred into 96-well tissue culture plates (Sarstedt, Germany). Passage 1 isolates of SARS-CoV-2 were thawed from −80° C freezer and titrated from 1:10 to 1:10−12 in U-shaped 96-well plates (Greiner, Germany) and pipetted into each corresponding well of the 96-well tissue culture plate. Plates were incubated at 36° C. Three days post infection, incubation was stopped by gently removing the supernatant, washing the cells three times with PBS and fixing cells in 1:1 ice-cold acetone-methanol. For easier optical evaluation, cells were dyed by crystal violet staining and tissue culture infectious dose of 70% (TCID70) and PFU were calculated (Ramakrishnan, 2016).
Whole genome sequencing and mutational analysis
Libraries were prepared according to the Ion AmpliSeq(tm) SARS-CoV-2 Research Panel (Thermofisher, USA), library construction and sequencing protocol with the Library Kit Plus (Thermo Fisher Scientific, Waltham, Massachusetts, USA; Cat. No. 4488990). The Amplicons were cleaned up with AMPure XP beads (Beckman Coulter, Germany) with a 1:1 ratio. The libraries were quantified using the Ion Library TaqMan™ Quantitation Kit (Cat. No. 4468802), normalizing, pooling and sequencing was performed using an Ion Torrent™ S5 Plus. Ion Torrent Suite software (v 5.12.2) of the Ion S5 sequencer was used to map the generated reads to a SARS-CoV-2 reference genome (Wuhan-Hu-1; GenBank accession numbers NC_045512 and MN908947.3), using TMAP software included in the Torrent Suite. The following plugins were used: Coverage Analysis (v5.10.0.3), Variant Caller (v.5.12.04) for mutation calls both with “Generic □S5/S5XL (510/520/530) □Somatic □Low Stringency” and “Generic□S5/S5XL (510/520/530) □Germ Line□Low Stringency” default parameters and COVID19AnnotateSnpEff (v.1.0.), a plugin specifically developed for SARS=CoV=2 that can predict the effect of a base substitution. No ultra-deep sequencing was performed and only mutations visible in the stated analysis methods were listed and rated.
FASTA files containing the raw reads were inspected for quality criteria (mapped, targeted, filtered reads, mean depth and uniformity) using Thermofisher Software. Multiple sequence alignments were performed using Unipro UGENE (23) as well as MEGA X (Kumar et al., 2018). The SARS-CoV-2 genomes were compared to the reference NC 045512.2-Wuhan-Hu-1. Viral genome assembly and screening for distinct mutations was performed online using nextstrain.org (https://github.com/nextstrain/ncov/blob/master/defaults/clades.tsv; https://clades.nextstrain.org/). The identification of pangolin lineages was carried out using Pangolin software, v.2.4.2. (https://pangolin.cog-uk.io/). The generated full-genome sequences are available at GISAID EpiCoV (https://gisaid.org/no.EPI_ISL_2106191-21061201). Additional sequences of frequent Austrian strains and prominent variants of concern, sequenced in the same study period, were retrieved from the GISAID EpiCoV database (24) to calculate a consensus tree.
Indels were coded using ‘simple indel coding’ (25) as implemented in 2matrix v.1.0 (70).
The best-fit model of nucleotide substitutions (TIM2+F+I) was selected under the Akaike and the Bayesian formation criteria using ModelFinder (71) as implemented in the PhyloSuite Software package (26). Phylogenomic inference was based on a Maximum Likelihood (ML). An ML tree with 5000 ultrafast bootstrap replicates was inferred in the IQ Tree plugin (27, 28) of PhyloSuite.
Bioinformatics CoV-Seq Workflow
The frequency of the various mutations and the homology to the most widespread variants of concern (Alpha-, Beta, Gamma- and Delta-variant) were investigated based on BAM files. Reads from CoV-Seq samples were demultiplexed by using in-house tools. Reads originating from human were filtered out by mapping against hg38 with bwa-mem 0.7.17 (29). All reads not mapping to human were trimmed for adapters und quality by using Cutadapt 3.2(30). The trimmed reads were mapped with bwa-mem 0.7.17 to the SARS-CoV-2 reference MN908947.3 from the NCBI. Mutations were called using breseq 0.35.5 (31). Graphics were created using pandas 1.2 for Python 3.
Statistics
Dichotomous data were evaluated by a chi-squared test or Fisher’s exacta in the case of small group size (n < 60) (Microsoft® Excel®, Microsoft 395 MSO, Windows 10). A two-sided significance level of p < 0.05 was used for determining statistical significance. After testing for distribution (Kolmogorov-Smirnov-test), non-parametric continuous independent variables were compared using Mann-Whitney-U test for each time point. Dependent non-parametric variables were compared using Wilcoxon-rank test.
Data Availability
The generated full-genome sequences are available at GISAID EpiCoV (https://gisaid.org/no.EPI_ISL_2106191-21061201).
AUTHOR CONTRIBUTION
Sissy Therese Sonnleitner: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology; Writing – original draft preparation. Martina Prelog: Writing – original draft preparation, editing; Stefanie Sonnleitner: Methodology, Writing – editing. Eva Hinterbichler: Investigation; Methodology; Hannah Halbfurter: Project administration. Dominik B. C. Kopecky: Visualization, Organization. Giovanni Almanzar: Methodology. Stephan Koblmüller: Conceptualization, Writing – editing, Funding acquisition. Christian Sturmbauer: Conceptualization, Writing – editing, Funding acquisition. Leonard Feist: Validation, Verification. Ralf Horres: Validation, Verification. Wilfried Posch: Validation, Supervision. Gernot Walder: Supervision.
ETHICAL APPROVAL
All clinical samples and data were collected for routine patient care and for public health interventions. The patient gave full written consent for the case to be attended and published. Ethical approval to use residual routinely taken serum samples for retrospective analyses was obtained by the Ethics Committee of the University Hospital Wuerzburg (no. 20201105_01).
COMPETING INTERESTS
The authors have no conflict of interest to declare.
DATA AVAILIBILITY
The generated full-genome sequences are available at GISAID EpiCoV (https://gisaid.org/no.EPI_ISL_2106191-21061201).
ACKNOWLEDGEMENT
We gratefully acknowledge the financial support of the Austrian Research Promotion Agency (FFG), Grant No. 889135. The authors report no potential conflict of interest. We also thank Ramona Polster, MSc, for excellent technical assistance, DDI Martin Lamprecht for valuable suggestions and all the kind colleagues from the routine diagnostics team of the Dr. Gernot Walder laboratory under the direction of BSc Viktoria Muehlmann and Mag. Hannes Mahl for great support.