We know the city where HIV first emerged
bbc.com
bbc.com
Here is the first report of simian AIDS: "Examination of the species-specific annual mortality rates of macaques at the center during the previous 4 yr showed a significant increase in deaths in 1980 and 1981" http://www.ncbi.nlm.nih.gov/pmc/articles/PMC393899/
Also, it is convenient they do not share the mcmc traces. Based on these settings I bet they were ugly:
>"For each data set, at least 3 MCMC chains of 250 million steps were computed. Parameters and trees were sampled every 50,000th step. Samples were combined with LogCombiner (77) and between 10 to 30% of each MCMC chain was discarded as burn-in. MCMC mixing was diagnosed using visual trace inspection and calculation of effective sample sizes in Tracer (77). We report the posterior mean and 95% Bayesian credible intervals for evolutionary parameters."
http://dx.doi.org/10.1126/science.1256739
Another thing, they write:
"Our estimated location of pandemic origin explains the observation that Kinshasa exhibits more contemporary HIV-1 genetic diversity than anywhere else"
This is totally circular. As shown in their figure S2, they used genetic diversity in a region as an indication of earlier presence.
I'm not sure we (that is, humans) realized that SIV was a thing until HIV prompted us to look in those places. It's the classic example of not knowing what to look for until you've found it.
>"Our estimated location of pandemic origin explains the observation that Kinshasa exhibits more contemporary HIV-1 genetic diversity than anywhere else" >This is totally circular. As shown in their figure S2, they used genetic diversity in a region as an indication of earlier presence.
It seems like they also used phylogenetic history. Diversity alone is not a sufficient to identify origination.
>A very high genetic diversity of HIV-1 has been reported, not only in Kinshasa and the north and south of the DRC (12, 13, 31, 32), but also in Brazzaville in the RC and, to a lesser extent, in the Mayombe area of RC near Pointe-Noire, all of which have been suggested as potential source locations of the pandemic (22, 33, 34). We therefore performed phylogeographic analyses of viruses collected in both the DRC and RC (table S1) and compared sequence sampling locations with phylogenetic history to formally test hypotheses concerning the location of ancestral viral lineages (30). Our analyses robustly place the spatial origin of the HIV-1 group M pandemic in Kinshasa [posterior probability (PP) = 0.99]
They could have used "consistent with" rather than "explains", though that may be a bit pedantic and debatable and we wouldn't be in the sorry state of science education/trust if all researchers would communicate like Neil deGrasse Tyson.
Perhaps I am misunderstanding their model, but look at figure S2. It shows the less closely related the sequences at each location, the more likely that was the originating region relative to the others. See how the leafs of the right tree corresponding to location B are more widely distributed (indicating diversity) than those for location A or C (which are more clustered within the tree)?
Obviously, just as humans are higher than monkeys on the Great Chain of Being, so must HIV be higher than SIV. ;)
See http://news.nationalgeographic.com/news/2006/05/060525-aids-... http://www.ncbi.nlm.nih.gov/m/pubmed/8026477/#
There is a video documentary that I watched in the last couple years that I was really hoping to find again and link to.
>"Phylogenetic analysis Phylogenetic relationships were estimated from comparisons of predicted protein sequences. Sequences were aligned using CLUSTAL (Higgins and Sharp, 1988, 1989). Evolutionary distances between all pairs of sequences were computed using Kimura's empirical method (Kimura, 1983; eqn. 4.8) to estimate the number of superimposed amino acid replacements; sites at which there was a gap in any sequence in the alignment were excluded from all comparisons. Phylogenetic relationships were estimated from these distances by the neighbor-joining method (Saitou and Nei, 1987). The reliability of branching orders was estimated by the bootstrap approach (Felsenstein, 1985). These methods were implemented using CLUSTAL V (Higgins et al., 1992).
Nucleotide accession numbers All sequences were submitted to GenBank and are available under accession numbers U03994-U04018."
It looks like those sequences are a mixture of the various proteins (env, pol, etc) so it would take some time to figure out which is which. I don't feel it like right now. It seems like getting the sequences from genbank and then aligning using default params with clustalw would be closest to their method. If someone wants to do it: http://www.ncbi.nlm.nih.gov/genbank/ http://www.genome.jp/tools/clustalw/
>"Barring the possibility that humans introduced multiple species-specific SIV lineages to the wild monkey populations of Bioko, the mainland and island SIVdrl variants must have been evolving independently since Bioko became isolated ~10,000 yr B.P., and perhaps longer given the high levels of genetic diversity seen within local SIV populations."
http://www.sciencemag.org/cgi/content/abstract/329/5998/1487
Interesting stuff. I hadn't thought of this until now, but from these few papers I don't think anyone has tried very hard to rule out the human to simian transmission idea.
That is them ruling out that it came from humans.
Several of the authors are also pretty hardcore Bayesian methods people in the field of phylogenetics specifically applied to the evolution of diseases. It's unlikely that their long MCMCs were due to some kind of coverup; it's possible that either the large number of sequences made convergence more difficult and / or they wanted to be absolutely certain of their results by running their analyses longer than usual. These sorts of phylogeographic methods (especially with an asymmetric movement model and lots of different localities) tend to have likelihood surfaces with many ridges, making MCMC quite difficult for all but the most trivial datasets.
They thinned 250 million steps to 500, then it sounds like they dropped the first 50-100 on top of that post-hoc (the 10-30% burn in gives an impression of manual tinkering).
Obviously this is going to bring up questions regarding convergence, so reviewers should have asked for the diagnostic charts. This would probably have been better suited to ABC:
https://en.wikipedia.org/wiki/Approximate_Bayesian_computati...
The standard number of samples for phylogenetic analyses is 1,000-10,000 samples, trending towards the smaller as the phylogenies get larger. Since they did 3 independent analyses and combined their samples, as long as they assessed convergence and ensured that all 3 chains were converged this is likely fine. Excluding burnin is also quite typical since the MCMC move operators for phylogenetics are not very good due to the high dimensionality of the search space. ABC is also challenging to work with due to the high dimensionality (though I admit I haven't worked much with these methods in phylogenetics)
Typically in this field it is not enough to say "you didn't run your MCMC chains for long enough" or anything like that. Of course the chains should be run for longer -- theoretically speaking they should be run for infinite amount of time, but there are routine disagreements over whether 100 or 1000 effective samples are sufficient for phylogenetics. But unless there is a serious model inadequacy that the authors haven't addressed there's typically no reason to nitpick about these sorts of things. The authors have covered a lot of their bases by trying a number of different models and priors and I don't see any reason to doubt this particular study based on that.
Would they be able to run these MCMC chains in a way that perhaps progressively renders results at greater and greater accuracy? Then they could get some initial "low resolution" results and they wouldn't necessarily have to wait for weeks, but the full data will eventually be available.
As far as any circular reasoning, they're just reiterating the claim in a different way. It's saying "Based on X data, we deduce Y. Y would explain why X happened."
Direct download if you don't want the Flash audio player: http://50.31.154.42/radiolab/radiolab111411.mp3?downloadId=5...
If anyone else is curious like I was.
After reading it I think that the exact city is the least interesting fact they have to offer.
Direct link to Jacques Pépin's "The Origin of AIDS" on amazon: http://www.amazon.com/dp/0521186374
Are they saying that the jump happened in a large metropolis ?
Or that Kinshasa was the first large city in which HIV became an epidemic ?
In many ways, whether you support e.g. more direct foreign investment areas in India and China or not, I'd prefer the conversations about costs to be as well formulated as this (I also suspect free trade and foreign investment in China and India has done a lot more good than HIV has done bad, even if you just count lives saved).
whuhh???
Just 'caus it had a CIA sponsored coup, doesn't mean things got any better.