regioneReloaded: evaluating the association of multiple genomic region sets.

MOTIVATION: Next-generation sequencing methods continue improving the annotation of genomes in part by determining the distribution of features such as epigenetic marks. Evaluating and interpreting the association between genomic regions and their features has become a common and challenging analysis in genomic and epigenomic studies. RESULTS: With regioneR we provided an R package allowing to assess the statistical significance of pairwise associations between genomic region sets using permutation tests. We now present the R package regioneReloaded that builds upon regioneR's statistical foundation and extends the functionality for the simultaneous analysis and visualization of the associations between multiple genomic region sets. Thus, we provide a novel discovery tool for the identification of significant associations that warrant to be tested for functional interdependence. AVAILABILITY AND IMPLEMENTATION: regioneReloaded is an R package released under an Artistic-2.0 License. The source code and documentation are freely available through Bioconductor: http://www.bioconductor.org/packages/regioneReloaded.

MacroH2A histone variants maintain nuclear organization and heterochromatin architecture.

Genetic loss-of-function studies on development, cancer and somatic cell reprogramming have suggested that the group of macroH2A histone variants might function through stabilizing the differentiated state by a yet unknown mechanism. Here, we present results demonstrating that macroH2A variants have a major function in maintaining nuclear organization and heterochromatin architecture. Specifically, we find that a substantial amount of macroH2A is associated with heterochromatic repeat sequences. We further identify macroH2A on sites of interstitial heterochromatin decorated by histone H3 trimethylated on K9 (H3K9me3). Loss of macroH2A leads to major defects in nuclear organization, including reduced nuclear circularity, disruption of nucleoli and a global loss of dense heterochromatin. Domains formed by DNA repeat sequences are disorganized, expanded and fragmented, and mildly re-expressed when depleted of macroH2A. At the molecular level, we find that macroH2A is required for the interaction of repeat sequences with the nucleostructural protein lamin B1. Taken together, our results argue that a major function of macroH2A histone variants is to link nucleosome composition to higher-order chromatin architecture.

Copy number profiling of adult relapsed B-cell precursor acute lymphoblastic leukemia reveals potential leukemia progression mechanisms.

The outcome of relapsed adult acute lymphoblastic leukemia (ALL) remains dismal despite new therapeutic approaches. Previous studies analyzing relapse samples have shown a high degree of heterogeneity regarding gene alterations without an evident relapse signature. Bone marrow or peripheral blood samples from 31 adult B-cell precursor ALL patients at first relapse, and 21 paired diagnostic samples were analyzed by multiplex ligation probe-dependent amplification (MLPA). Nineteen paired diagnostic and relapse samples of these 21 patients were also analyzed by SNP arrays. A trend to acquire homozygous CDKN2A/B deletions and a significant increase in the number of copy number alterations (CNA) was observed from diagnosis to first relapse. Evolution from an ancestral clone was the main pattern of clonal evolution. Relapse samples were extremely heterogeneous regarding CNA frequencies. However, CDKN2A/B, PAX5, ETV6, ATM, IKZF1, VPREB1, and TP53 deletions and duplications of 1q, 8q, 17q, 21, X/Y PAR1, and Xp were frequently detected at relapse. Duplications of genes involved in cell proliferation, drug resistance and stem cell homeostasis regulation, as well as deletions of KDM6A and STAG2 genes emerged as specific alterations at relapse. Genomics of relapsed adult B-cell precursor ALL is highly heterogeneous, although some recurrent lesions involved in essential pathways deregulation were frequently observed. Selective and simultaneous targeting of these deregulated pathways may improve the results of current salvage therapies.

regioneR: an R/Bioconductor package for the association analysis of genomic regions based on permutation tests.

MOTIVATION: Statistically assessing the relation between a set of genomic regions and other genomic features is a common challenging task in genomic and epigenomic analyses. Randomization based approaches implicitly take into account the complexity of the genome without the need of assuming an underlying statistical model. SUMMARY: regioneR is an R package that implements a permutation test framework specifically designed to work with genomic regions. In addition to the predefined randomization and evaluation strategies, regioneR is fully customizable allowing the use of custom strategies to adapt it to specific questions. Finally, it also implements a novel function to evaluate the local specificity of the detected association. AVAILABILITY AND IMPLEMENTATION: regioneR is an R package released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor (http://www.bioconductor.org/packages/regioneR). CONTACT: rmalinverni@carrerasresearch.org.

A Cbx8-containing polycomb complex facilitates the transition to gene activation during ES cell differentiation.

Polycomb proteins play an essential role in maintaining the repression of developmental genes in self-renewing embryonic stem cells. The exact mechanism allowing the derepression of polycomb target genes during cell differentiation remains unclear. Our project aimed to identify Cbx8 binding sites in differentiating mouse embryonic stem cells. Therefore, we used a genome-wide chromatin immunoprecipitation of endogenous Cbx8 coupled to direct massive parallel sequencing (ChIP-Seq). Our analysis identified 171 high confidence peaks. By crossing our data with previously published microarray analysis, we show that several differentiation genes transiently recruit Cbx8 during their early activation. Depletion of Cbx8 partially impairs the transcriptional activation of these genes. Both interaction analysis, as well as chromatin immunoprecipitation experiments support the idea that activating Cbx8 acts in the context of an intact PRC1 complex. Prolonged gene activation results in eviction of PRC1 despite persisting H3K27me3 and H2A ubiquitination. The composition of PRC1 is highly modular and changes when embryonic stem cells commit to differentiation. We further demonstrate that the exchange of Cbx7 for Cbx8 is required for the effective activation of differentiation genes. Taken together, our results establish a function for a Cbx8-containing complex in facilitating the transition from a Polycomb-repressed chromatin state to an active state. As this affects several key regulatory differentiation genes this mechanism is likely to contribute to the robust execution of differentiation programs.

A cellular model reflecting the phenotypic heterogeneity of mutant HRAS driven squamous cell carcinoma.

Squamous cell carcinomas have a range of histopathological manifestations. The parameters that determine this clinically observed heterogeneity are not fully understood. Here, we report the generation of a cell culture model that reflects part of this heterogeneity. We have used the catalytic subunit of human telomerase hTERT and large T to immortalize primary UV-unexposed keratinocytes. Then, mutant HRAS G12V has been introduced to transform these immortal keratinocytes. When injected into immunosuppressed mice, transformed cells grew as xenografts with distinct histopathological characteristics. We observed three major tissue architectures: solid, sarcomatoid and cystic growth types, which were primarily composed of pleomorphic and basaloid cells but in some cases displayed focal apocrine differentiation. We demonstrate that the cells generated represent different stages of skin cancerogenesis and as such can be used to identify novel tumor-promoting alterations such as the overexpression of the PADI2 oncogene in solid-type SCC. Importantly, the cultured cells maintain the characteristics from the xenograft they were derived from while being amenable to manipulation and analysis. The availability of cell lines representing different clinical manifestations opens a new tool to study the stochastic and deterministic factors that cause case-to-case heterogeneity despite departing from the same set of oncogenes and the same genetic background.

3D chromatin remodelling in the germ line modulates genome evolutionary plasticity.

Chromosome folding has profound impacts on gene regulation, whose evolutionary consequences are far from being understood. Here we explore the relationship between 3D chromatin remodelling in mouse germ cells and evolutionary changes in genome structure. Using a comprehensive integrative computational analysis, we (i) reconstruct seven ancestral rodent genomes analysing whole-genome sequences of 14 species representatives of the major phylogroups, (ii) detect lineage-specific chromosome rearrangements and (iii) identify the dynamics of the structural and epigenetic properties of evolutionary breakpoint regions (EBRs) throughout mouse spermatogenesis. Our results show that EBRs are devoid of programmed meiotic DNA double-strand breaks (DSBs) and meiotic cohesins in primary spermatocytes, but are associated in post-meiotic cells with sites of DNA damage and functional long-range interaction regions that recapitulate ancestral chromosomal configurations. Overall, we propose a model that integrates evolutionary genome reshuffling with DNA damage response mechanisms and the dynamic spatial genome organisation of germ cells.

MacroH2As regulate enhancer-promoter contacts affecting enhancer activity and sensitivity to inflammatory cytokines.

MacroH2A histone variants have a function in gene regulation that is poorly understood at the molecular level. We report that macroH2A1.2 and macroH2A2 modulate the transcriptional ground state of cancer cells and how they respond to inflammatory cytokines. Removal of macroH2A1.2 and macroH2A2 in hepatoblastoma cells affects the contact frequency of promoters and distal enhancers coinciding with changes in enhancer activity or preceding them in response to the cytokine tumor necrosis factor alpha. Although macroH2As regulate genes in both directions, they globally facilitate the nuclear factor κB (NF-κB)-mediated response. In contrast, macroH2As suppress the response to the pro-inflammatory cytokine interferon gamma. MacroH2A2 has a stronger contribution to gene repression than macroH2A1.2. Taken together, our results suggest that macroH2As have a role in regulating the response of cancer cells to inflammatory signals on the level of chromatin structure. This is likely relevant for the interaction of cancer cells with immune cells of their microenvironment.

Evolution of a histone variant involved in compartmental regulation of NAD metabolism.

NAD metabolism is essential for all forms of life. Compartmental regulation of NAD+ consumption, especially between the nucleus and the mitochondria, is required for energy homeostasis. However, how compartmental regulation evolved remains unclear. In the present study, we investigated the evolution of the macrodomain-containing histone variant macroH2A1.1, an integral chromatin component that limits nuclear NAD+ consumption by inhibiting poly(ADP-ribose) polymerase 1 in vertebrate cells. We found that macroH2A originated in premetazoan protists. The crystal structure of the macroH2A macrodomain from the protist Capsaspora owczarzaki allowed us to identify highly conserved principles of ligand binding and pinpoint key residue substitutions, selected for during the evolution of the vertebrate stem lineage. Metabolic characterization of the Capsaspora lifecycle suggested that the metabolic function of macroH2A was associated with nonproliferative stages. Taken together, we provide insight into the evolution of a chromatin element involved in compartmental NAD regulation, relevant for understanding its metabolism and potential therapeutic applications.

The Histone Variant MacroH2A1 Regulates Key Genes for Myogenic Cell Fusion in a Splice-Isoform Dependent Manner.

MacroH2A histone variants have functions in differentiation, somatic cell reprogramming and cancer. However, at present, it is not clear how macroH2As affect gene regulation to exert these functions. We have parted from the initial observation that loss of total macroH2A1 led to a change in the morphology of murine myotubes differentiated ex vivo. The fusion of myoblasts to myotubes is a key process in embryonic myogenesis and highly relevant for muscle regeneration after acute or chronic injury. We have focused on this physiological process, to investigate the functions of the two splice isoforms of macroH2A1. Individual perturbation of the two isoforms in myotubes forming in vitro from myogenic C2C12 cells showed an opposing phenotype, with macroH2A1.1 enhancing, and macroH2A1.2 reducing, fusion. Differential regulation of a subset of fusion-related genes encoding components of the extracellular matrix and cell surface receptors for adhesion correlated with these phenotypes. We describe, for the first time, splice isoform-specific phenotypes for the histone variant macroH2A1 in a physiologic process and provide evidence for a novel underlying molecular mechanism of gene regulation.

DNA methylation profile in chronic myelomonocytic leukemia associates with distinct clinical, biological and genetic features.

Chromosomal abnormalities are detected in 20-30% of patients with chronic myelomonocytic leukemia (CMML) and correlate with prognosis. On the mutation level, disruptive alterations are particularly frequent in chromatin regulatory genes. However, little is known about the consequential alterations in the epigenetic marking of the genome. Here, we report the analysis of genomic DNA methylation patterns of 64 CMML patients and 10 healthy controls, using a DNA methylation microarray focused on promoter regions. Differential methylation analysis between patients and controls allowed us to identify abnormalities in DNA methylation, including hypermethylation of specific genes and large genome regions with aberrant DNA methylation. Unsupervised hierarchical cluster analysis identified two main clusters that associated with the clinical, biological, and genetic features of patients. Group 1 was enriched in patients with adverse clinical and biological characteristics and poorer overall and progression-free survival. In addition, significant differences in DNA methylation were observed between patients with low risk and intermediate/high risk karyotypes and between TET2 mutant and wild type patients. Taken together, our results demonstrate that altered DNA methylation patterns reflect the CMML disease state and allow to identify patient groups with distinct clinical features.

Direct modulation of the bone marrow mesenchymal stromal cell compartment by azacitidine enhances healthy hematopoiesis.

Mesenchymal stromal cells (MSCs) are crucial components of the bone marrow (BM) microenvironment essential for regulating self-renewal, survival, and differentiation of hematopoietic stem/progenitor cells (HSPCs) in the stem cell niche. MSCs are functionally altered in myelodysplastic syndromes (MDS) and exhibit an altered methylome compared with MSCs from healthy controls, thus contributing to disease progression. To determine whether MSCs are amenable to epigenetic therapy and if this affects their function, we examined growth, differentiation, and HSPC-supporting capacity of ex vivo-expanded MSCs from MDS patients in comparison with age-matched healthy controls after direct treatment in vitro with the hypomethylating agent azacitidine (AZA). Strikingly, we find that AZA exerts a direct effect on healthy as well as MDS-derived MSCs such that they favor support of healthy over malignant clonal HSPC expansion in coculture experiments. RNA-sequencing analyses of MSCs identified stromal networks regulated by AZA. Notably, these comprise distinct molecular pathways crucial for HSPC support, foremost extracellular matrix molecules (including collagens) and interferon pathway components. Our study demonstrates that the hypomethylating agent AZA exerts its antileukemic activity in part through a direct effect on the HSPC-supporting BM niche and provides proof of concept for the therapeutic potential of epigenetic treatment of diseased MSCs. In addition, our comprehensive data set of AZA-sensitive gene networks represents a valuable framework to guide future development of targeted epigenetic niche therapy in myeloid malignancies such as MDS and acute myeloid leukemia.

Polycomb protein RING1A limits hematopoietic differentiation in myelodysplastic syndromes.

Genetic lesions affecting epigenetic regulators are frequent in myelodysplastic syndromes (MDS). Polycomb proteins are key epigenetic regulators of differentiation and stemness that act as two multimeric complexes termed polycomb repressive complexes 1 and 2, PRC1 and PRC2, respectively. While components and regulators of PRC2 such as ASXL1 and EZH2 are frequently mutated in MDS and AML, little is known about the role of PRC1. To analyze the role of PRC1, we have taken a functional approach testing PRC1 components in loss- and gain-of-function experiments that we found overexpressed in advanced MDS patients or dynamically expressed during normal hematopoiesis. This approach allowed us to identify the enzymatically active component RING1A as the key PRC1 component in hematopoietic stem cells and MDS. Specifically, we found that RING1A is expressed in CD34+ bone marrow progenitor cells and further overexpressed in high-risk MDS patients. Knockdown of RING1A in an MDS-derived AML cell line facilitated spontaneous and retinoic acid-induced differentiation. Similarly, inactivation of RING1A in primary CD34+ cells augmented erythroid differentiation. Treatment with a small compound RING1 inhibitor reduced the colony forming capacity of CD34+ cells from MDS patients and healthy controls. In MDS patients higher RING1A expression associated with an increased number of dysplastic lineages and blasts. Our data suggests that RING1A is deregulated in MDS and plays a role in the erythroid development defect.

MacroH2A1.1 regulates mitochondrial respiration by limiting nuclear NAD+ consumption.

Histone variants are structural components of eukaryotic chromatin that can replace replication-coupled histones in the nucleosome. The histone variant macroH2A1.1 contains a macrodomain capable of binding NAD+-derived metabolites. Here we report that macroH2A1.1 is rapidly induced during myogenic differentiation through a switch in alternative splicing, and that myotubes that lack macroH2A1.1 have a defect in mitochondrial respiratory capacity. We found that the metabolite-binding macrodomain was essential for sustained optimal mitochondrial function but dispensable for gene regulation. Through direct binding, macroH2A1.1 inhibits basal poly-ADP ribose polymerase 1 (PARP-1) activity and thus reduces nuclear NAD+ consumption. The resultant accumulation of the NAD+ precursor NMN allows for maintenance of mitochondrial NAD+ pools that are critical for respiration. Our data indicate that macroH2A1.1-containing chromatin regulates mitochondrial respiration by limiting nuclear NAD+ consumption and establishing a buffer of NAD+ precursors in differentiated cells.

Barcelona conference on epigenetics and cancer 2015: Coding and non-coding functions of the genome.

The Barcelona Conference on Epigenetics and Cancer (BCEC) entitled "Coding and Non-Coding functions of the Genome" took place October 29-30, 2015 in Barcelona. The 2015 BCEC was the third edition of a series of annual conferences jointly organized by 5 leading research centers in Barcelona together with B-Debate, an initiative of BioCat. Luciano Di Croce from the Center for Genomic Regulation and Marcus Buschbeck from the Josep Carreras Leukemia Research Institute put together the scientific program with a particular focus on the role of non-coding RNAs in enhancer regulation, epigenetic control by Polycomb complexes, histone variants, and nuclear organization. In one and a half days, 22 talks and 56 posters were presented to an audience of 215 participants.

Mammalian Comparative Genomics Reveals Genetic and Epigenetic Features Associated with Genome Reshuffling in Rodentia.

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.

A comparative gene index for the white sturgeon Acipenser transmontanus.

Sturgeons are archaic fishes phylogenetically distinct from Teleosts. They represent an important niche for aquaculture, particularly for the production of caviar and high quality fillets, while many natural populations in various world areas are today threatened by extinction. Knowledge of the sturgeon genome is limited, as it is the case of many other species of interest for fishery, aquaculture and conservation. Sequences from non-normalized libraries of skin and spleen of the American sturgeon (A. transmontanus) produced in our laboratories were analysed via a bioinformatic procedure, and compared to similar resources available for three Teleost species. Data collected during the analyses were stored in a database - the Sturgeon database (db) - that can be queried via a web interface. The Sturgeon db contains a total of 16,404 sequences from Acipenser transmontanus, Ictalurus punctatus, Salmo salar and Takifugu rubripes, each specie being represented by expressed sequence tags (ESTs) from skin and spleen. Data contained in the database are the results of a number of analyses that mostly focus on sequence annotation and intra- and inter-species comparison. Putative SNP sites, tandem repeats, and sequences matching known protein patterns and motifs were also identified. The Sturgeon db is by now the only online resource dedicated to the analysis of A. transmontanus EST sequences, and represents a starting point for the investigation of the genome of sturgeons from a physiological perspective; it will be used to identify polymorphic markers to study, for example, fish pathologies or to survey fish disease resistance, and to produce gene expression arrays. Introduction of sequences from other species in the analysis pipeline allowed inter-species comparisons of transcripts distribution in Gene Ontology categories, as well as orthologs identification, despite the high sturgeon phylogenetic distance from other fish species. As a result of the EST analysis procedure, 1058 sturgeon novel unigenes were identified.

Genome-wide transcriptional response of primary alveolar macrophages following infection with porcine reproductive and respiratory syndrome virus.

Porcine reproductive and respiratory syndrome is a major cause of economic loss for the swine industry worldwide. Porcine reproductive and respiratory syndrome virus (PRRSV) triggers weak and atypical innate immune responses, but key genes and mechanisms by which the virus interferes with the host innate immunity have not yet been elucidated. In this study, genes that control the response of the main target of PRRSV, porcine alveolar macrophages (PAMs), were profiled in vitro with a time-course experiment spanning the first round of virus replication. PAMs were obtained from six piglets and challenged with the Lelystad PRRSV strain, and gene expression was investigated using Affymetrix microarrays and real-time PCR. Of the 1409 differentially expressed transcripts identified by analysis of variance, two, five, 25, 16 and 100 differed from controls by a minimum of 1.5-fold at 1, 3, 6, 9 and 12 h post-infection (p.i.), respectively. A PRRSV infection effect was detectable between 3 and 6 h p.i., and was characterized by a consistent downregulation of gene expression, followed by the start of the host innate immune response at 9 h p.i. The expression of beta interferon 1 (IFN-beta), but not of IFN-alpha, was strongly upregulated, whilst few genes commonly expressed in response to viral infections and/or induced by interferons were found to be differentially expressed. A predominance of anti-apoptotic transcripts (e.g. interleukin-10), a shift towards a T-helper cell type 2 response and a weak upregulation of tumour necrosis factor-alpha expression were observed within 12 h p.i., reinforcing the hypotheses that PRRSV has developed sophisticated mechanisms to escape the host defence.

Gene expression profiling of porcine alveolar macrophages after antibody-mediated cross-linking of Sialoadhesin (Sn, Siglec-1).

Sialoadhesin (Sn) is the prototypic member of the Siglecs, a family of receptors mainly involved in cell-cell interactions. For several Siglecs, but not for Sn, intracellular signaling functions have been described. Because antibody-mediated cross-linking of surface transmembrane proteins is a powerful technique to investigate cell-molecular events, Sn expressed on porcine alveolar macrophages (PAM) was cross-linked with the antibody 41D3, and the expression profiles were compared with mock-treated macrophages by microarray analysis. Gene ontology analysis of 479 differentially expressed transcripts identified gene categories related to membrane localization, signal transduction, receptor and communication activities. Analyses of the human KEGG pathway database identified MAP kinase signaling, regulation of actin cytoskeleton, adipocytokine signaling, and wnt signaling as significantly altered pathways, supporting a role for Sn as intracellular signaling molecule. Real-time PCR of a subset of modulated genes confirmed these results and highlighted the reliability of a short-term cross-linking treatment for transcriptomic analysis of receptor functions.

The EADGENE Microarray Data Analysis Workshop (open access publication).

Microarray analyses have become an important tool in animal genomics. While their use is becoming widespread, there is still a lot of ongoing research regarding the analysis of microarray data. In the context of a European Network of Excellence, 31 researchers representing 14 research groups from 10 countries performed and discussed the statistical analyses of real and simulated 2-colour microarray data that were distributed among participants. The real data consisted of 48 microarrays from a disease challenge experiment in dairy cattle, while the simulated data consisted of 10 microarrays from a direct comparison of two treatments (dye-balanced). While there was broader agreement with regards to methods of microarray normalisation and significance testing, there were major differences with regards to quality control. The quality control approaches varied from none, through using statistical weights, to omitting a large number of spots or omitting entire slides. Surprisingly, these very different approaches gave quite similar results when applied to the simulated data, although not all participating groups analysed both real and simulated data. The workshop was very successful in facilitating interaction between scientists with a diverse background but a common interest in microarray analyses.