Transcriptional bursting from the HIV-1 promoter is a significant source of stochastic noise in HIV-1 gene expression

Stochastic Gene Expression in Proliferating Cells: Differing Noise Intensity in Single-Cell and Population Perspectives

Random fluctuations (noise) in gene expression can be studied from two complementary perspectives: following expression in a single cell over time or comparing expression between cells in a proliferating population at a given time. Here, we systematically investigated scenarios where both perspectives lead to different levels of noise in a given gene product. We first consider a stable protein, whose concentration is diluted by cellular growth, and the protein inhibits growth at high concentrations, establishing a positive feedback loop. For a stochastic model with molecular bursting of gene products, we analytically predict and contrast the steady-state distributions of protein concentration in both frameworks. Although positive feedback amplifies the noise in expression, this amplification is much higher in the population framework compared to following a single cell over time. We also study other processes that lead to different noise levels even in the absence of such dilution-based feedback. When considering randomness in the partitioning of molecules between daughters during mitosis, we find that in the single-cell perspective, the noise in protein concentration is independent of noise in the cell cycle duration. In contrast, partitioning noise is amplified in the population perspective by increasing randomness in cell-cycle time. Overall, our results show that the commonly used single-cell framework that does not account for proliferating cells can, in some cases, underestimate the noise in gene product levels. These results have important implications for studying the inter-cellular variation of different stress-related expression programs across cell types that are known to inhibit cellular growth.

The HIV-1 Transcriptional Program: From Initiation to Elongation Control

A large body of work in the last four decades has revealed the key pillars of HIV-1 transcription control at the initiation and elongation steps. Here, I provide a recount of this collective knowledge starting with the genomic elements (DNA and nascent TAR RNA stem-loop) and transcription factors (cellular and the viral transactivator Tat), and later transitioning to the assembly and regulation of transcription initiation and elongation complexes, and the role of chromatin structure. Compelling evidence support a core HIV-1 transcriptional program regulated by the sequential and concerted action of cellular transcription factors and Tat to promote initiation and sustain elongation, highlighting the efficiency of a small virus to take over its host to produce the high levels of transcription required for viral replication. I summarize new advances including the use of CRISPR-Cas9, genetic tools for acute factor depletion, and imaging to study transcriptional dynamics, bursting and the progression through the multiple phases of the transcriptional cycle. Finally, I describe current challenges to future major advances and discuss areas that deserve more attention to both bolster our basic knowledge of the core HIV-1 transcriptional program and open up new therapeutic opportunities.

N-terminal tagging of RNA Polymerase II shapes transcriptomes more than C-terminal alterations

RNA polymerase II (Pol II) has a C-terminal domain (CTD) that is unstructured, consisting of a large number of heptad repeats, and whose precise function remains unclear. Here, we investigate how altering the CTD's length and fusing it with protein tags affects transcriptional output on a genome-wide scale in mammalian cells at single-cell resolution. While transcription generally appears to occur in burst-like fashion, where RNA is predominantly made during short bursts of activity that are interspersed with periods of transcriptional silence, the CTD's role in shaping these dynamics seems gene-dependent; global patterns of bursting appear mostly robust to CTD alterations. Introducing protein tags with defined structures to the N terminus cause transcriptome-wide effects, however. We find the type of tag to dominate characteristics of the resulting transcriptomes. This is possibly due to Pol II-interacting factors, including non-coding RNAs, whose expression correlates with the tags. Proteins involved in liquid-liquid phase separation appear prominently.

Gram-positive bacteria are primed for surviving lethal doses of antibiotics and chemical stress

Antibiotic resistance kills millions worldwide yearly. However, a major contributor to recurrent infections lies in a small fraction of bacterial cells, known as persisters. These cells are not inherently antibiotic-resistant, yet they lead to increased antibiotic usage, raising the risk of developing resistant progenies. In a bacterial population, individual cells exhibit considerable fluctuations in their gene expression levels despite being cultivated under identical, stable conditions. This variability in cell-to-cell characteristics (phenotypic diversity) within an isogenic population enables persister cells to withstand antibiotic exposure by entering a non-dividing state. We recently showed the existence of "primed cells" in E. coli. Primed cells are dividing cells prepared for antibiotic stress before encountering it and are more prone to form persisters. They also pass their "prepared state" down for several generations through epigenetic memory. Here, we show that primed cells are common among distant bacterial lineages, allowing for survival against antibiotics and other chemical stress, and form in different growth phases. They are also responsible for increased persister levels in transition and stationary phases compared to the log phase. We tested and showed that the Gram-positive bacterium Bacillus megaterium, evolutionarily very distant from E. coli, forms primed cells and has a transient epigenetic memory that is maintained for 7 generations or more. We showed this using ciprofloxacin and the non-antibiotic chemical stress fluoride. It is well established that persister levels are higher in the stationary phase than in the log phase, and B. megaterium persisters levels are nearly identical from the early to late-log phase but are ~2-fold and ~4-fold higher in the transition and stationary phase, respectively. It was previously proposed that there are two distinct types of persisters: Type II forms in the log phase, while Type I forms in the stationary phase. However, we show that primed cells lead to increased persisters in the transition and stationary phase and found no evidence of Type I or II persisters with distant phenotypes. Overall, we have provided substantial evidence of the importance of primed cells and their transitory epigenetic memories to surviving stress.

Single-molecule RNA-FISH analysis reveals stochasticity in reactivation of latent HIV-1 regulated by Nuclear Orphan Receptors NR4A and cMYC

HIV-1 eradication strategies require complete reactivation of HIV-1 latent cells by Latency Reversing Agents (LRA). Current methods lack effectiveness due to incomplete proviral reactivation. We employed a single-molecule RNA-FISH (smRNA-FISH) and FISH-Quant analysis and found that proviral reactivation is highly variable from cell-to-cell, stochastic, and occurs in bursts and waves, with different kinetics in response to diverse LRAs. Approximately 1-5% of latent cells exhibited stochastic reactivation without LRAs. Through single-cell RNA-seq analysis, we identified NR4A3 and cMYC as extrinsic factors associated with stochastic HIV-1 reactivation. Concomitant with HIV-1 reactivation cMYC was downregulated and NR4A3 was upregulated in both latent cell lines and primary CD4 + T-cells from aviremic patients. By inhibiting cMYC using SN-38, an active metabolite of irinotecan, we induced NR4A3 and HIV-1 expression. Our results suggest that inherent stochasticity in proviral reactivation contributes to cell-to-cell variability, which could potentially be modulated by drugs targeting cMYC and NR4A3.

Phase separation reduces cell-to-cell variability of transcriptional bursting

Gene expression is a stochastic and noisy process often occurring in "bursts". Experiments have shown that the compartmentalization of proteins by liquid-liquid phase separation is conducive to reducing the noise of gene expression. Therefore, an important goal is to explore the role of bursts in phase separation noise reduction processes. We propose a coupled model that includes phase separation and a two-state gene expression process. Using the timescale separation method, we obtain approximate solutions for the expectation, variance, and noise strength of the dilute phase. We find that a higher burst frequency weakens the ability of noise reduction by phase separation, but as the burst size increases, this ability first increases and then decreases. This study provides a deeper understanding of phase separation to reduce noise in the stochastic gene expression with burst kinetics.

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.

CDK8 inhibitors antagonize HIV-1 reactivation and promote provirus latency in T cells

Latent HIV-1 provirus represents the barrier toward a cure for infection and is dependent upon the host RNA Polymerase (Pol) II machinery for reemergence. Here, we find that inhibitors of the RNA Pol II mediator kinases CDK8/19, Senexin A and BRD6989, inhibit induction of HIV-1 expression in response to latency-reversing agents and T cell signaling agonists. These inhibitors were found to impair recruitment of RNA Pol II to the HIV-1 LTR. Furthermore, HIV-1 expression in response to several latency reversal agents was impaired upon disruption of CDK8 by shRNA or gene knockout. However, the effects of CDK8 depletion did not entirely mimic CDK8/19 kinase inhibition suggesting that the mediator kinases are not functionally redundant. Additionally, treatment of CD4+ peripheral blood mononuclear cells isolated from people living with HIV-1 and who are receiving antiretroviral therapy with Senexin A inhibited induction of viral replication in response to T cell stimulation by PMA and ionomycin. These observations indicate that the mediator kinases, CDK8 and CDK19, play a significant role for regulation of HIV-1 transcription and that small molecule inhibitors of these enzymes may contribute to therapies designed to promote deep latency involving the durable suppression of provirus expression. IMPORTANCE A cure for HIV-1 infection will require novel therapies that can force elimination of cells that contain copies of the virus genome inserted into the cell chromosome, but which is shut off, or silenced. These are known as latently-infected cells, which represent the main reason why current treatment for HIV/AIDS cannot cure the infection because the virus in these cells is unaffected by current drugs. Our results indicate that chemical inhibitors of Cdk8 also inhibit the expression of latent HIV provirus. Cdk8 is an important enzyme that regulates the expression of genes in response to signals to which cells need to respond and which is produced by a gene that is frequently mutated in cancers. Our observations indicate that Cdk8 inhibitors may be employed in novel therapies to prevent expression from latent provirus, which might eventually enable infected individuals to cease treatment with antiretroviral drugs.

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.

Evaluating single-cell variability in proteasomal decay

Gene expression is a stochastic process that leads to variability in mRNA and protein abundances even within an isogenic population of cells grown in the same environment. This variation, often called gene-expression noise, has typically been attributed to transcriptional and translational processes while ignoring the contributions of protein decay variability across cells. Here we estimate the single-cell protein decay rates of two degron GFPs in Saccharomyces cerevisiae using time-lapse microscopy. We find substantial cell-to-cell variability in the decay rates of the degron GFPs. We evaluate cellular features that explain the variability in the proteasomal decay and find that the amount of 20s catalytic beta subunit of the proteasome marginally explains the observed variability in the degron GFP half-lives. We propose alternate hypotheses that might explain the observed variability in the decay of the two degron GFPs. Overall, our study highlights the importance of studying the kinetics of the decay process at single-cell resolution and that decay rates vary at the single-cell level, and that the decay process is stochastic. A complex model of decay dynamics must be included when modeling stochastic gene expression to estimate gene expression noise.

Analysis and design of single-cell experiments to harvest fluctuation information while rejecting measurement noise

Introduction: Despite continued technological improvements, measurement errors always reduce or distort the information that any real experiment can provide to quantify cellular dynamics. This problem is particularly serious for cell signaling studies to quantify heterogeneity in single-cell gene regulation, where important RNA and protein copy numbers are themselves subject to the inherently random fluctuations of biochemical reactions. Until now, it has not been clear how measurement noise should be managed in addition to other experiment design variables (e.g., sampling size, measurement times, or perturbation levels) to ensure that collected data will provide useful insights on signaling or gene expression mechanisms of interest. Methods: We propose a computational framework that takes explicit consideration of measurement errors to analyze single-cell observations, and we derive Fisher Information Matrix (FIM)-based criteria to quantify the information value of distorted experiments. Results and Discussion: We apply this framework to analyze multiple models in the context of simulated and experimental single-cell data for a reporter gene controlled by an HIV promoter. We show that the proposed approach quantitatively predicts how different types of measurement distortions affect the accuracy and precision of model identification, and we demonstrate that the effects of these distortions can be mitigated through explicit consideration during model inference. We conclude that this reformulation of the FIM could be used effectively to design single-cell experiments to optimally harvest fluctuation information while mitigating the effects of image distortion.

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.

Probing transient memory of cellular states using single-cell lineages

The inherent stochasticity in the gene product levels can drive single cells within an isoclonal population to different phenotypic states. The dynamic nature of this intercellular variation, where individual cells can transition between different states over time, makes it a particularly hard phenomenon to characterize. We reviewed recent progress in leveraging the classical Luria-Delbrück experiment to infer the transient heritability of the cellular states. Similar to the original experiment, individual cells were first grown into cell colonies, and then, the fraction of cells residing in different states was assayed for each colony. We discuss modeling approaches for capturing dynamic state transitions in a growing cell population and highlight formulas that identify the kinetics of state switching from the extent of colony-to-colony fluctuations. The utility of this method in identifying multi-generational memory of the both expression and phenotypic states is illustrated across diverse biological systems from cancer drug resistance, reactivation of human viruses, and cellular immune responses. In summary, this fluctuation-based methodology provides a powerful approach for elucidating cell-state transitions from a single time point measurement, which is particularly relevant in situations where measurements lead to cell death (as in single-cell RNA-seq or drug treatment) or cause an irreversible change in cell physiology.

Enhanced Transcriptional Strength of HIV-1 Subtype C Minimizes Gene Expression Noise and Confers Stability to the Viral Latent State

Stochastic fluctuations in gene expression emanating from the HIV-1 long terminal repeat (LTR), amplified by the Tat positive feedback circuit, determine the choice between viral infection fates: active transcription (ON) or transcriptional silence (OFF). The emergence of several transcription factor binding site (TFBS) variant strains in HIV-1 subtype C (HIV-1C), especially those containing the duplication of the NF-κB motif, mandates the evaluation of the effect of enhanced transcriptional strength on gene expression noise and its influence on viral fate selection switch. Using a panel of subgenomic LTR-variant strains containing different copy numbers of the NF-κB motif (ranging from 0 to 4), we used flow cytometry, mRNA quantification, and pharmacological perturbations to demonstrate an inverse correlation between promoter strength and gene expression noise in Jurkat T cells and primary CD4⁺ T cells. The inverse correlation is consistent in clonal cell populations at constant intracellular concentrations of Tat and when NF-κB levels were regulated pharmacologically. Further, we show that strong LTRs containing at least two copies of the NF-κB motif in the enhancer establish a more stable latent state and demonstrate more rapid latency reversal than weak LTRs containing fewer motifs. We also demonstrate a cooperative binding of NF-κB to the motif cluster in HIV-1C LTRs containing two, three, or four NF-κB motifs (Hill coefficient [H] = 2.61, 3.56, and 3.75, respectively). The present work alludes to a possible evolution of the HIV-1C LTR toward gaining transcriptional strength associated with attenuated gene expression noise with implications for viral latency. IMPORTANCE Over the past two consecutive decades, HIV-1 subtype C (HIV-1C) has been undergoing directional evolution toward augmenting the transcriptional strength of the long terminal repeat (LTR) by adding more copies of the existing transcription factor binding site (TFBS) by sequence duplication. Additionally, the duplicated elements are genetically diverse, suggesting broader-range signal receptivity by variant LTRs. The HIV-1 promoter is inherently noisy, and the stochastic fluctuations in gene expression of variant LTRs may influence the active transcription (ON)/transcriptional silence (OFF) latency decisions. The evolving NF-κB motif variations of HIV-1C offer a powerful opportunity to examine how the transcriptional strength of the LTR might influence gene expression noise. Our work here shows that the augmented transcriptional strength of the HIV-1C LTR leads to concomitantly reduced gene expression noise, consequently leading to stabler latency maintenance and rapid latency reversal. The present work offers a novel lead toward appreciating the molecular mechanisms governing HIV-1 latency.

An HIV stochastic model with cell-to-cell infection, B-cell immune response and distributed delay

In this study, a delayed HIV stochastic model with virus-to-cell infection, cell-to-cell transmission and B-cell immune response is proposed. We first transform the stochastic differential equation with distributed delay into a high-dimensional degenerate stochastic differential equation, and then theoretically analyze the dynamic behaviour of the degenerate model. The unique global solution of the model is given by rigorous analysis. By formulating suitable Lyapunov functions, the existence of the stationary Markov process is obtained if the stochastic B-cell-activated reproduction number is greater than one. We also use the law of large numbers theorem and the spectral radius analysis method to deduce that the virus can be cleared if the stochastic B-cell-inactivated reproduction number is less than one. Through uncertainty and sensitivity analysis, we obtain key parameters that determine the value of the stochastic B-cell-activated reproduction number. Numerically, we examine that low level noise can maintain the number of the virus and B-cell populations at a certain range, while high level noise is helpful for the elimination of the virus. Furthermore, the effect of the cell-to-cell infection on model behaviour, and the influence of the key parameters on the size of the stochastic B-cell-activated reproduction number are also investigated.

The effect of noise in an HIV infection model with cytotoxic T-lymphocyte impairment

The human immunodeficiency virus (HIV) interacts with the immune cells within the human body, where the environment is uncertain and noisy. Stochastic models can successfully encapsulate the effect of such a noisy environment compared to their deterministic counterparts. The human immune system is complex but well-coordinated with various immune cells like C D ⁴T cells, dendritic cells, and cytotoxic T-lymphocyte (CTL) cells, among many others. The CTL can kill the antigenic cells after its recognition. However, the efficacy of CTL in removing the infected C D ⁴T cells is progressively compromised in HIV-infected individuals. This paper considers a noise-induced HIV-immune cell interaction model with immune impairment. A multiplicative white noise is introduced in the infection rate parameter to represent the fluctuations around the average value of the rate parameter as a causative effect of the noise. We analyzed the deterministic and stochastic models and prescribed sufficient conditions for infection eradication and persistence. It is determined under what parametric restrictions the asymptotic solutions of the noise-induced system will be a limiting case of the deterministic solutions. Simulation results revealed that the solutions of the deterministic system either converge to a CTL-dominated interior equilibrium or a CTL-free immunodeficient equilibrium, depending on the initial values of the system. Stochastic analysis divulged that higher noise might be helpful in the infection removal process. The extinction time of infected C D ⁴T cells for some fixed immune impairment gradually decreases with increasing noise intensity and follows the power law.

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.

Persistent Correlation in Cellular Noise Determines Longevity of Viral Infections

The slowly decaying viral dynamics, even after 2-3 weeks from diagnosis, is one of the characteristics of COVID-19 infection that is still unexplored in theoretical and experimental studies. This long-lived characteristic of viral infections in the framework of inherent variations or noise present at the cellular level is often overlooked. Therefore, in this work, we aim to understand the effect of these variations by proposing a stochastic non-Markovian model that not only captures the coupled dynamics between the immune cells and the virus but also enables the study of the effect of fluctuations. Numerical simulations of our model reveal that the long-range temporal correlations in fluctuations dictate the long-lived dynamics of a viral infection and, in turn, also affect the rates of immune response. Furthermore, predictions of our model system are in agreement with the experimental viral load data of COVID-19 patients from various countries.

Activation domains can decouple the mean and noise of gene expression

Regulatory mechanisms set a gene's average level of expression, but a gene's expression constantly fluctuates around that average. These stochastic fluctuations, or expression noise, play a role in cell-fate transitions, bet hedging in microbes, and the development of chemotherapeutic resistance in cancer. An outstanding question is what regulatory mechanisms contribute to noise. Here, we demonstrate that, for a fixed mean level of expression, strong activation domains (ADs) at low abundance produce high expression noise, while weak ADs at high abundance generate lower expression noise. We conclude that differences in noise can be explained by the interplay between a TF's nuclear concentration and the strength of its AD's effect on mean expression, without invoking differences between classes of ADs. These results raise the possibility of engineering gene expression noise independently of mean levels in synthetic biology contexts and provide a potential mechanism for natural selection to tune the noisiness of gene expression.

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.

The nonequilibrium mechanism of noise-enhanced drug synergy in HIV latency reactivation

Noise-modulating chemicals can synergize with transcriptional activators in reactivating latent HIV to eliminate latent HIV reservoirs. To understand the underlying biomolecular mechanism, we investigate a previous two-gene-state model and identify two necessary conditions for the synergy: an assumption of the inhibition effect of transcription activators on noise enhancers; and frequent transitions to the gene non-transcription-permissive state. We then develop a loop-four-gene-state model with Tat transcription/translation and find that drug synergy is mainly determined by the magnitude and direction of energy input into the genetic regulatory kinetics of the HIV promoter. The inhibition effect of transcription activators is actually a phenomenon of energy dissipation in the nonequilibrium gene transition system. Overall, the loop-four-state model demonstrates that energy dissipation plays a crucial role in HIV latency reactivation, which might be useful for improving drug effects and identifying other synergies on lentivirus latency reactivation.

•

The inhibition of Activator on Noise enhancer is necessary for their synergy in reactivating HIV

•

The drug synergy is a nonequilibrium phenomenon in the gene regulatory system

•

The magnitude and direction of energy input determine the drug synergy

•

This nonequilibrium mechanism is general without regarding molecular details

Parallel analysis of transcription, integration, and sequence of single HIV-1 proviruses

HIV-1-infected cells that persist despite antiretroviral therapy (ART) are frequently considered "transcriptionally silent," but active viral gene expression may occur in some cells, challenging the concept of viral latency. Applying an assay for profiling the transcriptional activity and the chromosomal locations of individual proviruses, we describe a global genomic and epigenetic map of transcriptionally active and silent proviral species and evaluate their longitudinal evolution in persons receiving suppressive ART. Using genome-wide epigenetic reference data, we show that proviral transcriptional activity is associated with activating epigenetic chromatin features in linear proximity of integration sites and in their inter- and intrachromosomal contact regions. Transcriptionally active proviruses were actively selected against during prolonged ART; however, this pattern was violated by large clones of virally infected cells that may outcompete negative selection forces through elevated intrinsic proliferative activity. Our results suggest that transcriptionally active proviruses are dynamically evolving under selection pressure by host factors.

Emerging roles for endogenous retroviruses in immune epigenetic regulation

In recent years, there has been significant progress toward understanding the transcriptional networks underlying mammalian immune responses, fueled by advances in regulatory genomic technologies. Epigenomic studies profiling immune cells have generated detailed genome-wide maps of regulatory elements that will be key to deciphering the regulatory networks underlying cellular immune responses and autoimmune disorders. Unbiased analyses of these genomic maps have uncovered endogenous retroviruses as an unexpected ally in the regulation of human immune systems. Despite their parasitic origins, studies are finding an increasing number of examples of retroviral sequences having been co-opted for beneficial immune function and regulation by the host cell. Here, we review how endogenous retroviruses have given rise to numerous regulatory elements that shape the epigenetic landscape of host immune responses. We will discuss the implications of these elements on the function, dysfunction, and evolution of innate immunity.

Research priorities for an HIV cure: International AIDS Society Global Scientific Strategy 2021

Despite the success of antiretroviral therapy (ART) for people living with HIV, lifelong treatment is required and there is no cure. HIV can integrate in the host genome and persist for the life span of the infected cell. These latently infected cells are not recognized as foreign because they are largely transcriptionally silent, but contain replication-competent virus that drives resurgence of the infection once ART is stopped. With a combination of immune activators, neutralizing antibodies, and therapeutic vaccines, some nonhuman primate models have been cured, providing optimism for these approaches now being evaluated in human clinical trials. In vivo delivery of gene-editing tools to either target the virus, boost immunity or protect cells from infection, also holds promise for future HIV cure strategies. In this Review, we discuss advances related to HIV cure in the last 5 years, highlight remaining knowledge gaps and identify priority areas for research for the next 5 years.

PRC1 drives Polycomb-mediated gene repression by controlling transcription initiation and burst frequency

The Polycomb repressive system plays a fundamental role in controlling gene expression during mammalian development. To achieve this, Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) bind target genes and use histone modification-dependent feedback mechanisms to form Polycomb chromatin domains and repress transcription. The inter-relatedness of PRC1 and PRC2 activity at these sites has made it difficult to discover the specific components of Polycomb chromatin domains that drive gene repression and to understand mechanistically how this is achieved. Here, by exploiting rapid degron-based approaches and time-resolved genomics, we kinetically dissect Polycomb-mediated repression and discover that PRC1 functions independently of PRC2 to counteract RNA polymerase II binding and transcription initiation. Using single-cell gene expression analysis, we reveal that PRC1 acts uniformly within the cell population and that repression is achieved by controlling transcriptional burst frequency. These important new discoveries provide a mechanistic and conceptual framework for Polycomb-dependent transcriptional control.

Monitoring reactivation of latent HIV by label-free gradient light interference microscopy

Human immunodeficiency virus (HIV) can infect cells and take a quiescent and nonexpressive state called latency. In this study, we report insights provided by label-free, gradient light interference microscopy (GLIM) about the changes in dry mass, diameter, and dry mass density associated with infected cells that occur upon reactivation. We discovered that the mean cell dry mass and mean diameter of latently infected cells treated with reactivating drug, TNF-α, are higher for latent cells that reactivate than those of the cells that did not reactivate. Cells with mean dry mass and diameter less than approximately 10 pg and 8 μm, respectively, remain exclusively in the latent state. Also, cells with mean dry mass greater than approximately 28-30 pg and mean diameter greater than 11–12 μm have a higher probability of reactivating. This study is significant as it presents a new label-free approach to quantify latent reactivation of a virus in single cells.

•

GLIM imaging reveals differences between latent and reactivated HIV in JLat cells

•

Cells with reactivated HIV have higher dry mass and diameter

Noise suppression in stochastic genetic circuits using PID controllers

Inside individual cells, protein population counts are subject to molecular noise due to low copy numbers and the inherent probabilistic nature of biochemical processes. We investigate the effectiveness of proportional, integral and derivative (PID) based feedback controllers to suppress protein count fluctuations originating from two noise sources: bursty expression of the protein, and external disturbance in protein synthesis. Designs of biochemical reactions that function as PID controllers are discussed, with particular focus on individual controllers separately, and the corresponding closed-loop system is analyzed for stochastic controller realizations. Our results show that proportional controllers are effective in buffering protein copy number fluctuations from both noise sources, but this noise suppression comes at the cost of reduced static sensitivity of the output to the input signal. In contrast, integral feedback has no effect on the protein noise level from stochastic expression, but significantly minimizes the impact of external disturbances, particularly when the disturbance comes at low frequencies. Counter-intuitively, integral feedback is found to amplify external disturbances at intermediate frequencies. Next, we discuss the design of a coupled feedforward-feedback biochemical circuit that approximately functions as a derivate controller. Analysis using both analytical methods and Monte Carlo simulations reveals that this derivative controller effectively buffers output fluctuations from bursty stochastic expression, while maintaining the static input-output sensitivity of the open-loop system. In summary, this study provides a systematic stochastic analysis of biochemical controllers, and paves the way for their synthetic design and implementation to minimize deleterious fluctuations in gene product levels.

In the noisy cellular environment, biochemical species such as genes, RNAs and proteins that often occur at low molecular counts, are subject to considerable stochastic fluctuations in copy numbers over time. How cellular biochemical processes function reliably in the face of such randomness is an intriguing fundamental problem. Increasing evidence suggests that random fluctuations (noise) in protein copy numbers play important functional roles, such as driving genetically identical cells to different cell fates. Moreover, many disease states have been attributed to elevated noise levels in specific proteins. Here we systematically investigate design of biochemical systems that function as proportional, integral and derivative-based feedback controllers to suppress protein count fluctuations arising from bursty expression of the protein and external disturbance in protein synthesis. Our results show that different controllers are effective in buffering different noise components, and identify ranges of feedback gain for minimizing deleterious fluctuations in protein levels.

A near analytic solution of a stochastic immune response model considering variability in virus and T-cell dynamics

Biological processes at the cellular level are stochastic in nature, and the immune response system is no different. Therefore, models that attempt to explain this system need to also incorporate noise or fluctuations that can account for the observed variability. In this work, a stochastic model of the immune response system is presented in terms of the dynamics of T cells and virus particles. Making use of the Green's function and the Wilemski-Fixman approximation, this model is then solved to obtain the analytical expression for the joint probability density function of these variables in the early and late stages of infection. This is then also used to calculate the average level of virus particles in the system. Upon comparing the theoretically predicted average virus levels to those of COVID-19 patients, it is hypothesized that the long-lived dynamics that are characteristics of such viral infections are due to the long range correlations in the temporal fluctuations of the virions. This model, therefore, provides an insight into the effects of noise on viral dynamics.

TNF stimulation primarily modulates transcriptional burst size of NF-κB-regulated genes

Cell-to-cell heterogeneity is a feature of the tumor necrosis factor (TNF)-stimulated inflammatory response mediated by the transcription factor NF-κB, motivating an exploration of the underlying sources of this noise. Here, we combined single-transcript measurements with computational models to study transcriptional noise at six NF-κB-regulated inflammatory genes. In the basal state, NF-κB-target genes displayed an inverse correlation between mean and noise characteristic of transcriptional bursting. By analyzing transcript distributions with a bursting model, we found that TNF primarily activated transcription by increasing burst size while maintaining burst frequency for gene promoters with relatively high basal histone 3 acetylation (AcH3) that marks open chromatin environments. For promoters with lower basal AcH3 or when AcH3 was decreased with a small molecule drug, the contribution of burst frequency to TNF activation increased. Finally, we used a mathematical model to show that TNF positive feedback amplified gene expression noise resulting from burst size-mediated transcription, leading to a subset of cells with high TNF protein expression. Our results reveal potential sources of noise underlying intercellular heterogeneity in the TNF-mediated inflammatory response.

Synergistic Chromatin-Modifying Treatments Reactivate Latent HIV and Decrease Migration of Multiple Host-Cell Types

Upon infection of its host cell, human immunodeficiency virus (HIV) establishes a quiescent and non-productive state capable of spontaneous reactivation. Diverse cell types harboring the provirus form a latent reservoir, constituting a major obstacle to curing HIV. Here, we investigate the effects of latency reversal agents (LRAs) in an HIV-infected THP-1 monocyte cell line in vitro. We demonstrate that leading drug treatments synergize activation of the HIV long terminal repeat (LTR) promoter. We establish a latency model in THP-1 monocytes using a replication incompetent HIV reporter vector with functional Tat, and show that chromatin modifiers synergize with a potent transcriptional activator to enhance HIV reactivation, similar to T-cells. Furthermore, leading reactivation cocktails are shown to differentially affect latency reactivation and surface expression of chemokine receptor type 4 (CXCR4), leading to altered host cell migration. This study investigates the effect of chromatin-modifying LRA treatments on HIV latent reactivation and cell migration in monocytes. As previously reported in T-cells, epigenetic mechanisms in monocytes contribute to controlling the relationship between latent reactivation and cell migration. Ultimately, advanced “Shock and Kill” therapy needs to successfully target and account for all host cell types represented in a complex and composite latency milieu.

Correlation between external regulators governs the mean-noise relationship in stochastic gene expression

Gene transcription in single cells is inherently a probabilistic process. The relationship between variance ($ \sigma^{2} $) and mean expression ($ \mu $) is of paramount importance for investigations into the evolutionary origins and consequences of noise in gene expression. It is often formulated as $ \log \left({{{\sigma}^{2}}}/{{{\mu}^{2}}}\; \right) = \beta\log\mu+\log\alpha $, where $ \beta $ is a key parameter since its sign determines the qualitative dependence of noise on mean. We reveal that the sign of $ \beta $ is controlled completely by external regulation, but independent of promoter structure. Specifically, it is negative if regulators as stochastic variables are independent but positive if they are correlated. The essential mechanism revealed here can well interpret diverse experimental phenomena underlying expression noise. Our results imply that external regulation rather than promoter sequence governs the mean-noise relationship.

Dynamical behavior of a stochastic multigroup staged-progression HIV model with saturated incidence rate and higher-order perturbations

In this paper, we analyze a higher-order stochastically perturbed multigroup staged-progression model for the transmission of HIV with saturated incidence rate. We obtain sufficient conditions for the existence and uniqueness of an ergodic stationary distribution of positive solutions to the system by establishing a suitable stochastic Lyapunov function. In addition, we make up adequate conditions for complete eradication and wiping out the infectious disease. In a biological interpretation, the existence of a stationary distribution implies that the disease will prevail and persist in the long term. Finally, examples and numerical simulations are introduced to validate our theoretical results.

Mathematical analysis and simulation of a stochastic COVID-19 Lévy jump model with isolation strategy

This paper investigates the dynamics of a COVID-19 stochastic model with isolation strategy. The white noise as well as the Lévy jump perturbations are incorporated in all compartments of the suggested model. First, the existence and uniqueness of a global positive solution are proven. Next, the stochastic dynamic properties of the stochastic solution around the deterministic model equilibria are investigated. Finally, the theoretical results are reinforced by some numerical simulations.

A transient heritable memory regulates HIV reactivation from latency

Reactivation of human immunodeficiency virus 1 (HIV-1) from latently infected T cells is a critical barrier to cure patients. It remains unknown whether reactivation of individual latent cells occurs stochastically in response to latency reversal agents (LRAs) or is a deterministic outcome of an underlying cell state. To characterize these single-cell responses, we leverage the classical Luria-Delbrück fluctuation test where single cells are isolated from a clonal population and exposed to LRAs after colony expansion. Data show considerable colony-to-colony fluctuations with the fraction of reactivating cells following a skewed distribution. Modeling systematic measurements of fluctuations over time uncovers a transient heritable memory that regulates HIV-1 reactivation, where single cells are in an LRA-responsive state for a few weeks before switching back to an irresponsive state. These results have enormous implications for designing therapies to purge the latent reservoir and further utilize fluctuation-based assays to uncover hidden transient cellular states underlying phenotypic heterogeneity.

3 '-5 ' crosstalk contributes to transcriptional bursting

BACKGROUND

Transcription in mammalian cells is a complex stochastic process involving shuttling of polymerase between genes and phase-separated liquid condensates. It occurs in bursts, which results in vastly different numbers of an mRNA species in isogenic cell populations. Several factors contributing to transcriptional bursting have been identified, usually classified as intrinsic, in other words local to single genes, or extrinsic, relating to the macroscopic state of the cell. However, some possible contributors have not been explored yet. Here, we focus on processes at the 3 ' and 5 ' ends of a gene that enable reinitiation of transcription upon termination.

RESULTS

Using Bayesian methodology, we measure the transcriptional bursting in inducible transgenes, showing that perturbation of polymerase shuttling typically reduces burst size, increases burst frequency, and thus limits transcriptional noise. Analysis based on paired-end tag sequencing (PolII ChIA-PET) suggests that this effect is genome wide. The observed noise patterns are also reproduced by a generative model that captures major characteristics of the polymerase flux between the ends of a gene and a phase-separated compartment.

CONCLUSIONS

Interactions between the 3 ' and 5 ' ends of a gene, which facilitate polymerase recycling, are major contributors to transcriptional noise.

Identifying chromatin features that regulate gene expression distribution

Gene expression variability, differences in the number of mRNA per cell across a population of cells, is ubiquitous across diverse organisms with broad impacts on cellular phenotypes. The role of chromatin in regulating average gene expression has been extensively studied. However, what aspects of the chromatin contribute to gene expression variability is still underexplored. Here we addressed this problem by leveraging chromatin diversity and using a systematic investigation of randomly integrated expression reporters to identify what aspects of chromatin microenvironment contribute to gene expression variability. Using DNA barcoding and split-pool decoding, we created a large library of isogenic reporter clones and identified reporter integration sites in a massive and parallel manner. By mapping our measurements of reporter expression at different genomic loci with multiple epigenetic profiles including the enrichment of transcription factors and the distance to different chromatin states, we identified new factors that impact the regulation of gene expression distributions.

Microbial metabolic noise

From the time a cell was first placed under the microscope, it became apparent that identifying two clonal cells that "look" identical is extremely challenging. Since then, cell-to-cell differences in shape, size, and protein content have been carefully examined, informing us of the ultimate limits that hinder two cells from occupying an identical phenotypic state. Here, we present recent experimental and computational evidence that similar limits emerge also in cellular metabolism. These limits pertain to stochastic metabolic dynamics and, thus, cell-to-cell metabolic variability, including the resulting adapting benefits. We review these phenomena with a focus on microbial metabolism and conclude with a brief outlook on the potential relationship between metabolic noise and adaptive evolution. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering.

Uncoupling gene expression noise along the central dogma using genome engineered human cell lines

Eukaryotic protein synthesis is an inherently stochastic process. This stochasticity stems not only from variations in cell content between cells but also from thermodynamic fluctuations in a single cell. Ultimately, these inherently stochastic processes manifest as noise in gene expression, where even genetically identical cells in the same environment exhibit variation in their protein abundances. In order to elucidate the underlying sources that contribute to gene expression noise, we quantify the contribution of each step within the process of protein synthesis along the central dogma. We uncouple gene expression at the transcriptional, translational, and post-translational level using custom engineered circuits stably integrated in human cells using CRISPR. We provide a generalized framework to approximate intrinsic and extrinsic noise in a population of cells expressing an unbalanced two-reporter system. Our decomposition shows that the majority of intrinsic fluctuations stem from transcription and that coupling the two genes along the central dogma forces the fluctuations to propagate and accumulate along the same path, resulting in increased observed global correlation between the products.

Mixture distributions in a stochastic gene expression model with delayed feedback: a WKB approximation approach

Noise in gene expression can be substantively affected by the presence of production delay. Here we consider a mathematical model with bursty production of protein, a one-step production delay (the passage of which activates the protein), and feedback in the frequency of bursts. We specifically focus on examining the steady-state behaviour of the model in the slow-activation (i.e. large-delay) regime. Using a formal asymptotic approach, we derive an autonomous ordinary differential equation for the inactive protein that applies in the slow-activation regime. If the differential equation is monostable, the steady-state distribution of the inactive (active) protein is approximated by a single Gaussian (Poisson) mode located at the globally stable fixed point of the differential equation. If the differential equation is bistable (due to cooperative positive feedback), the steady-state distribution of the inactive (active) protein is approximated by a mixture of Gaussian (Poisson) modes located at the stable fixed points; the weights of the modes are determined from a WKB approximation to the stationary distribution. The asymptotic results are compared to numerical solutions of the chemical master equation.

Enhancement of gene expression noise from transcription factor binding to genomic decoy sites

The genome contains several high-affinity non-functional binding sites for transcription factors (TFs) creating a hidden and unexplored layer of gene regulation. We investigate the role of such “decoy sites” in controlling noise (random fluctuations) in the level of a TF that is synthesized in stochastic bursts. Prior studies have assumed that decoy-bound TFs are protected from degradation, and in this case decoys function to buffer noise. Relaxing this assumption to consider arbitrary degradation rates for both bound/unbound TF states, we find rich noise behaviors. For low-affinity decoys, noise in the level of unbound TF always monotonically decreases to the Poisson limit with increasing decoy numbers. In contrast, for high-affinity decoys, noise levels first increase with increasing decoy numbers, before decreasing back to the Poisson limit. Interestingly, while protection of bound TFs from degradation slows the time-scale of fluctuations in the unbound TF levels, the decay of bound TFs leads to faster fluctuations and smaller noise propagation to downstream target proteins. In summary, our analysis reveals stochastic dynamics emerging from nonspecific binding of TFs and highlights the dual role of decoys as attenuators or amplifiers of gene expression noise depending on their binding affinity and stability of the bound TF.

'Rinse and Replace': Boosting T Cell Turnover To Reduce HIV-1 Reservoirs

Latent HIV-1 persists indefinitely during antiretroviral therapy (ART) as an integrated silent genome in long-lived memory CD4+ T cells. In untreated infections, immune activation increases the turnover of intrinsically long-lived provirus-containing CD4+ T cells. Those are 'washed out' as a result of their activation, which when coupled to viral protein expression can facilitate local inflammation and recruitment of uninfected cells to activation sites, causing latently infected cells to compete for survival. De novo infection can counter this washout. During ART, inflammation and CD4+ T cell activation wane, resulting in reduced cell turnover and a persistent reservoir. We propose accelerating reservoir washout during ART by triggering sequential waves of polyclonal CD4+ T cell activation while simultaneously enhancing virus protein expression. Reservoir reduction as an adjunct to other therapies might achieve lifelong viral control.

HIV-1 Proviral Transcription and Latency in the New Era

Three decades of extensive work in the HIV field have revealed key viral and host cell factors controlling proviral transcription. Various models of transcriptional regulation have emerged based on the collective information from in vitro assays and work in both immortalized and primary cell-based models. Here, we provide a recount of the past and current literature, highlight key regulatory aspects, and further describe potential limitations of previous studies. We particularly delve into critical steps of HIV gene expression including the role of the integration site, nucleosome positioning and epigenomics, and the transition from initiation to pausing and pause release. We also discuss open questions in the field concerning the generality of previous regulatory models to the control of HIV transcription in patients under suppressive therapy, including the role of the heterogeneous integration landscape, clonal expansion, and bottlenecks to eradicate viral persistence. Finally, we propose that building upon previous discoveries and improved or yet-to-be discovered technologies will unravel molecular mechanisms of latency establishment and reactivation in a “new era”.

MicroRNAs organize intrinsic variation into stem cell states

Understanding how mammalian organisms achieve the full diversity of cell types in the adult organism is a central goal of developmental cell biology. Recent work has shown that some embryonic precursor cells can self-organize into developmental structures but the mechanisms of gene regulation that contribute to this process remain unknown. Here we show embryonic stem cells self-organize into distinct gene expression states that resemble developmental gene programs. We find that microRNAs, small noncoding regulators of gene expression, play a critical role in organizing fluctuations across gene networks to help achieve this organization into distinct expression states.

Pluripotent embryonic stem cells (ESCs) contain the potential to form a diverse array of cells with distinct gene expression states, namely the cells of the adult vertebrate. Classically, diversity has been attributed to cells sensing their position with respect to external morphogen gradients. However, an alternative is that diversity arises in part from cooption of fluctuations in the gene regulatory network. Here we find ESCs exhibit intrinsic heterogeneity in the absence of external gradients by forming interconverting cell states. States vary in developmental gene expression programs and display distinct activity of microRNAs (miRNAs). Notably, miRNAs act on neighborhoods of pluripotency genes to increase variation of target genes and cell states. Loss of miRNAs that vary across states reduces target variation and delays state transitions, suggesting variable miRNAs organize and propagate variation to promote state transitions. Together these findings provide insight into how a gene regulatory network can coopt variation intrinsic to cell systems to form robust gene expression states. Interactions between intrinsic heterogeneity and environmental signals may help achieve developmental outcomes.

Differences in Transcriptional Dynamics Between T-cells and Macrophages as Determined by a Three-State Mathematical Model

HIV-1 viral transcription persists in patients despite antiretroviral treatment, potentially due to intermittent HIV-1 LTR activation. While several mathematical models have been explored in the context of LTR-protein interactions, in this work for the first time HIV-1 LTR model featuring repressed, intermediate, and activated LTR states is integrated with generation of long (env) and short (TAR) RNAs and proteins (Tat, Pr55, and p24) in T-cells and macrophages using both cell lines and infected primary cells. This type of extended modeling framework allows us to compare and contrast behavior of these two cell types. We demonstrate that they exhibit unique LTR dynamics, which ultimately results in differences in the magnitude of viral products generated. One of the distinctive features of this work is that it relies on experimental data in reaction rate computations. Two RNA transcription rates from the activated promoter states are fit by comparison of experimental data to model predictions. Fitting to the data also provides estimates for the degradation/exit rates for long and short viral RNA. Our experimentally generated data is in reasonable agreement for the T-cell as well macrophage population and gives strong evidence in support of using the proposed integrated modeling paradigm. Sensitivity analysis performed using Latin hypercube sampling method confirms robustness of the model with respect to small parameter perturbations. Finally, incorporation of a transcription inhibitor (F07#13) into the governing equations demonstrates how the model can be used to assess drug efficacy. Collectively, our model indicates transcriptional differences between latently HIV-1 infected T-cells and macrophages and provides a novel platform to study various transcriptional dynamics leading to latency or activation in numerous cell types and physiological conditions.

Mammalian gene expression variability is explained by underlying cell state

Gene expression variability in mammalian systems plays an important role in physiological and pathophysiological conditions. This variability can come from differential regulation related to cell state (extrinsic) and allele-specific transcriptional bursting (intrinsic). Yet, the relative contribution of these two distinct sources is unknown. Here, we exploit the qualitative difference in the patterns of covariance between these two sources to quantify their relative contributions to expression variance in mammalian cells. Using multiplexed error robust RNA fluorescent in situ hybridization (MERFISH), we measured the multivariate gene expression distribution of 150 genes related to Ca2+ signaling coupled with the dynamic Ca2+ response of live cells to ATP. We show that after controlling for cellular phenotypic states such as size, cell cycle stage, and Ca2+ response to ATP, the remaining variability is effectively at the Poisson limit for most genes. These findings demonstrate that the majority of expression variability results from cell state differences and that the contribution of transcriptional bursting is relatively minimal.

Cell Size-Based Decision-Making of a Viral Gene Circuit

Latently infected T cells able to reinitiate viral propagation throughout the body remain a major barrier to curing HIV. Distinguishing between latently infected cells and uninfected cells will advance efforts for viral eradication. HIV decision-making between latency and active replication is stochastic, and drug cocktails that increase bursts of viral gene expression enhance reactivation from latency. Here, we show that a larger host-cell size provides a natural cellular mechanism for enhancing burst size of viral expression and is necessary to destabilize the latent state and bias viral decision-making. Latently infected Jurkat and primary CD4+ T cells reactivate exclusively in larger activated cells, while smaller cells remain silent. In addition, reactivation is cell-cycle dependent and can be modulated with cell-cycle-arresting compounds. Cell size and cell-cycle dependent decision-making of viral circuits may guide stochastic design strategies and applications in synthetic biology and may provide important determinants to advance diagnostics and therapies.

Bohn-Wippert et al. investigate reactivation of T cells latently infected with HIV. They discover that only larger cells exit latency, while smaller cells remain silent. Viral expression bursts are cell size and cell-cycle dependent, presenting dynamic cell states, capable of active control, as sources of viral fate determination.

Stochastic transcription in the p53-mediated response to DNA damage is modulated by burst frequency

Discontinuous transcription has been described for different mammalian cell lines and numerous promoters. However, our knowledge of how the activity of individual promoters is adjusted by dynamic signaling inputs from transcription factors is limited. To address this question, we characterized the activity of selected target genes that are regulated by pulsatile accumulation of the tumor suppressor p53 in response to ionizing radiation. We performed time-resolved measurements of gene expression at the single-cell level by smFISH and used the resulting data to inform a mathematical model of promoter activity. We found that p53 target promoters are regulated by frequency modulation of stochastic bursting and can be grouped along three archetypes of gene expression. The occurrence of these archetypes cannot solely be explained by nuclear p53 abundance or promoter binding of total p53. Instead, we provide evidence that the time-varying acetylation state of p53's C-terminal lysine residues is critical for gene-specific regulation of stochastic bursting.

Stable integrant-specific differences in bimodal HIV-1 expression patterns revealed by high-throughput analysis

HIV-1 gene expression is regulated by host and viral factors that interact with viral motifs and is influenced by proviral integration sites. Here, expression variation among integrants was followed for hundreds of individual proviral clones within polyclonal populations throughout successive rounds of virus and cultured cell replication, with limited findings using CD4+ cells from donor blood consistent with observations in immortalized cells. Tracking clonal behavior by proviral “zip codes” indicated that mutational inactivation during reverse transcription was rare, while clonal expansion and proviral expression states varied widely. By sorting for provirus expression using a GFP reporter in the nef open reading frame, distinct clone-specific variation in on/off proportions were observed that spanned three orders of magnitude. Tracking GFP phenotypes over time revealed that as cells divided, their progeny alternated between HIV transcriptional activity and non-activity. Despite these phenotypic oscillations, the overall GFP+ population within each clone was remarkably stable, with clones maintaining clone-specific equilibrium mixtures of GFP+ and GFP- cells. Integration sites were analyzed for correlations between genomic features and the epigenetic phenomena described here. Integrants inserted in the sense orientation of genes were more frequently found to be GFP negative than those in the antisense orientation, and clones with high GFP+ proportions were more distal to repressive H3K9me3 peaks than low GFP+ clones. Clones with low frequencies of GFP positivity appeared to expand more rapidly than clones for which most cells were GFP+, even though the tested proviruses were Vpr-. Thus, much of the increase in the GFP- population in these polyclonal pools over time reflected differential clonal expansion. Together, these results underscore the temporal and quantitative variability in HIV-1 gene expression among proviral clones that are conferred in the absence of metabolic or cell-type dependent variability, and shed light on cell-intrinsic layers of regulation that affect HIV-1 population dynamics.

Very few HIV-1 infected cells persist in patients for more than a couple days, but those that do pose life-long health risks. Strategies designed to eliminate these cells have been based on assumptions about what viral properties allow infected cell survival. However, such approaches for HIV-1 eradication have not yet shown therapeutic promise, possibly because many assumptions about virus persistence are based on studies involving a limited number of infected cell types, the averaged behavior of cells in diverse populations, or snapshot views. Here, we developed a high-throughput approach to study hundreds of distinct HIV-1 infected cells and their progeny over time in an unbiased way. This revealed that each virus established its own pattern of gene expression that, upon infected cell division, was stably transmitted to all progeny cells. Expression patterns consisted of alternating waves of activity and inactivity, with the extent of activity differing among infected cell families over a 1000-fold range. The dynamics and variability among infected cells and within complex populations that the work here revealed has not previously been evident, and may help establish more accurate correlates of persistent HIV-1 infection.

Discontinuous transcription

Numerous studies based on new single-cell and single-gene techniques show that individual genes can be transcribed in short bursts or pulses accompanied by changes in pulsing frequencies. Since so many examples of such discontinuous or fluctuating transcription have been found from prokaryotes to mammals, it now seems to be a common mode of gene expression. In this review we discuss the occurrence of the transcriptional fluctuations, the techniques used for their detection, their putative causes, kinetic characteristics, and probable physiological significance.