George Reid

Cellular decision‐making and environmental adaptation are dependent upon a heterogeneous response of gene expression to external cues. Heterogeneity arises in transcription from random switching between transcriptionally active and inactive states, resulting in bursts of RNA synthesis. Furthermore, the cellular state influences the competency of transcription, thereby globally affecting gene expression in a cell‐specific manner. We determined how external stimuli interplay with cellular state to modulate the kinetics of bursting. To this end, single‐cell dynamics of nascent transcripts were monitored at the endogenous estrogen‐responsive GREB1 locus. Stochastic modeling of gene expression implicated a two‐state promoter model in which the estrogen stimulus modulates the frequency of transcriptional bursting. The cellular state affects transcriptional dynamics by altering initiation and elongation kinetics and acts globally, as GREB1 alleles in the same cell correlate in their transcriptional output. Our results suggest that cellular state strongly affects the first step of the central dogma of gene expression, to promote heterogeneity in the transcriptional output of isogenic cells.

Cells detect and adapt to hypoxic and nutritional stress through immediate transcriptional, translational and metabolic responses. The environmental effects of ischemia on chromatin nanostructure were investigated using single molecule localization microscopy of DNA binding dyes and of acetylated histones, by the sensitivity of chromatin to digestion with DNAseI, and by fluorescence recovery after photobleaching (FRAP) of core and linker histones. Short-term oxygen and nutrient deprivation of the cardiomyocyte cell line HL-1 induces a previously undescribed chromatin architecture, consisting of large, chromatin-sparse voids interspersed between DNA-dense hollow helicoid structures 40–700 nm in dimension. The chromatin compaction is reversible, and upon restitution of normoxia and nutrients, chromatin transiently adopts a more open structure than in untreated cells. The compacted state of chromatin reduces transcription, while the open chromatin structure induced upon recovery provokes a transitory increase in transcription. Digestion of chromatin with DNAseI confirms that oxygen and nutrient deprivation induces compaction of chromatin. Chromatin compaction is associated with depletion of ATP and redistribution of the polyamine pool into the nucleus. FRAP demonstrates that core histones are not displaced from compacted chromatin; however, the mobility of linker histone H1 is considerably reduced, to an extent that far exceeds the difference in histone H1 mobility between heterochromatin and euchromatin. These studies exemplify the dynamic capacity of chromatin architecture to physically respond to environmental conditions, directly link cellular energy status to chromatin compaction and provide insight into the effect ischemia has on the nuclear architecture of cells.

Significance: Gene transcription is a highly regulated dynamic process. Delays in transcription have important consequences on dynamics of gene expression and consequently on downstream biological function. We model temporal dynamics of transcription using genome-wide time course data measuring transcriptional activity and mRNA concentration. We find a significant number of genes exhibit a long RNA processing delay between transcription termination and mRNA production. These long processing delays are more common for short genes, which would otherwise be expected to transcribe most rapidly. The distribution of intronic reads suggests that these delays are required for splicing to be completed. Understanding such delays is essential for understanding how a rapid cellular response is regulated.

The environmental effects of ischemia on chromatin nanostructure were evaluated using single molecule localisation microscopy (SMLM) of DNA binding dyes. Short-term oxygen and nutrient deprivation (OND) of the cardiomyocyte cell-line HL-1induces a previously undescribed chromatin architecture, consisting of large, chromatin sparse voids interspersed between DNA-dense hollow helicoid structures of the order of 40 to 700 nm in dimension. OND induced chromatin compaction is reversible, and upon restitution of normoxia and nutrients, chromatin transiently adopts a significantly more open structure than in untreated cells. We show that this compacted state of chromatin reduces transcription, while the open chromatin structure following recovery has a higher transcriptional rate than in untreated cells. Digestion of chromatin with DNAseI and DNA binding dye loading assays confirm that OND induces compaction of chromatin and a general redistribution of chromatin to the nuclear periphery. Mechanistically, chromatin compaction is associated with a depletion of intracellular ATP and a redistribution of the cellular polyamine pool into the nucleus. Additionally, Fluorescence Recovery After Photobleaching (FRAP) shows that core histones are not displaced from compacted chromatin and that the mobility of linker histone H1 is considerably reduced by OND treatment, to an extent that far exceeds the difference in histone H1 mobility between heterochromatin and euchromatin. These studies exemplify the dynamic capacity of chromatin architecture to physically respond to environmental conditions, directly link cellular energy status to chromatin compaction and provide insight into the effect ischemia has on the nuclear architecture of cells.

Primary cells enter replicative senescence after a limited number of cell divisions. This process needs to be considered in cell culture experiments, and it is particularly important for regenerative medicine. Replicative senescence is associated with reproducible changes in DNA methylation (DNAm) at specific sites in the genome. The mechanism that drives senescence-associated DNAm changes remains unknown - it may involve stochastic DNAm drift due to imperfect maintenance of epigenetic marks or it is directly regulated at specific sites in the genome. In this study, we analyzed the reorganization of nuclear architecture and DNAm changes during long-term culture of human fibroblasts and mesenchymal stromal cells (MSCs). We demonstrate that telomeres shorten and shift towards the nuclear center at later passages. In addition, DNAm profiles, either analyzed by MethylCap-seq or by 450k IlluminaBeadChip technology, revealed consistent senescence-associated hypermethylation in regions associated with H3K27me3, H3K4me3, and H3K4me1 histone marks, whereas hypomethylation was associated with chromatin containing H3K9me3 and lamina-associated domains (LADs). DNA hypermethylation was significantly enriched in the vicinity of genes that are either up- or downregulated at later passages. Furthermore, specific transcription factor binding motifs (e.g. EGR1, TFAP2A, and ETS1) were significantly enriched in differentially methylated regions and in the promoters of differentially expressed genes. Senescence-associated DNA hypermethylation occurs at specific sites in the genome and reflects functional changes in the course of replicative senescence. These results indicate that tightly regulated epigenetic modifications during long-term culture contribute to changes in nuclear organization and gene expression.

Gene transcription mediated by RNA polymerase II (pol-II) is a key step in gene expression. The dynamics of pol-II moving along the transcribed region influence the rate and timing of gene expression. In this work, we present a probabilistic model of transcription dynamics which is fitted to pol-II occupancy time course data measured using ChIP-Seq. The model can be used to estimate transcription speed and to infer the temporal pol-II activity profile at the gene promoter. Model parameters are estimated using either maximum likelihood estimation or via Bayesian inference using Markov chain Monte Carlo sampling. The Bayesian approach provides confidence intervals for parameter estimates and allows the use of priors that capture domain knowledge, e.g. the expected range of transcription speeds, based on previous experiments. The model describes the movement of pol-II down the gene body and can be used to identify the time of induction for transcriptionally engaged genes. By clustering the inferred promoter activity time profiles, we are able to determine which genes respond quickly to stimuli and group genes that share activity profiles and may therefore be co-regulated. We apply our methodology to biological data obtained using ChIP-seq to measure pol-II occupancy genome-wide when MCF-7 human breast cancer cells are treated with estradiol (E2). The transcription speeds we obtain agree with those obtained previously for smaller numbers of genes with the advantage that our approach can be applied genome-wide. We validate the biological significance of the pol-II promoter activity clusters by investigating cluster-specific transcription factor binding patterns and determining canonical pathway enrichment. We find that rapidly induced genes are enriched for both estrogen receptor alpha (ER) and FOXA1 binding in their proximal promoter regions. Cells express proteins in response to changes in their environment so as to maintain normal function. An initial step in the expression of proteins is transcription, which is mediated by RNA polymerase II (pol-II). To understand changes in transcription arising due to stimuli it is useful to model the dynamics of transcription. We present a probabilistic model of pol-II transcription dynamics that can be used to compute RNA transcription speed and infer the temporal pol-II activity at the gene promoter. The inferred promoter activity profile is used to determine genes that are responding in a coordinated manner to stimuli and are therefore potentially co-regulated. Model parameters are inferred using data from high-throughput sequencing assays, such as ChIP-Seq and GRO-Seq, and can therefore be applied genome-wide in an unbiased manner. We apply the method to pol-II ChIP-Seq time course data from breast cancer cells stimulated by estradiol in order to uncover the dynamics of early response genes in this system.

Identification of responsive genes to an extra-cellular cue enables characterization of pathophysiologically crucial biological processes. Deep sequencing technologies provide a powerful means to identify responsive genes, which creates a need for computational methods able to analyze dynamic and multi-level deep sequencing data. To answer this need we introduce here a data-driven algorithm, SPINLONG, which is designed to search for genes that match the user-defined hypotheses or models. SPINLONG is applicable to various experimental setups measuring several molecular markers in parallel. To demonstrate the SPINLONG approach, we analyzed ChIP-seq data reporting PolII, estrogen receptor ( ), H3K4me3 and H2A.Z occupancy at five time points in the MCF-7 breast cancer cell line after estradiol stimulus. We obtained 777 early responsive genes and compared the biological functions of the genes having binding within 20 kb of the transcription start site (TSS) to genes without such binding site. Our results show that the non-genomic action of via the MAPK pathway, instead of direct binding, may be responsible for early cell responses to activation. Our results also indicate that the responsive genes triggered by the genomic pathway are transcribed faster than those without binding sites. The survival analysis of the 777 responsive genes with 150 primary breast cancer tumors and in two independent validation cohorts indicated the ATAD3B gene, which does not have binding site within 20 kb of its TSS, to be significantly associated with poor patient survival. Cellular processes in mammalian cells are tightly regulated to ensure that the cells function properly as a part of an organism. Dysregulation of some of these processes, such as apoptosis, cell proliferation and growth, can lead to cancer. One of the most important regulation mechanisms for cellular processes is via activation of membrane receptors by extra-cellular stimulus. Such cues trigger signal cascades that lead to altered expression of a number of genes in the cell nucleus; a key challenge in biomedicine is to identify which genes respond to a specific stimulus. These so called response genes can be investigated on a whole-genome scale with genomic sequencing, which is a technology that can quantify protein binding to DNA or gene activation. Analysis of such whole-genome data, however, is challenging due to billions of data points measured in the experiments. Here we introduce a novel computational method, SPINLONG, which is a widely applicable novel computational method that integrates multiple levels of deep sequencing data to produce experimentally testable hypotheses. We applied SPINLONG to breast cancer data and found early responsive genes for estrogen receptor and analyzed their regulation. These analyses resulted in a gene whose high activity is associated with decreased breast cancer patient survival.

Estrogen receptor alpha (ERα) is a ligand dependent transcription factor that regulates the expression of target genes through interacting with cis‐acting estrogen response elements (EREs). However, only a minority of ERα binding sites are located within the proximal promoter regions of responsive genes. Here we report the characterization of an ERE located 9kbp upstream of the TSS of the cathepsin D gene (CTSD) that up‐regulates CTSD expression upon estrogen stimulation in MCF‐7 cells. Using ChIP, we show recruitment of ERα and phosphorylated PolII at the CTSD distal enhancer region. Moreover, we determine the kinetics of transient CpG methylation on the promoter region of CTSD and for the first time, at a distal enhancer element. We show that ERα is crucial for long‐distance regulation of CTSD expression involving a looping mechanism.