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

Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (277 seqs) #1814

Closed
FedeGueli opened this issue Mar 28, 2023 · 13 comments
Closed
Assignees
Milestone

Comments

@FedeGueli
Copy link
Contributor

I found a sublineage of BN.1.2.3 that likely emerged in Taiwan reaching there an high prevalence in February .
It was sampled both in Taiwan homeland and in travellers from Taiwan to Singapore,Usa, South Korea and Myanmar.
Then it started to circulate when it was sample little less than 100x times and showed there a growth comparable to XBB.1.5 (that doesnt seem to have started a wave there)

Gisaid query: C1471T,C1513T,C25572T, C29370T

Nr of total sequences: 188 (99 Japan)
Taiwanese sequences (27)
IDs
EPI_ISL_16259451, EPI_ISL_16537531, EPI_ISL_16676294,
EPI_ISL_16828717-16828720, EPI_ISL_16853149, EPI_ISL_17016390,
EPI_ISL_17040800, EPI_ISL_17040887, EPI_ISL_17053655,
EPI_ISL_17075667-17075668, EPI_ISL_17089525-17089526, EPI_ISL_17104729-17104730,
EPI_ISL_17104732, EPI_ISL_17104736, EPI_ISL_17172467,
EPI_ISL_17213326, EPI_ISL_17256547-17256548, EPI_ISL_17280157,
EPI_ISL_17289382, EPI_ISL_17289384

Defining mutations:
BN.1.2.3 (C1471T) >> S:K187E (A22121G) > C1513T >C25572T, N:T366I (C29370T)

Tree:
Schermata 2023-03-28 alle 13 59 18
https://nextstrain.org/fetch/genome.ucsc.edu/trash/ct/subtreeAuspice1_genome_1fb0c_2cf250.json?c=country&label=id:node_8037222

Its gain of fitness seems clear
Schermata 2023-03-28 alle 14 01 14
It looked competitive with XBB.1.5 although this plot is highly biased by its peculiar geographical spread
Schermata 2023-03-28 alle 14 02 09

its prevalence during 2023 in Taiwan was 27 sequences out of 306 globally foud searching for Taiwan/2023 on Gisaid.
but it was 18/128 in February maybe pointing that it managed to gain shares there.

@FedeGueli FedeGueli changed the title Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (188seqs) Mar 28, 2023
@xz-keg
Copy link
Contributor

xz-keg commented Mar 28, 2023

Samples uploads from Taiwan are mostly imported cases. Exclude them, then you get a much higher ‘local share“ for the variant.

When calculating growth rates shall be careful on this, as the imported/local ratio is not constant.

@FedeGueli
Copy link
Contributor Author

Samples uploads from Taiwan are mostly imported cases. Exclude them, then you get a much higher ‘local share“ for the variant.

def right @aviczhl2

@FedeGueli FedeGueli changed the title Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (188seqs) Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (208 seqs) Mar 31, 2023
@FedeGueli
Copy link
Contributor Author

This grew to 208. Last sequences coming in : from Singapore (2 local + 2 imported from TW) Taiwan and Japan.

Looking at samples from March only it is 4 out of 36 sequences. 2 out of 15 considering just the ones sampled in Taiwan

@FedeGueli FedeGueli changed the title Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (208 seqs) Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (251 seqs) Apr 10, 2023
@FedeGueli
Copy link
Contributor Author

251 sequences as today also sampled in New Zealand. Worth designating due its local relevance and international spread notwithstanding the fact is surely less transmissible than 486P carrying lineages

@FedeGueli FedeGueli changed the title Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (251 seqs) Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (274 seqs) Apr 16, 2023
@FedeGueli
Copy link
Contributor Author

274 sequences, 37 out of 274 from Taiwan or travellers from Taiwan.

The tree shows a lot of diversity and very little clustering both for countries and for timing ( i mean last sequences are all scattered in the tree)
and it is quite strikingly that the biggest branch of this lineage counts only 21 sequences
Screenshot 2023-04-16 alle 08 03 29
Screenshot 2023-04-16 alle 08 03 48
Screenshot 2023-04-16 alle 08 04 38

My guess these data point to a prolonged not very fast local circulation in an undersampled area likely Taiwan

@xz-keg
Copy link
Contributor

xz-keg commented Apr 16, 2023

Do you find any GenBase seqs on this tree? They are not on GISAID but they've been added to usher recently.

@FedeGueli FedeGueli changed the title Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (274 seqs) Taiwanese BN.1.2.3 sublineage defined by S:K187E and N:T366I (277 seqs) Apr 17, 2023
@InfrPopGen InfrPopGen self-assigned this Apr 21, 2023
InfrPopGen added a commit to InfrPopGen/pango-designation that referenced this issue Apr 21, 2023
…esignations, and 4 updated from BN.1.2.3
InfrPopGen added a commit that referenced this issue Apr 21, 2023
Added new lineage FR.1 from #1814 with 244 new sequence designations, and 4 updated from BN.1.2.3
@InfrPopGen InfrPopGen added this to the FR.1 milestone Apr 21, 2023
@InfrPopGen
Copy link
Contributor

Thanks for submitting. We've added lineage FR.1 with 244 newly designated sequences, and 4 updated. Defining mutations C29370T (N:T366I), C25572T (following C1513T).

@ryhisner
Copy link

ryhisner commented Jun 1, 2023

Has anyone noticed that this lineage has a frameshift in ORF7b? GISAID, Nextclade, and BLAST all call it a one-nucleotide insertion: ins27791T. This corresponds to the third nucleotide in ORF7b:12. (Nextclade interprets a few sequences as having 2, 3, 4, 5, or 6 inserted T's at the same position These don't seem especially closely related though, so I suspect those are misreads, probably due to the fact that five consecutive T's normally follow 27791.)

Notably, a stop codon doesn't arrive until 27851-27853, which, after the 1-nucleotide insertion, corresponds to ORF7b:33. So ORF7b's AA composition is changed dramatically, but the length of ORF7b is reduced from 43 AA to 32 AA, about a 25% reduction.

I believe the change in the AA composition of ORF7b should be as pictured below. (I did this manually, so there could be errors.) My first impression is that this insertion results in quite a few more polar residues than previously.
image

@FedeGueli
Copy link
Contributor Author

Great catch Ryan! @ryhisner

@ryhisner
Copy link

ryhisner commented Jun 1, 2023

Actually, I think the new ORF7b this creates is much more widespread than just FR.1. I've noticed a recent trend for frameshift-causing, two-nucleotide deletions in this same area of ORF7b, but I didn't look into them very closely. I believe they should create the same modified, somewhat truncated version of ORF7b as the one-nucleotide insertion in FR.1 does.
image

XBC.1* lineages have one of these 2-nuc deletions, for example (∆27792-27793). However, they have also independently occurred (in small numbers) in a number of different lineages and seem to be very much on the increase of late.
image

@ryhisner
Copy link

ryhisner commented Jun 2, 2023

I noticed that there's a synonymous ORF7b nuc mutation in all Omicron (C28007T) that makes a slight change in the sequence of the new version of ORF7b this creates. It's a small change: alanine to valine in the sixth residue after the insertion. XBC* does not have this mutation, so its sequence remains the same. I've put the sequences next to the old ORF7b sequence, with polar and charges residues in red and nonpolar residues in black.
image

@FedeGueli
Copy link
Contributor Author

Wow @ryhisner great work.
FR.1 is scoring very high in Collection42 by @alurqu , partially due the fact China is uploading fastly and a lot, but still its significance could be beyond the simply acknowledgement of something new

@FedeGueli
Copy link
Contributor Author

Look here what happened in XBC.1.1.2 @ryhisner it got a mutation right before the new stop codon:
Schermata 2023-06-13 alle 19 47 41
https://nextstrain.org/fetch/genome.ucsc.edu/trash/ct/subtreeAuspice1_genome_17adb_88fb00.json?c=userOrOld&label=id:node_6126821

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