Interplay between liver and blood stages of Plasmodium infection dictates malaria severity via γδ T cells and IL-17-promoted stress erythropoiesis.

Plasmodium replicates within the liver prior to reaching the bloodstream and infecting red blood cells. Because clinical manifestations of malaria only arise during the blood stage of infection, a perception exists that liver infection does not impact disease pathology. By developing a murine model where the liver and blood stages of infection are uncoupled, we showed that the integration of signals from both stages dictated mortality outcomes. This dichotomy relied on liver stage-dependent activation of Vγ4+ γδ T cells. Subsequent blood stage parasite loads dictated their cytokine profiles, where low parasite loads preferentially expanded IL-17-producing γδ T cells. IL-17 drove extra-medullary erythropoiesis and concomitant reticulocytosis, which protected mice from lethal experimental cerebral malaria (ECM). Adoptive transfer of erythroid precursors could rescue mice from ECM. Modeling of γδ T cell dynamics suggests that this protective mechanism may be key for the establishment of naturally acquired malaria immunity among frequently exposed individuals.

Understanding early HIV-1 rebound dynamics following antiretroviral therapy interruption: The importance of effector cell expansion.

Most people living with HIV-1 experience rapid viral rebound once antiretroviral therapy is interrupted; however, a small fraction remain in viral remission for an extended duration. Understanding the factors that determine whether viral rebound is likely after treatment interruption can enable the development of optimal treatment regimens and therapeutic interventions to potentially achieve a functional cure for HIV-1. We built upon the theoretical framework proposed by Conway and Perelson to construct dynamic models of virus-immune interactions to study factors that influence viral rebound dynamics. We evaluated these models using viral load data from 24 individuals following antiretroviral therapy interruption. The best-performing model accurately captures the heterogeneity of viral dynamics and highlights the importance of the effector cell expansion rate. Our results show that post-treatment controllers and non-controllers can be distinguished based on the effector cell expansion rate in our models. Furthermore, these results demonstrate the potential of using dynamic models incorporating an effector cell response to understand early viral rebound dynamics post-antiretroviral therapy interruption.

Modeling the emergence of viral resistance for SARS-CoV-2 during treatment with an anti-spike monoclonal antibody.

To mitigate the loss of lives during the COVID-19 pandemic, emergency use authorization was given to several anti-SARS-CoV-2 monoclonal antibody (mAb) therapies for the treatment of mild-to-moderate COVID-19 in patients with a high risk of progressing to severe disease. Monoclonal antibodies used to treat SARS-CoV-2 target the spike protein of the virus and block its ability to enter and infect target cells. Monoclonal antibody therapy can thus accelerate the decline in viral load and lower hospitalization rates among high-risk patients with variants susceptible to mAb therapy. However, viral resistance has been observed, in some cases leading to a transient viral rebound that can be as large as 3-4 orders of magnitude. As mAbs represent a proven treatment choice for SARS-CoV-2 and other viral infections, evaluation of treatment-emergent mAb resistance can help uncover underlying pathobiology of SARS-CoV-2 infection and may also help in the development of the next generation of mAb therapies. Although resistance can be expected, the large rebounds observed are much more difficult to explain. We hypothesize replenishment of target cells is necessary to generate the high transient viral rebound. Thus, we formulated two models with different mechanisms for target cell replenishment (homeostatic proliferation and return from an innate immune response antiviral state) and fit them to data from persons with SARS-CoV-2 treated with a mAb. We showed that both models can explain the emergence of resistant virus associated with high transient viral rebounds. We found that variations in the target cell supply rate and adaptive immunity parameters have a strong impact on the magnitude or observability of the viral rebound associated with the emergence of resistant virus. Both variations in target cell supply rate and adaptive immunity parameters may explain why only some individuals develop observable transient resistant viral rebound. Our study highlights the conditions that can lead to resistance and subsequent viral rebound in mAb treatments during acute infection.

Biphasic decay of intact SHIV genomes following initiation of antiretroviral therapy complicates analysis of interventions targeting the reservoir.

The latent reservoir for HIV-1 in resting CD4+ T cells persists despite antiretroviral therapy (ART) and precludes cure. Reservoir-targeting interventions are evaluated in ART-treated macaques infected with simian immunodeficiency virus (SIV) or simian-human immunodeficiency virus (SHIV). Efficacy is determined by reservoir measurements before and after the intervention. However, most proviruses persisting in the setting of ART are defective. In addition, intact HIV-1 and SIV genomes undergo complex, multiphasic decay observable when new infection events are blocked by ART. Intervention-induced elimination of latently infected cells must be distinguished from natural decay. Here, we address these issues for SHIV. We describe an intact proviral DNA assay that allows digital counting of SHIV genomes lacking common fatal defects. We show that intact SHIV genomes in circulating CD4+ T cells undergo biphasic decay during the first year of ART, with a rapid first phase (t1/2 = 30.1 d) and a slower second phase (t1/2 = 8.1 mo) that is still more rapid that the slow decay observed in people with HIV-1 on long-term ART (t1/2 = 3.7 y). In SHIV models, most interventions are tested during 2nd phase decay. Natural 2nd phase decay must be considered in evaluating interventions as most infected cells present at this time do not become part of the stable reservoir. In addition, for interventions tested during 2nd phase decay, a caveat is that the intervention may not be equally effective in people with HIV on long-term ART whose reservoirs are dominated by latently infected cells with a slower decay rate.

Modeling dynamics of acute HIV infection incorporating density-dependent cell death and multiplicity of infection.

Understanding the dynamics of acute HIV infection can offer valuable insights into the early stages of viral behavior, potentially helping uncover various aspects of HIV pathogenesis. The standard viral dynamics model explains HIV viral dynamics during acute infection reasonably well. However, the model makes simplifying assumptions, neglecting some aspects of HIV infection. For instance, in the standard model, target cells are infected by a single HIV virion. Yet, cellular multiplicity of infection (MOI) may have considerable effects in pathogenesis and viral evolution. Further, when using the standard model, we take constant infected cell death rates, simplifying the dynamic immune responses. Here, we use four models-1) the standard viral dynamics model, 2) an alternate model incorporating cellular MOI, 3) a model assuming density-dependent death rate of infected cells and 4) a model combining (2) and (3)-to investigate acute infection dynamics in 43 people living with HIV very early after HIV exposure. We find that all models qualitatively describe the data, but none of the tested models is by itself the best to capture different kinds of heterogeneity. Instead, different models describe differing features of the dynamics more accurately. For example, while the standard viral dynamics model may be the most parsimonious across study participants by the corrected Akaike Information Criterion (AICc), we find that viral peaks are better explained by a model allowing for cellular MOI, using a linear regression analysis as analyzed by R2. These results suggest that heterogeneity in within-host viral dynamics cannot be captured by a single model. Depending on the specific aspect of interest, a corresponding model should be employed.

How robust are estimates of key parameters in standard viral dynamic models?

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

Complex decay dynamics of HIV virions, intact and defective proviruses, and 2LTR circles following initiation of antiretroviral therapy.

In persons living with HIV-1 (PLWH) who start antiretroviral therapy (ART), plasma virus decays in a biphasic fashion to below the detection limit. The first phase reflects the short half-life (<1 d) of cells that produce most of the plasma virus. The second phase represents the slower turnover (t1/2 = 14 d) of another infected cell population, whose identity is unclear. Using the intact proviral DNA assay (IPDA) to distinguish intact and defective proviruses, we analyzed viral decay in 17 PLWH initiating ART. Circulating CD4+ T cells with intact proviruses include few of the rapidly decaying first-phase cells. Instead, this population initially decays more slowly (t1/2 = 12.9 d) in a process that largely represents death or exit from the circulation rather than transition to latency. This more protracted decay potentially allows for immune selection. After ∼3 mo, the decay slope changes, and CD4+ T cells with intact proviruses decay with a half-life of 19 mo, which is still shorter than that of the latently infected cells that persist on long-term ART. Two-long-terminal repeat (2LTR) circles decay with fast and slow phases paralleling intact proviruses, a finding that precludes their use as a simple marker of ongoing viral replication. Proviruses with defects at the 5' or 3' end of the genome show equivalent monophasic decay at rates that vary among individuals. Understanding these complex early decay processes is important for correct use of reservoir assays and may provide insights into properties of surviving cells that can constitute the stable latent reservoir.

In vivo kinetics of SARS-CoV-2 infection and its relationship with a person's infectiousness.

The within-host viral kinetics of SARS-CoV-2 infection and how they relate to a person's infectiousness are not well understood. This limits our ability to quantify the impact of interventions on viral transmission. Here, we develop viral dynamic models of SARS-CoV-2 infection and fit them to data to estimate key within-host parameters such as the infected cell half-life and the within-host reproductive number. We then develop a model linking viral load (VL) to infectiousness and show a person's infectiousness increases sublinearly with VL and that the logarithm of the VL in the upper respiratory tract is a better surrogate of infectiousness than the VL itself. Using data on VL and the predicted infectiousness, we further incorporated data on antigen and RT-PCR tests and compared their usefulness in detecting infection and preventing transmission. We found that RT-PCR tests perform better than antigen tests assuming equal testing frequency; however, more frequent antigen testing may perform equally well with RT-PCR tests at a lower cost but with many more false-negative tests. Overall, our models provide a quantitative framework for inferring the impact of therapeutics and vaccines that lower VL on the infectiousness of individuals and for evaluating rapid testing strategies.

Second-Phase Hepatitis C Plasma Viral Kinetics Directly Reflects Reduced Intrahepatic Burden of Hepatitis C Virus.

BACKGROUND: Mathematical models explain how antivirals control viral infections. Hepatitis C virus (HCV) treatment results in at least 2 phases of decline in viremia. The first phase reflects clearance of rapidly produced virions. The second phase is hypothesized to derive from loss of infected cells but has been challenging to prove. METHODS: Using single-cell methods, we quantified the number of hepatitis C virus (HCV)-infected hepatocytes in liver biopsies taken before and within 7 days of initiating direct-acting antivirals (DAAs) in a double-blinded randomized controlled trial testing 2 (sofosbuvir-velpatasvir) versus 3 (sofosbuvir-velpatasvir-voxilaprevir) DAAs. RESULTS: We employed thousands of intrahepatic measurements in 10 persons with chronic genotype 1a HCV infection: median proportion of infected hepatocytes declined from 11.3% (range, 1.3%-59%) to 0.6% (range, <0.3%-5.8%), a loss of 75%-95% infected hepatocytes. Plasma viremia correlated with numbers of HCV-infected hepatocytes (r = 0.77; P < .0001). Second-phase plasma dynamics and changes in infected hepatocytes were indistinct (P = .16), demonstrating that second-phase viral dynamics derive from loss of infected cells. DAAs led to a decline in intracellular HCV RNA and interferon-stimulated gene expression (P < .05 for both). CONCLUSIONS: We proved that second-phase viral dynamics reflect decay of intrahepatic burden of HCV, partly due to clearance of HCV RNA from hepatocytes. CLINICAL TRIALS REGISTRATION: NCT02938013.

Hepatitis B e Antigen-Negative Single Hepatocyte Analysis Shows Transcriptional Silencing and Slow Decay of Infected Cells With Treatment.

BACKGROUND: Nucleos(t)ide analogues (NUCs) rarely cure chronic hepatitis B (CHB) because they do not eliminate covalently closed circular deoxyribonucleic acid, the stable replication template. In hepatitis B e antigen (HBeAg)-positive CHB during NUCs, HBV-infected cells decline slowly and are transcriptionally silenced. Whether these occur in HBeAg-negative CHB is unknown. METHODS: Using paired liver biopsies separated by 2.7-3.7 years in 4 males with HIV and HBeAg-negative CHB at both biopsies and 1 male with HIV who underwent HBeAg seroconversion between biopsies, we quantified amounts of viral nucleic acids in hundreds of individual hepatocytes. RESULTS: In the 4 persistently HBeAg-negative participants, HBV-infected hepatocytes ranged from 6.2% to 17.7% (biopsy 1) and significantly declined in 3 of 4 by biopsy 2. In the HBeAg seroconverter, the proportion was 97.4% (biopsy 1) and declined to 81.9% at biopsy 2 (P < .05). We extrapolated that HBV eradication with NUCs would take >100 years. At biopsy 1 in the persistently HBeAg-negative participants, 23%-56.8% of infected hepatocytes were transcriptionally inactive-higher than we observed in HBeAg-positive CHB-and significantly declined in 1 of 4 at biopsy 2. CONCLUSIONS: In HBeAg-negative CHB on NUCs, the negligible decline in infected hepatocytes is similar to HBeAg-positive CHB, supporting the need for more potent therapeutics to achieve functional cure.

Variant-Specific Viral Kinetics in Acute COVID-19.

Understanding variant-specific differences in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral kinetics may explain differences in transmission efficiency and provide insights on pathogenesis and prevention. We evaluated SARS-CoV-2 kinetics from nasal swabs across multiple variants (Alpha, Delta, Epsilon, Gamma) in placebo recipients of the ACTIV-2/A5401 trial. Delta variant infection led to the highest maximum viral load and shortest time from symptom onset to viral load peak. There were no significant differences in time to viral clearance across the variants. Viral decline was biphasic with first- and second-phase decays having half-lives of 11 hours and 2.5 days, respectively, with differences among variants, especially in the second phase. These results suggest that while variant-specific differences in viral kinetics exist, post-peak viral load all variants appeared to be efficiently cleared by the host. Clinical Trials Registration. NCT04518410.

Longitudinal SARS-CoV-2 seroprevalence in Portugal and antibody maintenance 12 months after infection.

During the COVID-19 pandemic, Portugal has experienced three distinct SARS-CoV-2 infection waves. We previously documented the prevalence of SARS-CoV-2 immunity, measured by specific antibodies, in September 2020, 6 months after the initial moderate wave. Here, we show the seroprevalence changes 6 months later, up to the second week of March 2021, shortly following the third wave, which was one of the most severe in the world, and 2 months following the start of the vaccination campaign. A longitudinal epidemiological study was conducted, with a stratified quota sample of the Portuguese population. Serological testing was performed, including ELISA determination of antibody class and titers. The proportion of seropositives, which was 2.2% in September 2020, rose sharply to 17.3% (95% CI: 15.8-18.8%) in March 2021. Importantly, circulating IgG and IgA antibody levels were very stable 6 months after the initial determination and up to a year after initial infection, indicating long-lasting infection immunity against SARS-CoV-2. Moreover, vaccinated people had higher IgG levels from 3 weeks post-vaccination when compared with previously infected people at the same time post-infection.

Zika virus dynamics: Effects of inoculum dose, the innate immune response and viral interference.

Experimental Zika virus infection in non-human primates results in acute viral load dynamics that can be well-described by mathematical models. The inoculum dose that would be received in a natural infection setting is likely lower than the experimental infections and how this difference affects the viral dynamics and immune response is unclear. Here we study a dataset of experimental infection of non-human primates with a range of doses of Zika virus. We develop new models of infection incorporating both an innate immune response and viral interference with that response. We find that such a model explains the data better than models with no interaction between virus and the immune response. We also find that larger inoculum doses lead to faster dynamics of infection, but approximately the same total amount of viral production.

The dynamics of simian immunodeficiency virus after depletion of CD8+ cells.

Human immunodeficiency virus infection is still one of the most important causes of morbidity and mortality in the world, with a disproportionate human and economic burden especially in poorer countries. Despite many years of intense research, an aspect that still is not well understood is what (immune) mechanisms control the viral load during the prolonged asymptomatic stage of infection. Because CD8+ T cells have been implicated in this control by multiple lines of evidence, there has been a focus on understanding the potential mechanisms of action of this immune effector population. One type of experiment used to this end has been depleting these cells with monoclonal antibodies in the simian immunodeficiency virus-macaque model and then studying the effect of that depletion on the viral dynamics. Here we review what these experiments have told us. We emphasize modeling studies to interpret the changes in viral load observed in these experiments, including discussion of alternative models, assumptions and interpretations, as well as potential future experiments.

Prevalence of SARS-CoV-2 Antibodies after First 6 Months of COVID-19 Pandemic, Portugal.

In September 2020, we tested 13,398 persons in Portugal for antibodies against severe acute respiratory syndrome coronavirus 2 by using a quota sample stratified by age and population density. We found a seroprevalence of 2.2%, 3-4 times larger than the official number of cases at the end of the first wave of the pandemic.

HIV influences clustering and intracellular replication of hepatitis C virus.

HCV and HIV coinfection is common and HIV leads to increased HCV viraemia and accelerated disease progression. However, the biological basis of this interaction remains poorly understood and little is known about the impact of HIV on HCV replication at the cellular level. We analysed HCV RNA, based on single-cell laser-capture microdissection, in liver biopsies from monoinfected (n = 4) and HCV/HIV-coinfected (n = 5) participants. HCV RNA was assayed in 3200 hepatocytes with information of spatial position. We compared HCV RNA levels and clustering properties of infection between mono- and coinfected participants, and developed a mathematical model of infection. Although the median plasma HCV RNA level and the fraction of infected cells were comparable in monoinfected (7.0 log10 IU/mL and ~ 30%) and coinfected (7.3 log10 IU/mL and ~ 40%) participants, the median HCV RNA per infected hepatocyte in monoinfected (2.8IU) was significantly lower than in coinfected (8.2IU) participants (p = .03). Clustering of infected cells was more prominent in monoinfected participants (91% of samples) than in coinfected participants (~48%), p = .0045, suggesting that spatial spread may be influenced by HIV coinfection. Interestingly, when clustering does occur, the size of clusters is similar in both types of infection. A mathematical model of infection suggested that HIV allows higher intracellular accumulation of HCV RNA by impeding the export of HCV RNA. Our observations show that HIV coinfection impacts intracellular accumulation of HCV RNA and the clustering of HCV-infected cells, but to a less extent the fraction of HCV-infected cells.

Serum lipids and prostate cancer.

BACKGROUND: Conflicting results are found in the literature relating serum lipids levels and prostate cancer. Some results imply a relationship between them; others contradict this association. The purpose of this study was to investigate a possible association between serum lipids levels and prostate cancer, at time of diagnosis. METHODS: We measured serum levels of total cholesterol, HDL cholesterol, LDL cholesterol, and triglycerides in 237 patients submitted to a prostate biopsy, with PSA between 2 and 10 ng/ml. Patients without cancer at biopsy were used as controls, and the others were considered as cases. No information about lipid-lowering therapy, including statins, was available neither in cases nor in controls. Cases were divided into risk groups, according to the disease severity, based on staging. Lipids levels were compared between groups, using parametric and nonparametric tests. Logistic regression analysis and odds ratios were calculated. RESULTS: LDL and total cholesterol levels were lower in patients with cancer, with the difference being statistically significant for LDL cholesterol (p = 0.010) and borderline for total cholesterol (p = 0.050). No significant differences were found between the several risk groups. Odds ratios for low LDL cholesterol (<130 mg/dl) and low total cholesterol (<200 mg/dl), with prostate cancer as the outcome, were 1.983 and 1.703, respectively. There were no significant differences between cases and controls for the other lipids. CONCLUSION: Lower LDL cholesterol (<130 mg/dl) and lower total cholesterol (<200 mg/dl) serum levels seem to associate with prostate cancer, at time of diagnosis.

Probabilistic control of HIV latency and transactivation by the Tat gene circuit.

The reservoir of HIV latently infected cells is the major obstacle for eradication of HIV infection. The "shock-and-kill" strategy proposed earlier aims to reduce the reservoir by activating cells out of latency. While the intracellular HIV Tat gene circuit is known to play important roles in controlling latency and its transactivation in HIV-infected cells, the detailed control mechanisms are not well understood. Here we study the mechanism of probabilistic control of the latent and the transactivated cell phenotypes of HIV-infected cells. We reconstructed the probability landscape, which is the probability distribution of the Tat gene circuit states, by directly computing the exact solution of the underlying chemical master equation. Results show that the Tat circuit exhibits a clear bimodal probability landscape (i.e., there are two distinct probability peaks, one associated with the latent cell phenotype and the other with the transactivated cell phenotype). We explore potential modifications to reactions in the Tat gene circuit for more effective transactivation of latent cells (i.e., the shock-and-kill strategy). Our results suggest that enhancing Tat acetylation can dramatically increase Tat and viral production, while increasing the Tat-transactivation response binding affinity can transactivate latent cells more rapidly than other manipulations. Our results further explored the "block and lock" strategy toward a functional cure for HIV. Overall, our study demonstrates a general approach toward discovery of effective therapeutic strategies and druggable targets by examining control mechanisms of cell phenotype switching via exactly computed probability landscapes of reaction networks.

Superinfection and cure of infected cells as mechanisms for hepatitis C virus adaptation and persistence.

RNA viruses exist as a genetically diverse quasispecies with extraordinary ability to adapt to abrupt changes in the host environment. However, the molecular mechanisms that contribute to their rapid adaptation and persistence in vivo are not well studied. Here, we probe hepatitis C virus (HCV) persistence by analyzing clinical samples taken from subjects who were treated with a second-generation HCV protease inhibitor. Frequent longitudinal viral load determinations and large-scale single-genome sequence analyses revealed rapid antiviral resistance development, and surprisingly, dynamic turnover of dominant drug-resistant mutant populations long after treatment cessation. We fitted mathematical models to both the viral load and the viral sequencing data, and the results provided strong support for the critical roles that superinfection and cure of infected cells play in facilitating the rapid turnover and persistence of viral populations. More broadly, our results highlight the importance of considering viral dynamics and competition at the intracellular level in understanding rapid viral adaptation. Thus, we propose a theoretical framework integrating viral and molecular mechanisms to explain rapid viral evolution, resistance, and persistence despite antiviral treatment and host immune responses.