Impact Factor 5.091 | CiteScore 4.1
More on impact ›

ORIGINAL RESEARCH article

Front. Med., 26 November 2021 | https://doi.org/10.3389/fmed.2021.760170

Atypical Prolonged Viral Shedding With Intra-Host SARS-CoV-2 Evolution in a Mildly Affected Symptomatic Patient

  • 1Scientific Platform Pasteur—USP, São Paulo, Brazil
  • 2Institute of Biomedical Sciences, University of São Paulo, São Paulo, Brazil
  • 3Institut Pasteur, Paris, France

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is caused by a respiratory virus with a wide range of manifestations, varying from asymptomatic to fatal cases, with a generally short outcome. However, some individuals present long-term viral shedding. We monitored 38 individuals who were mildly affected by the SARS-CoV-2 infection. Out of the total studied population, three (7.9%) showed atypical events regarding the duration of positivity for viral RNA detection. In one of these atypical cases, a previously HIV-positive male patient presented a SARS-CoV-2 RNA shedding and subgenomic RNA (sgRNA) detected from the upper respiratory tract, respectively, for 232 and 224 days after the onset of the symptoms. The SARS-CoV-2 B.1.1.28 lineage, one of the most prevalent in Brazil in 2020, was identified in this patient in three serial samples. Interestingly, the genomic analyses performed throughout the infectious process showed an increase in the genetic diversity of the B.1.1.28 lineage within the host itself, with viral clearance occurring naturally, without any intervention measures to control the infection. Contrasting widely spread current knowledge, our results indicate that potentially infectious SARS-CoV-2 virus might be shed by much longer periods by some infected patients. This data call attention to better adapted non-pharmacological measures and clinical discharge of patients aiming at preventing the spread of SARS-CoV-2 to the population.

Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in Wuhan, Hubei province, China, in late December 2019 (13) causing highly transmissible respiratory infection and acute disease in humans. The SARS-CoV-2 (family Coronaviridae, genus Betacoronavirus) quickly spread over a large geographical area (4), impacting a global scale in terms of morbidity, mortality, and economic impact. The first diagnosed SARS-CoV-2 patient in Brazil occurred on February 26, 2020 (5). The virus rapidly started its community circulation process caused by multiple lineages in a short period in all Brazilian states (6), and distinct viral lineages were identified (68).

Since the SARS-CoV-2 emergence, hundreds of thousands of viral consensus genomes have been sequenced and quickly made available around the world (9). The rapid spread of SARS-CoV-2, initially with an apparent low diversity, suggested that the ancestor of the virus accumulated a genetic diversity allowing for the emergence of several phylogenetic lineages while spreading geographically across the world (9). Recently, genetically diverse viral populations that co-circulate within a single host began to be explored in some population groups infected with SARS-CoV-2 (1013).

Generally, SARS-CoV-2 causes acute respiratory syndrome with a rapid and broad clinical outcome. However, some studies had reported viral RNA (vRNA) and infectious particles for an atypical considered period (1416). Other studies that characterized the intra-host viral diversity in patients with persistent infection showed a strong selection bias, with the introduction of neutralizing antibody cocktails or several interventions to treat secondary symptoms to the virus action (11, 12). Here, we explore the clinical and molecular aspects of the prolonged viral shedding process in a patient with evidence of viral replication for extended period. An intra-host genomic evolution process was evidenced by the detection of subgenomic RNA (sgRNA) and by the accumulated viral mutations over the period.

Materials and Methods

Ethical Statement

The samples analyzed in this study were collected after the consent of the patient following the protocol approved by the Human Research Ethics Committee of the Biomedical Sciences Institute—University of São Paulo (CEPSH-ICB-USP) (no. 4.036.252).

Patients and Samples

This study is a molecular investigation of symptomatic, mildly affected patients, sampled between April 2020 and November 2020 in the city of São Paulo. In total, 721 patients with symptoms were screened for the detection of viral RNA of SARS-CoV-2. This first criterion of inclusion consisted in the presence of cold/flu symptoms, such as fever and/or respiratory symptoms (cough, breath shortness, and sore throat), following criteria recommended by the Brazilian Ministry of Health guidelines for Coronavirus Disease 2019 (COVID-19) diagnosis and treatment. Exclusion criteria included patients presenting symptoms of viral infections related to other agents than SARS-CoV-2, according to the same guidelines. Thirty-eight mild symptomatic positive to the SARS-CoV-2 detection out of 721 symptomatic patients were included in the present work. For inclusion in this group, the criteria were the acceptance to be followed weekly from the first episode of molecular positivity to the SARS-CoV-2 infection until they completed at least two or three consecutive episodes of negativity by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) in nasopharyngeal-oropharyngeal (NP-OP) swab samples used for the screening and viral molecular diagnosis. Blood of all patients was collected individually after the first episode of molecular positivity for further serological analyses.

Molecular Characterization

Nucleic acid extraction from all the NP-OP swabs was performed using the MagMAX™ Viral/Pathogen II (MVP II) Nucleic Acid Isolation Kit (Catalog number: A48383; Applied Biosystems, Waltham, MA, USA) and carried out according to the instructions of the manufacturer. Molecular detection of SARS-CoV-2 was performed using the AgPath-ID One-Step RT-PCR Reagents (Ambion, Austin, TX, USA). For the first screening, specific SARS-CoV-2 primers and probes were applied to the Envelope (E) gene, followed by the detection of the RNA-dependent RNA polymerase (RdRp) gene (17), as recommended by the WHO (18, 19). The screening study included some patients who were positive for the SARS-CoV-2 E gene but were considered as having an inconclusive diagnosis for COVID-19 due to negative results for the RdRp gene. They were followed up according to the same criteria used for confirmed cases (Supplementary Table 1). To detect sgRNA, we applied one of the Envelope primer and the probe in combination with other primer previously described (20). The RT-qPCR reactions consisted of a step of reverse transcription at 45°C for 10 min, enzyme activation at 95°C for 10 min, and 50 cycles at 95°C for 15 sec, and 60°C for 45 sec for hybridization and extension using QuantStudio™ 3 Real-Time PCR System (Thermo Fisher Scientific Inc., Waltham, MA, USA) to collect a fluorescence signal at the end of each cycle. We used the isolated virus SP02/BRA as the positive control, which was kindly provided by the Laboratory of Clinical Virology, Institute of Biomedical Sciences, University of São Paulo (21).

Serological Characterization

We assessed patients anti-SARS-CoV-2 antibody production performing an enzyme-linked immunosorbent assay (ELISA), using sera of positive patients to analyze the presence of specific IgA and IgG against the viral nucleocapsid C-terminal portion, which was kindly provided by Prof. L. C. S. Ferreira (Institute of Biomedical Sciences, University of São Paulo, São Paulo, Brazil).

MaxiSorp plates (Nalge Nunc International, Rochester, NY, USA) were coated with antigen (375 ng in 50 μl of 1 × phosphate-buffered saline [PBS] per well) and incubated overnight at 4°C. We used 1 × PBS plus 2.5% of heath inactivated Fetal Bovine Serum as a blocking buffer (PBS-FBS −200 μl/well), incubating for 1 h at room temperature (RT). After, the diluted sera (1:100 in a PBS-FBS) was added (50 μl/well) and incubated at RT for 2 h. Subsequently, secondary peroxidase-conjugated anti-IgA (1:16,000 in PBS-FBS; Sigma Aldrich, St Luis, MO, USA) or anti-IgG (1:8,000 in PBS-FBS; Sigma Aldrich, St Luis, MO, USA) antibodies were added to each well (75 μl/well) and incubated at 37°C for 50 min. To reveal the reaction, 100 μl/well of the reagent 3,3′,5,5′-TetraMethylBenzidine (TMB) (Thermo Fisher Scientific Inc., Waltham, MA, USA) were added and incubated in the dark for 10 min at RT. Then, the enzymatic reaction was stopped using 100 μl/well of 0.2 N H2SO4. Between the blocking steps and the enzymatic reaction, four washing steps were performed using 1 × PBS with 0.05% Tween 20. We used a Multiskan™ FC Microplate Photometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) to read the reaction, considering the difference between the optical density (OD) at 450 and 620 nm.

Receptor-Binding Inhibition Assay (Surrogate Virus Neutralization Test Kit, SVNT)

The ability of the patient to neutralize SARS-CoV-2 was assessed through a sVNT (Catalog number: L00847-A; GenScript, Piscataway, NJ, USA), following manufacturer instructions (22). Their sera were incubated with a recombinant receptor-binding domain conjugated with horseradish peroxidase (RBD-HRP) at the final dilution of 1:20 for 30 min at 37°C. Afterward, the combination was added to a plate pre-coated with human ACE2 (hACE2) and incubated for 15 min at 37°C. Unbound RBD-HRP was removed by four washes, using the provided washing solution. The colorimetric signal was developed on the enzymatic reaction of HRP with 100 μl with TMB. Then, 50 μl of provided stop solution was added, and the OD at 450 nm using Multiskan™ FC Microplate Photometer (Thermo Fisher Scientific Inc., Waltham, MA, USA). Percentage of inhibition (%) = (1—Sample OD value/Negative Control OD value) × 100. The sVNT inhibition was settled as positive when the percentage value resulted in inhibition ≥ 20% (98.9% sensitivity and 100% specificity), as previously validated (22).

Virus Isolation

Samples with positive RT-qPCR NP-OP were submitted to viral isolation. To do this, Vero CCL81 cells were seeded in 24-well plates and incubated overnight at 37°C, using a 5% CO2 incubator. The initial inoculum (passage 1) was prepared to dilute the clinical sample (1:5) in non-supplemented Dulbecco's Modified Eagle Medium (DMEM) low glucose media (LGC Biotecnologia, São Paulo, Brazil). The inoculum was then added to the monolayer cells, homogenized, and maintained for 1 h at 37°C. The inoculum was then removed, cells were washed two times with warm PBS and fresh DMEM low glucose media supplemented with FBS (2%) and penicillin/streptomycin (1%, 10,000 U/ml) was reloaded. The cell culture was observed for a cytopathic effect every day after the inoculation. Seventy-two-hour post-infection the cell culture supernatant was collected, centrifuged for the removal of debris (500 g, 10 min, 4°C) and stored at −80°C as the first viral passage. The same strategy was employed three times to obtain a third isolate passage. To confirm the isolation, we submitted the supernatants of each passage to RT-qPCR molecular assay to detect the Envelope gene, as mentioned above. Viral samples were considered isolated when the Ct value has dropped between passages one and three. We considered samples with Ct below 35 suitable for the next isolation passages. Samples in which the Ct value has not dropped were considered as non-isolated samples.

Sequencing and Viral Genome Assembly

Based on the Ct values of the NP-OP samples, the total RNA was extracted using the MagMAX™ Viral/Pathogen II (MVP II) Nucleic Acid Isolation Kit (Catalog number: A48383) (Applied Biosystems) and carried out according to the instructions of the manufacturer. The target whole viral genome library preparation was constructed using the QIAseqTM SARS-CoV-2 Primer Panel (Catalog number: 333896; Qiagen, Hilden, Germany). Sequencing was done at the Core Facility to Support Research—University of São Paulo (CEFAP-USP/GENIAL) using the Illumina MiSeq platform. Each sample was barcoded individually, which allowed the separation of reads for each one of them. Short unpaired reads and low-quality bases and reads were removed using Trimmomatic version 0.39 (LEADING:20 TRAILING:20 SLIDINGWINDOW:4:25 MINLEN:36) (23). Consensus genomes were assembled with paired end reads using Bowtie2 version 2.0.6 using default parameters (24). All the datasets analyzed during the current study are available from the corresponding author on reasonable request. The new sequence here characterized was deposited in GenBank under the accession number MW495017.

SARS-CoV-2 Lineage Nomenclature and Phylogenetic Analysis

The viral genome sequences obtained from the NP-OP samples and assembled in this study were submitted to the SARS-CoV-2 lineage assignment using the Phylogenetic Assignment of Named Global Outbreak LINeages web application following a methodology previously described (9). Phylogenetic tree of SARS-CoV-2 based on full-length, curated sequences was estimated using the Maximum Likelihood (ML) method implemented in IQ-TREE 1.5.5 (25) with automatic model selection by ModelFinder and using the Bayesian Information Criterion (BIC) (26). To characterize the sequences of the atypical patient studied, we retrieved all complete sequences in GISAID (June 28, 2021), collected in Brazil for SARS-CoV-2 lineage B.1.1.28 with high coverage and with collection date completed. In total, 1,303 sequences were retrieved (Supplementary Figure 1). Sequences were initially resampled, excluding all identical sequences (https://biopython.org/wiki/Sequence_Cleaner), and after separation by Brazilian states, the outliers were resampled considering the average of the low sampled Brazilian states + 2 x SD (n = 53 sequences; https://web.expasy.org/decrease_redundancy/). In the end, we worked with a final dataset containing 56 sequences, 53 from São Paulo state and 2 sequences from the atypical patient, and we included the reference genome as an outgroup (NC_045512.2). The robustness of the groupings observed was assessed using 1,000 non-parametric bootstrap replicates. ML tree was visualized and plotted using FigTree v.1.4.3 (27). All taxon labels for sequences used in this work are presented in the format: hCoV-19/local of isolation/strain name/date of isolation.

Statistical Analysis

The Shapiro–Wilk normality test was used in the present study. To compare the differences in the duration of positivity and symptoms according to each biological sex, the Mann–Whitney U test was applied, assuming a significance level of 0.05. All the statistical analysis and most of the graphs presented here were performed using R scripts and are available upon request.

Results

From February 2020 to the epidemiological week 45 (first week of November 2020), Brazil had intense SARS-CoV-2 dissemination characterized in the country as the first pandemic phase. In this period, we had monitored symptomatic individuals to characterize the SARS-CoV-2 infection status of the patients. Once positive on the molecular assay, we monitored the patients weekly by collecting NP-OP swabs and blood derivatives, such as whole blood, serum, and plasma. We found a positive correlation by the detection of Envelope and RdRp genes (r = 0.89) in the first sample collected in the molecular assay (Figure 1A). Thirty-eight mild symptomatic patients were followed weekly from the first episode of positivity until they completed at least two or three consecutive episodes of negativity by RT-qPCR in the NP-OP. The duration of the positivity after the onset of symptoms was not significantly different (p-value = 0.3531) between female (average 22.67 days ± SD 19.88) and male (average 33.34 ± SD 41.71) biological sexes (Figure 1B). The duration of symptoms in females (average 32 days ± SD 11.51) and males (average 34 days ± SD 18.14; Figure 1C) was similar (p-value = 0.8887). The length of time of the symptoms was weakly correlated with the length of the positivity (Figure 1D). Among the total number of surveyed patients, 7.9% (3/38) of them were RT-qPCR positive for an “atypical” period and were considered as “outliers” (Figure 1B), such as one female (71 days) and two male (81 and 232 days) individuals. The positive symptomatic patients presented a wide spectrum of mild-illness symptoms, the most frequent being coryza, ageusia and/or anosmia, cough, fever, headache, myalgia, nasal congestion, and presence of sputum (Figure 2).

FIGURE 1
www.frontiersin.org

Figure 1. Molecular characterization of the studied individuals during the first SARS-CoV-2 pandemic phase in Brazil. The different colors represent the female and male individuals. (A) Correlation of the SARS-CoV-2 detection to Envelope and RdRp genes. (B) Duration of positivity on Envelope molecular assay after the symptom onset. Boxplots represent the 75th percentile, median, 25th percentile, and the whiskers extend to the highest and lowest value in the 1.5 × interquartile range. (C) Duration of symptomatology since the onset of symptoms. Boxplots represent the 75th percentile, median, 25th percentile, and the whiskers extend to the highest and lowest value in the 1.5 × interquartile range. (D) Correlation of the duration of symptoms to the duration of positivity after the onset of symptoms. SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.

FIGURE 2
www.frontiersin.org

Figure 2. Frequency of each symptom of the Coronavirus Disease 2019 in the 38 SARS-CoV-2 positive patients. Normalized frequency considering the presence of the symptom with the lowest (n = 1) and highest (n = 26) value, to values between 0 and 1. SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.

One of those outliers presented a very peculiar picture regarding positivity over time. The male patient, 38-year-old, was previously diagnosed with HIV without AIDS as confirmed by negative molecular diagnosis for HIV infection and normal levels of T CD4+ and CD8+ lymphocytes both before and after the infection with SARS-CoV-2 (Supplementary Table 2). The patient did not present any other comorbidity, started the clinical symptoms compatible with the classic condition caused by COVID-19 on April 21, 2020, notably fever, asthenia, headache, cough, coryza, ageusia, and anosmia until May 11, 2020 (Figure 3A). In total, the symptoms lasted for 20 days. The patient was diagnosed as positive to the SARS-CoV-2 by RT-qPCR on the seventh day after the onset of symptoms (Figure 3B), with the lineage B.1.1.28, as ascertained by sequencing (Figure 4). Up to 232 days after the onset of symptoms, the patient was still positive in the RT-qPCR assay for the presence of vRNA (Figure 3B). Furthermore, two of the collected samples were also positive for sgRNA, the first (April 27, 2020) and the one before the last sample (December 1, 2020) (Supplementary Figure 2), suggesting at least 224 days of viral replicative activity. Before the last positivity collection point, the patient had four negative harvestings in three moments of the infectious process, the last being for 2 consecutive weeks in November 2020 (Figure 3B). None of the samples was positive for viral isolation in cell culture.

FIGURE 3
www.frontiersin.org

Figure 3. Timeline infection of the atypical studied patient with prolonged detection of viral RNA. (A) Schematic figure containing the diagnostic points, harvesting, and symptoms over time. (B) Cycle threshold to the Envelope gene of the SARS-CoV-2 according to the timeline of infection, as determined by RT-qPCR. The NP-OP swabs positive to the detection of subgenomic RNA are indicated by the dashed circle. (C) ELISA data showing the detection of IgA and IgG specific anti- SARS-CoV-2 nucleocapsid according to the timeline of infection. (D) Inhibition of the SARS-CoV-2/ACE2 ligation mediated by receptor-binding domain by neutralizing antibodies, according to the timeline of infection. SARS-CoV-2, severe acute respiratory syndrome coronavirus 2; NP-OP, nasopharyngeal-oropharyngeal; ELISA, enzyme-linked immunosorbent assay.

FIGURE 4
www.frontiersin.org

Figure 4. Maximum likelihood phylogenetic tree for SARS-CoV-2 lineage B.1.1.28 based on full-length genome sequences (n = 56). The sequence NC_045512.2 was used as an outgroup to reroot the tree and the nodes are labeled with bootstrap support values higher than 50. The two sequences obtained from the atypical patient are highlighted in yellow. SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.

Serological samples of this atypical patient were then tested in ELISA for the presence of anti-nucleocapsid (Nc) IgA and IgG antibody levels. We observed that specific anti-Nc IgA was detected from 3 weeks after the onset of symptoms (May 11, 2020) and persisted for the next two sera sample harvestings, (May 14, 2020 and May 27, 2020), remaining at basal levels until the last blood collection of the studied period (Figure 3C). On the other hand, anti-Nc IgG levels were detectable during the entire course of infection, with high OD values throughout the analyzed period, being the highest values observed in the first 2 months of the study, corresponding to the points with the greatest amount of vRNA (Figure 3B). Inhibition assays showed that total neutralizing antibody levels increased in the first 14–21 days of positivity for the molecular envelope gene assay, and around 85% of inhibition persisted during the sampling period (Figure 3D). In accordance, a perfect correlation was observed between the peaks of serum neutralizing antibodies and the decrease of the viral charge at the beginning of the asymptomatic period and the negative amplification of the envelope gene by RT-qPCR every since.

Among all NP-OP positive samples, the genomic monitoring was performed in three positive samples of this patient to perform genomic sequencing, one at the beginning (collection date = 27-Apr-2020), one at the middle (collection date = August 18, 2020), and one at the end (collection date = December 1, 2020) of the positivity period, all of them been chosen based on the combination of Ct values and period of infection (Figure 4). We obtained a complete genome for the first one, and for the other two samples only partial genomes. All sequences were characterized as lineage B.1.1.28. The last sequenced sample (collection date = December 1, 2020) had low coverage and was discarded from the phylogenetic analysis. The two sequences included in the phylogenetic studies were grouped in a monophyletic group when compared to other sequences isolated in São Paulo, Brazil for the same lineage and diversified between the first and the last atypical patient's harvestings, suggesting a within-host evolution of the SARS-CoV-2 (Figure 4).

Discussion

The clinical manifestations induced by SARS-CoV-2 vary from a broad spectrum of symptoms, ranging from asymptomatic or mild-to-severe cases of COVID-19 (2831). In autopsy investigations, the SARS-CoV-2 was identified in several organs, such as the lung, heart, and kidney (31). The viral infection has been associated with respiratory, gastrointestinal, hepatic, cardiac, renal, and neurological dysfunctions (30, 32), which leads to a multisystem inflammatory disease, acute respiratory distress syndrome, multiple organ failure, and death in severe cases (30, 32, 33). The dysfunctions induced by COVID-19 are manifested by signs and symptoms, which when well elucidated can be crucial to the clinical diagnosis and treatment (34). Initial descriptions of the average period of vRNA presence in samples of the patients who survived the infection are around 20 days, although other studies have identified atypical cases that can last up to 154 days (11, 12).

In the present study, the duration of positivity to SARS-CoV-2 by RT-qPCR was not statistically significantly different comparing biological sexes, with an average of 33.34 days for males and 22.67 days for females. However, one outlier was noted within the female group and two males. One male outlier patient remained detectable for SARS-CoV-2 for 232 days from the onset of symptoms. It is worth noting that a suggestive replicative viral activity observed in this particular patient was detectable at least until 224 days after the beginning of the infection, with lower Ct values even in advanced phases of the COVID-19 what may be of epidemiological significance in terms of transmission. This was clearly shown by the correlation between the detection of both sgRNA and the genomic envelope gene of SARS-CoV-2 (Supplementary Figure 2). Although the clinical and molecular conditions caused by viruses characterized as acute infections are generally quickly resolved, some of them can cause a long-lasting infection in mild (35) or severe conditions (36), and the impact in transmissibility and pathology still needs to be understood. Thus, it is very important to clarify whether a prolonged vRNA shedding is correlated to the clinical outcome of the infection, or yet to the condition of the patients, such as immunosuppressive states or occurrence of comorbidities.

The serological data of the atypical individual analyzed here above indicated that during the period of molecular vRNA detection, it was possible to detect specific IgA and IgG antibodies against SARS-CoV-2 Nc antigen. These data are compatible with the development of a humoral immune response raised against the virus very early after infection. Our results showed that both IgA and IgG anti-SARS-CoV-2 were already present in the first collected serum sample, 20 days after the onset of the symptoms. As with other viral infections (37, 38), the detection of serum IgA antibodies was limited in time, while the presence of IgG antibodies could be identified in the serum for long periods and present during the entire period of the study. Additionally, from the beginning of the sera collection, we identified the presence of total neutralizing antibodies (which can be of all immunoglobulin classes), with fluctuations during sampling. Despite an unrestrained innate immunity that can account for virus clearance at the beginning of acute infections (39), the specific responses observed here may reflect an early attempt of the immune system to control the viral infection by the induction of a robust neutralizing humoral immune response.

On top of the development of humoral and cellular immune responses, viruses escape from other immunological barriers imposed to clear the infection. Usually, to avoid an immune-mediated viral clearance, viruses frequently use a combination of several different strategies to subvert recognition by the immune system, such as (i) settling latent infections, (ii) replicating in immune-privileged sites, (iii) downregulating the expression of immune recognition signals on the surface of infected cells, or (iv) undergoing antigenic variation or else mechanisms for suppressing the immunological response (40). Possibly, a combination of one or more of these mechanisms could explain the prolonged viral shedding in the studied patient.

Nevertheless, the integration of the reverse-transcribed RNA into the genome of cultured human cells with the possibility of being expressed in tissue cells derived from humans seems to be an explanation, contradictorily proposed to clarify the persistent detection of SARS-CoV-2 RNA after COVID-19 recovery (41). This hypothesis does not seem to be the case of the outlier male patient studied here, since (i) we showed that sgRNA, a molecular marker for active SARS-CoV-2 replication (42, 43) was found during the long-term infection, (ii) there was a fluctuation in the levels of virus neutralization over time possibly characterizing the dynamics of the humoral immune response directed to the RBD antigen that binds to the human ACE2 receptor and yet, and (iii) the genetic divergence accumulated throughout the viral infection process (Figure 4) indicating an intra-host evolution over time. In addition, although the viral shedding was prolonged, the elimination of the viral infection was ultimately observed.

Conversely, since the outlier male patient studied here is seropositive for HIV since 2018, it could be hypothesized that this comorbidity would be impacting an appropriate immune response against SARS-CoV-2 through possible mechanisms of immunosuppression. It is worth noting that the differential lymphocytic markers to monitor HIV infection of this patient were consistently normal since the beginning of the antiretroviral therapy in 2019. Still, in-depth implications of HIV infection on the overall functioning of the immune system, even in treated individuals, are not fully comprehended (44), and data published elsewhere show that during the acute phase of HIV, there is a critical loss of memory CD4+ lymphocytes mainly in lymphoid tissues (45), which are critical to the maintenance of a fully competent immune system. The post-effects of this massive CD4+ loss may range to increased cell turnover and disrupted activation/differentiation and maturation of immune system cells, possibly due to indirect effects of the HIV replication (44, 46). Also, post-acute phase HIV causes a dramatic skewing of the lymphocyte population that is not fully recovered after effective antiretroviral treatment (47), possibly due to a decreased thymic functioning and HIV-induced lymph nodes architecture changes (48) that may impact innate and adaptive immune responses.

Thus, it remains to be concluded how the effects of HIV infection could be impacting positively or negatively the clearance of the SARS-CoV-2. In this sense, it is well recognized that an appropriate cell-mediated immune response mediated by the activation of CD4+ and CD8+ T cells producing class I interferon, i.e., interferon (IFN)-α and IFN-β, and the activation of B cells producing neutralizing antibodies are associated with a favorable clinical outcome of COVID-19 (4952). However, other data have correlated HIV treatments to the increased incidence of diabetes mellitus, hyperglycemia, and drug-altered metabolism (53), known as rather enhancing the SARS-CoV-2 replication (54). Altogether, our data suggest that if the previously existing HIV-positive status is somehow related to a prolonged infection, it may instead be favoring the emergence of new variants.

Finally, since the severity of the clinical condition of patients with COVID-19 or the resistance status to the infection is previously shown to be associated with a genetic background of the host (55, 56), it could also be hypothesized that the prolonged viral shedding is related to the host genetic environment. Still, genetic factors associated with the SARS-CoV-2 viruses or specific for certain viral lineages could certainly influence viral shedding. Other linked factors to the host, such as age, other comorbidities association, nutritional state, and previous exposure to different pathogens, could interfere with such a differential behavior of the SARS-CoV-2 during the infectious process and affect healing and transmission.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/genbank/, MW495017.

Ethics Statement

The studies involving human participants were reviewed and approved by Human Research Ethics Committee of the Biomedical Sciences Institute—University of São Paulo (CEPSH-ICB-USP) (n°. 4.036.252). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

MC, AV, LG, AC, and PM: conceptualization. MC and SM: data curation. MC, AV, CM, SA, SM, VB, SB, LO, YA, and AC: formal analysis. PM: funding acquisition. MC and PM: project administration. MC, AV, CM, SA, SM, VB, SB, LO, YA, LG, AC, and PM: investigation, methodology, resources, and writing—review and editing. MC: software. JP and PM: supervision. MC, AV, SM, and PM: validation and writing—original draft. MC: visualization. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the São Paulo Research Foundation (FAPESP—process #2017/27131-9), the Institut Pasteur (#DI2020-16/18), and the Cooperation and Cultural Action Services of the São Paulo French Consulate (#185BRA1184-−2020/SPPU). MC, AV, CM, SA, and VB received FAPESP Grants Nos. 2019/24518-5, 2019/20708-4, 2020/01487-4, 2017/23519-2, and 2018/11612-0, respectively. LO and YA received a Coordination for the Improvement of Higher Education Personnel (CAPES) Grant No. 88887.423542/2019-00 and 88887.474625/2020-00, respectively. JP is PI of a G4-grant from Institut Pasteur (#FUSP-3303-01). LG received a Brazilian National Council of Scientific and Technological Development (CNPq) Grant No. 152365/2019-2. AC was funded by Institut Pasteur/FUSP-Project #3303.

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.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors are indebted to Luis Carlos de S. Ferreira and Robert Andreata-Santos from the Institute of Biomedical Sciences, USP, for providing the SARS-CoV-2 Nc protein. We thank the Core Facility for Scientific Research—University of São Paulo (CEFAP-USP/GENIAL) for excellent technical assistance.

Supplementary Material

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

References

1. Organization WH. Pneumonia of Unknown Cause—China. (2020). Available online at: https://www.who.int/csr/don/05-january-2020-pneumonia-of-unkown-cause-china/en/

2. Zhou P, Yang X, Lou, Wang XG, Hu B, Zhang L, Zhang W, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. (2020) 579:270–3. doi: 10.1038/s41586-020-2012-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. (2020) 382:727–33. doi: 10.1056/NEJMoa2001017

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Li X, Wang W, Zhao X, Zai J, Zhao Q, Li Y, et al. Transmission dynamics and evolutionary history of 2019-nCoV. J Med Virol. (2020) 92:501–11. doi: 10.1002/jmv.25701

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Jesus JG de, Sacchi C, Candido D da S, Claro IM, Sales FCS, Manuli ER, et al. Importation and early local transmission of COVID-19 in Brazil, 2020. Rev Inst Med Trop Sáo Paulo. (2020) 62:1–5. doi: 10.1590/s1678-9946202062030

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Resende PC, Delatorre E, Graf T, Mir D, Motta FC, Appolinario L, et al. Genomic surveillance of SARS-CoV-2 reveals community transmission of a major lineage during the early pandemic phase in Brazil. bioRxiv. (2020) 6:1–18. Available at: http://biorxiv.org/content/early/2020/06/18/2020.06.17.158006.abstract

Google Scholar

7. Nascimento VA, Do, Corado A de LG, Nascimento FO Do, da Costa ÁKA, Duarte DCG, Luz SLB, Gonçalves LMF, et al. Genomic and phylogenetic characterisation of an imported case of SARS-CoV-2 in amazonas state, Brazil. Mem Inst Oswaldo Cruz. (2020) 115:1–6. doi: 10.1590/0074-02760200310

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Naveca FG, Nascimento V, de Souza VC, Corado A de L, Nascimento F, Silva G, et al. COVID-19 in Amazonas, Brazil, was driven by the persistence of endemic lineages and P.1 emergence. Nat Med. (2021) 21:1–9. doi: 10.1038/s41591-021-01378-7

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Rambaut A, Holmes EC, O'Toole Á, Hill V, McCrone JT, Ruis C, et al. 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

10. Siqueira JD, Goes LR, Alves BM, Carvalho PS de, Cicala C, Arthos J, et al. SARS-CoV-2 genomic analyses in cancer patients reveal elevated intrahost genetic diversity. Virus Evol. (2021) 7:13. doi: 10.1093/ve/veab013

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kemp SA, Collier DA, Datir RP, Ferreira IATM, Gayed S, Jahun A, et al. SARS-CoV-2 evolution during treatment of chronic infection. Nature. (2021) 592:277–82. doi: 10.1038/s41586-021-03291-y

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Choi B, Choudhary MC, Regan J, Sparks JA, Padera RF, Qiu X, et al. Persistence and evolution of SARS-CoV-2 in an Immunocompromised Host. N Engl J Med. (2020) 383:2291–3. doi: 10.1056/NEJMc2031364

PubMed Abstract | CrossRef Full Text | Google Scholar

13. San JE, Ngcapu S, Kanzi AM, Tegally H, Fonseca V, Giandhari J, et al. Transmission Dynamics of SARS-CoV-2 Within-Host diversity in two major hospital outbreaks in South Africa. Virus Evol. (2021) 7:41. doi: 10.1093/ve/veab041

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Beran A, Zink E, Mhanna M, Abugharbyeh A, Hanrahan J, Duggan J, et al. Transmissibility and viral replication of SARS-COV-2 in immunocompromised patients. J Med Virol. (2021) 93:4156–60. doi: 10.1002/jmv.26970

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tarhini H, Recoing A, Bridier-nahmias A, Rahi M, Lambert C, Martres P, 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. J Infect Dis. (2021) 223:1522–7. doi: 10.1093/infdis/jiab075

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Decker A, Welzel M, Laubner K, Grundmann S, Kochs G, Panning M, et al. Prolonged SARS-CoV-2 shedding and mild course of COVID-19 in a patient after recent heart transplantation. Am J Transplant. (2020) 20:3239–45. doi: 10.1111/ajt.16133

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Corman VM, Landt O, Kaiser M, Molenkamp R, Meijer A, Chu DKW, et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Eurosurveillance. (2020) 25:2000045. doi: 10.2807/1560-7917.ES.2020.25.3.2000045

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Corman V, Bleicker T, Brünink S, Drosten C, Landt O, Koopmans M, et al. Diagnostic detection of 2019-nCoV by real-time RT-PCR—Protocol and preliminary evaluation as of Jan 17, 2020. (2020). Available online at: https://virologie-ccm.charite.de/en/ (accessed September 14, 2021).

Google Scholar

19. Pan American Health Organization. Laboratory Guidelines for Detection and Diagnosis of the Novel Coronavirus (2019-nCoV) Infection. (2020). Available online at: https://www.paho.org/en/documents/laboratory-guidelines-detection-and-diagnosis-novel-coronavirus-2019-ncov-infection (accessed September 14, 2021).

20. 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. doi: 10.1038/s41586-020-2196-x

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Araujo DB, Machado RRG, Amgarten DE, Malta F de M, Araujo GG de, Monteiro CO. SARS-CoV-2 isolation from the first reported patients in Brazil and establishment of a coordinated task network. Mem Inst Oswaldo Cruz. (2020) 74:1–27. doi: 10.1590/0074-02760200342

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Tan CW, Chia WN, Qin X, Liu P, Chen MIC, Tiu C, et al. A SARS-CoV-2 surrogate virus neutralization test based on antibody-mediated blockage of ACE2–spike protein–protein interaction. Nat Biotechnol. (2020) 38:1073–8. doi: 10.1038/s41587-020-0631-z

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. (2014) 30:2114–20. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. (2012) 9:357–9. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Nguyen L-T, Schmidt HA, Haeseler A, von, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies Mol Biol Evol. (2015) 32:268–74. doi: 10.1093/molbev/msu300

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Kalyaanamoorthy S, Minh BQ, Wong TKF, Von Haeseler A, Jermiin LS. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat Methods. (2017) 14:587–9. doi: 10.1038/nmeth.4285

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Rambaut A. FigTree Version 1.4.3. (2016). Available online at: http://tree.bio.ed.ac.uk/software/figtree/

28. Guan W, Ni Z, Hu Y, Liang W, Ou C, He J, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. (2020) 382:1708–20. doi: 10.1056/NEJMoa2002032

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Arunachalam PS, Wimmers F, Mok CKP, Perera RAPM, Scott M, Hagan T, et al. Systems biological assessment of immunity to mild versus severe COVID-19 infection in humans. Science. (2020) 369:1210–20. doi: 10.1126/science.abc6261

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Duarte-Neto AN, Monteiro RAA, da Silva LFF, Malheiros DMAC, de Oliveira EP, Theodoro-Filho J, et al. Pulmonary and systemic involvement in COVID-19 patients assessed with ultrasound-guided minimally invasive autopsy. Histopathology. (2020) 77:186–97. doi: 10.1111/his.14160

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Duarte-Neto AN, Caldini EG, Gomes-Gouvêa MS, Kanamura CT, de Almeida Monteiro RA, Ferranti JF, et al. An autopsy study of the spectrum of severe COVID-19 in children: From SARS to different phenotypes of MIS-C. EClinicalMedicine. (2021) 35:100850. doi: 10.1016/j.eclinm.2021.100850

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Dolhnikoff M, Ferreira Ferranti J, de Almeida Monteiro RA, Duarte-Neto AN, Soares Gomes-Gouvêa M, Viu Degaspare N, et al. SARS-CoV-2 in cardiac tissue of a child with COVID-19-related multisystem inflammatory syndrome. Lancet Child Adolesc Heal. (2020) 4:790–4. doi: 10.1016/S2352-4642(20)30257-1

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Tay MZ, Poh CM, Rénia L, MacAry PA, Ng LFP. The trinity of COVID-19: immunity, inflammation and intervention. Nat Rev Immunol. (2020) 20:363–74. doi: 10.1038/s41577-020-0311-8

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zhou X, Menche J, Barabási AL, Sharma A. Human symptoms-disease network. Nat Commun. (2014) 5:1–10. doi: 10.1038/ncomms5212

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Oliveira DBL, Durigon GS, Mendes ÉA, Ladner JT, Andreata-Santos R, Araujo DB, et al. Persistence and intra-host genetic evolution of zika virus infection in symptomatic adults: a special view in the male reproductive system. Viruses. (2018) 10. doi: 10.3390/v10110615

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Oliveira DBL, Almeida FJ, Durigon EL, Mendes ÉA, Braconi CT, Marchetti I, et al. Prolonged Shedding of Zika Virus Associated with Congenital Infection. N Engl J Med. (2016) 15:1386–8. doi: 10.1056/NEJMc1500960

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Carter I, Fraser J, Cloonan M. Specific IgA antibody response in Ross River virus infection. Immunol Cell Biol. (1987) 65:511–3. doi: 10.1038/icb.1987.60

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Hsueh PR, Huang LM, Chen PJ, Kao CL, Yang PC. Chronological evolution of IgM, IgA, IgG and neutralisation antibodies after infection with SARS-associated coronavirus. Clin Microbiol Infect. (2004) 10:1062–6. doi: 10.1111/j.1469-0691.2004.01009.x

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Schuch V, Martins F, Ten Caten F, Araujo-Pereira M, Cunha MDP, El Khawanky N, et al. Systems Immunology of Flavivirus Infection. London: Academic Press. (2021).

Google Scholar

40. Borrow P. Mechanisms of viral clearance and persistence. J Viral Hepat. (1997) 4:16–24. doi: 10.1111/j.1365-2893.1997.tb00176.x

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Zhang L, Richards A, Barrasa MI, Hughes SH, Young RA, Jaenisch R. Reverse-transcribed SARS-CoV-2 RNA can integrate into the genome of cultured human cells and can be expressed in patient-derived tissues. Proc Natl Acad Sci. (2021) 118:e2105968118. doi: 10.1073/pnas.2105968118

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Avanzato VA, Matson MJ, Seifert SN, Pryce R, Williamson BN, Anzick SL, et al. Case Study: Prolonged Infectious SARS-CoV-2 Shedding from an Asymptomatic Immunocompromised Individual with Cancer. Cell. (2020) 183:1901–12. doi: 10.1016/j.cell.2020.10.049

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Romano M, Ruggiero A, Squeglia F, Maga G, Berisio R. A structural view of SARS-CoV-2 RNA replication machinery: rna synthesis, proofreading and final capping. Cells. (2020) 9:1267. doi: 10.3390/cells9051267

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Moir S, Chun TW, Fauci AS. Pathogenic mechanisms of HIV disease. Annu Rev Pathol Mech Dis. (2011) 6:223–48. doi: 10.1146/annurev-pathol-011110-130254

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Guadalupe M, Reay E, Sankaran S, Prindiville T, Flamm J, McNeil A, et al. Severe CD4+ T-Cell depletion in gut lymphoid tissue during primary human immunodeficiency virus type 1 infection and substantial delay in restoration following highly active antiretroviral therapy. J Virol. (2003) 77:11708. doi: 10.1128/JVI.77.21.11708-11717.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Lewin-Smith MR, Klassen MK, Frankel SS, Nelson AM. Pathology of human immunodeficiency virus infection: infectious conditions. Ann Diagn Pathol. (1998) 2:181–94. doi: 10.1016/S1092-9134(98)80006-3

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Hellerstein MK, Hoh RA, Hanley MB, Cesar D, Lee D, Neese RA, et al. Subpopulations of long-lived and short-lived T cells in advanced HIV-1 infection. J Clin Invest. (2003) 112:956. doi: 10.1172/JCI200317533

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Schacker TW, Nguyen PL, Beilman GJ, Wolinsky S, Larson M, Reilly C, et al. Collagen deposition in HIV-1 infected lymphatic tissues and T cell homeostasis. J Clin Invest. (2002) 110:1133. doi: 10.1172/JCI0216413

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Shah VK, Firmal P, Alam A, Ganguly D, Chattopadhyay S. Overview of Immune Response During SARS-CoV-2 Infection: Lessons From the Past. Front Immunol. (2020) 11:1949. doi: 10.3389/fimmu.2020.01949

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Park A, Iwasaki A. Type I and Type III Interferons—Induction, Signaling, Evasion, and Application to Combat COVID-19. Cell Host Microbe. (2020) 27:870–878. doi: 10.1016/j.chom.2020.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Sa Ribero M, Jouvenet N, Dreux M, Nisole S. Interplay between SARS-CoV-2 and the type I interferon response. PLoS Pathog. (2020) 16:e1008737. doi: 10.1371/journal.ppat.1008737

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Vanderheiden A, Ralfs P, Chirkova T, Upadhyay AA, Zimmerman MG, Bedoya S, et al. Type I and Type III interferons restrict SARS-CoV-2 infection of human airway epithelial cultures. J Virol. (2020) 94:e00985–20. doi: 10.1128/JVI.00985-20

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Montessori V, Press N, Harris M, Akagi L, Montaner JSG. Adverse effects of antiretroviral therapy for HIV infection. C Can Med Assoc J. (2004) 170:229.

PubMed Abstract | Google Scholar

54. Codo AC, Davanzo GG, Monteiro L de B, Souza GF de, Muraro SP, Virgilio-da-Silva JV, Prodonoff JS, Carregari VC, Junior CAO de B, Crunfli F, et al. Elevated Glucose Levels Favor SARS-CoV-2 infection and monocyte response through a hif-1α/glycolysis-dependent axis. Cell Metab. (2020) 32:437. doi: 10.1016/j.cmet.2020.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Alessio A, Claudio A, Salvatore B, Luca B, Monica B, Maurizio C, et al. Genomewide association study of severe covid-19 with respiratory failure. N Engl J Med. (2020) 383:1522–34. doi: 10.1056/NEJMoa2020283

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Zhang Q, Bastard P, Liu Z, Pen J, Le, Moncada-Velez M, Chen J, et al. . Inborn errors of type I IFN immunity in patients with life-threatening COVID-19 Science. (2020) 370:1–13. doi: 10.1126/science.abd4570

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: SARS-CoV-2, prolonged RNA viral shedding, immune response, coronavirus, intra-host evolution

Citation: Cunha MP, Vilela APP, Molina CV, Acuña SM, Muxel SM, Barroso VM, Baroni S, Gomes de Oliveira L, Angelo YS, Peron JPS, Góes LGB, Campos ACA and Minóprio P (2021) Atypical Prolonged Viral Shedding With Intra-Host SARS-CoV-2 Evolution in a Mildly Affected Symptomatic Patient. Front. Med. 8:760170. doi: 10.3389/fmed.2021.760170

Received: 17 August 2021; Accepted: 13 October 2021;
Published: 26 November 2021.

Edited by:

Sanjay Kumar, Armed Forces Medical College, Pune, India

Reviewed by:

Jayanthi S. Shastri, Brihanmumbai Municipal Corporation, India
Nivedita Gupta, Indian Council of Medical Research (ICMR), India

Copyright © 2021 Cunha, Vilela, Molina, Acuña, Muxel, Barroso, Baroni, Gomes de Oliveira, Angelo, Peron, Góes, Campos and Minóprio. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Marielton dos Passos Cunha, marieltondospassos@gmail.com; Paola Minóprio, paola.minoprio@pasteur.fr

These authors have contributed equally to this work