I am curious why the tree makes it look like the Thai strains have evolved from particular Wuhan strains when they don't share mutations?
-
-
Als antwoord op @pathogenomenick @trvrb
That will be because they are both from a later time than all the rest. So , all else being equal the coalescent model will put them together because it infers a small population.
1 antwoord 0 retweets 7 vind-ik-leuks -
I think it is premature to be trying to reconstruct time-trees.
1 antwoord 0 retweets 5 vind-ik-leuks -
Als antwoord op @arambaut @pathogenomenick
I wanted mainly to estimate TMRCA. I think the observation of early Dec ancestor fits with Poisson expectation given 5 zeros, 2 ones, 5 twos and 1 three mutation away from common ancestor and a genomic subs rate of 1.1 subs per site per month.
1 antwoord 0 retweets 3 vind-ik-leuks -
Als antwoord op @trvrb @pathogenomenick
What I mean is it is premature to try to estimate a single time tree. The BEAST estimate for TMRCA (depending on the rate you assume) is:pic.twitter.com/ntVKydbxLx
1 antwoord 1 retweet 12 vind-ik-leuks -
Als antwoord op @arambaut @pathogenomenick
Indeed a very good point. But it seems like we're going to need to assume some rate (or range of rates) for a while, while we lack temporal sampling range. I chose 4.59 x 10^-4 owing to https://elifesciences.org/articles/31257 . But a range here would be good.
2 antwoorden 1 retweet 8 vind-ik-leuks -
Right after a host switch the rate might well be different (higher). And the high fraction of non-syn mutations and residual clustering suggests this might be the case.
1 antwoord 0 retweets 3 vind-ik-leuks -
Als antwoord op @richardneher @trvrb en
I'm concerned about ~ 50% of the SNPs in the tree, so I think there's a lot of spurious signal at this stage. Mainly based on the following: 1. Ends that don't match 2. Indels that break ORFs 3. Ts/Tv rate that is too low (too many transversions) 4. Unusual Non-synonymous Tv's
1 antwoord 1 retweet 5 vind-ik-leuks -
Als antwoord op @K_G_Andersen @richardneher en
5. SNPs that are right next to each other. With this in mind, a couple of things: 1. Correct indels in 402120 2. Correct 5' end of 402127 3. After these corrections, I count 14 SNPs across the tree 4. Out of the 14 SNPs, 7 (50%) are highly suspicious
2 antwoorden 0 retweets 2 vind-ik-leuks -
Als antwoord op @K_G_Andersen @trvrb en
Yes. But all of these caveats suggest that the tmrca is very recent. If this virus keeps circulating, real SNPs should soon start dominating. Good thing they sequenced the entire genome. 1000 bases would have left us guessing even more wildly...
2 antwoorden 1 retweet 5 vind-ik-leuks
Agreed. And the genomes themselves are very complete, so really impressive work - we're talking about 1-2 potential sequencing errors throwing off a sequence, which is pretty insane given a 30kb genome!
-
-
Als antwoord op @K_G_Andersen @richardneher en
Based on the current data, I don't really think we can distinguish between H2H and people infected from the same source? We can rule out a widespread reservoir at this stage though - at least for the current cases.
1 antwoord 0 retweets 5 vind-ik-leuks -
Als antwoord op @K_G_Andersen @richardneher en
That would be my reading of it.
1 antwoord 0 retweets 1 vind-ik-leuk - Nog 1 antwoord
Nieuw gesprek -
Het laden lijkt wat langer te duren.
Twitter is mogelijk overbelast of ondervindt een tijdelijke onderbreking. Probeer het opnieuw of bekijk de Twitter-status voor meer informatie.