Predicting stimulation-dependent enhancer-promoter interactions from ChIP-Seq time course data.

We have developed a machine learning approach to predict stimulation-dependent enhancer-promoter interactions using evidence from changes in genomic protein occupancy over time. The occupancy of estrogen receptor alpha (ERα), RNA polymerase (Pol II) and histone marks H2AZ and H3K4me3 were measured over time using ChIP-Seq experiments in MCF7 cells stimulated with estrogen. A Bayesian classifier was developed which uses the correlation of temporal binding patterns at enhancers and promoters and genomic proximity as features to predict interactions. This method was trained using experimentally determined interactions from the same system and was shown to achieve much higher precision than predictions based on the genomic proximity of nearest ERα binding. We use the method to identify a genome-wide confident set of ERα target genes and their regulatory enhancers genome-wide. Validation with publicly available GRO-Seq data demonstrates that our predicted targets are much more likely to show early nascent transcription than predictions based on genomic ERα binding proximity alone.

A transient ischemic environment induces reversible compaction of chromatin.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Genome-wide modeling of transcription kinetics reveals patterns of RNA production delays.

Genes with similar transcriptional activation kinetics can display very different temporal mRNA profiles because of differences in transcription time, degradation rate, and RNA-processing kinetics. Recent studies have shown that a splicing-associated RNA production delay can be significant. To investigate this issue more generally, it is useful to develop methods applicable to genome-wide datasets. We introduce a joint model of transcriptional activation and mRNA accumulation that can be used for inference of transcription rate, RNA production delay, and degradation rate given data from high-throughput sequencing time course experiments. We combine a mechanistic differential equation model with a nonparametric statistical modeling approach allowing us to capture a broad range of activation kinetics, and we use Bayesian parameter estimation to quantify the uncertainty in estimates of the kinetic parameters. We apply the model to data from estrogen receptor α activation in the MCF-7 breast cancer cell line. We use RNA polymerase II ChIP-Seq time course data to characterize transcriptional activation and mRNA-Seq time course data to quantify mature transcripts. We find that 11% of genes with a good signal in the data display a delay of more than 20 min between completing transcription and mature mRNA production. The genes displaying these long delays are significantly more likely to be short. We also find a statistical association between high delay and late intron retention in pre-mRNA data, indicating significant splicing-associated production delays in many genes.

Replicative senescence is associated with nuclear reorganization and with DNA methylation at specific transcription factor binding sites.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Integrative analysis of deep sequencing data identifies estrogen receptor early response genes and links ATAD3B to poor survival in breast cancer.

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 α (ERα), H3K4me3 and H2A.Z occupancy at five time points in the MCF-7 breast cancer cell line after estradiol stimulus. We obtained 777 ERa early responsive genes and compared the biological functions of the genes having ERα 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 ERα via the MAPK pathway, instead of direct ERa binding, may be responsible for early cell responses to ERα activation. Our results also indicate that the ERα responsive genes triggered by the genomic pathway are transcribed faster than those without ERα binding sites. The survival analysis of the 777 ERα responsive genes with 150 primary breast cancer tumors and in two independent validation cohorts indicated the ATAD3B gene, which does not have ERα binding site within 20 kb of its TSS, to be significantly associated with poor patient survival.

The nonspecific lethal complex is a transcriptional regulator in Drosophila.

Here, we report the biochemical characterization of the nonspecific lethal (NSL) complex (NSL1, NSL2, NSL3, MCRS2, MBD-R2, and WDS) that associates with the histone acetyltransferase MOF in both Drosophila and mammals. Chromatin immunoprecipitation-Seq analysis revealed association of NSL1 and MCRS2 with the promoter regions of more than 4000 target genes, 70% of these being actively transcribed. This binding is functional, as depletion of MCRS2, MBD-R2, and NSL3 severely affects gene expression genome wide. The NSL complex members bind to their target promoters independently of MOF. However, depletion of MCRS2 affects MOF recruitment to promoters. NSL complex stability is interdependent and relies mainly on the presence of NSL1 and MCRS2. Tethering of NSL3 to a heterologous promoter leads to robust transcription activation and is sensitive to the levels of NSL1, MCRS2, and MOF. Taken together, we conclude that the NSL complex acts as a major transcriptional regulator in Drosophila.

Contrasting effects of Elg1-RFC and Ctf18-RFC inactivation in the absence of fully functional RFC in fission yeast.

Proliferating cell nuclear antigen loading onto DNA by replication factor C (RFC) is a key step in eukaryotic DNA replication and repair processes. In this study, the C-terminal domain (CTD) of the large subunit of fission yeast RFC is shown to be essential for its function in vivo. Cells carrying a temperature-sensitive mutation in the CTD, rfc1-44, arrest with incompletely replicated chromosomes, are sensitive to DNA damaging agents, are synthetically lethal with other DNA replication mutants, and can be suppressed by mutations in rfc5. To assess the contribution of the RFC-like complexes Elg1-RFC and Ctf18-RFC to the viability of rfc1-44, genes encoding the large subunits of these complexes have been deleted and overexpressed. Inactivation of Ctf18-RFC by the deletion of ctf18+, dcc1+ or ctf8+ is lethal in an rfc1-44 background showing that full Ctf18-RFC function is required in the absence of fully functional RFC. In contrast, rfc1-44 elg1Delta cells are viable and overproduction of Elg1 in rfc1-44 is lethal, suggesting that Elg1-RFC plays a negative role when RFC function is inhibited. Consistent with this, the deletion of elg1+ is shown to restore viability to rfc1-44 ctf18Delta cells.