Genotypic and Phylogenetic Insights on Prevention of the Spread of HIV-1 and Drug Resistance in "Real-World" Settings.

HIV continues to spread among vulnerable heterosexual (HET), Men-having-Sex with Men (MSM) and intravenous drug user (IDU) populations, influenced by a complex array of biological, behavioral and societal factors. Phylogenetics analyses of large sequence datasets from national drug resistance testing programs reveal the evolutionary interrelationships of viral strains implicated in the dynamic spread of HIV in different regional settings. Viral phylogenetics can be combined with demographic and behavioral information to gain insights on epidemiological processes shaping transmission networks at the population-level. Drug resistance testing programs also reveal emergent mutational pathways leading to resistance to the 23 antiretroviral drugs used in HIV-1 management in low-, middle- and high-income settings. This article describes how genotypic and phylogenetic information from Quebec and elsewhere provide critical information on HIV transmission and resistance, Cumulative findings can be used to optimize public health strategies to tackle the challenges of HIV in "real-world" settings.

1. Introduction

The HIV/AIDS pandemic has had a devastating global impact, eroding public health, social and economic infrastructures worldwide [1,2]. Whereas concentrated subtype B epidemics prevail among men having sex with men (MSM) and intravenous drug users (IDU) in the Western world [3,4], generalized heterosexual (HET) epidemics in Africa and Asia have diversified into an increasingly complex disease, with nine major viral subtypes and more than 70 circulating recombinant forms [1,2]. Migration and globalization patterns have contributed to the introduction and crossover of epidemics in America and Europe, wherein non-B viral subtypes now account for 20% and 60% of new infections in North America and Europe, respectively [5,6].

Outstanding advances in antiretroviral therapy (ART) have transformed HIV/AIDS from a deadly disease to a treatable and potentially preventable disease [7,8]. Since 2010, expanded ART coverage to 20.9 million of the 36.7 million persons living with HIV has resulted in 11% and 47%, declines in new infections among adults and children by respectively (www.unaids.org) [9,10,11]. The goal of ART has shifted from addressing the health benefits of the individual to control of the global HIV-1 pandemic [9,10,11]. The World Health Organization has advanced the “Getting to Zero” policy by 2030 with the “90-90-90” target for 2020, wherein 90% of HIV-infected persons are aware of their status, 90% of infected persons receive immediate ART and 90% of treated individuals achieve long-term viral suppression.

First-line ART in low- and middle-resource settings are largely limited to combination regimens containing two-nucleoside reverse transcriptase (RT) inhibitors (NRTIs) and a non-nucleoside reverse transcriptase inhibitor (NNRTI) [12,13,14]. In resource-rich settings, recent guidelines recommend the combination of two NRTIs with an integrase strand transfer inhibitor (INSTI). The superiority of triple combinations that include an INSTI over an NNRTI or a boosted protease inhibitor (PI) is based on improved tolerability, better drug-drug interactions dosing and higher genetic barrier to resistance [14,15].

HIV drug resistance remains an emerging threat to the long-term success of HIV treatment programs and epidemic control [16]. As more individuals are started on ART, treatment goals will require more robust and durable first-line regimens. Drug resistance testing programs can identify pre-existing and acquired resistance mutations, customizing and optimizing treatment options and clinical management over time [17,18,19,20,21]. Drug resistance testing is recommended at diagnosis, before ART initiation and at the time of virological failure, where possible.

While the role of ART in reducing HIV transmission is undisputed [22], there has been a resurgence of HIV epidemics among MSM and IDU in the Western world and growing rates of transmitted resistance in the Third world. This underscores the need for strong surveillance systems to understand drivers of regional epidemics and improve antiviral treatment and prevention strategies in “real-world” settings [3,23,24,25,26,27,28]. Genotypic drug resistance programs can also be used to provide salient information on the introduction and spread of viral subtypes in different at-risk groups at a population-level networks [27,29,30,31,32,33,34,35,36,37,38].

Phylogenetic analyses of sequence data from national genotyping programs provide an evidence-based framework to deduce the evolutionary interrelationship of viral strains and identify transmission clusters. Phylogenetic surveillance can clarify factors underlying HIV transmission and better inform public health policy [29,30,39,40,41,42,43]. A better understanding of the mechanisms underlying transmission networks may help balance and maximize the efficiency of Treatment-as Prevention strategies.

Mark A. Wainberg was a pioneer in the field of HIV drug resistance, defining signature resistance mutational patterns and related viral fitness costs implicated in the evolution of drug resistance in different HIV-1 viral subtypes. He was also instrumental in advocating for expanded access of the best antiretroviral drugs in Africa [14,15,44]. In this tribute, we provide an update on current knowledge on genotypic and phylogenetic inferences on HIV transmission and drug resistance in Quebec and elsewhere. Cumulative phylogenetic and genotypic findings provide important insights to inform public health interventions in “real-world” settings.

2. Results and Discussion

2.1. Phylogenetic Inferences on the Geographic Introduction and Spread of HIV

The high genetic variability of HIV-1 is a direct consequence of the high replicative rate of the virus, the high error rates of the viral reverse transcriptase (RT) enzyme and viral plasticity under selective host and drug pressure [17,45]. Although HIV continually evolves and adapts in infected hosts, there are stringent bottlenecks leading to transmission of a single monophyletic variant to most (>95%) newly-infected persons. Phylogenetic analysis can track the linkage of viral variants at a population-level, providing a molecular epidemiological framework for surveillance of HIV transmission dynamics. Moreover, phylogenetics can map the geographic introduction and domestic expansion of viral subtypes in different regional settings [29,31,34,46].

The global diversity of major non-B viral subtypes (A, C, D, F, G, H, J, K) and circulating recombinant forms (CRFs) is reflected in non-B subtype epidemic in Quebec (Figure 1). The non-B viral subtypes into Quebec (18% of the provincial epidemic) has arisen through immigration of individuals from francophone countries in Central and West Africa, at the epicenter of the global pandemic. Overall, subtype C, subtype CRF02_AG/G and subtype A/CRF01_AE represent 29%, 28% and 24% of new infections. Subtypes D, F and other CRFs contribute to 6.4%, 2.5% and 10% of the remaining non-B subtype epidemic. Phylogenetics reveal an infrequent onward domestic spread of non-B subtype strains within the province. Two major clusters, circled in Figure 1, depict a Quebec founder subtype A/B CRF (n = 40 members) strain and a subtype CRF01_AE (n = 21 members) cluster spreading among MSM.

Figure 1

Phylogenetic surveillance of the introduction and spread of non-B subtypes in Quebec [47]. First genotypes of new non-B subtypes (n = 1395) reflects the diversity of the African pandemic. Transmission clusters of 5+ members reveal onward spread among Men having Sex with Men (MSM, blue) and Heterosexual (HET, pink) groups. Large clusters (20+ members) are circled.

As illustrated, phylogenetic surveillance can provide important insights on the geographic entry, expansion and crossover of epidemics among HET, MSM and IDU populations [48,49,50,51,52,53]. This can assist in the design of public health interventions, tailored to at-risk vulnerable groups.

2.2. Phylogenetic Analysis of Transmission Clustering

Despite the remarkable global spread of HIV, the relative risk of acquiring infection per exposure is low; with risk estimates of 0.1%, 0.6% and 1.4% per contact for HET, IDU and MSM route of transmission, respectively [54,55]. Cumulative findings in multiple cohorts, across subtypes, show a bottleneck in HIV-1 transmission wherein a single founder establishes HIV infection in any given newly infected partner despite high diversity of viral quasispecies in infected source partners [56,57,58,59].

Drug resistance testing programs provide viral polymerase pol sequence datasets, covering a relatively large proportion of the infected populations at a state/provincial, national, or regional scale [27,28,29,31,35,43,49]. HIV-1 transmission chains are determined based on sequence similarity. Typically, a single (pol) gene sequence from national drug resistance testing program are obtained from each newly HIV-infected person in the population. Phylogenetic trees of viral sequences can be reconstructed using Neighbor-joining, Maximum Likelihood, Bayesian, or Gap approaches [33,49,60]. Transmission clustering of linked viral sequences is generally based on strong bootstrap support (with cutoffs generally exceeding 0.90), short within-cluster genetic distance (0.01–0.05 substitutions/site) or posterior probabilities [33,49].

Cumulative analyses of datasets with high genotypic coverage (>30% of the infected population) have implicated transmission clustering as a major driver in the spread of concentrated HIV-1 epidemics among MSM and IDU populations [29,31,34]. Unfortunately, drug resistance programs are not available or affordable in the majority of low- and middle-income countries ART [61,62]. As such, analysis of epidemics principally driven by viral spread in HET populations such as Africa, is far less understood.

2.3. Phylogenetic Surveillance of Transmission Dynamics of the Quebec MSM Epidemic

In Quebec, the provincial genotyping program has 27,487 pol sequences from 9785 persons, capturing driving forces implicated in the spread of subtype B epidemics in MSM (65%) and IDU (15%) and non-B subtype HET populations (20%) (2002–2016) [27]. Overall, IDU epidemics have declined in the post-ART era. In contrast, the high frequency of viral co-clustering (66%) among newly-diagnosed and ART-naïve MSM infer frequent re-transmissions among newly-infected persons, often unaware of their HIV status.

Genotypic analysis across multiple cohorts has shown the role of recent infection as a driver of onward spread of the provincial epidemic [27,29,30,39,41]. Phylogenetic linkage studies show an increasingly large role of large cluster networks implicated in viral spread, as well as drug-resistant sub-epidemics [27,29,30,39,41]. Half of the epidemic in MSM can be ascribed to “dead-end” transmissions or small clusters of 2–4 individuals expanding over median 4.75 month intervals. The remaining half of the epidemic can be ascribed to increasingly large cluster networks. Phylodynamics reveal that thirty viral strains have sustained the epidemic in MSM. These large cluster networks, involving 20 to 140 persons have risen from 13%, 25%, to 42% of new infections over the 2004–2007, 2008–2011 and 2012–2015 periods, respectively (Figure 2) [27].

Figure 2

Phylogenetic tree of 31 subtype B large cluster (20+ members) transmission networks fueling onward spread of the epidemic among Men having Sex with Men (MSM) (2004–2015).

Integration of phylogenetic, virological and behavioral information show large 20+ clusters outbreaks as compared to singleton transmissions are associated with primary/recent stage infection, younger populations (under 30 years of age) and having more than 5 partnerships prior to PHI (Odds-ratios of 3.7, 3.3 and 1.4, respectively (App. 2) [34].

Overall, 1.7% of viral species (n = 60) have contributed to large cluster networks (10+ members), sustaining the growth of the epidemic among MSM over the last decade. Although it is possible that this can arise partially because of random, stochastic events reflective of the low frequency at which the HIV-1 is transmitted, it may also indicate that there is a selective advantage for variants with certain genotypic and phenotypic features. Our recent findings used in vitro drug selections show a selection bias for large cluster viral variants showing higher replicative fitness under selective drug pressure [63]. Other lines of evidence to support selective transmission of viral species include the observations that (1) early virus populations are less genetically diverse than the source-virus quasispecies populations; (2) viruses present early in infection generally use the chemokine receptor type 5 (CCR5) co-receptor rather than the C-X-C chemokine receptor type 4 (CXCR4) variant for entry; and (3) early viruses are selectively resistant to type 1 interferons [64,65]. These findings indicate that signature sequence characteristics may provide a new opportunity to characterize which biological features of viruses increase their fitness for transmission.

2.4. Phylogenetic Inferences on Public Health Strategies to Control MSM and IDU Epidemics

Phylogenetic research has demonstrated that persons with recent infection drive transmission clustering among MSM and IDU groups; however, estimates are highly variable, representing 10 to 65% of all transmissions [26,27,63,66]. Discrepancies can be related to variations in HIV-1 prevalence (concentrated vs. diffuse epidemics), route of transmission (MSM vs. HET), sexual risk behaviors and depth of population sampling [29,41,49].

Studies in Quebec show a steady decline in IDU from 2002 to 2016, related in part to needle exchange programs, expanded access of ART, community-level declines in viral load. In contrast, rates of new infections among MSM have remained steady over this period. Phylogenetics reveal the decline in new infections associated with small cluster transmission networks have been offset by a relatively few persistent, self-staining and in some cases growing sub-epidemics (n = 30, cluster size 20–140) [27]. This has also been observed in the Netherlands MSM cohort [31].

Numerous groups have shown increased transmission efficiency and a disproportionate role of acute (0–6 weeks) and recent infection (first 6 months) in transmission dynamics [30,52,58]. Phylogenetics reveal how untimely diagnosis, particularly in younger populations can fuel transmission cascades. [52] These findings underscore how the inadequacies of current screening strategies can offset UN 90-90-90 treatment goals. The SPOT (2009–2016) and Actuel-sur-Rue (2012–2013) rapid testing sites in Montreal showed an improved ability to attract MSM (2% HIV+) [27,67]. The SPOT site revealed testing habits among MSM in the province with 50% of persons in Montreal reporting no prior test in the last year. Of note, testing frequency was inversely related to the number of reported partnerships in the three months prior to testing [27].

2.5. Genotypic Analysis Reveals Importance of Drug Regimen Selection in Long-Term HIV Management

The emergence of drug resistance can threaten the long-term benefit of ART regimens in resource-limited settings [14]. The evolution of drug resistance is multifactorial and is dependent on the genetic barrier, the fitness of the mutant variant, the level of resistance that a specific mutation confers and the frequency of that mutation within a given viral population (Figure 3) [17,68,69,70]. Drug resistance testing has enabled personalized strategies for the treatment of HIV-infected individuals in resource-rich settings. Baseline screening at diagnosis and prior to treatment initiation can identify transmitted drug resistance. Selection of treatment regimens can be guided by routine viral load (VL) testing and genotyping upon treatment failure (VL > 400 copies/mL).

Figure 3

Signature genotypic features of HIV-1 variants that are implicated in the genetic barrier to development of resistance (orange), fold-resistance (blue) and viral replicative fitness (green). These include mutations implicated in the resistance to the nucleoside reverse transcriptase (RT) inhibitors (NRTIs), including tenofovir (TDF), zidovudine/AZT (ZDV), stavudine (d4T), didanosine (ddI), non-nucleoside RT inhibitors (NNRTIs), protease inhibitors (PIs), and integrase inhibitors (INIs).

An understanding of the impact of viral subtype in acquired and transmitted drug resistance are essential to guide the clinical management of HIV/AIDS in low- and middle-income settings where drug resistance testing is largely unaffordable [17,20,71]. The first-line therapies available in low- and middle-income setting often includes regimens that are no longer recommended in resource-rich settings. The NRTI components in the first-line ART options in low-income setting, often combine zidovudine (ZDV/AZT), stavudine (d4T) or didanosine (ddI) with lamivudine (3TC), as compared with high-income settings where tenofovir (TDF), or abacavir (ABC) are combined with 3TC or FTC [12,13,14]. The use of d4T-based regimens has been shown in many settings to be particularly prone to toxicity, treatment failure and the development of drug resistance. Switches of regimen in high-income settings are based on genotypic analysis to tailor regimens to individual needs as compared to resource limited settings where shifts are chosen by drug supply.

HIV-1 genetic variation arises from repeated cycles of virus polymerization errors, recombination, APOBEC-mediated RNA editing and selective drug and immune pressure. The evolution of drug resistance is multifactorial, influenced by the genetic barrier to resistance, the replicative fitness and the level of resistance conferred by the acquired resistant variant (Figure 3). Viruses of the NNRTI drug class are replicatively fit with single point mutations conferring high level (>100-fold resistance). Select mutations, including Y181C/I, render group O and all strains of HIV-2 resistant to all drugs within the entire NNRTI class. Although the M184V/I point mutations confer high level resistance (>100 fold) to 3TC or FTC, an essential component most combination regimens, this mutation has high fitness costs, enhancing viral fidelity and resensitizing viruses to ZDV, d4T and TDF [72,73,74]. APOBEC-mediated G-to-A hypermutation—an ancient host defense mechanism responsible for lethal mutagenesis—facilitates the acquisition of select mutations against all drug classes [19,69] (Figure 4).

Figure 4

Impact of the choice of drugs in combination regimen on the development of drug resistance. The frequency of resistance in individuals failing first-line treatment, is influenced by the level (fold) resistance conferred by acquired mutations, the negative impact of resistance mutations on viral fitness.

Recent studies reveal the high frequency of resistance mutations in West/Central and sub-Saharan Africa failing first-line therapy with NNRTI-based regimens after a 2 to 4-years period of treatment [12,13]. Most patients harbored resistance to nevirapine (NVP) (97.5%) or efavirenz (EFV) (75%), acquiring NNRTI -associated resistance mutations, including K103N, Y181C and G190A. There was a facilitated development of V106M in subtype C infections, related in part to a signature natural polymorphism at codon 106 [75]. The higher frequency of E138A, V106M, Y181C and H221Y in subtype C limit the utility of etravirine and rilpivirine in salvage regimens [12]. K65R was more common in subtype C, related to the signature K-K-K codon motif 64–66 codon motif [76,77]. The relatively low frequency of K65R may be related to the high fitness costs. The Q151M nucleoside analogue mutation observed 10% of subtype C isolates, negligible in other subtypes [12,13].

Recent survey data from the World Health Organization have documented an increase in the frequency of transmitted drug resistance (TDR) to NNRTIs related to the scale-up of NVP- or EFV-based regimens (http://www.who.int/hiv/pub/drugresistance/hivdr-report-2017). A rise of transmitted resistance to NNRTIs above 10% will offset the benefit of “treatment-as-prevention” strategies to end the pandemic by 2030. Indeed, 10% of genotyped treatment-naïve individuals acquired a virus resistant to nevirapine or efavirenz in 6 of 11 surveyed countries in sub Saharan Africa and South and Central America.

The rise in TDR to NNRTIs may be more pronounced in MSM or IDU epidemics may also be related to transmission clustering [40,78]. As illustrated in Figure 5, transmission clustering has contributed to the selective spread of resistance to NNRTIs in Quebec. Indeed, the frequency of resistance to NNRTIs in treatment-naïve individuals has exceeded the frequency of resistance in NNRTI-experienced persons since 2009 following community-level declines in viral load associated with the rise in the use of TDF and emtricitabine (FTC) [40,79].

Figure 5

Transmitted resistance to various drug classes. Transmitted resistance to NNRTIs is related to viral fitness of viruses harboring K103N, G190A or Y181C point mutations. Transmitted resistance to NRTIs are related to viruses harboring M41L and revertants at codon 215.

Transmitted resistance to NNRTIs arise due a relatively small number of NNRTI-resistance point mutations (K130N, Y181C, G190A and V106M) responsible for most cases of high-level resistance. This suggests that inexpensive point-mutation assays to detect these mutations may be useful for pre-therapy screening in regions with high levels of TDR [71]. A reliable point-of-care genotypic resistance test could identify which patients should receive standard first-line therapy and which should receive a protease-inhibitor-containing regimen.

2.6. The Need for Better Drugs in Africa

Cumulative findings suggest that the use of drugs with a low genetic barrier (e.g., NNRTIs) and high toxicity (e.g., d4T) have limited the efficacy and durability of the benefit of ART in sustaining long-term viral suppression.

Recent clinical trials have documented the superiority of triple drug regimens that include an INSTI over boosted protease inhibitors (PIs) or NNRTIs [80,81,82,83,84]. Current guidelines have shifted to INSTI-based regimens as the recommended treatment options for ART-naïve populations in resource-rich settings (http://aidsinfo.nih.gov/guidelines). Recent studies have shown the potential utility and cos-effectiveness of dolutegravir (DTG)- and elvitegravir (EVG)-based regimens in low- and middle-income countries [85,86,87,88].

DTG, the newest agent in the INSTI class, has advantages over other INSTIs [89]. Indeed, DTG was the first drug to out-perform EFV over 48 weeks in fully powered randomized trials [81,82]. DTG appears to be impervious to the development of drug resistance in “real-world” settings and this may have important ramifications in the future control of the HIV epidemic. The non-occurrence of drug resistance against DTG has extended to the nucleoside compounds within triple drug regimens [80,82,83,90]. This is noteworthy given that the M184V mutation, conferring resistance to 3TC and FTC is often the first to emerge in the aftermath of treatment failure.

This underscores the potential benefit of DTG as a cornerstone of HIV therapy in all countries in the world including those in sub-Saharan Africa and other developing countries [91,92]. The high barrier to resistance will likely maintain the durability of first-line regimens in settings in which genotypic resistance testing is unavailable. Successful prevention of HIV mother-to-child transmission with dolutegravir-based combination antiretroviral therapy has been observed in a vertically infected pregnant woman harboring multiclass highly drug-resistant HIV-1 [93]. The global shift from EFV-based to DTG-based regimens has begun in settings such as South Africa [91,92,94]. The high barrier to emergent resistance against DTG and INSTIs may mitigate the development of transmitted resistance compared to growing rates of resistance to EFV and other NNRTIs, ranging from 8 to 10% in different regional settings [40,71,95,96,97].

3. Conclusions

Stagnant HIV incidence rates emphasizes the importance of tailoring HIV interventions for select MSM, IDU and HET groups in different regional settings. Phylogenetics can be used to track disease and elucidate dynamic patterns of transmission. Most studies show that persons with undiagnosed recent infection drive onward transmission cascades. The failure to identify new infections in a timely fashion with early treatment intervention may partially explain the persistence of large cluster outbreaks. Public health interventions need to adequately address current screening strategies and linkage to prevention and care, particularly among younger populations. Indeed, there can be substantial overlap between persons in need for testing and those eligible for pre-exposure prophylaxis.

Phylogenetic strategies are advancing the concept of monitoring transmission networks in “real time” [28,34]. A high threshold will be needed to avoid falsely linking viral isolates within “active” transmission chains. Using phylogenetics to understand transmission patterns requires careful consideration of ethics, confidentiality and privacy. Similarly, understanding the impact of migration requires careful attention to avoid further discrimination in this group.

In the era of expanded ART coverage, it also remains imperative to provide the best treatment options available in low- and middle-income settings given the limited access to viral load and drug resistance screening. Newer integrase inhibitors, including dolutegravir and bictegravir appear to be impervious to the development of resistance in the clinic and the laboratory. Incorporation of integrase inhibitors into first-line regimens may be beneficial and cost-effective for the long-term management of HIV/AIDS and realizing “Getting-to-Zero” objectives.