Using nearly full-genome HIV sequence data improves phylogeny reconstruction in a simulated epidemic
Yebra G., Hodcroft EB., Ragonnet-Cronin ML., Pillay D., Brown AJL., Fraser C., Kellam P., de Oliveira T., Dennis A., Hoppe A., Kityo C., Frampton D., Ssemwanga D., Tanser F., Keshani J., Lingappa J., Herbeck J., Wawer M., Essex M., Cohen MS., Paton N., Ratmann O., Kaleebu P., Hayes R., Fidler S., Quinn T., Novitsky V., Haywards A., Nastouli E., Morris S., Clark D., Kozlakidis Z.
AbstractHIV molecular epidemiology studies analyse viral pol gene sequences due to their availability, but whole genome sequencing allows to use other genes. We aimed to determine what gene(s) provide(s) the best approximation to the real phylogeny by analysing a simulated epidemic (created as part of the PANGEA_HIV project) with a known transmission tree. We sub-sampled a simulated dataset of 4662 sequences into different combinations of genes (gag-pol-env, gag-pol, gag, pol, env and partial pol) and sampling depths (100%, 60%, 20% and 5%), generating 100 replicates for each case. We built maximum-likelihood trees for each combination using RAxML (GTR + Γ), and compared their topologies to the corresponding true tree’s using CompareTree. The accuracy of the trees was significantly proportional to the length of the sequences used, with the gag-pol-env datasets showing the best performance and gag and partial pol sequences showing the worst. The lowest sampling depths (20% and 5%) greatly reduced the accuracy of tree reconstruction and showed high variability among replicates, especially when using the shortest gene datasets. In conclusion, using longer sequences derived from nearly whole genomes will improve the reliability of phylogenetic reconstruction. With low sample coverage, results can be highly variable, particularly when based on short sequences.