Histone variant macroH2A marks embryonic differentiation in vivo and acts as an epigenetic barrier to induced pluripotency

Multifunctional histone variants in genome function

Histones are integral components of eukaryotic chromatin that have a pivotal role in the organization and function of the genome. The dynamic regulation of chromatin involves the incorporation of histone variants, which can dramatically alter its structural and functional properties. Contrary to an earlier view that limited individual histone variants to specific genomic functions, new insights have revealed that histone variants exert multifaceted roles involving all aspects of genome function, from governing patterns of gene expression at precise genomic loci to participating in genome replication, repair and maintenance. This conceptual change has led to a new understanding of the intricate interplay between chromatin and DNA-dependent processes and how this connection translates into normal and abnormal cellular functions.

The Function of H2A Histone Variants and Their Roles in Diseases

Epigenetic regulation, which is characterized by reversible and heritable genetic alterations without changing DNA sequences, has recently been increasingly studied in diseases. Histone variant regulation is an essential component of epigenetic regulation. The substitution of canonical histones by histone variants profoundly alters the local chromatin structure and modulates DNA accessibility to regulatory factors, thereby exerting a pivotal influence on gene regulation and DNA damage repair. Histone H2A variants, mainly including H2A.Z, H2A.B, macroH2A, and H2A.X, are the most abundant identified variants among all histone variants with the greatest sequence diversity. Harboring varied chromatin occupancy and structures, histone H2A variants perform distinct functions in gene transcription and DNA damage repair. They are implicated in multiple pathophysiological mechanisms and the emergence of different illnesses. Cancer, embryonic development abnormalities, neurological diseases, metabolic diseases, and heart diseases have all been linked to histone H2A variant alterations. This review focuses on the functions of H2A histone variants in mammals, including H2A.Z, H2A.B, macroH2A, and H2A.X, and their current roles in various diseases.

Reprogramming Cell Identity: Past Lessons, Challenges, and Future Directions

Reprogramming is traditionally defined as the fate conversion of a cell to a stage of increased developmental potential. In its broader meaning, the reprogramming term is also applied to all forms of cell fate conversion that do not follow a developmental trajectory. Reprogramming is now a well-established field of research that gained rapid progress upon the advent of induced pluripotency. In this perspective, I reflect on the reprogramming lessons of the past, in the contributions to other fields of research and on the potential transformative future use of reprogrammed cells and of its cell derivatives.

Casting histone variants during mammalian reproduction

During mammalian reproduction, germ cell chromatin packaging is key to prepare parental genomes for fertilization and to initiate embryonic development. While chromatin modifications such as DNA methylation and histone post-translational modifications are well known to carry regulatory information, histone variants have received less attention in this context. Histone variants alter the stability, structure and function of nucleosomes and, as such, contribute to chromatin organization in germ cells. Here, we review histone variants expression dynamics during the production of male and female germ cells, and what is currently known about their parent-of-origin effects during reproduction. Finally, we discuss the apparent conundrum behind these important functions and their recent evolutionary diversification.

macroH2A2 antagonizes epigenetic programs of stemness in glioblastoma

Self-renewal is a crucial property of glioblastoma cells that is enabled by the choreographed functions of chromatin regulators and transcription factors. Identifying targetable epigenetic mechanisms of self-renewal could therefore represent an important step toward developing effective treatments for this universally lethal cancer. Here we uncover an epigenetic axis of self-renewal mediated by the histone variant macroH2A2. With omics and functional assays deploying patient-derived in vitro and in vivo models, we show that macroH2A2 shapes chromatin accessibility at enhancer elements to antagonize transcriptional programs of self-renewal. macroH2A2 also sensitizes cells to small molecule-mediated cell death via activation of a viral mimicry response. Consistent with these results, our analyses of clinical cohorts indicate that high transcriptional levels of this histone variant are associated with better prognosis of high-grade glioma patients. Our results reveal a targetable epigenetic mechanism of self-renewal controlled by macroH2A2 and suggest additional treatment approaches for glioblastoma patients.

The dynamics of chromatin states mediated by epigenetic modifications during somatic cell reprogramming

Somatic cell reprogramming (SCR) is the conversion of differentiated somatic cells into totipotent or pluripotent cells through a variety of methods. Somatic cell reprogramming also provides a platform to investigate the role of chromatin-based factors in establishing and maintaining totipotency or pluripotency, since high expression of totipotency- or pluripotency-related genes usually require an active chromatin state. Several studies in plants or mammals have recently shed light on the molecular mechanisms by which epigenetic modifications regulate the expression of totipotency or pluripotency genes by altering their chromatin states. In this review, we present a comprehensive overview of the dynamic changes in epigenetic modifications and chromatin states during reprogramming from somatic cells to totipotent or pluripotent cells. In addition, we illustrate the potential role of DNA methylation, histone modifications, histone variants, and chromatin remodeling during somatic cell reprogramming, which will pave the way to developing reliable strategies for efficient cellular reprogramming.

Chromatin Structure from Development to Ageing

Nuclear structure influences genome architecture, which contributes to determine patterns of gene expression. Global changes in chromatin dynamics are essential during development and differentiation, and are one of the hallmarks of ageing. This chapter describes the molecular dynamics of chromatin structure that occur during development and ageing. In the first part, we introduce general information about the nuclear lamina, the chromatin structure, and the 3D organization of the genome. Next, we detail the molecular hallmarks found during development and ageing, including the role of DNA and histone modifications, 3D genome dynamics, and changes in the nuclear lamina. Within the chapter we discuss the implications that genome structure has on the mechanisms that drive development and ageing, and the physiological consequences when these mechanisms fail.

Cleavage-embryo genes and transposable elements are regulated by histone variant H2A.X

During mammalian preimplantation development, stimulation of zygotic genome activation (ZGA) and transposable elements (TEs) shapes totipotency profiling. A rare mouse embryonic stem cells (mESCs) subpopulation is capable of transiently entering a state resembling 2-cell stage embryos, with subtypes of TEs expressed and ZGA genes transiently activated. In this study, we found that deletion of H2A.X in mESCs led to a significant upregulation of ZGA genes and misregulated TEs. ChIP-seq analysis indicated a direct association of H2A.X at the Dux locus for silencing the Dux gene and its downstream ZGA genes in mESCs. We also demonstrated that histone variant H2A.X is highly enriched in human cleavage embryos when ZGA genes and TEs are active. Therefore, we propose that H2A.X plays an important role in regulating ZGA genes and TEs to establish totipotency.

Cellular senescence or stemness: hypoxia flips the coin

Cellular senescence is a complex physiological state whose main feature is proliferative arrest. Cellular senescence can be considered the reverse of cell immortalization and continuous tumor growth. However, cellular senescence has many physiological functions beyond being a putative tumor suppressive trait. It remains unknown whether low levels of oxygen or hypoxia, which is a feature of every tissue in the organism, modulate cellular senescence, altering its capacity to suppress the limitation of proliferation. It has been observed that the lifespan of mammalian primary cells is increased under low oxygen conditions. Additionally, hypoxia promotes self-renewal and pluripotency maintenance in adult and embryonic stem cells (ESCs), induced pluripotent stem cells (iPSCs) and cancer stem cells (CSCs). In this study, we discuss the role of hypoxia facilitating senescence bypass during malignant transformation and acquisition of stemness properties, which all contribute to tumor development and cancer disease aggressiveness.

Islet Epigenetic Impacts on β-Cell Identity and Function

The development and maintenance of differentiation is vital to the function of mature cells. Terminal differentiation is achieved by locking in the expression of genes essential for the function of those cells. Gene expression and its memory through generations of cell division is controlled by transcription factors and a host of epigenetic marks. In type 2 diabetes, β cells have altered gene expression compared to controls, accompanied by altered chromatin marks. Mutations, diet, and environment can all disrupt the implementation and preservation of the distinctive β-cell transcriptional signature. Understanding of the full complement of genomic control in β cells is still nascent. This article describes the known effects of histone marks and variants, DNA methylation, how they are regulated in the β cell, and how they affect cell-fate specification, maintenance, and lineage propagation. © 2021 American Physiological Society. Compr Physiol 11:1-18, 2021.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

The diagnostic and prognostic value of H2AFY in hepatocellular carcinoma

Background

The potential correlation between H2AFY (also known as MacroH2A1) and the clinical characteristics of hepatocellular carcinoma (HCC) patients was analysed through gene expression profiles and clinical data in The Cancer Genome Atlas (TCGA) database, and the diagnostic and prognostic value of H2AFY in HCC was discussed.

Methods

The gene expression data of HCC and the corresponding clinical characteristics of HCC patients were downloaded from the TCGA database. The differences in H2AFY in normal liver tissues and HCC were analysed. The relationship between H2AFY and clinical characteristics was analysed by Wilcoxon signed-rank test, logistic regression and Kruskal-Wallis test. The Kaplan-Meier method and the Cox regression method were used to analyse the relationship between overall survival and clinical characteristics of the patients. An ROC curve was used to predict the diagnostic value of H2AFY in HCC. Gene set enrichment analysis (GSEA) was used to analyse the pathway enrichment of H2AFY.

Result

Compared with normal liver tissues, H2AFY was significantly highly expressed in HCC. H2AFY was positively correlated with the age, clinical stage, G stage (grade) and T stage (tumor stage) of liver cancer patients. Higher H2AFY expression predicted a poor prognosis in HCC patients. Cox regression analysis suggested that H2AFY was an independent risk factor for the prognosis of HCC patients. The ROC curve suggested that H2AFY had certain diagnostic value in HCC. GSEA suggested that H2AFY was correlated with lipid metabolism and a variety of tumour pathways.

Conclusion

Our study showed that H2AFY was significantly overexpressed in HCC. H2AFY may be a potential diagnostic and prognostic marker for HCC, and high expression of H2AFY predicts a poor prognosis in patients with HCC.

Solid tumours hijack the histone variant network

Cancer is a complex disease characterized by loss of cellular homeostasis through genetic and epigenetic alterations. Emerging evidence highlights a role for histone variants and their dedicated chaperones in cancer initiation and progression. Histone variants are involved in processes as diverse as maintenance of genome integrity, nuclear architecture and cell identity. On a molecular level, histone variants add a layer of complexity to the dynamic regulation of transcription, DNA replication and repair, and mitotic chromosome segregation. Because these functions are critical to ensure normal proliferation and maintenance of cellular fate, cancer cells are defined by their capacity to subvert them. Hijacking histone variants and their chaperones is emerging as a common means to disrupt homeostasis across a wide range of cancers, particularly solid tumours. Here we discuss histone variants and histone chaperones as tumour-promoting or tumour-suppressive players in the pathogenesis of cancer.

Redox and Epigenetics in Human Pluripotent Stem Cells Differentiation

Significance: Since their discovery, induced pluripotent stem cells (iPSCs) had generated considerable interest in the scientific community for their great potential in regenerative medicine, disease modeling, and cell-based therapeutic approach, due to their unique characteristics of self-renewal and pluripotency. Recent Advances: Technological advances in iPSC genome-wide epigenetic profiling led to the elucidation of the epigenetic control of cellular identity during nuclear reprogramming. Moreover, iPSC physiology and metabolism are tightly regulated by oxidation-reduction events that mainly occur during the respiratory chain. In theory, iPSC-derived differentiated cells would be ideal for stem cell transplantation as autologous cells from donors, as the risks of rejection are minimal. Critical Issues: However, iPSCs experience high oxidative stress that, in turn, confers a high risk of increased genomic instability, which is most often linked to DNA repair deficiencies. Genomic instability has to be assessed before iPSCs can be used in therapeutic designs. Future Directions: This review will particularly focus on the links between redox balance and epigenetic modifications-in particular based on the histone variant macroH2A1-that determine DNA damage response in iPSCs and derived differentiated cells, and that might be exploited to decrease the teratogenic potential on iPSC transplantation. Antioxid. Redox Signal. 34, 335-349.

Histone chaperone APLF level dictates the implantation of mouse embryos

Our recent findings demonstrated that the histone chaperone and DNA repair factor aprataxin and PNK-like factor (APLF) could regulate epithelial to mesenchymal transition (EMT) during the reprogramming of murine fibroblasts and in breast cancer metastasis. Therefore, we investigated the function of APLF in EMT associated with mouse development. Here, we show that APLF is predominantly enhanced in trophectoderm (TE) and lineages derived from TE in pre- and post-implantation embryos. Downregulation of APLF induced the hatching of embryos in vitro, with a significant increase in Cdh1 and Cdx2 expression. Aplf short hairpin RNA-microinjected embryos failed to implant in vivo Rescue experiments neutralized the knockdown effects of APLF both in vitro and in vivo Reduced expression of Snai2 and Tead4, and the gain in Cdh1 and sFlt1 (also known as Flt1) level, marked the differentiation of APLF-knocked down trophoblast stem cells that might contribute towards the impaired implantation of embryos. Hence, our findings suggest a novel role for APLF during implantation and post-implantation development of mouse embryos. We anticipate that APLF might contribute to the establishment of maternal-fetal connection, as its fine balance is required to achieve implantation and thereby attain proper pregnancy.

LSH mediates gene repression through macroH2A deposition

The human Immunodeficiency Centromeric Instability Facial Anomalies (ICF) 4 syndrome is a severe disease with increased mortality caused by mutation in the LSH gene. Although LSH belongs to a family of chromatin remodeling proteins, it remains unknown how LSH mediates its function on chromatin in vivo. Here, we use chemical-induced proximity to rapidly recruit LSH to an engineered locus and find that LSH specifically induces macroH2A1.2 and macroH2A2 deposition in an ATP-dependent manner. Tethering of LSH induces transcriptional repression and silencing is dependent on macroH2A deposition. Loss of LSH decreases macroH2A enrichment at repeat sequences and results in transcriptional reactivation. Likewise, reduction of macroH2A by siRNA interference mimicks transcriptional reactivation. ChIP-seq analysis confirmed that LSH is a major regulator of genome-wide macroH2A distribution. Tethering of ICF4 mutations fails to induce macroH2A deposition and ICF4 patient cells display reduced macroH2A deposition and transcriptional reactivation supporting a pathogenic role for altered marcoH2A deposition. We propose that LSH is a major chromatin modulator of the histone variant macroH2A and that its ability to insert marcoH2A into chromatin and transcriptionally silence is disturbed in the ICF4 syndrome.

The human ICF 4 syndrome is caused by mutation of the chromatin remodeller LSH. Here, the authors show that LSH depletion disrupts the ability of histone variant macroH2A to insert into chromatin and silence transcription.

Histone Variants: Guardians of Genome Integrity

Chromatin integrity is key for cell homeostasis and for preventing pathological development. Alterations in core chromatin components, histone proteins, recently came into the spotlight through the discovery of their driving role in cancer. Building on these findings, in this review, we discuss how histone variants and their associated chaperones safeguard genome stability and protect against tumorigenesis. Accumulating evidence supports the contribution of histone variants and their chaperones to the maintenance of chromosomal integrity and to various steps of the DNA damage response, including damaged chromatin dynamics, DNA damage repair, and damage-dependent transcription regulation. We present our current knowledge on these topics and review recent advances in deciphering how alterations in histone variant sequence, expression, and deposition into chromatin fuel oncogenic transformation by impacting cell proliferation and cell fate transitions. We also highlight open questions and upcoming challenges in this rapidly growing field.

An Insight into Reprogramming Barriers to iPSC Generation

Derivation of induced Pluripotent Stem Cells (iPSCs) by reprogramming somatic cells to a pluripotent state has revolutionized stem cell research. Ensuing this, various groups have used genetic and non-genetic approaches to generate iPSCs from numerous cell types. However, achieving a pluripotent state in most of the reprogramming studies is marred by serious limitations such as low reprogramming efficiency and slow kinetics. These limitations are mainly due to the presence of potent barriers that exist during reprogramming when a mature cell is coaxed to achieve a pluripotent state. Several studies have revealed that intrinsic factors such as non-optimal stoichiometry of reprogramming factors, specific signaling pathways, cellular senescence, pluripotency-inhibiting transcription factors and microRNAs act as a roadblock. In addition, the epigenetic state of somatic cells and specific epigenetic modifications that occur during reprogramming also remarkably impede the generation of iPSCs. In this review, we present a comprehensive overview of the barriers that inhibit reprogramming and the understanding of which will pave the way to develop safe strategies for efficient reprogramming.

Histone variants in skeletal myogenesis

Histone variants regulate chromatin accessibility and gene transcription. Given their distinct properties and functions, histone varint substitutions allow for profound alteration of nucleosomal architecture and local chromatin landscape. Skeletal myogenesis driven by the key transcription factor MyoD is characterized by precise temporal regulation of myogenic genes. Timed substitution of variants within the nucleosomes provides a powerful means to ensure sequential expression of myogenic genes. Indeed, growing evidence has shown H3.3, H2A.Z, macroH2A, and H1b to be critical for skeletal myogenesis. However, the relative importance of various histone variants and their associated chaperones in myogenesis is not fully appreciated. In this review, we summarize the role that histone variants play in altering chromatin landscape to ensure proper muscle differentiation. The temporal regulation and cross talk between histones variants and their chaperones in conjunction with other forms of epigenetic regulation could be critical to understanding myogenesis and their involvement in myopathies.

Coaching from the sidelines: the nuclear periphery in genome regulation

The genome is packaged and organized nonrandomly within the 3D space of the nucleus to promote efficient gene expression and to faithfully maintain silencing of heterochromatin. The genome is enclosed within the nucleus by the nuclear envelope membrane, which contains a set of proteins that actively participate in chromatin organization and gene regulation. Technological advances are providing views of genome organization at unprecedented resolution and are beginning to reveal the ways that cells co-opt the structures of the nuclear periphery for nuclear organization and gene regulation. These genome regulatory roles of proteins of the nuclear periphery have important influences on development, disease and ageing.

Zebrafish macroH2A variants have distinct embryo localization and function

Mouse and cell-based studies have shown that macroH2A histone variants predominantly associate with heterochromatin. Functional studies found that macroH2As are involved in gene repression, inhibiting the acquisition of pluripotency and preserving cell differentiation. However, only a few studies have analysed the role of macroH2A during early embryo development. We report the development of transgenic zebrafish lines expressing macroH2A isoforms (mH2A1 and mH2A2) fusion proteins (with GFP) under identified endogenous promoters. We found that mH2A1 and mH2A2 have different spatial and temporal expression patterns during embryonic development. mH2A1 is expressed mostly in the extraembryonic Yolk Syncytial Layer (YSL) starting before shield stage and decreasing once morphogenesis is completed. mH2A2 expression lags behind mH2A1, becoming evident at 24 hpf, within the whole body of the embryo proper. Our ChIP-seq analysis showed that mH2A1 and mH2A2 bind to different DNA regions, changing dramatically after gastrulation. We further analysed RNA-seq data and showed that there is not a general/unspecific repressing function of mH2A1 or mH2A2 associated with heterochromatin but a fine regulation depending on cell types and stage of development. mH2A1 downregulates DNA expression in specific cells and embryo stages and its effect is independent of heterochromatin formation but it is correlated with nucleus quiescence instead. Whereas mH2A2 DNA association correlates with upregulation of differentially expressed genes between 75% epiboly and 24 hpf stages. Our data provide information for underlying molecules that participate in crucial early developmental events, and open new venues to explore mH2A related mechanisms that involve cell proliferation, differentiation, migration and metabolism.

Retargeting of macroH2A following mitosis to cytogenetic-scale heterochromatic domains

How macroH2A, a histone variant involved in silencing gene expression, is inherited from parent to daughter cells is unclear. Using a combination of imaging, biochemical, and genomic approaches, Sato et al. describe how newly synthesized macroH2A is incorporated predominantly in the G1 phase of human mitosis, targeting heterochromatic regions.

The heritability of chromatin states through cell division is a potential contributor to the epigenetic maintenance of cellular memory of prior states. The macroH2A histone variant has properties of a regulator of epigenetic cell memory, including roles controlling gene silencing and cell differentiation. Its mechanisms of regional genomic targeting and maintenance through cell division are unknown. Here, we combined in vivo imaging with biochemical and genomic approaches to show that human macroH2A is incorporated into chromatin in the G1 phase of the cell cycle following DNA replication. The newly incorporated macroH2A retargets the same large heterochromatic domains where macroH2A was already enriched in the previous cell cycle. It remains heterotypic, targeting individual nucleosomes that do not already contain a macroH2A molecule. The pattern observed resembles that of a new deposition of centromeric histone variants during the cell cycle, indicating mechanistic similarities for macrodomain-scale regulation of epigenetic properties of the cell.

Histone variant MacroH2A1 is downregulated in prostate cancer and influences malignant cell phenotype

BACKGROUND

Prostate cancer (PCa), a major cause of cancer-related morbidity and mortality worldwide and mostly asymptomatic at earliest stages, is characterized by disruption of genetic and epigenetic balance. A better understanding of how those mechanisms orchestrate disease might improve diagnostic and prognostic tools, allowing for improvements in treatment efficacy. Replacement of canonical histones, an epigenetic mechanism, is highly conserved among species and altered expression of histones variants (e.g., MacroH2A1) has been associated with tumorigenesis. H2AFY gene encodes two isoforms of H2A histone variant MacroH2A1: MacroH2A1.1 and MacroH2A1.2. Specifically, MacroH2A1.1 isoform inhibits cell proliferation and promotes cellular differentiation. Because the contribution of this histone variant to carcinogenesis has been reported in several cancer types, but not for PCa, we aimed to investigate the contribution of MacroH2A1 for prostate carcinogenesis.

METHODS

MacroH2A1, MacroH2A1.1 and MacroH2A1.2 isoforms and the corresponding splicing regulators transcript levels were evaluated by RT-qPCR, in a tissue cohort composed by PCa, prostatic intraepithelial neoplasia (PIN) and normal prostate cases. Knockdown for MacroH2A1 and MacroH2A1.1 was performed through lentiviral transduction in DU145 cells, and MacroH2A1.1 overexpression was achieved in LNCaP cells by plasmid transfection, followed by functional assays. Biological and/or experimental replicates were performed when necessary, and specific statistical tests were applied to perform data analysis.

RESULTS

MacroH2A1.1 transcript levels were downregulated in PIN and primary PCa compared to normal prostate tissues. The same was found for QKI, a MacroH2A1.1's splicing regulator. Moreover, lower MacroH2A1.1 and QKI expression levels associated with less differentiated tumors (Gleason score ≥ 7). Interestingly, MacroH2A1.1, but more impressively DDX17 (AUC = 0.93; p < 0.0001) and QKI (AUC = 0.94; p < 0.0001), accurately discriminated cancerous from noncancerous prostate tissues. Furthermore, in PCa cell lines, total MacroH2A1 knockdown augmented malignant features, whereas MacroH2A1.1 overexpression impressively attenuated the malignant phenotype.

CONCLUSIONS

Overall, our data, derived from primary PCa tissues and cell lines, anticipate a tumor suppressive role for MacroH2A1, particularly for the MacroH2A1.1 isoform, in prostate carcinogenesis.

Histone variant macroH2A: from chromatin deposition to molecular function

The eukaryotic genome is regulated in the context of chromatin. Specialized histones, known as histone variants, incorporate into chromatin to replace their canonical counterparts and represent an important layer of regulation to diversify the structural characteristics and functional outputs of chromatin. MacroH2A is an unusual histone variant with a bulky C-terminal non-histone domain that distinguishes it from all other histones. It is a critical player in stabilizing differentiated cell identity by posing as a barrier to somatic cell reprogramming toward pluripotency and acts as a tumor suppressor in a wide range of cancers. MacroH2A histones are generally regarded as repressive variants that are enriched at the inactive X chromosome (Xi) and broad domains across autosomal chromatin. Recent studies have shed light on to how macroH2A influences transcriptional outputs within distinct genomic contexts and revealed new intriguing molecular functions of macroH2A variants beyond transcriptional regulation. Furthermore, the mechanisms of its mysterious chromatin deposition are beginning to be unraveled, facilitating our understanding of its complex regulation of genome function.

MacroH2A1 chromatin specification requires its docking domain and acetylation of H2B lysine 20

The histone variant macroH2A1 localizes to two functionally distinct chromatin subtypes marked by either H3K27me3 or H2B acetylations, where it is thought to directly regulate transcription. The recent finding, that macroH2A1 regulates mitochondrial respiration by globally dampening PARP activity, requires the field to re-evaluate which functions of macroH2A1 are due to global effects on cellular metabolism and which are direct effects determined by macroH2A1 chromatin localization. Here, we demonstrate macroH2A1 incorporation into H2B-acetylated chromatin requires a feature in its histone-fold domain, distinguishing this process from incorporation into H3K27me3-containing chromatin in which multiple features of macroH2A1 are sufficient for targeting. In addition, we identify H2BK20 acetylation as a critical modification required to target macroH2A1 to H2B-acetylated chromatin. Our findings have allowed us to definitively establish that macroH2A1's regulation of an important transcriptional program, the senescence-associated secretory phenotype (SASP), requires its accurate genomic localization.

The MacroH2A1.1 - PARP1 Axis at the Intersection Between Stress Response and Metabolism

The exchange of replication-coupled canonical histones by histone variants endows chromatin with specific features. The replacement of the canonical H2A histone for the histone variant macroH2A is one of the most remarkable epigenetic modifications. The three vertebrate macroH2A proteins have a unique tripartite structure consisting of H2A-like domain, unstructured linker, and macrodomain. Macrodomains are ancient globular folds that are able to bind nicotinamide adenine dinucleotide (NAD⁺) derived metabolites. Here, we will briefly describe the physiological relevance of the metabolite binding in the context of chromatin. In particular, we will focus on the macroH2A1.1 isoform that binds ADP-ribose and poly-ADP-ribose polymerase 1 (PARP1) enzyme, a cellular stress sensor. We will discuss the impact of this interaction in the context of cancer, senescence, cell stress and energy metabolism.

Histone variant macroH2A1 rewires carbohydrate and lipid metabolism of hepatocellular carcinoma cells towards cancer stem cells

Hepatocellular carcinomas (HCCs) contain a sub-population of cancer stem cells (CSCs) that are responsible for tumor relapse, metastasis, and chemoresistance. We recently showed that loss of macroH2A1, a variant of the histone H2A and an epigenetic regulator of stem-cell function, in HCC leads to CSC-like features such as resistance to chemotherapeutic agents and growth of large and relatively undifferentiated tumors in xenograft models. These HCC cells silenced for macroH2A1 also exhibited stem-like metabolic changes consistent with enhanced glycolysis. However, there is no consensus as to the metabolic characteristics of CSCs that render them adaptable to microenvironmental changes by conveniently shifting energy production source or by acquiring intermediate metabolic phenotypes. Here, we assessed long-term proliferation, energy metabolism, and central carbon metabolism in human hepatoma HepG2 cells depleted in macroH2A1. MacroH2A1-depleted HepG2 cells were insensitive to serum exhaustion and showed two distinct, but interdependent changes in glucose and lipid metabolism in CSCs: (1) massive upregulation of acetyl-coA that is transformed into enhanced lipid content and (2) increased activation of the pentose phosphate pathway, diverting glycolytic intermediates to provide precursors for nucleotide synthesis. Integration of metabolomic analyses with RNA-Seq data revealed a critical role for the Liver X Receptor pathway, whose inhibition resulted in attenuated CSCs-like features. These findings shed light on the metabolic phenotype of epigenetically modified CSC-like hepatic cells, and highlight a potential approach for selective therapeutic targeting.

MacroH2A histone variants limit chromatin plasticity through two distinct mechanisms

MacroH2A histone variants suppress tumor progression and act as epigenetic barriers to induced pluripotency. How they impart their influence on chromatin plasticity is not well understood. Here, we analyze how the different domains of macroH2A proteins contribute to chromatin structure and dynamics. By solving the crystal structure of the macrodomain of human macroH2A2 at 1.7 Å, we find that its putative binding pocket exhibits marked structural differences compared with the macroH2A1.1 isoform, rendering macroH2A2 unable to bind ADP-ribose. Quantitative binding assays show that this specificity is conserved among vertebrate macroH2A isoforms. We further find that macroH2A histones reduce the transient, PARP1-dependent chromatin relaxation that occurs in living cells upon DNA damage through two distinct mechanisms. First, macroH2A1.1 mediates an isoform-specific effect through its ability to suppress PARP1 activity. Second, the unstructured linker region exerts an additional repressive effect that is common to all macroH2A proteins. In the absence of DNA damage, the macroH2A linker is also sufficient for rescuing heterochromatin architecture in cells deficient for macroH2A.

ZFR coordinates crosstalk between RNA decay and transcription in innate immunity

Control of type I interferon production is crucial to combat infection while preventing deleterious inflammatory responses, but the extent of the contribution of post-transcriptional mechanisms to innate immune regulation is unclear. Here, we show that human zinc finger RNA-binding protein (ZFR) represses the interferon response by regulating alternative pre-mRNA splicing. ZFR expression is tightly controlled during macrophage development; monocytes express truncated ZFR isoforms, while macrophages induce full-length ZFR to modulate macrophage-specific alternative splicing. Interferon-stimulated genes are constitutively activated by ZFR depletion, and immunostimulation results in hyper-induction of interferon β (IFNβ/IFNB1). Through whole-genome analyses, we show that ZFR controls interferon signaling by preventing aberrant splicing and nonsense-mediated decay of histone variant macroH2A1/H2AFY mRNAs. Together, our data suggest that regulation of ZFR in macrophage differentiation guards against aberrant interferon responses and reveal a network of mRNA processing and decay that shapes the transcriptional response to infection.

Type I interferon signaling is critical for the control of infection. Here the authors show that zinc finger RNA-binding protein (ZFR) can control type I interferon responses, and that this control is itself regulated by distinct ZFR truncation patterns that differ between monocytes and macrophages.

Histone Variant MacroH2A1 Plays an Isoform-Specific Role in Suppressing Epithelial-Mesenchymal Transition

Epithelial-Mesenchymal Transition (EMT) is a biological program that plays key roles in various developmental and pathological processes. Although much work has been done on signaling pathways and transcription factors regulating EMT, the epigenetic regulation of EMT remains not well understood. Histone variants have been recognized as a key group of epigenetic regulators. Among them, macroH2A1 is involved in stem cell reprogramming and cancer progression. We postulated that macroH2A1 may play a role in EMT, a process involving reprogramming of cellular states. In this study, we demonstrate that expression of macroH2A1 is dramatically reduced during EMT induction in immortalized human mammary epithelial cells (HMLE). Moreover, ectopic expression of the macroH2A1.1 isoform, but not macroH2A1.2, can suppress EMT induction and reduce the stem-like cell population in HMLE. Interestingly, macroH2A1.1 overexpression cannot revert stable mesenchymal cells back to the epithelial state, suggesting a stage-specific role of macroH2A1.1 in EMT. We further pinpointed that the function of macroH2A1.1 in EMT suppression is dependent on its ability to bind the NAD+ metabolite PAR, in agreement with the inability to suppress EMT by macroH2A1.2, which lacks the PAR binding domain. Thus, our work discovered a previously unrecognized isoform-specific function of macroH2A1 in regulating EMT induction.

Histone MacroH2A1: A Chromatin Point of Intersection between Fasting, Senescence and Cellular Regeneration

Histone variants confer chromatin unique properties. They have specific genomic distribution, regulated by specific deposition and removal machineries. Histone variants, mostly of canonical histones H2A, H2B and H3, have important roles in early embryonic development, in lineage commitment of stem cells, in the converse process of somatic cell reprogramming to pluripotency and, in some cases, in the modulation of animal aging and life span. MacroH2A1 is a variant of histone H2A, present in two alternatively exon-spliced isoforms macroH2A1.1 and macroH2A1.2, regulating cell plasticity and proliferation, during pluripotency and tumorigenesis. Furthermore, macroH2A1 participates in the formation of senescence-associated heterochromatic foci (SAHF) in senescent cells, and multiple lines of evidence in genetically modified mice suggest that macroH2A1 integrates nutritional cues from the extracellular environment to transcriptional programs. Here, we review current molecular evidence based on next generation sequencing data, cell assays and in vivo models supporting different mechanisms that could mediate the function of macroH2A1 in health span and life span. We will further discuss context-dependent and isoform-specific functions. The aim of this review is to provide guidance to assess histone variant macroH2A1 potential as a therapeutic intervention point.

The histone variant macroH2A confers functional robustness to the intestinal stem cell compartment

A stem cell's epigenome directs cell fate during development, homeostasis, and regeneration. Epigenetic dysregulation can lead to inappropriate cell fate decisions, aberrant cell function, and even cancer. The histone variant macroH2A has been shown to influence gene expression, guide cell fate, and safeguard against genotoxic stress. Interestingly, mice lacking functional macroH2A histones (hereafter referred to as macroH2A DKO) are viable and fertile; yet suffer from increased perinatal death and reduced weight and size compared to wildtype (WT). Here, we ask whether the ostensible reduced vigor of macroH2A DKO mice extends to intestinal stem cell (ISC) function during homeostasis, regeneration, and oncogenesis. Lgr5-eGFP-IRES-CreERT2 or Hopx-CreERT2::Rosa26-LSL-tdTomato ISC reporter mice or the C57BL/6J-Apcmin/J murine intestinal adenoma model were bred into a macroH2A DKO or strain-matched WT background and assessed for ISC functionality, regeneration and tumorigenesis. High-dose (12Gy) whole-body γ-irradiation was used as an injury model. We show that macroH2A is dispensable for intestinal homeostasis and macroH2A DKO mice have similar numbers of active crypt-base columnar ISCs (CBCs). MacroH2A DKO intestine exhibits impaired regeneration following injury, despite having significantly more putative reserve ISCs. DKO reserve ISCs disproportionately undergo apoptosis compared to WT after DNA damage infliction. Interestingly, a macroH2A DKO background does not significantly increase tumorigenesis in the Apcmin model of intestinal adenoma. We conclude that macroH2A influences reserve ISC number and function during homeostasis and regeneration. These data suggest macroH2A enhances reserve ISC survival after DNA damage and thus confers functional robustness to the intestinal epithelium.

Gene Resistance to Transcriptional Reprogramming following Nuclear Transfer Is Directly Mediated by Multiple Chromatin-Repressive Pathways

Understanding the mechanism of resistance of genes to reactivation will help improve the success of nuclear reprogramming. Using mouse embryonic fibroblast nuclei with normal or reduced DNA methylation in combination with chromatin modifiers able to erase H3K9me3, H3K27me3, and H2AK119ub1 from transplanted nuclei, we reveal the basis for resistance of genes to transcriptional reprogramming by oocyte factors. A majority of genes is affected by more than one type of treatment, suggesting that resistance can require repression through multiple epigenetic mechanisms. We classify resistant genes according to their sensitivity to 11 chromatin modifier combinations, revealing the existence of synergistic as well as adverse effects of chromatin modifiers on removal of resistance. We further demonstrate that the chromatin modifier USP21 reduces resistance through its H2AK119 deubiquitylation activity. Finally, we provide evidence that H2A ubiquitylation also contributes to resistance to transcriptional reprogramming in mouse nuclear transfer embryos.

•

Identification of genes resistant to direct transcriptional reprogramming

•

Determination of resistant gene sensitivity to 11 chromatin modifier combinations

•

USP21 removes resistance through its H2AK119 deubiquitylation activity

•

USP21 improves the reprogramming of gene expression in two-cell-stage mouse embryos

Epigenetic regulation of somatic cell reprogramming

Pluripotent stem cells, having self-renewal capacities and multi-lineage differentiation abilities, offer great potential in disease modeling and therapeutic applications. The successful generation of induced pluripotent stem cells (iPSCs) by the Yamanaka group in 2006 is a milestone event in both reprogramming and stem cell research fields, which makes in vitro somatic cell reprogramming and personalized stem cell therapy feasible. During the past 10 years, several important progresses have been made in uncovering the molecular mechanisms involved in the reprogramming process, which shed light on improving the reprogramming efficiency and iPSC quality. Here, we briefly review the important progresses in the epigenetic regulation including histone and DNA modifications during somatic cell reprogramming.

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.

Multiple Roles of MYC in Integrating Regulatory Networks of Pluripotent Stem Cells

Pluripotent stem cells (PSCs) are defined by their self-renewal potential, which permits their unlimited propagation, and their pluripotency, being able to generate cell of the three embryonic lineages. These properties render PSCs a valuable tool for both basic and medical research. To induce and stabilize the pluripotent state, complex circuitries involving signaling pathways, transcription regulators and epigenetic mechanisms converge on a core transcriptional regulatory network of PSCs, thus determining their cell identity. Among the transcription factors, MYC represents a central hub, which modulates and integrates multiple mechanisms involved both in the maintenance of pluripotency and in cell reprogramming. Indeed, it instructs the PSC-specific cell cycle, metabolism and epigenetic landscape, contributes to limit exit from pluripotency and modulates signaling cascades affecting the PSC identity. Moreover, MYC extends its regulation on pluripotency by controlling PSC-specific non-coding RNAs. In this report, we review the MYC-controlled networks, which support the pluripotent state and discuss how their perturbation could affect cell identity. We further discuss recent finding demonstrating a central role of MYC in triggering epigenetic memory in PSCs, which depends on the establishment of a WNT-centered self-reinforcing circuit. Finally, we comment on the therapeutic implications of the role of MYC in affecting PSCs. Indeed, PSCs are used for both disease and cancer modeling and to derive cells for regenerative medicine. For these reasons, unraveling the MYC-mediated mechanism in those cells is fundamental to exploit their full potential and to identify therapeutic targets.

The characterization of macroH2A beyond vertebrates supports an ancestral origin and conserved role for histone variants in chromatin

Histone variants play a critical role in chromatin structure and epigenetic regulation. These "deviant" proteins have been historically considered as the evolutionary descendants of ancestral canonical histones, helping specialize the nucleosome structure during eukaryotic evolution. Such view is now challenged by 2 major observations: first, canonical histones present extremely unique features not shared with any other genes; second, histone variants are widespread across many eukaryotic groups. The present work further supports the ancestral nature of histone variants by providing the first in vivo characterization of a functional macroH2A histone (a variant long defined as a specific refinement of vertebrate chromatin) in a non-vertebrate organism (the mussel Mytilus) revealing its recruitment into heterochromatic fractions of actively proliferating tissues. Combined with in silico analyses of genomic data, these results provide evidence for the widespread presence of macroH2A in metazoan animals, as well as in the holozoan Capsaspora, supporting an evolutionary origin for this histone variant lineage before the radiation of Filozoans (including Filasterea, Choanoflagellata and Metazoa). Overall, the results presented in this work help configure a new evolutionary scenario in which histone variants, rather than modern "deviants" of canonical histones, would constitute ancient components of eukaryotic chromatin.

A Tale of Two Cities: How Xist and its partners localize to and silence the bicompartmental X

Sex chromosomal dosage compensation in mammals takes the form of X chromosome inactivation (XCI), driven by the non-coding RNA Xist. In contrast to dosage compensation systems of flies and worms, mammalian XCI has to restrict its function to the Xist-producing X chromosome, while leaving autosomes and active X untouched. The mechanisms behind the long-range yet cis-specific localization and silencing activities of Xist have long been enigmatic, but genomics, proteomics, super-resolution microscopy, and innovative genetic approaches have produced significant new insights in recent years. In this review, I summarize and integrate these findings with a particular focus on the redundant yet mutually reinforcing pathways that enable long-term transcriptional repression throughout the soma. This includes an exploration of concurrent epigenetic changes acting in parallel within two distinct compartments of the inactive X. I also examine how Polycomb repressive complexes 1 and 2 and macroH2A may bridge XCI establishment and maintenance. XCI is a remarkable phenomenon that operates across multiple scales, combining changes in nuclear architecture, chromosome topology, chromatin compaction, and nucleosome/nucleotide-level epigenetic cues. Learning how these pathways act in concert likely holds the answer to the riddle posed by Cattanach's and other autosomal translocations: What makes the X especially receptive to XCI?

Induced Pluripotent Stem Cells: Generation Strategy and Epigenetic Mystery behind Reprogramming

Possessing the ability of self-renewal with immortalization and potential for differentiation into different cell types, stem cells, particularly embryonic stem cells (ESC), have attracted significant attention since their discovery. As ESC research has played an essential role in developing our understanding of the mechanisms underlying reproduction, development, and cell (de)differentiation, significant efforts have been made in the biomedical study of ESC in recent decades. However, such studies of ESC have been hampered by the ethical issues and technological challenges surrounding them, therefore dramatically inhibiting the potential applications of ESC in basic biomedical studies and clinical medicine. Induced pluripotent stem cells (iPSCs), generated from the reprogrammed somatic cells, share similar characteristics including but not limited to the morphology and growth of ESC, self-renewal, and potential differentiation into various cell types. The discovery of the iPSC, unhindered by the aforementioned limitations of ESC, introduces a viable alternative to ESC. More importantly, the applications of iPSC in the development of disease models such as neurodegenerative disorders greatly enhance our understanding of the pathogenesis of such diseases and also facilitate the development of clinical therapeutic strategies using iPSC generated from patient somatic cells to avoid an immune rejection. In this review, we highlight the advances in iPSCs generation methods as well as the mechanisms behind their reprogramming. We also discuss future perspectives for the development of iPSC generation methods with higher efficiency and safety.

H3K9me3-Dependent Heterochromatin: Barrier to Cell Fate Changes

Establishing and maintaining cell identity depends on the proper regulation of gene expression, as specified by transcription factors and reinforced by epigenetic mechanisms. Among the epigenetic mechanisms, heterochromatin formation is crucial for the preservation of genome stability and the cell type-specific silencing of genes. The heterochromatin-associated histone mark H3K9me3, although traditionally associated with the noncoding portions of the genome, has emerged as a key player in repressing lineage-inappropriate genes and shielding them from activation by transcription factors. Here we describe the role of H3K9me3 heterochromatin in impeding the reprogramming of cell identity and the mechanisms by which H3K9me3 is reorganized during development and cell fate determination.

MacroH2A1 associates with nuclear lamina and maintains chromatin architecture in mouse liver cells

In the interphase nucleus, chromatin is organized into three-dimensional conformation to coordinate genome functions. The lamina-chromatin association is important to facilitate higher-order chromatin in mammalian cells, but its biological significances and molecular mechanisms remain poorly understood. One obstacle is that the list of lamina-associated proteins remains limited, presumably due to the inherent insolubility of lamina proteins. In this report, we identified 182 proteins associated with lamin B1 (a constitutive component of lamina) in mouse hepatocytes, by adopting virus-based proximity-dependent biotin identification. These proteins are functionally related to biological processes such as chromatin organization. As an example, we validated the association between lamin B1 and core histone macroH2A1, a histone associated with repressive chromatin. Furthermore, we mapped Lamina-associated domains (LADs) in mouse liver cells and found that boundaries of LADs are enriched for macroH2A. More interestingly, knocking-down of macroH2A1 resulted in the release of heterochromatin foci marked by histone lysine 9 trimethylation (H3K9me3) and the decondensation of global chromatin structure. However, down-regulation of lamin B1 led to redistribution of macroH2A1. Taken together, our data indicated that macroH2A1 is associated with lamina and is required to maintain chromatin architecture in mouse liver cells.

Elevated levels of ZAC1 disrupt neurogenesis and promote rapid in vivo reprogramming

The zinc finger transcription factor Zac1 is expressed in dividing progenitors of the nervous system with expression levels negatively controlled by genomic imprinting. To explore the consequences of elevated ZAC1 levels during neurogenesis we overexpressed it in the developing CNS. Increased levels of ZAC1 rapidly promoted upregulation of CDK inhibitors P57 and P27 followed by cell cycle exit. Surprisingly this was accompanied by stalled neuronal differentiation. Genome wide expression analysis of cortical cells overexpressing Zac1 revealed a decrease in neuronal gene expression and an increased expression of imprinted genes, factors regulating mesoderm formation as well as features of differentiated muscle. In addition, we observed a rapid induction of several genes regulating pluripotency. Taken together, our data suggests that expression levels of Zac1 need to be kept under strict control to avoid premature cell cycle exit, disrupted neurogenesis and aberrant expression of non-neuronal genes including pluripotency associated factors.

X chromosome reactivation in reprogramming and in development

Dramatic epigenetic changes take place during mammalian differentiation from the naïve pluripotent state including the silencing of one of the two X chromosomes in female cells through X chromosome inactivation. Conversely, reprogramming of somatic cells to naive pluripotency is coupled to X chromosome reactivation (XCR). Recent studies in the mouse system have shed light on the mechanisms of XCR by uncovering the timing and steps of XCR during reprogramming to induced pluripotent stem cells (iPSCs), allowing the generation of testable hypotheses during embryogenesis. In contrast, analyses of the X chromosome in human iPSCs have revealed important differences between mouse and human reprogramming processes that can partially be explained by the establishment of distinct pluripotent states and impact disease modeling and the application of human pluripotent stem cells. Here, we review recent literature on XCR as a readout and determinant of reprogramming to pluripotency.

Reprogramming barriers and enhancers: strategies to enhance the efficiency and kinetics of induced pluripotency

Induced pluripotent stem cells are powerful tools for disease modeling, drug screening, and cell transplantation therapies. These cells can be generated directly from somatic cells by ectopic expression of defined factors through a reprogramming process. However, pluripotent reprogramming is an inefficient process because of various defined and unidentified barriers. Recent studies dissecting the molecular mechanisms of reprogramming have methodically improved the quality, ease, and efficiency of reprogramming. Different strategies have been applied for enhancing reprogramming efficiency, including depletion/inhibition of barriers (p53, p21, p57, p16(Ink4a)/p19(Arf), Mbd3, etc.), overexpression of enhancing genes (e.g., FOXH1, C/EBP alpha, UTF1, and GLIS1), and administration of certain cytokines and small molecules. The current review provides an in-depth overview of the cutting-edge findings regarding distinct barriers of reprogramming to pluripotency and strategies to enhance reprogramming efficiency. By incorporating the mechanistic insights from these recent findings, a combined method of inhibition of roadblocks and application of enhancing factors may yield the most reliable and effective approach in pluripotent reprogramming.

Heterochromatin Protein 1β (HP1β) has distinct functions and distinct nuclear distribution in pluripotent versus differentiated cells

Background

Pluripotent embryonic stem cells (ESCs) have the unique ability to differentiate into every cell type and to self-renew. These characteristics correlate with a distinct nuclear architecture, epigenetic signatures enriched for active chromatin marks and hyperdynamic binding of structural chromatin proteins. Recently, several chromatin-related proteins have been shown to regulate ESC pluripotency and/or differentiation, yet the role of the major heterochromatin proteins in pluripotency is unknown.

Results

Here we identify Heterochromatin Protein 1β (HP1β) as an essential protein for proper differentiation, and, unexpectedly, for the maintenance of pluripotency in ESCs. In pluripotent and differentiated cells HP1β is differentially localized and differentially associated with chromatin. Deletion of HP1β, but not HP1α, in ESCs provokes a loss of the morphological and proliferative characteristics of embryonic pluripotent cells, reduces expression of pluripotency factors and causes aberrant differentiation. However, in differentiated cells, loss of HP1β has the opposite effect, perturbing maintenance of the differentiation state and facilitating reprogramming to an induced pluripotent state. Microscopy, biochemical fractionation and chromatin immunoprecipitation reveal a diffuse nucleoplasmic distribution, weak association with chromatin and high expression levels for HP1β in ESCs. The minor fraction of HP1β that is chromatin-bound in ESCs is enriched within exons, unlike the situation in differentiated cells, where it binds heterochromatic satellite repeats and chromocenters.

Conclusions

We demonstrate an unexpected duality in the role of HP1β: it is essential in ESCs for maintaining pluripotency, while it is required for proper differentiation in differentiated cells. Thus, HP1β function both depends on, and regulates, the pluripotent state.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0760-8) contains supplementary material, which is available to authorized users.