Stochastic holin expression can account for lysis time variation in the bacteriophage λ

The inherent stochastic nature of biochemical processes can drive differences in gene expression between otherwise identical cells. While cell-to-cell variability in gene expression has received much attention, randomness in timing of events has been less studied. We investigate event timing at the single-cell level in a simple system, the lytic pathway of the bacterial virus phage lambda. In individual cells, lysis occurs on average at 65 min, with an s.d. of 3.5 min. Interestingly, mutations in the lysis protein, holin, alter both the lysis time (LT) mean and variance. In our analysis, LT is formulated as the first-passage time (FPT) for cellular holin levels to cross a critical threshold. Exact analytical formulae for the FPT moments are derived for stochastic gene expression models. These formulae reveal how holin transcription and translation efficiencies independently modulate the LT mean and variation. Analytical expressions for the LT moments are used to evaluate previously published single-cell LT data for lambda phages with mutations in the holin sequence or its promoter. Our results show that stochastic holin expression is sufficient to account for the intercellular LT differences in both wild-type phages, and phage variants where holin transcription and the threshold for lysis have been experimentally altered. Finally, our analysis reveals regulatory motifs that enhance the robustness of lysis timing to cellular noise.