Drivers of adaptive evolution during chronic SARS-CoV-2 infections

Sheri Harari, Maayan Tahor, Natalie Rutsinsky, Suzy Meijer, Danielle Miller, Oryan Henig, Ora Halutz, Katia Levytskyi, Ronen Ben-Ami, Amos Adler, Yael Paran, Adi Stern
doi: https://doi.org/10.1101/2022.02.17.22270829

Abstract

In some immunocompromised patients with chronic SARS-CoV-2 infection, dramatic adaptive evolution occurs, with substitutions reminiscent of those in variants of concern (VOCs). Here, we searched for drivers of VOC-like emergence by consolidating sequencing results from a set of twenty-seven chronic infections. Most substitutions in this set reflected lineage-defining VOC mutations, yet a subset of mutations associated with successful global transmission was absent from chronic infections. The emergence of these mutations might dictate when variants from chronic infections can dramatically spread onwards. Next, we tested the ability to predict antibody-evasion mutations from patient- and viral-specific features, and found that viral rebound is strongly associated with the emergence of antibody-evasion. We found evidence for dynamic polymorphic viral populations in most patients, suggesting that a compromised immune system selects for antibody-evasion in particular niches in a patient’s body. We suggest that a trade-off exists between antibody-evasion and transmissibility that potentially constrains VOC emergence, and that monitoring chronic infections may be a means to predict future VOCs.

Introduction

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infections normally resolve clinically within a few days, and RNA shedding lasts a few days to a few weeks 1. However, more and more case reports are accumulating that document chronic infections spanning weeks to many months of infection 2,3. Notably, chronic infection should not be confused with long-COVID, where infection is cleared rapidly, yet symptoms persist 4; in cases of chronic infections, replicative virus is detected for lengthy periods of time. To date, all cases of chronic SRAS-CoV-2 infection were associated with a severely immunosuppressed condition, ranging from primary immunodeficiencies, patients following solid organ transplants, acquired immunodeficiency syndrome (AIDS), patients undergoing therapy with immunosuppressive drugs, and hematologic cancers patients. Presumably, these immune system disorders prevent clearance of the virus as compared to patients with an intact immune system, and thus the virus thrives for lengthy periods of time 5.

Longitudinal sequencing of some cases of chronic infection has uncovered striking mutational patterns of evolution, and revealed that the rate of evolution is much higher than that observed along transmission chains of acutely infected individuals (e.g., 6-8). Indeed, several recent reports have found that the intra-host variation during acute infections is quite limited 9-12. Moreover, the transmission bottleneck size was inferred to be quite low in these studies, as evidenced by little to no shared genetic diversity in transmission pairs. Overall, this suggests that adaptive evolution may be limited during acute infections, and that most of the diversity observed during SARS-CoV-2 circulation is due to neutral mutations fixed during the small transmission bottleneck. The low genetic diversity observed in acute infections contrasts the high genetic diversity observed in rare cases of chronic infections.

Much interest has arisen in cases of chronic infection since the first emergence of variants of concern (VOCs). One major hypothesis that came up when the Alpha VOC was first discovered, and more recently with the rise of the Omicron variant, was that these variants may have evolved due to selection pressure created during chronic infections 2,13. More specifically, the hypothesis posited that strong selection for antibody-evasion mutations may occur in immunocompromised individuals with chronic infection who are treated with convalescent plasma (CP) or monoclonal antibodies (mAbs), collectively denoted here as antibody-based treatments (ABT).

However, there does not seem to be a consistent pattern that describes the evolution of SARS-CoV-2 across all chronic infections. Whereas some cases display dramatic evolution in the spike (s) protein (e.g., 8), in other chronic infections relatively limited evolution is observed (e.g., 14). Here, we set out to consolidate the evolutionary patterns found across chronic infections, by re-analyzing previous reports, and by sequencing a cohort of six patients with chronic SARS-CoV-2 infection patients at the Tel Aviv Sourasky Medical Center 15,16. We explore the shared and unique mutational patterns that emerge in different chronically infected patients, and compare them to those observed in global data reflecting transmission chains in acute individuals. We focus on inferring the conditions that lead to dramatic evolution and the potential for the creation of new VOCs.

Results

Different patterns of evolution detected in chronically infected patients

We begin by defining criteria for a chronic infection. In clinical settings, a chronic infection is often defined as one with both prolonged shedding of viral RNA, coupled with evidence for infectious virus (either through virus isolation in tissue culture, or via detection of sub-genomic RNA). However, when surveying various studies reporting chronic infection, a lack of standardization was noted. We hence expanded our focus to include patients displaying high-viral-load shedding for 20 or more days, while mining the literature for all such cases that were accompanied by longitudinal whole-genome sequencing of the virus (Methods). The criterion of 20 days was based on a meta-analysis of the duration of viral shedding (defined as a positive nasopharyngeal PCR test) across thousands of patients diagnosed up till June 2020, which revealed that mean duration of upper-respiratory tract shedding was around 17 days, with a 95% confidence interval ranging from 15.5 to 18 days 17. Of note, shedding of replication-competent virus lasted a lot less than twenty days. Moreover, estimates of viral shedding are different in some of the more recently detected SARS-CoV-2 variants such as Delta and Omicron 18,19, yet as described below, our analysis focused on variants that were found in earlier stages of the pandemic.

Our search yielded a total of 21 case reports, all of which were diagnosed during 2020 or early 2021, and all of which were infected with viruses belonging to lineages that predated the Alpha variant (Table S2). In addition, six patients adhering to the above criterion were detected in Tel Aviv Sourasky Medical Center (TASMC), and all available samples were sequenced (Methods). Five TASMC patients suffered from hematological cancers. The sixth patient suffered from an autoimmune disorder, and was treated with a high dose of steroids. The six TASMC patients were all diagnosed in late 2020 or early 2021, with four patients infected with a virus from pre-Alpha lineages, and two patients infected with a virus from the Alpha lineage (Table S2).

Overall, when observing the entire set of 27 chronically infected patients, we observed that all were immunosuppressed, due to one or more of the following: hematological cancer, anti-B cell treatment, high dosage steroid treatment, or very low CD4 counts (acquired immunodeficiency disease, AIDS). We observed very different evolutionary outcomes across the range of patients examined, from dramatic evolution and antibody-evasion observed in some patients, to relatively static evolution in others (Table 1, Tables S1-S2).

View this table:
Table 1. Summary of all 27 patients with chronic SARS-CoV-2 infections.

The last six patients were sequenced herein. See Tables S1,S2 for more details.

Comparing patterns of evolution between chronic infections and global transmission chains

We searched for patterns of evolution across the total of twenty-seven patients with chronic infection, and compared this pattern with the pattern observed (a) under mostly neutral evolution, in the first approximately nine months of viral circulation 28,29, and (b) under presumable positive selection, which occurred in the lineages leading to the five currently defined VOCs (Alpha, Beta, Gamma, Delta, Omicron) (Fig. 1, Table S4). In each scenario, we searched for bins of mutations (along consecutive regions of 500 bases) enriched for mutations (p<0.05, binomial test, after correction for multiple testing; Methods).

Figure 1.
Figure 1. Number of substitutions observed along the SARS-CoV-2 genome, in bins of 500 nt.

The upper panel displays substitutions observed at any time-point of the twenty-seven chronic infections. The middle panel displays lineage-defining mutations of the five currently recognized VOCs. The lower panel displays substitutions observed globally during the first nine months of the pandemic, mostly before the emergence of VOCs. Asterisks and dots mark bins enriched for more substitutions (P<0.05 and P<0.1, respectively, Methods). The genomic positions are based on the Wuhan-Hu-1 reference genome (GenBank ID NC_045512), and the banner on the top shows a breakdown of ORF1a/b into individual proteins, and domains for the spike protein (see main text).

Figure 2.
Figure 2. A network of co-occurring substitutions across patients with chronic SARS-CoV-2 infection.

Each colored circle represents a locus, and an asterisk and dot represent P<0.05 and P<0.1, respectively.

During the first nine months of circulation, we noted that 65% of substitutions were non-synonymous, which is generally what we could expect under lack of both positive and purifying selection, and in line with reports suggesting incomplete purifying selection during the early stages of SARS-CoV-2 spread 30. During this time, we observed a relatively uniform distribution of substitutions across most of the genome, with some enrichment in ORF3a, ORF7, ORF8, and in N. This enrichment has been previously reported and may be due to more relaxed purifying selection in these regions, or higher mutation rates 28; adaptive evolution at these regions cannot be ruled out as well. We also noted some enrichment during the first nine months of the pandemic in the N-terminal domain (NTD) of the spike, which appeared to be mostly confined to the very first residues of the protein. These residues belong to the “NTD supersite” and are associated with immune evasion 31,32, which seems unlikely at this early stage of the pandemic, yet the NTD is also associated with structural rearrangements following ACE2 receptor binding33.

In general, the patterns obtained in chronic infections and in VOCs were very similar. The average proportion of non-synonymous substitutions was 78% and 82% in chronic infections and VOCs, respectively, generally suggestive of positive selection operating on both. The most striking similarity was observed in patterns along the Spike protein, and in particular at the regions which correspond to the NTD (genomic nucleotides 21,598 to 22,472) and the receptor binding domain (RBD) (genomic nucleotides 22,517 to 23,183). Several mutations at the RBD have been shown to enhance affinity to the ACE2 receptor and allow for better replication 34,35, whereas other mutations, both at RBD and NTD, are known to enhance antibody-evasion 31,36,37. Notably, chronic infections exhibited mutations both in the first NTD residues observed in the “neutral” pattern described above, but also in later positions. The most commonly observed substitutions in chronic infections were in the spike protein: E484K/Q and various deletions in the region spanning the NTD supersite, particularly amino-acids 140 through 145, all shown to confer antibody-evasion 38. Chronic infections shared the enrichment of ORF3a/ORF7a/ORF8 mutations with the “neutral” set, but lacked an enrichment across most of the N protein. Overall, it seems that mutations in chronic infections are predictive of those occurring in VOCs, as has been noted previously 2.

When focusing on the differences between VOCs and chronic infections, several intriguing differences emerged. First, four VOCs bear a three-amino acid deletion in the nsp6 gene (ORF1a: Δ3675-3677), an event not observed in our set of chronic infections. Next, in VOCs there is an enrichment in the region of the spike encompassing the S1/S2 boundary (positions 23,500-24,000 in Fig.1). This enrichment is primarily driven by S: P681H/R, a highly recurrent globally occurring mutation 39, surprisingly never observed in our chronic set. A recent analysis has analyzed the relationship between recurrent mutations and their success in driving clade expansion, i.e., onwards transmission 40. Notably, we observed that successful recurrent mutations were almost never observed in our chronic set, whereas non-successful recurrent mutations (S:E484K/Q and S:Δ144) were the ones most abundantly observed (Table 2). Overall, these results suggest that there may be a trade-off between antibody-evasion and transmissibility. This trade-off, if it exists, would not play a role in chronic infections, but would impact the ability of a variant created in a chronic infection to be transmitted onwards. Thus, only under specific conditions, a transmissible variant would emerge in chronic infections. Four of five VOCs independently acquired a mutation at or near the S1/S2 boundary (S:P681H/R or H655Y) suggesting that this may be a factor driving transmissibility. We note that Beta is an exception with no such mutations, yet this variant also displayed very limited global transmission.

View this table:
Table 2. Recurrent mutations observed along the SARS-CoV-2 phylogeny.

Clade success or lack of success is reported based on measurements of clade logistic growth, with + or – representing higher or lower than average growth, respectively 40. The last column marks the number of patients in the set of chronic infections herein where a substitution was observed.

We went on to examine co-occurring substitutions, defined as pairs of substitutions that appeared in two or more patients. We used a chi-square test to assess whether pairs of substitution occurred together more often than expected from their individual frequencies (Methods), as a measure of possible epistasis. Intriguingly, four pairs of substitutions across four different proteins emerged as significantly enriched, and formed a network of interactions: T30I in envelope, H125Y in the membrane glycoprotein, S13I in the spike protein, and T3058I in ORF1a (Fig. 4). This finding was intriguing on multiple fronts: first, envelope and membrane glycoprotein have generally remained very conserved throughout the entire pandemic, and specifically the two replacements found are at highly conserved sites (Table S1). The replacements in spike and ORF1a, on the other hand, have been observed a small number times in the global phylogeny. Notably, all the first three proteins form a part in the virion structure itself – however the functional meaning of this remains unclear. Other pairs of mutations found to co-occur were the three most common spike antibody-evasion mutations, yet these co-occurrences were not statistically significant. Larger cohorts of patients and further data will be required to determine the implications of these findings.

Figure 3.
Figure 3. Viral rebound predicts mutations in the spike protein associated with immune evasion.

Ct values are presented according to the day of infection (denoted as number of days post the first positive PCR test), with the dashed red horizontal line and shaded area representing a negative or borderline result, respectively. Blue dots represent samples that were sequenced. Only amino-acid replacements in the spike protein are shown, with predicted immune evasion mutations shown in bold (Table S1). Positive samples from bronchoalveolar lavage (BAL), endotracheal aspirates (ETA) or sputum are indicated in brown. Antibody-based anti-COVID19 treatments are represented by dashed vertical lines on the day of administration.

Figure 4.
Figure 4. Polymorphic populations observed across patients.

Each series of boxed lines represents a patient, and each line represents a sequenced time-point with time since infection on the right. The different open reading frames are color-coded. For each patient, only mutations relative to the first time point sequenced that appeared at a frequency ranging from 20% to 100% are shown.

Predictors of antibody-evasion substitutions

We went on to examine the variation observed across different patients. When examining medical background, the patients could be roughly classified into one of the following categories: hematological cancers, HIV/AIDS, organ transplantation, and auto-immune disorders (Table 1). The latter two categories were often treated with steroids. Some but not all of the hematological cancer patients and others were treated with antibodies targeting B-cells, presumably causing profound B-cell depletion. In line with this, most of the patients with confirmed B-cell depletion showed negative serology for SARS-CoV-2 at one or more time-points (Table S1). Some patients were treated with ABT against SARS-CoV-2, whereas others were not, and in some antibody-evasion mutations were detected while in others they were not. Finally, we found that while in some ABT patients, immune evasion mutations were detected, sometimes these mutations fixed before the treatment (see e.g., Fig. 1, P5).

We noted that many patients (four of the six patients sequenced herein) displayed an intriguing cycling pattern of viral load (VL) (reflected by cycle threshold Ct values), with very high Ct values reaching negative or borderline negative results at one or more stages of the infection, followed by rebound of the virus. In the four abovementioned patients, this rebound was accompanied by clinical evidence for disease, highly suggestive of active viral replication. Several different hypotheses could explain this pattern. First, the virus may have cleared and been followed by re-infection with another variant. Since this pattern can be ruled out using sequencing, such cases were excluded from our analysis (Methods). Second, the virus may cycle between different niches, such as upper and lower airways. Its re-emergence in the upper airways (nasopharynx) may be due to selective forces or genetic drift. When considering selective forces, viral rebound may occur due to the near clearance of the virus, either via ABT or by the self-immune system, followed by the emergence of a more fit variant with antibody-evasion properties.

We set out to perform logistic regression to predict an outcome of antibody-evasion (Methods). We first tested if B-cell depletion or steroid treatment predict antibody-evasion mutations sometime during the course of infection (Table S2), and found this was not a significant predictor (P=0.6, P=0.4, N.S., respectively). Moreover, the length of sequenced infection (i.e., the maximal time-point sequenced) was not predictive of antibody-evasion either (P=0.9, N.S.). Next, we went on to test whether ABT predicts antibody-evasion, and if antibody-evasion tends to be preceded by a drop in viral load (Table S3). Our results showed that ABT is not a significant predictor (P=0.17, NS), whereas a drop in viral load followed by viral rebound is (P=2.6E-05) (Fig. 1). These results suggest that in some cases, what may be driving immune escape is actually the weakened immune system of the patient, although the somewhat low P-value of ABT suggests it is possible that the treatment itself also plays a role in driving antibody-evasion as well. This result suggests that the adaptive immune response of the patients herein is not completely absent, and that a weak antibody response is what drives adaptive evolution in the virus.

Next, we went on to examine patterns of variation over time across the different patients. In many of the case reports, the authors noted the emergence and disappearance (and sometimes re-emergence) of particular substitutions (Fig. 4). When re-analyzing the data, we noted that this pattern of variation across time was observed in the majority of patients (Table S2). From an evolutionary point of view, it is quite unlikely for one or more substitutions to disappear from a given population, and since we observe this at very different loci across all patients, it is not likely that this is due to recurrent sequencing problems or due to biases of the viral polymerase. We and others have previously noted sequencing errors that occur predominantly when viral load is low, when errors that occur during reverse transcription or early PCR cycles are carried over to higher frequencies 10,11,41. However, this phenomenon most often leads to errors in intra-host variants segregating at relatively low frequency, and is less common at the consensus sequence level, which is defined here as mutations present at a frequency of 80% or higher. We thus conclude that the most likely explanation for the highly dynamic polymorphisms observed, is the existence of distinct populations residing in different niches, which may correspond to different infected organs, or niches within an infected organ.

Discussion

Tying the results obtained so far together, we suggest that viral rebound coupled with antibody-evasion may be enabled by partial penetration of antibodies to a specific niche. Accordingly, this partial penetration will prevent viral clearance on one hand, but on the other will promote selection for an antibody-evasion mutation. Notably, this is supported by the case of patient P1 in this study: during the time a negative nasopharyngeal swab was obtained, a BAL sample returned as positive, just prior to the emergence of the E484Q antibody evasion mutation. This suggests that the virus continued to thrive in lungs, possibly as a variant adapted to lungs, or as the same nasopharyngeal variant that migrated to the lungs. Alternatively, it is also possible that the immune evasion mutants continued to thrive in the nasopharyngeal swabs but were not picked up by the sample for PCR.

We go on to discuss the conditions that enable a variant in a chronic infections to keep spreading in the general population, and to potentially become a variant of concern. First, we searched for evidence of onwards transmission from a chronically infected patient in the cohort studied herein, yet found no such evidence in the publications or from epidemiological investigations of our TASMC patients. This is not necessarily surprising, as often such patients are bedridden, and in strict isolation both due to COVID-19 and due to their background conditions. However, this does raise the possibility that (a) chronic infections harbor variants that are non-infectious, and (b) chronic infections harbor non-transmissible virus, options we discuss below.

We first consider the option that chronic infections harbor non-infectious virus. In twelve of twenty-seven patients, evidence for viral replication was obtained, in the form of culturable virus, or in the form of PCR positive for sub-genomic RNA. In an additional eight patients, we noted either consistently high viral loads (Ct values around 20 or lower) or a negative difference of ten cycles between one PCR outcome and a later one (Table S6). In other words, in the latter scenario viral load increased dramatically across time in these patients, strongly suggestive of active viral replication. In line with this, we did not find evidence for mutations associated with defective virus, such as premature stop codons, frameshifts, or dramatic deletions outside of the spike protein (although in-frame deletions were noted in two patients, in ORF1a; Table S1). Thus, to summarize, while we cannot completely rule out the presence of non-infectious defective virus in chronic infections, we consider it is unlikely that all virus is non-infectious in most of the cases examined herein.

We next go on to discuss the probability that variants found in chronic infections are transmissible. In most of our patients (15/27) we observed either S:E484K/Q and/or deletions at S:Y144 in the NTD, both of which are highly recurrent mutations in the global phylogeny. However, globally circulating variants bearing these mutations were inferred to be less transmissible than other highly recurrent mutations (Table 2), and E484K/Q is notably absent in the once globally dominant variants Alpha and Delta. On the other hand, several mutations associated with higher transmissibility (e.g., S:P681R, ORF1a: Δ3675-3677) 40, were never observed in the chronic infections herein (Table 2).

Why if so were these mutations never observed? Evidently, in chronically infected patients we do not expect selection for mutations enhancing transmissibility, but rather expect selection to promote mutations that enhance viral replication (with immune evasion included in this latter category). However, there is often a link between better net viral replication and higher transmissibility 42-44, which should theoretically promote transmissibility-enhancing mutations also in in chronic infections. One possibility is that the sample size of patients in this study was too small. Given the large number of time we observed S:E484K/Q in chronic infections, it seems more likely that there is another explanation. Another possibility is that for example, S:P681H/R is under weaker positive selection in the conditions particular to immunocompromised patients. A recent study has suggested that part of the evolutionary advantage of S:P681H is its ability to evade the interferon response of the cell 45. Perhaps the particular conditions in immunocompromised prevent selection for interferon evasion: for example, some patients were treated with interferon inhibitors, although these are a minority of patients (Table S2); alternatively, the interferon response of immunocompromised patients is also likely weakened per se.

Finally, an additional possibility is that epistatic interactions among mutations prevent the emergence of some mutations, or in other words, two adaptive mutations may become deleterious when residing on the same genome. In line with this, the combination S: E484K/Q and S: P681H/R is not a common one, except for in Omicron (albeit E484A), where over 30 Spike amino-acid replacements have been observed. Indeed, several reports have suggested abundant epistatic interactions led to the emergence of Omicron 46, and these may allow the co-existence of S: E484A with S: P681H and others. It remains an enigma where Omicron first arose 47, yet the fact we see potentially epistatic interactions occurring in some chronic infections, coupled with the fact Omicron is characterized by many such pairs of mutations, suggest that lengthy chronic infections may allow an exploration of the SARS-CoV-2 fitness landscape, such that some of these chronic infections lead to the emergence of a highly transmissible variant.

Conclusions

To summarize, we have found that the overall patterns of mutations observed in chronic infections closely mirror the pattern observed in VOCs, with some notable exceptions. We find that when viral rebound is observed, this strongly suggests that antibody-evasion mutations may have emerged in a COVID-19 patients. Thus, viral rebound can be viewed as a warning signal that VOC-like mutation occurred in a patient, and extra caution may be warranted: genetic sequencing, isolation, and close monitoring of contacts may be crucial when viral rebound is observed. While these immune-evasion mutations are of concern, we suggest that most variants emerging in chronically infected patients lack the potential for dramatic onwards transmission, possibly due to an absence of key mutations. Future research is necessary to understand the precise factors determining when and if a variant generated in chronic infection becomes highly transmissible. Overall, we suggest that extensive monitoring of chronic infections can be used for forecasting the set of mutations in future VOCs.

Methods

TASMC virus genome sequencing

Leftover RNA from nasopharyngeal swabs were obtained from six chronically infected individuals, across several different time points for each patients (Table S1). Patient 2 sequences were available from our previous study16; patients 1 and 3 were also previously sequenced 15, but here we added on additional time-points. All patients were positive via COVID-19 RT-qPCR tests of SARS-CoV-2 for a period longer than 3 weeks, and all were immunocompromised. All relevant clinical data was collected and was summarized in Tables 1, S2. All samples underwent whole-genome SARS-CoV-2 using the V3 ARTIC protocol (https://artic.network/ncov-2019) and the samples were sequenced using Illumina Miseq 250-cycle V2 kits at the Technion Genomic Center (Israel). This study was approved by the TASMC Helsinki committee (approval No. 1042-20-TLV), and by the Tel Aviv University IRB (approval number 0004435-1), with an exemption from informed consent.

Determining genome consensus sequences and lineage determination

pTrimmer was used in order to trim the various primers that were used in the multiplex PCR (https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2854-x). The raw reads were mapped to the reference genome of SARS-CoV-2 (GenBank ID MN908947). Mapping and variant calling was performed using our AccuNGS pipeline 41. The consensus sequence of each sample was determined based on sites with coverage of at least 10X sites. As we and other have observed, the sequencing process may lead to erroneous mutations being observed, and this is particularly exacerbated at low viral loads 10,11,41. To this end, only mutations with a frequency of at least 80% were introduced into the consensus sequence (hereby referred to as substitutions or fixed mutations); loci with mixed populations at frequencies lower than 80% were considered ambiguous and marked by “N”.

For each patient, substitutions as compared to the reference that appeared in the first time point were marked as background substitutions, and were not considered for any further analyses. Only substitutions that were added in the following time points were considered. The “Pangolin” network was used to identify the consensus sequence lineage (https://cov-lineages.org/resources/pangolin.html). This allowed us to rule out re-infection in all patients but one: the last time-point of P6 from day 48 of infection was inferred to be due to re-infection (B.1.1.7 as compared to B.1.1.50 in all earlier time-points). The last time point from this patients was removed from the analysis, leaving a total of 37 days. Visual inspection of the results led to a suspected contamination at day 27 of P5, as several B.1.1.7 lineage defining mutations were observed on an otherwise B.1.1.50 background. This sample was omitted from the analysis.

Data availability

All consensus sequences for the patients sequenced herein were deposited in GISAID, and raw sequencing reads were uploaded to the NCBI Sequence Read Archive under submission number PRJNA803960 (Table S5).

Reanalysis of previous case reports

Previously published case reports were selected by mining the literature for the key-words “prolonged”, “persistent” or “chronic” “SARS-CoV-2 infection” up till Nov. 15 2021. While some studies validated the chronic nature of infection by performing culturing of the virus or by measuring sub-genomic RNA (sgRNA), several did not. Thus, when culture/sgRNA data was unavailable, we filtered out patients with infections of less than 20 days, based on an a meta-analysis that showed viral shedding was mostly limited to 20 days (see main text). We verified intermediate to high viral load (Ct values of 27 or lower) at some time exceeding day 20 as a proxy for ensuring that shedding was not of residual non-infectious virus. Only studies that performed longitudinal whole genome sequencing were maintained. In Tarhini. et al.25 patient 3 was removed from the analysis as the authors identified a co-infection event. Sequencing data for Jensen et al.22 was obtained from the authors, albeit with missing data – the first time point for patient D was missing. We thus relied on the report by the authors on S:E484 at time point 1, and only mutation S:E484K was considered in time point 2 for patient D.

When analyzing the data, and in line with the above, the focus of almost all our analyses, was on mutations with a frequency higher than 80%, not present in the first time point sequenced. When analyzing viral polymorphisms across time (Fig. 4) we also present mutations at lower frequencies. All relevant clinical data was collected and was summarized in Tables 1 and S2. Ct values for RT-qPCR tests were also extracted from all publications (Table S6). A test was considered borderline if the Ct ranged from 34 to 37, and was considered negative if reported as such, or was equal to 38 or higher. In one study 22, raw viral load measurements were presented in terms of number of copies per ml sample, and a negative result was reported when zero copies were detected.

Summary information of the patients

All available medical records including the patient’s background condition, treatments and anti SRAS-CoV-2 serology test are specified in Table S2. Our cohort composed of many immunocompromised conditions: seven-teen of the patients suffer from different kinds of lymphomas and leukemias, four of the patients undergo organ transplant, tree suffer from AIDS, two of the patients have an autoimmune disease and one suffer from renal disease. Regarding age groups: ten of the patients are classified in age range of 60-69, five belong to 70-79, and tree patients to 40-49 years old. In each of the age ranges of 0-9, 20-29, 30-39 and 50-59 two patients were added and only one patient is in age range of 80-89 years old.

Summary of substitution events

A summary of all substitution events was compiled into three tables. Table S1 lists every substitution and the patient/s where it was found. Data on the proportion of a given substitution in the general global viral samples was taken from http://cov-glue-viz.cvr.gla.ac.uk/, and counts of substitutions along a sampled phylogeny of approximately 3,500 sequences were taken from NextStrain (www.nextrain.org) 48 accessed on Jan 2, 2022. Each of the spike amino-acid replacements was queried to see if information exists on antibody-evasion for a given mutation, and was marked as enhancing antibody-evasion if a publication was available showing direct experimental evidence supporting this (Table S1).

Table S2 lists information (clinical and mutational, across all time-points) per patient, whereas Table S3 lists information per time point sequenced for each patient. When two samples from the same time point were available from different niches, we presented only the nasopharyngeal one, since our focus is primarily on this niche. Counts and information of substitutions per patient are listed in table S2 for both the spike region and the rest of the genome. Transient substitutions (that re-appeared at a certain time point) were counted only once.

Substitution events histogram

We counted the total number of substitutions observed at each position in the genome in three sets of sequences:

  1. our set of chronic infections

  2. the lineage defining mutation of the five VOCs (Alpha, Beta, Gamma, Beta, Omicron), as defined by https://covariants.org. Substitutions shared by two or more VOCs were omitted.

  3. a subset of sequences from NextStrain (www.nextrain.org)48 accessed on Jan 2, 2022. Sequences were filtered to only those earlier than 16/8/2020. The numbers of events per genome position were downloaded from the site.

A histogram was created for each of the three sets, with bin sizes of 500 bases. The probability of observing a given number of substitutions in a bin was calculated based on a binomial distribution with p=500/L (L being the reference genome length), n=the sum of all substitutions observed in a set, and k=observed number of substitutions in a given bin. Correction for multiple testing was performed using the false discovery rate (FDR) of Benjamini–Hochberg49.

Logistic regression

Logistic regression was performed twice to test for significant predictors of an antibody-evasion outcome. First, at the level of each patient, we tested whether severe B-cell depletion, steroid treatment, or the maximal time-point since infection that was sequenced. Severe B-cell depletion was defined as one of the following conditions: treatment with anti-B cell directed antibodies up to nine months prior to COVID-9 detection, or blood-test measurements of low B-cell markers. Second, at the level of time-point sequenced, we tested whether ABT (convalescent plasma, hyperimmune plasma, or monoclonal antibody/ies) or viral rebound predicted an outcome of an antibody-evasion mutation. Each sequenced time point was hence characterized as ABT-true if treatment was given up to seven days before the treatment, and viral rebound was defined as true if the COVID-19 PCR test just prior to the sequenced time point was borderline or negative. Binomial logistic regression was performed using the glm package in R. All the raw data is available in tables S2 and S3.

Co-occurring substitutions

We began by identifying all pairs of substitutions that occurred in two or more patients. A χ2 test was used in order to identify a deviation from expected frequencies, by comparing marginal probabilities with paired probabilities (number of degrees of freedom was 1). Multiple testing was accounted for via FDR 49.

Data Availability

All consensus sequences for the patients sequenced herein were deposited in GISAID, and raw sequencing reads were uploaded to the NCBI Sequence Read Archive under submission number PRJNA803960. Full details are available in the manuscript supplementary material.

Acknowledgements

We would like to thank Shay Fleishon, Neta Zuckerman, Tzachi Hagai, David Burstein, and Noam Harel for critical reading and discussions. This study was supported by an ERC starting grant 852223 (RNAVirFitness), and by an Israeli Science Foundation grant 3963/19. This study was supported in part by a fellowship to SH from the Edmond J. Safra Center for Bioinformatics at Tel-Aviv University. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

References

  1. ↵
    Fu, Y. et al. Risk factors for viral RNA shedding in COVID-19 patients. Eur Respir J 56, doi:10.1183/13993003.01190-2020 (2020).
    OpenUrlAbstract/FREE Full Text
  2. ↵
    Corey, L. et al. SARS-CoV-2 Variants in Patients with Immunosuppression. The New England journal of medicine 385, 562–566, doi:10.1056/NEJMsb2104756 (2021).
    OpenUrlCrossRefPubMed
  3. ↵
    Moran, E. et al. Persistent SARS-CoV-2 infection: the urgent need for access to treatment and trials. Lancet Infect Dis 21, 1345–1347, doi:10.1016/S1473-3099(21)00464-3 (2021).
    OpenUrlCrossRef
  4. ↵
    Sudre, C. H. et al. Attributes and predictors of long COVID. Nat Med 27, 626–631, doi:10.1038/s41591-021-01292-y (2021).
    OpenUrlCrossRefPubMed
  5. ↵
    Fung, M. & Babik, J. M. COVID-19 in Immunocompromised Hosts: What We Know So Far. Clin Infect Dis 72, 340–350, doi:10.1093/cid/ciaa863 (2021).
    OpenUrlCrossRef
  6. ↵
    Borges, V. et al. Long-Term Evolution of SARS-CoV-2 in an Immunocompromised Patient with Non-Hodgkin Lymphoma. mSphere 6, e0024421, doi:10.1128/mSphere.00244-21 (2021).
    OpenUrlCrossRef
  7. Kemp, S. A. et al. SARS-CoV-2 evolution during treatment of chronic infection. Nature 592, 277–282, doi:10.1038/s41586-021-03291-y (2021).
    OpenUrlCrossRef
  8. ↵
    Choi, B. et al. Persistence and Evolution of SARS-CoV-2 in an Immunocompromised Host. The New England journal of medicine 383, 2291–2293, doi:10.1056/NEJMc2031364 (2020).
    OpenUrlCrossRefPubMed
  9. ↵
    Lythgoe, K. A. et al. SARS-CoV-2 within-host diversity and transmission. Science 372, doi:10.1126/science.abg0821 (2021).
    OpenUrlAbstract/FREE Full Text
  10. ↵
    Tonkin-Hill, G. et al. Patterns of within-host genetic diversity in SARS-CoV-2. eLife 10, doi:10.7554/eLife.66857 (2021).
    OpenUrlCrossRef
  11. ↵
    Valesano, A. L. et al. Temporal dynamics of SARS-CoV-2 mutation accumulation within and across infected hosts. PLoS pathogens 17, e1009499, doi:10.1371/journal.ppat.1009499 (2021).
    OpenUrlCrossRef
  12. ↵
    Braun, K. M. et al. Acute SARS-CoV-2 infections harbor limited within-host diversity and transmit via tight transmission bottlenecks. PLoS pathogens 17, e1009849, doi:10.1371/journal.ppat.1009849 (2021).
    OpenUrlCrossRef
  13. ↵
    Rambaut, A. & al., e. Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations, <https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563> (2020).
  14. ↵
    Baang, J. H. et al. Prolonged Severe Acute Respiratory Syndrome Coronavirus 2 Replication in an Immunocompromised Patient. The Journal of infectious diseases 223, 23–27, doi:10.1093/infdis/jiaa666 (2021).
    OpenUrlCrossRefPubMed
  15. ↵
    Bronstein, Y. et al. Evolution of spike mutations following antibody treatment in two immunocompromised patients with persistent COVID-19 infection. Journal of medical virology, doi:10.1002/jmv.27445 (2021).
    OpenUrlCrossRef
  16. ↵
    Meijer, S. E. et al. Unraveling a Nosocomial Outbreak of COVID-19: The Role of Whole-Genome Sequence Analysis. Open Forum Infectious Diseases 8, doi:10.1093/ofid/ofab120 (2021).
    OpenUrlCrossRef
  17. ↵
    Cevik, M. et al. SARS-CoV-2, SARS-CoV, and MERS-CoV viral load dynamics, duration of viral shedding, and infectiousness: a systematic review and meta-analysis. Lancet Microbe 2, e13–e22, doi:10.1016/S2666-5247(20)30172-5 (2021).
    OpenUrlCrossRefPubMed
  18. ↵
    Hay, J. A. et al. Viral dynamics and duration of PCR positivity of the SARS-CoV-2 Omicron variant. medRxiv, 2022.2001.2013.22269257, doi:10.1101/2022.01.13.22269257 (2022).
    OpenUrlAbstract/FREE Full Text
  19. ↵
    Wang, Y. P. et al. Transmission, viral kinetics and clinical characteristics of the emergent SARS-CoV-2 Delta VOC in Guangzhou, China. Eclinicalmedicine 40, doi:ARTN 101129 10.1016/j.eclinm.2021.101129 (2021).
  20. Avanzato, V. A. et al. Case Study: Prolonged Infectious SARS-CoV-2 Shedding from an Asymptomatic Immunocompromised Individual with Cancer. Cell 183, 1901–1912 e1909, doi:10.1016/j.cell.2020.10.049 (2020).
    OpenUrlCrossRefPubMed
  21. Chen, L. et al. Emergence of Multiple SARS-CoV-2 Antibody Escape Variants in an Immunocompromised Host Undergoing Convalescent Plasma Treatment. mSphere 6, e0048021, doi:10.1128/mSphere.00480-21 (2021).
    OpenUrlCrossRef
  22. ↵
    Jensen, B. et al. Emergence of the E484K mutation in SARS-COV-2-infected immunocompromised patients treated with bamlanivimab in Germany. Lancet Reg Health Eur 8, 100164, doi:10.1016/j.lanepe.2021.100164 (2021).
    OpenUrlCrossRef
  23. Karim, F. et al. Persistent SARS-CoV-2 infection and intra-host evolution in association with advanced HIV infection. medRxiv, 2021.2006.2003.21258228, doi:10.1101/2021.06.03.21258228 (2021).
    OpenUrlAbstract/FREE Full Text
  24. Leung, W. F. et al. COVID-19 in an immunocompromised host: persistent shedding of viable SARS-CoV-2 and emergence of multiple mutations: a case report. Int J Infect Dis 114, 178–182, doi:10.1016/j.ijid.2021.10.045 (2022).
    OpenUrlCrossRef
  25. ↵
    Tarhini, H. et al. Long-Term Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infectiousness Among Three Immunocompromised Patients: From Prolonged Viral Shedding to SARS-CoV-2 Superinfection. The Journal of infectious diseases 223, 1522–1527, doi:10.1093/infdis/jiab075 (2021).
    OpenUrlCrossRef
  26. Truong, T. T. et al. Increased viral variants in children and young adults with impaired humoral immunity and persistent SARS-CoV-2 infection: A consecutive case series. EBioMedicine 67, 103355, doi:10.1016/j.ebiom.2021.103355 (2021).
    OpenUrlCrossRefPubMed
  27. Perez-Lago, L. et al. Different Within-Host Viral Evolution Dynamics in Severely Immunosuppressed Cases with Persistent SARS-CoV-2. Biomedicines 9, doi:10.3390/biomedicines9070808 (2021).
    OpenUrlCrossRef
  28. ↵
    Rochman, N. D. et al. Ongoing global and regional adaptive evolution of SARS-CoV-2. Proc Natl Acad Sci U S A 118, doi:10.1073/pnas.2104241118 (2021).
    OpenUrlAbstract/FREE Full Text
  29. ↵
    MacLean, O. A. et al. Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen. PLoS Biol 19, e3001115, doi:10.1371/journal.pbio.3001115 (2021).
    OpenUrlCrossRef
  30. ↵
    Ghafari, M. et al. Purifying selection determines the short-term time dependency of evolutionary rates in SARS-CoV-2 and pH1N1 influenza. Mol Biol Evol, doi:10.1093/molbev/msac009 (2022).
    OpenUrlCrossRef
  31. ↵
    Harvey, W. T. et al. SARS-CoV-2 variants, spike mutations and immune escape. Nat Rev Microbiol 19, 409–424, doi:10.1038/s41579-021-00573-0 (2021).
    OpenUrlCrossRefPubMed
  32. ↵
    Lok, S. M. An NTD supersite of attack. Cell host & microbe 29, 744–746, doi:10.1016/j.chom.2021.04.010 (2021).
    OpenUrlCrossRef
  33. ↵
    Benton, D. J. et al. Receptor binding and priming of the spike protein of SARS-CoV-2 for membrane fusion. Nature 588, 327–330, doi:10.1038/s41586-020-2772-0 (2020).
    OpenUrlCrossRefPubMed
  34. ↵
    Barton, M. I. et al. Effects of common mutations in the SARS-CoV-2 Spike RBD and its ligand, the human ACE2 receptor on binding affinity and kinetics. eLife 10, doi:10.7554/eLife.70658 (2021).
    OpenUrlCrossRef
  35. ↵
    Starr, T. N. et al. Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding. Cell 182, 1295–1310 e1220, doi:10.1016/j.cell.2020.08.012 (2020).
    OpenUrlCrossRefPubMed
  36. ↵
    Greaney, A. J. et al. Complete Mapping of Mutations to the SARS-CoV-2 Spike Receptor-Binding Domain that Escape Antibody Recognition. Cell host & microbe 29, 44–57 e49, doi:10.1016/j.chom.2020.11.007 (2021).
    OpenUrlCrossRefPubMed
  37. ↵
    Starr, T. N. et al. Prospective mapping of viral mutations that escape antibodies used to treat COVID-19. Science 371, 850-+, doi:10.1126/science.abf9302 (2021).
    OpenUrlAbstract/FREE Full Text
  38. ↵
    Chen, R. E. et al. Resistance of SARS-CoV-2 variants to neutralization by monoclonal and serum-derived polyclonal antibodies. Nat Med 27, 717–726, doi:10.1038/s41591-021-01294-w (2021).
    OpenUrlCrossRefPubMed
  39. ↵
    Tao, K. et al. The biological and clinical significance of emerging SARS-CoV-2 variants. Nat Rev Genet 22, 757–773, doi:10.1038/s41576-021-00408-x (2021).
    OpenUrlCrossRefPubMed
  40. ↵
    Kistler, K. E., Huddleston, J. & Bedford, T. Rapid and parallel adaptive mutations in spike S1 drive clade success in SARS-CoV-2. bioRxiv, 2021.2009.2011.459844, doi:10.1101/2021.09.11.459844 (2021).
    OpenUrlAbstract/FREE Full Text
  41. ↵
    Gelbart, M. et al. Drivers of within-host genetic diversity in acute infections of viruses. PLoS pathogens 16, e1009029, doi:10.1371/journal.ppat.1009029 (2020).
    OpenUrlCrossRefPubMed
  42. ↵
    Zhou, B. et al. SARS-CoV-2 spike D614G change enhances replication and transmission. Nature 592, 122–127, doi:10.1038/s41586-021-03361-1 (2021).
    OpenUrlCrossRefPubMed
  43. Lipsitch, M. & Moxon, E. R. Virulence and transmissibility of pathogens: What is the relationship? Trends in Microbiology 5, 31–37, doi:Doi 10.1016/S0966-842x(97)81772-6 (1997).
    OpenUrlCrossRefPubMedWeb of Science
  44. ↵
    Mlcochova, P. et al. SARS-CoV-2 B.1.617.2 Delta variant replication and immune evasion. Nature 599, 114–119, doi:10.1038/s41586-021-03944-y (2021).
    OpenUrlCrossRefPubMed
  45. ↵
    Lista, M. J. et al. The P681H mutation in the Spike glycoprotein confers Type I interferon resistance in the SARS-CoV-2 alpha (B.1.1.7) variant. bioRxiv, 2021.2011.2009.467693, doi:10.1101/2021.11.09.467693 (2021).
    OpenUrlAbstract/FREE Full Text
  46. ↵
    Martin, D. P. et al. Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function. bioRxiv, 2022.2001.2014.476382, doi:10.1101/2022.01.14.476382 (2022).
    OpenUrlAbstract/FREE Full Text
  47. ↵
    Wei, C. et al. Evidence for a mouse origin of the SARS-CoV-2 Omicron variant. J Genet Genomics, doi:10.1016/j.jgg.2021.12.003 (2021).
    OpenUrlCrossRef
  48. ↵
    Hadfield, J. et al. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics 34, 4121–4123, doi:10.1093/bioinformatics/bty407 (2018).
    OpenUrlCrossRefPubMed
  49. ↵
    Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B (Methodological), 289–300 (1995).
Back to top
Posted February 28, 2022.
Download PDF

Supplementary Material

Data/Code
Email
Drivers of adaptive evolution during chronic SARS-CoV-2 infections
Sheri Harari, Maayan Tahor, Natalie Rutsinsky, Suzy Meijer, Danielle Miller, Oryan Henig, Ora Halutz, Katia Levytskyi, Ronen Ben-Ami, Amos Adler, Yael Paran, Adi Stern
medRxiv 2022.02.17.22270829; doi: https://doi.org/10.1101/2022.02.17.22270829
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area