Stationary moments, distribution conjugation and phenotypic regions in stochastic gene transcription

Transcription is a pivotal step in gene expression yet a complex biochemical process. It occurs often in a bursty manner, leading to cell-to-cell variability important for the cells surviving in complex environments. Quantitative experiments of such stochastic transcription call for model analysis in a systematic rather than case-by-case fashion. Here we analyze two general yet biologically-reasonable classes of stochastic transcription models: the first class called promoter models considers that promoter structure is general, and the second class called queuing models considers that waiting-time distributions are general. For the former, we show that there are conjugate relationships between models with different transcription exits. This property well reveals the mechanic principle of stochastic transcription in complex cases of transcription factor regulation. For the latter, we establish an integral equation for the mRNA moment-generating function, through which stationary mRNA moments of any orders can be analytically derived. Finally, we analyze parametric regions for robust unimodality and bimodality. The overall analysis not only lays a foundation for quantitative analysis of stochastic transcription but also reveals how the upstream promoter kinetics impact the downstream expression dynamics.