Analysis of ChIP-seq Data in R/Bioconductor

The development of novel high-throughput sequencing methods for ChIP (chromatin immunoprecipitation) has provided a very powerful tool to study gene regulation in multiple conditions at unprecedented resolution and scale. Proactive quality-control and appropriate data analysis techniques are of critical importance to extract the most meaningful results from the data. Over the last years, an array of R/Bioconductor tools has been developed allowing researchers to process and analyze ChIP-seq data. This chapter provides an overview of the methods available to analyze ChIP-seq data based primarily on software packages from the open-source Bioconductor project. Protocols described in this chapter cover basic steps including data alignment, peak calling, quality control and data visualization, as well as more complex methods such as the identification of differentially bound regions and functional analyses to annotate regulatory regions. The steps in the data analysis process were demonstrated on publicly available data sets and will serve as a demonstration of the computational procedures routinely used for the analysis of ChIP-seq data in R/Bioconductor, from which readers can construct their own analysis pipelines.

Epileptic Seizure Prediction Using Big Data and Deep Learning: Toward a Mobile System

BACKGROUND

Seizure prediction can increase independence and allow preventative treatment for patients with epilepsy. We present a proof-of-concept for a seizure prediction system that is accurate, fully automated, patient-specific, and tunable to an individual's needs.

METHODS

Intracranial electroencephalography (iEEG) data of ten patients obtained from a seizure advisory system were analyzed as part of a pseudoprospective seizure prediction study. First, a deep learning classifier was trained to distinguish between preictal and interictal signals. Second, classifier performance was tested on held-out iEEG data from all patients and benchmarked against the performance of a random predictor. Third, the prediction system was tuned so sensitivity or time in warning could be prioritized by the patient. Finally, a demonstration of the feasibility of deployment of the prediction system onto an ultra-low power neuromorphic chip for autonomous operation on a wearable device is provided.

RESULTS

The prediction system achieved mean sensitivity of 69% and mean time in warning of 27%, significantly surpassing an equivalent random predictor for all patients by 42%.

CONCLUSION

This study demonstrates that deep learning in combination with neuromorphic hardware can provide the basis for a wearable, real-time, always-on, patient-specific seizure warning system with low power consumption and reliable long-term performance.

An efficient method for generation of bi-allelic null mutant mouse embryonic stem cells and its application for investigating epigenetic modifiers

Mouse embryonic stem (ES) cells are a popular model system to study biological processes, though uncovering recessive phenotypes requires inactivating both alleles. Building upon resources from the International Knockout Mouse Consortium (IKMC), we developed a targeting vector for second allele inactivation in conditional-ready IKMC 'knockout-first' ES cell lines. We applied our technology to several epigenetic regulators, recovering bi-allelic targeted clones with a high efficiency of 60% and used Flp recombinase to restore expression in two null cell lines to demonstrate how our system confirms causality through mutant phenotype reversion. We designed our strategy to select against re-targeting the 'knockout-first' allele and identify essential genes in ES cells, including the histone methyltransferase Setdb1. For confirmation, we exploited the flexibility of our system, enabling tamoxifen inducible conditional gene ablation while controlling for genetic background and tamoxifen effects. Setdb1 ablated ES cells exhibit severe growth inhibition, which is not rescued by exogenous Nanog expression or culturing in naive pluripotency '2i' media, suggesting that the self-renewal defect is mediated through pluripotency network independent pathways. Our strategy to generate null mutant mouse ES cells is applicable to thousands of genes and repurposes existing IKMC Intermediate Vectors.

Coupling shRNA screens with single-cell RNA-seq identifies a dual role for mTOR in reprogramming-induced senescence

Aarts et al. developed an innovative approach that integrates single-cell RNA sequencing with a shRNA screen in primary human fibroblasts expressing OCT4, SOX2, KLF4, and cMYC to investigate the mechanism of action of the identified candidates. This approach unveiled regulation of senescence as a novel way by which mechanistic target of rapamycin (mTOR) influences reprogramming.

Expression of the transcription factors OCT4, SOX2, KLF4, and cMYC (OSKM) reprograms somatic cells into induced pluripotent stem cells (iPSCs). Reprogramming is a slow and inefficient process, suggesting the presence of safeguarding mechanisms that counteract cell fate conversion. One such mechanism is senescence. To identify modulators of reprogramming-induced senescence, we performed a genome-wide shRNA screen in primary human fibroblasts expressing OSKM. In the screen, we identified novel mediators of OSKM-induced senescence and validated previously implicated genes such as CDKN1A. We developed an innovative approach that integrates single-cell RNA sequencing (scRNA-seq) with the shRNA screen to investigate the mechanism of action of the identified candidates. Our data unveiled regulation of senescence as a novel way by which mechanistic target of rapamycin (mTOR) influences reprogramming. On one hand, mTOR inhibition blunts the induction of cyclin-dependent kinase (CDK) inhibitors (CDKIs), including p16^INK4a, p21^CIP1, and p15^INK4b, preventing OSKM-induced senescence. On the other hand, inhibition of mTOR blunts the senescence-associated secretory phenotype (SASP), which itself favors reprogramming. These contrasting actions contribute to explain the complex effect that mTOR has on reprogramming. Overall, our study highlights the advantage of combining functional screens with scRNA-seq to accelerate the discovery of pathways controlling complex phenotypes.

Ret receptor tyrosine kinase sustains proliferation and tissue maturation in intestinal epithelia

Expression of the Ret receptor tyrosine kinase is a defining feature of enteric neurons. Its importance is underscored by the effects of its mutation in Hirschsprung disease, leading to absence of gut innervation and severe gastrointestinal symptoms. We report a new and physiologically significant site of Ret expression in the intestine: the intestinal epithelium. Experiments in Drosophila indicate that Ret is expressed both by enteric neurons and adult intestinal epithelial progenitors, which require Ret to sustain their proliferation. Mechanistically, Ret is engaged in a positive feedback loop with Wnt/Wingless signalling, modulated by Src and Fak kinases. We find that Ret is also expressed by the developing intestinal epithelium of mice, where its expression is maintained into the adult stage in a subset of enteroendocrine/enterochromaffin cells. Mouse organoid experiments point to an intrinsic role for Ret in promoting epithelial maturation and regulating Wnt signalling. Our findings reveal evolutionary conservation of the positive Ret/Wnt signalling feedback in both developmental and homeostatic contexts. They also suggest an epithelial contribution to Ret loss‐of‐function disorders such as Hirschsprung disease.

Primary Palliative Care Education: A Pilot Survey

The demand for palliative services is outpacing the availability of specialist palliative care clinicians. One strategy to fill this gap is to improve "primary palliative care" skills and knowledge of all clinicians who care for seriously ill patients. Previous educational efforts have shown mixed results, and one possible explanation is unrecognized discordance of educational goals between those offering education and potential primary palliative care learners. The article describes the results and feasibility of a needs assessment survey comparing interest in palliative care education topics and settings among both palliative care specialists (PCS) and nonpalliative care specialists (NPCS). This is the first attempt to measure the perceived importance of primary palliative care topics and preferences about learning settings from the perspectives of both NPCS and PCS. The results suggest substantial areas of both concordant and discordant opinions with respect to educational topics and learning settings. Such data are essential to guide primary palliative care educational efforts. Future work will be needed to determine whether these results are consistent across diverse health systems and what variables influence educational preferences.

A high-resolution map of transcriptional repression

Turning genes on and off is essential for development and homeostasis, yet little is known about the sequence and causal role of chromatin state changes during the repression of active genes. This is surprising, as defective gene silencing underlies developmental abnormalities and disease. Here we delineate the sequence and functional contribution of transcriptional repression mechanisms at high temporal resolution. Inducible entry of the NuRD-interacting transcriptional regulator Ikaros into mouse pre-B cell nuclei triggered immediate binding to target gene promoters. Rapid RNAP2 eviction, transcriptional shutdown, nucleosome invasion, and reduced transcriptional activator binding required chromatin remodeling by NuRD-associated Mi2beta/CHD4, but were independent of HDAC activity. Histone deacetylation occurred after transcriptional repression. Nevertheless, HDAC activity contributed to stable gene silencing. Hence, high resolution mapping of transcriptional repression reveals complex and interdependent mechanisms that underpin rapid transitions between transcriptional states, and elucidates the temporal order, functional role and mechanistic separation of NuRD-associated enzymatic activities.

DOI:http://dx.doi.org/10.7554/eLife.22767.001

SWI/SNF regulates a transcriptional program that induces senescence to prevent liver cancer

Here, Tordella et. al identified senescence regulators relevant to cancer by screening an shRNA library targeting genes deleted in hepatocellular carcinoma (HCC). They show that knockdown of the SWI/SNF component ARID1B prevents oncogene-induced senescence and cooperates with RAS to induce liver tumors, and their results provide new insights into the mechanisms by which epigenetic regulators can affect tumor progression.

Oncogene-induced senescence (OIS) is a potent tumor suppressor mechanism. To identify senescence regulators relevant to cancer, we screened an shRNA library targeting genes deleted in hepatocellular carcinoma (HCC). Here, we describe how knockdown of the SWI/SNF component ARID1B prevents OIS and cooperates with RAS to induce liver tumors. ARID1B controls p16^INK4a and p21^CIP1a transcription but also regulates DNA damage, oxidative stress, and p53 induction, suggesting that SWI/SNF uses additional mechanisms to regulate senescence. To systematically identify SWI/SNF targets regulating senescence, we carried out a focused shRNA screen. We discovered several new senescence regulators, including ENTPD7, an enzyme that hydrolyses nucleotides. ENTPD7 affects oxidative stress, DNA damage, and senescence. Importantly, expression of ENTPD7 or inhibition of nucleotide synthesis in ARID1B-depleted cells results in re-establishment of senescence. Our results identify novel mechanisms by which epigenetic regulators can affect tumor progression and suggest that prosenescence therapies could be employed against SWI/SNF-mutated cancers.

MicroRNAs of the miR-290-295 Family Maintain Bivalency in Mouse Embryonic Stem Cells

Numerous developmentally regulated genes in mouse embryonic stem cells (ESCs) are marked by both active (H3K4me3)- and polycomb group (PcG)-mediated repressive (H3K27me3) histone modifications. This bivalent state is thought to be important for transcriptional poising, but the mechanisms that regulate bivalent genes and the bivalent state remain incompletely understood. Examining the contribution of microRNAs (miRNAs) to the regulation of bivalent genes, we found that the miRNA biogenesis enzyme DICER was required for the binding of the PRC2 core components EZH2 and SUZ12, and for the presence of the PRC2-mediated histone modification H3K27me3 at many bivalent genes. Genes that lost bivalency were preferentially upregulated at the mRNA and protein levels. Finally, reconstituting Dicer-deficient ESCs with ESC miRNAs restored bivalent gene repression and PRC2 binding at formerly bivalent genes. Therefore, miRNAs regulate bivalent genes and the bivalent state itself.

p53 Enables metabolic fitness and self-renewal of nephron progenitor cells

Contrary to its classic role in restraining cell proliferation, we demonstrate here a divergent function of p53 in the maintenance of self-renewal of the nephron progenitor pool in the embryonic mouse kidney. Nephron endowment is regulated by progenitor availability and differentiation potential. Conditional deletion of p53 in nephron progenitor cells (Six2Cre(+);p53(fl/fl)) induces progressive depletion of Cited1(+)/Six2(+) self-renewing progenitors and loss of cap mesenchyme (CM) integrity. The Six2(p53-null) CM is disorganized, with interspersed stromal cells and an absence of a distinct CM-epithelia and CM-stroma interface. Impaired cell adhesion and epithelialization are indicated by decreased E-cadherin and NCAM expression and by ineffective differentiation in response to Wnt induction. The Six2Cre(+);p53(fl/fl) cap has 30% fewer Six2(GFP(+)) cells. Apoptotic index is unchanged, whereas proliferation index is significantly reduced in accordance with cell cycle analysis showing disproportionately fewer Six2Cre(+);p53(fl/fl) cells in the S and G2/M phases compared with Six2Cre(+);p53(+/+) cells. Mutant kidneys are hypoplastic with fewer generations of nascent nephrons. A significant increase in mean arterial pressure is observed in early adulthood in both germline and conditional Six2(p53-null) mice, linking p53-mediated defects in kidney development to hypertension. RNA-Seq analyses of FACS-isolated wild-type and Six2(GFP(+)) CM cells revealed that the top downregulated genes in Six2Cre(+);p53(fl/fl) CM belong to glucose metabolism and adhesion and/or migration pathways. Mutant cells exhibit a ∼ 50% decrease in ATP levels and a 30% decrease in levels of reactive oxygen species, indicating energy metabolism dysfunction. In summary, our data indicate a novel role for p53 in enabling the metabolic fitness and self-renewal of nephron progenitors.

mTOR regulates MAPKAPK2 translation to control the senescence-associated secretory phenotype

Senescent cells secrete a combination of factors collectively known as the senescence-associated secretory phenotype (SASP). The SASP reinforces senescence and activates an immune surveillance response, but it can also show pro-tumorigenic properties and contribute to age-related pathologies. In a drug screen to find new SASP regulators, we uncovered the mTOR inhibitor rapamycin as a potent SASP suppressor. Here we report a mechanism by which mTOR controls the SASP by differentially regulating the translation of the MK2 (also known as MAPKAPK2) kinase through 4EBP1. In turn, MAPKAPK2 phosphorylates the RNA-binding protein ZFP36L1 during senescence, inhibiting its ability to degrade the transcripts of numerous SASP components. Consequently, mTOR inhibition or constitutive activation of ZFP36L1 impairs the non-cell-autonomous effects of senescent cells in both tumour-suppressive and tumour-promoting contexts. Altogether, our results place regulation of the SASP as a key mechanism by which mTOR could influence cancer, age-related diseases and immune responses.

Jarid2 Coordinates Nanog Expression and PCP/Wnt Signaling Required for Efficient ESC Differentiation and Early Embryo Development

Jarid2 is part of the Polycomb Repressor complex 2 (PRC2) responsible for genome-wide H3K27me3 deposition. Unlike other PRC2-deficient embryonic stem cells (ESCs), however, Jarid2-deficient ESCs show a severe differentiation block, altered colony morphology, and distinctive patterns of deregulated gene expression. Here, we show that Jarid2^−/− ESCs express constitutively high levels of Nanog but reduced PCP signaling components Wnt9a, Prickle1, and Fzd2 and lowered β-catenin activity. Depletion of Wnt9a/Prickle1/Fzd2 from wild-type ESCs or overexpression of Nanog largely phenocopies these cellular defects. Co-culture of Jarid2^−/− with wild-type ESCs restores variable Nanog expression and β-catenin activity and can partially rescue the differentiation block of mutant cells. In addition, we show that ESCs lacking Jarid2 or Wnt9a/Prickle1/Fzd2 or overexpressing Nanog induce multiple ICM formation when injected into normal E3.5 blastocysts. These data describe a previously unrecognized role for Jarid2 in regulating a core pluripotency and Wnt/PCP signaling circuit that is important for ESC differentiation and for pre-implantation development.

•

ESCs lacking Jarid2 show constitutive Nanog expression

•

ESCs lacking Jarid2 have reduced PCP/Wnt signaling

•

Co-culture of Jarid2-null and WT ESCs restores differentiation capability

•

Jarid2-null ESCs form more than one ICM upon injection to E3.5 mouse blastocysts

5-hydroxymethylcytosine marks promoters in colon that resist DNA hypermethylation in cancer

BACKGROUND

The discovery of cytosine hydroxymethylation (5hmC) as a mechanism that potentially controls DNA methylation changes typical of neoplasia prompted us to investigate its behaviour in colon cancer. 5hmC is globally reduced in proliferating cells such as colon tumours and the gut crypt progenitors, from which tumours can arise.

RESULTS

Here, we show that colorectal tumours and cancer cells express Ten-Eleven-Translocation (TET) transcripts at levels similar to normal tissues. Genome-wide analyses show that promoters marked by 5hmC in normal tissue, and those identified as TET2 targets in colorectal cancer cells, are resistant to methylation gain in cancer. In vitro studies of TET2 in cancer cells confirm that these promoters are resistant to methylation gain independently of sustained TET2 expression. We also find that a considerable number of the methylation gain-resistant promoters marked by 5hmC in normal colon overlap with those that are marked with poised bivalent histone modifications in embryonic stem cells.

CONCLUSIONS

Together our results indicate that promoters that acquire 5hmC upon normal colon differentiation are innately resistant to neoplastic hypermethylation by mechanisms that do not require high levels of 5hmC in tumours. Our study highlights the potential of cytosine modifications as biomarkers of cancerous cell proliferation.

Spatial enhancer clustering and regulation of enhancer-proximal genes by cohesin

In addition to mediating sister chromatid cohesion during the cell cycle, the cohesin complex associates with CTCF and with active gene regulatory elements to form long-range interactions between its binding sites. Genome-wide chromosome conformation capture had shown that cohesin's main role in interphase genome organization is in mediating interactions within architectural chromosome compartments, rather than specifying compartments per se. However, it remains unclear how cohesin-mediated interactions contribute to the regulation of gene expression. We have found that the binding of CTCF and cohesin is highly enriched at enhancers and in particular at enhancer arrays or "super-enhancers" in mouse thymocytes. Using local and global chromosome conformation capture, we demonstrate that enhancer elements associate not just in linear sequence, but also in 3D, and that spatial enhancer clustering is facilitated by cohesin. The conditional deletion of cohesin from noncycling thymocytes preserved enhancer position, H3K27ac, H4K4me1, and enhancer transcription, but weakened interactions between enhancers. Interestingly, ∼ 50% of deregulated genes reside in the vicinity of enhancer elements, suggesting that cohesin regulates gene expression through spatial clustering of enhancer elements. We propose a model for cohesin-dependent gene regulation in which spatial clustering of enhancer elements acts as a unified mechanism for both enhancer-promoter "connections" and "insulation."

Initiation and maintenance of pluripotency gene expression in the absence of cohesin

Cohesin is implicated in establishing and maintaining pluripotency. Whether this is because of essential cohesin functions in the cell cycle or in gene regulation is unknown. Here we tested cohesin's contribution to reprogramming in systems that reactivate the expression of pluripotency genes in the absence of proliferation (embryonic stem [ES] cell heterokaryons) or DNA replication (nuclear transfer). Contrary to expectations, cohesin depletion enhanced the ability of ES cells to initiate somatic cell reprogramming in heterokaryons. This was explained by increased c-Myc (Myc) expression in cohesin-depleted ES cells, which promoted DNA replication-dependent reprogramming of somatic fusion partners. In contrast, cohesin-depleted somatic cells were poorly reprogrammed in heterokaryons, due in part to defective DNA replication. Pluripotency gene induction was rescued by Myc, which restored DNA replication, and by nuclear transfer, where reprogramming does not require DNA replication. These results redefine cohesin's role in pluripotency and reveal a novel function for Myc in promoting the replication-dependent reprogramming of somatic nuclei.

Salt-inducible kinase 2 regulates mitotic progression and transcription in prostate cancer

Salt-inducible kinase 2 (SIK2) is a multifunctional kinase of the AMPK family that plays a role in CREB1-mediated gene transcription and was recently reported to have therapeutic potential in ovarian cancer. The expression of this kinase was investigated in prostate cancer clinical specimens. Interestingly, auto-antibodies against SIK2 were increased in the plasma of patients with aggressive disease. Examination of SIK2 in prostate cancer cells found that it functions both as a positive regulator of cell-cycle progression and a negative regulator of CREB1 activity. Knockdown of SIK2 inhibited cell growth, delayed cell-cycle progression, induced cell death, and enhanced CREB1 activity. Expression of a kinase-dead mutant of SIK2 also inhibited cell growth, induced cell death, and enhanced CREB1 activity. Treatment with a small-molecule SIK2 inhibitor (ARN-3236), currently in preclinical development, also led to enhanced CREB1 activity in a dose- and time-dependent manner. Because CREB1 is a transcription factor and proto-oncogene, it was posited that the effects of SIK2 on cell proliferation and viability might be mediated by changes in gene expression. To test this, gene expression array profiling was performed and while SIK2 knockdown or overexpression of the kinase-dead mutant affected established CREB1 target genes; the overlap with transcripts regulated by forskolin (FSK), the adenylate cyclase/CREB1 pathway activator, was incomplete.

IMPLICATIONS

This study demonstrates that targeting SIK2 genetically or therapeutically will have pleiotropic effects on cell-cycle progression and transcription factor activation, which should be accounted for when characterizing SIK2 inhibitors.