No abstract

Stochasticity in gene expression can result in fluctuations in gene product levels. Recent experiments indicated that feedback regulation plays an important role in controlling the noise in gene expression. A quantitative understanding of the feedback effect on gene expression requires analysis of the corresponding stochastic model. However, for stochastic models of gene expression with general regulation functions, exact analytical results for gene product distributions have not been given so far. Here, we propose a technique to solve a generalized ON-OFF model of stochastic gene expression with arbitrary (positive or negative, linear or nonlinear) feedbacks including posttranscriptional or posttranslational regulation. The obtained results, which generalize results obtained previously, provide new insights into the role of feedback in regulating gene expression. The proposed analytical framework can easily be extended to analysis of more complex models of stochastic gene expression.

<abstract> <p>While transcription often occurs in a bursty manner, various possible regulations can lead to complex promoter patterns such as promoter cycles, giving rise to an important question: How do promoter kinetics shape transcriptional bursting kinetics? Here we introduce and analyze a general model of the promoter cycle consisting of multi-OFF states and multi-ON states, focusing on the effects of multi-ON mechanisms on transcriptional bursting kinetics. The derived analytical results indicate that burst size follows a mixed geometric distribution rather than a single geometric distribution assumed in previous studies, and ON and OFF times obey their own mixed exponential distributions. In addition, we find that the multi-ON mechanism can lead to bimodal burst-size distribution, antagonistic timing of ON and OFF, and diverse burst frequencies, each further contributing to cell-to-cell variability in the mRNA expression level. These results not only reveal essential features of transcriptional bursting kinetics patterns shaped by multi-state mechanisms but also can be used to the inferences of transcriptional bursting kinetics and promoter structure based on experimental data.</p> </abstract>

Bimodality of gene expression, as a mechanism generating phenotypic diversity, enhances the survival of cells in a fluctuating environment. Previous studies have shown that the bimodality of gene expression is...

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.