Dan J. Woodcock

Dynamic Analysis of Stochastic Transcription Cycles

Cycling of gene expression in individual cells is controlled by dynamic chromatin remodeling.

Physiological levels of TNFα stimulation induce stochastic dynamics of NF-κB responses in single living cells

Nuclear factor kappa B (NF-κB) signalling is activated by cellular stress and inflammation and regulates cytokine expression. We applied single-cell imaging to investigate dynamic responses to different doses of tumour necrosis factor alpha (TNFα). Lower doses activated fewer cells and those responding showed an increasingly variable delay in the initial NF-κB nuclear translocation and associated IκBα degradation. Robust 100 minute nuclear:cytoplasmic NF-κB oscillations were observed over a wide range of TNFα concentrations. The result is supported by computational analyses, which identified a limit cycle in the system with a stable 100 minute period over a range of stimuli, and indicated no co-operativity in the pathway activation. These results suggest that a stochastic threshold controls functional all-or-nothing responses in individual cells. Deterministic and stochastic models simulated the experimentally observed activation threshold and gave rise to new predictions about the structure of the system and open the way for better mechanistic understanding of physiological TNFα activation of inflammatory responses in cells and tissues.

Quantifying intrinsic and extrinsic noise in gene transcription using the linear noise approximation: An application to single cell data

A central challenge in computational modeling of dynamic biological systems is parameter inference from experimental time course measurements. However, one would not only like to infer kinetic parameters but also study their variability from cell to cell. Here we focus on the case where single-cell fluorescent protein imaging time series data are available for a population of cells. Based on van Kampen?s linear noise approximation, we derive a dynamic state space model for molecular populations which is then extended to a hierarchical model. This model has potential to address the sources of variability relevant to single-cell data, namely, intrinsic noise due to the stochastic nature of the birth and death processes involved in reactions and extrinsic noise arising from the cell-to-cell variation of kinetic parameters. In order to infer such a model from experimental data, one must also quantify the measurement process where one has to allow for nonmeasurable molecular species as well as measurement noise of unknown level and variance. The availability of multiple single-cell time series data here provides a unique testbed to fit such a model and quantify these different sources of variation from experimental data.

Evidence for Phenotypic Plasticity among Multihost Campylobacter jejuni and C. coli Lineages, Obtained Using Ribosomal Multilocus Sequence Typing and Raman Spectroscopy

ABSTRACT Closely related bacterial isolates can display divergent phenotypes. This can limit the usefulness of phylogenetic studies for understanding bacterial ecology and evolution. Here, we compare phenotyping based on Raman spectrometric analysis of cellular composition to phylogenetic classification by ribosomal multilocus sequence typing (rMLST) in 108 isolates of the zoonotic pathogens Campylobacter jejuni and C. coli. Automatic relevance determination (ARD) was used to identify informative peaks in the Raman spectra that could be used to distinguish strains in taxonomic and host source groups (species, clade, clonal complex, and isolate source/host). Phenotypic characterization based on Raman spectra showed a degree of agreement with genotypic classification using rMLST, with segregation accuracy between species (83.95%), clade (in C. coli, 98.41%), and, to some extent, clonal complex (86.89% C. jejuni ST-21 and ST-45 complexes) being achieved. This confirmed the utility of Raman spectroscopy for lineage classification and the correlation between genotypic and phenotypic classification. In parallel analysis, relatively distantly related isolates (different clonal complexes) were assigned the correct host origin irrespective of the clonal origin (74.07 to 96.97% accuracy) based upon different Raman peaks. This suggests that the phenotypic characteristics, from which the phenotypic signal is derived, are not fixed by clonal descent but are influenced by the host environment and change as strains move between hosts.

A hierarchical model of transcriptional dynamics allows robust estimation of transcription rates in populations of single cells with variable gene copy number

Motivation: cis-regulatory DNA sequence elements, such as enhancers and silencers, function to control the spatial and temporal expression of their target genes. Although the overall levels of gene expression in large cell populations seem to be precisely controlled, transcription of individual genes in single cells is extremely variable in real time. It is, therefore, important to understand how these cis-regulatory elements function to dynamically control transcription at single-cell resolution. Recently, statistical methods have been proposed to back calculate the rates involved in mRNA transcription using parameter estimation of a mathematical model of transcription and translation. However, a major complication in these approaches is that some of the parameters, particularly those corresponding to the gene copy number and transcription rate, cannot be distinguished; therefore, these methods cannot be used when the copy number is unknown. Results: Here, we develop a hierarchical Bayesian model to estimate biokinetic parameters from live cell enhancer–promoter reporter measurements performed on a population of single cells. This allows us to investigate transcriptional dynamics when the copy number is variable across the population. We validate our method using synthetic data and then apply it to quantify the function of two known developmental enhancers in real time and in single cells. Availability: Supporting information is submitted with the article. Contact: d.j.woodcock@warwick.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Genomic plasticity and rapid host switching can promote the evolution of generalism: a case study in the zoonotic pathogen Campylobacter

Horizontal gene transfer accelerates bacterial adaptation to novel environments, allowing selection to act on genes that have evolved in multiple genetic backgrounds. This can lead to ecological specialization. However, little is known about how zoonotic bacteria maintain the ability to colonize multiple hosts whilst competing with specialists in the same niche. Here we develop a stochastic evolutionary model and show how genetic transfer of host segregating alleles, distributed as predicted for niche specifying genes, and the opportunity for host transition could interact to promote the emergence of host generalist lineages of the zoonotic bacterium Campylobacter. Using a modelling approach we show that increasing levels of homologous recombination enhance the efficiency with which selection can fix combinations of beneficial alleles, speeding adaptation. We then show how these predictions change in a multi-host system, with low levels of recombination, consistent with real r/m estimates, increasing the standing variation in the population, allowing a more effective response to changes in the selective landscape. Our analysis explains how observed gradients of host specialism and generalism can evolve in a multihost system through the transfer of ecologically important loci among coexisting strains.

Lineage-specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co‐evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large‐scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co‐evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.

The long non‐coding RNA Paupar promotes KAP1‐dependent chromatin changes and regulates olfactory bulb neurogenesis

Many long non‐coding RNAs (lncRNAs) are expressed during central nervous system (CNS) development, yet their in vivo roles and mechanisms of action remain poorly understood. Paupar, a CNS‐expressed lncRNA, controls neuroblastoma cell growth by binding and modulating the activity of transcriptional regulatory elements in a genome‐wide manner. We show here that the Paupar lncRNA directly binds KAP1, an essential epigenetic regulatory protein, and thereby regulates the expression of shared target genes important for proliferation and neuronal differentiation. Paupar promotes KAP1 chromatin occupancy and H3K9me3 deposition at a subset of distal targets, through the formation of a ribonucleoprotein complex containing Paupar, KAP1 and the PAX6 transcription factor. Paupar‐KAP1 genome‐wide co‐occupancy reveals a fourfold enrichment of overlap between Paupar and KAP1 bound sequences, the majority of which also appear to associate with PAX6. Furthermore, both Paupar and Kap1 loss‐of‐function in vivo disrupt olfactory bulb neurogenesis. These observations provide important conceptual insights into the trans‐acting modes of lncRNA‐mediated epigenetic regulation and the mechanisms of KAP1 genomic recruitment, and identify Paupar and Kap1 as regulators of neurogenesis in vivo.

Temperature regulates NF-κB dynamics and function through timing of A20 transcription

Significance Inflammation is often accompanied by temperature change, but little is known about the role of temperature in the inflammatory response. We show that physiologically relevant temperature changes significantly perturb NF-κB dynamics following TNFα stimulation in single cells. Using experimentation informed by mathematical modeling, we found that these changes were mediated, at least in part, through the key feedback gene TNFAIP3/A20. Curtailing A20 expression removed temperature sensitivity across the fever range (37 °C to 40 °C). Gene expression was generally unaffected between these temperatures, although a select set of NF-κB−regulated genes was up-regulated at early time points. These genes were predominantly involved in inflammation, signaling, and cell fate. The cellular response to inflammation may therefore be mechanistically and functionally regulated by temperature. NF-κB signaling plays a pivotal role in control of the inflammatory response. We investigated how the dynamics and function of NF-κB were affected by temperature within the mammalian physiological range (34 °C to 40 °C). An increase in temperature led to an increase in NF-κB nuclear/cytoplasmic oscillation frequency following Tumor Necrosis Factor alpha (TNFα) stimulation. Mathematical modeling suggested that this temperature sensitivity might be due to an A20-dependent mechanism, and A20 silencing removed the sensitivity to increased temperature. The timing of the early response of a key set of NF-κB target genes showed strong temperature dependence. The cytokine-induced expression of many (but not all) later genes was insensitive to temperature change (suggesting that they might be functionally temperature-compensated). Moreover, a set of temperature- and TNFα-regulated genes were implicated in NF-κB cross-talk with key cell-fate–controlling pathways. In conclusion, NF-κB dynamics and target gene expression are modulated by temperature and can accurately transmit multidimensional information to control inflammation.

Paupar LncRNA Promotes KAP1 Dependent Chromatin Changes And Regulates Subventricular Zone Neurogenesis

Many long non-coding RNAs (lncRNAs) are expressed during central nervous system (CNS) development, yet their in vivo roles and molecular mechanisms of action remain poorly understood. Paupar, a CNS expressed lncRNA, controls neuroblastoma cell growth by binding and modulating the activity of genome-wide transcriptional regulatory elements. We show here that Paupar transcript directly binds KAP1, an essential epigenetic regulatory protein, and thereby regulates the expression of shared target genes important for proliferation and neuronal differentiation. Paupar promotes KAP1 chromatin occupancy and H3K9me3 deposition at a subset of distal targets, through formation of a DNA binding ribonucleoprotein complex containing Paupar, KAP1 and the PAX6 transcription factor. Paupar-KAP1 genome-wide co-occupancy reveals a 4-fold enrichment of overlap between Paupar and KAP1 bound sequences. Furthermore, both Paupar and Kap1 loss of function in vivo accelerates lineage progression in the mouse postnatal subventricular zone (SVZ) stem cell niche and disrupts olfactory bulb neurogenesis. These observations provide important conceptual insights into the trans-acting modes of lncRNA-mediated epigenetic regulation, the mechanisms of KAP1 genomic recruitment and identify Paupar and Kap1 as regulators of SVZ neurogenesis.