To characterize HIV-1 transmission dynamics in regions where the burden of HIV-1 is greatest, the “Phylogenetics and Networks for Generalised HIV Epidemics in Africa” consortium (PANGEA-HIV) is sequencing full-genome viral isolates from across sub-Saharan Africa. We report the first 3,985 PANGEA-HIV consensus sequences from four cohort sites (Rakai Community Cohort Study, n = 2,833; MRC/UVRI Uganda, n = 701; Mochudi Prevention Project, n = 359; Africa Health Research Institute Resistance Cohort, n = 92). Next-generation sequencing success rates varied: more than 80% of the viral genome from the gag to the nef genes could be determined for all sequences from South Africa, 75% of sequences from Mochudi, 60% of sequences from MRC/UVRI Uganda, and 22% of sequences from Rakai. Partial sequencing failure was primarily associated with low viral load, increased for amplicons closer to the 3′ end of the genome, was not associated with subtype diversity except HIV-1 subtype D, and remained significantly associated with sampling location after controlling for other factors. We assessed the impact of the missing data patterns in PANGEA-HIV sequences on phylogeny reconstruction in simulations. We found a threshold in terms of taxon sampling below which the patchy distribution of missing characters in next-generation sequences (NGS) has an excess negative impact on the accuracy of HIV-1 phylogeny reconstruction, which is attributable to tree reconstruction artifacts that accumulate when branches in viral trees are long. The large number of PANGEA-HIV sequences provides unprecedented opportunities for evaluating HIV-1 transmission dynamics across sub-Saharan Africa and identifying prevention opportunities. Molecular epidemiological analyses of these data must proceed cautiously because sequence sampling remains below the identified threshold and a considerable negative impact of missing characters on phylogeny reconstruction is expected.
Get full access to this article
View all access options for this article.
References
1.
BrennerBG, RogerM, RoutyJP, MoisiD, NtemgwaM, MatteC, et al.: High rates of forward transmission events after acute/early HIV-1 infection. J Infect Dis, 2007; 195:951–959.
2.
OsterAM, DorellCG, MenaLA, ThomasPE, ToledoCA, HeffelfingerJD: HIV risk among young African American men who have sex with men: A case-control study in Mississippi. Am J Public Health, 2011; 101:137–143.
3.
VolzE, IonidesE, Romero-SeversonE, BrandtMG, MokotoffE, KoopmanJ: HIV-1 transmission during early infection in men who have sex with men: A phylodynamic analysis. PLoS Med, 2013; 10:e1001568.
4.
RatmannO, van SighemA, BezemerD, GavryushkinaA, JuurriansS, WensingAM, et al.: Sources of HIV infection among men having sex with men and implications for prevention. Sci Transl Med, 2016; 8:320ra2.
5.
PoonAF, GustafsonR, DalyP, ZerrL, DemlowSE, WongJ, et al.: Near real-time monitoring of HIV transmission hotspots from routine HIV genotyping: An implementation case study. Lancet HIV, 2016; 3:e231–e238.
6.
PillayD, HerbeckJ, CohenMS, de OliveiraT, FraserC, RatmannO, et al.: PANGEA-HIV: Phylogenetics for generalised epidemics in Africa. Lancet Infect Dis, 2015; 15:259–261.
7.
HPTN 071 (PopART) Phylogenetics Protocol Team.. HPTN 071-2 Phylogenetics in HPTN 071: An ancillary study to “Population Effects of Antiretroviral Therapy to Reduce HIV Transmission (PopART): A cluster-randomized trial of the impact of a combination prevention package on population-level HIV incidence in Zambia and South Africa”. 2015. Available at https://www.hptn.org/sites/default/files/2016-05/HPTN?20071-2_Phylogenetics?20Ancillary?20Protocol_v?201.0_15Jan2015.pdf Accessed June17, 2017.
8.
FraserC: Bridging the evolution and epidemiology of HIV in Europe 2014. Available at http://cordis.europa.eu/project/rcn/185401_en.html Accessed June17, 2017.
9.
KozlakidisZ: ICONIC: Infection response through virus genomics 2014. Available at www.hra.nhs.uk/news/research-summaries/infection-response-through-virus-genomics-iconic/ Accessed June17, 2017.
10.
GallA, FernsB, MorrisC, WatsonS, CottenM, RobinsonM, et al.: Universal amplification, next-generation sequencing, and assembly of HIV-1 genomes. J Clin Microbiol, 2012; 50:3838–3844.
11.
GrabowskiMK, LesslerJ, BazaaleJM, ChangLW, NabukaluD, WawerM, et al.: Migration, HIV-infection and access to combination HIV prevention in Rakai district, Uganda. CROI 20162016.
12.
NovitskyV, BussmannH, LoganA, MoyoS, van WidenfeltE, OkuiL, et al.: Phylogenetic relatedness of circulating HIV-1C variants in Mochudi, Botswana. PLoS One, 2013; 8:e80589.
13.
NovitskyV, BussmannH, OkuiL, LoganA, MoyoS, van WidenfeltE, et al.: Estimated age and gender profile of individuals missed by a home-based HIV testing and counselling campaign in a Botswana community. J Int AIDS Soc, 2015; 18:19918.
14.
AsikiG, MpendoJ, AbaasaA, AgabaC, NanvubyaA, NielsenL, et al.: HIV and syphilis prevalence and associated risk factors among fishing communities of Lake Victoria, Uganda. Sex Trans Infect, 2011; 87:511–515.
15.
SeeleyJ, Nakiyingi-MiiroJ, KamaliA, MpendoJ, AsikiG, AbaasaA, et al.: High HIV incidence and socio-behavioral risk patterns in fishing communities on the shores of Lake Victoria, Uganda. Sex Trans Dis, 2012; 39:433–439.
16.
AsikiG, MurphyG, Nakiyingi-MiiroJ, SeeleyJ, NsubugaRN, KarabarindeA, et al.: The general population cohort in rural south-western Uganda: A platform for communicable and non-communicable disease studies. Int J Epidemiol, 2013; 42:129–141.
17.
VandepitteJ, BukenyaJ, WeissHA, NakubulwaS, FrancisSC, HughesP, et al.: HIV and other sexually transmitted infections in a cohort of women involved in high-risk sexual behavior in Kampala, Uganda. Sex Trans Dis, 2011; 38:316–323.
18.
DanaviahS, ManasaJ, WilkinsonE, PillayS, SibisiZ, MsweliS, et al.: Near Full Length HIV-1 Sequencing to Understand HIV Phylodynamics in Africa in Real Time. CROI, Seattle, Washington, 2015.
19.
RokasA, WilliamsBL, KingN, CarrollSB: Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature, 2003; 425:798–804.
PhilippeH, SnellEA, BaptesteE, LopezP, HollandPW, CasaneD: Phylogenomics of eukaryotes: Impact of missing data on large alignments. Mol Biol Evol, 2004; 21:1740–1752.
22.
DelsucF, BrinkmannH, PhilippeH: Phylogenomics and the reconstruction of the tree of life. Nat Rev Genet, 2005; 6:361–375.
23.
RegierJC, ShultzJW, ZwickA, HusseyA, BallB, WetzerR, et al.: Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences. Nature, 2010; 463:1079–1083.
24.
WiensJJ, MorrillMC: Missing data in phylogenetic analysis: Reconciling results from simulations and empirical data. Syst Biol, 2011; 60:719–731.
25.
PhilippeH, BrinkmannH, LavrovDV, LittlewoodDT, ManuelM, WorheideG, et al.: Resolving difficult phylogenetic questions: Why more sequences are not enough. PLoS Biol, 2011; 9:e1000602.
26.
NovitskyV, MoyoS, LeiQ, DeGruttolaV, EssexM: Importance of viral sequence length and number of variable and informative sites in analysis of HIV clustering. AIDS Res Human Retroviruses, 2015; 31:531–542.
27.
YebraG, HodcroftEB, Ragonnet-CroninML, PillayD, Brown AJ; PANGEA HIVConsortium, et al.: Using nearly full-genome HIV sequence data improves phylogeny reconstruction in a simulated epidemic. Sci Rep, 2016; 6:39489.
WilkinsonE, RasmussenDA, RatmannO, StadlerT, EngelbrechtS, de OliveiraT: Origin, imports and exports of HIV-1 subtype C in South Africa: A historical perspective. Infect Genet Evol. 2016; 46:200–208.
30.
GrabowskiMK, LesslerJ, ReddAD, KagaayiJ, LaeyendeckerO, NdyanaboA, et al.: The role of viral introductions in sustaining community-based HIV epidemics in rural Uganda: Evidence from spatial clustering, phylogenetics, and egocentric transmission models. PLoS Med, 2014; 11:e1001610.
PybusOG, RambautA: Evolutionary analysis of the dynamics of viral infectious disease. Nat Rev Genet, 2009; 10:540–550.
33.
WiensJJ: Does adding characters with missing data increase or decrease phylogenetic accuracy?. Syst Biol, 1998; 47:625–640.
34.
RoureB, BaurainD, PhilippeH: Impact of missing data on phylogenies inferred from empirical phylogenomic data sets. Mol Biol Evol, 2013; 30:197–214.
35.
RatmannO, HodcroftE, PicklesM, CoriA, HallM, LycettS, et al.: Phylogenetic tools for generalised HIV epidemics: Findings from the PANGEA-HIV methods comparison. Mol Biol Evol, 2017; 34:185–203.
36.
GallA, MorrisC, KellamP, BerryN: Complete genome sequence of the WHO International Standard for HIV-1 RNA determined by deep sequencing. Genome Announc, 2014; 2:e01254–13.
37.
WymantC, BlanquartF, GallA, BakkerM, BezemerD, CroucherNJ, et al.: Easy and accurate reconstruction of whole HIV genomes from short-read sequence data. bioRxiv 2016.
38.
HuntM, GallA, OngSH, BrenerJ, FernsB, GoulderP, et al.: IVA: Accurate de novo assembly of RNA virus genomes. Bioinformatics, 2015; 31:2374–2376.
39.
KuikenC, FoleyB, LeitnerT, ApetreiV, HahnB, MizrachiI, et al.: HIV Sequence Compendium 2012. In: Theoretical Biology and Biophysics Group (LANL, ed.) NM, LA-UR 12-246532012.
40.
StasinopoulosDM, RigbyRA: Generalized additive models for location scale and shape (GAMLSS) in R. J Stat Softw, 2007; 23.
41.
RambautA, GrasslyNC: Seq-Gen: An application for the Monte Carlo simulation of DNA sequence evolution along phylogenetic trees. Comput Appl Biosci, 1997; 13:235–238.
42.
NguyenLT, SchmidtHA, von HaeselerA, MinhBQ: IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol, 2015; 32:268–274.
43.
ChernomorO, von HaeselerA, MinhBQ: Terrace aware data structure for phylogenomic inference from supermatrices. Syst Biol, 2016; 65:997–1008.
44.
GuindonS, DufayardJF, LefortV, AnisimovaM, HordijkW, GascuelO: New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst Biol, 2010; 59:307–321.
45.
StamatakisA: RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics, 2014; 30:1312–1313.
46.
PriceMN, DehalPS, ArkinAP: FastTree 2—Approximately maximum-likelihood trees for large alignments. PLoS One, 2010; 5:e9490.
47.
ToTH, JungM, LycettS, GascuelO: Fast dating using least-squares criteria and algorithms. Syst Biol, 2015; 65:82–97.
48.
SandA, HoltMK, JohansenJ, FagerbergR, BrodalGS, PedersenCN, et al.: Algorithms for computing the triplet and quartet distances for binary and general trees. Biology (Basel), 2013; 2:1189–1209.
49.
KendallM, ColijnC: Mapping phylogenetic trees to reveal distinct patterns of evolution. Mol Biol Evol, 2016; 33:2735–2743.
50.
SteelM, PennyD: Distributions of tree comparison metrics—Some new results. Syst Biol, 1993; 42:126–141.
51.
EshlemanSH, HudelsonSE, ReddAD, WangL, DebesR, ChenYQ, et al.: Analysis of genetic linkage of HIV from couples enrolled in the HIV Prevention Trials Network 052 trial. J Infect Dis, 2011; 204:1918–1926.
52.
CampbellMS, MullinsJI, HughesJP, CelumC, WongKG, RaugiDN, et al.: Viral linkage in HIV-1 seroconverters and their partners in an HIV-1 prevention clinical trial. PLoS One, 2011; 6:e16986.
53.
ChangLW, GrabowskiMK, SsekubuguR, NalugodaF, KigoziG, NantumeB, et al.: Heterogeneity of the HIV epidemic in agrarian, trading, and fishing communities in Rakai, Uganda: An observational epidemiological study. Lancet HIV, 2016; 3:e388–e396.
54.
YebraG, Ragonnet-CroninM, SsemwangaD, ParryCM, LogueCH, CanePA, et al.: Analysis of the history and spread of HIV-1 in Uganda using phylodynamics. J Gen Virol, 2015; 96(Pt 7):1890–1898.
NovitskyV, KühnertD, MoyoS, WidenfeltE, OkuiL, EssexM: Phylodynamic analysis of HIV sub-epidemics in Mochudi, Botswana. Epidemics, 2015; 13:44–55.
57.
JeffroyO, BrinkmannH, DelsucF, PhilippeH: Phylogenomics: The beginning of incongruence?. Trends Genet, 2006; 22:225–231.
58.
HennMR, BoutwellCL, CharleboisP, LennonNJ, PowerKA, MacalaladAR, et al.: Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection. PLoS Pathog, 2012; 8:e1002529.
59.
ZaniniF, BrodinJ, TheboL, LanzC, BrattG, AlbertJ, et al.: Population genomics of intrapatient HIV-1 evolution. Elife, 2015; 4:e11282.
60.
HedskogC, MildM, JernbergJ, SherwoodE, BrattG, LeitnerT, et al.: Dynamics of HIV-1 quasispecies during antiviral treatment dissected using ultra-deep pyrosequencing. PLoS One, 2010; 5:e11345.
61.
BuzonMJ, CodonerFM, FrostSD, PouC, PuertasMC, MassanellaM, et al.: Deep molecular characterization of HIV-1 dynamics under suppressive HAART. PLoS Pathog, 2011; 7:e1002314.
62.
FischerW, GanusovVV, GiorgiEE, HraberPT, KeeleBF, LeitnerT, et al.: Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing. PLoS One, 2010; 5:e12303.
63.
ReddAD, QuinnTC, TobianAA: Frequency and implications of HIV superinfection. Lancet Infect Dis, 2013; 13:622–628.
64.
BrummeCJ, PoonAF: Promises and pitfalls of Illumina sequencing for HIV resistance genotyping. Virus Res, 2016.
65.
HueS, HassanAS, NabweraH, SandersEJ, PillayD, BerkleyJA, et al.: HIV type 1 in a rural coastal town in Kenya shows multiple introductions with many subtypes and much recombination. AIDS Res Human Retroviruses, 2012; 28:220–224.
66.
WertheimJO, FourmentM, Kosakovsky PondSL: Inconsistencies in estimating the age of HIV-1 subtypes due to heterotachy. Mol Biol Evol, 2012; 29:451–456.
67.
AlizonS, FraserC: Within-host and between-host evolutionary rates across the HIV-1 genome. Retrovirology, 2013; 10:49.
68.
van der KuylAC, BerkhoutB: The biased nucleotide composition of the HIV genome: A constant factor in a highly variable virus. Retrovirology, 2012; 9:92.
69.
LamersS, BarbierA, RatmannO, FraserC, RoseR, LaeyendeckerO, et al.: HIV-1 sequence data coverage in Central East Africa from 1959–2013. AIDS Res Hum Retroviruses. 2016; 32:904–908.
70.
BaurainD, BrinkmannH, PhilippeH: Lack of resolution in the animal phylogeny: Closely spaced cladogeneses or undetected systematic errors?. Mol Biol Evol, 2007; 24:6–9.
71.
TanG, MuffatoM, LedergerberC, HerreroJ, GoldmanN, GilM, et al.: Current methods for automated filtering of multiple sequence alignments frequently worsen single-gene phylogenetic inference. Syst Biol, 2015; 64:778–791.
72.
DidelotX, GardyJ, ColijnC: Bayesian inference of infectious disease transmission from whole-genome sequence data. Mol Biol Evol, 2014; 31:1869–1879.
73.
CornelissenM, GallA, VinkM, ZorgdragerF, BinterS, EdwardsSet al.: From clinical sample to complete genome: Comparing methods for the extraction of HIV-1 RNA for high-throughput deep sequencing. Virus Res., 2016.
74.
NovitskyV, Zahralban-SteeleM, McLaneMF, MoyoS, van WidenfeltE, GaseitsiweS, et al.: Long-range HIV genotyping using viral RNA and proviral DNA for analysis of HIV drug resistance and HIV clustering. J Clin Microbiol, 2015; 53:2581–2592.
75.
BattyEM, WongTH, TrebesA, ArgoudK, AttarM, BuckD, et al.: A modified RNA-Seq approach for whole genome sequencing of RNA viruses from faecal and blood samples. PLoS One, 2013; 8:e66129.
76.
MalboeufCM, YangX, CharleboisP, QuJ, BerlinAM, CasaliM, et al.: Complete viral RNA genome sequencing of ultra-low copy samples by sequence-independent amplification. Nucleic Acids Res, 2013; 41:e13.
77.
BonsallD, AnsariMA, IpC, TrebesA, BrownA, KlenermanP, et al.: ve-SEQ: Robust, unbiased enrichment for streamlined detection and whole-genome sequencing of HCV and other highly diverse pathogens. F1000Res, 2015; 4:1062.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
