Tracking members of the simian immunodeficiency virus deltaB670 quasispecies population in vivo at single-cell resolution

Large Variations in HIV-1 Viral Load Explained by Shifting-Mosaic Metapopulation Dynamics

The viral population of HIV-1, like many pathogens that cause systemic infection, is structured and differentiated within the body. The dynamics of cellular immune trafficking through the blood and within compartments of the body has also received wide attention. Despite these advances, mathematical models, which are widely used to interpret and predict viral and immune dynamics in infection, typically treat the infected host as a well-mixed homogeneous environment. Here, we present mathematical, analytical, and computational results that demonstrate that consideration of the spatial structure of the viral population within the host radically alters predictions of previous models. We study the dynamics of virus replication and cytotoxic T lymphocytes (CTLs) within a metapopulation of spatially segregated patches, representing T cell areas connected by circulating blood and lymph. The dynamics of the system depend critically on the interaction between CTLs and infected cells at the within-patch level. We show that for a wide range of parameters, the system admits an unexpected outcome called the shifting-mosaic steady state. In this state, the whole body’s viral population is stable over time, but the equilibrium results from an underlying, highly dynamic process of local infection and clearance within T-cell centers. Notably, and in contrast to previous models, this new model can explain the large differences in set-point viral load (SPVL) observed between patients and their distribution, as well as the relatively low proportion of cells infected at any one time, and alters the predicted determinants of viral load variation.

A novel metapopulation model of HIV suggests that within-host infections are characterized by a highly dynamic process of localized infection followed by clearance within T cell centers.

When a person is infected with HIV, the initial peak level of virus in the blood is usually very high before a lower, relatively stable level is reached and maintained for the duration of the chronic infection. This stable level is known as the set-point viral load (SPVL) and is associated with severity of infection. SPVL is also highly variable among patients, ranging from 100 to a million copies of the virus per mL of blood. The replicative capacity of the infecting virus and the strength of the immune response both influence SPVL. However, standard mathematical models show that variation in these two factors cannot easily reproduce the observed distribution of SPVL among patients. Standard models typically treat infected individuals as well-mixed systems, but in reality viral replication is localised in T-cell centres, or patches, found in secondary lymphoid tissue. To account for this population structure, we developed a carefully parameterised metapopulation model. We find the system can reach a steady state at which the viral load in the blood is relatively stable, representing SPVL, but surprisingly, the patches are highly dynamic, characterised by bursts of infection followed by elimination of virus due to localised host immune responses. Significantly, this model can reproduce the wide distribution of SPVLs found among infected individuals for realistic distributions of viral replicative capacity and strength of immune response. Our model can also be used in the future to understand other aspects of chronic HIV infection.

Increased APOBEC3G expression is associated with extensive G-to-A hypermutation in viral DNA in rhesus macaque brain during lentiviral infection

APOBEC3G restricts retrovirus replication through inducing guanosine-to-adenosine (G-to-A) hypermutations in viral DNA. Its role in brain "intrinsic immunity" has not been elucidated nor has it been convincingly demonstrated which brain cell compartments produce this defense factor in human immunodeficiency virus (HIV) infection, acquired immunodeficiency syndrome (AIDS), and antiretroviral therapy. Here, we investigated by immunohistochemistry and in situ hybridization the cell-specific regulation of APOBEC3G in rhesus macaque brains during infection with simian immunodeficiency virus (SIV) clone deltaB670, a primate model of HIV disease. We found that APOBEC3G protein and mRNA were exclusively expressed by some perivascular macrophages throughout the brain of noninfected and asymptomatic SIV-infected monkeys. Depending on virus burden, APOBEC3G was induced in microglia/macrophage-derived cells and T lymphocytes in late-stage SIV infection. Intracellularly, APOBEC3G was found in the cytoplasm and/or in the nucleus. APOBEC3G-positive cells were in close proximity to SIV gag-positive cells or were SIV-copositive. Induction of APOBEC3G was accompanied by G-to-A hypermutations in the gag and pol regions of retroviral DNA isolated from brain sections of AIDS-symptomatic monkeys. Although brain-directed treatment with antiretroviral 6-chloro-2',3'-dideoxyguanosine suppressed brain SIV burden, encephalitis and reduced cerebral APOBEC3G synthesis hypermutations were still detectable. Upregulation of APOBEC3G may restrict spread of SIV in the brain and thus limit brain damage during lentiviral infection.

Genetic analysis of simian immunodeficiency virus expressed in milk and selectively transmitted through breastfeeding

To develop effective intervention strategies that prevent breast milk transmission of human immunodeficiency virus (HIV), we must understand the specific viral properties and mechanisms responsible for infant infection. We have used lactating rhesus macaques infected with a pathogenic simian immunodeficiency virus (SIV) stock to analyze the viral genotypes expressed in plasma and milk throughout the disease course and to identify those variants ultimately transmitted to infants through breastfeeding. In these studies we observed mother-to-infant transmission of SIV/Delta(B670) by eight females during the chronic phase of disease, and we analyzed by heteroduplex tracking assays and sequence analysis the distribution and fluctuations in viral genotypes expressed. Each female expressed multiple V1 envelope genotypes in milk near the time of transmission, while a single genotype was found in each of the infants. Variants transmitted to infants were not expressed throughout the maternal disease course but were only detected near the time of transmission. The emergence of the transmitted genotype in the dam typically occurred in plasma before milk and was coincident with increased milk viral loads. Transmitted genotypes tended to be longer and more glycosylated and had a less negative charge over the V1 region compared to viral genotypes expressed in milk but not transmitted. These observations demonstrate that specific viral genotypes are selectively transmitted to infants through breastfeeding and support the hypothesis that transmission occurs as genotypes adapt for efficient expression in milk.

Progressive selection for neurovirulent genotypes in the brain of SIV-infected macaques

OBJECTIVE

To compare the viral genotypes present in RNA from brain and peripheral blood mononuclear cells (PBMC) and DNA from brain during acute, asymptomatic and late stages of SIV infection of macaques.

METHODS

Eighteen pigtailed macaques were intravenously inoculated with SIV. At 10, 21 and 56 days postinoculation, six were euthanized and the severity of encephalitis was assessed by microscopic examination. DNA and RNA were isolated from brain and PBMC, and the V1 region of env was amplified by the polymerase chain reaction and sequenced from over 800 different clones.

RESULTS

Similar genotypes were detected in RNA from brain and PBMC at 10 days postinoculation, suggesting an unrestricted exchange of virus between the periphery and the brain during acute infection. There was a progressive increase in the percentage of neurovirulent genotypes in brain RNA from acute (14% of all genotypes detected in brain RNA) to early asymptomatic (45%), to late asymptomatic (52%) and to terminal (95%) infection. Fewer different genotypes were found in brain RNA than in PBMC RNA from macaques euthanized during early asymptomatic (2.5 and 5 different genotypes, respectively; P = 0.007), late asymptomatic (2 and 5 different genotypes, respectively; P = 0.003) and terminal (2 and 4 different genotypes, respectively; P < 0.001) infection.

CONCLUSION

These data demonstrate that the almost exclusive replication of neurovirulent genotypes in the brain seen at late-stage infection is a progressive process that begins early in infection and continues to late stage disease.

Spatial models of virus-immune dynamics

To date, the majority of theoretical models describing the dynamics of infectious diseases in vivo are based on the assumption of well-mixed virus and cell populations. Because many infections take place in solid tissues, spatially structured models represent an important step forward in understanding what happens when the assumption of well-mixed populations is relaxed. Here, we explore models of virus and virus-immune dynamics where dispersal of virus and immune effector cells was constrained to occur locally. The stability properties of our spatial virus-immune dynamics models remained robust under almost all biologically plausible dispersal schemes, regardless of their complexity. The various spatial dynamics were compared to the basic non-spatial dynamics and important differences were identified: When space was assumed to be homogeneous, the dynamics generated by non-spatial and spatially structured models differed substantially at the peak of the infection. Thus, non-spatial models may lead to systematic errors in the estimates of parameters underlying acute infection dynamics. When space was assumed to be heterogeneous, spatial coupling not only changed the equilibrium properties of the uncoupled populations but also equalized the dynamics and thereby reduced the likelihood of dynamic elimination of the infection. In line with experimental and clinical observations, long-lasting oscillation periods were virtually absent. When source-sink dynamics were considered, the long-term outcome of the infection depended critically on the degree of spatial coupling. The infection collapsed when emigration from source sites became too large. Finally, we discuss the implications of spatially structured models on medical treatment of infectious diseases, and note that a huge gap exists in data accurately describing infection dynamics in solid tissues.

Roles of substrate availability and infection of resting and activated CD4+ T cells in transmission and acute simian immunodeficiency virus infection

In studies of sexual mucosal transmission and early stages of simian immunodeficiency virus (SIV) and HIV infections, productive infection predominates in CD4(+) T cell populations, in both ostensibly resting and activated cells. The surprising ability of SIV and HIV to replicate in resting cells in vivo, in contrast to propagation of infection in vitro, suggested a model in which during the early stages of infection these viruses exploit the greater availability of resting cells to maintain unbroken chains of transmission from an infected resting cell to another resting cell nearby. Because immune activation in response to infection provides more activated CD4(+) T cells, these viruses take advantage of the greater efficiency of virus production and spread in activated cells for propagation and dissemination of infection. In this article, we report the results of experimental tests of this model, including visualization at the light microscopic level and direct analysis of virus production by cells in tissues. Analysis of tissues of rhesus macaques inoculated intravaginally or i.v. with SIV supports the proposed roles of target cell availability, susceptibility, and virus production by infected resting and activated CD4(+) T cells in mucosal transmission and early infection, and points to a potential role for topical anti-inflammatory agents in moderating the initial propagation of infection.

Transmission, acute HIV-1 infection and the quest for strategies to prevent infection

By the acute stage of HIV-1 infection, the immune system already faces daunting challenges. Research on mucosal barriers and the events immediately after heterosexual transmission that precede this acute stage could facilitate the development of effective microbicides and vaccines.

Role of microglial cells in selective replication of simian immunodeficiency virus genotypes in the brain

An accelerated, consistent macaque simian immunodeficiency virus (SIV) model in which over 90% of pigtailed macaques (Macaca nemestrina) coinoculated with SIV/17E-Fr and SIV/DeltaB670 developed encephalitis was used to determine whether central nervous system (CNS) lesions are associated with the replication of specific genotypes in the brain and, more specifically, in the microglia. Ten of 11 inoculated macaques had severe (n = 3), moderate (n = 5), or mild (n = 2) encephalitis at 3 months postinoculation. To compare actively replicating viral genotypes in the CNS and in microglia with those in the periphery, the V1 region of the SIV envelope gene was amplified and sequenced from RNA extracted from basal ganglia, from microglial cells isolated from the brain, and from peripheral blood mononuclear cells (PBMC) isolated from blood at the time of death. To distinguish between actively replicating with latent viral genotypes in the CNS, viral genotypes in RNA and DNA from basal ganglia were compared. Two macrophage-tropic, neurovirulent viruses, SIV/17E-Fr and SIV/DeltaB670 Cl-2, predominated in the brain RNA of macaques with encephalitis, comprising 95% of the genotypes detected. The same two viral genotypes were present at the same frequencies in microglial cell RNA, suggesting that microglia are pivotal in the selective replication of neurovirulent viruses. There was a significantly greater number of viral genotypes in DNA than there were in RNA in the brain (P = 0.004), including those of both the macrophage- and lymphocyte-tropic viral strains. Furthermore, significantly fewer viral genotypes were detected in brain RNA than in PBMC RNA at the time of death (P = 0.004) and the viral strain that predominated in the brain frequently was different from that which predominated in the PBMC of the same animal. These data suggest that many viral genotypes enter the brain, but only a limited subset of macrophage-tropic, neurovirulent viruses replicate terminally in the brains of macaques with encephalitis. They further suggest that the selection of macrophage-tropic, neurovirulent viruses occurs not at the level of the blood-brain barrier but at a stage after virus entry and that microglial cells may play an important role in that selection process.

Spatiotemporal Dynamics of HIV Propagation

Although viral propagation is a localized process, mathematical models of viral replication kinetics have been limited to systems of ordinary differential equations describing spatially averaged behavior. In this paper, we introduce a cellular automaton model of viral propagation based on the known biophysical properties of HIV. In particular, we include the competition between viral lability and Brownian motion. The model predicts three testable effects not present in previous descriptions. First, we find a profound dependence of viral infectivity on cell concentration; virion instability decreases infectivity more than 100-fold under typical experimental conditions, resulting in misleading estimates of the number of infectious particles. Second, we find that, in a large parameter regime, infection extinguishes itself due to insufficient target cell replenishment. Finally, we find that propagation is limited by viral stability at low cell density and by geometry at high cell density. The geometry-limited regime can be modulated by downregulation of CD4. These different properties are analysed quantitatively and compared with previous experimental results.

Recovery of replication-competent HIV type 1-infected circulating monocytes from individuals receiving antiretroviral therapy

The affect of antiretroviral therapy (ART) on HIV-1 recovery from blood monocytes was determined in purified peripheral blood monocyte-derived macrophage (MDM) cultures from HIV-1-infected subjects with undetectable plasma viremia or active viral replication. Additionally, the association between replication-competent HIV-1-infected MDM and neurocognitive status was examined. Fifty-two individual with previous AIDS-defining illnesses receiving nucleoside analogues with and without protease inhibitors or no ART were followed for up to 1.5 years. Detection of plasma viremia significantly correlated with the occurrence of infected monocytes. Viral replication was detected in less than 10% of the MDM cultures from 23 individuals receiving effective antiretroviral therapy. In contrast, approximately 50% of the MDM cultures from 29 individuals with active viral replication and evidence of decreased immune function, including all individuals with neurocognitive impairment, produced detectable virus indicating that a lack of adequate ART results in increased abundance of replication-competent blood monocytes. Proviral DNA levels were a minimum of 13-fold higher in MDM from subjects with active viral replication. The infrequent detection of viral DNA in cultures from individuals receiving effective ART suggested low levels of circulating monocytes harboring replication-incompetent virus. These studies demonstrate that HIV-infected individuals on ART with breakthrough viremia have significantly higher levels of circulating infected monocytes, the precursors of tissue macrophages.

Pathogenesis of SIV pneumonia: selective replication of viral genotypes in the lung

Lymphocytic interstitial pneumonia of HIV-infected individuals and SIV pneumonia of macaques are both characterized by diffuse infiltration of the lungs with lymphocytes, plasma cells, and macrophages. This study was undertaken to determine whether there are specific, macrophage-tropic genotypes that selectively replicate in the lung of macaques with SIV pneumonia, as in SIV encephalitis. Using a rapid, reproducible SIV/macaque model of AIDS, 11 pig-tailed macaques were intravenously inoculated with an immunosuppressive viral strain, SIV/DeltaB670, and a macrophage-tropic molecule clone, SIV/17E-Fr, and euthanized at 3 months postinoculation. All 11 macaques had severe (6 macaques) or moderate (5 macaques) pneumonia. To identify the viral genotypes that were replicating in the lung parenchyma, bronchoalveolar lavage (BAL) cells, and peripheral blood mononuclear cells (PBMC) of each macaque, RNA was isolated and the SIV env V1 region was amplified, cloned, and sequenced. Lung homogenates and BAL cells contained a more limited repertoire of viral genotypes than PBMC. SIV/17E-Fr was the major genotype in the lungs of 5 macaques and in BAL cells of 6 macaques. The remainder of the macaques had SIV/17E-Fr and the macrophage-tropic strains of SIV/DeltaB670 clones 2 and 12. In contrast, SIV/17E-Fr was the predominant strain in the PBMC of only 3 of 11 macaques. The viral strain that predominated in PBMC was rarely the strain that predominated in the lungs (only 3 of 11 macaques). The severity of pulmonary lesions did not correlate with the levels of viral RNA in lung homogenates or in plasma. However, when only SIV/17E-Fr was expressed in the lung, the viral load in the lung was significantly higher (P = 0.016) than when SIV/DeltaB670 was present alone or in combination with SIV/17E-Fr. These data suggest that SIV pneumonia is associated with selective replication of specific macrophage-tropic genotypes in the lung and that SIV/17E-Fr has a selective advantage for replication in the lung.

Variation in simian immunodeficiency virus env V1 region in simian AIDS-associated lymphoma

Genetic variation of SIV env during the course of infection provides a large population pool that is continually shaped by selective forces in vivo and may influence the development of clinical disease. SAIDS-associated lymphoma (SAL) in the SIV-infected macaque is typically a clonal or oligoclonal mass of B cell origin, extranodal in anatomic distribution, in which SIV is restricted largely to infiltrating macrophages. To explore the degree of genetic variation in SIV env represented in SAL, a 480-bp DNA fragment containing the V1 region was PCR amplified from seven cases of SAL and from a nonneoplastic lymph node of an SIV-infected macaque. The nucleotide sequence of the V1 region was determined from at least 10 clones from multiple independent amplification reactions of each tissue. Overall, the degree of V1 variability within lymphomas was found not to be restricted but to resemble the heterogeneity reported in SIV-infected lymphoid and other tissues. V1 variation in the nonneoplastic lymph node was unexpectedly limited, perhaps related to the unusual disease condition associated with SAIDS in that animal. Unlike observations from SIV-infected tissues of animals without neoplastic disease, no increase was detected in the number of O- or N-linked glycosylation sites in the V1 regions isolated from lymphomas as compared with the original inoculum. These findings suggest that, within the microenvironment of the lymphoma, the immune evasion conferred by increased glycosylation may offer little selective advantage.

Transition between stochastic evolution and deterministic evolution in the presence of selection: general theory and application to virology

We present here a self-contained analytic review of the role of stochastic factors acting on a virus population. We develop a simple one-locus, two-allele model of a haploid population of constant size including the factors of random drift, purifying selection, and random mutation. We consider different virological experiments: accumulation and reversion of deleterious mutations, competition between mutant and wild-type viruses, gene fixation, mutation frequencies at the steady state, divergence of two populations split from one population, and genetic turnover within a single population. In the first part of the review, we present all principal results in qualitative terms and illustrate them with examples obtained by computer simulation. In the second part, we derive the results formally from a diffusion equation of the Wright-Fisher type and boundary conditions, all derived from the first principles for the virus population model. We show that the leading factors and observable behavior of evolution differ significantly in three broad intervals of population size, N. The "neutral limit" is reached when N is smaller than the inverse selection coefficient. When N is larger than the inverse mutation rate per base, selection dominates and evolution is "almost" deterministic. If the selection coefficient is much larger than the mutation rate, there exists a broad interval of population sizes, in which weakly diverse populations are almost neutral while highly diverse populations are controlled by selection pressure. We discuss in detail the application of our results to human immunodeficiency virus population in vivo, sampling effects, and limitations of the model.

Search for the mechanism of genetic variation in the pro gene of human immunodeficiency virus

To study the mechanism of evolution of the human immunodeficiency virus (HIV) protease gene (pro), we analyzed a database of 213 pro sequences isolated from 11 HIV type 1-infected patients who had not been treated with protease inhibitors. Variation in pro is restricted to rare variable bases which are highly diverse and differ in location among individuals; an average variable base appears in about 16% of individuals. The average intrapatient distance per individual variable site, 27%, is similar for synonymous and nonsynonymous sites, although synonymous sites are twice as abundant. The latter observation excludes selection for diversity as an important, permanently acting factor in the evolution of pro and leaves purifying selection as the only kind of selection. Based on this, we developed a model of evolution, both within individuals and along the transmission chain, which explains variable sites as slightly deleterious mutants slowly reverting to the better-fit variant during individual infection. In the case of a single-source transmission, genetic bottlenecks at the moment of transmission effectively suppress selection, allowing mutants to accumulate along the transmission chain to high levels. However, even very rare coinfections from independent sources are, as we show, able to counteract the bottleneck effect. Therefore, there are two possible explanations for the high mutant frequency. First, the frequency of coinfection in the natural host population may be quite low. Alternatively, a strong variation of the best-adapted sequence between individuals could be caused by a combination of an immune response present in early infection and coselection.

Linkage disequilibrium test implies a large effective population number for HIV in vivo

The effective size of the HIV population in vivo, although critically important for the prediction of appearance of drug-resistant variants, is currently unknown. To address this issue, we have developed a simple virus population model, within which the relative importance of stochastic factors and purifying selection for genetic evolution differs over, at least, three broad intervals of the effective population size, with approximate boundaries given by the inverse selection coefficient and the inverse mutation rate per base per cycle. Random drift and selection dominate the smallest (stochastic) and largest (deterministic) population intervals, respectively. In the intermediate (selection-drift) interval, random drift controls weakly diverse populations, whereas strongly diverse populations are controlled by selection. To estimate the effective size of the HIV population in vivo, we tested 200 pro sequences isolated from 11 HIV-infected patients for the presence of a linkage disequilibrium effect which must exist only in small populations. This analysis demonstrated a steady-state virus population of 10(5) infected cells or more, which is either in or at the border of the deterministic regime with respect to evolution of separate bases.

Effect of population patchiness and migration rates on the adaptation and divergence of vesicular stomatitis virus quasispecies populations

The effect of migration among different isolated virus quasispecies populations on their adaptation and diversity was analysed through experimental evolution. An in vitro cell system was employed to simulate migration of vesicular stomatitis virus between isolated homogeneous host cell populations. The results clearly demonstrated a positive correlation between the migration rate and the magnitude of the mean fitness reached by the virus quasispecies populations. The results also showed, although less clearly, that fitness differences among quasispecies decreased with the magnitude of migration. These results are in close agreement with predictions of standard population genetics theory. These results can be explained in terms of the spread of beneficial mutations, originating in a single isolated quasispecies, through the entire system formed by the different quasispecies populations contained in different host cell populations.

Population biology of HIV-1 infection: viral and CD4+ T cell demographics and dynamics in lymphatic tissues

Human immunodeficiency virus-1 (HIV-1) is usually transmitted through sexual contact and in the very early stages of infection establishes a persistent infection in lymphatic tissues (LT). Virus is produced and stored at this site in a dynamic process that slowly depletes the immune system of CD4+ T cells, setting the stage for AIDS. In this review, I describe the changes in viral and CD4+ T cell populations in LT over the course of infection and after treatment. I present recent evidence that productively infected CD4+ T cells play an important role in establishing persistent infection from the onset, and that the LT are the major reservoir where virus is produced and stored on follicular dendritic cells (FDCs). I discuss the methods used to define the size of viral and CD4+ T cell populations in LT and the nature of virus-host cell interactions in vivo. These experimental approaches have identified populations of latently and chronically infected cells in which virus can elude host defenses, perpetuate infection, and escape eradication by highly active antiretroviral treatment (HAART). I discuss the dramatic impact of HAART on suppressing virus production, reducing the pool of stored virus, and restoring CD4+ T cell populations. I discuss the contributions of thymopoiesis and other renewal mechanisms, lymphatic homeostasis and trafficking to these changes in CD4+ T cell populations in LT, and conclude with a model of immune depletion and repopulation based on the limited regenerative capacity of the adult and the uncompensated losses of productively infected cells that treatment stems. The prediction of this model is that immune regeneration will be slow, variable, and partial. It is nonetheless encouraging to know that even in late stages of infection, control of active replication of HIV-1 provides an opportunity for the immune system to recover from the injuries inflicted by infection.