Viral genetic variation accounts for a third of variability in HIV-1 set-point viral load in Europe

Using viral diversity to identify HIV-1 variants under HLA-dependent selection in a systematic viral genome-wide screen

The pathogenesis of HIV-1 infection is governed by a highly dynamic, time-dependent interaction between the host and the viral genome. In this study, we developed a novel systematic approach to assess the host-virus interaction, using average pairwise viral diversity as a proxy for time since infection, and applied this method to nearly whole viral genome sequences (n = 4,464), human leukocyte antigen (HLA) genotyping data (n = 1,044), and viral RNA load (VL) measurements during the untreated chronic phase (n = 829) of Swiss HIV Cohort Study participants. Our systematic genome-wide screen revealed for 98 HLA/viral-variant pairs a signature of immune-driven selection in the form of an HLA-dependent effect of infection time on the presence of HIV amino acid variants. Of these pairs, 12 were found to have an effect on VL. Furthermore, 28/58 pairs were validated by time-to-event analyses and 48/92 by computational HLA-epitope predictions. Our diversity-based approach allows a powerful and systematic investigation of the interaction between the virus and cellular immunity, revealing a notable subset of such interaction effects. From an evolutionary perspective, these observations underscore the complexity of HLA-mediated selection pressures on the virus that shape viral evolution and pathogenesis.

Self-reported neurocognitive complaints in the Swiss HIV Cohort Study: a viral genome-wide association study

People with HIV may report neurocognitive complaints, with or without associated neurocognitive impairment, varying between individuals and populations. While the HIV genome could play a major role, large systematic viral genome-wide screens to date are lacking. The Swiss HIV Cohort Study biannually enquires neurocognitive complaints. We quantified broad-sense heritability estimates using partial 'pol' sequences from the Swiss HIV Cohort Study resistance database and performed a viral near full-length genome-wide association study for the longitudinal area under the curve of neurocognitive complaints. We performed all analysis (i) restricted to HIV Subtype B and (ii) including all HIV subtypes. From 8547 people with HIV with neurocognitive complaints, we obtained 6966 partial 'pol' sequences and 2334 near full-length HIV sequences. Broad-sense heritability estimates for presence of memory loss complaints ranged between 1% and 17% (Subtype B restricted 1-22%) and increased with the stringency of the phylogenetic distance thresholds. The genome-wide association study revealed one amino acid (Env L641E), after adjusting for multiple testing, positively associated with memory loss complaints (P = 4.3 * 10-6). Other identified mutations, while insignificant after adjusting for multiple testing, were reported in other smaller studies (Tat T64N, Env *291S). We present the first HIV genome-wide association study analysis of neurocognitive complaints and report a first estimate for the heritability of neurocognitive complaints through HIV. Moreover, we could identify one mutation significantly associated with the presence of memory loss complaints. Our findings indicate that neurocognitive complaints are polygenetic and highlight advantages of a whole genome approach for pathogenicity determination.

Accelerating Bayesian inference of dependency between mixed-type biological traits

Inferring dependencies between mixed-type biological traits while accounting for evolutionary relationships between specimens is of great scientific interest yet remains infeasible when trait and specimen counts grow large. The state-of-the-art approach uses a phylogenetic multivariate probit model to accommodate binary and continuous traits via a latent variable framework, and utilizes an efficient bouncy particle sampler (BPS) to tackle the computational bottleneck-integrating many latent variables from a high-dimensional truncated normal distribution. This approach breaks down as the number of specimens grows and fails to reliably characterize conditional dependencies between traits. Here, we propose an inference pipeline for phylogenetic probit models that greatly outperforms BPS. The novelty lies in 1) a combination of the recent Zigzag Hamiltonian Monte Carlo (Zigzag-HMC) with linear-time gradient evaluations and 2) a joint sampling scheme for highly correlated latent variables and correlation matrix elements. In an application exploring HIV-1 evolution from 535 viruses, the inference requires joint sampling from an 11,235-dimensional truncated normal and a 24-dimensional covariance matrix. Our method yields a 5-fold speedup compared to BPS and makes it possible to learn partial correlations between candidate viral mutations and virulence. Computational speedup now enables us to tackle even larger problems: we study the evolution of influenza H1N1 glycosylations on around 900 viruses. For broader applicability, we extend the phylogenetic probit model to incorporate categorical traits, and demonstrate its use to study Aquilegia flower and pollinator co-evolution.

A metagenome-wide association study of HIV disease progression in HIV controllers

Some HIV controllers experience immunologic progression with CD4⁺ T cell decline. We aimed to identify genetic factors associated with CD4⁺ T cell lost in HIV controllers. A total of 561 HIV controllers were included, 442 and 119 from the International HIV controllers Study Cohort and the Swiss HIV Cohort Study, respectively. No SNP or gene was associated with the long-term non-progressor HIV spontaneous control phenotype in the individual GWAS or in the meta-analysis. However, SNPs previously associated with natural HIV control linked to HLA-B (rs2395029 [p = 0.005; OR = 1.70], rs59440261 [p = 0.003; OR = 1.78]), MICA (rs112243036 [p = 0.011; OR = 1.45]), and PSORS1C1 loci (rs3815087 [p = 0.017; OR = 1.39]) showed nominal association with this phenotype. Genetic factors associated with the long-term HIV controllers without risk of immunologic progression are those previously related to the overall HIV controller phenotype.

Deep-sequencing of viral genomes from a large and diverse cohort of treatment-naive HIV-infected persons shows associations between intrahost genetic diversity and viral load

BACKGROUND

Infection with human immunodeficiency virus type 1 (HIV) typically results from transmission of a small and genetically uniform viral population. Following transmission, the virus population becomes more diverse because of recombination and acquired mutations through genetic drift and selection. Viral intrahost genetic diversity remains a major obstacle to the cure of HIV; however, the association between intrahost diversity and disease progression markers has not been investigated in large and diverse cohorts for which the majority of the genome has been deep-sequenced. Viral load (VL) is a key progression marker and understanding of its relationship to viral intrahost genetic diversity could help design future strategies for HIV monitoring and treatment.

METHODS

We analysed deep-sequenced viral genomes from 2,650 treatment-naive HIV-infected persons to measure the intrahost genetic diversity of 2,447 genomic codon positions as calculated by Shannon entropy. We tested for associations between VL and amino acid (AA) entropy accounting for sex, age, race, duration of infection, and HIV population structure.

RESULTS

We confirmed that the intrahost genetic diversity is highest in the env gene. Furthermore, we showed that mean Shannon entropy is significantly associated with VL, especially in infections of >24 months duration. We identified 16 significant associations between VL (p-value<2.0x10-5) and Shannon entropy at AA positions which in our association analysis explained 13% of the variance in VL. Finally, equivalent analysis based on variation in HIV consensus sequences explained only 2% of VL variance.

CONCLUSIONS

Our results elucidate that viral intrahost genetic diversity is associated with VL and could be used as a better disease progression marker than HIV consensus sequence variants, especially in infections of longer duration. We emphasize that viral intrahost diversity should be considered when studying viral genomes and infection outcomes.

TRIAL REGISTRATION

Samples included in this study were derived from participants who consented in the clinical trial, START (NCT00867048) (23), run by the International Network for Strategic Initiatives in Global HIV Trials (INSIGHT). All the participant sites are listed here: http://www.insight-trials.org/start/my_phpscript/participating.php?by=site.

Scalable Bayesian phylogenetics

Recent advances in Bayesian phylogenetics offer substantial computational savings to accommodate increased genomic sampling that challenges traditional inference methods. In this review, we begin with a brief summary of the Bayesian phylogenetic framework, and then conceptualize a variety of methods to improve posterior approximations via Markov chain Monte Carlo (MCMC) sampling. Specifically, we discuss methods to improve the speed of likelihood calculations, reduce MCMC burn-in, and generate better MCMC proposals. We apply several of these techniques to study the evolution of HIV virulence along a 1536-tip phylogeny and estimate the internal node heights of a 1000-tip SARS-CoV-2 phylogenetic tree in order to illustrate the speed-up of such analyses using current state-of-the-art approaches. We conclude our review with a discussion of promising alternatives to MCMC that approximate the phylogenetic posterior. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

A phylogeny-aware GWAS framework to correct for heritable pathogen effects on infectious disease traits

Infectious diseases are particularly challenging for genome-wide association studies (GWAS) because genetic effects from two organisms (pathogen and host) can influence a trait. Traditional GWAS assume individual samples are independent observations. However, pathogen effects on a trait can be heritable from donor to recipient in transmission chains. Thus, residuals in GWAS association tests for host genetic effects may not be independent due to shared pathogen ancestry. We propose a new method to estimate and remove heritable pathogen effects on a trait based on the pathogen phylogeny prior to host GWAS, thus restoring independence of samples. In simulations, we show this additional step can increase GWAS power to detect truly associated host variants when pathogen effects are highly heritable, with strong phylogenetic correlations. We applied our framework to data from two different host-pathogen systems, HIV in humans and X. arboricola in A. thaliana. In both systems, the heritability and thus phylogenetic correlations turn out to be low enough such that qualitative results of GWAS do not change when accounting for the pathogen shared ancestry through a correction step. This means that previous GWAS results applied to these two systems should not be biased due to shared pathogen ancestry. In summary, our framework provides additional information on the evolutionary dynamics of traits in pathogen populations and may improve GWAS if pathogen effects are highly phylogenetically correlated amongst individuals in a cohort.

Attenuated HIV-1 Nef But Not Vpu Function in a Cohort of Rwandan Long-Term Survivors

HIV-1 accessory proteins Nef and Vpu enhance viral pathogenesis through partially overlapping immune evasion activities. Attenuated Nef or Vpu functions have been reported in individuals who display slower disease progression, but few studies have assessed the relative impact of these proteins in non-B HIV-1 subtypes or examined paired proteins from the same individuals. Here, we examined the sequence and function of matched Nef and Vpu clones isolated from 29 long-term survivors (LTS) from Rwanda living with HIV-1 subtype A and compared our results to those of 104 Nef and 62 Vpu clones isolated from individuals living with chronic untreated HIV-1 subtype A from the same geographic area. Nef and vpu coding regions were amplified from plasma HIV RNA and cloned. The function of one intact, phylogenetically-validated Nef and Vpu clone per individual was then quantified by flow cytometry following transient expression in an immortalized CD4+ T-cell line. We measured the ability of each Nef clone to downregulate CD4 and HLA class I, and of each Vpu clone to downregulate CD4 and Tetherin, from the cell surface. Results were normalized to reference clones (Nef-SF2 and Vpu-NL4.3). We observed that Nef-mediated CD4 and HLA downregulation functions were lower in LTS compared to the control cohort (Mann-Whitney p=0.03 and p<0.0001, respectively). Moreover, we found a positive correlation between Nef-mediated CD4 downregulation function and plasma viral load in LTS and controls (Spearman ρ= 0.59, p=0.03 and ρ=0.30, p=0.005, respectively). In contrast, Vpu-mediated functions were similar between groups and did not correlate with clinical markers. Further analyses identified polymorphisms at Nef codon 184 and Vpu codons 60-62 that were associated with function, which were confirmed through mutagenesis. Overall, our results support attenuated function of Nef, but not Vpu, as a contributor to slower disease progression in this cohort of long-term survivors with HIV-1 subtype A.

The interplay between replication capacity of HIV-1 and surrogate markers of disease

BACKGROUND

HIV-1 replication capacity (RC) of transmitted/founder viruses may influence the further course of HIV-1 infection.

METHODS

Replication capacities (RCs) of 355 whole genome primary HIV-1 isolates derived from samples acquired during acute and recent primary HIV-1 infection (PHI) were determined using a novel high throughput infection assay in primary cells. The RCs were used to elucidate potential factors that could be associated with RC during PHI.

RESULTS

Increased RC was found to be associated with increased set point viral load (VL), and significant differences in RCs among 13 different HIV-1 subtypes were discerned. Notably, we observed an increase in RCs for primary HIV-1 isolates of HIV-1 subtype B over a 17-year period. Associations were not observed between RC and CD4 count at sample date of RC measurement, CD4 recovery after initiation of antiretroviral treatment (ART), CD4 decline in untreated individuals, and acute retroviral syndrome severity scores.

DISCUSSION

These findings highlight that RCs of primary HIV-1 isolates acquired during the acute and recent phase of infection are more associated with viral factors, i.e., set point VL, than with host factors. Furthermore, we observed a temporal increase in RC for HIV-1 subtype B viruses over a period of 17 years.

Phylogenetic Estimation of the Viral Fitness Landscape of HIV-1 Set-Point Viral Load

Set-point viral load (SPVL), a common measure of human immunodeficiency virus (HIV)-1 virulence, is partially determined by viral genotype. Epidemiological evidence suggests that this viral property has been under stabilising selection, with a typical optimum for the virus between 10⁴ and 10⁵ copies of viral RNA per ml. Here we aimed to detect transmission fitness differences between viruses from individuals with different SPVLs directly from phylogenetic trees inferred from whole-genome sequences. We used the local branching index (LBI) as a proxy for transmission fitness. We found that LBI is more sensitive to differences in infectiousness than to differences in the duration of the infectious state. By analysing subtype-B samples from the Bridging the Evolution and Epidemiology of HIV in Europe project, we inferred a significant positive relationship between SPVL and LBI up to approximately 10⁵ copies/ml, with some evidence for a peak around this value of SPVL. This is evidence of selection against low values of SPVL in HIV-1 subtype-B strains, likely related to lower infectiousness, and perhaps a peak in the transmission fitness in the expected range of SPVL. The less prominent signatures of selection against higher SPVL could be explained by an inherent limit of the method or the deployment of antiretroviral therapy.

A highly virulent variant of HIV-1 circulating in the Netherlands

We discovered a highly virulent variant of subtype-B HIV-1 in the Netherlands. One hundred nine individuals with this variant had a 0.54 to 0.74 log₁₀ increase (i.e., a ~3.5-fold to 5.5-fold increase) in viral load compared with, and exhibited CD4 cell decline twice as fast as, 6604 individuals with other subtype-B strains. Without treatment, advanced HIV-CD4 cell counts below 350 cells per cubic millimeter, with long-term clinical consequences-is expected to be reached, on average, 9 months after diagnosis for individuals in their thirties with this variant. Age, sex, suspected mode of transmission, and place of birth for the aforementioned 109 individuals were typical for HIV-positive people in the Netherlands, which suggests that the increased virulence is attributable to the viral strain. Genetic sequence analysis suggests that this variant arose in the 1990s from de novo mutation, not recombination, with increased transmissibility and an unfamiliar molecular mechanism of virulence.

Cohort profile: the Netherlands Cohort Study on Acute HIV infection (NOVA), a prospective cohort study of people with acute or early HIV infection who immediately initiate HIV treatment

PURPOSE

Initiation of combination antiretroviral therapy (cART) during acute or early HIV-infection (AEHI) limits the size of the viral reservoir and preserves immune function. This renders individuals who started cART during AEHI promising participants in HIV-cure trials. Therefore, we established a multicentre prospective cohort study in the Netherlands that enrols people with AEHI. In anticipation of future cure trials, we will longitudinally investigate the properties of the viral reservoir size and HIV-specific immune responses among cohort participants.

PARTICIPANTS

Participants immediately initiate intensified cART: dolutegravir, emtricitabine/tenofovir and darunavir/ritonavir (DRV/r). After 4 weeks, once baseline resistance data are available, DRV/r is discontinued. Three study groups are assembled based on the preparedness of individuals to participate in the extensiveness of sampling. Participants accepting immediate treatment and follow-up but declining additional sampling are included in study group 1 ('standard') and routine diagnostic procedures are performed. Participants willing to undergo blood, leukapheresis and semen sampling are included in study group 2 ('less invasive'). In study group 3 ('extended'), additional tissue (gut-associated lymphoid tissue, peripheral lymph node) and cerebrospinal fluid sampling are performed.

FINDINGS TO DATE

Between 2015 and 2020, 140 individuals with AEHI have been enrolled at nine study sites. At enrolment, median age was 36 (IQR 28-47) years, and 134 (95.7%) participants were men. Distribution of Fiebig stages was as follows: Fiebig I, 3 (2.1%); II, 20 (14.3%); III, 7 (5.0%); IV, 49 (35.0%); V, 39 (27.9%); VI, 22 (15.7%). Median plasma HIV RNA was 5.9 (IQR 4.7-6.7) log10 copies/mL and CD4 count 510 (IQR 370-700) cells/mm3. Median time from cART initiation to viral suppression was 8.0 (IQR 4.0-16.0) weeks.

FUTURE PLANS

The Netherlands Cohort Study on Acute HIV infection remains open for participant enrolment and for additional sites to join the network. This cohort provides a unique nationwide platform for conducting future in-depth virological, immunological, host genetic and interventional studies investigating HIV-cure strategies.

Human immunotypes impose selection on viral genotypes through viral epitope specificity

BACKGROUND

Understanding the genetic interplay between human hosts and infectious pathogens is crucial for how we interpret virulence factors. Here, we tested for associations between HIV and host genetics, and interactive genetic effects on viral load (VL) in HIV+ ART-naive clinical trial participants.

METHODS

HIV genomes were sequenced and the encoded amino acid (AA) variants were associated with VL, human single nucleotide polymorphisms (SNPs) and imputed HLA alleles, using generalized linear models with Bonferroni correction.

RESULTS

Human (388,501 SNPs) and HIV (3,010 variants) genetic data was available for 2,122 persons. Four HIV variants were associated with VL (p-values<1.66×10 -5). Twelve HIV variants were associated with a range of 1-512 human SNPs (p-value<4.28×10 -11). We found 46 associations between HLA alleles and HIV variants (p-values<1.29×10 -7). We found HIV variants and immunotypes when analyzed separately, were associated with lower VL, whereas the opposite was true when analyzed in concert. Epitope binding prediction showed HLA alleles to be weaker binders of associated HIV AA variants relative to alternative variants on the same position.

CONCLUSIONS

Our results show the importance of immunotype specificity on viral antigenic determinants, and the identified genetic interplay puts emphasis that viral and human genetics should be studied in the context of each other.

Heritability of the HIV-1 reservoir size and decay under long-term suppressive ART

The HIV-1 reservoir is the major hurdle to curing HIV-1. However, the impact of the viral genome on the HIV-1 reservoir, i.e. its heritability, remains unknown. We investigate the heritability of the HIV-1 reservoir size and its long-term decay by analyzing the distribution of those traits on viral phylogenies from both partial-pol and viral near full-length genome sequences. We use a unique nationwide cohort of 610 well-characterized HIV-1 subtype-B infected individuals on suppressive ART for a median of 5.4 years. We find that a moderate but significant fraction of the HIV-1 reservoir size 1.5 years after the initiation of ART is explained by genetic factors. At the same time, we find more tentative evidence for the heritability of the long-term HIV-1 reservoir decay. Our findings indicate that viral genetic factors contribute to the HIV-1 reservoir size and hence the infecting HIV-1 strain may affect individual patients’ hurdle towards a cure.

The HIV reservoir is a major hurdle for a cure of HIV, but the factors determining its size and dynamics remain unclear. Here the authors show in a large cohort of 610 HIV-1 infected individuals, who are on suppressive ART for a median of 5.4 years, that viral genetic factors contribute substantially to the HIV-1 reservoir size.

Single molecule, near full-length genome sequencing of dengue virus

Current methods for dengue virus (DENV) genome amplification, amplify parts of the genome in at least 5 overlapping segments and then combine the output to characterize a full genome. This process is laborious, costly and requires at least 10 primers per serotype, thus increasing the likelihood of PCR bias. We introduce an assay to amplify near full-length dengue virus genomes as intact molecules, sequence these amplicons with third generation “nanopore” technology without fragmenting and use the sequence data to differentiate within-host viral variants with a bioinformatics tool (Nano-Q). The new assay successfully generated near full-length amplicons from DENV serotypes 1, 2 and 3 samples which were sequenced with nanopore technology. Consensus DENV sequences generated by nanopore sequencing had over 99.5% pairwise sequence similarity to Illumina generated counterparts provided the coverage was > 100 with both platforms. Maximum likelihood phylogenetic trees generated from nanopore consensus sequences were able to reproduce the exact trees made from Illumina sequencing with a conservative 99% bootstrapping threshold (after 1000 replicates and 10% burn-in). Pairwise genetic distances of within host variants identified from the Nano-Q tool were less than that of between host variants, thus enabling the phylogenetic segregation of variants from the same host.

Bimodal distribution and set point HBV DNA viral loads in chronic infection: retrospective analysis of cohorts from the UK and South Africa

Hepatitis B virus (HBV) viral load (VL) is used as a biomarker to assess risk of disease progression, and to determine eligibility for treatment. While there is a well recognised association between VL and the expression of the viral e-antigen protein, the distributions of VL at a population level are not well described. We here present cross-sectional, observational HBV VL data from two large population cohorts in the UK and in South Africa, demonstrating a consistent bimodal distribution. The right skewed distribution and low median viral loads are different from the left-skew and higher viraemia in seen in HIV and hepatitis C virus (HCV) cohorts in the same settings. Using longitudinal data, we present evidence for a stable 'set-point' VL in peripheral blood during chronic HBV infection. These results are important to underpin improved understanding of HBV biology, to inform approaches to viral sequencing, and to plan public health interventions.

D. de Azevedo S, Delatorre E, Bello G. Few amino acid signatures distinguish HIV-1 subtype B pandemic and non-pandemic strains

The Human Immunodeficiency Virus Type I (HIV-1) subtype B comprises approximately 10% of all HIV infections in the world. The HIV-1 subtype B epidemic comprehends a pandemic variant (named B_PANDEMIC) disseminated worldwide and non-pandemic variants (named B_CAR) that are mostly restricted to the Caribbean. The goal of this work was the identification of amino acid signatures (AAs) characteristic to the B_CAR and B_PANDEMIC variants. To this end, we analyzed HIV-1 subtype B full-length (n = 486) and partial (n = 814) genomic sequences from the Americas classified within the B_CAR and B_PANDEMIC clades and reconstructed the sequences of their most recent common ancestors (MRCA). Analysis of contemporary HIV-1 sequences revealed 13 AAs between B_CAR and B_PANDEMIC variants (four on Gag, three on Pol, three on Rev, and one in Vif, Vpu, and Tat) of which only two (one on Gag and one on Pol) were traced to the MRCA. All AAs correspond to polymorphic sites located outside essential functional proteins domains, except the AAs in Tat. The absence of stringent AAs inherited from their ancestors between modern B_CAR and B_PANDEMIC variants support that ecological factors, rather than viral determinants, were the main driving force behind the successful spread of the B_PANDEMIC strain.

A Comprehensive Genomics Solution for HIV Surveillance and Clinical Monitoring in Low-Income Settings

Viral genetic sequencing can be used to monitor the spread of HIV drug resistance, identify appropriate antiretroviral regimes, and characterize transmission dynamics. Despite decreasing costs, next-generation sequencing (NGS) is still prohibitively costly for routine use in generalized HIV epidemics in low- and middle-income countries. Here, we present veSEQ-HIV, a high-throughput, cost-effective NGS sequencing method and computational pipeline tailored specifically to HIV, which can be performed using leftover blood drawn for routine CD4 cell count testing. This method overcomes several major technical challenges that have prevented HIV sequencing from being used routinely in public health efforts; it is fast, robust, and cost-efficient, and generates full genomic sequences of diverse strains of HIV without bias. The complete veSEQ-HIV pipeline provides viral load estimates and quantitative summaries of drug resistance mutations; it also exploits information on within-host viral diversity to construct directed transmission networks. We evaluated the method's performance using 1,620 plasma samples collected from individuals attending 10 large urban clinics in Zambia as part of the HPTN 071-2 study (PopART Phylogenetics). Whole HIV genomes were recovered from 91% of samples with a viral load of >1,000 copies/ml. The cost of the assay (30 GBP per sample) compares favorably with existing VL and HIV genotyping tests, proving an affordable option for combining HIV clinical monitoring with molecular epidemiology and drug resistance surveillance in low-income settings.

Bimodal distribution and set point HBV DNA viral loads in chronic infection: retrospective analysis of cohorts from the UK and South Africa

Hepatitis B virus (HBV) viral load (VL) is used as a biomarker to assess risk of disease progression, and to determine eligibility for treatment. While there is a well recognised association between VL and the expression of the viral e-antigen (HBeAg) protein, the precise determinants and distribution of VL at a population level are not well described. We here report the distribution of HBV VL in two large cross-sectional population cohorts in the UK and in South Africa, demonstrating a consistent bimodal distribution. The right skewed distribution and low median viral loads are significantly different from the left-skew and higher viraemia in seen in comparable HIV and hepatitis C virus (HCV) cohorts. Using longitudinal data, we present evidence for a stable 'set-point' VL in peripheral blood during chronic HBV infection. These results are important to underpin improved understanding of HBV biology and to plan public health interventions.

A Systematic Phylogenetic Approach to Study the Interaction of HIV-1 With Coinfections, Noncommunicable Diseases, and Opportunistic Diseases

To systematically test whether coinfections spread along the HIV-1 transmission network and whether similarities in HIV-1 genomes predict AIDS-defining illnesses and comorbidities, we analyzed the distribution of these variables on the HIV phylogeny of the densely sampled Swiss HIV Cohort Study. By combining different statistical methods, we could detect, quantify, and explain the clustering of diseases. Infectious conditions such as hepatitis C, but also Kaposi sarcoma, clustered significantly, suggesting transmission of these infections along the HIV-1 transmission network. The clustering of patients with neurocognitive complaints could not be completely explained by the clustering of patients with similar demographic risk factors, which suggests a potential impact of viral genetics. In summary, the consistent and robust signal for coinfections and comorbidities highlights the strong interaction of HIV-1 and other infections and shows the potential of combining phylogenetic methods to identify disease traits that are likely to be related to virus genetic factors.

Natural selection favoring more transmissible HIV detected in United States molecular transmission network

HIV molecular epidemiology can identify clusters of individuals with elevated rates of HIV transmission. These variable transmission rates are primarily driven by host risk behavior; however, the effect of viral traits on variable transmission rates is poorly understood. Viral load, the concentration of HIV in blood, is a heritable viral trait that influences HIV infectiousness and disease progression. Here, we reconstruct HIV genetic transmission clusters using data from the United States National HIV Surveillance System and report that viruses in clusters, inferred to be frequently transmitted, have higher viral loads at diagnosis. Further, viral load is higher in people in larger clusters and with increased network connectivity, suggesting that HIV in the United States is experiencing natural selection to be more infectious and virulent. We also observe a concurrent increase in viral load at diagnosis over the last decade. This evolutionary trajectory may be slowed by prevention strategies prioritized toward rapidly growing transmission clusters.

New HIV-1 circulating recombinant form 94: from phylogenetic detection of a large transmission cluster to prevention in the age of geosocial-networking apps in France, 2013 to 2017

BackgroundEnding the HIV pandemic must involve new tools to rapidly identify and control local outbreaks and prevent the emergence of recombinant strains with epidemiological advantages.AimThis observational study aimed to investigate in France a cluster of HIV-1 cases related to a new circulating recombinant form (CRF). The confirmation this CRF's novelty as well as measures to control its spread are presented.MethodsPhylogenetic analyses of HIV sequences routinely generated for drug resistance genotyping before 2018 in French laboratories were employed to detect the transmission chain. The CRF involved was characterised by almost full-length viral sequencing for six cases. Cases' clinical data were reviewed. Where possible, epidemiological information was collected with a questionnaire.ResultsThe transmission cluster comprised 49 cases, mostly diagnosed in 2016-2017 (n = 37). All were infected with a new CRF, CRF94_cpx. The molecular proximity of this CRF to X4 strains and the high median viraemia, exceeding 5.0 log₁₀ copies/mL, at diagnosis, even in chronic infection, raise concerns of enhanced virulence. Overall, 41 cases were diagnosed in the Ile-de-France region and 45 were men who have sex with men. Among 24 cases with available information, 20 reported finding partners through a geosocial networking app. Prevention activities in the area and population affected were undertaken.ConclusionWe advocate the systematic use of routinely generated HIV molecular data by a dedicated reactive network, to improve and accelerate targeted prevention interventions. Geosocial networking apps can play a role in the spread of outbreaks, but could also deliver local targeted preventive alerts.

Automatic generation of evolutionary hypotheses using mixed Gaussian phylogenetic models

Phylogenetic comparative methods are widely used to understand and quantify the evolution of phenotypic traits, based on phylogenetic trees and trait measurements of extant species. Such analyses depend crucially on the underlying model. Gaussian phylogenetic models like Brownian motion and Ornstein-Uhlenbeck processes are the workhorses of modeling continuous-trait evolution. However, these models fit poorly to big trees, because they neglect the heterogeneity of the evolutionary process in different lineages of the tree. Previous works have addressed this issue by introducing shifts in the evolutionary model occurring at inferred points in the tree. However, for computational reasons, in all current implementations, these shifts are "intramodel," meaning that they allow jumps in 1 or 2 model parameters, keeping all other parameters "global" for the entire tree. There is no biological reason to restrict a shift to a single model parameter or, even, to a single type of model. Mixed Gaussian phylogenetic models (MGPMs) incorporate the idea of jointly inferring different types of Gaussian models associated with different parts of the tree. Here, we propose an approximate maximum-likelihood method for fitting MGPMs to comparative data comprising possibly incomplete measurements for several traits from extant and extinct phylogenetically linked species. We applied the method to the largest published tree of mammal species with body- and brain-mass measurements, showing strong statistical support for an MGPM with 12 distinct evolutionary regimes. Based on this result, we state a hypothesis for the evolution of the brain-body-mass allometry over the past 160 million y.

Dissecting HIV Virulence: Heritability of Setpoint Viral Load, CD4+ T-Cell Decline, and Per-Parasite Pathogenicity

Pathogen strains may differ in virulence because they attain different loads in their hosts, or because they induce different disease-causing mechanisms independent of their load. In evolutionary ecology, the latter is referred to as “per-parasite pathogenicity”. Using viral load and CD4+ T-cell measures from 2014 HIV-1 subtype B-infected individuals enrolled in the Swiss HIV Cohort Study, we investigated if virulence—measured as the rate of decline of CD4+ T cells—and per-parasite pathogenicity are heritable from donor to recipient. We estimated heritability by donor–recipient regressions applied to 196 previously identified transmission pairs, and by phylogenetic mixed models applied to a phylogenetic tree inferred from HIV pol sequences. Regressing the CD4+ T-cell declines and per-parasite pathogenicities of the transmission pairs did not yield heritability estimates significantly different from zero. With the phylogenetic mixed model, however, our best estimate for the heritability of the CD4+ T-cell decline is 17% (5–30%), and that of the per-parasite pathogenicity is 17% (4–29%). Further, we confirm that the set-point viral load is heritable, and estimate a heritability of 29% (12–46%). Interestingly, the pattern of evolution of all these traits differs significantly from neutrality, and is most consistent with stabilizing selection for the set-point viral load, and with directional selection for the CD4+ T-cell decline and the per-parasite pathogenicity. Our analysis shows that the viral genotype affects virulence mainly by modulating the per-parasite pathogenicity, while the indirect effect via the set-point viral load is minor.

The phylogenomics of evolving virus virulence

How virulence evolves after a virus jumps to a new host species is central to disease emergence. Our current understanding of virulence evolution is based on insights drawn from two perspectives that have developed largely independently: long-standing evolutionary theory based on limited real data examples that often lack a genomic basis, and experimental studies of virulence-determining mutations using cell culture or animal models. A more comprehensive understanding of virulence mutations and their evolution can be achieved by bridging the gap between these two research pathways through the phylogenomic analysis of virus genome sequence data as a guide to experimental study.

Easy and accurate reconstruction of whole HIV genomes from short-read sequence data with shiver

Studying the evolution of viruses and their molecular epidemiology relies on accurate viral sequence data, so that small differences between similar viruses can be meaningfully interpreted. Despite its higher throughput and more detailed minority variant data, next-generation sequencing has yet to be widely adopted for HIV. The difficulty of accurately reconstructing the consensus sequence of a quasispecies from reads (short fragments of DNA) in the presence of large between- and within-host diversity, including frequent indels, may have presented a barrier. In particular, mapping (aligning) reads to a reference sequence leads to biased loss of information; this bias can distort epidemiological and evolutionary conclusions. De novo assembly avoids this bias by aligning the reads to themselves, producing a set of sequences called contigs. However contigs provide only a partial summary of the reads, misassembly may result in their having an incorrect structure, and no information is available at parts of the genome where contigs could not be assembled. To address these problems we developed the tool shiver to pre-process reads for quality and contamination, then map them to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences. Run with two commands per sample, it can easily be used for large heterogeneous data sets. We used shiver to reconstruct the consensus sequence and minority variant information from paired-end short-read whole-genome data produced with the Illumina platform, for sixty-five existing publicly available samples and fifty new samples. We show the systematic superiority of mapping to shiver's constructed reference compared with mapping the same reads to the closest of 3,249 real references: median values of 13 bases called differently and more accurately, 0 bases called differently and less accurately, and 205 bases of missing sequence recovered. We also successfully applied shiver to whole-genome samples of Hepatitis C Virus and Respiratory Syncytial Virus. shiver is publicly available from https://github.com/ChrisHIV/shiver.

The evolution of subtype B HIV-1 tat in the Netherlands during 1985-2012

For the production of viral genomic RNA, HIV-1 is dependent on an early viral protein, Tat, which is required for high-level transcription. The quantity of viral RNA detectable in blood of HIV-1 infected individuals varies dramatically, and a factor involved could be the efficiency of Tat protein variants to stimulate RNA transcription. HIV-1 virulence, measured by set-point viral load, has been observed to increase over time in the Netherlands and elsewhere. Investigation of tat gene evolution in clinical isolates could discover a role of Tat in this changing virulence. A dataset of 291 Dutch HIV-1 subtype B tat genes, derived from full-length HIV-1 genome sequences from samples obtained between 1985-2012, was used to analyse the evolution of Tat. Twenty-two patient-derived tat genes, and the control Tat_HXB2 were analysed for their capacity to stimulate expression of an LTR-luciferase reporter gene construct in diverse cell lines, as well as for their ability to complement a tat-defective HIV-1_LAI clone. Analysis of 291 historical tat sequences from the Netherlands showed ample amino acid (aa) variation between isolates, although no specific mutations were selected for over time. Of note, however, the encoded protein varied its length over the years through the loss or gain of stop codons in the second exon. In transmission clusters, a selection against the shorter Tat86 ORF was apparent in favour of the more common Tat101 version, likely due to negative selection against Tat86 itself, although random drift, transmission bottlenecks, or linkage to other variants could also explain the observation. There was no correlation between Tat length and set-point viral load; however, the number of non-intermediate variants in our study was small. In addition, variation in the length of Tat did not significantly change its capacity to stimulate transcription. From 1985 till 2012, variation in the length of the HIV-1 subtype B tat gene is increasingly found in the Dutch epidemic. However, as Tat proteins did not differ significantly in their capacity to stimulate transcription elongation in vitro, the increased HIV-1 virulence seen in recent years could not be linked to an evolving viral Tat protein.

Viral Characteristics Associated with the Clinical Nonprogressor Phenotype Are Inherited by Viruses from a Cluster of HIV-1 Elite Controllers

A small group of HIV-1-infected individuals, called long-term nonprogressors (LTNPs), and in particular a subgroup of LTNPs, elite controllers (LTNP-ECs), display permanent control of viral replication and lack of clinical progression. This control is the result of a complex interaction of host, immune, and viral factors. We identified, by phylogenetic analysis, a cluster of LTNP-ECs infected with very similar low-replication HIV-1 viruses, suggesting the contribution of common viral features to the clinical LTNP-EC phenotype. HIV-1 envelope (Env) glycoprotein mediates signaling and promotes HIV-1 fusion, entry, and infection, being a key factor of viral fitness in vitro, cytopathicity, and infection progression in vivo Therefore, we isolated full-length env genes from viruses of these patients and from chronically infected control individuals. Functional characterization of the initial events of the viral infection showed that Envs from the LTNP-ECs were ineffective in the binding to CD4 and in the key triggering of actin/tubulin-cytoskeleton modifications compared to Envs from chronic patients. The viral properties of the cluster viruses result in a defective viral fusion, entry, and infection, and these properties were inherited by every virus of the cluster. Therefore, inefficient HIV-1 Env functions and signaling defects may contribute to the low viral replication capacity and transmissibility of the cluster viruses, suggesting a direct role in the LTNP-EC phenotype of these individuals. These results highlight the important role of viral characteristics in the LTNP-EC clinical phenotype. These Env viral properties were common to all the cluster viruses and thus support the heritability of the viral characteristics.IMPORTANCE HIV-1 long-term nonprogressor elite controller patients, due to their permanent control of viral replication, have been the object of numerous studies to identify the factors responsible for this clinical phenotype. In this work, we analyzed the viral characteristics of the envelopes of viruses from a phylogenetic cluster of LTNP-EC patients. These envelopes showed ineffective binding to CD4 and the subsequent signaling activity to modify actin/tubulin cytoskeletons, which result in low fusion and deficient entry and infection capacities. These Env viral characteristics could explain the nonprogressor clinical phenotype of these patients. In addition, these inefficient env viral properties were present in all viruses of the cluster, supporting the heritability of the viral phenotype.

Reconciling Pasteur and Darwin to control infectious diseases

The continual emergence of new pathogens and the increased spread of antibiotic resistance in bacterial populations remind us that microbes are living entities that evolve at rates that impact public health interventions. Following the historical thread of the works of Pasteur and Darwin shows how reconciling clinical microbiology, ecology, and evolution can be instrumental to understanding pathology, developing new therapies, and prolonging the efficiency of existing ones.

A Practical Guide to Estimating the Heritability of Pathogen Traits

Pathogen traits, such as the virulence of an infection, can vary significantly between patients. A major challenge is to measure the extent to which genetic differences between infecting strains explain the observed variation of the trait. This is quantified by the trait’s broad-sense heritability, H². A recent discrepancy between estimates of the heritability of HIV-virulence has opened a debate on the estimators’ accuracy. Here, we show that the discrepancy originates from model limitations and important lifecycle differences between sexually reproducing organisms and transmittable pathogens. In particular, current quantitative genetics methods, such as donor–recipient regression of surveyed serodiscordant couples and the phylogenetic mixed model (PMM), are prone to underestimate H², because they neglect or do not fit to the loss of resemblance between transmission partners caused by within-host evolution. In a phylogenetic analysis of 8,483 HIV patients from the United Kingdom, we show that the phenotypic correlation between transmission partners decays with the amount of within-host evolution of the virus. We reproduce this pattern in toy-model simulations and show that a phylogenetic Ornstein–Uhlenbeck model (POUMM) outperforms the PMM in capturing this correlation pattern and in quantifying H². In particular, we show that POUMM outperforms PMM even in simulations without selection—as it captures the mentioned correlation pattern—which has not been appreciated until now. By cross-validating the POUMM estimates with ANOVA on closest phylogenetic pairs, we obtain H² ≈ 0.2, meaning ∼20% of the variation in HIV-virulence is explained by the virus genome both for European and African data.