Addressing Ethical Challenges in US-Based HIV Phylogenetic Research.

In recent years, phylogenetic analysis of HIV sequence data has been used in research studies to investigate transmission patterns between individuals and groups, including analysis of data from HIV prevention clinical trials, in molecular epidemiology, and in public health surveillance programs. Phylogenetic analysis can provide valuable information to inform HIV prevention efforts, but it also has risks, including stigma and marginalization of groups, or potential identification of HIV transmission between individuals. In response to these concerns, an interdisciplinary working group was assembled to address ethical challenges in US-based HIV phylogenetic research. The working group developed recommendations regarding (1) study design; (2) data security, access, and sharing; (3) legal issues; (4) community engagement; and (5) communication and dissemination. The working group also identified areas for future research and scholarship to promote ethical conduct of HIV phylogenetic research.

HIV Diversity and Genetic Compartmentalization in Blood and Testes during Suppressive Antiretroviral Therapy.

HIV's ability to persist during suppressive antiretroviral therapy is the main barrier to cure. Immune-privileged tissues, such as the testes, may constitute distinctive sites of HIV persistence, but this has been challenging to study in humans. We analyzed the proviral burden and genetics in the blood and testes of 10 individuals on suppressive therapy who underwent elective gender-affirming surgery. HIV DNA levels in matched blood and testes were quantified by quantitative PCR, and subgenomic proviral sequences (nef region) were characterized from single templates. HIV diversity, compartmentalization, and immune escape burden were assessed using genetic and phylogenetic approaches. Diverse proviruses were recovered from the blood (396 sequences; 354 nef-intact sequences) and testes (326 sequences; 309 nef-intact sequences) of all participants. Notably, the frequency of identical HIV sequences varied markedly between and within individuals. Nevertheless, proviral loads, within-host unique HIV sequence diversity, and the immune escape burden correlated positively between blood and testes. When all intact nef sequences were evaluated, 60% of participants exhibited significant blood-testis genetic compartmentalization, but none did so when the evaluation was restricted to unique sequences per site, suggesting that compartmentalization, when present, is attributable to the clonal expansion of HIV-infected cells. Our observations confirm the testes as a site of HIV persistence and suggest that individuals with larger and more diverse blood reservoirs will have larger and more diverse testis reservoirs. Furthermore, while the testis microenvironment may not be sufficiently unique to facilitate the seeding of unique viral populations therein, differential clonal expansion dynamics may be at play, which may complicate HIV eradication.IMPORTANCE Two key questions in HIV reservoir biology are whether immune-privileged tissues, such as the testes, harbor distinctive proviral populations during suppressive therapy and, if so, by what mechanism. While our results indicated that blood-testis HIV genetic compartmentalization was reasonably common (60%), it was always attributable to differential frequencies of identical HIV sequences between sites. No blood-tissue data set retained evidence of compartmentalization when only unique HIV sequences per site were considered; moreover, HIV immune escape mutation burdens were highly concordant between sites. We conclude that the principal mechanism by which blood and testis reservoirs differ is not via seeding of divergent HIV sequences therein but, rather, via differential clonal expansion of latently infected cells. Thus, while viral diversity and escape-related barriers to HIV eradication are of a broadly similar magnitude across the blood and testes, clonal expansion represents a challenge. The results support individualized analysis of within-host reservoir diversity to inform curative approaches.

node.dating: dating ancestors in phylogenetic trees in R.

Summary: Phylogenetic trees encode the evolutionary distances between species or populations. With sufficient information, these evolutionary distances can be rescaled over time to provide estimates of the dates of the most recent ancestors of the species. Here we present the R program node.dating, divergence-time analysis software, which uses a maximum-likelihood method to estimate the dates of the internal nodes of a phylogenetic tree. Availability and Implementation: node.dating is available as a part of the R v3.30 package ape v4.0 (cran.r-project.org). node.dating is also available in the GitHub repository: https://github.com/brj1/node.dating , along with supplementary software and tests. Contact: brj1@sfu.ca. Supplementary information: Supplementary data are available at Bioinformatics online.

Promises and pitfalls of Illumina sequencing for HIV resistance genotyping.

Genetic sequencing ("genotyping") plays a critical role in the modern clinical management of HIV infection. This virus evolves rapidly within patients because of its error-prone reverse transcriptase and short generation time. Consequently, HIV variants with mutations that confer resistance to one or more antiretroviral drugs can emerge during sub-optimal treatment. There are now multiple HIV drug resistance interpretation algorithms that take the region of the HIV genome encoding the major drug targets as inputs; expert use of these algorithms can significantly improve to clinical outcomes in HIV treatment. Next-generation sequencing has the potential to revolutionize HIV resistance genotyping by lowering the threshold that rare but clinically significant HIV variants can be detected reproducibly, and by conferring improved cost-effectiveness in high-throughput scenarios. In this review, we discuss the relative merits and challenges of deploying the Illumina MiSeq instrument for clinical HIV genotyping.

Near real-time monitoring of HIV transmission hotspots from routine HIV genotyping: an implementation case study.

BACKGROUND: HIV evolves rapidly and therefore infections with similar genetic sequences are likely linked by recent transmission events. Clusters of related infections can represent subpopulations with high rates of transmission. We describe the implementation of an automated near real-time system to monitor and characterise HIV transmission hotspots in British Columbia, Canada. METHODS: In this implementation case study, we applied a monitoring system to the British Columbia drug treatment database, which holds more than 32 000 anonymised HIV genotypes for nearly 9000 residents of British Columbia living with HIV. On average, five to six new HIV genotypes are deposited in the database every day, which triggers an automated reanalysis of the entire database. We extracted clusters of five or more individuals with short phylogenetic distances between their respective HIV sequences. The system generated monthly reports of the growth and characteristics of clusters that were distributed to public health officers. FINDINGS: In June, 2014, the monitoring system detected the expansion of a cluster by 11 new cases during 3 months, including eight cases with transmitted drug resistance. This cluster generally comprised young men who have sex with men. The subsequent report precipitated an enhanced public health follow-up to ensure linkage to care and treatment initiation in the affected subpopulation. Of the nine cases associated with this follow-up, all had already been linked to care and five cases had started treatment. Subsequent to the follow-up, three additional cases started treatment and most cases achieved suppressed viral loads. During the next 12 months, we detected 12 new cases in this cluster with reduction in the onward transmission of drug resistance. INTERPRETATION: Our findings show the first application of an automated phylogenetic system monitoring a clinical database to detect a recent HIV outbreak and support the ensuing public health response. By making secondary use of routinely collected HIV genotypes, this approach is cost-effective, attains near real-time monitoring of new cases, and can be implemented in all settings in which HIV genotyping is the standard of care. FUNDING: BC Centre for Excellence in HIV/AIDS, the Canadian Institutes for Health Research, the Genome Canada-CIHR Partnership in Genomics and Personalized Health, and the US National Institute on Drug Abuse.

Phylodynamic Inference with Kernel ABC and Its Application to HIV Epidemiology.

The shapes of phylogenetic trees relating virus populations are determined by the adaptation of viruses within each host, and by the transmission of viruses among hosts. Phylodynamic inference attempts to reverse this flow of information, estimating parameters of these processes from the shape of a virus phylogeny reconstructed from a sample of genetic sequences from the epidemic. A key challenge to phylodynamic inference is quantifying the similarity between two trees in an efficient and comprehensive way. In this study, I demonstrate that a new distance measure, based on a subset tree kernel function from computational linguistics, confers a significant improvement over previous measures of tree shape for classifying trees generated under different epidemiological scenarios. Next, I incorporate this kernel-based distance measure into an approximate Bayesian computation (ABC) framework for phylodynamic inference. ABC bypasses the need for an analytical solution of model likelihood, as it only requires the ability to simulate data from the model. I validate this "kernel-ABC" method for phylodynamic inference by estimating parameters from data simulated under a simple epidemiological model. Results indicate that kernel-ABC attained greater accuracy for parameters associated with virus transmission than leading software on the same data sets. Finally, I apply the kernel-ABC framework to study a recent outbreak of a recombinant HIV subtype in China. Kernel-ABC provides a versatile framework for phylodynamic inference because it can fit a broader range of models than methods that rely on the computation of exact likelihoods.

The impact of clinical, demographic and risk factors on rates of HIV transmission: a population-based phylogenetic analysis in British Columbia, Canada.

BACKGROUND: The diversification of human immunodeficiency virus (HIV) is shaped by its transmission history. We therefore used a population based province wide HIV drug resistance database in British Columbia (BC), Canada, to evaluate the impact of clinical, demographic, and behavioral factors on rates of HIV transmission. METHODS: We reconstructed molecular phylogenies from 27,296 anonymized bulk HIV pol sequences representing 7747 individuals in BC-about half the estimated HIV prevalence in BC. Infections were grouped into clusters based on phylogenetic distances, as a proxy for variation in transmission rates. Rates of cluster expansion were reconstructed from estimated dates of HIV seroconversion. RESULTS: Our criteria grouped 4431 individuals into 744 clusters largely separated with respect to risk factors, including large established clusters predominated by injection drug users and more-recently emerging clusters comprising men who have sex with men. The mean log10 viral load of an individual's phylogenetic neighborhood (composed of 5 other individuals with shortest phylogenetic distances) increased their odds of appearing in a cluster by >2-fold per log10 viruses per milliliter. CONCLUSIONS: Hotspots of ongoing HIV transmission can be characterized in near real time by the secondary analysis of HIV resistance genotypes, providing an important potential resource for targeting public health initiatives for HIV prevention.

Theoretical and experimental assessment of degenerate primer tagging in ultra-deep applications of next-generation sequencing.

Primer IDs (pIDs) are random oligonucleotide tags used in next-generation sequencing to identify sequences that originate from the same template. These tags are produced by degenerate primers during the reverse transcription of RNA molecules into cDNA. The use of pIDs helps to track the number of RNA molecules carried through amplification and sequencing, and allows resolution of inconsistencies between reads sharing a pID. Three potential issues complicate the above applications. First, multiple cDNAs may share a pID by chance; we found that while preventing any cDNAs from sharing a pID may be unfeasible, it is still practical to limit the number of these collisions. Secondly, a pID must be observed in at least three sequences to allow error correction; as such, pIDs observed only one or two times must be rejected. If the sequencing product contains copies from a high number of RT templates but produces few reads, our findings indicate that rejecting such pIDs will discard a great deal of data. Thirdly, the use of pIDs could influence amplification and sequencing. We examined the effects of several intrinsic and extrinsic factors on sequencing reads at both the individual and ensemble level.

An evaluation of phylogenetic methods for reconstructing transmitted HIV variants using longitudinal clonal HIV sequence data.

UNLABELLED: A population of human immunodeficiency virus (HIV) within a host often descends from a single transmitted/founder virus. The high mutation rate of HIV, coupled with long delays between infection and diagnosis, make isolating and characterizing this strain a challenge. In theory, ancestral reconstruction could be used to recover this strain from sequences sampled in chronic infection; however, the accuracy of phylogenetic techniques in this context is unknown. To evaluate the accuracy of these methods, we applied ancestral reconstruction to a large panel of published longitudinal clonal and/or single-genome-amplification HIV sequence data sets with at least one intrapatient sequence set sampled within 6 months of infection or seroconversion (n = 19,486 sequences, median [interquartile range] = 49 [20 to 86] sequences/set). The consensus of the earliest sequences was used as the best possible estimate of the transmitted/founder. These sequences were compared to ancestral reconstructions from sequences sampled at later time points using both phylogenetic and phylogeny-naive methods. Overall, phylogenetic methods conferred a 16% improvement in reproducing the consensus of early sequences, compared to phylogeny-naive methods. This relative advantage increased with intrapatient sequence diversity (P < 10(-5)) and the time elapsed between the earliest and subsequent samples (P < 10(-5)). However, neither approach performed well for reconstructing ancestral indel variation, especially within indel-rich regions of the HIV genome. Although further improvements are needed, our results indicate that phylogenetic methods for ancestral reconstruction significantly outperform phylogeny-naive alternatives, and we identify experimental conditions and study designs that can enhance accuracy of transmitted/founder virus reconstruction. IMPORTANCE: When HIV is transmitted into a new host, most of the viruses fail to infect host cells. Consequently, an HIV infection tends to be descended from a single "founder" virus. A priority target for the vaccine research, these transmitted/founder viruses are difficult to isolate since newly infected individuals are often unaware of their status for months or years, by which time the virus population has evolved substantially. Here, we report on the potential use of evolutionary methods to reconstruct the genetic sequence of the transmitted/founder virus from its descendants at later stages of an infection. These methods can recover this ancestral sequence with an overall error rate of about 2.3%-about 15% more information than if we had ignored the evolutionary relationships among viruses. Although there is no substitute for sampling infections at earlier points in time, these methods can provide useful information about the genetic makeup of transmitted/founder HIV.

Use of cellular HIV DNA to predict virologic response to maraviroc: performance of population-based and deep sequencing.

BACKGROUND: A tropism test is required before administration of the antiretroviral drug maraviroc. However, plasma RNA testing is not possible in patients with undetectable plasma viral loads. Here we assess genotypic testing of cellular human immunodeficiency virus (HIV) DNA from peripheral blood mononuclear cells (PBMCs) to predict virologic responses in treatment-experienced patients beginning maraviroc-containing regimens. METHODS: PBMC samples from 181 maraviroc recipients at study entry in MOTIVATE or A4001029 (51% R5 by original Trofile). The V3 loop was amplified in triplicate from cellular HIV DNA, and matching plasma RNA (n = 156). Sequencing was performed using standard population-based methods and next-generation deep sequencing, with tropism assessment as previously defined. RESULTS: Genotypic DNA-based tropism testing from the cellular compartment had 78%-81% sensitivity relative to RNA-based Trofile at the same time point. Cell-based genotypic tropism methods and plasma-based phenotypic and genotypic methods were predictive of virologic response. However, when classifications were discordant, the outcomes favored the plasma predictions over the DNA ones. CONCLUSIONS: Genotypic determination of HIV tropism can be performed using cell-derived viral DNA, and is a predictor of virologic success on maraviroc in therapy-experienced patients. However, the PBMC compartment appears to be a suboptimal predictor compared to plasma.

Dates of HIV infection can be estimated for seroprevalent patients by coalescent analysis of serial next-generation sequencing data.

OBJECTIVE: To reconstruct dates of HIV infection by the coalescent analysis of longitudinal next-generation sequencing (NGS) data. DESIGN: The coalescent predicts the time that has elapsed since the most recent common ancestor (MRCA) of a population. Because HIV tends to undergo severe bottlenecks upon transmission, the MRCA may be a good predictor of the time of infection. NGS provides an efficient means for performing large-scale clonal sequencing of HIV populations within patients, and the ideal raw material for coalescent analysis. METHODS: Baseline and follow-up plasma samples were obtained from 19 individuals enrolled into the Montréal Primary HIV Infection cohort. Dates of infection were initially estimated at baseline from nongenetic data (clinical and serological markers and patient questionnaires). HIV RNA was extracted and seven regions of the genome were amplified, subjected to parallel-tagged 454 pyrosequencing, and analyzed using the software package BEAST. RESULTS: Mean estimates of the time to the MRCA per patient were significantly correlated with nongenetic estimates (Spearman's ρ = 0.65, P = 4.4 × 10(-3)). The median absolute difference between coalescent and nongenetic date estimates was smallest (median 29.4 days) for highly variable regions of the HIV genome such as env V3, and greater (median 114.9 days) for more conserved regions such as pol. CONCLUSION: This application of NGS represents an important advancement, not only because accurate estimates of dates of infection can be derived retrospectively from archived specimens, but also because each analysis is patient-specific and, therefore, robust to variation in rates of HIV evolution.

Immune-driven recombination and loss of control after HIV superinfection.

After acute HIV infection, CD8(+) T cells are able to control viral replication to a set point. This control is often lost after superinfection, although the mechanism behind this remains unclear. In this study, we illustrate in an HLA-B27(+) subject that loss of viral control after HIV superinfection coincides with rapid recombination events within two narrow regions of Gag and Env. Screening for CD8(+) T cell responses revealed that each of these recombination sites (approximately 50 aa) encompassed distinct regions containing two immunodominant CD8 epitopes (B27-KK10 in Gag and Cw1-CL9 in Env). Viral escape and the subsequent development of variant-specific de novo CD8(+) T cell responses against both epitopes were illustrative of the significant immune selection pressures exerted by both responses. Comprehensive analysis of the kinetics of CD8 responses and viral evolution indicated that the recombination events quickly facilitated viral escape from both dominant WT- and variant-specific responses. These data suggest that the ability of a superinfecting strain of HIV to overcome preexisting immune control may be related to its ability to rapidly recombine in critical regions under immune selection pressure. These data also support a role for cellular immune pressures in driving the selection of new recombinant forms of HIV.

Herpes simplex virus type 2 acquisition during recent HIV infection does not influence plasma HIV levels.

We assessed the effect of herpes simplex virus type 2 (HSV-2) acquisition on the plasma HIV RNA and CD4 cell levels among individuals with primary HIV infection using a retrospective cohort analysis. We studied 119 adult, antiretroviral-naive, recently HIV-infected men with a negative HSV-2-specific enzyme immunoassay (EIA) result at enrollment. HSV-2 acquisition was determined by seroconversion on HSV-2 EIA, confirmed by Western blot analysis. Ten men acquired HSV-2 infection a median of 1.3 years after HIV infection (HSV-2 incidence rate of 7.4 per 100 person-years of follow-up). The median time of follow-up after acquiring HSV-2 infection was 303 days. All men except 1 were asymptomatic during HSV-2 acquisition, and only 1 HSV-2 seroconverter, who was asymptomatic, had a transient increase in blood HIV load (0.5 log10 copies/mL over 11 days). The HSV-2 incidence rate was high in our cohort of recently HIV-infected individuals; however, HSV-2 acquisition did not significantly change the plasma HIV dynamics and CD4 cell levels.

An evolutionary-network model reveals stratified interactions in the V3 loop of the HIV-1 envelope.

The third variable loop (V3) of the human immunodeficiency virus type 1 (HIV-1) envelope is a principal determinant of antibody neutralization and progression to AIDS. Although it is undoubtedly an important target for vaccine research, extensive genetic variation in V3 remains an obstacle to the development of an effective vaccine. Comparative methods that exploit the abundance of sequence data can detect interactions between residues of rapidly evolving proteins such as the HIV-1 envelope, revealing biological constraints on their variability. However, previous studies have relied implicitly on two biologically unrealistic assumptions: (1) that founder effects in the evolutionary history of the sequences can be ignored, and; (2) that statistical associations between residues occur exclusively in pairs. We show that comparative methods that neglect the evolutionary history of extant sequences are susceptible to a high rate of false positives (20%-40%). Therefore, we propose a new method to detect interactions that relaxes both of these assumptions. First, we reconstruct the evolutionary history of extant sequences by maximum likelihood, shifting focus from extant sequence variation to the underlying substitution events. Second, we analyze the joint distribution of substitution events among positions in the sequence as a Bayesian graphical model, in which each branch in the phylogeny is a unit of observation. We perform extensive validation of our models using both simulations and a control case of known interactions in HIV-1 protease, and apply this method to detect interactions within V3 from a sample of 1,154 HIV-1 envelope sequences. Our method greatly reduces the number of false positives due to founder effects, while capturing several higher-order interactions among V3 residues. By mapping these interactions to a structural model of the V3 loop, we find that the loop is stratified into distinct evolutionary clusters. We extend our model to detect interactions between the V3 and C4 domains of the HIV-1 envelope, and account for the uncertainty in mapping substitutions to the tree with a parametric bootstrap.

Evolutionary interactions between N-linked glycosylation sites in the HIV-1 envelope.

The addition of asparagine (N)-linked polysaccharide chains (i.e., glycans) to the gp120 and gp41 glycoproteins of human immunodeficiency virus type 1 (HIV-1) envelope is not only required for correct protein folding, but also may provide protection against neutralizing antibodies as a "glycan shield." As a result, strong host-specific selection is frequently associated with codon positions where nonsynonymous substitutions can create or disrupt potential N-linked glycosylation sites (PNGSs). Moreover, empirical data suggest that the individual contribution of PNGSs to the neutralization sensitivity or infectivity of HIV-1 may be critically dependent on the presence or absence of other PNGSs in the envelope sequence. Here we evaluate how glycan-glycan interactions have shaped the evolution of HIV-1 envelope sequences by analyzing the distribution of PNGSs in a large-sequence alignment. Using a "covarion"-type phylogenetic model, we find that the rates at which individual PNGSs are gained or lost vary significantly over time, suggesting that the selective advantage of having a PNGS may depend on the presence or absence of other PNGSs in the sequence. Consequently, we identify specific interactions between PNGSs in the alignment using a new paired-character phylogenetic model of evolution, and a Bayesian graphical model. Despite the fundamental differences between these two methods, several interactions are jointly identified by both. Mapping these interactions onto a structural model of HIV-1 gp120 reveals that negative (exclusive) interactions occur significantly more often between colocalized glycans, while positive (inclusive) interactions are restricted to more distant glycans. Our results imply that the adaptive repertoire of alternative configurations in the HIV-1 glycan shield is limited by functional interactions between the N-linked glycans. This represents a potential vulnerability of rapidly evolving HIV-1 populations that may provide useful glycan-based targets for neutralizing antibodies.

Functional origins of fitness effect-sizes of compensatory mutations in the DNA bacteriophage phiX174.

Epistasis is an important and poorly understood aspect of mutations and strongly influences the evolutionary impact of genetic variation on adaptation and fitness. Although recent studies have begun to characterize the distribution of epistatic effects between mutations affecting fitness, there is currently a lack of empirical information on the underlying biological causes of these epistatic interactions. What are the functional constraints that determine the effectiveness of a compensatory mutation at restoring fitness? We have measured the effect-sizes of 52 compensatory mutations affecting nine different deleterious mutations in the major capsid and spike proteins of the DNA bacteriophage phiX174. On average, an experimentally detectable compensatory mutation recovers about two-thirds of the fitness cost of the preceding deleterious mutation. Variation in fitness effect-sizes is only weakly associated with measures of the distance separating the deleterious and compensatory mutations in the amino acid sequence or the folded protein structure. However, there is a strong association of fitness effect-size with the correlation in the effects of the mutations on the biochemical properties of amino acids. A compensatory mutation has the largest effect-size, on average, when both the compensatory and deleterious mutations have radical effects on the overall biochemical make-up of the amino acids. By examining the relative contributions of specific biochemical properties to variation in fitness effect-size, we find that the area and charge of amino acids have a major influence, which suggests that the complexity of the amino acid phenotype is simplified by selection into a reduced number of phenotypic components.