RNA-seq is now the technology of choice for genome-wide differential gene expression experiments, but it is not clear how many biological replicates are needed to ensure valid biological interpretation of the results or which statistical tools are best for analyzing the data. An RNA-seq experiment with 48 biological replicates in each of two conditions was performed to answer these questions and provide guidelines for experimental design. With three biological replicates, nine of the 11 tools evaluated found only 20%-40% of the significantly differentially expressed (SDE) genes identified with the full set of 42 clean replicates. This rises to >85% for the subset of SDE genes changing in expression by more than fourfold. To achieve >85% for all SDE genes regardless of fold change requires more than 20 biological replicates. The same nine tools successfully control their false discovery rate at ≲5% for all numbers of replicates, while the remaining two tools fail to control their FDR adequately, particularly for low numbers of replicates. For future RNA-seq experiments, these results suggest that at least six biological replicates should be used, rising to at least 12 when it is important to identify SDE genes for all fold changes. If fewer than 12 replicates are used, a superior combination of true positive and false positive performances makes edgeR and DESeq2 the leading tools. For higher replicate numbers, minimizing false positives is more important and DESeq marginally outperforms the other tools.

The positioning of nucleosomes within the coding regions of eukaryotic genes is aligned with respect to transcriptional start sites. This organization is likely to influence many genetic processes, requiring access to the underlying DNA. Here we show that the combined action of Isw1 and Chd1 nucleosome spacing enzymes is required to maintain this organization. In the absence of these enzymes regular positioning of the majority of nucleosomes is lost. Exceptions include the region upstream of the promoter, the +1 nucleosome and a subset of locations distributed throughout coding regions where other factors are likely to be involved. These observations indicated that ATP-dependent remodeling enzymes are responsible for directing the positioning of the majority of nucleosomes within the Saccharomyces cerevisiae genome.

We present a genomewide, high-resolution study of nucleolar-associated chromatin using comparative genome hybridization, deep sequencing, and photoactivation microscopy. We show specific regions from most chromosomes associate with nucleoli. These regions are AT-rich, enriched in repressed genes, and overlap with nuclear lamina-associated loci.

Oxygen is essential for the life of most multicellular organisms. Cells possess enzymes called molecular dioxygenases that depend on oxygen for activity. A subclass of molecular dioxygenases is the histone demethylase enzymes, which are characterized by the presence of a Jumanji-C (JmjC) domain. Hypoxia can alter chromatin, but whether this is a direct effect on JmjC-histone demethylases or due to other mechanisms is unknown. Here, we report that hypoxia induces a rapid and hypoxia-inducible factor–independent induction of histone methylation in a range of human cultured cells. Genomic locations of histone-3 lysine-4 trimethylation (H3K4me3) and H3K36me3 after a brief exposure of cultured cells to hypoxia predict the cell’s transcriptional response several hours later. We show that inactivation of one of the JmjC-containing enzymes, lysine demethylase 5A (KDM5A), mimics hypoxia-induced cellular responses. These results demonstrate that oxygen sensing by chromatin occurs via JmjC-histone demethylase inhibition.

Within the genomes of metazoans, nucleosomes are highly organised adjacent to the binding sites for a subset of transcription factors. Here we have sought to investigate which chromatin remodelling enzymes are responsible for this. We find that the ATP-dependent chromatin remodelling enzyme SNF2H plays a major role organising arrays of nucleosomes adjacent to the binding sites for the architectural transcription factor CTCF sites and acts to promote CTCF binding. At many other factor binding sites SNF2H and the related enzyme SNF2L contribute to nucleosome organisation. The action of SNF2H at CTCF sites is functionally important as depletion of CTCF or SNF2H affects transcription of a common group of genes. This suggests that chromatin remodelling ATPase’s most closely related to the Drosophila ISWI protein contribute to the function of many human gene regulatory elements.

Throughout their lifetime, cells are subject to extrinsic and intrinsic mutational processes leaving behind characteristic signatures in the genome. DNA mismatch repair (MMR) deficiency leads to hypermutation and is found in different cancer types. Although it is possible to associate mutational signatures extracted from human cancers with possible mutational processes, the exact causation is often unknown. Here, we use genome sequencing of and knockouts to reveal the mutational patterns linked to MMR deficiency and their dependency on endogenous replication errors and errors caused by deletion of the polymerase ε subunit Signature extraction from 215 human colorectal and 289 gastric adenocarcinomas revealed three MMR-associated signatures, one of which closely resembles the MMR spectrum and strongly discriminates microsatellite stable and unstable tumors (AUC = 98%). A characteristic difference between human and MMR deficiency is the lack of elevated levels of NG > NTG mutations in likely caused by the absence of cytosine (CpG) methylation in worms The other two human MMR signatures may reflect the interaction between MMR deficiency and other mutagenic processes, but their exact cause remains unknown. In summary, combining information from genetically defined models and cancer samples allows for better aligning mutational signatures to causal mutagenic processes.

Motivation: High-throughput RNA sequencing (RNA-seq) is now the standard method to determine differential gene expression. Identifying differentially expressed genes crucially depends on estimates of read-count variability. These estimates are typically based on statistical models such as the negative binomial distribution, which is employed by the tools edgeR, DESeq and cuffdiff. Until now, the validity of these models has usually been tested on either low-replicate RNA-seq data or simulations.Results: A 48-replicate RNA-seq experiment in yeast was performed and data tested against theoretical models. The observed gene read counts were consistent with both log-normal and negative binomial distributions, while the mean-variance relation followed the line of constant dispersion parameter of ∼0.01. The high-replicate data also allowed for strict quality control and screening of ‘bad’ replicates, which can drastically affect the gene read-count distribution.Availability and implementation: RNA-seq data have been submitted to ENA archive with project ID PRJEB5348.Contact: g.j.barton@dundee.ac.uk

The small ubiquitin-like modifier 2 (SUMO-2) is required for survival when cells are exposed to treatments that induce proteotoxic stress by causing the accumulation of misfolded proteins.Exposure of cells to heat shock or other forms of proteotoxic stress induces the conjugation of SUMO-2 to proteins in the nucleus. Here, we investigated the chromatin landscape of SUMO-2 modifications in response to heat stress. Through chromatin immunoprecipitation assays coupled to high-throughput DNA sequencing and with mRNA sequencing, we showed that in response to heat shock, SUMO-2 accumulated at nucleosome-depleted, active DNA-regulatory elements, which represented binding sites for large protein complexes and were predominantly associated with active genes. However, SUMO did not act as a direct transcriptional repressor or activator of these genes during heat shock. Instead, integration of our results with published proteomics data on heat shock-induced SUMO-2 substrates supports a model in which the conjugation of SUMO-2 to proteins acts as an acute stress response that is required for the stability of protein complexes involved in gene expression and posttranscriptional modification of mRNA. We showed that the conjugation of SUMO-2 to chromatin-associated proteins is an integral component of the proteotoxic stress response, and propose that SUMO-2 fulfills its essential role in cell survival by contributing to the maintenance of protein complex homeostasis.