[Paper Review] Abortive Initiation as a Bottleneck for Transcription in the Early Drosophila Embryo

Alexander S. Serov, Alexander Levine+1 more|arXiv (Cornell University)|Jan 21, 2017

Genomics and Chromatin Dynamics33 references3 citations

TL;DR

This study identifies abortive transcription initiation as the primary bottleneck limiting transcription rates in early Drosophila embryos, explaining why observed maximal transcription rates (40% of theoretical maximum) are lower than predicted by the standard k-TASEP model. By incorporating experimentally measured abortive initiation kinetics into a modified k-TASEP model with extended particles, the authors achieve quantitative agreement with live MS2 imaging data, shifting the transcription bottleneck from the gene body to the promoter region and implying that synonymous codon optimization in the gene body cannot enhance transcription rates.

ABSTRACT

Gene transcription is a critical step in gene expression. The currently accepted physical model of transcription predicts the existence of a physical limit on the maximal rate of transcription, which does not depend on the transcribed gene. This limit appears as a result of polymerase "traffic jams" forming in the bulk of the 1D DNA chain at high polymerase concentrations. Recent experiments have, for the first time, allowed one to access live gene expression dynamics in the Drosophila fly embryo in vivo under the conditions of heavy polymerase load and test the predictions of the model. Our analysis of the data shows that the maximal rate of transcription is indeed the same for the Hunchback, Snail and Knirps gap genes, and modified gene constructs in nuclear cycles 13, 14, but the experimentally observed value of the maximal transcription rate corresponds to only 40 % of the one predicted by this model. We argue that such a decrease must be due to a slower polymerase elongation rate in the vicinity of the promoter region. This effectively shifts the bottleneck of transcription from the bulk to the promoter region of the gene. We suggest a quantitative explanation of the difference by taking into account abortive transcription initiation. Our calculations based on the independently measured abortive initiation constant in vitro confirm this hypothesis and find quantitative agreement with MS2 fluorescence live imaging data in the early fruit fly embryo. If our explanation is correct, then the transcription rate cannot be increased by replacing "slow codons" in the bulk with synonymous codons, and experimental efforts must be focused on the promoter region instead. This study extends our understanding of transcriptional regulation, re-examines physical constraints on the kinetics of transcription and re-evaluates the validity of the standard TASEP model of transcription.

Motivation & Objective

• To determine whether physical constraints on transcription limit developmental gene expression in early Drosophila embryos.
• To test whether the standard k-TASEP model accurately predicts maximal transcription rates under high polymerase load.
• To investigate why experimentally observed transcription rates are only 40% of the theoretical maximum predicted by the MacDonald-Gibbs-Pipkin model.
• To identify the true location of the transcription bottleneck—whether in the gene body or near the promoter—using live imaging and theoretical modeling.
• To evaluate the impact of abortive initiation on transcriptional kinetics and assess its implications for experimental strategies to enhance gene expression.

Proposed method

• Adaptation of the totally asymmetric simple exclusion process (TASEP) model to include extended polymerase footprints and inhomogeneous elongation rates.
• Incorporation of a slow, rate-limiting first site (promoter region) to model abortive initiation, with the initiation time constant τ derived from in vitro measurements.
• Use of a modified k-TASEP model with quenched disorder to simulate polymerase traffic under high-coverage conditions in nuclear cycles 13 and 14.
• Comparison of simulated transcription dynamics (polymerase density, mRNA production rate) with live MS2 fluorescence imaging data from Drosophila embryos.
• Analysis of non-steady-state dynamics and distribution widths of steady-state polymerase numbers to detect noise sources and validate model assumptions.
• Evaluation of whether the system reaches a true steady-state by examining temporal profiles of transcription rates and polymerase density.

Experimental results

Research questions

• RQ1 Does the observed transcription rate in early Drosophila embryos saturate the physical limit predicted by the k-TASEP model under high polymerase load?
• RQ2 Why is the experimentally observed maximal transcription rate only 40% of the theoretical maximum predicted by the standard k-TASEP model?
• RQ3 Is the transcription bottleneck located in the gene body due to polymerase traffic jams, or is it shifted to the promoter region due to abortive initiation?
• RQ4 Can the discrepancy between theory and experiment be quantitatively explained by incorporating abortive initiation kinetics into the k-TASEP framework?
• RQ5 What are the implications of a promoter-proximal bottleneck for strategies aimed at increasing transcriptional output?

Key findings

• The experimentally observed maximal transcription rate for Hunchback, Snail, and Knirps genes is consistently 40% of the theoretical maximum predicted by the standard k-TASEP model.
• The discrepancy between theory and experiment is quantitatively explained by incorporating abortive initiation with an experimentally measured initiation time constant τ into the k-TASEP model.
• The model with a slow first site (promoter region) and extended polymerase footprints reproduces the live MS2 imaging data with high fidelity, confirming the promoter as the primary bottleneck.
• The model predicts that synonymous codon replacement in the gene body cannot increase transcription rates, as the bottleneck is not in elongation but in initiation.
• Non-steady-state dynamics, such as gradual polymerase density decline after peak production, cannot be explained by the standard or abortive k-TASEP models and require further investigation.
• The wider-than-expected experimental distributions of steady-state polymerase numbers suggest additional noise sources beyond those modeled, indicating a need for further noise analysis.