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

Dating divergent sequences/tips #28

Open
ekankaka opened this issue May 22, 2024 · 1 comment
Open

Dating divergent sequences/tips #28

ekankaka opened this issue May 22, 2024 · 1 comment

Comments

@ekankaka
Copy link

For HIV infection, diversity generally increases with time. But based on the viral load trajectory in natural HIV infection, early viremic infection also seems to have a lot of evolution (which reduces and increases again in late chronic infection). If true, divergent sequences with censored/missing dates would have equal chance of mapping to early or late viremic infection. How does bayroot handle this?

@ArtPoon
Copy link
Contributor

ArtPoon commented May 27, 2024

The short answer is that bayroot does not handle variation in the rate of diversification at all.

You're making the assumption that rate of diversification is roughly proportional to viral load. This is not necessarily a bad assumption but it is complicated. For example, neutral evolution in terms of substitution (divergence, not diversification) should be independent of population size (mutation input is $N \times \mu$, probability of fixation is $1/N$, so substitution rate becomes $u$).

It is conceivable that there is a burst of diversification at an early stage of infection associated with viremia. However, this is exponential growth from a fairly homogeneous population. Have we ever seen a surge of diversity associated with viremia in HIV RNA? Off the top of my head, I can't think of an example of this, although I assume these samples would be quite difficult to obtain.

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

2 participants