Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Feb 23;10(1):e0236621.
doi: 10.1128/spectrum.02366-21. Epub 2022 Feb 23.

Spread of Gamma (P.1) Sub-Lineages Carrying Spike Mutations Close to the Furin Cleavage Site and Deletions in the N-Terminal Domain Drives Ongoing Transmission of SARS-CoV-2 in Amazonas, Brazil

Felipe Gomes Naveca  1   2 , Valdinete Nascimento  1 , Victor Souza  1 , André de Lima Corado  1 , Fernanda Nascimento  1 , George Silva  1   3 , Matilde Contreras Mejía  1 , Maria Júlia Brandão  1 , Ágatha Costa  1 , Débora Duarte  1 , Karina Pessoa  1 , Michele Jesus  4 , Luciana Gonçalves  1   5 , Cristiano Fernandes  5 , Tirza Mattos  6 , Ligia Abdalla  7 , João Hugo Santos  8 , Alex Martins  7 , Fabiola Mendonça Chui  7 , Fernando Fonseca Val  9 , Gisely Cardoso de Melo  7   9 , Mariana Simão Xavier  9   10 , Vanderson de Souza Sampaio  5   9 , Maria Paula Mourão  7   9 , Marcus Vinícius Lacerda  9   11 , Érika Lopes Rocha Batista  12 , Alessandro Leonardo Álvares Magalhães  12 , Nathânia Dábilla  13 , Lucas Carlos Gomes Pereira  14 , Fernando Vinhal  14 , Fabio Miyajima  15 , Fernando Braga Stehling Dias  15 , Eduardo Ruback Dos Santos  16 , Danilo Coêlho  17 , Matheus Ferraz  17 , Roberto Lins  17 , Gabriel Luz Wallau  18 , Edson Delatorre  19 , Tiago Gräf  20 , Marilda Mendonça Siqueira  21 , Paola Cristina Resende  21 , Gonzalo Bello  22 ; Fiocruz COVID-19 Genomic Surveillance Network
Affiliations

Spread of Gamma (P.1) Sub-Lineages Carrying Spike Mutations Close to the Furin Cleavage Site and Deletions in the N-Terminal Domain Drives Ongoing Transmission of SARS-CoV-2 in Amazonas, Brazil

Felipe Gomes Naveca et al. Microbiol Spectr. .

Abstract

The Amazonas was one of the most heavily affected Brazilian states by the COVID-19 epidemic. Despite a large number of infected people, particularly during the second wave associated with the spread of the Variant of Concern (VOC) Gamma (lineage P.1), SARS-CoV-2 continues to circulate in the Amazonas. To understand how SARS-CoV-2 persisted in a human population with a high immunity barrier, we generated 1,188 SARS-CoV-2 whole-genome sequences from individuals diagnosed in the Amazonas state from 1st January to 6th July 2021, of which 38 were vaccine breakthrough infections. Our study reveals a sharp increase in the relative prevalence of Gamma plus (P.1+) variants, designated Pango Lineages P.1.3 to P.1.6, harboring two types of additional Spike changes: deletions in the N-terminal (NTD) domain (particularly Δ144 or Δ141-144) associated with resistance to anti-NTD neutralizing antibodies or mutations at the S1/S2 junction (N679K or P681H) that probably enhance the binding affinity to the furin cleavage site, as suggested by our molecular dynamics simulations. As lineages P.1.4 (S:N679K) and P.1.6 (S:P681H) expanded (Re > 1) from March to July 2021, the lineage P.1 declined (Re < 1) and the median Ct value of SARS-CoV-2 positive cases in Amazonas significantly decreases. Still, we did not find an increased incidence of P.1+ variants among breakthrough cases of fully vaccinated patients (71%) in comparison to unvaccinated individuals (93%). This evidence supports that the ongoing endemic transmission of SARS-CoV-2 in the Amazonas is driven by the spread of new local Gamma/P.1 sublineages that are more transmissible, although not more efficient to evade vaccine-elicited immunity than the parental VOC. Finally, as SARS-CoV-2 continues to spread in human populations with a declining density of susceptible hosts, the risk of selecting more infectious variants or antibody evasion mutations is expected to increase. IMPORTANCE The continuous evolution of SARS-CoV-2 is an expected phenomenon that will continue to happen due to the high number of cases worldwide. The present study analyzed how a Variant of Concern (VOC) could still circulate in a population hardly affected by two COVID-19 waves and with vaccination in progress. Our results showed that the answer behind that was a new generation of Gamma-like viruses, which emerged locally carrying mutations that made it more transmissible and more capable of spreading, partially evading prior immunity triggered by natural infections or vaccines. With thousands of new cases daily, the current pandemics scenario suggests that SARS-CoV-2 will continue to evolve and efforts to reduce the number of infected subjects, including global equitable access to COVID-19 vaccines, are mandatory. Thus, until the end of pandemics, the SARS-CoV-2 genomic surveillance will be an essential tool to better understand the drivers of the viral evolutionary process.

Keywords: Brazil; COVID-19; SARS-CoV-2; coronavirus; variant gamma; virus evolution.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflict of interest.

Figures

FIG 1
FIG 1
Temporal distribution and genetic diversity of SARS-CoV-2-positive samples from the Amazonas state during 2021. (A) Graph depicting the temporal evolution of SARI cases and SARI deaths based on the date of symptom onset (source, http://info.gripe.fiocruz.br) as a proxy for the COVID-19 epidemic curve in Amazonas state, along with the number of SARS-CoV-2 whole-genome sequences generated between January and July 2021. (B) Relative frequency of different P.1 lineage variants among SARS-CoV-2 positive cases sequenced in the Amazonas. (C) Relative frequency of different P.1 lineage variants among SARS-CoV-2 P.1 Brazilian sequences detected outside the Amazonas states.
FIG 2
FIG 2
Maximum likelihood phylogenetic tree of P.1 Amazonian sequences and P.1+NTDdel, P.1+N679K, and P.1+P681H sequences detected outside Amazonas. Tips were colored according to the S mutations as indicated in the legend. Major P.1 sublineages carrying additional mutations/deletions in the S protein were highlighted with colored boxes. The aLRT support values are indicated in key branches, and branch lengths are drawn to scale with the bar indicating nucleotide substitutions per site.
FIG 3
FIG 3
Temporal structure and phylogeographic reconstruction of the P.1+NTDdel, P.1+N679K, and P.1+P681H clades. (A) Root-to-tip regression of genetic divergence against dates of sample collection. P.1 sequences were colored green, while each P.1 subclade carrying deletions or additional mutations in S protein was colored following the legend. Time-resolved maximum clade credibility phylogenies of each P.1 subclade defined in the ML analysis: (B) P.1.4; (C) P.1.6; (D) P.1.7; (E) minor clades P.1.3, P.1.5, P.1.8 and P.1+Δ141-144AM. Tips and branches colors indicate the sampling state and the most probable inferred state of the nodes, respectively, as indicated in the legend for each tree. Bayesian posterior probabilities are indicated in key branches. All horizontal branch lengths are time-scaled, and the tree was automatically rooted under the assumption of the molecular clock model.
FIG 4
FIG 4
Epidemic trajectories of major SARS-CoV-2 lineages circulating in Amazonas in 2021. Graphs depicting the temporal variation in Re (median and 95% HPD) of Amazonian lineages P.1, P.1.4 and P.1.6 were estimated using the BDSKY approach.
FIG 5
FIG 5
Estimation by RT–PCR of viral load in the upper respiratory tract of SARS-CoV-2 infected patients in Amazonas. (A) Graph depicting the relative prevalence of P.1+ lineages estimated from whole-genome sequencing (dashed gray line) and the Ct (mean and 95% Confidence Interval) among SARS-CoV-2 positive cases (solid line) in Amazonas between March and July 2021. (B) Comparison of Ct values from March to July 2021. Horizontal bars represent Ct medians and IQR. Two-sided P values for the nonparametric Mann–Whitney test are shown for each group. Two-sided P-values <0.05 were considered statistically significant.
FIG 6
FIG 6
Binding between furin enzyme and the structural motif 679NSPRRARS686 of the SARS-CoV-2 Spike protein. (A) Representation of furin as electrostatic potential surface, showing the negative charge distribution around the loop substrate (represented in licorice cyan). Arrows indicate the native position of variant mutations N679K and P681H/R. (B) Free energy surface landscape of dissociation of modeled loops from furin enzyme depicted through the potential of mean force (PMF) as a function of the chosen CV. Dashed lines represent the lowest energy basin for the dissociation of the peptides.
FIG 7
FIG 7
Relative prevalence of P.1 and P.1+ lineages among unvaccinated and fully vaccinated groups. All breakthrough cases were detected among individuals fully vaccinated with CoronaVac. The P-value for the Chi–Square test is shown. P-values <0.05 were considered statistically significant.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Fundação de Vigilância em Saúde do Amazonas. 2021. Boletim diaria COVID-19 no Amazonas 15–06-2021. https://www.fvs.am.gov.br/noticias_view/6264. Accessed July 03, 2021.
    1. Naveca FG, Nascimento V, Souza VCd, Corado AL, Nascimento F, Silva G, Costa Á, Duarte D, Pessoa K, Mejía M, Brandão MJ, Jesus M, Gonçalves L, Costa CFd, Sampaio V, Barros D, Silva M, Mattos T, Pontes G, Abdalla L, Santos JH, Arantes I, Dezordi FZ, Siqueira MM, Wallau GL, Resende PC, Delatorre E, Gräf T, Bello G. 2021. COVID-19 in Amazonas, Brazil, was driven by the persistence of endemic lineages and P.1 emergence. Nat Med 27:1230–1238. doi:10.1038/s41591-021-01378-7. - DOI - PubMed
    1. Faria NR, Mellan TA, Whittaker C, Claro IM, Candido DDS, Mishra S, Crispim MAE, Sales FCS, Hawryluk I, McCrone JT, Hulswit RJG, Franco LAM, Ramundo MS, de Jesus JG, Andrade PS, Coletti TM, Ferreira GM, Silva CAM, Manuli ER, Pereira RHM, Peixoto PS, Kraemer MUG, Gaburo N, Jr, Camilo CDC, Hoeltgebaum H, Souza WM, Rocha EC, de Souza LM, de Pinho MC, Araujo LJT, Malta FSV, de Lima AB, Silva JDP, Zauli DAG, Ferreira ACS, Schnekenberg RP, Laydon DJ, Walker PGT, Schluter HM, Dos Santos ALP, Vidal MS, Del Caro VS, Filho RMF, Dos Santos HM, Aguiar RS, Proenca-Modena JL, Nelson B, Hay JA, Monod M, Miscouridou X, et al. . 2021. Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil. Science 372:815–821. doi:10.1126/science.abh2644. - DOI - PMC - PubMed
    1. Fujino T, Nomoto H, Kutsuna S, Ujiie M, Suzuki T, Sato R, Fujimoto T, Kuroda M, Wakita T, Ohmagari N. 2021. Novel SARS-CoV-2 variant identified in travelers from Brazil to Japan. Emerg Infect Dis 27:1243–1245. - PMC - PubMed
    1. Governo do Estado do Amazonas, Fundação de Vigilância em Saúde do Amazonas. 2021. Vacinômetro - COVID-19. https://www.fvs.am.gov.br/indicadorSalaSituacao_view/75/2. Accessed July 03, 2021.

Publication types

MeSH terms