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

Samples Unassigned #451

Closed
huipingis opened this issue May 17, 2022 · 10 comments
Closed

Samples Unassigned #451

huipingis opened this issue May 17, 2022 · 10 comments

Comments

@huipingis
Copy link

Dear Sir/Madam,

Recently we got results on COVID genotyping, and found some samples unassigned, see below:

NEG6_13052022 | Unassigned | 0,090909091 | PUSHER-v1.8v4.0.6 | pass | Ambiguous_content:0.15 | Omicron (Unassigned) | Usher placements: BA.2(10/11) BA.2.11(1/11); scorpio replaced lineage inference BA.2

POS1_13052022 | Unassigned | 0 | PUSHER-v1.8v4.0.6 | pass | Ambiguous_content:0.12 | Probable Omicron (Unassigned) | Usher placements: BA.2.3.3(2/2); scorpio replaced lineage inference BA.2.3.3

pOS2_13052022 | Unassigned | 0,333333333 | PUSHER-v1.8v4.0.6 | pass | Ambiguous_content:0.1 | Probable Omicron (Unassigned) | Usher placements: BA.2(1/3) BA.2.3.3(2/3); scorpio replaced lineage inference BA.2.3.3

I am wondering why these samples were unassigned for lineages? Thanks.

Best regards,
Chen Huiping
University Hospital of Iceland

@Hideaki615
Copy link

Hideaki615 commented May 18, 2022

I've got the same issue when I updated my local pangolin_3.1.20 to 4.0.6.

Then it could be avoided by downgrading constellations-0.1.10 to 0.1.9, by executing a command like below.

$ conda activate pangolin
$ pip install git+https://github.com/cov-lineages/[email protected]

Adding, I don't know how much it would be related, but I found the output of scorpio call is decreasing in the latest version, like below.

sample-01
-- constellations-0.1.9
---- scorpio call: Alt alleles 37; Ref alleles 1; Amb alleles 20; Oth alleles 1
-- constellations-0.1.10
---- scorpio call: Alt alleles 21; Ref alleles 0; Amb alleles 12; Oth alleles 0

So the latest one may have underestimated alleles and therefore made the sample "Unassigned", I guess...

@huipingis
Copy link
Author

Thanks. The 3 samples I showed: 1 negative control and 2 positive controls. I speculate that they were unassigned due to contamination during PCR and sequencing. It´s clear the negative control turned positive because of contamination. The other 2 positive controls (lineage is B) were unassigned probably due to contamination too, because I saw some reads from other strains by alignment of reads. I wonder if the contamination is one of reasons for lineage un-assignment? Thanks.

Best regards,
Chen

@huipingis
Copy link
Author

huipingis commented May 18, 2022 via email

@aineniamh
Copy link
Member

Unassigned is often a result of conflicting signals or reference calls that scorpio thinks would make a lineage call unreliable. If you suspect your sequences are from contaminated samples, this would exactly be why scorpio would mask the lineage call.

@huipingis
Copy link
Author

Thanks, Áine.

I wonder if the contamination can cause a wrong call, e.g., from B to BA.3 by using pangolin?

Regards,
Chen

@huipingis
Copy link
Author

huipingis commented May 18, 2022 via email

@aineniamh
Copy link
Member

Hi @huipingis, this can absolutely happen as conflicting signals can confuse the inference. This is why we curate the assignments and check this, reporting Unassigned if the sequence doesn't meet the scorpio checks.

@eddieimada
Copy link

Hi @aineniamh,
Does this means that if scorpio cannot confidently assign a lineage, the lineage call from UShER should not be trusted and discarded?

@huipingis
Copy link
Author

huipingis commented May 19, 2022 via email

@aineniamh
Copy link
Member

Hi all, in the latest release of pangolin (4.1), we now have introduced an updated process.

We are confident in the calls by usher, so now bypass scorpio in usher mode and do not overwrite the usher assignment (so now default behaviour is to take the usher assignment as is, but still report the scorpio call). In pangoLEARN/ fast mode we still overwrite assignments that do not explicitly agree with scorpio to elimitate false positive calls with the pangoLEARN method.

Hopefully this will resolve these Unassigned samples, but bear in mind that they would have been unassigned because they did not pass the scorpio checks.

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

No branches or pull requests

4 participants