Stochastic pausing at latent HIV-1 promoters generates transcriptional bursting

The cell biology of HIV-1 latency and rebound

Transcriptionally latent forms of replication-competent proviruses, present primarily in a small subset of memory CD4+ T cells, pose the primary barrier to a cure for HIV-1 infection because they are the source of the viral rebound that almost inevitably follows the interruption of antiretroviral therapy. Over the last 30 years, many of the factors essential for initiating HIV-1 transcription have been identified in studies performed using transformed cell lines, such as the Jurkat T-cell model. However, as highlighted in this review, several poorly understood mechanisms still need to be elucidated, including the molecular basis for promoter-proximal pausing of the transcribing complex and the detailed mechanism of the delivery of P-TEFb from 7SK snRNP. Furthermore, the central paradox of HIV-1 transcription remains unsolved: how are the initial rounds of transcription achieved in the absence of Tat? A critical limitation of the transformed cell models is that they do not recapitulate the transitions between active effector cells and quiescent memory T cells. Therefore, investigation of the molecular mechanisms of HIV-1 latency reversal and LRA efficacy in a proper physiological context requires the utilization of primary cell models. Recent mechanistic studies of HIV-1 transcription using latently infected cells recovered from donors and ex vivo cellular models of viral latency have demonstrated that the primary blocks to HIV-1 transcription in memory CD4+ T cells are restrictive epigenetic features at the proviral promoter, the cytoplasmic sequestration of key transcription initiation factors such as NFAT and NF-κB, and the vanishingly low expression of the cellular transcription elongation factor P-TEFb. One of the foremost schemes to eliminate the residual reservoir is to deliberately reactivate latent HIV-1 proviruses to enable clearance of persisting latently infected cells-the "Shock and Kill" strategy. For "Shock and Kill" to become efficient, effective, non-toxic latency-reversing agents (LRAs) must be discovered. Since multiple restrictions limit viral reactivation in primary cells, understanding the T-cell signaling mechanisms that are essential for stimulating P-TEFb biogenesis, initiation factor activation, and reversing the proviral epigenetic restrictions have become a prerequisite for the development of more effective LRAs.

HIV-Tocky system to visualize proviral expression dynamics

Determinants of HIV-1 latency establishment are yet to be elucidated. HIV reservoir comprises a rare fraction of infected cells that can survive host and virus-mediated killing. In vitro reporter models so far offered a feasible means to inspect this population, but with limited capabilities to dissect provirus silencing dynamics. Here, we describe a new HIV reporter model, HIV-Timer of cell kinetics and activity (HIV-Tocky) with dual fluorescence spontaneous shifting to reveal provirus silencing and reactivation dynamics. This unique feature allows, for the first time, identifying two latent populations: a directly latent, and a recently silenced subset, with the latter having integration features suggestive of stable latency. Our proposed model can help address the heterogeneous nature of HIV reservoirs and offers new possibilities for evaluating eradication strategies.

Release of P-TEFb from the Super Elongation Complex promotes HIV-1 latency reversal

The persistence of HIV-1 in long-lived latent reservoirs during suppressive antiretroviral therapy (ART) remains one of the principal barriers to a functional cure. Blocks to transcriptional elongation play a central role in maintaining the latent state, and several latency reversal strategies focus on the release of positive transcription elongation factor b (P-TEFb) from sequestration by negative regulatory complexes, such as the 7SK complex and BRD4. Another major cellular reservoir of P-TEFb is in Super Elongation Complexes (SECs), which play broad regulatory roles in host gene expression. Still, it is unknown if the release of P-TEFb from SECs is a viable latency reversal strategy. Here, we demonstrate that the SEC is not required for HIV-1 replication in primary CD4+ T cells and that a small molecular inhibitor of the P-TEFb/SEC interaction (termed KL-2) increases viral transcription. KL-2 acts synergistically with other latency reversing agents (LRAs) to reactivate viral transcription in several cell line models of latency in a manner that is, at least in part, dependent on the viral Tat protein. Finally, we demonstrate that KL-2 enhances viral reactivation in peripheral blood mononuclear cells (PBMCs) from people living with HIV on suppressive ART, most notably in combination with inhibitor of apoptosis protein antagonists (IAPi). Taken together, these results suggest that the release of P-TEFb from cellular SECs may be a novel route for HIV-1 latency reactivation.

HIV-1 Proviral Genome Engineering with CRISPR-Cas9 for Mechanistic Studies

HIV-1 latency remains a barrier to a functional cure because of the ability of virtually silent yet inducible proviruses within reservoir cells to transcriptionally reactivate upon cell stimulation. HIV-1 reactivation occurs through the sequential action of host transcription factors (TFs) during the "host phase" and the viral TF Tat during the "viral phase", which together facilitate the positive feedback loop required for exponential transcription, replication, and pathogenesis. The sequential action of these TFs poses a challenge to precisely delineate the contributions of the host and viral phases of the transcriptional program to guide future mechanistic and therapeutic studies. To address this limitation, we devised a genome engineering approach to mutate tat and create a genetically matched pair of Jurkat T cell clones harboring HIV-1 at the same integration site with and without Tat expression. By comparing the transcriptional profile of both clones, the transition point between the host and viral phases was defined, providing a system that enables the temporal mechanistic interrogation of HIV-1 transcription prior to and after Tat synthesis. Importantly, this CRISPR method is broadly applicable to knockout individual viral proteins or genomic regulatory elements to delineate their contributions to various aspects of the viral life cycle and ultimately may facilitate therapeutic approaches in our race towards achieving a functional cure.

Time will tell: comparing timescales to gain insight into transcriptional bursting

Recent imaging studies have captured the dynamics of regulatory events of transcription inside living cells. These events include transcription factor (TF) DNA binding, chromatin remodeling and modification, enhancer-promoter (E-P) proximity, cluster formation, and preinitiation complex (PIC) assembly. Together, these molecular events culminate in stochastic bursts of RNA synthesis, but their kinetic relationship remains largely unclear. In this review, we compare the timescales of upstream regulatory steps (input) with the kinetics of transcriptional bursting (output) to generate mechanistic models of transcription dynamics in single cells. We highlight open questions and potential technical advances to guide future endeavors toward a quantitative and kinetic understanding of transcription regulation.

Tissue-specific RNA Polymerase II promoter-proximal pause release and burst kinetics in a Drosophila embryonic patterning network

BACKGROUND

Formation of tissue-specific transcriptional programs underlies multicellular development, including dorsoventral (DV) patterning of the Drosophila embryo. This involves interactions between transcriptional enhancers and promoters in a chromatin context, but how the chromatin landscape influences transcription is not fully understood.

RESULTS

Here we comprehensively resolve differential transcriptional and chromatin states during Drosophila DV patterning. We find that RNA Polymerase II pausing is established at DV promoters prior to zygotic genome activation (ZGA), that pausing persists irrespective of cell fate, but that release into productive elongation is tightly regulated and accompanied by tissue-specific P-TEFb recruitment. DV enhancers acquire distinct tissue-specific chromatin states through CBP-mediated histone acetylation that predict the transcriptional output of target genes, whereas promoter states are more tissue-invariant. Transcriptome-wide inference of burst kinetics in different cell types revealed that while DV genes are generally characterized by a high burst size, either burst size or frequency can differ between tissues.

CONCLUSIONS

The data suggest that pausing is established by pioneer transcription factors prior to ZGA and that release from pausing is imparted by enhancer chromatin state to regulate bursting in a tissue-specific manner in the early embryo. Our results uncover how developmental patterning is orchestrated by tissue-specific bursts of transcription from Pol II primed promoters in response to enhancer regulatory cues.

HIV-1 Gag co-localizes with euchromatin histone marks at the nuclear periphery

IMPORTANCE

The traditional view of retrovirus assembly posits that packaging of gRNA by HIV-1 Gag occurs in the cytoplasm or at the plasma membrane. However, our previous studies showing that HIV-1 Gag enters the nucleus and binds to USvRNA at transcription sites suggest that gRNA selection may occur in the nucleus. In the present study, we observed that HIV-1 Gag trafficked to the nucleus and co-localized with USvRNA within 8 hours of expression. In infected T cells (J-Lat 10.6) reactivated from latency and in a HeLa cell line stably expressing an inducible Rev-dependent HIV-1 construct, we found that Gag preferentially localized with euchromatin histone marks associated with enhancer and promoter regions near the nuclear periphery, which is the favored site HIV-1 integration. These observations support the innovative hypothesis that HIV-1 Gag associates with euchromatin-associated histones to localize to active transcription sites, promoting capture of newly synthesized gRNA for packaging.

4D nucleome equation predicts gene expression controlled by long-range enhancer-promoter interaction

Recent experimental evidence strongly supports that three-dimensional (3D) long-range enhancer-promoter (E-P) interactions have important influences on gene-expression dynamics, but it is unclear how the interaction information is translated into gene expression over time (4D). To address this question, we developed a general theoretical framework (named as a 4D nucleome equation), which integrates E-P interactions on chromatin and biochemical reactions of gene transcription. With this equation, we first present the distribution of mRNA counts as a function of the E-P genomic distance and then reveal a power-law scaling of the expression level in this distance. Interestingly, we find that long-range E-P interactions can induce bimodal and trimodal mRNA distributions. The 4D nucleome equation also allows for model selection and parameter inference. When this equation is applied to the mouse embryonic stem cell smRNA-FISH data and the E-P genomic-distance data, the predicted E-P contact probability and mRNA distribution are in good agreement with experimental results. Further statistical inference indicates that the E-P interactions prefer to modulate the mRNA level by controlling promoter activation and transcription initiation rates. Our model and results provide quantitative insights into both spatiotemporal gene-expression determinants (i.e., long-range E-P interactions) and cellular fates during development.

Mathematical Models of HIV-1 Dynamics, Transcription, and Latency

HIV-1 latency is a major barrier to curing infections with antiretroviral therapy and, consequently, to eliminating the disease globally. The establishment, maintenance, and potential clearance of latent infection are complex dynamic processes and can be best described with the help of mathematical models followed by experimental validation. Here, we review the use of viral dynamics models for HIV-1, with a focus on applications to the latent reservoir. Such models have been used to explain the multi-phasic decay of viral load during antiretroviral therapy, the early seeding of the latent reservoir during acute infection and the limited inflow during treatment, the dynamics of viral blips, and the phenomenon of post-treatment control. Finally, we discuss that mathematical models have been used to predict the efficacy of potential HIV-1 cure strategies, such as latency-reversing agents, early treatment initiation, or gene therapies, and to provide guidance for designing trials of these novel interventions.

Gene Activity as the Predictive Indicator for Transcriptional Bursting Dynamics

Transcription commonly occurs in bursts, with alternating productive (ON) and quiescent (OFF) periods, governing mRNA production rates. Yet, how transcription is regulated through bursting dynamics remains unresolved. In this study, we conduct real-time measurements of endogenous transcriptional bursting with single-mRNA sensitivity. Leveraging the diverse transcriptional activities in early fly embryos, we uncover stringent relationships between bursting parameters. Specifically, we find that the durations of ON and OFF periods are linked. Regardless of the developmental stage or body-axis position, gene activity levels predict the average ON and OFF periods of individual alleles. Lowly transcribing alleles predominantly modulate OFF durations (burst frequency), while highly transcribing alleles primarily tune ON durations (burst size). Importantly, these relationships persist even under perturbation of cis-regulatory elements or trans-factors. This suggests a novel mechanistic constraint governing bursting dynamics rather than a modular control of distinct parameters by distinct regulatory processes. Our study provides a foundation for future investigations into the molecular mechanisms underpinning spatiotemporal transcriptional control.

Transcriptional Stochasticity as a Key Aspect of HIV-1 Latency

This review summarizes current advances in the role of transcriptional stochasticity in HIV-1 latency, which were possible in a large part due to the development of single-cell approaches. HIV-1 transcription proceeds in bursts of RNA production, which stem from the stochastic switching of the viral promoter between ON and OFF states. This switching is caused by random binding dynamics of transcription factors and nucleosomes to the viral promoter and occurs at several time scales from minutes to hours. Transcriptional bursts are mainly controlled by the core transcription factors TBP, SP1 and NF-κb, the chromatin status of the viral promoter and RNA polymerase II pausing. In particular, spontaneous variability in the promoter chromatin creates heterogeneity in the response to activators such as TNF-α, which is then amplified by the Tat feedback loop to generate high and low viral transcriptional states. This phenomenon is likely at the basis of the partial and stochastic response of latent T cells from HIV-1 patients to latency-reversing agents, which is a barrier for the development of shock-and-kill strategies of viral eradication. A detailed understanding of the transcriptional stochasticity of HIV-1 and the possibility to precisely model this phenomenon will be important assets to develop more effective therapeutic strategies.

BurstDECONV: a signal deconvolution method to uncover mechanisms of transcriptional bursting in live cells

Monitoring transcription in living cells gives access to the dynamics of this complex fundamental process. It reveals that transcription is discontinuous, whereby active periods (bursts) are separated by one or several types of inactive periods of distinct lifetimes. However, decoding temporal fluctuations arising from live imaging and inferring the distinct transcriptional steps eliciting them is a challenge. We present BurstDECONV, a novel statistical inference method that deconvolves signal traces into individual transcription initiation events. We use the distribution of waiting times between successive polymerase initiation events to identify mechanistic features of transcription such as the number of rate-limiting steps and their kinetics. Comparison of our method to alternative methods emphasizes its advantages in terms of precision and flexibility. Unique features such as the direct determination of the number of promoter states and the simultaneous analysis of several potential transcription models make BurstDECONV an ideal analytic framework for live cell transcription imaging experiments. Using simulated realistic data, we found that our method is robust with regards to noise or suboptimal experimental designs. To show its generality, we applied it to different biological contexts such as Drosophila embryos or human cells.

Bayesian analysis dissects kinetic modulation during non-stationary gene expression

Labelling of nascent stem loops with fluorescent proteins has fostered the visualization of transcription in living cells. Quantitative analysis of recorded fluorescence traces can shed light on kinetic transcription parameters and regulatory mechanisms. However, existing methods typically focus on steady state dynamics. Here, we combine a stochastic process transcription model with a hierarchical Bayesian method to infer global as well locally shared parameters for groups of cells and recover unobserved quantities such as initiation times and polymerase loading of the gene. We apply our approach to the cyclic response of the yeast CUP1 locus to heavy metal stress. Within the previously described slow cycle of transcriptional activity on the scale of minutes, we discover fast time-modulated bursting on the scale of seconds. Model comparison suggests that slow oscillations of transcriptional output are regulated by the amplitude of the bursts. Several polymerases may initiate during a burst.

Modulation of transcription burst amplitude underpins dosage compensation in the Drosophila embryo

Dosage compensation, the balancing of X-linked gene expression between sexes and to the autosomes, is critical to an organism's fitness and survival. In Drosophila, dosage compensation involves hypertranscription of the male X chromosome. Here, we use quantitative live imaging and modeling at single-cell resolution to study X chromosome dosage compensation in Drosophila. We show that the four X chromosome genes studied undergo transcriptional bursting in male and female embryos. Mechanistically, our data reveal that transcriptional upregulation of male X chromosome genes is primarily mediated by a higher RNA polymerase II initiation rate and burst amplitude across the expression domain. In contrast, burst frequency is spatially modulated in nuclei within the expression domain in response to different transcription factor concentrations to tune the transcriptional response. Together, these data show how the local and global regulation of distinct burst parameters can establish the complex transcriptional outputs underpinning developmental patterning.

Inferring transcriptional bursting kinetics from single-cell snapshot data using a generalized telegraph model

Gene expression has inherent stochasticity resulting from transcription's burst manners. Single-cell snapshot data can be exploited to rigorously infer transcriptional burst kinetics, using mathematical models as blueprints. The classical telegraph model (CTM) has been widely used to explain transcriptional bursting with Markovian assumptions. However, growing evidence suggests that the gene-state dwell times are generally non-exponential, as gene-state switching is a multi-step process in organisms. Therefore, interpretable non-Markovian mathematical models and efficient statistical inference methods are urgently required in investigating transcriptional burst kinetics. We develop an interpretable and tractable model, the generalized telegraph model (GTM), to characterize transcriptional bursting that allows arbitrary dwell-time distributions, rather than exponential distributions, to be incorporated into the ON and OFF switching process. Based on the GTM, we propose an inference method for transcriptional bursting kinetics using an approximate Bayesian computation framework. This method demonstrates an efficient and scalable estimation of burst frequency and burst size on synthetic data. Further, the application of inference to genome-wide data from mouse embryonic fibroblasts reveals that GTM would estimate lower burst frequency and higher burst size than those estimated by CTM. In conclusion, the GTM and the corresponding inference method are effective tools to infer dynamic transcriptional bursting from static single-cell snapshot data.

Enhancing HIV-1 latency reversal through regulating the elongating RNA Pol II pause-release by a small-molecule disruptor of PAF1C

The polymerase-associated factor 1 complex (PAF1C) is a key, post-initiation transcriptional regulator of both promoter-proximal pausing and productive elongation catalyzed by RNA Pol II and is also involved in transcriptional repression of viral gene expression during human immunodeficiency virus-1 (HIV-1) latency. Using a molecular docking-based compound screen in silico and global sequencing-based candidate evaluation in vivo, we identified a first-in-class, small-molecule inhibitor of PAF1C (iPAF1C) that disrupts PAF1 chromatin occupancy and induces global release of promoter-proximal paused RNA Pol II into gene bodies. Transcriptomic analysis revealed that iPAF1C treatment mimics acute PAF1 subunit depletion and impairs RNA Pol II pausing at heat shock-down-regulated genes. Furthermore, iPAF1C enhances the activity of diverse HIV-1 latency reversal agents both in cell line latency models and in primary cells from persons living with HIV-1. In sum, this study demonstrates that efficient disruption of PAF1C by a first-in-class, small-molecule inhibitor may have therapeutic potential for improving current HIV-1 latency reversal strategies.

HIV-1 Gag colocalizes with euchromatin histone marks at the nuclear periphery

The retroviral Gag protein of human immunodeficiency virus type 1 (HIV-1) plays a central role in the selection of unspliced viral genomic RNA for packaging into new virions. Previously, we demonstrated that full-length HIV-1 Gag undergoes nuclear trafficking where it associates with unspliced viral RNA (vRNA) at transcription sites. To further explore the kinetics of HIV-1 Gag nuclear localization, we used biochemical and imaging techniques to examine the timing of HIV-1 entry into the nucleus. We also aimed to determine more precisely Gag's subnuclear distribution to test the hypothesis that Gag would be associated with euchromatin, the transcriptionally active region of the nucleus. We observed that HIV-1 Gag localized to the nucleus shortly after its synthesis in the cytoplasm, suggesting that nuclear trafficking was not strictly concentration-dependent. Furthermore, we found that HIV-1 Gag preferentially localized to the transcriptionally active euchromatin fraction compared to the heterochromatin-rich region in a latently-infected CD4+ T cell line (J-Lat 10.6) treated with latency-reversal agents. Interestingly, HIV-1 Gag was more closely associated with transcriptionally-active histone markers near the nuclear periphery, where the HIV-1 provirus was previously shown to integrate. Although the precise function of Gag's association with histones in transcriptionally-active chromatin remains uncertain, together with previous reports, this finding is consistent with a potential role for euchromatin-associated Gag molecules to select newly transcribed unspliced vRNA during the initial stage of virion assembly.

Importance

The traditional view of retroviral assembly posits that HIV-1 Gag selection of unspliced vRNA begins in the cytoplasm. However, our previous studies demonstrated that HIV-1 Gag enters the nucleus and binds to unspliced HIV-1 RNA at transcription sites, suggesting that genomic RNA selection may occur in the nucleus. In the present study, we observed nuclear entry of HIV-1 Gag and co-localization with unspliced viral RNA within 8 hours post-expression. In CD4+ T cells (J-Lat 10.6) treated with latency reversal agents, as well as a HeLa cell line stably expressing an inducible Rev-dependent provirus, we found that HIV-1 Gag preferentially localized with histone marks associated with enhancer and promoter regions of transcriptionally active euchromatin near the nuclear periphery, which favors HIV-1 proviral integration sites. These observations support the hypothesis that HIV-1 Gag hijacks euchromatin-associated histones to localize to active transcription sites, promoting capture of newly synthesized genomic RNA for packaging.

Genome-wide inference reveals that feedback regulations constrain promoter-dependent transcriptional burst kinetics

Gene expression in mammalian cells is highly variable and episodic, resulting in a series of discontinuous bursts of mRNAs. A challenge is to understand how static promoter architecture and dynamic feedback regulations dictate bursting on a genome-wide scale. Although single-cell RNA sequencing (scRNA-seq) provides an opportunity to address this challenge, effective analytical methods are scarce. We developed an interpretable and scalable inference framework, which combined experimental data with a mechanistic model to infer transcriptional burst kinetics (sizes and frequencies) and feedback regulations. Applying this framework to scRNA-seq data generated from embryonic mouse fibroblast cells, we found Simpson's paradoxes, i.e. genome-wide burst kinetics exhibit different characteristics in two cases without and with distinguishing feedback regulations. We also showed that feedbacks differently modulate burst frequencies and sizes and conceal the effects of transcription start site distributions on burst kinetics. Notably, only in the presence of positive feedback, TATA genes are expressed with high burst frequencies and enhancer-promoter interactions mainly modulate burst frequencies. The developed inference method provided a flexible and efficient way to investigate transcriptional burst kinetics and the obtained results would be helpful for understanding cell development and fate decision.

Shock and kill within the CNS: A promising HIV eradication approach?

The most studied HIV eradication approach is the "shock and kill" strategy, which aims to reactivate the latent reservoir by latency reversing agents (LRAs) and allowing elimination of these cells by immune-mediated clearance or viral cytopathic effects. The CNS is an anatomic compartment in which (persistent) HIV plays an important role in HIV-associated neurocognitive disorder. Restriction of the CNS by the blood-brain barrier is important for maintenance of homeostasis of the CNS microenvironment, which includes CNS-specific cell types, expression of transcription factors, and altered immune surveillance. Within the CNS predominantly myeloid cells such as microglia and perivascular macrophages are thought to be a reservoir of persistent HIV infection. Nevertheless, infection of T cells and astrocytes might also impact HIV infection in the CNS. Genetic adaptation to this microenvironment results in genetically distinct, compartmentalized viral populations with differences in transcription profiles. Because of these differences in transcription profiles, LRAs might have different effects within the CNS as compared with the periphery. Moreover, reactivation of HIV in the brain and elimination of cells within the CNS might be complex and could have detrimental consequences. Finally, independent of activity on latent HIV, LRAs themselves can have adverse neurologic effects. We provide an extensive overview of the current knowledge on compartmentalized (persistent) HIV infection in the CNS and on the "shock and kill" strategy. Subsequently, we reflect on the impact and promise of the "shock and kill" strategy on the elimination of persistent HIV in the CNS.

Viro-fluidics: Real-time analysis of virus production kinetics at the single-cell level

Real-time visualization and quantification of viruses released by a cell are crucial to further decipher infection processes. Kinetics studies at the single-cell level will circumvent the limitations of bulk assays with asynchronous virus replication. We have implemented a "viro-fluidic" method, which combines microfluidics and virology at single-cell and single-virus resolutions. As an experimental model, we used standard cell lines producing fluorescent HIV-like particles (VLPs). First, to scale the strategy to the single-cell level, we validated a sensitive flow virometry system to detect VLPs in low concentration samples (≥10⁴ VLPs/mL). Then, this system was coupled to a single-cell trapping device to monitor in real-time the VLPs released, one at a time, from single cells under cell culture conditions. Our results revealed an average production rate of 50 VLPs/h/cell similar to the rate estimated for the same cells grown in population. Thus, the virus-producing capacities of the trapped cells were preserved and its real-time monitoring was accurate. Moreover, single-cell analysis revealed a release of VLPs with stochastic bursts with typical time intervals of few minutes, revealing the existence of limiting step(s) in the virus biogenesis process. Our tools can be applied to other pathogens or to extracellular vesicles to elucidate the dissemination mechanisms of these biological nanoparticles.

A transcriptional cycling model recapitulates chromatin-dependent features of noisy inducible transcription

Activation of gene expression in response to environmental cues results in substantial phenotypic heterogeneity between cells that can impact a wide range of outcomes including differentiation, viral activation, and drug resistance. An important source of gene expression noise is transcriptional bursting, or the process by which transcripts are produced during infrequent bursts of promoter activity. Chromatin accessibility impacts transcriptional bursting by regulating the assembly of transcription factor and polymerase complexes on promoters, suggesting that the effect of an activating signal on transcriptional noise will depend on the initial chromatin state at the promoter. To explore this possibility, we simulated transcriptional activation using a transcriptional cycling model with three promoter states that represent chromatin remodeling, polymerase binding and pause release. We initiated this model over a large parameter range representing target genes with different chromatin environments, and found that, upon increasing the polymerase pause release rate to activate transcription, changes in gene expression noise varied significantly across initial promoter states. This model captured phenotypic differences in activation of latent HIV viruses integrated at different chromatin locations and mediated by the transcription factor NF-κB. Activating transcription in the model via increasing one or more of the transcript production rates, as occurs following NF-κB activation, reproduced experimentally measured transcript distributions for four different latent HIV viruses, as well as the bimodal pattern of HIV protein expression that leads to a subset of reactivated virus. Importantly, the parameter ‘activation path’ differentially affected gene expression noise, and ultimately viral activation, in line with experimental observations. This work demonstrates how upstream signaling pathways can be connected to biological processes that underlie transcriptional bursting, resulting in target gene-specific noise profiles following stimulation of a single upstream pathway.

Many genes are transcribed in infrequent bursts of mRNA production through a process called transcriptional bursting, which contributes to variability in responses between cells. Heterogeneity in cell responses can have important biological impacts, such as whether a cell supports viral replication or responds to a drug, and thus there is an effort to describe this process with mathematical models to predict biological outcomes. Previous models described bursting as a transition between an “OFF” state or an “ON” state, an elegant and simple mathematical representation of complex molecular mechanisms, but one which failed to capture how upstream activation signals affected bursting. To address this, we added an additional promoter state to better reflect biological mechanisms underlying bursting. By fitting this model to variable activation of quiescent HIV infections in T cells, we showed that our model more accurately described viral expression variability across cells in response to an upstream stimulus. Our work highlights how mathematical models can be further developed to understand complex biological mechanisms and suggests ways to connect transcriptional bursting to upstream activation pathways.

Defective HIV-1 genomes and their potential impact on HIV pathogenesis

Defective HIV-1 proviruses represent a population of viral genomes that are selected for by immune pressures, and clonally expanded to dominate the persistent HIV-1 proviral genome landscape. There are examples of RNA and protein expression from these compromised genomes which are generated by a variety of mechanisms. Despite the evidence that these proviruses are transcribed and translated, their role in HIV pathogenesis has not been fully explored. The potential for these genomes to participate in immune stimulation is particularly relevant considering the accumulation of cells harboring these defective proviruses over the course of antiretroviral therapy in people living with HIV. The expression of defective proviruses in different cells and tissues could drive innate sensing mechanisms and inflammation. They may also alter antiviral T cell responses and myeloid cell functions that directly contribute to HIV-1 associated chronic comorbidities. Understanding the impact of these defective proviruses needs to be considered as we advance cure strategies that focus on targeting the diverse population of HIV-1 proviral genomes.

Infectious RNA: Human Immunodeficiency Virus (HIV) Biology, Therapeutic Intervention, and the Quest for a Vaccine

Different mechanisms mediate the toxicity of RNA. Genomic retroviral mRNA hijacks infected host cell factors to enable virus replication. The viral genomic RNA of the human immunodeficiency virus (HIV) encompasses nine genes encoding in less than 10 kb all proteins needed for replication in susceptible host cells. To do so, the genomic RNA undergoes complex alternative splicing to facilitate the synthesis of the structural, accessory, and regulatory proteins. However, HIV strongly relies on the host cell machinery recruiting cellular factors to complete its replication cycle. Antiretroviral therapy (ART) targets different steps in the cycle, preventing disease progression to the acquired immunodeficiency syndrome (AIDS). The comprehension of the host immune system interaction with the virus has fostered the development of a variety of vaccine platforms. Despite encouraging provisional results in vaccine trials, no effective vaccine has been developed, yet. However, novel promising vaccine platforms are currently under investigation.

Silent transcription intervals and translational bursting lead to diverse phenotypic switching

For complex process of gene expression, we use theoretical analysis and stochastic simulations to study the phenotypic diversity induced by silent transcription intervals and translational bursting.

A Two-Color Haploid Genetic Screen Identifies Novel Host Factors Involved in HIV-1 Latency

To identify novel host factors as putative targets to reverse HIV-1 latency, we performed an insertional mutagenesis genetic screen in a latent HIV-1 infected pseudohaploid KBM7 cell line (Hap-Lat). Following mutagenesis, insertions were mapped to the genome, and bioinformatic analysis resulted in the identification of 69 candidate host genes involved in maintaining HIV-1 latency. A select set of candidate genes was functionally validated using short hairpin RNA (shRNA)-mediated depletion in latent HIV-1 infected J-Lat A2 and 11.1 T cell lines. We confirmed ADK, CHD9, CMSS1, EVI2B, EXOSC8, FAM19A, GRIK5, IRF2BP2, NF1, and USP15 as novel host factors involved in the maintenance of HIV-1 latency. Chromatin immunoprecipitation assays indicated that CHD9, a chromodomain helicase DNA-binding protein, maintains HIV-1 latency via direct association with the HIV-1 5′ long terminal repeat (LTR), and its depletion results in increased histone acetylation at the HIV-1 promoter, concomitant with HIV-1 latency reversal. FDA-approved inhibitors 5-iodotubercidin, trametinib, and topiramate, targeting ADK, NF1, and GRIK5, respectively, were characterized for their latency reversal potential. While 5-iodotubercidin exhibited significant cytotoxicity in both J-Lat and primary CD4⁺ T cells, trametinib reversed latency in J-Lat cells but not in latent HIV-1 infected primary CD4⁺ T cells. Importantly, topiramate reversed latency in cell line models, in latently infected primary CD4⁺ T cells, and crucially in CD4⁺ T cells from three people living with HIV-1 (PLWH) under suppressive antiretroviral therapy, without inducing T cell activation or significant toxicity. Thus, using an adaptation of a haploid forward genetic screen, we identified novel and druggable host factors contributing to HIV-1 latency.

Fighting HIV-1 Persistence: At the Crossroads of "Shoc-K and B-Lock"

Despite the success of highly active antiretroviral therapy (HAART), integrated HIV-1 proviral DNA cannot be eradicated from an infected individual. HAART is not able to eliminate latently infected cells that remain invisible to the immune system. Viral sanctuaries in specific tissues and immune-privileged sites may cause residual viral replication that contributes to HIV-1 persistence. The “Shock or Kick, and Kill” approach uses latency reversing agents (LRAs) in the presence of HAART, followed by cell-killing due to viral cytopathic effects and immune-mediated clearance. Different LRAs may be required for the in vivo reactivation of HIV-1 in different CD4⁺ T cell reservoirs, leading to the activation of cellular transcription factors acting on the integrated proviral HIV-1 LTR. An important requirement for LRA drugs is the reactivation of viral transcription and replication without causing a generalized immune activation. Toll-like receptors, RIG-I like receptors, and STING agonists have emerged recently as a new class of LRAs that augment selective apoptosis in reactivated T lymphocytes. The challenge is to extend in vitro observations to HIV-1 positive patients. Further studies are also needed to overcome the mechanisms that protect latently infected cells from reactivation and/or elimination by the immune system. The Block and Lock alternative strategy aims at using latency promoting/inducing agents (LPAs/LIAs) to block the ability of latent proviruses to reactivate transcription in order to achieve a long term lock down of potential residual virus replication. The Shock and Kill and the Block and Lock approaches may not be only alternative to each other, but, if combined together (one after the other), or given all at once [namely “Shoc-K(kill) and B(block)-Lock”], they may represent a better approach to a functional cure.

Clonal Expansion of Infected CD4+ T Cells in People Living with HIV

HIV infection is not curable with current antiretroviral therapy (ART) because a small fraction of CD4+ T cells infected prior to ART initiation persists. Understanding the nature of this latent reservoir and how it is created is essential to development of potentially curative strategies. The discovery that a large fraction of the persistently infected cells in individuals on suppressive ART are members of large clones greatly changed our view of the reservoir and how it arises. Rather than being the products of infection of resting cells, as was once thought, HIV persistence is largely or entirely a consequence of infection of cells that are either expanding or are destined to expand, primarily due to antigen-driven activation. Although most of the clones carry defective proviruses, some carry intact infectious proviruses; these clones comprise the majority of the reservoir. A large majority of both the defective and the intact infectious proviruses in clones of infected cells are transcriptionally silent; however, a small fraction expresses a few copies of unspliced HIV RNA. A much smaller fraction is responsible for production of low levels of infectious virus, which can rekindle infection when ART is stopped. Further understanding of the reservoir will be needed to clarify the mechanism(s) by which provirus expression is controlled in the clones of cells that constitute the reservoir.

Distinguishing between models of mammalian gene expression: telegraph-like models versus mechanistic models

Two-state models (telegraph-like models) have a successful history of predicting distributions of cellular and nascent mRNA numbers that can well fit experimental data. These models exclude key rate limiting steps, and hence it is unclear why they are able to accurately predict the number distributions. To answer this question, here we compare these models to a novel stochastic mechanistic model of transcription in mammalian cells that presents a unified description of transcriptional factor, polymerase and mature mRNA dynamics. We show that there is a large region of parameter space where the first, second and third moments of the distributions of the waiting times between two consecutively produced transcripts (nascent or mature) of two-state and mechanistic models exactly match. In this region: (i) one can uniquely express the two-state model parameters in terms of those of the mechanistic model, (ii) the models are practically indistinguishable by comparison of their transcript numbers distributions, and (iii) they are distinguishable from the shape of their waiting time distributions. Our results clarify the relationship between different gene expression models and identify a means to select between them from experimental data.

Quantitative imaging of transcription in living Drosophila embryos reveals the impact of core promoter motifs on promoter state dynamics

Genes are expressed in stochastic transcriptional bursts linked to alternating active and inactive promoter states. A major challenge in transcription is understanding how promoter composition dictates bursting, particularly in multicellular organisms. We investigate two key Drosophila developmental promoter motifs, the TATA box (TATA) and the Initiator (INR). Using live imaging in Drosophila embryos and new computational methods, we demonstrate that bursting occurs on multiple timescales ranging from seconds to minutes. TATA-containing promoters and INR-containing promoters exhibit distinct dynamics, with one or two separate rate-limiting steps respectively. A TATA box is associated with long active states, high rates of polymerase initiation, and short-lived, infrequent inactive states. In contrast, the INR motif leads to two inactive states, one of which relates to promoter-proximal polymerase pausing. Surprisingly, the model suggests pausing is not obligatory, but occurs stochastically for a subset of polymerases. Overall, our results provide a rationale for promoter switching during zygotic genome activation.