Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation to analyze extracellular viral markers.

Chronic infection with hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. Thus, quantifying and understanding the dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, such study requires repeated liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because liver biopsy is potentially morbid and not common during hepatitis B treatment. We here aimed to develop a noninvasive method for quantifying cccDNA in the liver using surrogate markers in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations, integrates experimental data from in vitro and in vivo investigations. By applying this model, we roughly predicted the amount and dynamics of intrahepatic cccDNA within a certain range using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The noninvasive quantification of cccDNA using our proposed method holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

Age-related changes in the hematopoietic stem cell pool revealed via quantifying the balance of symmetric and asymmetric divisions.

Hematopoietic stem cells (HSCs) are somatic stem cells that continuously generate lifelong supply of blood cells through a balance of symmetric and asymmetric divisions. It is well established that the HSC pool increases with age. However, not much is known about the underlying cause for these observed changes. Here, using a novel method combining single-cell ex vivo HSC expansion with mathematical modeling, we quantify HSC division types (stem cell-stem cell (S-S) division, stem cell-progenitor cell (S-P) division, and progenitor cell-progenitor cell (P-P) division) as a function of the aging process. Our time-series experiments reveal how changes in these three modes of division can explain the increase in HSC numbers with age. Contrary to the popular notion that HSCs divide predominantly through S-P divisions, we show that S-S divisions are predominant throughout the lifespan of the animal, thereby expanding the HSC pool. We, therefore, provide a novel mathematical model-based experimental validation for reflecting HSC dynamics in vivo.

Anti-HTLV-1 immunity combined with proviral load as predictive biomarkers for adult T-cell leukemia-lymphoma.

Human T-cell leukemia virus type 1 (HTLV-1) establishes chronic infection in humans and induces a T-cell malignancy called adult T-cell leukemia-lymphoma (ATL) and several inflammatory diseases such as HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). Persistent HTLV-1 infection is established under the pressure of host immunity, and therefore the immune response against HTLV-1 is thought to reflect the status of the disease it causes. Indeed, it is known that cellular immunity against viral antigens is suppressed in ATL patients compared to HAM/TSP patients. In this study, we show that profiling the humoral immunity to several HTLV-1 antigens, such as Gag, Env, and Tax, and measuring proviral load are useful tools for classifying disease status and predicting disease development. Using targeted sequencing, we found that several carriers whom this profiling method predicted to be at high risk for developing ATL indeed harbored driver mutations of ATL. The clonality of HTLV-1-infected cells in those carriers was still polyclonal; it is consistent with an early stage of leukemogenesis. Furthermore, this study revealed significance of anti-Gag proteins to predict high risk group in HTLV-1 carriers. Consistent with this finding, anti-Gag cytotoxic T lymphocytes (CTLs) were increased in patients who received hematopoietic stem cell transplantation and achieved remission state, indicating the significance of anti-Gag CTLs for disease control. Our findings suggest that our strategy that combines anti-HTLV-1 antibodies and proviral load may be useful for prediction of the development of HTLV-1-associated diseases.

Multiscale modeling of HBV infection integrating intra- and intercellular viral propagation for analyzing extracellular viral markers.

Chronic infection of hepatitis B virus (HBV) is caused by the persistence of closed circular DNA (cccDNA) in the nucleus of infected hepatocytes. Despite available therapeutic anti-HBV agents, eliminating the cccDNA remains challenging. The quantifying and understanding dynamics of cccDNA are essential for developing effective treatment strategies and new drugs. However, it requires a liver biopsy to measure the intrahepatic cccDNA, which is basically not accepted because of the ethical aspect. We here aimed to develop a non-invasive method for quantifying cccDNA in the liver using surrogate markers present in peripheral blood. We constructed a multiscale mathematical model that explicitly incorporates both intracellular and intercellular HBV infection processes. The model, based on age-structured partial differential equations (PDEs), integrates experimental data from in vitro and in vivo investigations. By applying this model, we successfully predicted the amount and dynamics of intrahepatic cccDNA using specific viral markers in serum samples, including HBV DNA, HBsAg, HBeAg, and HBcrAg. Our study represents a significant step towards advancing the understanding of chronic HBV infection. The non-invasive quantification of cccDNA using our proposed methodology holds promise for improving clinical analyses and treatment strategies. By comprehensively describing the interactions of all components involved in HBV infection, our multiscale mathematical model provides a valuable framework for further research and the development of targeted interventions.

Smoking enhances the expression of angiotensin-converting enzyme 2 involved in the efficiency of severe acute respiratory syndrome coronavirus 2 infection.

Smoking is one of the risk factors most closely related to the severity of coronavirus disease 2019 (COVID-19). However, the relationship between smoking history and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infectivity is unknown. In this study, we evaluated the ACE2 expression level in the lungs of current smokers, ex-smokers, and nonsmokers. The ACE2 expression level of ex-smokers who smoked cigarettes until recently (cessation period shorter than 6 months) was higher than that of nonsmokers and ex-smokers with a long history of nonsmoking (cessation period longer than 6 months). We also showed that the efficiency of SARS-CoV-2 infection was enhanced in a manner dependent on the angiotensin-converting enzyme 2 (ACE2) expression level. Using RNA-seq analysis on the lungs of smokers, we identified that the expression of inflammatory signaling genes was correlated with ACE2 expression. Notably, with increasing duration of smoking cessation among ex-smokers, not only ACE2 expression level but also the expression levels of inflammatory signaling genes decreased. These results indicated that smoking enhances the expression levels of ACE2 and inflammatory signaling genes. Our data suggest that the efficiency of SARS-CoV-2 infection is enhanced by smoking-mediated upregulation of ACE2 expression level.

[Digital transformation of COVID-19 research].

In a current life sciences research, we are in an era in which advanced technology emerging and utilize big data. Data-driven approaches such as machine learnings play an important role to analyze these datasets. However, limited clinical (time-course) datasets are available for infectious diseases, cancer, and other diseases. Especially in the case of emerging infectious disease outbreaks, clinical data obtained from a limited number of cases must be used to develop treatment strategies and public health policies. This means that many clinical data are not big data, which often makes the application of data-driven approaches difficult. In this paper, we mainly apply a mathematical model-based approach to the clinical data of COVID-19 and discuss how biologically important information can be extracted from the limited data and how they can benefit society.

Incomplete antiviral treatment may induce longer durations of viral shedding during SARS-CoV-2 infection.

The duration of viral shedding is determined by a balance between de novo infection and removal of infected cells. That is, if infection is completely blocked with antiviral drugs (100% inhibition), the duration of viral shedding is minimal and is determined by the length of virus production. However, some mathematical models predict that if infected individuals are treated with antiviral drugs with efficacy below 100%, viral shedding may last longer than without treatment because further de novo infections are driven by entry of the virus into partially protected, uninfected cells at a slower rate. Using a simple mathematical model, we quantified SARS-CoV-2 infection dynamics in non-human primates and characterized the kinetics of viral shedding. We counterintuitively found that treatments initiated early, such as 0.5 d after virus inoculation, with intermediate to relatively high efficacy (30-70% inhibition of virus replication) yield a prolonged duration of viral shedding (by about 6.0 d) compared with no treatment.

Should a viral genome stay in the host cell or leave? A quantitative dynamics study of how hepatitis C virus deals with this dilemma.

Virus proliferation involves gene replication inside infected cells and transmission to new target cells. Once positive-strand RNA virus has infected a cell, the viral genome serves as a template for copying ("stay-strategy") or is packaged into a progeny virion that will be released extracellularly ("leave-strategy"). The balance between genome replication and virion release determines virus production and transmission efficacy. The ensuing trade-off has not yet been well characterized. In this study, we use hepatitis C virus (HCV) as a model system to study the balance of the two strategies. Combining viral infection cell culture assays with mathematical modeling, we characterize the dynamics of two different HCV strains (JFH-1, a clinical isolate, and Jc1-n, a laboratory strain), which have different viral release characteristics. We found that 0.63% and 1.70% of JFH-1 and Jc1-n intracellular viral RNAs, respectively, are used for producing and releasing progeny virions. Analysis of the Malthusian parameter of the HCV genome (i.e., initial proliferation rate) and the number of de novo infections (i.e., initial transmissibility) suggests that the leave-strategy provides a higher level of initial transmission for Jc1-n, whereas, in contrast, the stay-strategy provides a higher initial proliferation rate for JFH-1. Thus, theoretical-experimental analysis of viral dynamics enables us to better understand the proliferation strategies of viruses, which contributes to the efficient control of virus transmission. Ours is the first study to analyze the stay-leave trade-off during the viral life cycle and the significance of the replication-release switching mechanism for viral proliferation.

The machinery for endocytosis of epidermal growth factor receptor coordinates the transport of incoming hepatitis B virus to the endosomal network.

Sodium taurocholate cotransporting polypeptide (NTCP) is expressed at the surface of human hepatocytes and functions as an entry receptor of hepatitis B virus (HBV). Recently, we have reported that epidermal growth factor receptor (EGFR) is involved in NTCP-mediated viral internalization during the cell entry process. Here, we analyzed which function of EGFR is essential for mediating HBV internalization. In contrast to the reported crucial function of EGFR-downstream signaling for the entry of hepatitis C virus (HCV), blockade of EGFR-downstream signaling proteins, including mitogen-activated protein kinase (MAPK), phosphoinositide 3-kinase (PI3K), and signal transducer and activator of transcription (STAT), had no or only minor effects on HBV infection. Instead, deficiency of EGFR endocytosis resulting from either a deleterious mutation in EGFR or genetic knockdown of endocytosis adaptor molecules abrogated internalization of HBV via NTCP and prevented viral infection. EGFR activation triggered a time-dependent relocalization of HBV preS1 to the early and late endosomes and to lysosomes in concert with EGFR transport. Suppression of EGFR ubiquitination by site-directed mutagenesis or by knocking down two EGFR-sorting molecules, signal-transducing adaptor molecule (STAM) and lysosomal protein transmembrane 4ß (LAPTM4B), suggested that EGFR transport to the late endosome is critical for efficient HBV infection. Cumulatively, these results support the idea that the EGFR endocytosis/sorting machinery drives the translocation of NTCP-bound HBV from the cell surface to the endosomal network, which eventually enables productive viral infection.

Revealing uninfected and infected target cell dynamics from peripheral blood data in highly and less pathogenic simian/human immunodeficiency virus infected Rhesus macaque.

Since chimeric simian and human immunodeficiency viruses (SHIVs) used here, that is, SHIV-#64 and -KS661 utilize both CCR5 and CXCR4 chemokine receptors, they have broad target cell properties. A highly pathogenic SHIV strain, SHIV-KS661, causes an infection that systemically depletes the CD4+ T cells of Rhesus macaques, while a less pathogenic strain, SHIV-#64, does not cause severe symptoms in the macaques. In our previous studies, we established in vitro quantification system for virus infection dynamics, and concluded that SHIV-KS661 effectively produces infectious virions compared with SHIV-#64 in the HSC-F cell culture. However, in vivo dynamics of SHIV infection have not been well understood. To quantify SHIV-#64 and -KS661 infection dynamics in Rhesus macaques, we developed a novel approach and analyzed total CD4+ T cells and viral load in peripheral blood, and reproduced the expected dynamics for the uninfected and infected CD4+ T cells in silico. Using our previous cell culture experimental datasets, we revealed that an infection rate constant is different between SHIV-#64 and -KS661, but the viral production rate and the death rate are similar for the both strains. Thus, here, we assumed these relations in our in vivo data and carried out the data fitting. We performed Bayesian estimation for the whole dataset using MCMC sampling, and simultaneously fitted our novel model to total CD4+ T cells and viral load of SHIV-#64 and -KS661 infection. Our analyses explained that the Malthusian parameter (i.e., fitness of virus infection) and the basic reproduction number (i.e., potential of virus infection) for SHIV-KS661 are significantly higher than those of SHIV-#64. In addition, we demonstrated that the number of uninfected CD4+ T cells in SHIV-KS661 infected Rhesus macaques decreases to the significantly lower value than that before the inoculation several days earlier compared with SHIV-#64 infection. Taken together, the differences between SHIV-#64 and -KS661 infection before the peak viral load might determine the subsequent destiny, that is, whether the systemic CD4+ T cell depletion occurs or the host immune response develop.

Epidermal growth factor receptor is a host-entry cofactor triggering hepatitis B virus internalization.

Sodium taurocholate cotransporting polypeptide (NTCP) is a host cell receptor required for hepatitis B virus (HBV) entry. However, the susceptibility of NTCP-expressing cells to HBV is diverse depending on the culture condition. Stimulation with epidermal growth factor (EGF) was found to potentiate cell susceptibility to HBV infection. Here, we show that EGF receptor (EGFR) plays a critical role in HBV virion internalization. In EGFR-knockdown cells, HBV or its preS1-specific fluorescence peptide attached to the cell surface, but its internalization was attenuated. PreS1 internalization and HBV infection could be rescued by complementation with functional EGFR. Interestingly, the HBV/preS1-NTCP complex at the cell surface was internalized concomitant with the endocytotic relocalization of EGFR. Molecular interaction between NTCP and EGFR was documented by immunoprecipitation assay. Upon dissociation from functional EGFR, NTCP no longer functioned to support viral infection, as demonstrated by either (i) the introduction of NTCP point mutation that disrupted its interaction with EGFR, (ii) the detrimental effect of decoy peptide interrupting the NTCP-EGFR interaction, or (iii) the pharmacological inactivation of EGFR. Together, these data support the crucial role of EGFR in mediating HBV-NTCP internalization into susceptible cells. EGFR thus provides a yet unidentified missing link from the cell-surface HBV-NTCP attachment to the viral invasion beyond the host cell membrane.

Human-Specific Adaptations in Vpu Conferring Anti-tetherin Activity Are Critical for Efficient Early HIV-1 Replication In Vivo.

The HIV-1-encoded accessory protein Vpu exerts several immunomodulatory functions, including counteraction of the host restriction factor tetherin, downmodulation of CD4, and inhibition of NF-κB activity to facilitate HIV-1 infection. However, the relative contribution of individual Vpu functions to HIV-1 infection in vivo remained unclear. Here, we used a humanized mouse model and HIV-1 strains with selective mutations in vpu to demonstrate that the anti-tetherin activity of Vpu is a prerequisite for efficient viral spread during the early phase of infection. Mathematical modeling and gain-of-function mutations in SIVcpz, the simian precursor of pandemic HIV-1, corroborate this finding. Blockage of interferon signaling combined with transcriptome analyses revealed that basal tetherin levels are sufficient to control viral replication. These results establish tetherin as a key effector of the intrinsic immune defense against HIV-1, and they demonstrate that Vpu-mediated tetherin antagonism is critical for efficient viral spread during the initial phase of HIV-1 replication.

Sporadic on/off switching of HTLV-1 Tax expression is crucial to maintain the whole population of virus-induced leukemic cells.

Viruses causing chronic infection artfully manipulate infected cells to enable viral persistence in vivo under the pressure of immunity. Human T-cell leukemia virus type 1 (HTLV-1) establishes persistent infection mainly in CD4+ T cells in vivo and induces leukemia in this subset. HTLV-1-encoded Tax is a critical transactivator of viral replication and a potent oncoprotein, but its significance in pathogenesis remains obscure due to its very low level of expression in vivo. Here, we show that Tax is expressed in a minor fraction of leukemic cells at any given time, and importantly, its expression spontaneously switches between on and off states. Live cell imaging revealed that the average duration of one episode of Tax expression is ∼19 hours. Knockdown of Tax rapidly induced apoptosis in most cells, indicating that Tax is critical for maintaining the population, even if its short-term expression is limited to a small subpopulation. Single-cell analysis and computational simulation suggest that transient Tax expression triggers antiapoptotic machinery, and this effect continues even after Tax expression is diminished; this activation of the antiapoptotic machinery is the critical event for maintaining the population. In addition, Tax is induced by various cytotoxic stresses and also promotes HTLV-1 replication. Thus, it seems that Tax protects infected cells from apoptosis and increases the chance of viral transmission at a critical moment. Keeping the expression of Tax minimal but inducible on demand is, therefore, a fundamental strategy of HTLV-1 to promote persistent infection and leukemogenesis.

HIV-1 competition experiments in humanized mice show that APOBEC3H imposes selective pressure and promotes virus adaptation.

APOBEC3 (A3) family proteins are DNA cytosine deaminases recognized for contributing to HIV-1 restriction and mutation. Prior studies have demonstrated that A3D, A3F, and A3G enzymes elicit a robust anti-HIV-1 effect in cell cultures and in humanized mouse models. Human A3H is polymorphic and can be categorized into three phenotypes: stable, intermediate, and unstable. However, the anti-viral effect of endogenous A3H in vivo has yet to be examined. Here we utilize a hematopoietic stem cell-transplanted humanized mouse model and demonstrate that stable A3H robustly affects HIV-1 fitness in vivo. In contrast, the selection pressure mediated by intermediate A3H is relaxed. Intriguingly, viral genomic RNA sequencing reveled that HIV-1 frequently adapts to better counteract stable A3H during replication in humanized mice. Molecular phylogenetic analyses and mathematical modeling suggest that stable A3H may be a critical factor in human-to-human viral transmission. Taken together, this study provides evidence that stable variants of A3H impose selective pressure on HIV-1.

A heart-brain-kidney network controls adaptation to cardiac stress through tissue macrophage activation.

Heart failure is a complex clinical syndrome characterized by insufficient cardiac function. In addition to abnormalities intrinsic to the heart, dysfunction of other organs and dysregulation of systemic factors greatly affect the development and consequences of heart failure. Here we show that the heart and kidneys function cooperatively in generating an adaptive response to cardiac pressure overload. In mice subjected to pressure overload in the heart, sympathetic nerve activation led to activation of renal collecting-duct (CD) epithelial cells. Cell-cell interactions among activated CD cells, tissue macrophages and endothelial cells within the kidney led to secretion of the cytokine CSF2, which in turn stimulated cardiac-resident Ly6Clo macrophages, which are essential for the myocardial adaptive response to pressure overload. The renal response to cardiac pressure overload was disrupted by renal sympathetic denervation, adrenergic ß2-receptor blockade or CD-cell-specific deficiency of the transcription factor KLF5. Moreover, we identified amphiregulin as an essential cardioprotective mediator produced by cardiac Ly6Clo macrophages. Our results demonstrate a dynamic interplay between the heart, brain and kidneys that is necessary for adaptation to cardiac stress, and they highlight the homeostatic functions of tissue macrophages and the sympathetic nervous system.

Duration of SHIV production by infected cells is not exponentially distributed: Implications for estimates of infection parameters and antiviral efficacy.

The duration of the eclipse phase, from cell infection to the production and release of the first virion progeny, immediately followed by the virus-production phase, from the first to the last virion progeny, are important steps in a viral infection, by setting the pace of infection progression and modulating the response to antiviral therapy. Using a mathematical model (MM) and data for the infection of HSC-F cells with SHIV in vitro, we reconfirm our earlier finding that the eclipse phase duration follows a fat-tailed distribution, lasting 19 h (18-20 h). Most importantly, for the first time, we show that the virus-producing phase duration, which lasts 11 h (9.8-12 h), follows a normal-like distribution, and not an exponential distribution as is typically assumed. We explore the significance of this finding and its impact on analysis of plasma viral load decays in HIV patients under antiviral therapy. We find that incorrect assumptions about the eclipse and virus-producing phase distributions can lead to an overestimation of antiviral efficacy. Additionally, our predictions for the rate of plasma HIV decay under integrase inhibitor therapy offer an opportunity to confirm whether HIV production duration in vivo also follows a normal distribution, as demonstrated here for SHIV infections in vitro.

Dynamics of HIV infection in lymphoid tissue network.

Human immunodeficiency virus (HIV) is a fast replicating ribonucleic acid virus, which can easily mutate in order to escape the effects of drug administration. Hence, understanding the basic mechanisms underlying HIV persistence in the body is essential in the development of new therapies that could eradicate HIV infection. Lymphoid tissues are the primary sites of HIV infection. Despite the recent progress in real-time monitoring technology, HIV infection dynamics in a whole body is unknown. Mathematical modeling and simulations provide speculations on global behavior of HIV infection in the lymphatic system. We propose a new mathematical model that describes the spread of HIV infection throughout the lymphoid tissue network. In order to represent the volume difference between lymphoid tissues, we propose the proportionality of several kinetic parameters to the lymphoid tissues' volume distribution. Under this assumption, we perform extensive numerical computations in order to simulate the spread of HIV infection in the lymphoid tissue network. Numerical computations simulate single drug treatments of an HIV infection. One of the important biological speculations derived from this study is a drug saturation effect generated by lymphoid network connection. This implies that a portion of reservoir lymphoid tissues to which drug is not sufficiently delivered would inhibit HIV eradication despite of extensive drug injection.

HIV-1 Vpr accelerates viral replication during acute infection by exploitation of proliferating CD4+ T cells in vivo.

The precise role of viral protein R (Vpr), an HIV-1-encoded protein, during HIV-1 infection and its contribution to the development of AIDS remain unclear. Previous reports have shown that Vpr has the ability to cause G2 cell cycle arrest and apoptosis in HIV-1-infected cells in vitro. In addition, vpr is highly conserved in transmitted/founder HIV-1s and in all primate lentiviruses, which are evolutionarily related to HIV-1. Although these findings suggest an important role of Vpr in HIV-1 pathogenesis, its direct evidence in vivo has not been shown. Here, by using a human hematopoietic stem cell-transplanted humanized mouse model, we demonstrated that Vpr causes G2 cell cycle arrest and apoptosis predominantly in proliferating CCR5(+) CD4(+) T cells, which mainly consist of regulatory CD4(+) T cells (Tregs), resulting in Treg depletion and enhanced virus production during acute infection. The Vpr-dependent enhancement of virus replication and Treg depletion is observed in CCR5-tropic but not CXCR4-tropic HIV-1-infected mice, suggesting that these effects are dependent on the coreceptor usage by HIV-1. Immune activation was observed in CCR5-tropic wild-type but not in vpr-deficient HIV-1-infected humanized mice. When humanized mice were treated with denileukin diftitox (DD), to deplete Tregs, DD-treated humanized mice showed massive activation/proliferation of memory T cells compared to the untreated group. This activation/proliferation enhanced CCR5 expression in memory CD4(+) T cells and rendered them more susceptible to CCR5-tropic wild-type HIV-1 infection than to vpr-deficient virus. Taken together, these results suggest that Vpr takes advantage of proliferating CCR5(+) CD4(+) T cells for enhancing viremia of CCR5-tropic HIV-1. Because Tregs exist in a higher cycling state than other T cell subsets, Tregs appear to be more vulnerable to exploitation by Vpr during acute HIV-1 infection.

A race between tumor immunoescape and genome maintenance selects for optimum levels of (epi)genetic instability.

The human immune system functions to provide continuous body-wide surveillance to detect and eliminate foreign agents such as bacteria and viruses as well as the body's own cells that undergo malignant transformation. To counteract this surveillance, tumor cells evolve mechanisms to evade elimination by the immune system; this tumor immunoescape leads to continuous tumor expansion, albeit potentially with a different composition of the tumor cell population ("immunoediting"). Tumor immunoescape and immunoediting are products of an evolutionary process and are hence driven by mutation and selection. Higher mutation rates allow cells to more rapidly acquire new phenotypes that help evade the immune system, but also harbor the risk of an inability to maintain essential genome structure and functions, thereby leading to an error catastrophe. In this paper, we designed a novel mathematical framework, based upon the quasispecies model, to study the effects of tumor immunoediting and the evolution of (epi)genetic instability on the abundance of tumor and immune system cells. We found that there exists an optimum number of tumor variants and an optimum magnitude of mutation rates that maximize tumor progression despite an active immune response. Our findings provide insights into the dynamics of tumorigenesis during immune system attacks and help guide the choice of treatment strategies that best inhibit diverse tumor cell populations.

Vpu augments the initial burst phase of HIV-1 propagation and downregulates BST2 and CD4 in humanized mice.

While human cells express potent antiviral proteins as part of the host defense repertoire, viruses have evolved their own arsenal of proteins to antagonize them. BST2 was identified as an inhibitory cellular protein of HIV-1 replication, which tethers virions to the cell surface to prevent their release. On the other hand, the HIV-1 accessory protein, Vpu, has the ability to downregulate and counteract BST2. Vpu also possesses the ability to downmodulate cellular CD4 and SLAMF6 molecules expressed on infected cells. However, the role of Vpu in HIV-1 infection in vivo remains unclear. Here, using a human hematopoietic stem cell-transplanted humanized mouse model, we demonstrate that Vpu contributes to the efficient spread of HIV-1 in vivo during the acute phase of infection. Although Vpu did not affect viral cytopathicity, target cell preference, and the level of viral protein expression, the amount of cell-free virions in vpu-deficient HIV-1-infected mice was profoundly lower than that in wild-type HIV-1-infected mice. We provide a novel insight suggesting that Vpu concomitantly downregulates BST2 and CD4, but not SLAMF6, from the surface of infected cells. Furthermore, we show evidence suggesting that BST2 and CD4 impair the production of cell-free infectious virions but do not associate with the efficiency of cell-to-cell HIV-1 transmission. Taken together, our findings suggest that Vpu downmodulates BST2 and CD4 in infected cells and augments the initial burst of HIV-1 replication in vivo. This is the first report demonstrating the role of Vpu in HIV-1 infection in an in vivo model.