German BA.2/BA.1/BA.2 recombinant with S:K147E, S:R346K, S:460K, S:493 reversion, S:H1101Y [55 seq as of 2022-07-22] #823
Lähete asiasta #823: Fedeueli spotted this one and mentioned it in #773
Description:
Recombinant between: BA.2/BA.1/BA.2
Earliest sequence: 1st of June 2022
Most recent sequence: 23th of June 2022
Country: Germany
Likely breakpoint: between 13195 and 15240 and between 21618 and 21762 [S:27 to S:67]
Private Mutations: C4927T, T5386G (like BA.1), C7834T, C12049T,
A220001G/S:K147E, G22599A/S:R346K, G22992/S:S477N, C24863T/S:H1101Y,
A27507G .
(20 days ago) "Thanks for analysing these sequences! great job.
edited:
i would add that this sublineage has the S:H69- S:V70- 69/70del
and the S:493Q reverted to wildtype,
both important features.
I am still not completely convinced that it is a simple recombinant
event or something popping up from or next to the "Omicrons" source.
silcn: " apart from the fact this has only been found in Germany, another argument against a link to the Omicron source is the presence of 9866T in the BA.2 section. A combined BA.1/BA.2 chronic infection could certainly produce something like this".
" Good point thx, silcn".
JS, from BA.1 there is also the 5386G that was highlighted by C198.. in the other issue".
(19 days ago) Shay: "In Israel we got a sample (will soon be uploaded to GISAID) of a sample clustering with this branch of a passenger coming from Greece (no contact to Germany)".
"If this is a real recombinant, it means that it's a recombinant +
saltation. Because then we consider the 69/70 del not as part of BA.1
section (As the C21762T mutation of BA.1 is missing). And also it has
many nonBA.2/BA.1 mutations which are highly converged in Omicrons 2nd
gen/stepwise :
S:K147E
S:R346K
S:N460K
and the highly converged reversion in 493".
FedeGueli:
"Thanks shay ,
to me it is still contested if it is a plain recombinant , a
recombinant in a chronic case or something coming from the entourage of
original Omicron source.
The silcn
argument on Germany is weakened by the fact that BA.5+S:346I was barely
sequenced in South Africa and then popped up and emerged fastly in
Germany, representing a brand new branch of BA.5 distinct from the main
designated or even proposed ones
"More:
And i was just reading casually this paper yesterday" :
"Potential evidence of a single recombination event
involving BA.1, BA.2 and BA3 was identified by 3SEQ (P = 0.03), GARD
(delta c-AIC = 20) and RDP5 (GENECONV P = 0.027; RDP P = 0.006) within
the NTD encoding region of spike. The most likely breakpoint locations
for this recombination event were 21690 for the 5′ breakpoint (high
likelihood interval between 15716 and 21761) and 22198 for the 3′
breakpoint (high likelihood interval between 22197 and 22774)." (Tulio
de Oliveira et Al:
https://www.nature.com/articles/s41586-022-04411-y)
the original BA.1/BA.2/BA.3 breakpoint seems to overlap with the
breakpoints of this one described here. Sincerely i dunno know if
significant or not, very far from me i would say to say anything
relevant, but maybe someone else could take a 2nd look at this".
Silcn: "BA.5+S:346I doesn't come directly from the Omicron source though. There have been at least 5 (maybe more, see #455, #561) separate emergences from the Omicron source and in every case most of the early sequences came from South Africa".
FedeGueli: "Yes true i mean if a very transmissible lineage is barely detected in the country where it likely emerged i think a much lower transmissible one can fly under the radar easily. Not saying it is wrong what have you said but that is not the strongest argument vs BA.6 hypothesis, i think. Also as you i think a new recomb in a chronic host different from Omicron source is the most likely explanation".
Josette: "Thanks for all the comments! It gave me some extra insight.
I think it is more complicated than a normal recombinant with 2 breaking
points. I missed the 69/70 del, since it is not in Sc2rf, but that
would make it 2 BA.1 mutation clusters instead of one.
Together with the T5386G BA.1 mutation, which is shown by Sc2rf, and I
did notice, but didn't highlight as BA.1 (just edited it in the first
part of the issue) and the 493 reversion (A23040A) that really makes it
very complicated.
To summarize the nucleotides:
1-4321 BA.2
T5386G BA.1
8393-13195 BA.2
15240-21618 BA.1
C21762C BA.2
21765-21770 BA.1 (69/70del)
21846-22882 BA.2
A23040A S:493Q reversion, no BA.1/2
23048-29510 BA.2
Plus of course several private mutations."
Josette: " changed the title
German BA.2/BA.1/BA.2 recombinant with extra S:K147E, S:R346K, S:S477N, S:H1101Y and 5 other private mutations
German BA.2/BA.1/BA.2 recombinant with extra S:K147E, S:R346K, S:H1101Y and 5 other private mutations
19 days ago
Silcn: " also possible that 69/70del occurred independently rather than being BA.1-derived. Wouldn't be too surprising given the other additional spike mutations".
(18 days ago)
FedeGueli: "Thanks, Josette! I suggest to edit your main comment so it would be clearer for everyone !"
FedeGueli: " Silcn, that is possible, but in a landscape of multiple breakpoints and a key reversion for other Omicron lineages i dont know how much we can be sure of that or not".
Silcn: "Josette, also interesting that you mention BA.3, because something very similar happened there: it contains S:G446S from BA.1 in an otherwise ancestral BA.2-derived section. Again, impossible to know if it arose independently or though recombination. Still doesn't change my opinion that this lineage probably has nothing to do with the Omicron source (the C9866T is convincing enough for me). The S:R493Q reversion doesn't make things any more complicated than they already are: it is extremely common in saltation-derived BA.2 sequences/lineages with multiple Spike mutations."
Shay: " As I mentioned yesterday, a first Israeli sample (traveler from Greece) :
(15 days ago)
Ulrich: "We have two cases in Austria as well, currently uploading to Github".
FedeGueli: " Ulrich, thx a lot! To the highest Authorities: C.R. C.R A.H. IPG : in my view this has to be granted an early designation, whatever as recombinant or something else, still very different from everything else circulating (or better saying that it has a unique mix of mutations actually on the rise) .
Ulrich: "What I find remarkable about this variant is that it appears to be convergent evolution with BA.2.75 (147, 460, 493) and with regards to the 493 reversion even BA.5. It seems obvious the virus has identified the next surfaces to optimise."
Shay : "Ulrich This is why it might be a combination of recombination+saltation in a chronic patient. The mutations converged in this one are indeed more related to variants that arose outside SA and the parallel BA.1-5. And the recombination is much cleaner than could be seen in cases where presumably a wide range of variants been co-evolved like might be in the original chronic patient ( like BA.3 probably coming from the same body where BA.1 and BA.2 evolved and XB probably coming from the same body were B.1.631 and B.1.634 evolved)."
"If anyone needs -here is a comparison of the mutational profile of this variant , BA.1 and BA.2
PDI 823 .xlsx "
FedeGueli: " Thx, Shay!"
Shay: "For CovSpectrum query i would replace to :
[4-of:C4927T,C7834T,C12049T,A22001G,G22599A,T22942G,C24863T,A27507G]
searching 4 out of the 8 unique mutations in this lineage (not present in BA.1 or BA.2)."
Josette: Here's a link for Shays 's query. ( Kts. linkistä)
(11 days ago)
Josette: "We found another one. Also I tried to illustrate a bit what the mutations are like":
(9 days ago)
Fede Gueli: " 23 sequences as today
Early signs of growth advantage vs BA.5 limited to Germany".
(6 days ago)
Shay: "27 samples as of yesterday".
Ulrich: "With the search term Spike_K147E, Spike_N460K, Spike_G339D, Spike_H69del currently a total of 60 cases are found. We mostly sequence only spike and selected regions in Austria, so these samples are only found with a term that does not demand too many mutations beyond spike. Could we maybe adapt the search term to account for that? 2 samples were also full-genome sequenced already. Really looking like growth advantage. Plus the spread across Germany and beyond says it all...https://user-images.githubusercontent.com/109073004/180005968-d957dfe5-c284-47fc-aa58-67036b0ec2a8.png
FedeGueli: "Ulrich, i would add Spike_H1101Y to your query".
C.Roemer. "FedeGueli, I would recommend using advanced query with NextcladePangoLineage:BA.5* instead of just two mutations which do not immediately show which lineage you compare to. I've edited your growth advantage query
Josedtte, I think your sc2rf plot is confusing as it has Delta BA.2 as parents rather than BA.1 and BA.2. Can you explain why you did that and maybe change it to BA.1/BA.2 if this is a mistake? It's confusing because the title of the issue suggests this is a BA.2/BA.1/BA.2 recombinant.
Here's a quick summary of how I think about this recombinant
Spike is mostly BA.2* but with some notable changes:
- Missing: T19I, 24-26del, 27S not there because that part is from BA.1* parent
- Donated: 69-70del convergent with many VOCs
- Additional S:R346K, convergent with BA.1.1
- Additional S:K147E, S:N460K and reverted S:493 convergent with BA.2.75
I've added a BA.1/BA.2 sc2rf pic to the first post in this issue to keep things simple for those reading from top to bottom, Josette
I tried to find which branch the BA.2* donor is from but couldn't find anything interesting in Germany/Europe. For the BA.1 part, there are no private mutations so if anything one can go by exclusion principle.
corneliusroemer
changed the title
German BA.2/BA.1/BA.2 recombinant with extra S:K147E, S:R346K, S:H1101Y and 5 other private mutations
German BA.2/BA.1/BA.2 recombinant with extra S:K147E, S:R346K, S:460K, S:493 reversion, S:H1101Y and 4 other private mutations
6 days ago
Josette Shoenma: "Thanks, Cornelius Roemer . Somehow I forgot to include it. Couldn't have done my analysis without it. Here is the one I made from the first 16 sequences. It won't be very different..."
corneliusroemer
changed the title
German BA.2/BA.1/BA.2 recombinant with extra S:K147E, S:R346K, S:460K, S:493 reversion, S:H1101Y and 4 other private mutations
German BA.2/BA.1/BA.2 recombinant with S:K147E, S:R346K, S:460K, S:493 reversion, S:H1101Y [55 seq as of 2022-07-22]
5 days ago
Corneliusroemer:
@chrisruis @InfrPopGen @thomasppeacock
There are now at least 55 sequences of this recombinant, 37 have been submitted since Tuesday! I would strongly recommend quick designation as the criteria are satisified, the Spike profile is very interesting and there appears to be growth, too.
Here are recent EPI_ISLs of complete sequences (3 are Spike only) so you can designate even if these sequences aren't yet in the Usher tree (maybe excluded due to long branch? @AngieHinrichs
|
Thanks @corneliusroemer
-- yep, they are filtered out due to having too many reversions
relative to BA.2 placement by nextclade, but I've been manually
exempting the EPI_ISL's from Josette's initial list and EPI_ISL_13719505
from the cov-spectrum query as of a few days ago. They're in the tree,
on a branch with the label proposed823 (Pango lineage assigned by
UShER). .... I will exempt the rest from your list, so they should be added in tomorrow's build. CorneliusRoemer: " Thanks AngieHinrichs , My list was created using manual AA substitution queries in GISAID - not super reliable. So not surprised you catch a few more."
InfrPopGen commented 4 days ago
|
Inga kommentarer:
Skicka en kommentar