Immuno-epidemiological modeling of HIV-1 predicts high heritability of the set-point virus load, while selection for CTL escape dominates virulence evolution

Mathematical methods for scaling from within-host to population-scale in infectious disease systems

Mathematical modellers model infectious disease dynamics at different scales. Within-host models represent the spread of pathogens inside an individual, whilst between-host models track transmission between individuals. However, pathogen dynamics at one scale affect those at another. This has led to the development of multiscale models that connect within-host and between-host dynamics. In this article, we systematically review the literature on multiscale infectious disease modelling according to PRISMA guidelines, dividing previously published models into five categories governing their methodological approaches (Garira (2017)), explaining their benefits and limitations. We provide a primer on developing multiscale models of infectious diseases.

Reconstructing multi-strain pathogen interactions from cross-sectional survey data via statistical network inference

Infectious diseases often involve multiple pathogen species or multiple strains of the same pathogen. As such, knowledge of how different pathogens interact is key to understand and predict the outcome of interventions targeting only a subset of species or strains involved in disease. Population-level data may be useful to infer pathogen strain interactions, but most previously used inference methods only consider uniform interactions between all strains or focus on marginal pairwise interactions. As such, these methods are prone to bias induced by indirect interactions through other strains. Here, we evaluated statistical network inference for reconstructing heterogeneous interactions from cross-sectional surveys detecting joint presence/absence patterns of pathogen strains within hosts. We applied various network models to simulated survey data, representing endemic infection states of multiple pathogen strains with potential interactions in acquisition or clearance of infection. Satisfactory performance was demonstrated by the estimators converging to the true interactions. Accurate reconstruction of interaction networks was achieved by regularization or penalization for sample size. Although performance deteriorated in the presence of host heterogeneity, this was overcome by correcting for individual-level risk factors. Our work demonstrates how statistical network inference could prove useful for detecting multi-strain pathogen interactions and may have applications beyond epidemiology.

Key Factors and Parameter Ranges for Immune Control of Equine Infectious Anemia Virus Infection

Equine Infectious Anemia Virus (EIAV) is an important infection in equids, and its similarity to HIV creates hope for a potential vaccine. We analyze a within-host model of EIAV infection with antibody and cytotoxic T lymphocyte (CTL) responses. In this model, the stability of the biologically relevant endemic equilibrium, characterized by the coexistence of long-term antibody and CTL levels, relies upon a balance between CTL and antibody growth rates, which is needed to ensure persistent CTL levels. We determine the model parameter ranges at which CTL and antibody proliferation rates are simultaneously most influential in leading the system towards coexistence and can be used to derive a mathematical relationship between CTL and antibody production rates to explore the bifurcation curve that leads to coexistence. We employ Latin hypercube sampling and least squares to find the parameter ranges that equally divide the endemic and boundary equilibria. We then examine this relationship numerically via a local sensitivity analysis of the parameters. Our analysis is consistent with previous results showing that an intervention (such as a vaccine) intended to control a persistent viral infection with both immune responses should moderate the antibody response to allow for stimulation of the CTL response. Finally, we show that the CTL production rate can entirely determine the long-term outcome, regardless of the effect of other parameters, and we provide the conditions for this result in terms of the identified ranges for all model parameters.

Patient characteristics and determinants of CD4 at diagnosis of HIV in Mexico from 2008 to 2017: a 10-year population-based study

BACKGROUND

In 2007-2012 the Mexican government launched the National HIV program and there was a major change in HIV policies implemented in 2013-2018, when efforts focused on prevention, increase in early diagnosis and timely treatment. Still, late HIV diagnosis is a major concern in Mexico due to its association with the development of AIDS development and mortality. Thus, the objectives of this study were to identify the determinants of late HIV diagnosis (i.e. CD4 count less than 200 cells/mm3) in Mexico from 2008 to 2017 and to evaluate the impact of the 2013-2017 National HIV program.

METHODS

Using patient level data from the SALVAR database, which includes 64% of the population receiving HIV care in Mexico, an adjusted logistic model was conducted. Main study outcomes were HIV late diagnosis which was defined as CD4 count less than 200 cells/mm3 at diagnosis.

RESULTS

The study included 106,830 individuals newly diagnosed with HIV and treated in Mexican public health facilities between 2008 and 2017 (mean age: 33 years old, 80% male). HIV late diagnosis decreased from 45 to 43% (P < 0.001) between 2008 and 2012 and 2013-2017 (i.e. before and after the implementation of the 2013-2017 policy). Multivariable logistic regressions indicated that being diagnosed between 2013 and 2017 (odds ratio [OR] = 0.96 [95% Confidence interval [CI] [0.93, 0.98]) or in health facilities specialized in HIV care (OR = 0.64 [95% CI 0.60, 0.69]) was associated with early diagnosis. Being male, older than 29 years old, diagnosed in Central East, the South region of Mexico or in high-marginalized locality increased the odds of a late diagnosis.

CONCLUSIONS

The results of this study indicate that the 2013-2017 National HIV program in Mexico has been marginally successful in decreasing the proportion of individuals with late HIV diagnosis in Mexico. We identified several predictors of late diagnosis which could help establishing health policies. The main determinants for late diagnosis were being male, older than 29 years old, and being diagnosed in a Hospital or National Institute.

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.

Linked within-host and between-host models and data for infectious diseases: a systematic review

The observed dynamics of infectious diseases are driven by processes across multiple scales. Here we focus on two: within-host, that is, how an infection progresses inside a single individual (for instance viral and immune dynamics), and between-host, that is, how the infection is transmitted between multiple individuals of a host population. The dynamics of each of these may be influenced by the other, particularly across evolutionary time. Thus understanding each of these scales, and the links between them, is necessary for a holistic understanding of the spread of infectious diseases. One approach to combining these scales is through mathematical modeling. We conducted a systematic review of the published literature on multi-scale mathematical models of disease transmission (as defined by combining within-host and between-host scales) to determine the extent to which mathematical models are being used to understand across-scale transmission, and the extent to which these models are being confronted with data. Following the PRISMA guidelines for systematic reviews, we identified 24 of 197 qualifying papers across 30 years that include both linked models at the within and between host scales and that used data to parameterize/calibrate models. We find that the approach that incorporates both modeling with data is under-utilized, if increasing. This highlights the need for better communication and collaboration between modelers and empiricists to build well-calibrated models that both improve understanding and may be used for prediction.

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.

Evolution of HIV-1 within untreated individuals and at the population scale in Uganda

HIV-1 undergoes multiple rounds of error-prone replication between transmission events, resulting in diverse viral populations within and among individuals. In addition, the virus experiences different selective pressures at multiple levels: during the course of infection, at transmission, and among individuals. Disentangling how these evolutionary forces shape the evolution of the virus at the population scale is important for understanding pathogenesis, how drug- and immune-escape variants are likely to spread in populations, and the development of preventive vaccines. To address this, we deep-sequenced two regions of the HIV-1 genome (p24 and gp41) from 34 longitudinally-sampled untreated individuals from Rakai District in Uganda, infected with subtypes A, D, and inter-subtype recombinants. This dataset substantially increases the availability of HIV-1 sequence data that spans multiple years of untreated infection, in particular for different geographical regions and viral subtypes. In line with previous studies, we estimated an approximately five-fold faster rate of evolution at the within-host compared to the population scale for both synonymous and nonsynonymous substitutions, and for all subtypes. We determined the extent to which this mismatch in evolutionary rates can be explained by the evolution of the virus towards population-level consensus, or the transmission of viruses similar to those that establish infection within individuals. Our findings indicate that both processes are likely to be important.

The impact of HIV-1 within-host evolution on transmission dynamics

The adaptive potential of HIV-1 is a vital mechanism to evade host immune responses and antiviral treatment. However, high evolutionary rates during persistent infection can impair transmission efficiency and alter disease progression in the new host, resulting in a delicate trade-off between within-host virulence and between-host infectiousness. This trade-off is visible in the disparity in evolutionary rates at within-host and between-host levels, and preferential transmission of ancestral donor viruses. Understanding the impact of within-host evolution for epidemiological studies is essential for the design of preventive and therapeutic measures. Herein, we review recent theoretical and experimental work that generated new insights into the complex link between within-host evolution and between-host fitness, revealing temporal and selective processes underlying the structure and dynamics of HIV-1 transmission.

How Germinal Centers Evolve Broadly Neutralizing Antibodies: the Breadth of the Follicular Helper T Cell Response

Many HIV-1-infected patients evolve broadly neutralizing antibodies (bnAbs). This evolutionary process typically takes several years and is poorly understood as selection taking place in germinal centers occurs on the basis of antibody affinity. B cells with the highest-affinity receptors tend to acquire the most antigen from the follicular dendritic cell (FDC) network and present the highest density of cognate peptides to follicular helper T (Tfh) cells, which provide survival signals to the B cell. bnAbs are therefore expected to evolve only when the B cell lineage evolving breadth is consistently capturing and presenting more peptides to Tfh cells than other lineages of more specific B cells. Here we develop mathematical models of Tfh cells in germinal centers to explicitly define the mechanisms of selection in this complex evolutionary process. Our results suggest that broadly reactive B cells presenting a high density of peptides bound to major histocompatibility complex class II molecules (pMHC) are readily outcompeted by B cells responding to lineages of HIV-1 that transiently dominate the within host viral population. Conversely, if broadly reactive B cells acquire a large variety of several HIV-1 proteins from the FDC network and present a high diversity of several pMHC, they can be rescued by a large fraction of the Tfh cell repertoire in the germinal center. Under such circumstances the evolution of bnAbs is much more consistent. Increasing either the magnitude of the Tfh cell response or the breadth of the Tfh cell repertoire markedly facilitates the evolution of bnAbs. Because both the magnitude and breadth can be increased by vaccination with several HIV-1 proteins, this calls for experimental testing.IMPORTANCE Many HIV-infected patients slowly evolve antibodies that can neutralize a large variety of viruses. Such broadly neutralizing antibodies (bnAbs) could in the future become therapeutic agents. bnAbs appear very late, and patients are typically not protected by them. At the moment, we fail to understand why this takes so long and how the immune system selects for broadly neutralizing capacity. Typically, antibodies are selected based on affinity and not on breadth. We developed mathematical models to study two different mechanisms by which the immune system can select for broadly neutralizing capacity. One of these is based upon the repertoire of different follicular helper T (Tfh) cells in germinal centers. We suggest that broadly reactive B cells may interact with a larger fraction of this repertoire and demonstrate that this would select for bnAbs. Intriguingly, this suggests that broadening the Tfh cell repertoire by vaccination may speed up the evolution of bnAbs.

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

HIV-1 set-point viral load-the approximately stable value of viraemia in the first years of chronic infection-is a strong predictor of clinical outcome and is highly variable across infected individuals. To better understand HIV-1 pathogenesis and the evolution of the viral population, we must quantify the heritability of set-point viral load, which is the fraction of variation in this phenotype attributable to viral genetic variation. However, current estimates of heritability vary widely, from 6% to 59%. Here we used a dataset of 2,028 seroconverters infected between 1985 and 2013 from 5 European countries (Belgium, Switzerland, France, the Netherlands and the United Kingdom) and estimated the heritability of set-point viral load at 31% (CI 15%-43%). Specifically, heritability was measured using models of character evolution describing how viral load evolves on the phylogeny of whole-genome viral sequences. In contrast to previous studies, (i) we measured viral loads using standardized assays on a sample collected in a strict time window of 6 to 24 months after infection, from which the viral genome was also sequenced; (ii) we compared 2 models of character evolution, the classical "Brownian motion" model and another model ("Ornstein-Uhlenbeck") that includes stabilising selection on viral load; (iii) we controlled for covariates, including age and sex, which may inflate estimates of heritability; and (iv) we developed a goodness of fit test based on the correlation of viral loads in cherries of the phylogenetic tree, showing that both models of character evolution fit the data well. An overall heritability of 31% (CI 15%-43%) is consistent with other studies based on regression of viral load in donor-recipient pairs. Thus, about a third of variation in HIV-1 virulence is attributable to viral genetic variation.

Estimating the Respective Contributions of Human and Viral Genetic Variation to HIV Control

We evaluated the fraction of variation in HIV-1 set point viral load attributable to viral or human genetic factors by using joint host/pathogen genetic data from 541 HIV infected individuals. We show that viral genetic diversity explains 29% of the variation in viral load while host factors explain 8.4%. Using a joint model including both host and viral effects, we estimate a total of 30% heritability, indicating that most of the host effects are reflected in viral sequence variation.

Viral loads of Human Immunodeficiency Virus infections are correlated between the donor and the recipient of the transmission pair. Similarly, human genetic factors may modulate viral load. In this study we estimate the extents to which viral load is heritable either via the viral genotype (from donor to recipient) or via the host’s Human Leukocyte Antigen (HLA) genotype. We find that a major fraction of inter individual variability is explained by the similarity of the viral genotypes, and that human genetic variation in the HLA region provide little additional explanatory power.

Effect of the Latent Reservoir on the Evolution of HIV at the Within- and Between-Host Levels

The existence of long-lived reservoirs of latently infected CD4+ T cells is the major barrier to curing HIV, and has been extensively studied in this light. However, the effect of these reservoirs on the evolutionary dynamics of the virus has received little attention. Here, we present a within-host quasispecies model that incorporates a long-lived reservoir, which we then nest into an epidemiological model of HIV dynamics. For biologically plausible parameter values, we find that the presence of a latent reservoir can severely delay evolutionary dynamics within a single host, with longer delays associated with larger relative reservoir sizes and/or homeostatic proliferation of cells within the reservoir. These delays can fundamentally change the dynamics of the virus at the epidemiological scale. In particular, the delay in within-host evolutionary dynamics can be sufficient for the virus to evolve intermediate viral loads consistent with maximising transmission, as is observed, and not the very high viral loads that previous models have predicted, an effect that can be further enhanced if viruses similar to those that initiate infection are preferentially transmitted. These results depend strongly on within-host characteristics such as the relative reservoir size, with the evolution of intermediate viral loads observed only when the within-host dynamics are sufficiently delayed. In conclusion, we argue that the latent reservoir has important, and hitherto under-appreciated, roles in both within- and between-host viral evolution.

The impact of antiretroviral therapy on population-level virulence evolution of HIV-1

In HIV-infected patients, an individual's set point viral load (SPVL) strongly predicts disease progression. Some think that SPVL is evolving, indicating that the virulence of the virus may be changing, but the data are not consistent. In addition, the widespread use of antiretroviral therapy (ART) has the potential to drive virulence evolution. We develop a simple deterministic model designed to answer the following questions: what are the expected patterns of virulence change in the initial decades of an epidemic? Could administration of ART drive changes in virulence evolution and, what is the potential size and direction of this effect? We find that even without ART we would not expect monotonic changes in average virulence. Transient decreases in virulence following the peak of an epidemic are not necessarily indicative of eventual evolution to avirulence. In the short term, we would expect widespread ART to cause limited downward pressure on virulence. In the long term, the direction of the effect is determined by a threshold condition, which we define. We conclude that, given the surpassing benefits of ART to the individual and in reducing onward transmission, virulence evolution considerations need have little bearing on how we treat.

A trend towards increasing viral load in newly diagnosed HIV-infected inpatients in southeast China

Peripheral blood viral load is an important indicator of viral production and clearance. Previous studies have suggested that viral load might predict the rate of decrease in CD4+ cell count and progression to AIDS and death. Here, we conducted a retrospective analysis of the trends in HIV-1 viral load in southeast China. Among inpatients newly diagnosed with HIV infection, we found that viral load has increased over the past decade from 4·20 log10 copies/ml in 2002 to 6·61 log10 copies/ml in 2014, with a mean increase of 0·19 log10 copies/ml each year. However, the CD4+ cell count was stable and insensitive to changes in viral load. Thus, increasing viral load appears to be an emerging trend in newly diagnosed HIV-infected inpatients.

Impact of pre-adapted HIV transmission

Human leukocyte antigen class I (HLA)-restricted CD8(+) T lymphocyte (CTL) responses are crucial to HIV-1 control. Although HIV can evade these responses, the longer-term impact of viral escape mutants remains unclear, as these variants can also reduce intrinsic viral fitness. To address this, we here developed a metric to determine the degree of HIV adaptation to an HLA profile. We demonstrate that transmission of viruses that are pre-adapted to the HLA molecules expressed in the recipient is associated with impaired immunogenicity, elevated viral load and accelerated CD4(+) T cell decline. Furthermore, the extent of pre-adaptation among circulating viruses explains much of the variation in outcomes attributed to the expression of certain HLA alleles. Thus, viral pre-adaptation exploits 'holes' in the immune response. Accounting for these holes may be key for vaccine strategies seeking to elicit functional responses from viral variants, and to HIV cure strategies that require broad CTL responses to achieve successful eradication of HIV reservoirs.

The impact of host genetic variation on infection with HIV-1

The outcome after infection with the human immunodeficiency virus type 1 (HIV-1) is a complex phenotype determined by interactions among the pathogen, the human host and the surrounding environment. An impact of host genetic variation on HIV-1 susceptibility was identified early in the pandemic, with a major role attributed to the genes encoding class I human leukocyte antigens (HLA) and the chemokine receptor CCR5. Studies using genome-wide data sets have underscored the strength of these associations relative to variants located throughout the rest of the genome. However, the extent to which additional polymorphisms influence HIV-1 disease progression, and how much of the variability in outcome can be attributed to host genetics, remain largely unclear. Here we discuss findings concerning the functional impact of associated variants, outline methods for quantifying the host genetic component and examine how available genome-wide data sets may be leveraged to discover gene variants that affect the outcome of HIV-1 infection.