Clinical, Serological, Whole Genome Sequence Analyses to Confirm SARS-CoV-2 Reinfection in Patients From Mumbai, India

Background: SARS-CoV-2 infection may not provide long lasting post-infection immunity. While hundreds of reinfections have reported only a few have been confirmed. Whole genome sequencing (WGS) of the viral isolates from the different episodes is mandatory to establish reinfection.

Methods: Nasopharyngeal (NP), oropharyngeal (OP) and whole blood (WB) samples were collected from paired samples of four individuals who were suspected of SARS-CoV-2 reinfection based on distinct clinical episodes and RT-PCR tests. Details from their case record files and investigations were documented. RNA was extracted from the NP and OP samples and subjected to WGS, and the nucleotide and amino acid sequences were subjected to genome and protein-based functional annotation analyses. Serial serology was performed for Anti-N IgG, Anti- S1 RBD IgG, and sVNT (surrogate virus neutralizing test).

Findings: Three patients were more symptomatic with lower Ct values and longer duration of illness. Seroconversion was detected soon after the second episode in three patients. WGS generated a genome coverage ranging from 80.07 to 99.7%. Phylogenetic analysis revealed sequences belonged to G, GR and “Other” clades. A total of 42mutations were identified in all the samples, consisting of 22 non-synonymous, 17 synonymous, two in upstream, and one in downstream regions of the SARS-CoV-2 genome. Comparative genomic and protein-based annotation analyses revealed differences in the presence and absence of specific mutations in the virus sequences from the two episodes in all four paired samples.

Interpretation: Based on the criteria of genome variations identified by whole genome sequencing and supported by clinical presentation, molecular and serological tests, we were able to confirm reinfections in two patients, provide weak evidence of reinfection in the third patient and unable to rule out a prolonged infection in the fourth. This study emphasizes the importance of detailed analyses of clinical and serological information as well as the virus's genomic variations while assessing cases of SARS-CoV-2 reinfection.

Introduction

In December 2019, a novel coronavirus (n-CoV-19) sparked an outbreak in Wuhan, China. This virus was subsequently named SARS-CoV-2 and the disease COVID-19. On 11th March 2020, there were 1,18,000 cases in 114 countries with 4,291 deaths and the World Health Organization (WHO) declared that COVID-19 was a pandemic (1).

In August, the first report of reinfection by a phylogenetically distinct strain of SARS-CoV-2 was confirmed in Hong Kong (2) and subsequently Nevada reported a confirmed reinfection in USA (3). While there have been many reports of putative reinfections based on RT-PCR positivity, this has been confounded by prolonged shedding of viral RNA in the absence of replication competent virus (4) which can continue to cause RT-PCR positivity for up to at least 83 days (5). Nevertheless, the samples from the two episodes can be sequenced and genomic analysis may demonstrate genetic variation that can't be explained by short term in vivo evolution, which when combined with epidemiological and clinical evidence, may confirm reinfection (2, 3).

The present study was undertaken using samples collected from individuals tested for SARS-COV-2 as standard of care either for contact tracing or diagnostic purposes in symptomatic individuals. We report a case series of four individuals who had asymptomatic or mild RT-PCR proven COVID-19 followed by a second symptomatic RT-PCR positive episode with lower Ct values and varying degrees of increased clinical severity in the second episode.

Materials and Methods

Study Design and Participants

We identified four individuals who had tested RT-PCR positive for SARS-CoV-2 between April to June 2020 and who tested RT-PCR positive for SARS-CoV-2 once again between July to September after presenting with symptoms suggestive of COVID-19. Based on the RT-PCR results and clinical presentation of the patients, we suspected reinfection with SARS-CoV-2. Upon confirmation of the RT-PCR findings, whole genome sequencing was performed on the stored paired samples. Clinical findings and investigations were retrieved from their case records. Blood samples were collected prior to and after the second episode for anti-SARS-CoV-2 serology including anti-N, anti-S1 RBD, sVNT (surrogate virus neutralization test). The study was approved by the Institutional Review Board of Kasturba Hospital of Infectious Diseases; IRB number 015/2020. The patients provided written informed consent.

Procedures

Sample Collection

Nasopharyngeal (NP) and oropharyngeal (OP) samples for SARS-CoV-2 RT-PCR were collected, aliquoted and stored for future use as detailed in the Supplementary Table 1. Phlebotomy was performed and blood was collected in dipotassium EDTA tubes for anti-SARS-CoV-2 serology at time points between the first and second episode, early in the second episode and a longitudinal sample as described in Table 1.

TABLE 1

Table 1. Clinical course, RT-PCR, and serology.

RT-PCR

One of the aliquots was used for RNA extraction and tested by multiplex real time RT-PCR TaqPath™ COVID19 RTPCR kit for the qualitative detection of nucleic acid of SARS-CoV-2 from Applied Biosystems. Additional details of RT-PCR testing are described in Supplementary Table 1.

Serology

Anti-N protein IgG antibodies were tested by qualitative ARCHITECT chemiluminescence microparticle immunoassay (Abbott Laboratories, USA). Anti-S1 RBD antibodies were tested using SARS-CoV-2 Total antibody test on Atellica IM analyzer (Siemens, Germany). Neutralizing antibodies were tested by SARS-CoV-2 Surrogate Virus Neutralization test (GenScript USA, Inc).

Whole Genome Sequencing

Extracted RNA from all four paired stored samples was transported at −80°C for whole genome sequencing. Sample preparation, sequencing, and data analysis was performed by previously published protocols (6). Briefly, double-stranded cDNA was synthesized from 50 ng of total RNA for all the SARS-CoV-2 positive samples. The first strand of cDNA was synthesized using Superscript IV followed by RNA digestion with RNase H for second strand synthesis using DNA Polymerase I Large fragment (Klenow fragment). One hundred nanograms of purified double-stranded cDNA for both pools of ARTIC tiling PCR primers (V3 Primer pools) were taken forward. Post-amplification, pool 1 and 2 amplicons were pooled and purified using 1x AMpure beads (AMPure XP, Beckman Coulter, Cat. No. A63881). Further, 200 ng of each purified sample of multiplexed PCR amplicons obtained was taken for library preparation using Oxford Nanopore Technology (ONT) as per document no. PTC_9096_V109_REVf_06fEB2020. This included End Repair/dA tailing, Native Barcode Ligation, and Adapter Ligation of the PCR amplicons. One hundred nanograms of the pooled and purified library was sequenced using ONT's MinION Mk1B platform.

Phylogenetic and Comparative genomic analysis

Samples were base called and demultiplexed using Guppy basecaller (https://community.nanoporetech.com). Reads having phead quality score <7 were discarded to filter the low-quality reads. The resulting fastq files were normalized by read length (300–500) and reads were aligned using Minimap2 (v2.17) (7) to the reference (MN908947.3). Variants were called using Nanopolish (8) from the aligned reads and further creating consensus fasta using bcftools (v1.8) (9). Assembled fasta files from the SARS CoV-2 were aligned using CLC workbench and a UPGMA tree was constructed using default parameters. A secondary tree was generated after downloading whole genome sequences from VIPR (10) database from India submitted during the period from March 2020 to June 2020. Phylogenetic Analysis was done on all the compiled datasets using Vipr.

Lineage Analysis

Further, the assembled SARS-CoV-2 genomes were assigned lineages using the package Phylogenetic Assignment of Named Global Outbreak LINeages (PANGOLIN) (11).

Protein-Based Annotation

In order to categorize the specific amino acid variants present, the genomes were annotated by SnpEff version 4.5 (12). NC_045512 was taken as the reference genome of SARS-CoV-2 (13). The synonymous variants were filtered out from the analysis. The global frequency data for these 12 unique missense variations present across the four pairs was taken from cov-GLUE database which lists amino acid changes observed in GISAID SARS-CoV-2 sequences (14, 15). Total number of GISAID sequences retrieved at the time of analysis was 82,927, out of which 75,734 passed the exclusion criteria of CoV-GLUE.

Role of the Funding Source

The funder of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Results

A timeline summary of the clinical presentation during the two episodes, RT-PCR testing and serology are provided in Figure 1.

FIGURE 1

Figure 1. Timeline of infections in the four patients along with their serological and RT PCR investigations.

Clinical Analysis Reveals Increased Severity in the Second SARS-CoV-2 Episode

The four patients included in the study were assigned the IDs of Patient A, Patient B, Patient D and Patient E and their follow up samples were suffixed with f/u after each ID. Patient A (aged 27, male), B (aged 31, male) and D (aged 24, female) had no history of pre-existing illnesses or immunodeficiency. They were all directly involved in the clinical care of COVID-19 patients. Patient E (aged 51, female) was a controlled hypertensive, had no history of other pre-existing illnesses or immunodeficiency and worked as a technician in a COVID-19 diagnostic laboratory. Patients A and B were tested as part of a contact tracing exercise. Patient A developed sore throat and rhinitis 2 days after testing positive and recovered completely in 2 days. Patient B remained asymptomatic. Counting from their first positive RT-PCR tests A and B tested RT-PCR negative 4 and 3 days later, respectively. On day 64 and 62, respectively, they both developed COVID-19 like symptoms. Patient A tested RT-PCR positive on day 65 (1 day after symptom onset) and patient B on day 64 (2 days after symptoms onset). Patient A had fever, cough, myalgia and fatigue that lasted a week while Patient B had myalgia that lasted 2 days. Patient D's first episode was symptomatic and she tested RT-PCR positive a day after symptom onset. Symptoms included sore throat, rhinitis, and myalgia and lasted 5 days. Counting from the first positive RT-PCR, Patient D developed symptoms compatible with COVID-19 on day 52 and 2 days later on Day 54 tested positive by RT-PCR. Symptoms included sore throat, rhinitis, cough, fever, myalgia, and fatigue. Most symptoms resolved in 3 weeks but fatigue persisted for over a month. Patients A, B, and D were hospitalized during both episodes for isolation and monitoring. All three had normal respiratory rates, pulse oximetry and chest X-rays during both episodes. During the first episode, Patient E developed cough, fever, myalgia, and tested RT-PCR positive 2 days after symptoms onset. Fever remitted in 5 days but fatigue persisted. Counting from the first positive test, on day 3, RT-PCR was negative. On day 136, Patient E developed symptoms compatible with COVID-19 and 3 days later on day 139, tested RT-PCR positive. Symptoms included fever, cough, breathlessness, myalgia, nausea, and abdominal pain. Fever lasted 8 days, but breathlessness on exertion and fatigue persisted for more than 6 weeks. She was hospitalized for isolation and monitoring in the first episode but was managed as an outpatient during the second episode. Her respiratory rate and pulse oximetry were normal during both episodes but a HRCT of the chest during the second episode demonstrated pneumonia and pulmonary fibrosis. In all four patients, the second episode was more symptomatic and lasted longer in duration. All four reported that their second episodes were subjectively worse.

RT-PCR samples were collected within 3 days of symptom onset for all patients during both episodes. Patient A's sample was collected 2 days before and 1 day after symptom onset in the first and second episodes, respectively. Patient B was asymptomatic during the first episode and the sample was collected 2 days after symptom onset in the second episode. Patient D's samples were collected 1 day and 2 days after symptom onset in the first and second episodes, respectively. Patient E's samples were collected 2 days and 3 days after symptom onset in the first and second episodes, respectively. Similar time points of sample collection for the first and second episodes for the patients along with harmonized RT-PCR sample collection, processing and testing methodology allowed us to compare Ct values despite the short window for RT-PCR positivity in some COVID-19 patients. Patients A, D, and E had lower Ct values in the second episode compared to the first. Patient B's Ct values were higher during the second episode. Details of Ct values are presented in Table 1.

Seroconversion Detected After the Second Episode

Three serological tests performed, anti-N IgG, anti-S1 RBD IgG, and neutralizing antibodies by sVNT. Counting from the first positive RT-PCR test, on day 47 Patients A and B were both negative for anti-N IgG antibodies. Their plasma samples drawn on day 47 were not stored for additional tests (which became available later). On day 69 both patients had already developed symptoms for the second time and serological sampling was repeated. Patient A became symptomatic 5 days prior and RT-PCR positive 4 days prior to serological sampling. Patient A's sample was sVNT was positive but anti-N and anti-S1 RBD IgG were both negative. Patient B became symptomatic 7 days prior and RT-PCR positive 5 days prior to serological sampling. All three serological tests were negative on day 69. A third sample was drawn for both A and B on day 124. All three serological tests were positive for Patient A. Patient B was positive by sVNT but negative for anti-N and anti-S1 RBD IgG. Counting from the first positive RT-PCR, on day 21 Patient D was negative for all three antibodies. On day 55, just 3 days after symptom onset and 1 day after RT-PCR positivity in the second episode, Patient D was positive for all three serological tests. A longitudinal sample collected on day 73 was more strongly positive for all three tests. Counting from the first positive RT-PCR test, Patient E tested negative for all three antibodies on day 137 (1 day after symptom onset in the second episode). On day 153 (17 days after symptom onset in the second episode) Patient E was positive for all three antibodies.

Genome Analysis Reveals Clade Change and/or Distinct Mutations in the Virus Populations Between Episodes

Genome sequencing generated genome coverage of 80.07–99.7% (Table 2).The assembled genomes were curated and taken for further analysis. Phylogenetic tree analysis of the eight sequences, along with 160 complete viral genome sequences submitted from India in GISAID between the months of May to September 2020 because both phases of samples used for the study has been collected in this duration, revealed two samples (Patients A and B) sub-clustered together with their f/u samples respectively while samples Patient D and E and their f/u sequences clustered in different clades (Figure 2).

TABLE 2

Table 2. Clade, lineage of patients with reinfections (n = 4).

FIGURE 2

Figure 2. Circular Phylogram generated using UPMGA on MEGAX. A total of 160 sequences were used in the analysis. Each patient sample pairs are colored. Patient A and f/u jade green, Patient B and f/u olive green, Patient D and f/u orange, and Patient E and f/u red. Sequences downloaded from the public database are colored in purple.

Clade based analysis revealed that two of eight sequences belonged to the G clade while one sequence belonged to clade GR while the remaining five sequences categorized under “Other” category. Further, analysis of lineage by PANGOLIN revealed distribution of the eight with variations of B lineages including B, B.1, B.1.80, and B.1.1.32 (Table 2).

The samples from the first and second episode of infection of the four patients are predominantly from the SARS-CoV-2 clade 19A and 20A. The clades from the first and second episode, respectively, were 20A and 19A in Patient A, 20B and 20B in patient B, 19A and 20B in Patient D, and 19A and 20B in Patient D.

Mutation analysis of the samples revealed distinct mutations in all the samples (Table 2). Interestingly, we observed a higher number of mutations in the follow-up samples except Pair-B, which had 10 mutations in first infection compared to three in the follow-up. Pair-E had the highest number of 13 mutations in the follow-up sample compared to two in the first sample, followed by Pair-D with 10 mutations in follow-up and one in the first sample and lastly by Pair-A with two in follow-up and one in the first sample. A total of 42 (Figure 3) mutations were observed in our sample set of four patients. Twenty-two non-synonymous, 17 synonymous, and 2 upstream UTR and 1 downstream UTR mutation is observed. Interestingly the non-synonymous mutation P323L in the nsp 12 RNA-dependent RNA polymerase gene has been reported to be concurrently present with D614G mutation in the spike protein, is observed in all patient samples, whereas D614G mutation was observed only in four samples (16, 17). In the nsp3 region, part of the replicase complex, two synonymous mutations F924F, N1123N, and one non-synonymous mutation A1812D observed in mild cases of COVID-19 (18) were observed in Patient E, Patient B, and Patient B f/u samples, respectively.

FIGURE 3

Figure 3. Heatmap with 42 overall mutations (unique set of 25 mutations). Purple, mustard yellow, and sky blue colors show the presence of mutation in samples while gray color shows the absence of mutation in samples.

To evaluate amino-acid alterations, we performed protein-based annotation of the 22 non-synonymous mutations found from our genome analysis of the four pair of samples (Figure 4). It was observed that Pair 1, i.e., Patient A shows minor variations, with common ones occurring within Nsp12. With respect to the other patients, interestingly, we found heterogeneity within mutations in both episodes. For instance, in Patient B, the mutations within Spike protein (D614G, Q677H) in the first episode were missing in the followup sample. Similarly, in Patients D and E, we found presence of additional mutations in samples of followup. Interestingly, in re-infection cases, a higher number of mutations were found in non-structural proteins, including nsp1, nsp2, nsp3, nsp5, nsp6, and nsp12, and nsp 14. Further, we also performed correlations of these mutations with viral genomes from world-wide populations (~1,44,426) to understand their relative frequency (Figure 2). While P323L mutation within nsp12 was found in all samples without exception, other frequent mutations showed abrupt patterns. In particular, D614G mutation within the Spike protein was consistently present in both infections in Patient E but was present only in one of the episodes in Patients B and D.

FIGURE 4

Figure 4. Mapping of amino-acid substitutions within SARS-CoV-2 genome of four pairs of samples. The upper plot demonstrates the seven proteins in different colors that harbor 12 non-synonymous mutations shown in dots. The Y-axis shows the four pair of patient samples. The blue and red dot indicates the presence of the mutation in the first and second episode of infections respectively. The lower plot shows the frequency of that particular mutation in 82,927 genomes deposited in GISAID.

Sequence Submission

All SARS-CoV-2 sequences from eight patients were submitted to GISAID under the accession number EPI_ISL_528419 and EPI_ISL_528420 for patient A, A_f/u, EPI_ISL_528421, and EPI_ISL_528422 for patient B, B_f/u, EPI_ISL_528425, and EPI_ISL_528426 for patient D, D_f/u, EPI_ISL_801538, and EPI_ISL_676509 for patient E, E_f/u.

Discussion

Clinically SARS-CoV-2 infection can present with or without symptoms and severity has been categorized into four types ranging from asymptomatic to critical illness based on symptoms, clinical findings, chest imaging and blood gases as presented in Supplementary Figure 1 (19). New immunological evidence is enriching our knowledge of the immune response to infection (20) and duration of immunity following infection (21). Emerging evidence suggests Ct values and viral loads at the time of diagnosis maybe implicated in pathogenesis and disease severity (22). A handful of confirmed SARS-CoV-2 reinfection have been published on the basis of genome variation observed in the viruses between the two episodes with varying clinical manifestations between the episodes (2, 3, 23, 24). The European Center for Disease Control and Prevention (ECDC) (25) and United States Center for Disease Control and Prevention (US CDC) (26) have considered multiple criteria to investigate a case of suspected reinfection.

On the basis of these criteria, we discuss our patients and confirm or reject a case as SARS-CoV-2 reinfection. As per the US CDC, SARS-CoV-2 reinfection should be considered in individuals with COVID-19 like symptoms and a positive RT-PCR for SARS-CoV-2 with a Ct value <33 at least 45 days after the first positive RT-PCR. There should not be an obvious alternative etiology for the symptomatic second episode. Paired samples from the two episodes should undergo genomic testing that includes evaluation of single nucleotide variations (SNV) and clades to distinguish between viral persistence within host evolution vs. reinfections. In patients meeting the above criteria, genomic testing revealing differing clades as defined in Nextstrain (27) and GISAID of SARS-CoV-2 between the first and second infection is considered the best evidence of SARS-CoV-2 reinfection. More than two nucleotide differences per month in consensus between sequences that meet quality metrics is considered moderate evidence. The US CDC also recommends serial serological testing.

Accordingly our present study evaluates clinical, RT-PCR, genomic and serological information to evaluate reinfections in four patients who presented with repeat episodes of SARS-CoV-2 infections. Of the four patients in the study, Patients A, D, and E had COVID-19 like symptoms during both first episodes and second episode and did not have an obvious alternate etiology for their COVID-19 like symptoms. Their symptoms were also accompanied by a positive RT-PCR for COVID-19 over 45 days from the first positive RT-PCR. Interestingly, Patients A, D, and E had increased clinical severity and lower Ct values in the second episode. All three had Ct values not exceeding 23. Such Ct values correlate with active viral replication and positively correlate with virus culture positivity (28). Analysis of whole genome sequence data generated from the samples of both episodes of Patients A, D and E revealed that the two paired samples clustered in different clades and belonged to different lineages.

Patient A's paired samples contained viruses from different clades but were separated by a single mutation. Moreover, the sample from the second episode had low Ct values (23 in confirmatory gene) and the clinical picture strongly suggested active SARS-CoV-2 infection. Crucially, Patient A was positive for neutralizing antibodies just 5 and 4 days after symptom onset and RT-PCR positivity during the second episode. While WGS showed a single distinct mutation in consensus sequences, the clinical picture, low Ct values, difference in clade and presence of neutralizing antibodies within 5 days of symptom onset supports reinfection. It should be noted that Patient A's first sample genome coverage was 80.37 and in the second episode was 83.01. This could have resulted in detection of fewer mutations. Despite the clade change, clinical picture, lower Ct values, and nAb positivity, with the caveat of genomic coverage and based on the CDC criteria for defining reinfection, we determined the evidence as weak evidence for assigning the second episode of Patient A as a reinfection.

Patient B was asymptomatic in the first episode and but had a symptomatic second episode about 2 months later with myalgia and malaise. The Ct value from samples for RT-PCR was 33 in the first episode but 36 in the second episode. The genome analysis of the paired samples of this patient further showed no clade or lineage difference. However, mutation analysis revealed difference in mutations observed including the presence of the D614G mutation only in the sample from the first episode. There were addition/deletion of both synonymous and non-synonymous mutations between the samples of the two episodes as was observed in the functional protein annotation analysis. Most of the mutations were found in the spike protein, the region most likely to undergo mutations to escape immune pressure during prolonged infections. Three synonymous and two non-synonymous mutations occurred in the spike region. Additionally, in the second episode, 7 and 5 days after symptom onset and RT-PCR positivity all three antibody tests (anti-N, anti-S1 RBD, and sVNT) were negative. All these analyses put together make it difficult to differentiate between a prolonged infection and a reinfection in Patient B.

Both patient D and E had symptoms compatible with COVID-19 during both episodes and the clinical picture was strongly suggestive of COVID-19. Both had lower Ct values in the second episode suggestive of active viral replication. Additionally, during the second episode Patient E had radiological evidence of acute pulmonary infection (pneumonitis) superimposed on COVID-19 pulmonary sequelae (pulmonary fibrosis). Paired samples from both Patient D and E contained viruses from different clades and had distinct mutations exceeding the cut off requiring >2 distinct mutations per month between consensus sequences clearly confirming SARS-CoV-2 reinfection.

In the present study, we found priming of immunity in the first episode leading to a boosting effect following the second episode by production of neutralizing antibodies early in the second episode. Analysis of the serological profiles of all the patients failed to reveal seroconversion after the first episode but during the second episode, neutralizing antibodies were detected 5 and 3 days after symptom onset as seen in Patients A and D, respectively. Further, longitudinal samples of these patients revealed increasing titers of neutralizing antibodies. In the case of Patient E, seroconversion was not detected early in the second episode but was observed two and a half weeks after symptoms onset. While most individuals do seroconvert following SARS-CoV-2 infection, some individuals do not seroconvert (20). It is possible that the patient sin our study had failure of humoral immunity which may explain the absence of detectable antibodies. It is possible that the absence of seroconversion predisposed them to reinfection.

While our study found that the second episode was more symptomatic with a longer duration of illness, our study was not designed to identify reasons for increased severity in the second episode. Nevertheless, we hypothesize a few possible reasons for the observed increased severity in the second episode.

Some evidence from animal studies suggests that increased inoculum size or a higher infecting dose may result in increased clinical severity (29). Owing to their status as health care workers caring for COVID-19 patients or handling their samples all four patients had an occupational risk of exposure. It is possible that the participants in our study were exposed to a larger infecting dose in the second episode as compared the primary infection. Another aspect to consider is the impact of mutations in the viral genome. Recent detection of SARS-CoV-2 variants has raised important questions about the impact of S gene mutations and deletions on increased transmissibility, ACE-2 receptor affinity, viral loads, immune escape, and severity. S variants of SARS-CoV-2 have been associated with significantly lower median Ct values suggesting that changes in the S protein RBD may result in increased viral loads (30). While our sample size and absence of viral culture studies does not allow us to make determinations about the impact of S gene mutations and deletions on clinical severity and viral load, it is possible that mutations at the Spike gene may explain lower Ct values and increased severity in the second episodes.

Some experimental in vitro studies suggest the possibility of antibody dependent enhancement of SARS-CoV-2 (31, 32) which has also been observed in other coronaviruses. It is possible immune enhancement may have increased the severity of the second episode.

Taken altogether, our present study provides a level of evidence classified by US CDC as best evidence of reinfection in two patients (Patients D and E), weak evidence with possible reinfection in one patient (Patient A), and we were unable to differentiate between prolonged infection and reinfection in the case of Patient B. Our study adds to the growing body of evidence of SARS-CoV-2 reinfections and demonstrates the value of serial serological data in supporting reinfection claims. Our study highlights that SARS-CoV-2 reinfections do occur, and individuals who have recovered from SARS-CoV-2 infection should continue to take infection prevention precautions.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Kasturba Hospital of Infectious Diseases; IRB number 015/2020. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JS conceptualized and designed the study. JS and LP identified the study participants. SP and SA collected and compiled data from different sources. NC and MP performed RNA extraction, aliquoting, and RT-PCR. RP, VA, JSV, AK, RM, and SF performed genome sequencing. RP, VA, JSV, AK, RM, SF, LT, SS, SC, and CS performed genomic and linage analyses. SP, SS, and JS drafted and revised the manuscript. JS, SS, and AA provided resources and participated in overall supervision. All authors contributed to data interpretation, critically reviewed the manuscript, provided contributions to tables, figures and text in the manuscript, and approved the final manuscript for submission.

Funding

This work was supported by the Municipal Corporation of Greater Mumbai's funding for Kasturba Hospital's COVID-19 diagnostic laboratory; International Centre for Genetic Engineering and Biotechnology core funding and Council of Scientific and Industrial Research (CSIR) Institute of Genomics and Integrative Biology from CSIR project (MLP-2005) and Foundation Botnar project (CLP-0031).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dean, BYL Nair Hospital and Municipal Corporation of Greater Mumbai (MCGM) for permission to publish. We thank Dr. Unnati Desai, Associate Professor, Pulmonary Medicine, BYL Nair Hospital for facilitating access to clinical information. The authors acknowledge Dr. Camilla Rodrigues for sample transfer and Dr. Gagandeep Kang through THSTI for facilitating sequencing and for their involvement in sample transfer for sequencing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.631769/full#supplementary-material

References

1. WHO. Director-General's Opening Remarks at the Media Briefing on COVID-19. (2020). Available online at: https://www.who.int/director-general/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19—11-march-2020 (accessed September 4, 2020).

2. To KK, Hung IF, Ip JD, Chu AW, Chan WM, Tam AR, et al. COVID-19 re-infection by a phylogenetically distinct SARS-coronavirus-2 strain confirmed by whole genome sequencing. Clin Infect Dis. (2020). doi: 10.1093/cid/ciaa1275. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Tillett R, Sevinsky J, Hartley P, Kerwin H, Crawford N, Gorzalski A, et al. Genomic evidence for reinfection with SARS-CoV-2: a case study. Lancet Infect Dis. (2020) 21:52–8. doi: 10.1016/S1473-3099(20)30764-7

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wölfel R, Corman VM, Guggemos W, Seilmaier M, Zange S, Müller MA, et al. Virological assessment of hospitalized patients with COVID-2019. Nature. (2020) 581:465–9. doi: 10.1038/s41586-020-2196-x

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Li N, Wang X, Lv T. Prolonged SARS-CoV-2 RNA shedding: not a rare phenomenon. J Med Virol. (2020) 92:2286–7. doi: 10.1002/jmv.25952

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Kumar P, Pandey R, Sharma P, Dhar M, Vivekanand A, Uppili B, et al. Integrated genomic view of SARS-CoV-2 in India. Wellcome Open Res. (2020) 5:184. doi: 10.12688/wellcomeopenres.16119.1

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. (2018) 34:3094–100. doi: 10.1093/bioinformatics/bty191

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Loman N, Quick J, Simpson J. A complete bacterial genome assembled de novo using only nanopore sequencing data. Nat Methods. (2015) 12:733–5. doi: 10.1038/nmeth.3444

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. (2011) 27:2987–93. doi: 10.1093/bioinformatics/btr509

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pickett BE, Sadat EL, Zhang Y, Noronha JM, Squires RB, Hunt V, et al. ViPR: an open bioinformatics database and analysis resource for virology research. Nucleic Acids Res. (2012) 40:D593–8. doi: 10.1093/nar/gkr859

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Rambaut A, Holmes EC, O'Toole Á, Hill V, McCrone JT, Ruis C, et al. A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat Microbiol. (2020) 5:1403–7. doi: 10.1038/s41564-020-0770-5

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. (2012) 6:80–92. doi: 10.4161/fly.19695

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang C, Liu Z, Chen Z, Huang X, Xu M, He T, et al. The establishment of reference sequence for SARS-CoV-2 and variation analysis. J Med Virol. (2020) 92:667–74. doi: 10.1002/jmv.25762

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Singer J, Gifford R, Cotten M, Robertson D. CoV-GLUE: a web application for tracking SARS-CoV-2 genomic variation. Preprints. (2020) 2020:2020060225. doi: 10.20944/preprints202006.0225.v1

CrossRef Full Text | Google Scholar

15. Shu Y, McCauley J. GISA ID: global initiative on sharing all in?uenza data – from vision to reality. Euro Surveill. (2017) 22:30494. doi: 10.2807/1560-7917.ES.2017.22.13.30494

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Ahn DG, Choi JK, Taylor DR, Oh JW. Biochemical characterization of a recombinant SARS coronavirus nsp12 RNA-dependent RNA polymerase capable of copying viral RNA templates. Arch Virol. (2012) 157:2095–104. doi: 10.1007/s00705-012-1404-x

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Toyoshima Y, Nemoto K, Matsumoto S, Nakamura Y, Kiyotani K. SARS-CoV-2 genomic variations associated with mortality rate of COVID-19. J Hum Genet. (2020) 65:1075–82. doi: 10.1038/s10038-020-0808-9

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Biswas SK, Mudi SR. RNA-dependent RNA polymerase and spike protein mutant variants of SARS-CoV-2 predominate in severely affected COVID-19 patients. Preprints. (2020) 2020:2020070251. doi: 10.20944/preprints202007.0251.v1

CrossRef Full Text | Google Scholar

19. Directorate General of Health Services. MoHFW, Government of India Clinical Management Protocol: COVID-19 Version 4.0. (2020). Available online at: https://www.mohfw.gov.in/pdf/ClinicalManagementProtocolforCOVID19dated27062020.pdf (accessed October 4, 2020).

20. Wajnberg A, Amanat F, Firpo A, Altman DR, Bailey MJ, Mansour M, et al. Robust neutralizing antibodies to SARS-CoV-2 infection persist for months. Science. (2020) 370:1227–30. doi: 10.1126/science.abd7728

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Seow J, Graham C, Merrick B, Acors S, Pickering S, Steel KJA, et al. Longitudinal observation and decline of neutralizing antibody responses in the three months following SARS-CoV-2 infection in humans. Nat Microbiol. (2020) 5:1598–607. doi: 10.1038/s41564-020-00813-8

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Rao SN, Manissero D, Steele VR, Pareja J. A systematic review of the clinical utility of cycle threshold values in the context of COVID-19. Infect Dis Ther. (2020) 9:587. doi: 10.1007/s40121-020-00328-z

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Van Elslande J, Vermeersch P, Vandervoort K, Wawina-Bokalanga T, Vanmechelen B, Wollants E, et al. Symptomatic SARS-CoV-2 reinfection by a phylogenetically distinct strain. Clin Infect Dis. (2020). doi: 10.1093/cid/ciaa1330. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Torres DA, Ribeiro LDCB, Riello APFL, Horovitz DDG, Pinto LFR, Croda J. Reinfection of COVID-19 after 3 months with a distinct and more aggressive clinical presentation: case report. J Med Virol. (2020) 93:1857–9. doi: 10.1002/jmv.2663723

PubMed Abstract | CrossRef Full Text | Google Scholar

25. European CDC. Reinfection with SARS-CoV-2: Considerations for Public Health Response. (2020). Available online at: https://www.ecdc.europa.eu/sites/default/files/documents/Re-infection-and-viral-shedding-threat-assessment-brief.pdf (accessed October 4, 2020).

26. US CDC. Common Investigation Protocol for Investigating Suspected SARS-CoV-2 Reinfection. (2020). Available online at: https://www.cdc.gov/coronavirus/2019-ncov/php/reinfection.html (accessed October 31, 2020).

27. Hadfield J, Megill C, Bell SM, Huddleston J, Potter B, Callender C, et al. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics. (2018) 34:4121–3. doi: 10.1093/bioinformatics/bty407

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Basile K, McPhie K, Carter I, Alderson S, Rahman H, Donovan L, et al. Cell-based culture of SARS-CoV-2 informs infectivity and safe de-isolation assessments during COVID-19. Clin Infect Dis. (2020). doi: 10.1093/cid/ciaa1579. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Ryan KA, Bewley KR, Fotheringham SA, Slack GS, Brown P, Hall Y, et al. Dose-dependent response to infection with SARS-CoV-2 in the ferret model and evidence of protective immunity. Nat Commun. (2021) 12:81. doi: 10.1038/s41467-020-20439-y

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kidd M, Richter A, Best A, Mirza J, Percival B, Mayhew M, et al. S-variant SARS-CoV-2 lineage B1.1.7 is associated with significantly higher viral loads in samples tested by ThermoFisher TaqPath RT-qPCR. J Infect Dis. (2020). doi: 10.1101/2020.12.24.20248834. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Wu F, Yan R, Liu M, Liu Z, Wang Y, Luan D, et al. Antibody-dependent enhancement (ADE) of SARS-CoV-2 infection in recovered COVID-19 patients: studies based on cellular and structural biology analysis medRxiv. (2020). doi: 10.1101/2020.10.08.20209114

CrossRef Full Text | Google Scholar

32. Yip MS, Leung NH, Cheung CY, Li PH, Lee HH, Daëron M, et al. Antibody-dependent infection of human macrophages by severe acute respiratory syndrome coronavirus. Virol J. (2014) 11:82. doi: 10.1186/1743-422X-11-82

PubMed Abstract | CrossRef Full Text