1
|
A de novo approach to inferring within-host fitness effects during untreated HIV-1 infection. PLoS Pathog 2020; 16:e1008171. [PMID: 32492061 PMCID: PMC7295245 DOI: 10.1371/journal.ppat.1008171] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Revised: 06/15/2020] [Accepted: 05/11/2020] [Indexed: 12/15/2022] Open
Abstract
In the absence of effective antiviral therapy, HIV-1 evolves in response to the within-host environment, of which the immune system is an important aspect. During the earliest stages of infection, this process of evolution is very rapid, driven by a small number of CTL escape mutations. As the infection progresses, immune escape variants evolve under reduced magnitudes of selection, while competition between an increasing number of polymorphic alleles (i.e., clonal interference) makes it difficult to quantify the magnitude of selection acting upon specific variant alleles. To tackle this complex problem, we developed a novel multi-locus inference method to evaluate the role of selection during the chronic stage of within-host infection. We applied this method to targeted sequence data from the p24 and gp41 regions of HIV-1 collected from 34 patients with long-term untreated HIV-1 infection. We identify a broad distribution of beneficial fitness effects during infection, with a small number of variants evolving under strong selection and very many variants evolving under weaker selection. The uniquely large number of infections analysed granted a previously unparalleled statistical power to identify loci at which selection could be inferred to act with statistical confidence. Our model makes no prior assumptions about the nature of alleles under selection, such that any synonymous or non-synonymous variant may be inferred to evolve under selection. However, the majority of variants inferred with confidence to be under selection were non-synonymous in nature, and in most cases were have previously been associated with either CTL escape in p24 or neutralising antibody escape in gp41. We also identified a putative new CTL escape site (residue 286 in gag), and a region of gp41 (including residues 644, 648, 655 in env) likely to be associated with immune escape. Sites inferred to be under selection in multiple hosts have high within-host and between-host diversity although not all sites with high between-host diversity were inferred to be under selection at the within-host level. Our identification of selection at sites associated with resistance to broadly neutralising antibodies (bNAbs) highlights the need to fully understand the role of selection in untreated individuals when designing bNAb based therapies.
Collapse
|
2
|
Garcia V, Bonhoeffer S, Fu F. Cancer-induced immunosuppression can enable effectiveness of immunotherapy through bistability generation: A mathematical and computational examination. J Theor Biol 2020; 492:110185. [PMID: 32035826 PMCID: PMC7079339 DOI: 10.1016/j.jtbi.2020.110185] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2019] [Revised: 01/09/2020] [Accepted: 02/03/2020] [Indexed: 12/22/2022]
Abstract
The presence of an immunological barrier in cancer- immune system interaction (CISI) is consistent with the bistability patterns in that system. In CISI models, bistability patterns are consistent with immunosuppressive effects dominating immunoproliferative effects. Bistability could be harnessed to devise effective combination immunotherapy approaches.
Cancer immunotherapies rely on how interactions between cancer and immune system cells are constituted. The more essential to the emergence of the dynamical behavior of cancer growth these interactions are, the more effectively they may be used as mechanisms for interventions. Mathematical modeling can help unearth such connections, and help explain how they shape the dynamics of cancer growth. Here, we explored whether there exist simple, consistent properties of cancer-immune system interaction (CISI) models that might be harnessed to devise effective immunotherapy approaches. We did this for a family of three related models of increasing complexity. To this end, we developed a base model of CISI, which captures some essential features of the more complex models built on it. We find that the base model and its derivates can plausibly reproduce biological behavior that is consistent with the notion of an immunological barrier. This behavior is also in accord with situations in which the suppressive effects exerted by cancer cells on immune cells dominate their proliferative effects. Under these circumstances, the model family may display a pattern of bistability, where two distinct, stable states (a cancer-free, and a full-grown cancer state) are possible. Increasing the effectiveness of immune-caused cancer cell killing may remove the basis for bistability, and abruptly tip the dynamics of the system into a cancer-free state. Additionally, in combination with the administration of immune effector cells, modifications in cancer cell killing may be harnessed for immunotherapy without the need for resolving the bistability. We use these ideas to test immunotherapeutic interventions in silico in a stochastic version of the base model. This bistability-reliant approach to cancer interventions might offer advantages over those that comprise gradual declines in cancer cell numbers.
Collapse
Affiliation(s)
- Victor Garcia
- Institute of Applied Simulation, Zurich University of Applied Sciences, Einsiedlerstrasse 31a, 8820 Wädenswil, Switzerland; ETH Zurich, Universitätstrasse 16, 8092 Zürich, Switzerland; Institute for Social and Preventive Medicine, University of Bern, Finkenhubelweg 11, 3012 Bern, Switzerland; Department of Biology, Stanford University, 371 Serra Mall, Stanford CA 94305, USA.
| | | | - Feng Fu
- Department of Mathematics, Dartmouth College, 27 N. Main Street, 6188 Kemeny Hall, Hanover, NH 03755-3551, USA; ETH Zurich, Universitätstrasse 16, 8092 Zürich, Switzerland
| |
Collapse
|
3
|
Garcia V, Glassberg EC, Harpak A, Feldman MW. Clonal interference can cause wavelet-like oscillations of multilocus linkage disequilibrium. J R Soc Interface 2019; 15:rsif.2017.0921. [PMID: 29563246 DOI: 10.1098/rsif.2017.0921] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2017] [Accepted: 02/23/2018] [Indexed: 11/12/2022] Open
Abstract
Within-host adaptation of pathogens such as human immunodeficiency virus (HIV) often occurs at more than two loci. Multiple beneficial mutations may arise simultaneously on different genetic backgrounds and interfere, affecting each other's fixation trajectories. Here, we explore how these evolutionary dynamics are mirrored in multilocus linkage disequilibrium (MLD), a measure of multi-way associations between alleles. In the parameter regime corresponding to HIV, we show that deterministic early infection models induce MLD to oscillate over time in a wavelet-like fashion. We find that the frequency of these oscillations is proportional to the rate of adaptation. This signature is robust to drift, but can be eroded by high variation in fitness effects of beneficial mutations. Our findings suggest that MLD oscillations could be used as a signature of interference among multiple equally advantageous mutations and may aid the interpretation of MLD in data.
Collapse
Affiliation(s)
- Victor Garcia
- Department of Biology, Stanford University, 371 Serra Mall, Stanford, CA 94305, USA
| | - Emily C Glassberg
- Department of Biology, Stanford University, 371 Serra Mall, Stanford, CA 94305, USA
| | - Arbel Harpak
- Department of Biology, Stanford University, 371 Serra Mall, Stanford, CA 94305, USA
| | - Marcus W Feldman
- Department of Biology, Stanford University, 371 Serra Mall, Stanford, CA 94305, USA
| |
Collapse
|
4
|
Baral S, Raja R, Sen P, Dixit NM. Towards multiscale modeling of the CD8 + T cell response to viral infections. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2019; 11:e1446. [PMID: 30811096 PMCID: PMC6614031 DOI: 10.1002/wsbm.1446] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/24/2018] [Revised: 01/23/2019] [Accepted: 01/28/2019] [Indexed: 12/22/2022]
Abstract
The CD8+ T cell response is critical to the control of viral infections. Yet, defining the CD8+ T cell response to viral infections quantitatively has been a challenge. Following antigen recognition, which triggers an intracellular signaling cascade, CD8+ T cells can differentiate into effector cells, which proliferate rapidly and destroy infected cells. When the infection is cleared, they leave behind memory cells for quick recall following a second challenge. If the infection persists, the cells may become exhausted, retaining minimal control of the infection while preventing severe immunopathology. These activation, proliferation and differentiation processes as well as the mounting of the effector response are intrinsically multiscale and collective phenomena. Remarkable experimental advances in the recent years, especially at the single cell level, have enabled a quantitative characterization of several underlying processes. Simultaneously, sophisticated mathematical models have begun to be constructed that describe these multiscale phenomena, bringing us closer to a comprehensive description of the CD8+ T cell response to viral infections. Here, we review the advances made and summarize the challenges and opportunities ahead. This article is categorized under: Analytical and Computational Methods > Computational Methods Biological Mechanisms > Cell Fates Biological Mechanisms > Cell Signaling Models of Systems Properties and Processes > Mechanistic Models.
Collapse
Affiliation(s)
- Subhasish Baral
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, India
| | - Rubesh Raja
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, India
| | - Pramita Sen
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, India
| | - Narendra M Dixit
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, India.,Centre for Biosystems Science and Engineering, Indian Institute of Science, Bangalore, India
| |
Collapse
|
5
|
Abstract
The interplay between immune response and HIV is intensely studied via mathematical modeling, with significant insights but few direct answers. In this short review, we highlight advances and knowledge gaps across different aspects of immunity. In particular, we identify the innate immune response and its role in priming the adaptive response as ripe for modeling. The latter have been the focus of most modeling studies, but we also synthesize key outstanding questions regarding effector mechanisms of cellular immunity and development of broadly neutralizing antibodies. Thus far, most modeling studies aimed to infer general immune mechanisms; we foresee that significant progress will be made next by detailed quantitative fitting of models to data, and prediction of immune responses.
Collapse
Affiliation(s)
- Jessica M Conway
- Department of Mathematics and Center for Infectious Disease Dynamics, Pennsylvania State University, University Park PA 16802, USA
| | - Ruy M Ribeiro
- Laboratorio de Biomatematica, Faculdade de Medicina da Universidade de Lisboa, Portugal and Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| |
Collapse
|
6
|
Yang Y, Ganusov VV. Kinetics of HIV-Specific CTL Responses Plays a Minimal Role in Determining HIV Escape Dynamics. Front Immunol 2018; 9:140. [PMID: 29472921 PMCID: PMC5810297 DOI: 10.3389/fimmu.2018.00140] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2017] [Accepted: 01/16/2018] [Indexed: 11/13/2022] Open
Abstract
Cytotoxic T lymphocytes (CTLs) have been suggested to play an important role in controlling human immunodeficiency virus (HIV-1 or simply HIV) infection. HIV, due to its high mutation rate, can evade recognition of T cell responses by generating escape variants that cannot be recognized by HIV-specific CTLs. Although HIV escape from CTL responses has been well documented, factors contributing to the timing and the rate of viral escape from T cells have not been fully elucidated. Fitness costs associated with escape and magnitude of the epitope-specific T cell response are generally considered to be the key in determining timing of HIV escape. Several previous analyses generally ignored the kinetics of T cell responses in predicting viral escape by either considering constant or maximal T cell response; several studies also considered escape from different T cell responses to be independent. Here, we focus our analysis on data from two patients from a recent study with relatively frequent measurements of both virus sequences and HIV-specific T cell response to determine impact of CTL kinetics on viral escape. In contrast with our expectation, we found that including temporal dynamics of epitope-specific T cell response did not improve the quality of fit of different models to escape data. We also found that for well-sampled escape data, the estimates of the model parameters including T cell killing efficacy did not strongly depend on the underlying model for escapes: models assuming independent, sequential, or concurrent escapes from multiple CTL responses gave similar estimates for CTL killing efficacy. Interestingly, the model assuming sequential escapes (i.e., escapes occurring along a defined pathway) was unable to accurately describe data on escapes occurring rapidly within a short-time window, suggesting that some of model assumptions must be violated for such escapes. Our results thus suggest that the current sparse measurements of temporal CTL dynamics in blood bear little quantitative information to improve predictions of HIV escape kinetics. More frequent measurements using more sensitive techniques and sampling in secondary lymphoid tissues may allow to better understand whether and how CTL kinetics impacts viral escape.
Collapse
Affiliation(s)
- Yiding Yang
- Department of Microbiology, University of Tennessee, Knoxville, TN, United States
| | - Vitaly V. Ganusov
- Department of Microbiology, University of Tennessee, Knoxville, TN, United States
- National Institute for Mathematical and Biological Synthesis, University of Tennessee, Knoxville, TN, United States
- Department of Mathematics, University of Tennessee, Knoxville, TN, United States
| |
Collapse
|
7
|
Yaseen MM, Abuharfeil NM, Alqudah MA, Yaseen MM. Mechanisms and Factors That Drive Extensive Human Immunodeficiency Virus Type-1 Hypervariability: An Overview. Viral Immunol 2017; 30:708-726. [PMID: 29064351 DOI: 10.1089/vim.2017.0065] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The extensive hypervariability of human immunodeficiency virus type-1 (HIV-1) populations represents a major barrier against the success of currently available antiretroviral therapy. Moreover, it is still the most important obstacle that faces the development of an effective preventive vaccine against this infectious virus. Indeed, several factors can drive such hypervariability within and between HIV-1 patients. These factors include: first, the very low fidelity nature of HIV-1 reverse transcriptase; second, the extremely high HIV-1 replication rate; and third, the high genomic recombination rate that the virus has. All these factors together with the APOBEC3 proteins family and the immune and antiviral drugs pressures drive the extensive hypervariability of HIV-1 populations. Studying these factors and the mechanisms that drive such hypervariability will provide valuable insights that may guide the development of effective therapeutic and preventive strategies against HIV-1 infection in the near future. To this end, in this review, we summarized recent advances in this area of HIV-1 research.
Collapse
Affiliation(s)
- Mahmoud Mohammad Yaseen
- 1 Department of Medical Laboratory Sciences, College of Applied Medical Sciences, Jordan University of Science and Technology , Irbid, Jordan
| | - Nizar Mohammad Abuharfeil
- 2 Department of Applied Biological Sciences, College of Science and Arts, Jordan University of Science and Technology , Irbid, Jordan
| | - Mohammad Ali Alqudah
- 3 Department of Clinical Pharmacy, College of Pharmacy, Jordan University of Science and Technology , Irbid, Jordan
| | - Mohammad Mahmoud Yaseen
- 4 Department of Public Health, College of Medicine, Jordan University of Science and Technology , Irbid, Jordan
| |
Collapse
|
8
|
Garcia V, Feldman MW. Within-Epitope Interactions Can Bias CTL Escape Estimation in Early HIV Infection. Front Immunol 2017; 8:423. [PMID: 28507544 PMCID: PMC5410659 DOI: 10.3389/fimmu.2017.00423] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2016] [Accepted: 03/27/2017] [Indexed: 01/03/2023] Open
Abstract
As human immunodeficiency virus (HIV) begins to replicate within hosts, immune responses are elicited against it. Escape mutations in viral epitopes—immunogenic peptide parts presented on the surface of infected cells—allow HIV to partially evade these responses, and thus rapidly go to fixation. The faster they go to fixation, i.e., the higher their escape rate, the larger the selective pressure exerted by the immune system is assumed to be. This relation underpins the rationale for using escapes to assess the strength of immune responses. However, escape rate estimates are often obtained by employing an aggregation procedure, where several mutations that affect the same epitope are aggregated into a single, composite epitope mutation. The aggregation procedure thus rests upon the assumption that all within-epitope mutations have indistinguishable effects on immune recognition. In this study, we investigate how violation of this assumption affects escape rate estimates. To this end, we extend a previously developed simulation model of HIV that accounts for mutation, selection, and recombination to include different distributions of fitness effects (DFEs) and inter-mutational genomic distances. We use this discrete time Wright–Fisher based model to simulate early within-host evolution of HIV for DFEs and apply standard estimation methods to infer the escape rates. We then compare true with estimated escape rate values. We also compare escape rate values obtained by applying the aggregation procedure with values estimated without use of that procedure. We find that across the DFEs analyzed, the aggregation procedure alters the detectability of escape mutations: large-effect mutations are overrepresented while small-effect mutations are concealed. The effect of the aggregation procedure is similar to extracting the largest-effect mutation appearing within an epitope. Furthermore, the more pronounced the over-exponential decay of the DFEs, the more severely true escape rates are underestimated. We conclude that the aggregation procedure has two main consequences. On the one hand, it leads to a misrepresentation of the DFE of fixed mutations. On the other hand, it conceals within-epitope interactions that may generate irregularities in mutation frequency trajectories that are thus left unexplained.
Collapse
Affiliation(s)
- Victor Garcia
- Department of Biology, Stanford University, Stanford, CA, USA
| | | |
Collapse
|
9
|
Garcia V, Feldman MW, Regoes RR. Investigating the Consequences of Interference between Multiple CD8+ T Cell Escape Mutations in Early HIV Infection. PLoS Comput Biol 2016; 12:e1004721. [PMID: 26829720 PMCID: PMC4735108 DOI: 10.1371/journal.pcbi.1004721] [Citation(s) in RCA: 182] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2015] [Accepted: 12/18/2015] [Indexed: 12/14/2022] Open
Abstract
During early human immunodeficiency virus (HIV) infection multiple CD8+ T cell responses are elicited almost simultaneously. These responses exert strong selective pressures on different parts of HIV’s genome, and select for mutations that escape recognition and are thus beneficial to the virus. Some studies reveal that the later these escape mutations emerge, the more slowly they go to fixation. This pattern of escape rate decrease(ERD) can arise by distinct mechanisms. In particular, in large populations with high beneficial mutation rates interference among different escape strains –an effect that can emerge in evolution with asexual reproduction and results in delayed fixation times of beneficial mutations compared to sexual reproduction– could significantly impact the escape rates of mutations. In this paper, we investigated how interference between these concurrent escape mutations affects their escape rates in systems with multiple epitopes, and whether it could be a source of the ERD pattern. To address these issues, we developed a multilocus Wright-Fisher model of HIV dynamics with selection, mutation and recombination, serving as a null-model for interference. We also derived an interference-free null model assuming initial neutral evolution before immune response elicitation. We found that interference between several equally selectively advantageous mutations can generate the observed ERD pattern. We also found that the number of loci, as well as recombination rates substantially affect ERD. These effects can be explained by the underexponential decline of escape rates over time. Lastly, we found that the observed ERD pattern in HIV infected individuals is consistent with both independent, interference-free mutations as well as interference effects. Our results confirm that interference effects should be considered when analyzing HIV escape mutations. The challenge in estimating escape rates and mutation-associated selective coefficients posed by interference effects cannot simply be overcome by improved sampling frequencies or sizes. This problem is a consequence of the fundamental shortcomings of current estimation techniques under interference regimes. Hence, accounting for the stochastic nature of competition between mutations demands novel estimation methodologies based on the analysis of HIV strains, rather than mutation frequencies. Within-host evolution of human immunodeficiency virus (HIV) is shaped by strong immune responses mounted against the virus. Multiple CD8+ T cell populations, each recognizing a specific part of an HIV protein, simultaneously suppress HIV growth. Escape mutations that arise in HIV genome regions coding for these virion protein parts, impair CD8+ T cell recognition and are consequently strongly selected. The emergence and rise of these escape mutations exhibits an intriguing temporal pattern: the earlier an escape mutation arises, the faster it goes to fixation. This pattern is termed escape rate decrease (ERD). In this paper, we tested computationally whether interference, i.e. the coexistence of multiple genetically distinct HIV strains engaged in competitive interaction within the host, could be a possible source of ERD. As an alternative, we also mathematically derived the temporal pattern of escapes under interference-free conditions, and compared this with data. We found that interference between multiple beneficial mutations could generate ERD. However, ERD does not imply the presence of interference. Thus, more detailed data is required to unambiguously determine whether interference effects influence ERD generation. Nevertheless, interference should be considered when studying the within-host evolution of HIV. Ignoring its effects on population dynamics can severely underestimate the protective capacity of CD8+ T cells.
Collapse
Affiliation(s)
- Victor Garcia
- Institute of Integrative Biology, ETH Zurich, Zurich, Switzerland
- * E-mail:
| | - Marcus W. Feldman
- Department of Biological Sciences, Stanford University, Stanford, California, United States of America
| | - Roland R. Regoes
- Institute of Integrative Biology, ETH Zurich, Zurich, Switzerland
| |
Collapse
|
10
|
Leviyang S, Ganusov VV. Broad CTL Response in Early HIV Infection Drives Multiple Concurrent CTL Escapes. PLoS Comput Biol 2015; 11:e1004492. [PMID: 26506433 PMCID: PMC4624722 DOI: 10.1371/journal.pcbi.1004492] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2014] [Accepted: 08/06/2015] [Indexed: 12/15/2022] Open
Abstract
Recent studies have highlighted the ability of HIV to escape from cytotoxic T lymphocyte (CTL) responses that concurrently target multiple viral epitopes. Yet, the viral dynamics involved in such escape are incompletely understood. Previous analyses have made several strong assumptions regarding HIV escape from CTL responses such as independent or non-concurrent escape from individual CTL responses. Using experimental data from evolution of HIV half genomes in four patients we observe concurrent viral escape from multiple CTL responses during early infection (first 100 days of infection), providing confirmation of a recent result found in a study of one HIV-infected patient. We show that current methods of estimating CTL escape rates, based on the assumption of independent escapes, are biased and perform poorly when CTL escape proceeds concurrently at multiple epitopes. We propose a new method for analyzing longitudinal sequence data to estimate the rate of CTL escape across multiple epitopes; this method involves few parameters and performs well in simulation studies. By applying our novel method to experimental data, we find that concurrent multiple escapes occur at rates between 0.03 and 0.4 day−1, a relatively broad range that reflects uncertainty due to sparse sampling and wide ranges of parameter values. However, we show that concurrent escape at rates 0.1–0.2 day−1 across multiple epitopes is consistent with our patient datasets. Since the early 1990s, cytotoxic T lymphocytes (CTLs) have been known to play an important role in HIV infection with CTLs targeting HIV epitopes and, in turn, HIV escapes arising through mutations in the targeted epitopes. Over the past decade, studies have shown that CTL responses concurrently target multiple HIV epitopes, yet the effect of concurrent responses on HIV dynamics and evolution is not well understood. Through an analysis of patient datasets and a novel statistical method, we show that during early HIV infection concurrent CTL responses drive concurrent HIV escapes at multiple epitopes with significant pressure, suggesting a complex picture in which HIV simultaneously explores multiple mutational pathways to escape from broad and potent CTL response.
Collapse
Affiliation(s)
- Sivan Leviyang
- Department of Mathematics and Statistics, Georgetown University, Washington, DC, United States of America
- * E-mail:
| | - Vitaly V. Ganusov
- Department of Microbiology, University of Tennessee, Knoxville, Tennessee, United States of America
| |
Collapse
|