Evolution during primary HIV infection does not require adaptive immune selection.

Modern HIV research depends crucially on both viral sequencing and population measurements. To directly link mechanistic biological processes and evolutionary dynamics during HIV infection, we developed multiple within-host phylodynamic models of HIV primary infection for comparative validation against viral load and evolutionary dynamics data. The optimal model of primary infection required no positive selection, suggesting that the host adaptive immune system reduces viral load but surprisingly does not drive observed viral evolution. Rather, the fitness (infectivity) of mutant variants is drawn from an exponential distribution in which most variants are slightly less infectious than their parents (nearly neutral evolution). This distribution was not largely different from either in vivo fitness distributions recorded beyond primary infection or in vitro distributions that are observed without adaptive immunity, suggesting the intrinsic viral fitness distribution may drive evolution. Simulated phylogenetic trees also agree with independent data and illuminate how phylogenetic inference must consider viral and immune-cell population dynamics to gain accurate mechanistic insights.

Relationship between CD4 T cell turnover, cellular differentiation and HIV persistence during ART.

The precise role of CD4 T cell turnover in maintaining HIV persistence during antiretroviral therapy (ART) has not yet been well characterized. In resting CD4 T cell subpopulations from 24 HIV-infected ART-suppressed and 6 HIV-uninfected individuals, we directly measured cellular turnover by heavy water labeling, HIV reservoir size by integrated HIV-DNA (intDNA) and cell-associated HIV-RNA (caRNA), and HIV reservoir clonality by proviral integration site sequencing. Compared to HIV-negatives, ART-suppressed individuals had similar fractional replacement rates in all subpopulations, but lower absolute proliferation rates of all subpopulations other than effector memory (TEM) cells, and lower plasma IL-7 levels (p = 0.0004). Median CD4 T cell half-lives decreased with cell differentiation from naïve to TEM cells (3 years to 3 months, p<0.001). TEM had the fastest replacement rates, were most highly enriched for intDNA and caRNA, and contained the most clonal proviral expansion. Clonal proviruses detected in less mature subpopulations were more expanded in TEM, suggesting that they were maintained through cell differentiation. Earlier ART initiation was associated with lower levels of intDNA, caRNA and fractional replacement rates. In conclusion, circulating integrated HIV proviruses appear to be maintained both by slow turnover of immature CD4 subpopulations, and by clonal expansion as well as cell differentiation into effector cells with faster replacement rates.

Optimizing clinical dosing of combination broadly neutralizing antibodies for HIV prevention.

Broadly neutralizing antibodies (bNAbs) are promising agents to prevent HIV infection and achieve HIV remission without antiretroviral therapy (ART). As with ART, bNAb combinations are likely needed to cover HIV's extensive diversity. Not all bNAbs are identical in terms of their breadth, potency, and in vivo longevity (half-life). Given these differences, it is important to optimally select the composition, or dose ratio, of combination bNAb therapies for future clinical studies. We developed a model that synthesizes 1) pharmacokinetics, 2) potency against a wide HIV diversity, 3) interaction models for how drugs work together, and 4) correlates that translate in vitro potency to clinical protection. We found optimization requires drug-specific balances between potency, longevity, and interaction type. As an example, tradeoffs between longevity and potency are shown by comparing a combination therapy to a bi-specific antibody (a single protein merging both bNAbs) that takes the better potency but the worse longevity of the two components. Then, we illustrate a realistic dose ratio optimization of a triple combination of VRC07, 3BNC117, and 10-1074 bNAbs. We apply protection estimates derived from both a non-human primate (NHP) challenge study meta-analysis and the human antibody mediated prevention (AMP) trials. In both cases, we find a 2:1:1 dose emphasizing VRC07 is nearly optimal. Our approach can be immediately applied to optimize the next generation of combination antibody prevention and cure studies.

A SARS-CoV-2 vaccine candidate would likely match all currently circulating variants.

The magnitude of the COVID-19 pandemic underscores the urgency for a safe and effective vaccine. Many vaccine candidates focus on the Spike protein, as it is targeted by neutralizing antibodies and plays a key role in viral entry. Here we investigate the diversity seen in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequences and compare it to the sequence on which most vaccine candidates are based. Using 18,514 sequences, we perform phylogenetic, population genetics, and structural bioinformatics analyses. We find limited diversity across SARS-CoV-2 genomes: Only 11 sites show polymorphisms in >5% of sequences; yet two mutations, including the D614G mutation in Spike, have already become consensus. Because SARS-CoV-2 is being transmitted more rapidly than it evolves, the viral population is becoming more homogeneous, with a median of seven nucleotide substitutions between genomes. There is evidence of purifying selection but little evidence of diversifying selection, with substitution rates comparable across structural versus nonstructural genes. Finally, the Wuhan-Hu-1 reference sequence for the Spike protein, which is the basis for different vaccine candidates, matches optimized vaccine inserts, being identical to an ancestral sequence and one mutation away from the consensus. While the rapid spread of the D614G mutation warrants further study, our results indicate that drift and bottleneck events can explain the minimal diversity found among SARS-CoV-2 sequences. These findings suggest that a single vaccine candidate should be efficacious against currently circulating lineages.

Reliability of Self-Sampling for Accurate Assessment of Respiratory Virus Viral and Immunologic Kinetics.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic demonstrates the need for accurate and convenient approaches to diagnose and therapeutically monitor respiratory viral infections. We demonstrated that self-sampling with mid-nasal foam swabs is well-tolerated and provides quantitative viral output concordant with flocked swabs. Using longitudinal home-based self-sampling, we demonstrate that nasal cytokine levels correlate and cluster according to immune cell of origin. Periods of stable viral loads are followed by rapid elimination, which could be coupled with cytokine expansion and contraction. Nasal foam swab self-sampling at home provides a precise, mechanistic readout of respiratory virus shedding and local immune responses.

Dynamics of HIV DNA reservoir seeding in a cohort of superinfected Kenyan women.

A reservoir of HIV-infected cells that persists despite suppressive antiretroviral therapy (ART) is the source of viral rebound upon ART cessation and the major barrier to a cure. Understanding reservoir seeding dynamics will help identify the best timing for HIV cure strategies. Here we characterize reservoir seeding using longitudinal samples from before and after ART initiation in individuals who sequentially became infected with genetically distinct HIV variants (superinfected). We previously identified cases of superinfection in a cohort of Kenyan women, and the dates of both initial infection and superinfection were determined. Six women, superinfected 0.2-5.2 years after initial infection, were subsequently treated with ART 5.4-18.0 years after initial infection. We performed next-generation sequencing of HIV gag and env RNA from plasma collected during acute infection as well as every ~2 years thereafter until ART initiation, and of HIV DNA from PBMCs collected 0.9-4.8 years after viral suppression on ART. We assessed phylogenetic relationships between HIV DNA reservoir sequences and longitudinal plasma RNA sequences prior to ART, to determine proportions of initial and superinfecting variants in the reservoir. The proportions of initial and superinfection lineage variants present in the HIV DNA reservoir were most similar to the proportions present in HIV RNA immediately prior to ART initiation. Phylogenetic analysis confirmed that the majority of HIV DNA reservoir sequences had the smallest pairwise distance to RNA sequences from timepoints closest to ART initiation. Our data suggest that while reservoir cells are created throughout pre-ART infection, the majority of HIV-infected cells that persist during ART entered the reservoir near the time of ART initiation. We estimate the half-life of pre-ART DNA reservoir sequences to be ~25 months, which is shorter than estimated reservoir decay rates during suppressive ART, implying continual decay and reseeding of the reservoir up to the point of ART initiation.

Molecular dating and viral load growth rates suggested that the eclipse phase lasted about a week in HIV-1 infected adults in East Africa and Thailand.

Most HIV-1 infected individuals do not know their infection dates. Precise infection timing is crucial information for studies that document transmission networks or drug levels at infection. To improve infection timing, we used the prospective RV217 cohort where the window when plasma viremia becomes detectable is narrow: the last negative visit occurred a median of four days before the first detectable HIV-1 viremia with an RNA test, referred below as diagnosis. We sequenced 1,280 HIV-1 genomes from 39 participants at a median of 4, 32 and 170 days post-diagnosis. HIV-1 infections were dated by using sequence-based methods and a viral load regression method. Bayesian coalescent and viral load regression estimated that infections occurred a median of 6 days prior to diagnosis (IQR: 9-3 and 11-4 days prior, respectively). Poisson-Fitter, which analyzes the distribution of hamming distances among sequences, estimated a median of 7 days prior to diagnosis (IQR: 15-4 days) based on sequences sampled 4 days post-diagnosis, but it did not yield plausible results using sequences sampled at 32 days. Fourteen participants reported a high-risk exposure event at a median of 8 days prior to diagnosis (IQR: 12 to 6 days prior). These different methods concurred that HIV-1 infection occurred about a week before detectable viremia, corresponding to 20 days (IQR: 34-15 days) before peak viral load. Together, our methods comparison helps define a framework for future dating studies in early HIV-1 infection.

Distinguishing Nanoparticle Aggregation from Viscosity Changes in MPS/MSB Detection of Biomarkers.

Magnetic particle spectroscopy (MPS) in the Brownian relaxation regime, also termed magnetic spectroscopy of Brownian motion (MSB), can detect and quantitate very low, sub-nanomolar concentrations of molecular biomarkers. MPS/MSB uses the harmonics of the magnetization induced by a small, low-frequency oscillating magnetic field to provide quantitative information about the magnetic nanoparticles' (mNPs') microenvironment. A key application uses antibody-coated mNPs to produce biomarker-mediated aggregation that can be detected using MPS/MSB. However, relaxation changes can also be caused by viscosity changes. To address this challenge, we propose a metric that can distinguish between aggregation and viscosity. Viscosity changes scale the MPS/MSB harmonic ratios with a constant multiplier across all applied field frequencies. The change in viscosity is exactly equal to the multiplier with generality, avoiding the need to understand the signal explicitly. This simple scaling relationship is violated when particles aggregate. Instead, a separate multiplier must be used for each frequency. The standard deviation of the multipliers over frequency defines a metric isolating viscosity (zero standard deviation) from aggregation (non-zero standard deviation). It increases monotonically with biomarker concentration. We modeled aggregation and simulated the MPS/MSB signal changes resulting from aggregation and viscosity changes. MPS/MSB signal changes were also measured experimentally using 100 nm iron-oxide mNPs in solutions with different viscosities (modulated by glycerol concentration) and with different levels of aggregation (modulated by concanavalin A linker concentrations). Experimental and simulation results confirmed that viscosity changes produced small changes in the standard deviation and aggregation produced larger values of standard deviation. This work overcomes a key barrier to using MPS/MSB to detect biomarkers in vivo with variable tissue viscosity.

Mathematical modeling to reveal breakthrough mechanisms in the HIV Antibody Mediated Prevention (AMP) trials.

The ongoing Antibody Mediated Prevention (AMP) trials will uncover whether passive infusion of the broadly neutralizing antibody (bNAb) VRC01 can protect against HIV acquisition. Previous statistical simulations indicate these trials may be partially protective. In that case, it will be crucial to identify the mechanism of breakthrough infections. To that end, we developed a mathematical modeling framework to simulate the AMP trials and infer the breakthrough mechanisms using measurable trial outcomes. This framework combines viral dynamics with antibody pharmacokinetics and pharmacodynamics, and will be generally applicable to forthcoming bNAb prevention trials. We fit our model to human viral load data (RV217). Then, we incorporated VRC01 neutralization using serum pharmacokinetics (HVTN 104) and in vitro pharmacodynamics (LANL CATNAP database). We systematically explored trial outcomes by reducing in vivo potency and varying the distribution of sensitivity to VRC01 in circulating strains. We found trial outcomes could be used in a clinical trial regression model (CTRM) to reveal whether partially protective trials were caused by large fractions of VRC01-resistant (IC50>50 µg/mL) circulating strains or rather a global reduction in VRC01 potency against all strains. The former mechanism suggests the need to enhance neutralizing antibody breadth; the latter suggests the need to enhance VRC01 delivery and/or in vivo binding. We will apply the clinical trial regression model to data from the completed trials to help optimize future approaches for passive delivery of anti-HIV neutralizing antibodies.

Autologous Stem Cell Transplantation Disrupts Adaptive Immune Responses during Rebound Simian/Human Immunodeficiency Virus Viremia.

Primary HIV-1 infection induces a virus-specific adaptive/cytolytic immune response that impacts the plasma viral load set point and the rate of progression to AIDS. Combination antiretroviral therapy (cART) suppresses plasma viremia to undetectable levels that rebound upon cART treatment interruption. Following cART withdrawal, the memory component of the virus-specific adaptive immune response may improve viral control compared to primary infection. Here, using primary infection and treatment interruption data from macaques infected with simian/human immunodeficiency virus (SHIV), we observe a lower peak viral load but an unchanged viral set point during viral rebound. The addition of an autologous stem cell transplant before cART withdrawal alters viral dynamics: we found a higher rebound set point but similar peak viral loads compared to the primary infection. Mathematical modeling of the data that accounts for fundamental immune parameters achieves excellent fit to heterogeneous viral loads. Analysis of model output suggests that the rapid memory immune response following treatment interruption does not ultimately lead to better viral containment. Transplantation decreases the durability of the adaptive immune response following cART withdrawal and viral rebound. Our model's results highlight the impact of the endogenous adaptive immune response during primary SHIV infection. Moreover, because we capture adaptive immune memory and the impact of transplantation, this model will provide insight into further studies of cure strategies inspired by the Berlin patient.IMPORTANCE HIV patients who interrupt combination antiretroviral therapy (cART) eventually experience viral rebound, the return of viral loads to pretreatment levels. However, the "Berlin patient" remained free of HIV rebound over a decade after stopping cART. His cure is attributed to leukemia treatment that included an HIV-resistant stem cell transplant. Inspired by this case, we studied the impact of stem cell transplantation in a macaque simian/HIV (SHIV) system. Using a mechanistic mathematical model, we found that while primary infection generates an adaptive immune memory response, stem cell transplantation disrupts this learned immunity. The results have implications for HIV cure regimens based on stem cell transplantation.

Approaches for modeling magnetic nanoparticle dynamics.

Magnetic nanoparticles are useful biological probes as well as therapeutic agents. Several approaches have been used to model nanoparticle magnetization dynamics for both Brownian as well as Neel rotation. Magnetizations are often of interest and can be compared with experimental results. Here we summarize these approaches, including the Stoner-Wohlfarth approach and stochastic approaches including thermal fluctuations. Non-equilibrium-related temperature effects can be described by a distribution function approach (Fokker-Planck equation) or a stochastic differential equation (Langevin equation). Approximate models in several regimes can be derived from these general approaches to simplify implementation.

Magnetic nanoparticle sensing: decoupling the magnetization from the excitation field.

Remote sensing of magnetic nanoparticles has exciting applications for magnetic nanoparticle hyperthermia and molecular detection. We introduce, simulate, and experimentally demonstrate an innovation-a sensing coil that is geometrically decoupled from the excitation field-for magnetic nanoparticle spectroscopy that increases the flexibility and capabilities of remote detection. The decoupling enhances the sensitivity absolutely; to small amounts of nanoparticles, and relatively; to small changes in the nanoparticle dynamics. We adapt a previous spectroscopic method that measures the relaxation time of nanoparticles and demonstrate a new measurement of nanoparticle temperature that could potentially be used concurrently during hyperthermia.

Mild HIV-specific selective forces overlaying natural CD4+ T cell dynamics explain the clonality and decay dynamics of HIV reservoir cells.

The latent reservoir of HIV persists for decades in people living with HIV (PWH) on antiretroviral therapy (ART). To determine if persistence arises from the natural dynamics of memory CD4+ T cells harboring HIV, we compared the clonal dynamics of HIV proviruses to that of memory CD4+ T cell receptors (TCRß) from the same PWH and from HIV-seronegative people. We show that clonal dominance of HIV proviruses and antigen-specific CD4+ T cells are similar but that the field's understanding of the persistence of the less clonally dominant reservoir is significantly limited by undersampling. We demonstrate that increasing reservoir clonality over time and differential decay of intact and defective proviruses cannot be explained by mCD4+ T cell kinetics alone. Finally, we develop a stochastic model of TCRß and proviruses that recapitulates experimental observations and suggests that HIV-specific negative selection mediates approximately 6% of intact and 2% of defective proviral clearance. Thus, HIV persistence is mostly, but not entirely, driven by natural mCD4+ T cell kinetics.

Quantification of magnetic nanoparticles with low frequency magnetic fields: compensating for relaxation effects.

Quantifying the number of nanoparticles present in tissue is central to many in vivo and in vitro applications. Magnetic nanoparticles can be detected with high sensitivity both in vivo and in vitro using the harmonics of their magnetization produced in a sinusoidal magnetic field. However, relaxation effects damp the magnetic harmonics rendering them of limited use in quantification. We show that an accurate measure of the number of nanoparticles can be made by correcting for relaxation effects. Correction for relaxation reduced errors of 50% for larger nanoparticles in high relaxation environments to 2%. The result is a method of nanoparticle quantification suitable for in vivo and in vitro applications including histopathology assays, quantitative imaging, drug delivery and thermal therapy preparation.

Impact of misclassified defective proviruses on HIV reservoir measurements.

Most proviruses persisting in people living with HIV (PWH) on antiretroviral therapy (ART) are defective. However, rarer intact proviruses almost always reinitiate viral rebound if ART stops. Therefore, assessing therapies to prevent viral rebound hinges on specifically quantifying intact proviruses. We evaluated the same samples from 10 male PWH on ART using the two-probe intact proviral DNA assay (IPDA) and near full length (nfl) Q4PCR. Both assays admitted similar ratios of intact to total HIV DNA, but IPDA found ~40-fold more intact proviruses. Neither assay suggested defective proviruses decay over 10 years. However, the mean intact half-lives were different: 108 months for IPDA and 65 months for Q4PCR. To reconcile this difference, we modeled additional longitudinal IPDA data and showed that decelerating intact decay could arise from very long-lived intact proviruses and/or misclassified defective proviruses: slowly decaying defective proviruses that are intact in IPDA probe locations (estimated up to 5%, in agreement with sequence library based predictions). The model also demonstrates how misclassification can lead to underestimated efficacy of therapies that exclusively reduce intact proviruses. We conclude that sensitive multi-probe assays combined with specific nfl-verified assays would be optimal to document absolute and changing levels of intact HIV proviruses.

Estimating the contribution of CD4 T cell subset proliferation and differentiation to HIV persistence.

Persistence of HIV in people living with HIV (PWH) on suppressive antiretroviral therapy (ART) has been linked to physiological mechanisms of CD4+ T cells. Here, in the same 37 male PWH on ART we measure longitudinal kinetics of HIV DNA and cell turnover rates in five CD4 cell subsets: naïve (TN), stem-cell- (TSCM), central- (TCM), transitional- (TTM), and effector-memory (TEM). HIV decreases in TTM and TEM but not in less-differentiated subsets. Cell turnover is ~10 times faster than HIV clearance in memory subsets, implying that cellular proliferation consistently creates HIV DNA. The optimal mathematical model for these integrated data sets posits HIV DNA also passages between CD4 cell subsets via cellular differentiation. Estimates are heterogeneous, but in an average participant's year ~10 (in TN and TSCM) and ~104 (in TCM, TTM, TEM) proviruses are generated by proliferation while ~103 proviruses passage via cell differentiation (per million CD4). In simulations, therapies blocking proliferation and/or enhancing differentiation could reduce HIV DNA by 1-2 logs over 3 years. In summary, HIV exploits cellular proliferation and differentiation to persist during ART but clears faster in more proliferative/differentiated CD4 cell subsets and the same physiological mechanisms sustaining HIV might be temporarily modified to reduce it.

High monoclonal neutralization titers reduced breakthrough HIV-1 viral loads in the Antibody Mediated Prevention trials.

The Antibody Mediated Prevention (AMP) trials (NCT02716675 and NCT02568215) demonstrated that passive administration of the broadly neutralizing monoclonal antibody VRC01 could prevent some HIV-1 acquisition events. Here, we use mathematical modeling in a post hoc analysis to demonstrate that VRC01 influenced viral loads in AMP participants who acquired HIV. Instantaneous inhibitory potential (IIP), which integrates VRC01 serum concentration and VRC01 sensitivity of acquired viruses in terms of both IC50 and IC80, follows a dose-response relationship with first positive viral load (p = 0.03), which is particularly strong above a threshold of IIP = 1.6 (r = -0.6, p = 2e-4). Mathematical modeling reveals that VRC01 activity predicted from in vitro IC80s and serum VRC01 concentrations overestimates in vivo neutralization by 600-fold (95% CI: 300-1200). The trained model projects that even if future therapeutic HIV trials of combination monoclonal antibodies do not always prevent acquisition, reductions in viremia and reservoir size could be expected.

Multi-scale modelling reveals that early super-spreader events are a likely contributor to novel variant predominance.

The emergence of new SARS-CoV-2 variants of concern (VOC) has hampered international efforts to contain the COVID-19 pandemic. VOCs have been characterized to varying degrees by higher transmissibility, worse infection outcomes and evasion of vaccine and infection-induced immunologic memory. VOCs are hypothesized to have originated from animal reservoirs, communities in regions with low surveillance and/or single individuals with poor immunologic control of the virus. Yet, the factors dictating which variants ultimately predominate remain incompletely characterized. Here we present a multi-scale model of SARS-CoV-2 dynamics that describes population spread through individuals whose viral loads and numbers of contacts (drawn from an over-dispersed distribution) are both time-varying. This framework allows us to explore how super-spreader events (SSE) (defined as greater than five secondary infections per day) contribute to variant emergence. We find stochasticity remains a powerful determinant of predominance. Variants that predominate are more likely to be associated with higher infectiousness, an SSE early after variant emergence and ongoing decline of the current dominant variant. Additionally, our simulations reveal that most new highly infectious variants that infect one or a few individuals do not achieve permanence in the population. Consequently, interventions that reduce super-spreading may delay or mitigate emergence of VOCs.

Improving vaccination coverage and offering vaccine to all school-age children allowed uninterrupted in-person schooling in King County, WA: Modeling analysis.

The rapid spread of highly transmissible SARS-CoV-2 variants combined with slowing pace of vaccination in Fall 2021 created uncertainty around the future trajectory of the epidemic in King County, Washington, USA. We analyzed the benefits of offering vaccination to children ages 5-11 and expanding the overall vaccination coverage using mathematical modeling. We adapted a mathematical model of SARS-CoV-2 transmission, calibrated to data from King County, Washington, to simulate scenarios of vaccinating children aged 5-11 with different starting dates and different proportions of physical interactions (PPI) in schools being restored. Dynamic social distancing was implemented in response to changes in weekly hospitalizations. Reduction of hospitalizations and estimated time under additional social distancing measures are reported over the 2021-2022 school year. In the scenario with 85% vaccination coverage of 12+ year-olds, offering early vaccination to children aged 5-11 with 75% PPI was predicted to prevent 756 (median, IQR 301-1434) hospitalizations cutting youth hospitalizations in half compared to no vaccination and largely reducing the need for additional social distancing measures over the school year. If, in addition, 90% overall vaccination coverage was reached, 60% of remaining hospitalizations would be averted and the need for increased social distancing would almost certainly be avoided. Our work suggests that uninterrupted in-person schooling in King County was partly possible because reasonable precaution measures were taken at schools to reduce infectious contacts. Rapid vaccination of all school-aged children provides meaningful reduction of the COVID-19 health burden over this school year but only if implemented early. It remains critical to vaccinate as many people as possible to limit the morbidity and mortality associated with future epidemic waves.

Substantial uneven proliferation of CD4+ T cells during recovery from acute HIV infection is sufficient to explain the observed expanded clones in the HIV reservoir.

The HIV reservoir is a population of 1-10 million anatomically dispersed, latently infected memory CD4+ T cells in which HIV DNA is quiescently integrated into human chromosomal DNA. When antiretroviral therapy (ART) is stopped and HIV replication initiates in one of these cells, systemic viral spread resumes, rekindling progression to AIDS. Therefore, HIV latency prevents cure. The detection of many populations of identical HIV sequences at unique integration sites implicates CD4+ T cell proliferation as the critical driver of reservoir sustainment after a prolonged period of effective ART. Initial reservoir formation occurs during the first week of primary infection usually before ART is started. While empirical data indicates that both de novo infection and cellular proliferation generate latently infected cells during early untreated infection, it is not known which of these mechanisms is predominant. We developed a mathematical model that recapitulates the profound depletion and brisk recovery of CD4+ T cells, reservoir creation, and viral load trajectory during primary HIV infection. We extended the model to stochastically simulate individual HIV reservoir clones. This model predicts the first detection of HIV infected clones approximately 5 weeks after infection as has recently been shown in vivo and suggests that substantial, uneven proliferation among clones during the recovery from CD4+ lymphopenia is the most plausible explanation for the observed clonal reservoir distribution during the first year of infection.