Senescence-associated reprogramming promotes cancer stemness.

Cellular senescence is a stress-responsive cell-cycle arrest program that terminates the further expansion of (pre-)malignant cells. Key signalling components of the senescence machinery, such as p16INK4a, p21CIP1 and p53, as well as trimethylation of lysine 9 at histone H3 (H3K9me3), also operate as critical regulators of stem-cell functions (which are collectively termed 'stemness'). In cancer cells, a gain of stemness may have profound implications for tumour aggressiveness and clinical outcome. Here we investigated whether chemotherapy-induced senescence could change stem-cell-related properties of malignant cells. Gene expression and functional analyses comparing senescent and non-senescent B-cell lymphomas from Eµ-Myc transgenic mice revealed substantial upregulation of an adult tissue stem-cell signature, activated Wnt signalling, and distinct stem-cell markers in senescence. Using genetically switchable models of senescence targeting H3K9me3 or p53 to mimic spontaneous escape from the arrested condition, we found that cells released from senescence re-entered the cell cycle with strongly enhanced and Wnt-dependent clonogenic growth potential compared to virtually identical populations that had been equally exposed to chemotherapy but had never been senescent. In vivo, these previously senescent cells presented with a much higher tumour initiation potential. Notably, the temporary enforcement of senescence in p53-regulatable models of acute lymphoblastic leukaemia and acute myeloid leukaemia was found to reprogram non-stem bulk leukaemia cells into self-renewing, leukaemia-initiating stem cells. Our data, which are further supported by consistent results in human cancer cell lines and primary samples of human haematological malignancies, reveal that senescence-associated stemness is an unexpected, cell-autonomous feature that exerts its detrimental, highly aggressive growth potential upon escape from cell-cycle blockade, and is enriched in relapse tumours. These findings have profound implications for cancer therapy, and provide new mechanistic insights into the plasticity of cancer cells.

Complexity against current cancer research: Are we on the wrong track?

Cancer genetics has led to major discoveries, including protooncogene and tumor-suppressor concepts, and cancer genomics generated concepts like driver and passenger genes, revealed tumor heterogeneity and clonal evolution. Reconstructing trajectories of tumorigenesis using spatial and single-cell genomics is possible. Patient stratification and prognostic parameters have been improved. Yet, despite these advances, successful translation into targeted therapies has been scarce and mostly limited to kinase inhibitors. Here, we argue that current cancer research may be on the wrong track, by considering cancer more as a "monogenic" disease, trying to extract common information from thousands of patients, while not properly considering complexity and individual diversity. We suggest to empower a systems cancer approach which reconstructs the information network that has been altered by the tumorigenic events, to analyze hierarchies and predict (druggable) key nodes that could interfere with/block the aberrant information transfer. We also argue that the interindividual variability between patients of similar cohorts is too high to extract common polygenic network information from large numbers of patients and argue in favor of an individualized approach. The analysis we propose would require a structured multinational and multidisciplinary effort, in which clinicians, and cancer, developmental, cell and computational biologists together with mathematicians and informaticians develop dynamic regulatory networks which integrate the entire information transfer in and between cells and organs in (patho)physiological conditions, revealing hierarchies and available drugs to interfere with key regulators. Based on this blueprint, the altered information transfer in individual cancers could be modeled and possible targeted (combo)therapies proposed.

Reconstructed cell fate-regulatory programs in stem cells reveal hierarchies and key factors of neurogenesis.

Cell lineages, which shape the body architecture and specify cell functions, derive from the integration of a plethora of cell intrinsic and extrinsic signals. These signals trigger a multiplicity of decisions at several levels to modulate the activity of dynamic gene regulatory networks (GRNs), which ensure both general and cell-specific functions within a given lineage, thereby establishing cell fates. Significant knowledge about these events and the involved key drivers comes from homogeneous cell differentiation models. Even a single chemical trigger, such as the morphogen all-trans retinoic acid (RA), can induce the complex network of gene-regulatory decisions that matures a stem/precursor cell to a particular step within a given lineage. Here we have dissected the GRNs involved in the RA-induced neuronal or endodermal cell fate specification by integrating dynamic RXRA binding, chromatin accessibility, epigenetic promoter epigenetic status, and the transcriptional activity inferred from RNA polymerase II mapping and transcription profiling. Our data reveal how RA induces a network of transcription factors (TFs), which direct the temporal organization of cognate GRNs, thereby driving neuronal/endodermal cell fate specification. Modeling signal transduction propagation using the reconstructed GRNs indicated critical TFs for neuronal cell fate specification, which were confirmed by CRISPR/Cas9-mediated genome editing. Overall, this study demonstrates that a systems view of cell fate specification combined with computational signal transduction models provides the necessary insight in cellular plasticity for cell fate engineering. The present integrated approach can be used to monitor the in vitro capacity of (engineered) cells/tissues to establish cell lineages for regenerative medicine.

The KDM5 family is required for activation of pro-proliferative cell cycle genes during adipocyte differentiation.

The KDM5 family of histone demethylases removes the H3K4 tri-methylation (H3K4me3) mark frequently found at promoter regions of actively transcribed genes and is therefore generally considered to contribute to corepression. In this study, we show that knockdown (KD) of all expressed members of the KDM5 family in white and brown preadipocytes leads to deregulated gene expression and blocks differentiation to mature adipocytes. KDM5 KD leads to a considerable increase in H3K4me3 at promoter regions; however, these changes in H3K4me3 have a limited effect on gene expression per se. By contrast, genome-wide analyses demonstrate that KDM5A is strongly enriched at KDM5-activated promoters, which generally have high levels of H3K4me3 and are associated with highly expressed genes. We show that KDM5-activated genes include a large set of cell cycle regulators and that the KDM5s are necessary for mitotic clonal expansion in 3T3-L1 cells, indicating that KDM5 KD may interfere with differentiation in part by impairing proliferation. Notably, the demethylase activity of KDM5A is required for activation of at least a subset of pro-proliferative cell cycle genes. In conclusion, the KDM5 family acts as dual modulators of gene expression in preadipocytes and is required for early stage differentiation and activation of pro-proliferative cell cycle genes.

Methylation specifies distinct estrogen-induced binding site repertoires of CBP to chromatin.

Multiple signaling pathways ultimately modulate the epigenetic information embedded in the chromatin of gene promoters by recruiting epigenetic enzymes. We found that, in estrogen-regulated gene programming, the acetyltransferase CREB-binding protein (CBP) is specifically and exclusively methylated by the coactivator-associated arginine methyltransferase (CARM1) in vivo. CARM1-dependent CBP methylation and p160 coactivators were required for estrogen-induced recruitment to chromatin targets. Notably, methylation increased the histone acetyltransferase (HAT) activity of CBP and stimulated its autoacetylation. Comparative genome-wide chromatin immunoprecipitation sequencing (ChIP-seq) studies revealed a variety of patterns by which p160, CBP, and methyl-CBP (meCBP) are recruited (or not) by estrogen to chromatin targets. Moreover, significant target gene-specific variation in the recruitment of (1) the p160 RAC3 protein, (2) the fraction of a given meCBP species within the total CBP, and (3) the relative recruitment of different meCBP species suggests the existence of a target gene-specific "fingerprint" for coregulator recruitment. Crossing ChIP-seq and transcriptomics profiles revealed the existence of meCBP "hubs" within the network of estrogen-regulated genes. Together, our data provide evidence for an unprecedented mechanism by which CARM1-dependent CBP methylation results in gene-selective association of estrogen-recruited meCBP species with different HAT activities and specifies distinct target gene hubs, thus diversifying estrogen receptor programming.

The inactive X chromosome is epigenetically unstable and transcriptionally labile in breast cancer.

Disappearance of the Barr body is considered a hallmark of cancer, although whether this corresponds to genetic loss or to epigenetic instability and transcriptional reactivation is unclear. Here we show that breast tumors and cell lines frequently display major epigenetic instability of the inactive X chromosome, with highly abnormal 3D nuclear organization and global perturbations of heterochromatin, including gain of euchromatic marks and aberrant distributions of repressive marks such as H3K27me3 and promoter DNA methylation. Genome-wide profiling of chromatin and transcription reveal modified epigenomic landscapes in cancer cells and a significant degree of aberrant gene activity from the inactive X chromosome, including several genes involved in cancer promotion. We demonstrate that many of these genes are aberrantly reactivated in primary breast tumors, and we further demonstrate that epigenetic instability of the inactive X can lead to perturbed dosage of X-linked factors. Taken together, our study provides the first integrated analysis of the inactive X chromosome in the context of breast cancer and establishes that epigenetic erosion of the inactive X can lead to the disappearance of the Barr body in breast cancer cells. This work offers new insights and opens up the possibility of exploiting the inactive X chromosome as an epigenetic biomarker at the molecular and cytological levels in cancer.

Development of biotin-retinoid conjugates as chemical probes for analysis of retinoid function.

Herein, we report the rational design, synthesis and biological evaluation of conjugates consisting of the synthetic retinoid Am580 and biotin connected via a linker moiety. We found that the linking substructure between the retinoid part and the biotin part is critical for retaining the biological activity. Conjugate 4 with a shorter linker showed similar potency to endogenous retinoid ATRA (1) and the parent compound Am580 (2) for neural differentiation of mouse embryotic carcinoma P19 cells, and showed the same pattern of induction of gene expression. It is expected to be useful as a probe for investigations of retinoid function. The design rationale and structure-activity relationship of the linker moiety are expected to be helpful for developing biotin conjugates of other nuclear receptor ligands.

Epimetheus - a multi-profile normalizer for epigenomic sequencing data.

BACKGROUND: Exponentially increasing numbers of NGS-based epigenomic datasets in public repositories like GEO constitute an enormous source of information that is invaluable for integrative and comparative studies of gene regulatory mechanisms. One of today's challenges for such studies is to identify functionally informative local and global patterns of chromatin states in order to describe the regulatory impact of the epigenome in normal cell physiology and in case of pathological aberrations. Critically, the most preferred Chromatin ImmunoPrecipitation-Sequencing (ChIP-Seq) is inherently prone to significant variability between assays, which poses significant challenge on comparative studies. One challenge concerns data normalization to adjust sequencing depth variation. RESULTS: Currently existing tools either apply linear scaling corrections and/or are restricted to specific genomic regions, which can be prone to biases. To overcome these restrictions without any external biases, we developed Epimetheus, a genome-wide quantile-based multi-profile normalization tool for histone modification data and related datasets. CONCLUSIONS: Epimetheus has been successfully used to normalize epigenomics data in previous studies on X inactivation in breast cancer and in integrative studies of neuronal cell fate acquisition and tumorigenic transformation; Epimetheus is freely available to the scientific community.

LOGIQA: a database dedicated to long-range genome interactions quality assessment.

BACKGROUND: Proximity ligation-mediated methods are essential to study the impact of three-dimensional chromatin organization on gene programming. Albeit significant progress has been made in the development of computational tools that assess long-range chromatin interactions, next to nothing is known about the quality of the generated datasets. METHOD: We have developed LOGIQA ( www.ngs-qc.org/logiqa ), a database hosting quality scores for long-range genome interaction assays, accessible through a user-friendly web-based environment. RESULTS: Currently, LOGIQA harbors QC scores for >900 datasets, which provides a global view of their relative quality and reveals the impact of genome size, coverage and other technical aspects. LOGIQA provides a user-friendly dataset query panel and a genome viewer to assess local genome-interaction maps at different resolution and quality-assessment conditions. CONCLUSIONS: LOGIQA is the first database hosting quality scores dedicated to long-range chromatin interaction assays, which in addition provides a platform for visualizing genome interactions made available by the scientific community.

Genome-wide studies of nuclear receptors in cell fate decisions.

Nuclear receptors (NRs) are important mediators of the information encoded in the chemical structure of its corresponding ligand, as they interpret such information in the context of the cell identity and physiological status and convert it into sequential transcription regulatory events. At the cell level this can result in temporally coordinated processes such as cell fate transitions, which comprise the regulation of a plethora of gene programs including among others regulation of cell proliferation, metabolism and specific functionalities that are acquired by the differentiated cell. While both the early steps of nuclear receptor function and their impact on animal/organ physiology is rather well understood, little is known about the dynamic gene networks that ultimately cause a particular (cell) physiological phenomenon induced by the cognate NR ligand/hormone. Thanks to advances in massive parallel sequencing and bioinformatics analyses of genome-wide data sets, time has come for the development of NR systems biology. Indeed it is now possible to integrate global transcription factor binding, epigenetic chromatin histone and DNA modification patterns with transcriptomes and 3-dimensional chromatin structures, extract decision points in temporal studies and decipher the temporal control of gene networks that are the ultimate genetic readouts of NR ligand-induced physiological phenomena. In this review we will summarize the chronology of the development of increasingly larger data sets for NR action, with a particular focus on studies performed with the RAR/RXR nuclear receptor family, and discuss the present attempts to integrate a multitude of genome-wide data sets in the ultimate context of the temporal 3-dimensional chromatin structure.

Thioether analogues of disulfide-bridged cyclic peptides targeting death receptor 5: conformational analysis, dimerisation and consequences for receptor activation.

Cyclic peptides containing redox-stable thioether bridges might provide a useful alternative to disulfide-bridged bioactive peptides. We report the effect of replacing the disulfide bridge with a lanthionine linkage in a 16-mer cyclic peptide that binds to death receptor 5 (DR5, TRAIL-R2). Upon covalent oligomerisation, the disulfide-bridged peptide has previously shown similar behaviour to that of TNF-related apoptosis inducing ligand (TRAIL), by selectively triggering the DR5 cell death pathway. The structural and biological properties of the DR5-binding peptide and its desulfurised analogue were compared. Surface plasmon resonance (SPR) data suggest that these peptides bind DR5 with comparable affinities. The same holds true for dimeric versions of these peptides: the thioether is able to induce DR5-mediated apoptosis of BJAB lymphoma and tumorigenic BJELR cells, albeit to a slightly lower extent compared to its disulfide homologue. NMR analysis revealed subtle variation in the conformations of the two peptides and suggests that the thioether peptide is slightly less folded than its disulfide homologue. These observations could account for the different capability of the two dimers to cluster DR5 receptors on the cell surface and to trigger apoptosis. Nevertheless, our results suggest that the thioether peptide is a potential candidate for evaluation in animal models.

Recruitment of RXR by homotetrameric RARalpha fusion proteins is essential for transformation.

While formation of higher-order oncogenic transcriptional complexes is critical for RARalpha fusion proteins in acute promyelocytic leukemia, the essential components and their roles in mediating transformation are still largely unknown. To this end, the present study demonstrates that homodimerization is not sufficient for RARalpha fusion-mediated transformation, which requires higher-order homotetramerization. Surprisingly, intrinsic homo-oligomeric DNA binding by the fusion proteins is also dispensable. Importantly, higher-order RXR/RARalpha fusion hetero-oligomeric complexes that aberrantly recruit transcriptional corepressors to downstream targets are essential for transformation. Intervention of RXR-dependent pathways by panRXR-agonists or RXRalpha shRNAs suppresses RARalpha fusion-mediated transformation. Taken together, these results define the oncogenic threshold for self-association and reveal the pathological significance of higher-order RARalpha fusion/RXR hetero-oligomeric complexes and their potential value as a therapeutic target.

A quality control system for profiles obtained by ChIP sequencing.

The absence of a quality control (QC) system is a major weakness for the comparative analysis of genome-wide profiles generated by next-generation sequencing (NGS). This concerns particularly genome binding/occupancy profiling assays like chromatin immunoprecipitation (ChIP-seq) but also related enrichment-based studies like methylated DNA immunoprecipitation/methylated DNA binding domain sequencing, global run on sequencing or RNA-seq. Importantly, QC assessment may significantly improve multidimensional comparisons that have great promise for extracting information from combinatorial analyses of the global profiles established for chromatin modifications, the bindings of epigenetic and chromatin-modifying enzymes/machineries, RNA polymerases and transcription factors and total, nascent or ribosome-bound RNAs. Here we present an approach that associates global and local QC indicators to ChIP-seq data sets as well as to a variety of enrichment-based studies by NGS. This QC system was used to certify >5600 publicly available data sets, hosted in a database for data mining and comparative QC analyses.

Integrative genomics to dissect retinoid functions.

Retinoids and rexinoids, as all other ligands of the nuclear receptor (NR) family, act as ligand-regulated trans-acting transcription factors that bind to cis-acting DNA regulatory elements in the promoter regions of target genes (for reviews see [12, 22, 23, 26, 36]). Ligand binding modulates the communication functions of the receptor with the intracellular environment, which essentially entails receptor-protein and receptor-DNA or receptor-chromatin interactions. In this communication network, the receptor simultaneously serves as both intracellular sensor and regulator of cell/organ functions. Receptors are "intelligent" mediators of the information encoded in the chemical structure of a nuclear receptor ligand, as they interpret this information in the context of cellular identity and cell-physiological status and convert it into a dynamic chain of receptor-protein and receptor-DNA interactions. To process input and output information, they are composed of a modular structure with several domains that have evolved to exert particular molecular recognition functions. As detailed in other chapters in this volume, the main functional domains are the DNA-binding (DBD) and ligand-binding (LBD) [5-7, 38, 56, 71]. The LBD serves as a dual input-output information processor. Inputs, such as ligand binding or receptor phosphorylations, induce allosteric changes in receptor surfaces that serve as docking sites for outputs, such as subunits of transcription and epigenetic machineries or enzyme complexes. The complexity of input and output signals and their interdependencies is far from being understood.

Single-tube linear DNA amplification (LinDA) for robust ChIP-seq.

Genome-wide profiling of transcription factors based on massive parallel sequencing of immunoprecipitated chromatin (ChIP-seq) requires nanogram amounts of DNA. Here we describe a high-fidelity, single-tube linear DNA amplification method (LinDA) for ChIP-seq and reChIP-seq with picogram DNA amounts obtained from a few thousand cells. This amplification technology will facilitate global analyses of transcription-factor binding and chromatin with very small cell populations, such as stem or cancer-initiating cells.