Comments (19)
Hi @andersonbrito, thanks for letting us know. I've updated the classifications for B.1.526, B.1.526.1 and B.1.526.2 from your genomes, the update should appear in release 1.1.14
from pango-designation.
@andersonbrito Can you send the alignment to me offline?
from pango-designation.
Seems like the Auspice tree may be problematic:
There are some conflicting homoplasies
from pango-designation.
So the two unassigned ones that split these branches are probably the problematic sequences. Also that long one looks odd.
from pango-designation.
FWIW S:L5F (C21575T) is in the Problematic Sites list as highly homoplasic, so we mask it out when building our trees.
Problematic Sites initial report on virological (there have been several updates, but 21575 was there from the beginning): https://virological.org/t/issues-with-sars-cov-2-sequencing-data/473
Problematic Sites VCF file: https://raw.githubusercontent.com/W-L/ProblematicSites_SARS-CoV2/master/problematic_sites_sarsCov2.vcf
(the 7th column contains either 'caution' or 'mask' -- we mask only the positions with 'mask')
from pango-designation.
I think it could be pertinent here because it is associated with the Q1011H and G1946S - I suspect all three need to be pushed back to the common ancestor of B.1.526 and B.1.526.2 with the two unassigned ones (the thin branches) being the problems (presumably they are missing the corresponding pairs). Homoplasic sites can just be neutral sites that are flip-flopping and thus actually provide useful fine scale information (I don't know the pattern of L5F - will look into it).
from pango-designation.
Doing this would make B.1.526.2 less distinguishable from B.1.526 and would probably warrant the withdrawal of B.1.526.2 as a sub lineage. It is possible the S477N should be pushed back to and then either the top clade has a 477S reversion or this part of the tree needs to be re-rooted with the bottome clade as more of an outgroup. This would probable help make more sense for the root-to-tip plot too:
from pango-designation.
Also C21575T (L5F) has occurred within patients suggesting that it is a true homoplasy so likely phylogenetically informative.
from pango-designation.
I think these are the two that are causing the problems:
from pango-designation.
That last one shares many of the same mutations as the stem of the B.1.526.1 lineage:
from pango-designation.
Possibly a recombinant but more likely an artifactual mosaic. @andersonbrito - worth checking these two out for mixtures or contamination issues.
from pango-designation.
from pango-designation.
@rambaut, I have just sent you an email with alignment and tree files.
from pango-designation.
@rambaut: I'm not sure what version of pangolin is being run by GISAID, but here's the latest I see:
on GISAID currently B.1.526.1
is a clear sister lineage to B.1.526
. From https://cov-lineages.org/lineages/lineage_B.1.526.html I believe that B.1.526
is supposed to capture S:D253G and its descendants. However, the lineage for B.1.526.1
is sister to D253G.
Take a look at: https://nextstrain.org/groups/blab/ncov/ny/B.1.526?c=gt-S_253
I have noticed that the homoplasies do complicate structure within B.1.526, but I believe that 253G is consistently partitioning these lineages in question as sister lineages.
from pango-designation.
In this link is the phylogeny with bootstrap values (PDF view here) , mentioned in my previous message.
Purple = B.1.526
Blue = B.1.526.1
Light Blue = B.1.526.2
from pango-designation.
There is something odd with B.1.526 and sublineages: B.1.526 often shows up as a sister clade of B.1.526.1 and B.1.526.2. Additionally, B.1.526.3 is not being assigned. I would expect this structure: (B.1.526(B.1.526.1, B.1.526.2, B.1.526.3)). But what we observe is: (B.1.526, B.1.526.1, B.1.526.2), with B.1.526.3 not being assigned properly.
It seems that what is being called B.1.526 in the link above is actually B.1.526.3. With that, I would expect the designation B.1.526 only assigned to genomes that lack the signatures found in the sublineages .1, .2, and .3. But so far, B.1.526 is presented as a sister group of all the other sublineages, while it should be a parental lineage (see below).
from pango-designation.
Looking at the designated sequences, B.1.526, B.1.526.1 and B.1.526.2 on the surface look very similar from an epi point of view (https://raw.githubusercontent.com/cov-lineages/pango-designation/master/lineages.csv), they've all got sequences from CT little other locations desigated.
Currently the only sequences designated B.1.526.3 are
Luxembourg/LNS9461736/2021,B.1.526.3
Luxembourg/LNS4948053/2021,B.1.526.3
Luxembourg/LNS5785410/2021,B.1.526.3
Belgium/CHUNamur13191626/2021,B.1.526.3
Luxembourg/LNS6053410/2021,B.1.526.3
Belgium/ULG-12466/2021,B.1.526.3
Luxembourg/LNS4693077/2021,B.1.526.3
Belgium/CHUNamur13191713/2021,B.1.526.3
Belgium/CHUNamur13191658/2021,B.1.526.3
Luxembourg/LNS1421682/2021,B.1.526.3
Any of the assignments that are beyond what's in the designation list already can be included if they are >95% complete to help with these assignments.
The question for the .1 and .2 sublineages here may be whether they should be sublineages or just merged back into the parent if it seems they're not robustly being distinguished in a phylogeny.
from pango-designation.
Merged in sublineages to B.1.526 in commit ec2335a.
from pango-designation.
Release tagged: https://github.com/cov-lineages/pango-designation/releases/tag/v1.2.13
from pango-designation.
Related Issues (20)
- JN.1+S:F456L on C23683T branch(132 seqs, 7 countries) HOT 4
- KP.3+N:204P,Orf9b:P3H(117 seqs, 7 countries, 1 with S:A475V) HOT 5
- KP.1.1+S:S31del(46 seqs, 11 countries) HOT 6
- JN.1+S:F456L on C4795T branch (114 seqs, 14 countries) HOT 5
- JN.1 FViRT (346T+456V) with S:G181R, Orf3a:Q213K and PLPro_T75I from Singapore (8 sequences) HOT 1
- JN.1.11.1 + C5467T got FLiRT (27 samples, 8 countries ) HOT 3
- JN.1 + S:F456L + S:V407i (37 seqs, Spain) HOT 13
- JN.1+S:F456L on C6145T branch (79 seqs, 9 countries) HOT 1
- Jn.1.1.5 FViRT lineage with S:L176F (13, mostly US)
- JN.1.50(JN.1 + S:A67V, S:L249F, S:V445P, S:F456L) + S:R346T(6 seq, 3 countries) HOT 7
- JN.1+S:F456L via T22930A (85 seqs, 14 countries, 26 with S:R346T) HOT 5
- Another JN.1 + FLiRT lineage emerged in the Iberian Peninsula (19) HOT 3
- JN.1.1 with a S:R346T sublineage (175 samples) HOT 1
- JN.1.22 with Orf1a:G445V (38 sequences, 2 countries) HOT 1
- JN.1.18.1 + S:F456L (FLiRT) with S:P272H (14 sequences 8+ countries)
- JN.1.16 got S:T572I (23 mostly Spain) HOT 2
- JN.1+S:F456L on C9979T branch (18 seqs, 2 countries)
- JN.1+S:R346T on C2395T branch (35 seqs, 4 countries)
- JN.1.4.5 + S:R346T that further gained S:T732I (52 seqs, 7 countries)
- JN.1.4 with Orf1a:D4117N, S:T22N and S:R346T (37 seqs)
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from pango-designation.