Ha Youn Lee

Modeling Sequence Evolution in Acute HIV-1 Infection

We describe a mathematical model and Monte Carlo (MC) simulation of viral evolution during acute infection. We consider both synchronous and asynchronous processes of viral infection of new target cells. The model enables an assessment of the expected sequence diversity in new HIV-1 infections originating from a single transmitted viral strain, estimation of the most recent common ancestor (MRCA) of the transmitted viral lineage, and estimation of the time to coalesce back to the MRCA. We also calculate the probability of the MRCA being the transmitted virus or an evolved variant. Excluding insertions and deletions, we assume HIV-1 evolves by base substitution without selection pressure during the earliest phase of HIV-1 infection prior to the immune response. Unlike phylogenetic methods that follow a lineage backwards to coalescence, we compare the observed data to a model of the diversification of a viral population forward in time. To illustrate the application of these methods, we provide detailed comparisons of the model and simulations results to 306 envelope sequences obtained from eight newly infected subjects at a single time point. The data from 68 patients were in good agreement with model predictions, and hence compatible with a single-strain infection evolving under no selection pressure. The diversity of the samples from the other two patients was too great to be explained by the model, suggesting multiple HIV-1-strains were transmitted. The model can also be applied to longitudinal patient data to estimate within-host viral evolutionary parameters.

Transfusion transmission of retroviruses: human T-lymphotropic virus types I and II compared with human immunodeficiency virus type 1.

Background: The incidence of transfusion transmission of human T‐ lymphotropic virus type I (HTLV‐I) and HTLV type II (HTLV‐II) has not been compared directly or to that of human immunodeficiency virus type 1 (HIV‐1). The effects of refrigerator storage of the blood component on infectivity of the viruses needs definition.

Short-Lived Infected Cells Support Virus Replication in Sooty Mangabeys Naturally Infected with Simian Immunodeficiency Virus: Implications for AIDS Pathogenesis

ABSTRACT Sooty mangabeys (SMs) naturally infected with simian immunodeficiency virus (SIV) do not develop AIDS despite high levels of virus replication. At present, the mechanisms underlying this disease resistance are poorly understood. Here we tested the hypothesis that SIV-infected SMs avoid immunodeficiency as a result of virus replication occurring in infected cells that live significantly longer than human immunodeficiency virus (HIV)-infected human cells. To this end, we treated six SIV-infected SMs with potent antiretroviral therapy (ART) and longitudinally measured the decline in plasma viremia. We applied the same mathematical models used in HIV-infected individuals and observed that SMs naturally infected with SIV also present a two-phase decay of viremia following ART, with the bulk (92 to 99%) of virus replication sustained by short-lived cells (average life span, 1.06 days), and only 1 to 8% occurring in longer-lived cells. In addition, we observed that ART had a limited impact on CD4+ T cells and the prevailing level of T-cell activation and proliferation in SIV-infected SMs. Collectively, these results suggest that in SIV-infected SMs, similar to HIV type 1-infected humans, short-lived activated CD4+ T cells, rather than macrophages, are the main source of virus production. These findings indicate that a short in vivo life span of infected cells is a common feature of both pathogenic and nonpathogenic primate lentivirus infections and support a model for AIDS pathogenesis whereby the direct killing of infected cells by HIV is not the main determinant of disease progression.

Dynamic Correlation between Intrahost HIV-1 Quasispecies Evolution and Disease Progression

Quantifying the dynamics of intrahost HIV-1 sequence evolution is one means of uncovering information about the interaction between HIV-1 and the host immune system. In the chronic phase of infection, common dynamics of sequence divergence and diversity have been reported. We developed an HIV-1 sequence evolution model that simulated the effects of mutation and fitness of sequence variants. The amount of evolution was described by the distance from the founder strain, and fitness was described by the number of offspring a parent sequence produces. Analysis of the model suggested that the previously observed saturation of divergence and decrease of diversity in later stages of infection can be explained by a decrease in the proportion of offspring that are mutants as the distance from the founder strain increases rather than due to an increase of viral fitness. The prediction of the model was examined by performing phylogenetic analysis to estimate the change in the rate of evolution during infection. In agreement with our modeling, in 13 out of 15 patients (followed for 3–12 years) we found that the rate of intrahost HIV-1 evolution was not constant but rather slowed down at a rate correlated with the rate of CD4+ T-cell decline. The correlation between the dynamics of the evolutionary rate and the rate of CD4+ T-cell decline, coupled with our HIV-1 sequence evolution model, explains previously conflicting observations of the relationships between the rate of HIV-1 quasispecies evolution and disease progression.

Unequal evolutionary rates in the human immunodeficiency virus type 1 (HIV-1) pandemic: the evolutionary rate of HIV-1 slows down when the epidemic rate increases.

ABSTRACT HIV-1 sequences in intravenous drug user (IDU) networks are highly homogenous even after several years, while this is not observed in most sexual epidemics. To address this disparity, we examined the human immunodeficiency virus type 1 (HIV-1) evolutionary rate on the population level for IDU and heterosexual transmissions. All available HIV-1 env V3 sequences from IDU outbreaks and heterosexual epidemics with known sampling dates were collected from the Los Alamos HIV sequence database. Evolutionary rates were calculated using phylogenetic trees with a t test root optimization of dated samples. The evolutionary rate of HIV-1 subtype A1 was found to be 8.4 times lower in fast spread among IDUs in the former Soviet Union (FSU) than in slow spread among heterosexual individuals in Africa. Mixed epidemics (IDU and heterosexual) showed intermediate evolutionary rates, indicating a combination of fast- and slow-spread patterns. Hence, if transmissions occur repeatedly during the initial stage of host infection, before selective pressures of the immune system have much impact, the rate of HIV-1 evolution on the population level will decrease. Conversely, in slow spread, where HIV-1 evolves under the pressure of the immune system before a donor infects a recipient, the virus evolution at the population level will increase. Epidemiological modeling confirmed that the evolutionary rate of HIV-1 depends on the rate of spread and predicted that the HIV-1 evolutionary rate in a fast-spreading epidemic, e.g., for IDUs in the FSU, will increase as the population becomes saturated with infections and the virus starts to spread to other risk groups.

Short-lived infected cells support virus replication in naturally SIV-infected sooty 1 mangabeys: implications for AIDS pathogenesis 2

ABSTRACT Sooty mangabeys (SMs) naturally infected with simian immunodeficiency virus (SIV) do not develop AIDS despite high levels of virus replication. At present, the mechanisms underlying this disease resistance are poorly understood. Here we tested the hypothesis that SIV-infected SMs avoid immunodeficiency as a result of virus replication occurring in infected cells that live significantly longer than human immunodeficiency virus (HIV)-infected human cells. To this end, we treated six SIV-infected SMs with potent antiretroviral therapy (ART) and longitudinally measured the decline in plasma viremia. We applied the same mathematical models used in HIV-infected individuals and observed that SMs naturally infected with SIV also present a two-phase decay of viremia following ART, with the bulk (92 to 99%) of virus replication sustained by short-lived cells (average life span, 1.06 days), and only 1 to 8% occurring in longer-lived cells. In addition, we observed that ART had a limited impact on CD4+ T cells and the prevailing level of T-cell activation and proliferation in SIV-infected SMs. Collectively, these results suggest that in SIV-infected SMs, similar to HIV type 1-infected humans, short-lived activated CD4+ T cells, rather than macrophages, are the main source of virus production. These findings indicate that a short in vivo life span of infected cells is a common feature of both pathogenic and nonpathogenic primate lentivirus infections and support a model for AIDS pathogenesis whereby the direct killing of infected cells by HIV is not the main determinant of disease progression.

Designing a genome-based HIV incidence assay with high sensitivity and specificity.

Objective:Considerable inaccuracy in estimates of HIV incidence has been a serious obstacle to the development of efficient HIV/AIDS prevention and interventions. Accurately distinguishing recent or incident infections from chronic infections enables one to monitor epidemics and evaluate the impact of HIV prevention/intervention trials. However, serological testing has not been able to realize these promises due to a number of critical limitations. Our study is to design a novel scheme of identifying incident infections in a highly accurate manner, based on the characteristics of HIV gene diversification within an infected individual. Methods:We perform a comprehensive meta-analysis on 5596 full envelope HIV genes generated by single genome amplification-direct sequencing from 182 incident and 43 chronic cases. We devise a binary classification test based on the tail characteristics of the Hamming distance distribution of sequences. Results:We identify a clear signature of incident infections, the presence of closely related strains in the sampled HIV envelope gene sequences in each HIV-infected patient, in both single-variant and multivariant transmissions. The sequence similarity used as a biomarker is found to have high specificity and sensitivity, greater than 95%, and is robust to viral and host-specific factors such as the clade of the viral strain, viral load, and the length and location of sequences in the HIV envelope gene. Conclusion:Because of rapid and continuing improvements in sequencing technology and cost, sequence-based incidence assays hold great promise as a means of quantifying HIV incidence from a single blood test.

Treatment response in acute/early infection versus advanced AIDS: equivalent first and second phases of HIV RNA decline

Objective:Compare the initial phases of virologic decay when acute/early and advanced HIV-infected adults are administered the same treatment regimen. Design:Mathematical modeling of a previously completed prospective treatment pilot study involving treatment-naive patients with early and advanced immunosuppression. Methods:We analyzed data from a treatment protocol in which 18 individuals with acute or recent HIV-1 seroconversion and six patients with advanced AIDS were administered the same four-drug antiretroviral regimen. Initial treatment responses were compared by fitting a mathematical model to frequent viral load measurements in order to calculate the first and second phase kinetics of viral clearance, and also by comparing viral load suppression over 24 weeks. Patients were also comprehensively compared in terms of protease inhibitor drug levels, HIV-specific immune responses at baseline, and the presence of drug resistance-conferring mutations. Results:There was no statistically meaningful difference in first phase clearance of comparable high-level viremia in the two groups, whether protease inhibitor levels were inserted into the model or 100% antiviral drug effectiveness was assumed. In contrast, acute/early patients had inferior sustained responses than advanced patients, reflecting erratic adherence. Conclusions:Despite many years of intervening immune destruction, the initial virologic decay on therapy appears to be the same at the extremes of the HIV disease spectrum.

Mathematical Modeling of Ultradeep Sequencing Data Reveals that Acute CD8+ T-Lymphocyte Responses Exert Strong Selective Pressure in Simian Immunodeficiency Virus-Infected Macaques but Still Fail To Clear Founder Epitope Sequences

ABSTRACT The prominent role of antiviral cytotoxic CD8+ T-lymphocytes (CD8-TL) in containing the acute viremia of human and simian immunodeficiency viruses (HIV-1 and SIV) has rationalized the development of T-cell-based vaccines. However, the presence of escape mutations in the acute stage of infection has raised a concern that accelerated escape from vaccine-induced CD8-TL responses might undermine vaccine efficacy. We reanalyzed previously published data of 101,822 viral genomes of three CD8-TL epitopes, Nef103-111RM9 (RM9), Tat28-35SL8 (SL8), and Gag181-189CM9 (CM9), sampled by ultradeep pyrosequencing from eight macaques. Multiple epitope variants appeared during the resolution of acute viremia, followed by the predominance of a single mutant epitope. By fitting a mathematical model, we estimated the first acute escape rate as 0.36 day−1 within escape-prone epitopes, RM9 and SL8, and the chronic escape rate as 0.014 day−1 within the CM9 epitope. Our estimate of SIV acute escape rates was found to be comparable to very early HIV-1 escape rates. The timing of the first escape was more highly correlated with the timing of the peak CD8-TL response than with the magnitude of the CD8-TL response. The transmitted epitope decayed more than 400 times faster during the acute viral decline stage than predicted by a neutral evolution model. However, the founder epitope persisted as a minor population even at the viral set point; in contrast, the majority of acute escape epitopes were completely cleared. Our results suggest that a reservoir of SIV infection is preferentially formed by virus with the transmitted epitope.

SPMM: estimating infection duration of multivariant HIV-1 infections

MOTIVATION Illustrating how HIV-1 is transmitted and how it evolves in the following weeks is an important step for developing effective vaccination and prevention strategies. It is currently possible through DNA sequencing to account for the diverse array of viral strains within an infected individual. This provides an unprecedented opportunity to pinpoint when each patient was infected and which viruses were transmitted. RESULTS Here we develop a mathematical tool for early HIV-1 evolution within a subject whose infection originates either from a single or multiple viral variants. The shifted Poisson mixture model (SPMM) provides a quantitative guideline for segregating viral lineages, which in turn enables us to assess when a subject was infected. The infection duration estimated by SPMM showed a statistically significant linear relationship with that by Fiebig laboratory staging (P = 0.00059) among 37 acutely infected subjects. Our tool provides a functional approach to understanding early genetic diversity, one of the most important parameters for deciphering HIV-1 transmission and predicting the rate of disease progression. AVAILABILITY AND IMPLEMENTATION SPMM, webserver, is available at http://www.hayounlee.org/web-tools.html. CONTACT hayoun@usc.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.