Rapid spread and high impact of the Variant of Concern P.1 in the largest city of Brazil ======================================================================================== * Gabriela Rodrigues Barbosa * Luiz Vinicius Leão Moreira * Alberto Fernando Oliveira Justo * Ana Helena Perosa * Ana Paula Cunha Chaves * Mariana Sardinha Bueno * Luciano Kleber de Souza Luna * Danielle Dias Conte * Joseane Mayara Almeida Carvalho * Janesly Prates * Patricia Sousa Dantas * Klinger Soares Faico-Filho * Clarice Camargo * Paola Cristina Resende * Marilda Mendonça Siqueira * Nancy Bellei ## Abstract First in Manaus in the Brazilian Northern region, the Variant of Concern P.1 traveled 3800 kilometers southeast to endanger Sao Paulo contributing to the collapse of the health system. Here, we show evidence of how fast the VOC P.1 has spread in the most populated city in South America. The COVID-19 has reached over 131 million cases around the world, with more than 2.8 million of death, according to the World Health Organization (1). As the SARS-CoV-2 cases continue to emerge globally, variants of concern (VOC) and variants of interest (VOI) have been described (2). The variants share important points of mutations at the receptor-binding domain (RBD) of the Spike protein, which might increase the transmissibility of COVID-19 and promote escape from neutralizing antibodies (3). Brazil is currently the epicenter of COVID-19, with more than 13 million confirmed cases of April 2021(1). Two variants, VOC P.1 and VOI P.2, evolved from lineage B.1.1.28, have taken over the scene since late 2020 in the country (4,5). The VOC P.1 was first detected in January 2021 in Japanese travelers returning from Manaus and was responsible for the second wave in Amazonas in late November 2020 (5). In October 2020, the VOI P.2 was reported in Rio de Janeiro and was estimated to have emerged in late August 2020 (5,6). Both variants have been associated also with reinfection cases (6). In a recent investigation (data not published) of the circulating variants in São Paulo city, during the first week of March, 64.4% of samples were identified as P.1. The investigation of lineages not only contributes to epidemiological surveillance but also provides a better comprehension of the spread and circulation of variants, allowing the association to clinical outcomes and response to vaccines (7). In this sense, we aimed to investigate the spread of P.1 and P.2 variants in samples from hospitalized patients (HP) and healthcare workers (HCW) attended in a university hospital in São Paulo city. This study was conducted in compliance with institutional guidelines, approved by the Ethics Committee of São Paulo Federal University (CEP/UNIFESP n. 29407720.4.0000.5505). From March 1st to March 15th, 427 nasopharyngeal samples were collected from 245 HP and 125 from HCW outpatients (25.5% and 23.2% of positivity rate, respectively). We then selected 60 samples with Ct value ≤30 (38 samples from HP, and 22 from HCW). All HCW presented only mild symptoms and did not need hospitalization. Of the 60 selected samples, 52 whole genome sequences were generated (30 from HP and 22 from HCW) following the sequencing protocol using the Illumina MiSeq platform and the analysis pipeline described by Resende et al (8). The SARS-CoV-2 lineages were classified by the PANGO lineages nomenclature (9). Genome sequences generated have been deposited at the EpiCoV database on GISAID ([https://www.gisaid.org/](https://www.gisaid.org/)) under accession numbers EPI\_ISL\_1464630 to EPI\_ISL\_1464677. Of the 52 sequenced samples, 44 (84.4%) were identified as VOC P.1; 5 (9.2%) as VOI P.2; 1 (1,9%) as B.1.1.7, and 2 (3,8%) B.1.1.28 (Figure). The most notable variants circulating in the second wave, including B.1.1.7 (detected first in the United Kingdom) and B.1.1.351 (detected first in South Africa), and P.1, are related to an increase of transmissibility (2,10). Interestingly, the P.1 variant was first identified in the State of Amazonas, about 3,800 kilometers apart from São Paulo (5). It is evident that the P.1 variant prevailed during the first two weeks of March, showing a regular distribution among HP and HCW with no difference in terms of age, sex, vaccination, and outcome (Table). From the first to the second weeks of March, we observed a higher frequency of P.1 (78.6% and 91.7%, respectively). In this survey, only one sample from a HP was identified as VOC B.1.1.7. The other two samples were identified as B.1.1.28, a widely spread lineage during the first wave in Brazil. ![Figure1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/14/2021.04.10.21255111/F1.medium.gif) [Figure1](http://medrxiv.org/content/early/2021/04/14/2021.04.10.21255111/F1) View this table: [Table](http://medrxiv.org/content/early/2021/04/14/2021.04.10.21255111/T1) Table Characteristics of individuals identified with the VOC P.1 There is a broad discussion about whether the available vaccines against SARS-CoV-2 will be less effective at preventing infection with the emerging variants (10). In this work, 14 samples (26.9%) of the 52 sequenced samples were from individuals that had received at least one dose of vaccine, ChAdOx1-S/nCoV-19 (n=2) or SINOVAC (n=26). Although they were vaccinated, they could not be considered immunized, regarding the days after vaccination. Among the hospitalized patients, 19 (63%) were admitted to the intensive care unit, from which nine were discharged and ten died. Comparing the RT-PCR Ct values of all attended patients since the first wave, we did not observe any difference in the Ct mean values with those of P.1 (data not shown). May 2020 registered the peak of number of positive cases with a Ct mean of 23.6. Now, as of April 2021, we are facing a rise in the number of cases. However, the Ct mean was 24.9, which may indicate that the spread of P.1 does not contribute to an actual increase in the viral load. There is still a need for more epidemiologic surveys to assure the role of the VOCs in transmission and escape to neutralizing antibodies. Our findings emphasize that the P.1 variant has spread widely throughout the country. Despite all actions of interventions such as the use of masks, physical distancing, flight travel reductions, and the currently established lockdown in São Paulo, the frequency rates of P.1 increased significantly in two weeks, evidencing its fast spread. ## Data Availability The data that support the findings of this study are available in GISAID at [https://www.gisaid.org/](https://www.gisaid.org/) reference number EPI\_ISL\_1464630 to EPI\_ISL\_1464677. ## Disclosure statement No potential conflict of interest was reported by the author(s). ## Acknowledgements The authors wish to thank all the health care workers and scientists who have worked hard to deal with this pandemic threat, in special the sequencing team from Laboratory of Respiratory Viruses and Measles, Oswaldo Cruz Institute and the Fiocruz COVID-19 Genomic Surveillance Network ([http://www.genomahcov.fiocruz.br/](http://www.genomahcov.fiocruz.br/)). This work was supported by grant 2020/11719-0, São Paulo Research Foundation (FAPESP) and Coordination for the Improvement of Higher Education Personnel CAPES. * Received April 10, 2021. * Revision received April 10, 2021. * Accepted April 14, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.WHO. WHO Coronavirus Disease (COVID-19) Dashboard [Internet]. Who.int. 2021 [cited 2021 Apr 7]. Available from: [https://covid19.who.int/](https://covid19.who.int/) 2. 2.Centers for Disease Control and Prevention (US). SARS-CoV-2 Variants of Concern [Internet]. Centers for Disease Control and Prevention (US). 2021 [cited 2021 Apr 2]. Available from: [https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/variant-surveillance/variant-info.html](https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/variant-surveillance/variant-info.html) 3. 3.Resende PC. The ongoing evolution of variants of concern and interest of SARS-CoV-2 in Brazil revealed by convergent indels in the amino (N)-terminal domain of the Spike protein -SARS-CoV-2 coronavirus / nCoV-2019 Genomic Epidemiology - Virological. VirologicalOrg [Internet]. 2021;1–37. Available from: [https://virological.org/t/the-ongoing-evolution-of-variants-of-concern-and-interest-of-sars-cov-2-in-brazil-revealed-by-convergent-indels-in-the-amino-n-terminal-domain-of-the-spike-protein/659](https://virological.org/t/the-ongoing-evolution-of-variants-of-concern-and-interest-of-sars-cov-2-in-brazil-revealed-by-convergent-indels-in-the-amino-n-terminal-domain-of-the-spike-protein/659) 4. 4.Voloch CM, da Silva Francisco R, de Almeida LGP, Cardoso CC, Brustolini OJ, Gerber AL, et al. Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil. J Virol. 2021;(March). 5. 5.Naveca F, Souza V, Costa C. COVID-19 epidemic in the Brazilian state of Amazonas was driven by long-term persistence of endemic SARS-CoV-2 lineages and the recent emergence of the new Variant of Concern P. 1. Res Sq Sq [Internet]. 2021;(February):1–21. Available from: [https://www.researchsquare.com/article/rs-275494/v1](https://www.researchsquare.com/article/rs-275494/v1) 6. 6.Resende PC, Bezerra JF, Teixeira de Vasconcelos RH, Arantes I, Appolinario L, Mendonça AC, et al. Spike E484K mutation in the first SARS-CoV-2 reinfection case confirmed in Brazil, 2020 - SARS-CoV-2 coronavirus / nCoV-2019 Genomic Epidemiology - Virological [Internet]. Virological. 2021. p. 1–8. Available from: [https://virological.org/t/spike-e484k-mutation-in-the-first-sars-cov-2-reinfection-case-confirmed-in-brazil-2020/584](https://virological.org/t/spike-e484k-mutation-in-the-first-sars-cov-2-reinfection-case-confirmed-in-brazil-2020/584) 7. 7.Walensky RP, Walke HT, Fauci AS. SARS-CoV-2 Variants of Concern in the United States-Challenges and Opportunities. JAMA - J Am Med Assoc. 2021;325(11):1037–8. 8. 8.Resende PC, Motta FC, Roy S, Appolinario L, Fabri A, Xavier J, et al. SARS-CoV-2 genomes recovered by long amplicon tiling multiplex approach using nanopore sequencing and applicable to other sequencing platforms. bioRxiv. 2020;1–11. 9. 9.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 [Internet]. 2020;5(11):1403–7. Available from: [http://dx.doi.org/10.1038/s41564-020-0770-5](http://dx.doi.org/10.1038/s41564-020-0770-5) 10. 10.Wang P, Nair MS, Liu L, Iketani S, Luo Y, Guo Y, et al. Antibody Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7. Nature. 2021;