Despite advances in HIV prevention over the past 30 years, an estimated 2.5 million persons were newly infected in 2011, bringing the number of people living with HIV worldwide to 34 million.1 The HIV burden continues to be greatest in sub-Saharan Africa; this region accounts for 75% of all HIV infections and the highest adult prevalence at nearly 5% overall,1 although this average mask extremes reported in subpopulations, some exceeding 50% infected.2 The use of antiretroviral therapy (ART) to reduce viral loads and associated transmission risk among serodiscordant couples (Treatment-as-Prevention, TasP) has garnered excitement as a means to curb the spread of the virus.3 A recent population-based cohort study conducted in a high prevalence region in KwaZulu-Natal found that the risk for HIV acquisition was lowest in areas with the highest ART coverage—providing ecological evidence for real-world effectiveness of TasP.4 As ART coverage has increased in Africa in the last decade, more HIV-infected individuals have been treated, and life expectancy of infected individuals has increased.5 Meanwhile, new transmissions continue, and therefore, overall prevalence of HIV can be expected to increase.6
To sustain an ongoing ART scale-up, and the potential widespread implementation of TasP in the future, expanded financial and public health resources will be required.7,8 However, the most biologically effective and financially efficient way to implement and evaluate prevention measures at the population level is unclear.9,10 Comprehensive knowledge about local epidemics will be required for successful prevention campaigns, including basic data about population demographics, transmission risk groups and viral subtypes, and complex estimates about transmission dynamics, social and sexual mixing networks, and patterns of geographic spread. Importantly, prospective information about success or failure of interventions is essential.
The tools of HIV phylogenetics and molecular epidemiology can be used to understand local transmission dynamics and assist in the design and evaluation of prevention trials. Since early in the epidemic, these approaches have been used to track HIV origin and geographic spread11–14 and in forensic studies evaluating small transmission chains.15–19 New developments, primarily led by the increased availability of viral sequences over the past 20 years, have allowed fine-scale transmission dynamics at the community, regional, and country level to be uncovered.20 The opportunity for such approaches has been facilitated by: (1) the increasing availability of large HIV sequence databases, driven by the routine provision of gene sequence–based drug resistance testing, with sequences linked to epidemiological data (temporal, clinical, demographic, behavioral, or geographic); (2) rapid advances in high-throughput sequencing technology and decreases in sequencing costs; and (3) theoretical and methodological advances in studies of viral transmission, using phylogenetics or genetic network analysis together with linking epidemiological and population genetic models. These advances provide a framework to identify individual traits or stages of infection that are associated with high relative infectiousness; the results of such studies can answer questions not easily resolved with standard epidemiological approaches. In effect, molecular epidemiology tools can help identify traits associated with ongoing HIV transmission (“who” is transmitting the virus?) rather than behaviors or demographic characteristics associated with high rates of infection (who is currently infected?). These analyses are increasing deployed for concentrated epidemics where HIV sequences are routinely available. However, the essential next step is to apply these tools to HIV prevention in generalized epidemics and resource-limited settings where the impact would be greatest. Furthermore, the development of such approaches for HIV can provide a model for new approaches to the transmission dynamics of other pathogens.
SCOPE OF PHYLOGENETIC ANALYSIS TO INFORM MOLECULAR EPIDEMIOLOGY
Molecular epidemiology is the use of genetic data to inform disease etiology and distribution. In HIV, the term can encompass disparate approaches. For this review, we divide the use of genetic analysis to inform HIV epidemiology into 3 general categories: molecular epidemiology, phylodynamics, and phylogeography (Table 1). Although there is generally overlap across these categories within studies, the main questions associated with each category are distinct. Molecular epidemiology allows understanding of the risk factors for HIV transmission and epidemic spread. Phylodynamics reconstructs epidemic history and quantifies epidemic growth or decline, using viral genealogies and explicit population genetic models. Phylogeography describes the distribution of subtype diversity, estimates the impact of human migration on viral spread, and places historical and risk factor data into geographic context, to identify hubs of transmission.
As for any approach to elucidate HIV transmission patterns, gene sequences can only be generated for those already diagnosed and sampled, and therefore, the contribution to transmission from those unsampled and undiagnosed individuals must be considered. Nevertheless, with a well-designed population sampling and sequencing strategy, with linkage to some key demographic, epidemiological, and clinical data, questions from all 3 categories can be addressed for any given population (using the same sequence dataset).
PHYLOGENETIC STUDIES OF TRANSMISSION IN CONCENTRATED EPIDEMICS
Molecular epidemiological tools have been used to address a diverse array of research questions over the past 2 decades, predominantly in settings with concentrated epidemics. These approaches have advanced to describe transmission dynamics at the local, regional, and national level through the post hoc utilization of large sequence datasets linked to epidemiological data. These sequence datasets are largely a consequence of routine antiretroviral drug resistance testing that accompanies initial HIV diagnosis and contains partial HIV pol sequences (typically full protease and partial reverse transcriptase). The pol region has been shown to have sufficient variability to allow for phylogenetic reconstruction to the same extent as the more variable gag and env regions.48 Studies focused on transmission networks in concentrated epidemics have most commonly been conducted in Europe or North America, where large pol sequence datasets linked with epidemiological data exist (Table 2).
Basic Methodological Approach
The basic method to identify HIV transmission dynamics (traits associated with ongoing transmission) is straightforward (Fig. 1): (1) HIV gene sequences are used to reconstruct phylogenetic trees or viral transmission networks; (2) transmission clusters are identified from the networks or phylogenies using ad hoc thresholds for inclusion (eg, minimum pairwise genetic distance between viral gene sequences, or by statistical robustness of the node that defines a phylogenetic clade); and (3) individual traits are evaluated for the strength of their association with cluster membership and used to infer underlying determinants of transmission.63,64 Cluster inclusion thresholds are ad hoc in the sense that there is no widely accepted threshold definition or determining convention. Variation in threshold definitions and approaches can strongly impact the size and number of inferred clusters in a sequence data set, yet many studies do perform sensitivity analyses to quantitatively assess the impact of this variation on cluster identification in their respective data sets. A recent study of transmission dynamics in Brighton, United Kingdom, took this basic methodological approach 1 step further: clinical data were used to determine the most likely transmitter within transmission pairs or clusters, which allowed for likely onward transmission events (and not simply cluster membership based on pairwise genetic distance) to be evaluated for association with individual traits.54
What is the Contribution of Acute and Early Infection to Ongoing Transmissions?
Acute HIV infection is generally defined as the time after HIV acquisition but before seroconversion.65 Early HIV infection represents the first few months of infection, sometimes using a staging system reported by Fiebig et al.66 Most recent results suggest that HIV transmission is increased during acute and early HIV infection, for a period of uncertain duration that may range from weeks to months.67 Multiple phylogenetic studies have focused on acute and early (primary) HIV infection to assess the contribution of acute infection transmissions to overall HIV incidence in a population.54,59,61,68–72 The majority of these studies analyzed only sequences from patients diagnosed during acute or early HIV infection and found variable rates of cluster formation, ranging from 13%71 to greater than 50%61,68 of the study population. These results suggest that acute and early infections are responsible for a disproportionate number of onward transmissions, relative to chronic infections. This might be expected, as: (1) individuals with newly acquired infections remain sexually active73,74; (2) such individuals are unlikely to know their HIV status75; (3) HIV viral load is exceptionally high for weeks to months after infection76; and (4) the transmitted/founder viruses that establish successful infections may have transmission advantages.77–79 However, the proportion of transmission events attributed to acute or primary infection is dependent on definitions used for these stages of infection; more rigorous definitions of early-stage infection, and focusing on individuals with estimated infection dates69 or known dates of seroconversion,70 can help to resolve this issue. It is essential to understand the contribution of early-stage infection to onward transmission, as a high frequency may compromise TasP strategies.80
Other Factors Linked to Transmissions
Trends or patterns associated with onward transmission have been described in greater detail by incorporating sequences from chronically infected patients, thereby increasing the study population sampling density. Using molecular epidemiology tools, transmission network characteristics have been described by various factors, such as race/ancestry and ethnicity,60,62,81,82 transmission risk group,60,62 HIV subtype,49,83,84 transmitted drug resistance (TDR) mutations,53,85,86 and transmission cluster growth.51,52,87 A detailed longitudinal study of men who have sex with men (MSM) sought to find the most likely transmitter for new infections and assess risk factors at the time of transmission.54 Even with 75% of all diagnosed HIV-infected individuals from the local clinic providing pol sequences, the most likely transmitter could be identified in only 25% of those recently infected. This implies that even high sampling fractions of local epidemics may not completely reveal underlying transmission linkage because of, for instance, extra-community or undiagnosed infections as major sources of new transmissions.
Large-Scale Analyses Can Delineate Sub-epidemics
Large-scale analyses of pol sequences have been used to assess the potential influence of viral subtype and regional scale on HIV transmission dynamics. The United Kingdom HIV Drug Resistance Database (UKRDB) (containing >85,000 pol sequences, http://www.hivrdb.org.uk/) has provided data for several transmission studies. Using phylodynamics and a relaxed molecular clock approach, transmission dynamics were reconstructed among MSM56,58 and heterosexuals.57 Among MSM, 6 large transmission clusters were identified, representing separate introductions of subtype B into the United Kingdom in the 1980s.56 Additional analysis within the large MSM clusters (reconstructed with ∼2000 individual sequences) indicated that 25% of transmissions likely occurred within 6 months of infection, with most clusters arising over periods of 3–4 years—an episodic epidemic with multiple clusters of transmission.58 However, in heterosexuals, where the epidemic is dominated by non-B subtypes, much slower transmission dynamics were found.57 This demonstrates that phylogenetic approaches actually reveal the differing dynamics between different risk groups.
The Swiss HIV Cohort includes large repositories of pol sequences that have been used to assess changes in subtype B and non-B transmission clusters over time. Largely independent epidemics of MSM and heterosexuals/injection drug users were noted in the subtype B clusters, but heterosexuals alone did not dominate any of the clusters.52 Over time, the contribution of injection drug use to the heterosexual epidemic notably decreased. The effect of migration was investigated by analyzing non-B subtypes sampled both inside and outside Switzerland. Less than 25% of the non-B subtypes sampled in Switzerland were found in clusters with other Swiss sequences, suggesting that most non-B infections in the country could not be prevented through national prevention measures targeting individuals of only Swiss origin.84
Transmission Network Analysis
The framework of social network analysis provides an alternative method to understand HIV transmission dynamics, which can be supplemental to the explicitly evolutionary approach of phylogenetic analysis. Transmission clusters are reconstructed similar to contact tracing but based on genetic distance metrics. Transmission network parameters among MSM were estimated using phylodynamics of more than 14,000 UKRDB sequences (representing ∼60% of UK MSM).88 Using an inferred network distribution and associated parameter values, the HIV epidemic was characterized by preferential association, which predicted that the epidemic would persist even under conditions of poor overall transmission following a randomly distributed intervention. Cluster growth also differed across the MSM and heterosexual transmission risk groups. However, the relationship between phylogenetic clusters and sexual networks is not direct (they are not one and the same), and the use of phylogenies to understand the underlying transmission network is complex and requires further elucidation.89
PHYLOGENETICS AND TRANSMISSION DYNAMICS IN GENERALIZED EPIDEMICS OF SUB-SAHARAN AFRICA: WHAT QUESTIONS CAN BE ADDRESSED?
HIV transmission dynamic studies in concentrated epidemics largely involve the post hoc use of HIV drug resistance screening datasets whose coverage can reach a substantial proportion of the HIV-infected population. Large datasets of this type are rarely found in Africa—despite sub-Saharan Africa accounting for two thirds of the global HIV infections, only 24% of sequences deposited in the LANL HIV Sequence Database (http://www.hiv.lanl.gov, which receives all HIV sequences deposited in GenBank) are derived from the region (Fig. 2). Notably, the Southern African Treatment and Resistance Network (SATuRN) has a growing database currently with >7000 HIV sequences, albeit sampled from a very large HIV-infected population.90 Nonetheless, informative phylogenetic analyses of generalized epidemics will still require significant de novo sequencing effort. To enable detailed molecular epidemiological studies to inform prevention, clinical, demographic, and behavioral data must be linked to each sequence. This represents a major logistical challenge, given that current sequence databases are generally based on opportunistic approaches to sequence acquisition. Several important questions surrounding HIV epidemiology and prevention can be addressed in the region with available molecular epidemiology tools, particularly when linked to traditional epidemiologic data.
Historical Pattern of Growth or Decline for a Given HIV Epidemic
Molecular epidemiology and phylogeography can be used to understand historical patterns in epidemic origin, spread, and growth over time and space. The high genetic diversity of HIV-1 Group M has given rise to 9 genetically divergent subtypes (A–D, F–H, and J–K), intersubtype recombinants, and circulating recombinant forms.40 Sub-Saharan Africa has the highest HIV genetic diversity and is also characterized by distinct geographical subtype distributions that have remained relatively stable over the past decade.41 Surveillance of genetic diversity has historically supported tracking global epidemiology42 and public health strategies to slow further viral spread.43
The integration of time-stamped sequence (those with known dates of sampling) data with phylogenetics, coalescent models, and molecular clock models allows inferences to be made on the timing of epidemic origin and spread in Africa.13,35–39,45,46,91,92 These analyses can contribute to the design of intervention strategies through better understanding of epidemic growth potential or factors contributing to historical geographic spread of subtypes or recombinants. Phylogeographic approaches showed that HIV-1C in Zimbabwe expanded through multiple introductions originating in southern Africa and localized exponential growth in the 1980s corresponding to demographic and political change.45 Travel accessibility and infrastructure are felt to be critical in epidemic spread, particularly the rapid growth of HIV-A and HIV-D into east Africa35 and the expansion of HIV-1C into east Africa36 and Angola.37
Geographic Source of Local Epidemics or Outbreaks
HIV prevalence in Africa is heterogeneously distributed within countries and often communities; epidemics can be overlapping sub-epidemics defined by geography, time, and a complex interplay of local epidemic drivers. This variation in HIV spread through populations requires that prevention efforts be tailored to characteristics of local epidemics.1,93 For example, the geographic clustering of HIV along roadways in KwaZulu-Natal indicated the need for more intensified interventions within these communities; this study also indicated that many infections were imported from outside local communities.94 High transport connectivity and mobile populations may explain the hyperendemic outbreaks experienced in eastern and southern Africa.95 Understanding the degree of transmission that occurs from outside communities may have a substantial impact on the design and success of targeted prevention efforts such as TasP.
Communities with high transport connectivity may have local HIV epidemics supported by a significant proportion of transmissions from outside communities. Among 153 HIV pol sequences from patients in 1 community in rural coastal Kenya, multiple subtypes and significant recombination were documented.32 In a phylogeographic analysis, many of these sequences were related to different regions in Africa, suggesting multiple introductions into this community, likely reflecting its extensive transport links. In the Rakai district in Uganda, 14,595 individuals in 46 communities underwent extensive HIV surveillance, including spatial and phylogenetic transmission linkage analysis.33 Of 189 HIV-incident cases, an estimated 39% of new cases were infected by household partners, and many new cases that were infected by an extra-household partner were from outside the community. The high degree of external HIV introductions into these communities suggests that the ability of test-and-treat strategies to reduce HIV transmissions may be difficult to measure in small populations unless external introductions are reduced, or unless geographic TasP coverage is sufficiently widespread that migration no longer becomes a relevant problem.
Results from phylogeographic studies must be interpreted in the context of the viral sequences analyzed. Studies based on convenience samples with limited temporal or spatial scales could lead to erroneous inferences about transmission rates between populations (eg, between the studied community and the extra-community). Although this potential limitation can be addressed to some extent with extensive sampling of the community of interest, the use of simulation to validate inferences can be useful.96
Tracking the Transmission and Evolution of Drug Resistance Mutations
The decreases in morbidity and mortality following ART provision97,98 have led to widespread ART rollout over the past decade.99 With increased ART access in Africa, the prevalence of TDR has duly increased.100,101 High levels of acquired drug resistance, often unrecognized when laboratory monitoring for virologic failure is unavailable, and the failure of ART to prevent secondary transmission contribute to increasing TDR in the population.102 Simulations suggest that TDR could have a significant impact on mortality,103 thus highlighting the need for ongoing attention for tracking and preventing TDR. Furthermore, with TasP programs expanding, and World Health Organization guidelines104 recommending treatment initiation at CD4 lymphocyte counts <500 cells per milliliter, then it is likely that an increasing proportion of new infections will be with resistant viruses, despite an overall drop in incidence.105
Although TDR prevalence in sub-Saharan Africa is moderate at an estimated 5.7%,102 it is projected to have increased by 14%–29% per year in southern and east Africa following ART rollout.100 Phylogenetic analyses have been used in TDR prevalence studies primarily to characterize the extensive genetic diversity in specific regions21–23,106,107 and the dominance of HIV-1C in southern Africa.24,26,108–110 Few studies have evaluated transmission lineages of drug-resistant strains through transmission cluster analysis to support or refute epidemiological linkages.21,24,111 Transmission clusters combined with participant life histories revealed a high degree of sexual partner mixing in Ugandan fishing communities and uncovered clusters sharing similar TDR mutations.107 Only 1 study incorporated antiretroviral history in the probable transmitting partners to further characterize linkages among individuals with TDR.111
In contrast to resource-rich settings, drug resistance testing is rarely performed at entry to clinical care, thereby limiting the number of sequences available for phylogenetics. However, sequences will become increasingly available as drug resistance monitoring strategies continue to expand and through ongoing or pending TasP protocols. Phylogenetics could be used to evaluate trends in genetic diversity, further the understanding of sexual networks or transmission clusters harboring resistant strains, track the evolution of drug resistance on a population level, and to help assess the effect of TDR on transmissibility by subtype.
Understanding Transmission Patterns to Help Design Targeted Prevention Measures
Targeting epidemic drivers, or core groups, for enhanced prevention may increase the effectiveness of an intervention but requires detailed epidemiological understanding of viral transmission in the community. This is challenging because HIV transmission dynamics of generalized African epidemics are largely unknown.112–114 Molecular epidemiological approaches can help uncover local HIV epidemic drivers by contributing the links between overlapping sub-epidemics that are characterized by geography, time, and social/sexual interaction.
Characterizing HIV subtypes and including linkage analysis within epidemiological studies can shed light on transmission patterns between high-risk groups and the general population.25,29,34,44 Among MSM in Senegal, HIV phylogenetics revealed different subtype distributions compared with the general population.44 Most of these men reported sex with women, thus likely contributing to bridging between these groups that may modify the subtype distribution. In coastal Kenya, local subtypes among MSM predominated, including those frequent among female sex workers, indicating an epidemic of local origin and confirmation of observed behavioral links between MSM and the general population.25 By combining partnership histories and phylogenetic analysis among Ugandan female sex workers, partial sexual networks and multiple infections were observed, confirming high-risk networks.34 Among HIV concordant heterosexual couples in Senegal, most couples had phylogenetically linked sequences. When combined with interview data, the male partner was often the most likely index, implying concurrency-associated transmission among stable partners.29
In the design of prevention measures, these phylogenetic tools could be expanded to other questions; currently unknown is whether epidemic growth in local epidemics of sub-Saharan Africa is driven primarily by those with high viremia or consistent low-level transmission from chronically infected individuals in concurrent partnerships. Determining epidemic drivers allows for targeted prevention in a more cost-effective manner. For example, Avahan, the India AIDS Initiative, focused prevention on groups at high risk for transmission and acquisition in India (sex workers, their clients, and injection drug users; identified without phylogenetic analyses) leading to more than 100,000 estimated HIV infections prevented in the general population between 2003 and 2008.115,116 The success of these types of approaches, however, depends on having an in-depth understanding of local epidemic drivers in the HIV epidemics of sub-Saharan Africa.
Assessing the Impact of an Intervention
Molecular epidemiology approaches also hold promise in evaluating HIV prevention interventions. At the individual level, phylogenetics has been used to confirm or refute transmission linkages among seroconverters and their partners in heterosexual serodiscordant partnership trials.30,31 By determining genetic linkages between enrolled partners, the primary efficacy of the intervention can be better assessed (if transmission occurred between partners despite the intervention or did the transmission arise from an outside partnership). The Partners in Prevention HSV/HIV Transmission study (PiP) assessed the efficacy of genital herpes suppression in reducing HIV transmission among serodiscordant couples in east and southern Africa.30 Nearly 27% of couples were found to be unlinked through phylogenetic analyses, showing that a substantial number of transmissions occurred through outside partnerships. In HPTN 052, pol and env sequences among index–partner pairs and controls were also evaluated with phylogenetic methods.31 Similar to PiP, 24% of the index cases were not linked to their partner. There was a strong association between linked transmission and the delayed ART initiation study arm: 28 of 29 linked transmission events were in the delayed arm. The association between early ART initiation and transmission reduction became stronger when only the linked events are included in the analysis; this emphasizes the importance of genetic linkage analysis to assess seroconversion events in prevention studies.
Several combination prevention trials that will assess the impact of TasP on a population level are either planned or are currently in progress in sub-Saharan Africa. Most of these trials have planned or are considering integration of molecular epidemiology analyses (Table 3).
At the population level, analysis of genetic data can potentially supplement standard approaches to evaluate the impact of an intervention. Comparative phylogenetic analyses of a baseline trial population and the population over the course of a trial can reveal the emergence or disappearance of clusters associated with particular traits. This approach was used to assess the impact of targeted hepatitis B vaccination in the Netherlands and showed that resulting decreases in hepatitis B virus incidence were due largely to declines in intravenous drug or heterosexual (but not MSM) risk groups.118 However, a follow-up study with increased sample size (n = 894 versus n = 85) suggested that reduced hepatitis B virus transmission was in fact because of reduced incidence in MSM, highlighting the importance of sample size and extended sampling periods for studies of this type.119 Alternatively, gene sequence data can be used to reconstruct HIV transmission networks rather than phylogenies, and cluster size distributions (CSDs) can be compared over the course of a prevention trial or intervention.88,120,121 In theory, CSDs will be dominated by larger clusters in populations where epidemic drivers (individuals with relatively high infectiousness) persist. Changes in CSD can reflect an intervention's impact on particular subgroups or the overall transmission patterns. These population-level approaches may be most suitable for trials in which clear incidence outcomes are equivocal or in which there are clear decreases in incidence but the underlying cause is unknown. Although the methods can be applied to a trial with any targeting strategy, clinical and demographic data from sequenced individuals are required.
Phylodynamics is the use of pathogen sequences to reconstruct epidemic history using viral genealogies and explicit population genetic models.120,122,123 Despite great methodological potential and scientific interest, phylodynamics has to date been rarely used for impact evaluation. This is likely related to a lack of consensus about the interpretation of the estimated parameter Ne, [nominally the effective population size, a quantity proportional to the number of infected individuals, and estimated by coalescent approaches within the product Ne × τ, where τ is the mean (viral) generation time]. Estimates of both Ne and τ can be strongly affected by epidemic stage, transmission dynamics, or population sampling,11,124,125 and come with large variances. There is poor resolution of population size changes in the recent past (∼5 years).125 Additionally, simulation studies have shown that it may be difficult to disentangle the effects of changing incidence and changing transmission networks on phylodynamic parameters125,126; information on transmission network structure might be required for accurate parameter estimation. On a positive note, there is a growing base of modeling approaches to understand the relationships between epidemic models, phylogenetics, and transmission networks, which could be used to better understand how transmission network structure affects phylogenetic trees, and to model outcomes of specific prevention trial designs.89,120,127–131 Of particular interest is the incorporation of stochastic birth–death processes into phylodynamic estimation of epidemic parameters in lieu of standard coalescent models, allowing for more realistic assumptions about changes in epidemic size through time.132–136
CHALLENGES AND SOLUTIONS FOR IMPLEMENTING PHYLOGENETIC STUDIES IN GENERALIZED EPIDEMICS
Phylogenetic analyses of HIV transmission and other epidemiological questions hold great promise to further our understanding of generalized epidemics and inform prevention efforts. However, we must consider how differences between concentrated and generalized epidemics will affect the design and implementation of such studies. Below we note several key challenges that must be met, followed by potential solutions.
HIV Transmission Networks Will be Difficult to Detect With Limited Population Sampling
As seen in phylogenetic studies of concentrated HIV epidemics, a large sampling fraction (eg, >25% of infected individuals in a community) is needed to identify transmission pairs or clusters. This result is seen empirically54 and in simulations.127,128 For phylogenetic studies in generalized HIV epidemics, especially in regions where prevalence can exceed 10%, a substantial number of individuals will need to be sampled and sequenced.
Potential Solution 1
Population sampling that is biased toward sequencing of incident cases can both decrease the required sample fraction and increase the probability of identifying transmission linkages. This “targeted” sampling strategy contrasts with the opportunistic sampling strategy generally found in standard phylogenetic studies. The approach is suited for studies that seek to identify HIV transmissions rather than reconstruct viral evolutionary history. As such, phylogenetic studies of transmission will be most informative and efficient when epidemiological questions are not simply overlaid onto ad hoc phylogenies reconstructed from randomly sampled individuals in a population (the standard approach for phylogenetic studies in, eg, systematics, biogeography, and phylodynamics).
Potential Solution 2
Population sampling conducted over multiple periods, with the initial sample completed before the intervention, increases the probability of identifying transmission pairs or clusters with phylogenetic analysis. Transmission studies in concentrated epidemics generally involve post hoc use of sequence datasets that allows for retrospective analyses of epidemic history or transmission dynamics. The utility of HIV phylogenetics to inform prevention trials cannot be fully realized, however, based solely on retrospective analyses. The sampling frequency required to improve transmission cluster detection is unclear and will likely vary according to local incidence rates.
HIV Sequencing Requirements May Be Extensive
As large sample fractions will be required for informative phylogenetic analyses of generalized epidemics, significant de novo sequencing effort will be necessary, even with targeted sampling of incident infections. This will require extensive technical capability and ability for large volumes of sequences to be generated relatively rapidly. Except for the SATuRN database,90 there are no standing HIV sequence databases in the region that can be readily used for molecular epidemiological studies, in contrast to resource-rich settings that have a larger proportion of sequences deposited compared with their epidemic size (Fig. 2).
Developing a sequencing pipeline is necessary for large-scale HIV phylogenetic and molecular epidemiology projects. Five features of this pipeline are required: (1) the pipeline must work on clinical samples; (2) it must scale across multiple diverse HIV genomes (different subtypes), based on a universal primer set, or alternative methods of genome enrichment that produce equivalent amplification rates across diverse samples; (3) it must scale from hundreds to thousands of samples with ∼75% sequencing success rate across a wide range of genome copy number (viral loads of the individual samples); (4) it must produce accurate consensus sequences with no manual editing; (5) it must detect and accurately quantify minority variants; and (6) it must maximize the informative phylogenetic signal, by extending sequencing from 1 gene (typically the partial pol gene sequenced to around 1 kilobase in length), to the whole viral genome (9.8 kilobase in length).137 A component of this solution may be the development in Africa of the capacity for high-throughput full-genome HIV sequencing.
The development of a pan-African sequence database (analogous to the LANL HIV Sequence Database) will provide an important resource for future studies. The few phylogenetic studies of HIV transmission in Africa to date33,138 suggest that extra-community HIV transmission sources are common. The potential role of a large pan-African database would be to provide sets of African outgroup sequences to identify extra-community sources and to clarify their impact versus undiagnosed local sources. The database would also be useful for studying the spread of intersubtype recombinants and for characterizing the diversity of regional epidemics in Africa, that is, for phylodynamics and phylogeography of the African HIV epidemic.
Methodological Challenges: Integrating Molecular Epidemiology Into Phylodynamics
Molecular epidemiology and phylodynamics have both been areas of active methodological development, as evidenced by the articles referenced above. Nonetheless, current studies in molecular epidemiology that define risk factors for transmission do not make full use of the data available and do not adequately account for the uncertainty and arbitrariness inherent in clustering.
A preferable approach would be to integrate the estimation of transmission risk factors (and other statistics of molecular epidemiology) directly into a phylodynamic inference framework, so that all the data available could be used and arbitrary clustering would not be a prerequisite step. Some authors have begun this process,139 but further methodological development and validation is needed. Additionally, the development of consistent quantitative definitions of transmission clusters, for example, based on tree shape characteristics, such as average branch lengths or nodal support, can make the identification of clusters more rigorous140; this includes assessing the statistical significance of trait clusters by simulation procedures similar to those used to examine gene flow among populations.141
Ethical Challenges: Risks to Individual Privacy and Stigmatization
Ethical challenges in studies involving transmission dynamics in HIV epidemics extend beyond those faced by randomized HIV control trials142 and apply to both concentrated and generalized epidemics. Phylogenetic studies conducted at the village or small community level may involve collection of HIV sequences and individual clinical and demographic data and in some studies may include geographic location data. The risk of individual identification could result in the loss of privacy and, in locations where HIV transmission or reckless exposure is a criminal offense, prosecution. Additionally, an important goal in phylogenetic studies of transmission is to identify traits associated with individuals and groups responsible for onward HIV transmission. Thus, there is the potential for stigmatization of individuals linked with or have common features of transmission network members, either underlying (eg, socioeconomic or demographic group) or proximate (eg, injection–drug use or sexual practice).63,64 In contrast to the stigma related to HIV infection, this challenge will include HIV-negative individuals as well.
Although sampling from generalized epidemics in regions of high prevalence might make such identifications unlikely, principles and governance on patient identifiable data will be necessary. Data from phylogenetic studies of transmission should be reported in ways where individuals cannot be identified. For example, the UKRDB and the Swiss HIV Cohort Study have adopted the strategy where only a minority (eg, 10%) of the sequences collected will be released, with these sequences chosen at random. This includes location data; one approach is to reduce the resolution of location data by including a set number of individuals (eg, 200) in the location set, such that individual identification by location is not possible. These data security strategies will also be addressed in the patient consent process and the tight restriction needed on data release must be recognized by funding agencies and scientific journals.
Opportunities to implement phylogenetic methods at the inception of HIV prevention studies should not be lost. Progress in computational and analytic techniques for reconstructing HIV phylogenies is ongoing. The costs associated with HIV gene or genome sequencing will continue to decrease; rapid, high-throughput sequencing will produce more sequences and larger databases. These advances make it all but certain that sequences or specimens collected in broad reaching studies will eventually be sequenced and used for phylogenetic analyses. However, planning for such analyses from the beginning will maximize their usefulness and the likelihood that phylogenetic analyses can be used in impact evaluations. Additionally, implementing these analyses prospectively will help in identifying hidden subpopulations or core-transmitter groups and in monitoring the spread of TDR, especially among rapidly transmitting networks or clusters.
The authors thank Ward Cates and Nancy Padian for their careful review of earlier versions of the manuscript and Vladimir Novitsky and Mary Kate Grabowski for their helpful discussions. This work stemmed from discussions held during meetings sponsored by the National Institutes of Health (NIH) HIV Prevention Trials Network (February 21–22, 2012) and the Bill and Melinda Gates Foundation (BMGF) (November 1–2, 2012). The authors thank participants of these meetings and specifically David Burns at the NIH and Gina Dallabeta at the BMGF. The authors would also like to acknowledge members of the Phylogenetic and Networks for Generalized HIV Epidemics in Africa (PANGEA), a consortium sponsored by the BMGF.
1. UNAIDS. Joint United Nations Programme on HIV/AIDS. UNAIDS World AIDS Day Report, 2012. Geneva, Switzerland; 2012. Available at: http://www.unaids.org/en/media/unaids/contentassets/documents/epidemiology/2012/gr2012/20121120_UNAIDS_Global_Report_2012_en.pdf
. Accessed July 2, 2013.
2. Karim QA, Kharsany AB, Frohlich JA, et al.. Stabilizing HIV prevalence masks high HIV incidence rates amongst rural and urban women in KwaZulu-Natal, South Africa. Int J Epidemiol. 2011;40:922–930.
3. Cohen MS, Chen YQ, McCauley M, et al.. Prevention of HIV-1
infection with early antiretroviral therapy. N Engl J Med. 2011;365:493–505.
4. Tanser F. High coverage of ART associated with decline in risk of HIV acquisition in rural KwaZulu-Natal, South Africa. Science. 2013;339:966–971.
5. Bor J, Herbst AJ, Newell ML, et al.. Increases in adult life expectancy in rural South Africa: valuing the scale-up of HIV treatment. Science. 2013;339:961–965.
6. Zaidi J, Grapsa E, Tanser F, et al.. Dramatic increases in HIV prevalence after scale-up of antiretroviral treatment: a longitudinal population-based HIV surveillance study in rural Kwazulu-Natal. AIDS. 2013;27:2301–2305.
7. Meyer-Rath G, Over M. HIV treatment as prevention: modelling the cost of antiretroviral treatment—state of the art and future directions. PLoS Med. 2012;9:e1001247.
8. Barnighausen T, Salomon JA, Sangrujee N. HIV treatment as prevention: issues in economic evaluation. PLoS Med. 2012;9:e1001263.
9. Garnett GP, Becker S, Bertozzi S. Treatment as prevention: translating efficacy trial results to population effectiveness. Curr Opin HIV AIDS. 2012;7:157–163.
10. Eaton JW, Johnson LF, Salomon JA, et al.. HIV treatment as prevention: systematic comparison of mathematical models of the potential impact of antiretroviral therapy on HIV incidence in South Africa. PLoS Med. 2012;9:e1001245.
11. Holmes EC, Nee S, Rambaut A, et al.. Revealing the history of infectious disease epidemics through phylogenetic
trees. Philos Trans R Soc Lond B Biol Sci. 1995;349:33–40.
12. Rambaut A, Robertson DL, Pybus OG, et al.. Human immunodeficiency virus. Phylogeny and the origin of HIV-1
. Nature. 2001;410:1047–1048.
13. Gilbert MT, Rambaut A, Wlasiuk G, et al.. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci U S A. 2007;104:18566–18570.
14. Worobey M, Gemmel M, Teuwen DE, et al.. Direct evidence of extensive diversity of HIV-1
in Kinshasa by 1960. Nature. 2008;455:661–664.
15. Ou CY, Ciesielski CA, Myers G, et al.. Molecular epidemiology
of HIV transmission in a dental practice. Science. 1992;256:1165–1171.
16. Scaduto DI, Brown JM, Haaland WC, et al.. Source identification in two criminal cases using phylogenetic
analysis of HIV-1
DNA sequences. Proc Natl Acad Sci U S A. 2010;107:21242–21247.
17. Leitner T, Escanilla D, Franzen C, et al.. Accurate reconstruction of a known HIV-1
transmission history by phylogenetic
tree analysis. Proc Natl Acad Sci U S A. 1996;93:10864–10869.
18. Lemey P, Van Dooren S, Van Laethem K, et al.. Molecular testing of multiple HIV-1
transmissions in a criminal case. AIDS. 2005;19:1649–1658.
19. Pistello M, Del Santo B, Buttò S, et al.. Genetic and phylogenetic
analyses of HIV-1
corroborate the transmission link hypothesis. J Clin Virol. 2004;30:11–18.
20. Brenner BG, Wainberg MA. Future of phylogeny in HIV prevention. J Acquir Immune Defic Syndr. 2013;63(suppl 2):S248–S254.
21. Bartolo I, Casanovas J, Bastos R, et al.. HIV-1
genetic diversity and transmitted drug resistance in health care settings in Maputo, Mozambique. J Acquir Immune Defic Syndr. 2009;51:323–331.
22. Yaotse DA, Nicole V, Roch NF, et al.. Genetic characterization of HIV-1
strains in Togo reveals a high genetic complexity and genotypic drug-resistance mutations in ARV naive patients. Infect Genet Evol. 2009;9:646–652.
23. Caron M, Lekana-Douki SE, Makuwa M, et al.. Prevalence, genetic diversity and antiretroviral drugs resistance-associated mutations among untreated HIV-1
-infected pregnant women in Gabon, central Africa. BMC Infect Dis. 2012;12:64.
24. Deho L, Walwema R, Cappelletti A, et al.. Subtype assignment and phylogenetic
analysis of HIV type 1 strains in patients from Swaziland. AIDS Res Hum Retroviruses. 2008;24:323–325.
25. Tovanabutra S, Sanders EJ, Graham SM, et al.. Evaluation of HIV type 1 strains in men having sex with men and in female sex workers in Mombasa, Kenya. AIDS Res Hum Retroviruses. 2010;26:123–131.
26. Jacobs GB, Laten A, van Rensburg EJ, et al.. Phylogenetic
diversity and low level antiretroviral resistance mutations in HIV type 1 treatment-naive patients from Cape Town, South Africa. AIDS Res Hum Retroviruses. 2008;24:1009–1012.
27. Redd AD, Collinson-Streng A, Martens C, et al.. Identification of HIV superinfection in seroconcordant couples in Rakai, Uganda, by use of next-generation deep sequencing. J Clin Microbiol. 2011;49:2859–2867.
28. Kraft CS, Basu D, Hawkins PA, et al.. Timing and source of subtype-C HIV-1
superinfection in the newly infected partner of Zambian couples with disparate viruses. Retrovirology. 2012;9:22.
29. Jennes W, Kyongo JK, Vanhommerig E, et al.. Molecular epidemiology
transmission in a cohort of HIV-1
concordant heterosexual couples from Dakar, Senegal. PLoS One. 2012;7:e37402.
30. Campbell MS, Mullins JI, Hughes JP, et al.. Viral linkage in HIV-1
seroconverters and their partners in an HIV-1
prevention clinical trial. PLoS One. 2011;6:e16986.
31. Eshleman SH, Hudelson SE, Redd AD, 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.
32. Hue S, Hassan AS, Nabwera H, et al.. HIV type 1 in a rural coastal town in Kenya shows multiple introductions with many subtypes and much recombination. AIDS Res Hum Retroviruses. 2012;28:220–224.
33. Grabowski MK, Lessler J, Redd AD, 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.
34. Ssemwanga D, Ndembi N, Lyagoba F, et al.. HIV type 1 subtype distribution, multiple infections, sexual networks, and partnership histories in female sex workers in Kampala, Uganda. AIDS Res Hum Retroviruses. 2012;28:357–365.
35. Gray RR, Tatem AJ, Lamers S, et al.. Spatial phylodynamics of HIV-1
epidemic emergence in east Africa. AIDS. 2009;23:F9–F17.
36. Delatorre EO, Bello G. Phylodynamics of HIV-1
subtype C epidemic in east Africa. PLoS One. 2012;7:e41904.
37. Afonso JM, Morgado MG, Bello G. Evidence of multiple introductions of HIV-1
subtype C in Angola. Infect Genet Evol. 2012;12:1458–1465.
38. Faria NR, Suchard MA, Abecasis A, et al.. Phylodynamics of the HIV-1
CRF02_AG clade in Cameroon. Infect Genet Evol. 2012;12:453–460.
39. de Silva TI, van Tienen C, Onyango C, et al.. Population dynamics of HIV-2 in rural West Africa: comparison with HIV-1
and ongoing transmission at the heart of the epidemic. AIDS. 2013;27:125–134.
40. Robertson DL, Anderson JP, Bradac JA, et al.. HIV-1
nomenclature proposal. Science. 2000;288:55–56.
41. Lihana RW, Ssemwanga D, Abimiku A, et al.. Update on HIV-1
diversity in Africa: a decade in review. AIDS Rev. 2012;14:83–100.
42. Hemelaar J, Gouws E, Ghys PD, et al.. Global trends in molecular epidemiology
during 2000-2007. AIDS. 2011;25:679–689.
43. Salemi M. Toward a robust monitoring of HIV subtypes distribution worldwide. AIDS. 2011;25:713–714.
44. Ndiaye HD, Toure-Kane C, Vidal N, et al.. Surprisingly high prevalence of subtype C and specific HIV-1
subtype/CRF distribution in men having sex with men in Senegal. J Acquir Immune Defic Syndr. 2009;52:249–252.
45. Dalai SC, de Oliveira T, Harkins GW, et al.. Evolution and molecular epidemiology
of subtype C HIV-1
in Zimbabwe. AIDS. 2009;23:2523–2532.
46. Esbjornsson J, Mild M, Mansson F, et al.. HIV-1 molecular epidemiology
in Guinea-Bissau, West Africa: origin, demography and migrations. PLoS One. 2011;6:e17025.
[computer program]. Version 1.4. 2009. Available at: http://tree.bio.ed.ac.uk/software/figtree/
. Accessed July 15, 2014.
48. Hue S, Clewley JP, Cane PA, et al.. HIV-1
pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy. AIDS. 2004;18:719–728.
49. Chalmet K, Staelens D, Blot S, et al.. Epidemiological study of phylogenetic
transmission clusters in a local HIV-1
epidemic reveals distinct differences between subtype B and non-B infections. BMC Infect Dis. 2010;10:262.
50. González-Alba JM, Holguín Á, Garcia R, et al.. Molecular surveillance of HIV-1
in Madrid, Spain: a phylogeographic analysis. J Virol. 2011;85:10755–10763.
51. Ambrosioni J, Junier T, Delhumeau C, et al.. Impact of highly active antiretroviral therapy on the molecular epidemiology
of newly diagnosed HIV infections. AIDS. 2012;26:2079–2086.
52. Kouyos RD, von Wyl V, Yerly S, et al.. Molecular epidemiology
reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. J Infect Dis. 2010;201:1488–1497.
53. Yerly S, Junier T, Gayet-Ageron A, et al.. The impact of transmission clusters on primary drug resistance in newly diagnosed HIV-1
infection. AIDS. 2009;23:1415–1423.
54. Fisher M, Pao D, Brown AE, et al.. Determinants of HIV-1
transmission in men who have sex with men: a combined clinical, epidemiological and phylogenetic
approach. AIDS. 2010;24:1739–1747.
55. Gifford RJ, de Oliveira T, Rambaut A, et al.. Phylogenetic
surveillance of viral genetic diversity and the evolving molecular epidemiology
of human immunodeficiency virus type 1. J Virol. 2007;81:13050–13056.
56. Hue S, Pillay D, Clewley JP, et al.. Genetic analysis reveals the complex structure of HIV-1
transmission within defined risk groups. Proc Natl Acad Sci U S A. 2005;102:4425–4429.
57. Hughes GJ, Fearnhill E, Dunn D, et al.. Molecular phylodynamics of the heterosexual HIV epidemic in the United Kingdom. PLoS Pathog. 2009;5:e1000590.
58. Lewis F, Hughes GJ, Rambaut A, et al.. Episodic sexual transmission of HIV revealed by molecular phylodynamics. PLoS Med. 2008;5:e50.
59. Pao D, Fisher M, Hue S, et al.. Transmission of HIV-1
during primary infection: relationship to sexual risk and sexually transmitted infections. AIDS. 2005;19:85–90.
60. Aldous JL, Pond SK, Poon A, et al.. Characterizing HIV transmission networks
across the United States. Clin Infect Dis. 2012;55:1135–1143.
61. Brenner BG, Roger M, Routy JP, et al.. High rates of forward transmission events after acute/early HIV-1
infection. J Infect Dis. 2007;195:951–959.
62. Dennis AM, Hue S, Hurt CB, et al.. Phylogenetic
insights into regional HIV transmission. AIDS. 2012;26:1813–1822.
63. Boerma JT, Weir SS. Integrating demographic and epidemiological approaches to research on HIV/AIDS: the proximate-determinants framework. J Infect Dis. 2005;191(suppl 1):S61–S67.
64. Lewis JJ, Donnelly CA, Mare P, et al.. Evaluating the proximate determinants framework for HIV infection in rural Zimbabwe. Sex Transm Infect. 2007;83(suppl 1):i61–i69.
65. Pilcher CD, Eron JJ Jr, Galvin S, et al.. Acute HIV revisited: new opportunities for treatment and prevention. J Clin Invest. 2004;113:937–945.
66. Fiebig EW, Wright DJ, Rawal BD, et al.. Dynamics of HIV viremia and antibody seroconversion in plasma donors: implications for diagnosis and staging of primary HIV infection. AIDS. 2003;17:1871–1879.
67. Powers KA, Ghani AC, Miller WC, et al.. The role of acute and early HIV infection in the spread of HIV and implications for transmission prevention strategies in Lilongwe, Malawi: a modelling study. Lancet. 2011;378:256–268.
68. Brenner BG, Roger M, Stephens D, et al.. Transmission clustering drives the onward spread of the HIV epidemic among men who have sex with men in Quebec. J Infect Dis. 2011;204:1115–1119.
69. Brown AE, Gifford RJ, Clewley JP, et al.. Phylogenetic
reconstruction of transmission events from individuals with acute HIV infection: toward more-rigorous epidemiological definitions. J Infect Dis. 2009;199:427–431.
70. Chibo D, Kaye M, Birch C. HIV transmissions during seroconversion contribute significantly to new infections in men who have sex with men in Australia. AIDS Res Hum Retroviruses. 2012;28:460–464.
71. Frange P, Meyer L, Deveau C, et al.. Recent HIV-1
infection contributes to the viral diffusion over the French territory with a recent increasing frequency. PLoS One. 2012;7:e31695.
72. Bezemer D, van Sighem A, Lukashov VV, et al.. Transmission networks
among men having sex with men in the Netherlands. AIDS. 2010;24:271–282.
73. Weinhardt LS, Kelly JA, Brondino MJ, et al.. HIV transmission risk behavior among men and women living with HIV in 4 cities in the United States. J Acquir Immune Defic Syndr. 2004;36:1057–1066.
74. Avants SK, Warburton LA, Hawkins KA, et al.. Continuation of high-risk behavior by HIV-positive drug users. Treatment implications. J Subst Abuse Treat. 2000;19:15–22.
75. Daar ES. Clinical presentation and diagnosis of primary HIV-1
infection. Curr Opin HIV AIDS. 2008;3:10–15.
76. Vergis EN, Mellors JW. Natural history of HIV-1
infection. Infect Dis Clin North Am. 2000;14:809–825, v–vi.
77. Lythgoe KA, Fraser C. New insights into the evolutionary rate of HIV-1
at the within-host and epidemiological levels. Proc Biol Sci. 2012;279:3367–3375.
78. Alizon S, Fraser C. Within-host and between-host evolutionary rates across the HIV-1
genome. Retrovirology. 2013;10:49.
79. Parrish NF, Gao F, Li H, et al.. Phenotypic properties of transmitted founder HIV-1
. Proc Natl Acad Sci U S A. 2013;110:6626–6633.
80. Cohen MS, Dye C, Fraser C, et al.. HIV treatment as prevention: debate and commentary—will early infection compromise treatment-as-prevention strategies? PLoS Med. 2012;9:e1001232.
81. Oster AM, Pieniazek D, Zhang X, et al.. Demographic but not geographic insularity in HIV transmission among young black MSM. AIDS. 2011;25:2157–2165.
82. Kramer MA, Cornelissen M, Paraskevis D, et al.. HIV transmission patterns among The Netherlands, Suriname, and The Netherlands Antilles: a molecular epidemiological study. AIDS Res Hum Retroviruses. 2010;27:123–130.
83. Callegaro A, Svicher V, Alteri C, et al.. Epidemiological network analysis in HIV-1
B infected patients diagnosed in Italy between 2000 and 2008. Infect Genet Evol. 2011;11:624–632.
84. von Wyl V, Kouyos RD, Yerly S, et al.. The role of migration and domestic transmission in the spread of HIV-1
non-B subtypes in Switzerland. J Infect Dis. 2011;204:1095–1103.
85. Kaye M, Chibo D, Birch C. Phylogenetic
investigation of transmission pathways of drug-resistant HIV-1
utilizing pol sequences derived from resistance genotyping. J Acquir Immune Defic Syndr. 2008;49:9–16.
86. Hue S, Gifford RJ, Dunn D, et al.; on Behalf of the UK Collaborative Group on HIV Drug Resistance. Demonstration of Sustained drug-resistant human immunodeficiency virus type 1 lineages circulating among treatment-naive individuals. J Virol. 2009;83:2645–2654.
87. Ragonnet-Cronin M, Ofner-Agostini M, Merks H, et al.. Longitudinal phylogenetic
surveillance identifies distinct patterns of cluster dynamics. J Acquir Immune Defic Syndr. 2010;55:102–108.
88. Leigh Brown AJ, Lycett SJ, Weinert L, et al.. Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J Infect Dis. 2011;204:1463–1469.
89. Robinson K, Fyson N, Cohen T, et al.. How the dynamics and structure of sexual contact networks shape pathogen phylogenies. PLoS Comput Biol. 2013;9:e1003105.
90. de Oliveira T, Shafer RW, Seebregts C. Public database for HIV drug resistance in southern Africa. Nature. 2010;464:673.
91. Bello G, Afonso JM, Morgado MG. Phylodynamics of HIV-1
subtype F1 in Angola, Brazil and Romania. Infect Genet Evol. 2012;12:1079–1086.
92. Lemey P, Pybus OG, Rambaut A, et al.. The molecular population genetics of HIV-1
group O. Genetics. 2004;167:1059–1068.
93. Hankins CA, de Zalduondo BO. Combination prevention: a deeper understanding of effective HIV prevention. AIDS. 2010;24:S70–S80.
94. Tanser F, Barnighausen T, Cooke GS, et al.. Localized spatial clustering of HIV infections in a widely disseminated rural South African epidemic. Int J Epidemiol. 2009;38:1008–1016.
95. Tatem AJ, Hemelaar J, Gray RR, et al.. Spatial accessibility and the spread of HIV-1
subtypes and recombinants. AIDS. 2012;26:2351–2360.
96. Carnegie NB, Wang R, Novitsky V, et al.. Linkage of viral sequences among HIV-infected village residents in Botswana: estimation of linkage rates in the presence of missing data. PLoS Comput Biol. 2014;10:e1003430.
97. Hogg RS, Yip B, Chan KJ, et al.. Rates of disease progression by baseline CD4 cell count and viral load after initiating triple-drug therapy. JAMA. 2001;286:2568–2577.
98. Egger M, May M, Chene G, et al.. Prognosis of HIV-1
-infected patients starting highly active antiretroviral therapy: a collaborative analysis of prospective studies. Lancet. 2002;360:119.
99. Gilks CF, Crowley S, Ekpini R, et al.. The WHO public-health approach to antiretroviral treatment against HIV in resource-limited settings. Lancet. 2006;368:505–510.
100. Gupta RK, Jordan MR, Sultan BJ, et al.. Global trends in antiretroviral resistance in treatment-naive individuals with HIV after rollout of antiretroviral treatment in resource-limited settings: a global collaborative study and meta-regression analysis. Lancet. 2012;380:1250–1258.
101. Hamers RL, Sigaloff KC, Kityo C, et al.. HIV-1
drug resistance in antiretroviral-naive patients in sub-Saharan Africa
. Lancet Infect Dis. 2013;13:196–197.
102. Stadeli KM, Richman DD. Rates of emergence of HIV drug resistance in resource-limited settings: a systematic review. Antivir Ther. 2013;18:115–123.
103. Cambiano V, Bertagnolio S, Jordan MR, et al.. Transmission of drug resistant HIV and its potential impact on mortality and treatment outcomes in resource-limited settings. J Infect Dis. 2013;207(suppl 2):S57–S62.
104. World Health Organization. Consolidated guidelines on general HIV care and the use of antiretroviral drugs for treating and preventing HIV infection: recommendations for a public health approach. 2013. Available at: http://www.who.int/hiv/pub/guidelines/arv2013/download/en/
. Accessed July 15, 2014.
105. Cambiano V, Bertagnolio S, Jordan MR, et al.. Predicted levels of HIV drug resistance in South Africa: potential impact of expanding diagnosis, retention and eligibility criteria for antiretroviral therapy initiation. AIDS. 2014;28(suppl 1):S15–S23.
106. Sigaloff KC, Mandaliya K, Hamers RL, et al.. Short communication: high prevalence of transmitted antiretroviral drug resistance among newly HIV type 1 diagnosed adults in Mombasa, Kenya. AIDS Res Hum Retroviruses. 2012;28:1033–1037.
107. Nazziwa J, Njai HF, Ndembi N, et al.. Short communication: HIV type 1 transmitted drug resistance and evidence of transmission clusters among recently infected antiretroviral-naive individuals from Ugandan fishing communities of Lake Victoria. AIDS Res Hum Retroviruses. 2013;29:788–795.
108. Bessong PO, Mphahlele J, Choge IA, et al.. Resistance mutational analysis of HIV type 1 subtype C among rural South African drug-naive patients prior to large-scale availability of antiretrovirals. AIDS Res Hum Retroviruses. 2006;22:1306–1312.
109. Bussmann H, de la Hoz Gomez F, Roels TH, et al.. Prevalence of transmitted HIV drug resistance in Botswana: lessons learned from the HIVDR-Threshold Survey conducted among women presenting for routine antenatal care as part of the 2007 national sentinel survey. AIDS Res Hum Retroviruses. 2011;27:365–372.
110. Nwobegahay J, Bessong P, Masebe T, et al.. Prevalence of drug-resistant mutations in newly diagnosed drug-naive HIV-1
-infected individuals in a treatment site in the Waterberg district, Limpopo province. S Afr Med J. 2011;101:335–337.
111. Price MA, Wallis CL, Lakhi S, et al.. Transmitted HIV type 1 drug resistance among individuals with recent HIV infection in East and Southern Africa. AIDS Res Hum Retroviruses. 2011;27:5–12.
112. Chemaitelly H, Shelton JD, Hallett TB, et al.. Only a fraction of new HIV infections occur within identifiable stable discordant couples in sub-Saharan Africa
. AIDS. 2013;27:251–260.
113. Morris M. Barking up the wrong evidence tree. Comment on Lurie and Rosenthal, “Concurrent partnerships as a driver of the HIV epidemic in sub-Saharan Africa
? The evidence is limited”. AIDS Behav. 2010;14:31–33; discussion 34–37.
114. Lurie M, Rosenthal S, Williams B. Concurrency driving the African HIV epidemics: where is the evidence? Lancet. 2009;374:1420; author reply 1420–1421.
115. Ng M, Gakidou E, Levin-Rector A, et al.. Assessment of population-level effect of Avahan, an HIV-prevention initiative in India. Lancet. 2011;378:1643–1652.
116. Boily MC, Pickles M, Lowndes CM, et al.. Positive impact of a large-scale HIV prevention program among female sex workers and clients in Karnataka state, India. AIDS. 2013;27:1449–1460.
117. Boily MC, Masse B, Alsallaq R, et al.. HIV treatment as prevention: considerations in the design, conduct, and analysis of cluster randomized controlled trials of combination HIV prevention. PLoS Med. 2012;9:e1001250.
118. van Houdt R, Bruisten SM, Koedijk FD, et al.. Molecular epidemiology
of acute hepatitis B in the Netherlands in 2004: nationwide survey. J Med Virol. 2007;79:895–901.
119. Hahne S, van Houdt R, Koedijk F, et al.. Selective hepatitis B virus vaccination has reduced hepatitis B virus transmission in the Netherlands. PLoS One. 2013;8:e67866.
120. Volz EM, Kosakovsky Pond SL, Ward MJ, et al.. Phylodynamics of infectious disease epidemics. Genetics. 2009;183:1421–1430.
121. Wertheim JO, Kosakovsky Pond SL, Little SJ, et al.. Using HIV transmission networks
to investigate community effects in HIV prevention trials. PLoS One. 2011;6:e27775.
122. Grenfell BT, Pybus OG, Gog JR, et al.. Unifying the epidemiological and evolutionary dynamics of pathogens. Science. 2004;303:327–332.
123. Drummond AJ, Rambaut A, Shapiro B, et al.. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:1185–1192.
124. Frost SD, Volz EM. Viral phylodynamics and the search for an “effective number of infections”. Philos Trans R Soc Lond B Biol Sci. 2010;365:1879–1890.
125. de Silva E, Ferguson NM, Fraser C. Inferring pandemic growth rates from sequence data. J R Soc Interface. 2012;9:1797–1808.
126. O'Dea EB, Wilke CO. Contact heterogeneity and phylodynamics: how contact networks shape parasite evolutionary trees. Interdiscip Perspect Infect Dis. 2011;2011:238743.
127. Volz EM, Koopman JS, Ward MJ, et al.. Simple epidemiological dynamics explain phylogenetic
clustering of HIV from patients with recent infection. PLoS Comput Biol. 2012;8:e1002552.
128. Frost SD, Volz EM. Modelling tree shape and structure in viral phylodynamics. Philos Trans R Soc Lond B Biol Sci. 2013;368:20120208.
129. Leventhal GE, Kouyos R, Stadler T, et al.. Inferring epidemic contact structure from phylogenetic
trees. PLoS Comput Biol. 2012;8:e1002413.
130. Magiorkinis G, Sypsa V, Magiorkinis E, et al.. Integrating phylodynamics and epidemiology to estimate transmission diversity in viral epidemics. PLoS Comput Biol. 2013;9:e1002876.
131. Ratmann O, Donker G, Meijer A, et al.. Phylodynamic inference and model assessment with approximate Bayesian computation: influenza as a case study. PLoS Comput Biol. 2012;8:e1002835.
132. Volz EM. Complex population dynamics and the coalescent under neutrality. Genetics. 2012;190:187–201.
133. Rasmussen DA, Ratmann O, Koelle K. Inference for nonlinear epidemiological models using genealogies and time series. PLoS Comput Biol. 2011;7:e1002136.
134. Leventhal GE, Gunthard HF, Bonhoeffer S, et al.. Using an epidemiological model for phylogenetic
inference reveals density dependence in HIV transmission. Mol Biol Evol. 2014;31:6–17.
135. Kuhnert D, Stadler T, Vaughan TG, et al.. Simultaneous reconstruction of evolutionary history and epidemiological dynamics from viral sequences with the birth-death SIR model. J R Soc Interface. 2014;11:20131106.
136. Stadler T. On incomplete sampling under birth-death models and connections to the sampling-based coalescent. J Theor Biol. 2009;261:58–66.
137. Gall A, Ferns B, Morris C, et al.. Universal amplification, next-generation sequencing, and assembly of HIV-1
genomes. J Clin Microbiol. 2012;50:3838–3844.
138. Novitsky V, Wang R, Logan A, et al.. The role of HIV-1
viral linkage and viral load in treatment-as-prevention studies. Paper presented at: HIV Dynamics and Evolution April 25, 2012; Asheville, NC.
139. Stadler T, Kuhnert D, Bonhoeffer S, et al.. Birth-death skyline plot reveals temporal changes of epidemic spread in HIV and hepatitis C virus (HCV). Proc Natl Acad Sci U S A. 2013;110:228–233.
140. Chevenet F, Jung M, Peeters M, et al.. Searching for virus phylotypes. Bioinformatics. 2013;29:561–570.
141. Slatkin M, Maddison WP. A cladistic measure of gene flow inferred from the phylogenies of alleles. Genetics. 1989;123:603–613.
142. Osrin D, Azad K, Fernandez A, et al.. Ethical challenges in cluster randomized controlled trials: experiences from public health interventions in Africa and Asia. Bull World Health Organ. 2009;87:772–779.