Comprehensive characterization of the epigenetic landscape in Multiple Myeloma.

Background: Human multiple myeloma (MM) cell lines (HMCLs) have been widely used to understand the molecular processes that drive MM biology. Epigenetic modifications are involved in MM development, progression, and drug resistance. A comprehensive characterization of the epigenetic landscape of MM would advance our understanding of MM pathophysiology and may attempt to identify new therapeutic targets. Methods: We performed chromatin immunoprecipitation sequencing to analyze histone mark changes (H3K4me1, H3K4me3, H3K9me3, H3K27ac, H3K27me3 and H3K36me3) on 16 HMCLs. Results: Differential analysis of histone modification profiles highlighted links between histone modifications and cytogenetic abnormalities or recurrent mutations. Using histone modifications associated to enhancer regions, we identified super-enhancers (SE) associated with genes involved in MM biology. We also identified promoters of genes enriched in H3K9me3 and H3K27me3 repressive marks associated to potential tumor suppressor functions. The prognostic value of genes associated with repressive domains and SE was used to build two distinct scores identifying high-risk MM patients in two independent cohorts (CoMMpass cohort; n = 674 and Montpellier cohort; n = 69). Finally, we explored H3K4me3 marks comparing drug-resistant and -sensitive HMCLs to identify regions involved in drug resistance. From these data, we developed epigenetic biomarkers based on the H3K4me3 modification predicting MM cell response to lenalidomide and histone deacetylase inhibitors (HDACi). Conclusions: The epigenetic landscape of MM cells represents a unique resource for future biological studies. Furthermore, risk-scores based on SE and repressive regions together with epigenetic biomarkers of drug response could represent new tools for precision medicine in MM.

SETDB1/NSD-dependent H3K9me3/H3K36me3 dual heterochromatin maintains gene expression profiles by bookmarking poised enhancers.

Gene silencing by heterochromatin plays a crucial role in cell identity. Here, we characterize the localization, the biogenesis, and the function of an atypical heterochromatin, which is simultaneously enriched in the typical H3K9me3 mark and in H3K36me3, a histone mark usually associated with gene expression. We identified thousands of dual regions in mouse embryonic stem (ES) cells that rely on the histone methyltransferases SET domain bifurcated 1 (SETDB1) and nuclear set domain (NSD)-containing proteins to generate H3K9me3 and H3K36me3, respectively. Upon SETDB1 removal, dual domains lose both marks, gain signatures of active enhancers, and come into contact with upregulated genes, suggesting that it might be an important pathway by which genes are controlled by heterochromatin. In differentiated tissues, a subset of these dual domains is destabilized and becomes enriched in active enhancer marks, providing a mechanistic insight into the involvement of heterochromatin in the maintenance of cell identity.

Regulation of single-cell genome organization into TADs and chromatin nanodomains.

The genome folds into a hierarchy of three-dimensional structures within the nucleus. At the sub-megabase scale, chromosomes form topologically associating domains (TADs)1-4. However, how TADs fold in single cells is elusive. Here, we reveal TAD features inaccessible to cell population analysis by using super-resolution microscopy. TAD structures and physical insulation associated with their borders are variable between individual cells, yet chromatin intermingling is enriched within TADs compared to adjacent TADs in most cells. The spatial segregation of TADs is further exacerbated during cell differentiation. Favored interactions within TADs are regulated by cohesin and CTCF through distinct mechanisms: cohesin generates chromatin contacts and intermingling while CTCF prevents inter-TAD contacts. Furthermore, TADs are subdivided into discrete nanodomains, which persist in cells depleted of CTCF or cohesin, whereas disruption of nucleosome contacts alters their structural organization. Altogether, these results provide a physical basis for the folding of individual chromosomes at the nanoscale.

4D Genome Rewiring during Oncogene-Induced and Replicative Senescence.

To understand the role of the extensive senescence-associated 3D genome reorganization, we generated genome-wide chromatin interaction maps, epigenome, replication-timing, whole-genome bisulfite sequencing, and gene expression profiles from cells entering replicative senescence (RS) or upon oncogene-induced senescence (OIS). We identify senescence-associated heterochromatin domains (SAHDs). Differential intra- versus inter-SAHD interactions lead to the formation of senescence-associated heterochromatin foci (SAHFs) in OIS but not in RS. This OIS-specific configuration brings active genes located in genomic regions adjacent to SAHDs in close spatial proximity and favors their expression. We also identify DNMT1 as a factor that induces SAHFs by promoting HMGA2 expression. Upon DNMT1 depletion, OIS cells transition to a 3D genome conformation akin to that of cells in replicative senescence. These data show how multi-omics and imaging can identify critical features of RS and OIS and discover determinants of acute senescence and SAHF formation.

Microscopy-Based Chromosome Conformation Capture Enables Simultaneous Visualization of Genome Organization and Transcription in Intact Organisms.

Eukaryotic chromosomes are organized in multiple scales, from nucleosomes to chromosome territories. Recently, genome-wide methods identified an intermediate level of chromosome organization, topologically associating domains (TADs), that play key roles in transcriptional regulation. However, these methods cannot directly examine the interplay between transcriptional activation and chromosome architecture while maintaining spatial information. Here we present a multiplexed, sequential imaging approach (Hi-M) that permits simultaneous detection of chromosome organization and transcription in single nuclei. This allowed us to unveil the changes in 3D chromatin organization occurring upon transcriptional activation and homologous chromosome unpairing during awakening of the zygotic genome in intact Drosophila embryos. Excitingly, the ability of Hi-M to explore the multi-scale chromosome architecture with spatial resolution at different stages of development or during the cell cycle will be key to understanding the mechanisms and consequences of the 4D organization of the genome.

Polycomb-Dependent Chromatin Looping Contributes to Gene Silencing during Drosophila Development.

Interphase chromatin is organized into topologically associating domains (TADs). Within TADs, chromatin looping interactions are formed between DNA regulatory elements, but their functional importance for the establishment of the 3D genome organization and gene regulation during development is unclear. Using high-resolution Hi-C experiments, we analyze higher order 3D chromatin organization during Drosophila embryogenesis and identify active and repressive chromatin loops that are established with different kinetics and depend on distinct factors: Zelda-dependent active loops are formed before the midblastula transition between transcribed genes over long distances. Repressive loops within polycomb domains are formed after the midblastula transition between polycomb response elements by the action of GAGA factor and polycomb proteins. Perturbation of PRE function by CRISPR/Cas9 genome engineering affects polycomb domain formation and destabilizes polycomb-mediated silencing. Preventing loop formation without removal of polycomb components also decreases silencing efficiency, suggesting that chromatin architecture can play instructive roles in gene regulation during development. VIDEO ABSTRACT.

TADs are 3D structural units of higher-order chromosome organization in Drosophila.

Deciphering the rules of genome folding in the cell nucleus is essential to understand its functions. Recent chromosome conformation capture (Hi-C) studies have revealed that the genome is partitioned into topologically associating domains (TADs), which demarcate functional epigenetic domains defined by combinations of specific chromatin marks. However, whether TADs are true physical units in each cell nucleus or whether they reflect statistical frequencies of measured interactions within cell populations is unclear. Using a combination of Hi-C, three-dimensional (3D) fluorescent in situ hybridization, super-resolution microscopy, and polymer modeling, we provide an integrative view of chromatin folding in Drosophila. We observed that repressed TADs form a succession of discrete nanocompartments, interspersed by less condensed active regions. Single-cell analysis revealed a consistent TAD-based physical compartmentalization of the chromatin fiber, with some degree of heterogeneity in intra-TAD conformations and in cis and trans inter-TAD contact events. These results indicate that TADs are fundamental 3D genome units that engage in dynamic higher-order inter-TAD connections. This domain-based architecture is likely to play a major role in regulatory transactions during DNA-dependent processes.

Multiscale 3D Genome Rewiring during Mouse Neural Development.

Chromosome conformation capture technologies have revealed important insights into genome folding. Yet, how spatial genome architecture is related to gene expression and cell fate remains unclear. We comprehensively mapped 3D chromatin organization during mouse neural differentiation in vitro and in vivo, generating the highest-resolution Hi-C maps available to date. We found that transcription is correlated with chromatin insulation and long-range interactions, but dCas9-mediated activation is insufficient for creating TAD boundaries de novo. Additionally, we discovered long-range contacts between gene bodies of exon-rich, active genes in all cell types. During neural differentiation, contacts between active TADs become less pronounced while inactive TADs interact more strongly. An extensive Polycomb network in stem cells is disrupted, while dynamic interactions between neural transcription factors appear in vivo. Finally, cell type-specific enhancer-promoter contacts are established concomitant to gene expression. This work shows that multiple factors influence the dynamics of chromatin interactions in development.

ERCC1-XPF cooperates with CTCF and cohesin to facilitate the developmental silencing of imprinted genes.

Inborn defects in DNA repair are associated with complex developmental disorders whose causal mechanisms are poorly understood. Using an in vivo biotinylation tagging approach in mice, we show that the nucleotide excision repair (NER) structure-specific endonuclease ERCC1-XPF complex interacts with the insulator binding protein CTCF, the cohesin subunits SMC1A and SMC3 and with MBD2; the factors co-localize with ATRX at the promoters and control regions (ICRs) of imprinted genes during postnatal hepatic development. Loss of Ercc1 or exposure to MMC triggers the localization of CTCF to heterochromatin, the dissociation of the CTCF-cohesin complex and ATRX from promoters and ICRs, altered histone marks and the aberrant developmental expression of imprinted genes without altering DNA methylation. We propose that ERCC1-XPF cooperates with CTCF and cohesin to facilitate the developmental silencing of imprinted genes and that persistent DNA damage triggers chromatin changes that affect gene expression programs associated with NER disorders.

GATA-1 Inhibits PU.1 Gene via DNA and Histone H3K9 Methylation of Its Distal Enhancer in Erythroleukemia.

GATA-1 and PU.1 are two important hematopoietic transcription factors that mutually inhibit each other in progenitor cells to guide entrance into the erythroid or myeloid lineage, respectively. PU.1 controls its own expression during myelopoiesis by binding to the distal URE enhancer, whose deletion leads to acute myeloid leukemia (AML). We herein present evidence that GATA-1 binds to the PU.1 gene and inhibits its expression in human AML-erythroleukemias (EL). Furthermore, GATA-1 together with DNA methyl Transferase I (DNMT1) mediate repression of the PU.1 gene through the URE. Repression of the PU.1 gene involves both DNA methylation at the URE and its histone H3 lysine-K9 methylation and deacetylation as well as the H3K27 methylation at additional DNA elements and the promoter. The GATA-1-mediated inhibition of PU.1 gene transcription in human AML-EL mediated through the URE represents important mechanism that contributes to PU.1 downregulation and leukemogenesis that is sensitive to DNA demethylation therapy.

GATA1 and PU.1 Bind to Ribosomal Protein Genes in Erythroid Cells: Implications for Ribosomopathies.

The clear connection between ribosome biogenesis dysfunction and specific hematopoiesis-related disorders prompted us to examine the role of critical lineage-specific transcription factors in the transcriptional regulation of ribosomal protein (RP) genes during terminal erythroid differentiation. By applying EMSA and ChIP methodologies in mouse erythroleukemia cells we show that GATA1 and PU.1 bind in vitro and in vivo the proximal promoter region of the RPS19 gene which is frequently mutated in Diamond-Blackfan Anemia. Moreover, ChIPseq data analysis also demonstrates that several RP genes are enriched as potential GATA1 and PU.1 gene targets in mouse and human erythroid cells, with GATA1 binding showing an association with higher ribosomal protein gene expression levels during terminal erythroid differentiation in human and mouse. Our results suggest that RP gene expression and hence balanced ribosome biosynthesis may be specifically and selectively regulated by lineage specific transcription factors during hematopoiesis, a finding which may be clinically relevant to ribosomopathies.

Genomic and proteomic analysis of transcription factor TFII-I reveals insight into the response to cellular stress.

The ubiquitously expressed transcription factor TFII-I exerts both positive and negative effects on transcription. Using biotinylation tagging technology and high-throughput sequencing, we determined sites of chromatin interactions for TFII-I in the human erythroleukemia cell line K562. This analysis revealed that TFII-I binds upstream of the transcription start site of expressed genes, both upstream and downstream of the transcription start site of repressed genes, and downstream of RNA polymerase II peaks at the ATF3 and other stress responsive genes. At the ATF3 gene, TFII-I binds immediately downstream of a Pol II peak located 5 kb upstream of exon 1. Induction of ATF3 expression increases transcription throughout the ATF3 gene locus which requires TFII-I and correlates with increased association of Pol II and Elongin A. Pull-down assays demonstrated that TFII-I interacts with Elongin A. Partial depletion of TFII-I expression caused a reduction in the association of Elongin A with and transcription of the DNMT1 and EFR3A genes without a decrease in Pol II recruitment. The data reveal different interaction patterns of TFII-I at active, repressed, or inducible genes, identify novel TFII-I interacting proteins, implicate TFII-I in the regulation of transcription elongation and provide insight into the role of TFII-I during the response to cellular stress.

Metal and metallothionein concentrations in Paracentrotus lividus from Amvrakikos gulf (Ionian Sea-Greece).

Concentrations of Cd, Cu, Cr, Ni, Zn, Fe and metallothioneins (MTs) were measured in the gonads of Paracentrotus lividus from Amvrakikos gulf (Ionian Sea, Greece). Three natural populations were selected; two of them, growing inside the gulf (Agios Thomas and Koronisia), presented higher density and smaller body size than the population living in a coastal area just outside the gulf (Mytikas). Metal and MT levels were not elevated, with the exception of Zn, showing high values related to the reproduction stage of the sea urchins. Significant geographical variations were recorded in the concentrations of Cu, Zn, Cd, Cr and MTs. The highest mean and maximum values of Cu, Zn and MTs were recorded in Agios Thomas while Cd and Cr were higher in Mytikas population. Copper, Zn, Fe and MT concentrations were negatively correlated to the sea urchin body size, while a positive concentration-size relationship was observed for Cd. Although all studied populations grow in a low metal level marine environment, urchins with smaller body size living in a food limited marine environment showed higher gonadosomatic index, metal concentrations and MT levels in their gonads (Agios Thomas and Koronisia) than larger specimens growing in a food unlimited area (Mytikas).

Differential microRNA profiles and their functional implications in different immunogenetic subsets of chronic lymphocytic leukemia.

Critical processes of B-cell physiology, including immune signaling through the B-cell receptor (BcR) and/or Toll-like receptors (TLRs), are targeted by microRNAs. With this in mind and also given the important role of BcR and TLR signaling and microRNAs in chronic lymphocytic leukemia (CLL), we investigated whether microRNAs could be implicated in shaping the behavior of CLL clones with distinct BcR and TLR molecular and functional profiles. To this end, we examined 79 CLL cases for the expression of 33 microRNAs, selected on the following criteria: (a) deregulated in CLL versus normal B-cells; (b) differentially expressed in CLL subgroups with distinct clinicobiological features; and, (c) if meeting (a) + (b), having predicted targets in the immune signaling pathways. Significant upregulation of miR-150, miR-29c, miR-143 and miR-223 and downregulation of miR-15a was found in mutated versus unmutated CLL, with miR-15a showing the highest fold difference. Comparison of two major subsets with distinct stereotyped BcRs and signaling signatures, namely subset 1 [IGHV1/5/7-IGKV1(D)-39, unmutated, bad prognosis] versus subset 4 [IGHV4-34/IGKV2-30, mutated, good prognosis] revealed differences in the expression of miR-150, miR-29b, miR-29c and miR-101, all down-regulated in subset 1. We were also able to link these distinct microRNA profiles with cellular phenotypes, importantly showing that, in subset 1, miR-101 downregulation is associated with overexpression of the enhancer of zeste homolog 2 (EZH2) protein, which has been associated with clinical aggressiveness in other B-cell lymphomas. In conclusion, specific miRNAs differentially expressed among CLL subgroups with distinct BcR and/or TLR signaling may modulate the biological and clinical behavior of the CLL clones.

GATA-1 genome-wide occupancy associates with distinct epigenetic profiles in mouse fetal liver erythropoiesis.

We report the genomic occupancy profiles of the key hematopoietic transcription factor GATA-1 in pro-erythroblasts and mature erythroid cells fractionated from day E12.5 mouse fetal liver cells. Integration of GATA-1 occupancy profiles with available genome-wide transcription factor and epigenetic profiles assayed in fetal liver cells enabled as to evaluate GATA-1 involvement in modulating local chromatin structure of target genes during erythroid differentiation. Our results suggest that GATA-1 associates preferentially with changes of specific epigenetic modifications, such as H4K16, H3K27 acetylation and H3K4 di-methylation. Furthermore, we used random forest (RF) non-linear regression to predict changes in the expression levels of GATA-1 target genes based on the genomic features available for pro-erythroblasts and mature fetal liver-derived erythroid cells. Remarkably, our prediction model explained a high proportion of 62% of variation in gene expression. Hierarchical clustering of the proximity values calculated by the RF model produced a clear separation of upregulated versus downregulated genes and a further separation of downregulated genes in two distinct groups. Thus, our study of GATA-1 genome-wide occupancy profiles in mouse primary erythroid cells and their integration with global epigenetic marks reveals three clusters of GATA-1 gene targets that are associated with specific epigenetic signatures and functional characteristics.

Accurate microRNA Target Prediction Using Detailed Binding Site Accessibility and Machine Learning on Proteomics Data.

MicroRNAs (miRNAs) are a class of small regulatory genes regulating gene expression by targeting messenger RNA. Though computational methods for miRNA target prediction are the prevailing means to analyze their function, they still miss a large fraction of the targeted genes and additionally predict a large number of false positives. Here we introduce a novel algorithm called DIANA-microT-ANN which combines multiple novel target site features through an artificial neural network (ANN) and is trained using recently published high-throughput data measuring the change of protein levels after miRNA overexpression, providing positive and negative targeting examples. The features characterizing each miRNA recognition element include binding structure, conservation level, and a specific profile of structural accessibility. The ANN is trained to integrate the features of each recognition element along the 3'untranslated region into a targeting score, reproducing the relative repression fold change of the protein. Tested on two different sets the algorithm outperforms other widely used algorithms and also predicts a significant number of unique and reliable targets not predicted by the other methods. For 542 human miRNAs DIANA-microT-ANN predicts 120000 targets not provided by TargetScan 5.0. The algorithm is freely available at http://microrna.gr/microT-ANN.

The DIANA-mirExTra web server: from gene expression data to microRNA function.

BACKGROUND: High-throughput gene expression experiments are widely used to identify the role of genes involved in biological conditions of interest. MicroRNAs (miRNA) are regulatory molecules that have been functionally associated with several developmental programs and their deregulation with diverse diseases including cancer. METHODOLOGY/PRINCIPAL FINDINGS: Although miRNA expression levels may not be routinely measured in high-throughput experiments, a possible involvement of miRNAs in the deregulation of gene expression can be computationally predicted and quantified through analysis of overrepresented motifs in the deregulated genes 3' untranslated region (3'UTR) sequences. Here, we introduce a user-friendly web-server, DIANA-mirExTra (www.microrna.gr/mirextra) that allows the comparison of frequencies of miRNA associated motifs between sets of genes that can lead to the identification of miRNAs responsible for the deregulation of large numbers of genes. To this end, we have investigated different approaches and measures, and have practically implemented them on experimental data. CONCLUSIONS/SIGNIFICANCE: On several datasets of miRNA overexpression and repression experiments, our proposed approaches have successfully identified the deregulated miRNA. Beyond the prediction of miRNAs responsible for the deregulation of transcripts, the web-server provides extensive links to DIANA-mirPath, a functional analysis tool incorporating miRNA targets in biological pathways. Additionally, in case information about miRNA expression changes is provided, the results can be filtered to display the analysis for miRNAs of interest only.

Accurate microRNA target prediction correlates with protein repression levels.

BACKGROUND: MicroRNAs are small endogenously expressed non-coding RNA molecules that regulate target gene expression through translation repression or messenger RNA degradation. MicroRNA regulation is performed through pairing of the microRNA to sites in the messenger RNA of protein coding genes. Since experimental identification of miRNA target genes poses difficulties, computational microRNA target prediction is one of the key means in deciphering the role of microRNAs in development and disease. RESULTS: DIANA-microT 3.0 is an algorithm for microRNA target prediction which is based on several parameters calculated individually for each microRNA and combines conserved and non-conserved microRNA recognition elements into a final prediction score, which correlates with protein production fold change. Specifically, for each predicted interaction the program reports a signal to noise ratio and a precision score which can be used as an indication of the false positive rate of the prediction. CONCLUSION: Recently, several computational target prediction programs were benchmarked based on a set of microRNA target genes identified by the pSILAC method. In this assessment DIANA-microT 3.0 was found to achieve the highest precision among the most widely used microRNA target prediction programs reaching approximately 66%. The DIANA-microT 3.0 prediction results are available online in a user friendly web server at http://www.microrna.gr/microT.

Mung bean nuclease mapping of RNAs 3' end.

A method is described that allows an accurate mapping of 3' ends of RNAs. In this method a labeled DNA probe, containing the presumed 3' end of the RNA under analysis is allowed to anneals to the RNA itself. Mung-bean nuclease is then used to digest single strands of both RNA and DNA. Electrophoretic fractionation of "protected" undigested, labeled DNA is than performed using a sequence reaction of a known DNA as length marker. This procedure was applied to the analysis of both a polyA RNA (Interleukin 10 mRNA) and non polyA RNAs (sea urchin 18S and 26S rRNAs). This method might be potentially relevant for the evaluation of the role of posttrascriptional control of IL-10 in the pathogenesis of the immune and inflammatory mediated diseases associated to ageing. This might allow to develop new strategies to approach to the diagnosis and therapy of age related diseases.