  A fundamental question of cell biology is how cells control the number of
organelles. The processes of organelle biogenesis, namely de novo synthesis,
fission, fusion, and decay, are inherently stochastic, producing cell-to-cell
variability in organelle abundance. In addition, experiments suggest that the
synthesis of some organelles can be bursty. We thus ask how bursty synthesis
impacts intracellular organelle number distribution. We develop an organelle
biogenesis model with bursty de novo synthesis by considering geometrically
distributed burst sizes. We analytically solve the model in biologically
relevant limits and provide exact expressions for the steady-state organelle
number distributions and their means and variances. We also present approximate
solutions for the whole model, complementing with exact stochastic simulations.
We show that bursts generally increase the noise in organelle numbers,
producing distinct signatures in noise profiles depending on different
mechanisms of organelle biogenesis. We also find different shapes of organelle
number distributions, including bimodal distributions in some parameter
regimes. Notably, bursty synthesis broadens the parameter regime of observing
bimodality compared to the `non-bursty' case. Together, our framework utilizes
number fluctuations to elucidate the role of bursty synthesis in producing
organelle number heterogeneity in cells.
