Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Increasing prevalence of BA.5.3.1 sublineage in Amazonas, Brazil, with S:K444T, S:N460K and S:Y144del #1302

Closed
akograf opened this issue Nov 8, 2022 · 5 comments
Assignees
Milestone

Comments

@akograf
Copy link

akograf commented Nov 8, 2022

By Rede Genômica Fiocruz

#Description:
We found a highly supported clade of BA.5.3.1 sequences in Amazonas, Brazil, harboring S:K444T, S:N460K and the deletions S:Y144. Such Amazonian BA.5.3.1 clade is independent of other BA.5.3.1 descendent lineages carrying the same mutations in the Spike (like lineages BE.1.* or BQ.1.*). Genomic surveillance in Amazonas shows that it is increasing fast, from around 25% of sequenced genomes in the first half of September to approximately 94% in the second half of October.

#Non-synonymous mutations outside the Spike:
ORF1a: V84del
ORF1a: M85del

#Synonymous mutations outside the Spike:
A27508C
A27749G
T28693C

Earliest sequence: 2022-09-08
Most recent sequence: 2022-10-20
Countries circulating: Brazil

250 Genomes
Data Availability
GISAID Identifier: EPI_SET_221108bx
doi: 10.55876/gis8.221108bx

Evidence
BA531_tree
Linhagens AM 2022% SET_OUT

@FedeGueli
Copy link
Contributor

It has the same Nuc mutation for S:N460K as BQ* family T22942A
but misses both the defining mutations of BQ outside the spike: Orf9b:P10F and orf1b:ORF1b:Y264H

the two mutations outside the spike are:
ORF7a:T39P A27508C
ORF7a:K119R A27749G
while the other is synonymous : T28693C

first sequence of this lineage should be EPI_ISL_14853495 from Brazil which had S:N460K and not S:444T
collected on 16/07/2022

But before that there were 298 sequence of BA.5.3.1 with A27508C,A27749G,T28693C
but without S:N460K and S:444T .
Although it is more likely convergent evolution especially after seeing this first sequene lacks genuinely the S:K444T mutation, being July the same time when BQ family emerged from Nigeria Ghana region i would not exclude this could have emerged after a recombinant event between the brazilian A27508C,A27749G,T28693C lineage and one virus from the rising african BQ. family.

@akograf
Copy link
Author

akograf commented Nov 9, 2022

Thanks for you comment @FedeGueli

We don’t believe this BA.5.3.1 can be a recombination with BQ.1 because:

  1. There is only one BQ.1 genome detected in Amazonas and it was detected in end of October. Thus, there was very few opportunities of a recombination have occurred.
  2. Sequence EPI_ISL_15603705, from July and located at the base of the clade, has S:N460K+S:Y144del, while BQ.1 doesn’t have Y144del. Recombination should had occurred in a very specific portion of the Spike protein.

Furthermore, we think that a new lineage should be designated only for the sequences harboring S:K444T, S:N460K and S:Y144del. Not including the ancestors with the incomplete set of mutations. This is the variant actually increasing in prevalence and responsible for the resurgence of covid19 in Amazonas.

Finally, we found an issue with the mutations A27508C and A27749G. They only appear in sequences with a 240 bp gap at the positions 27509 to 27748 (see the images below from my BA.5.3.1 alignment). This is observed in sequences from several countries, and I suspect it can be a sequencing issue. Has someone else already observed it?

Screen Shot 2022-11-09 at 11 25 21

Screen Shot 2022-11-09 at 11 28 56

@FedeGueli
Copy link
Contributor

Thx @akograf agree on everything. i also noticed that two nuc mutations arent always shown by Covspectrum as AA mutations , dont know why , as you said some artifact or sequencing issue at work there

@InfrPopGen InfrPopGen self-assigned this Nov 12, 2022
InfrPopGen added a commit that referenced this issue Nov 12, 2022
Added new lineage BE.9 from #1302 with 233 new sequence designations, and 0 updated designations
@InfrPopGen InfrPopGen added this to the BE.9 milestone Nov 12, 2022
@InfrPopGen
Copy link
Contributor

Thanks for submitting. We've added lineage BE.9 with 233 newly designated sequences, and 0 updated designations. Defining mutation A22893C (S:K444T) (following T22942A (S:N460K)).
The current usher tree has A27508C (ORF7a:T39P); then C27509A (ORF7a:P39Q) (this being on same node as nt:A27510G (ORF7a:Y40N), nt:T27511A, and nt:C27513T), in the path from BA.5.3.1 to A22893C (S:K444T).

@JosetteSchoenma
Copy link

CovSpectrum query. I exluded 2 mutations to exclude several BQ.1 and BW.1.
https://cov-spectrum.org/explore/World/AllSamples/Past2M/variants?variantQuery=S%3AK444T+%26+S%3AN460K+%26+S%3AY144-+%26+ORF1a%3AM85-++%26+%21ORF1b%3AT2537I+%26+%21ORF1a%3AL3829F&
GISAID: Spike_K444T, Spike_N460K, Spike_Y144del, NSP1_M85del, -NSP16_T140I, -NSP6_L260F
This query now catches 255, 250 from Brazil, 2 from Denmark, 1 from England, 1 from South Korea and 1 from the USA.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

5 participants