Introduction

SARS-CoV-2 is described as a causative agent of a pneumonia disease that was firstly reported in late 2019 in Wuhan, Hubei province, central China [1]. The fast spread of the virus worldwide led World Health Organization (WHO) to characterize Coronavirus Disease (COVID-19) as a pandemic on March 11, 2020. Since then, there have been more than 210 million confirmed cases and more than 4 million deaths caused by COVID-19 globally (https://covid19.who.int/).

A suspected case of reinfection may be described when a person has symptoms compatible with COVID-19 and a positive test result for SARS-CoV-2 following a period ≥ 45 days from the first infection for which prolonged shedding of viral RNA and infection by another agent has been ruled out [2]. Furthermore, to confirm a reinfection case, SARS-CoV-2 complete genomic sequencing for both primary and secondary infections indicating that they belong to different genetic clades or lineages should be used to dismissing recurrence or prolonged viral excretion[2, 3].

Reinfection cases might happen in convalescent or vaccinated individuals due to the lower neutralizing capacity of SARS-CoV-2 antibodies and partial immunity against new variants[4]. However, cases of reinfection have been reported for both VOC[5,6,7] and non-VOC strains[8,9,10]. In the State of Rio Grande do Sul (RS), South of Brazil, the lineages that were circulating during May to November of 2020 were, most of them, classified as B.1.1, B.1.1.33, B.1.1.248, and B.1.1.49[11, 12]. In late November, it was possible to see an increase of P.2 (Zeta) in the region[12]. One case of P.1 (Gamma) in late November was described in RS; however, only in late January of 2021 was evidenced the circulation of Gamma variant, coinciding with the peak of SARS-CoV-2 new cases, hospitalizations, and deaths[5].

Although many recent articles have highlighted the importance of possible Spike epitopes between former lineages and new variants as a potential cause of an increase in the frequency of reinfections[7, 13, 14], in this study, we reported two SARS-CoV-2 reinfection cases caused by non-variant-related lineages in South Brazil. Both patients had a first RT-qPCR positive result in July of 2020 and a second in September of 2020, with 45 and 61 days between both infections for patients 1 and 2, respectively. Samples were sequenced and the lineage responsible for infections were B.1.91, B.1.1.33, and B.1.1.28.

Material and methods

SARS-CoV-2 detection by RT-qPCR

For SARS-CoV-2 detection, naso-oropharyngeal swab samples were received at Laboratório de Microbiologia Molecular of Universidade Feevale. Viral RNA was extracted with the commercial MagMAX™ CORE Nucleic Acid Purification Kit (Applied biosystems™, Thermo Fisher Scientific, Waltham, MA, USA) kit, using the automated equipment KingFisher™ Duo Prime (Thermo Fisher Scientific™). Real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) for SARS-CoV-2 selective for the envelope (E) gene was performed according to the Charite Institute, Berlin, Germany, protocols[15] and using AgPath-ID One-Step RT-PCR Reagents (Thermo Fisher Scientific™). Assays were performed adding positive and negative controls. Both template control (NTC) and negative controls diluted in Milli-Q ultrapure RNase/DNase-free water were included to confirm lack of contamination among samples. The genome of the positive control used in all assays was deposited in the GISAID genomic database under access number EPI_ISL_2229839.

SARS-CoV-2 genome sequencing and phylogenetic analysis

Viral genome sequencing and phylogenetic analysis were performed as previously described by da Silva et al. (2021)[5]. Briefly, after naso-oropharyngeal swab samples extraction, as described above, reverse transcription reaction was carried out by SuperScript IV reverse transcriptase kit (Thermo Fisher Scientific, Waltham, MA, USA). Viral genome library preparation was carried out using the QIAseq SARS-CoV-2 Primer Panel paired for library enrichment and QIAseq FX DNA Library UDI Kit, according to the manufacturer’s instructions (QIAGEN, Hilden, Germany). Illumina MiSeq platform (Foster City, CA, USA), using MiSeq Reagent Kit v3 (600-cycle), was used. After sequencing, Geneious Prime™ suite was used for genome assembly and editing. Briefly, FASTQ reads were imported, trimmed to remove low-quality sequences and primers used for library generation (BBDuk 37.25), and mapped against the reference sequence hCoV-19/Wuhan/WIV04/2019 (EPI_ISL_402124) available in the EpiCoV database from GISAID[16]. Sequences generated in this study were aligned with 97 Brazilian SARS-CoV-2 complete genomes and the reference sequence (> 29 kb) that were retrieved from the GISAID database. Alignment using the reference sequence from Wuhan as the outgroup was performed using the Clustal Omega software. The Maximum likelihood phylogenetic analysis under the general time reversible model, allowing for a proportion of invariable sites and substitution rates, was inferred empirically in IQ-TREE v2.1.2 web server[17], especially for maximum-likelihood (ML, applying 200 replicates and 1000 bootstraps.

Results

Patient 1 is a 47-year-old woman living in Igrejinha city in the southernmost state of Brazil. In the 1st infection episode (hCoV-19/Brazil/LMM16242/2020, GISAID access: EPI_ISL_2928341), she had symptoms that started on July 15, 2020, and naso-oropharyngeal swab samples collection for COVID-19 diagnostic by RT-qPCR was on July 20, 2020. In the reinfection (hCoV-19/Brazil/LMM25714/2020, GISAID access: EPI_ISL_2928147), the first symptoms appeared on August 31, 2020, and the sample was collected on September 3, 2020. The interval between the two positive RT-qPCR for SARS-CoV-2 was 45 days, and results were 20.15 and 16.77 CT (cycle threshold) values for E gene, for infection and reinfection, respectively. She was symptomatic in both events, being that, in the 1st episode, we did not have access to the full record of symptoms, and in the reinfection, she reported myalgia, dry cough, loss of appetite, and dyspnea.

Patient 2 is a 28-year-old woman living in São Leopoldo and working in Novo Hamburgo as a nursing technique, RS, Brazil. In the 1st infection (hCoV-19/Brazil/LMM14996/2020, GISAID access: EPI_ISL_2928141), she had symptoms that started on July 10, 2020, and naso-oropharyngeal swab sample collection for COVID-19 diagnostic was on July 13, 2020. In reinfection (hCoV-19/Brazil/LMM27281/2020, GISAID access: EPI_ISL_2928344), the first symptoms appeared on September 4, 2020, and the sample was collected on September 12, 2020. The result of RT-qPCR was 14.91 and 15.43 CT values for E gene, for infection and reinfection, respectively. The interval between the two positive results for SARS-CoV-2 was 61 days, and she was symptomatic in both events. In the 1st episode of COVID-19, she reported experiencing flu symptoms, and, in the reinfection, she described body pain, tiredness, sore throat, headache, and runny nose. The patient’s data were described in Table 1.

Table 1 Timeline of SARS-CoV-2 reinfection cases in two patients from Southern Brazil

Whole-genome sequences for SARS-CoV-2 were covered in high quality for analyzed samples being two distinct lineages present in samples from the same person in different infections episodes demonstrating reinfection for both patients. According to Pangolin COVID-19 Lineage Assigner tool[18], sequences of patient 1 were classified as B.1.91 (hCoV-19/Brazil/LMM16242/2020) in infection and B.1.1.33 (hCoV-19/Brazil/LMM25714/2020) in reinfection (GISAID access: EPI_ISL_2928341 and EPI_ISL_2928147, respectively). Sequences of patient 2 were classified as B.1.1.33 (hCoV-19/Brazil/LMM14996/2020) in infection and B.1.1.28 (hCoV-19/Brazil/LMM27281/2020) in reinfection (GISAID access: EPI_ISL_2928141 and EPI_ISL_2928344, respectively). Our phylogenetic analysis also confirmed these results clustering the four samples with their respective branches of B.1.91, B.1.1.28, and B.1.1.33 (Fig. 1).

Fig. 1
figure 1

Phylogenetic tree of SARS-CoV-2 complete genomes with the two cases of reinfection

Discussion

This study reported two reinfection cases of SARS-CoV-2 by non-VOC strains in the South of Brazil. Although reinfection cases of SARS-CoV-2 by VOC strains[5,6,7] have been linked to lower neutralizing capacity antibodies and partial immunity against new variants[4], a limited and transitory immunity induced by primary infection might be responsible for the occurrence of secondary one reflecting the virus ability to evade the previous immune responses[7] also in reinfection cases caused by non-VOC strains, as were described in others studies[8,9,10].

In our study, both patients had mild to moderate symptomatic cases in primary and secondary infection. Even though there are cases of reinfection with an increase in severity[5, 6], more cases have been showing less severity, suggesting partial immune protection from disease[19]. The four samples presented CT values for the E gene between 14 and 20. As already described, when CT values are higher than 24, the SARS-CoV-2 infectivity might be significantly reduced[20]. Therefore, reinfection cases described in our study contributed to onward transmission of the virus if isolation measures were not correctly applied.

There is an additional concern when a case of reinfection is linked to a health professional, and it is the history of patient 2. This is due to the fact that transmission within the hospital environment can affect, in addition to professionals, patients who are already debilitated. Therefore, it is important to efficiently detect and manage healthcare-associated infections preventing transmission and complications[21].

Genomic surveillance of SARS-CoV-2 is responsible for reporting crucial issues to support epidemiological investigations such as the circulating lineages, the introduction of new sources of virus from other geographical areas, and the evolution of the virus in response to applied control measures[21]. The three lineages sequenced (B.1.91, B.1.1.28, and B.1.1.33) in our samples were descendants of B.1 lineage. Besides reinfection cases described above, SARS-CoV-2 lineage B has been related with most Brazilian sequences[22], and also with coinfection cases that might be another relevant factor to the complex interaction between immune response systems and virus mutations[12].

Moreover, B.1.1.28 lineage described herein is related as an ancestral of the Variants of Interest (VOI) P.2 (Zeta) (alias for B.1.1.28.2)[23] and P.1 (Gamma) (alias for B.1.1.28.1)[24]. Another relevant issue is that B.1.1.28 has also been related to the emergence of putative novel SARS-CoV-2 variants circulating in Southern Brazil[12, 25]. Since this lineage originated VOIs, VOCs, and proposed new variants, it is a need to continually evaluate and track its circulation in the population. Besides genomic surveillance of SARS-CoV-2, there is an evident need for clinical and epidemiological studies to enlighten how previous exposure through natural infection and vaccination protects against infection and reinfection with existing and newly emerging variants[4].

In accordance with our findings, we highlight the importance of continuous genomic monitoring of SARS-CoV-2 in samples from patients suspected of reinfection, especially because reinfection cases are not just linked to VOC that might have a significant reduction in virus neutralization by antibodies generated during infection or vaccination[4, 26]. Therefore, the correct surveillance of these cases can help to elucidate the durability of the immune response against existing and new variants and even the need for vaccine boosters in the population. Also, it may assist in decisions planning the next measures that must be applied to respond to the pandemic, evaluating the specifics and necessities of each country or locality.