Mohammad M. Karimi

Vitamin C induces Tet-dependent DNA demethylation and a blastocyst-like state in ES cells

DNA methylation is a heritable epigenetic modification involved in gene silencing, imprinting, and the suppression of retrotransposons. Global DNA demethylation occurs in the early embryo and the germ line, and may be mediated by Tet (ten eleven translocation) enzymes, which convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). Tet enzymes have been studied extensively in mouse embryonic stem (ES) cells, which are generally cultured in the absence of vitamin C, a potential cofactor for Fe(ii) 2-oxoglutarate dioxygenase enzymes such as Tet enzymes. Here we report that addition of vitamin C to mouse ES cells promotes Tet activity, leading to a rapid and global increase in 5hmC. This is followed by DNA demethylation of many gene promoters and upregulation of demethylated germline genes. Tet1 binding is enriched near the transcription start site of genes affected by vitamin C treatment. Importantly, vitamin C, but not other antioxidants, enhances the activity of recombinant Tet1 in a biochemical assay, and the vitamin-C-induced changes in 5hmC and 5mC are entirely suppressed in Tet1 and Tet2 double knockout ES cells. Vitamin C has a stronger effect on regions that gain methylation in cultured ES cells compared to blastocysts, and in vivo are methylated only after implantation. In contrast, imprinted regions and intracisternal A particle retroelements, which are resistant to demethylation in the early embryo, are resistant to vitamin-C-induced DNA demethylation. Collectively, the results of this study establish vitamin C as a direct regulator of Tet activity and DNA methylation fidelity in ES cells.DNA methylation is a heritable epigenetic modification involved in gene silencing, imprinting, and the suppression of retrotransposons1. Global DNA demethylation occurs in the early embryo and the germline2,3 and may be mediated by Tet (ten-eleven-translocation) enzymes4–6, which convert 5-methylcytosine (mC) to 5-hydroxymethylcytosine (hmC)7. Tet enzymes have been extensively studied in mouse embryonic stem cells (ESCs)8–12, which are generally cultured in the absence of Vitamin C (VitC), a potential co-factor for Fe(II) 2-oxoglutarate dioxygenase enzymes like Tets. Here we report that addition of VitC to ESCs promotes Tet activity leading to a rapid and global increase in hmC. This is followed by DNA demethylation of numerous gene promoters and up-regulation of demethylated germline genes. Tet1 binding is enriched near the transcription start site (TSS) of genes affected by VitC treatment. Importantly, VitC, but not other antioxidants, enhances the activity of recombinant Tet1 in a biochemical assay and the VitC-induced changes in hmC and mC are entirely suppressed in Tet1/2 double knockout (Tet DKO) ESCs. VitC has the strongest effects on regions that gain methylation in cultured ESCs compared to blastocysts and in vivo are methylated only after implantation. In contrast, imprinted regions and intracisternal A-particle (IAP) retroelements, which are resistant to demethylation in the early embryo2,13, are resistant to VitC-induced DNA demethylation. Collectively, this study establishes VitC as a direct regulator of Tet activity and DNA methylation fidelity in ESCs.

Setdb1 is required for germline development and silencing of H3K9me3-marked endogenous retroviruses in primordial germ cells

Transcription of endogenous retroviruses (ERVs) is inhibited by de novo DNA methylation during gametogenesis, a process initiated after birth in oocytes and at approximately embryonic day 15.5 (E15.5) in prospermatogonia. Earlier in germline development, the genome, including most retrotransposons, is progressively demethylated. Young ERVK and ERV1 elements, however, retain intermediate methylation levels. As DNA methylation reaches a low point in E13.5 primordial germ cells (PGCs) of both sexes, we determined whether retrotransposons are marked by H3K9me3 and H3K27me3 using a recently developed low-input ChIP-seq (chromatin immunoprecipitation [ChIP] combined with deep sequencing) method. Although these repressive histone modifications are found predominantly on distinct genomic regions in E13.5 PGCs, they concurrently mark partially methylated long terminal repeats (LTRs) and LINE1 elements. Germline-specific conditional knockout of the H3K9 methyltransferase SETDB1 yields a decrease of both marks and DNA methylation at H3K9me3-enriched retrotransposon families. Strikingly, Setdb1 knockout E13.5 PGCs show concomitant derepression of many marked ERVs, including intracisternal A particle (IAP), ETn, and ERVK10C elements, and ERV-proximal genes, a subset in a sex-dependent manner. Furthermore, Setdb1 deficiency is associated with a reduced number of male E13.5 PGCs and postnatal hypogonadism in both sexes. Taken together, these observations reveal that SETDB1 is an essential guardian against proviral expression prior to the onset of de novo DNA methylation in the germline.

An ultra-low-input native ChIP-seq protocol for genome-wide profiling of rare cell populations

Combined chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) has enabled genome-wide epigenetic profiling of numerous cell lines and tissue types. A major limitation of ChIP-seq, however, is the large number of cells required to generate high-quality data sets, precluding the study of rare cell populations. Here, we present an ultra-low-input micrococcal nuclease-based native ChIP (ULI-NChIP) and sequencing method to generate genome-wide histone mark profiles with high resolution from as few as 10(3) cells. We demonstrate that ULI-NChIP-seq generates high-quality maps of covalent histone marks from 10(3) to 10(6) embryonic stem cells. Subsequently, we show that ULI-NChIP-seq H3K27me3 profiles generated from E13.5 primordial germ cells isolated from single male and female embryos show high similarity to recent data sets generated using 50-180 × more material. Finally, we identify sexually dimorphic H3K27me3 enrichment at specific genic promoters, thereby illustrating the utility of this method for generating high-quality and -complexity libraries from rare cell populations.

Retrotransposon-Induced Heterochromatin Spreading in the Mouse Revealed by Insertional Polymorphisms

The “arms race” relationship between transposable elements (TEs) and their host has promoted a series of epigenetic silencing mechanisms directed against TEs. Retrotransposons, a class of TEs, are often located in repressed regions and are thought to induce heterochromatin formation and spreading. However, direct evidence for TE–induced local heterochromatin in mammals is surprisingly scarce. To examine this phenomenon, we chose two mouse embryonic stem (ES) cell lines that possess insertionally polymorphic retrotransposons (IAP, ETn/MusD, and LINE elements) at specific loci in one cell line but not the other. Employing ChIP-seq data for these cell lines, we show that IAP elements robustly induce H3K9me3 and H4K20me3 marks in flanking genomic DNA. In contrast, such heterochromatin is not induced by LINE copies and only by a minority of polymorphic ETn/MusD copies. DNA methylation is independent of the presence of IAP copies, since it is present in flanking regions of both full and empty sites. Finally, such spreading into genes appears to be rare, since the transcriptional start sites of very few genes are less than one Kb from an IAP. However, the B3galtl gene is subject to transcriptional silencing via IAP-induced heterochromatin. Hence, although rare, IAP-induced local heterochromatin spreading into nearby genes may influence expression and, in turn, host fitness.

Distinct roles of KAP1, HP1 and G9a/GLP in silencing of the two-cell-specific retrotransposon MERVL in mouse ES cells

BackgroundIn mouse embryonic stem cells (mESCs), transcriptional silencing of numerous class I and II endogenous retroviruses (ERVs), including IAP, ETn and MMERVK10C, is dependent upon the H3K9 methyltransferase (KMTase) SETDB1/ESET and its binding partner KAP1/TRIM28. In contrast, the H3K9 KMTases G9a and GLP and HP1 proteins are dispensable for this process. Intriguingly, MERVL retroelements are actively transcribed exclusively in the two-cell (2C) embryo, but the molecular basis of silencing of these class III ERVs at later developmental stages has not been systematically addressed.ResultsHere, we characterized the roles of these chromatin factors in MERVL silencing in mESCs. While MMERVK10C and IAP ERVs are bound by SETDB1 and KAP1 and are induced following their deletion, MERVL ERVs show relatively low levels of SETDB1 and KAP1 binding and are upregulated exclusively following KAP1 depletion, indicating that KAP1 influences MERVL expression independent of SETDB1. In contrast to class I and class II ERVs, MERVL and MERVL LTR-driven genic transcripts are also upregulated following depletion of G9a or GLP, and G9a binds directly to these ERVs. Consistent with a direct role for H3K9me2 in MERVL repression, these elements are highly enriched for G9a-dependent H3K9me2, and catalytically active G9a is required for silencing of MERVL LTR-driven transcripts. MERVL is also derepressed in HP1α and HP1β KO ESCs. However, like KAP1, HP1α and HP1β are only modestly enriched at MERVL relative to IAP LTRs. Intriguingly, as recently shown for KAP1, RYBP, LSD1 and G9a-deficient mESCs, many genes normally expressed in the 2C embryo are also induced in HP1 KO mESCs, revealing that aberrant expression of a subset of 2C-specific genes is a common feature in each of these KO lines.ConclusionsOur results indicate that G9a and GLP, which are not required for silencing of class I and II ERVs, are recruited to MERVL elements and play a direct role in silencing of these class III ERVs, dependent upon G9a catalytic activity. In contrast, induction of MERVL expression in KAP1, HP1α and HP1β KO ESCs may occur predominantly as a consequence of indirect effects, in association with activation of a subset of 2C-specific genes.

Distinct isoform of FABP7 revealed by screening for retroelement-activated genes in diffuse large B-cell lymphoma

Significance Sequences derived from transposable elements (TEs) are abundant in the human genome and can influence gene expression. In normal cells, most TEs are silenced by epigenetic mechanisms such as DNA methylation but, in cancer, normally dormant TEs can become active. We hypothesized that cancer-specific release of epigenetic suppression of TEs could result in gene expression perturbations, which could promote oncogenesis. Using a bioinformatics method, we identified many genes expressed in diffuse large B-cell lymphoma (DLBCL) via activation of TE promoters. Further analysis of one, FABP7, showed it was expressed in some DLBCL samples through use of a TE promoter. The TE-driven FABP7 transcript encodes a novel isoform of the protein, which is required for optimal DLBCL cell line proliferation. Remnants of ancient transposable elements (TEs) are abundant in mammalian genomes. These sequences harbor multiple regulatory motifs and hence are capable of influencing expression of host genes. In response to environmental changes, TEs are known to be released from epigenetic repression and to become transcriptionally active. Such activation could also lead to lineage-inappropriate activation of oncogenes, as one study described in Hodgkin lymphoma. However, little further evidence for this mechanism in other cancers has been reported. Here, we reanalyzed whole transcriptome data from a large cohort of patients with diffuse large B-cell lymphoma (DLBCL) compared with normal B-cell centroblasts to detect genes ectopically expressed through activation of TE promoters. We have identified 98 such TE-gene chimeric transcripts that were exclusively expressed in primary DLBCL cases and confirmed several in DLBCL-derived cell lines. We further characterized a TE-gene chimeric transcript involving a fatty acid-binding protein gene (LTR2-FABP7), normally expressed in brain, that was ectopically expressed in a subset of DLBCL patients through the use of an endogenous retroviral LTR promoter of the LTR2 family. The LTR2-FABP7 chimeric transcript encodes a novel chimeric isoform of the protein with characteristics distinct from native FABP7. In vitro studies reveal a dependency for DLBCL cell line proliferation and growth on LTR2-FABP7 chimeric protein expression. Taken together, these data demonstrate the significance of TEs as regulators of aberrant gene expression in cancer and suggest that LTR2-FABP7 may contribute to the pathogenesis of DLBCL in a subgroup of patients.

hnRNP K coordinates transcriptional silencing by SETDB1 in embryonic stem cells.

Retrotransposition of endogenous retroviruses (ERVs) poses a substantial threat to genome stability. Transcriptional silencing of a subset of these parasitic elements in early mouse embryonic and germ cell development is dependent upon the lysine methyltransferase SETDB1, which deposits H3K9 trimethylation (H3K9me3) and the co-repressor KAP1, which binds SETDB1 when SUMOylated. Here we identified the transcription co-factor hnRNP K as a novel binding partner of the SETDB1/KAP1 complex in mouse embryonic stem cells (mESCs) and show that hnRNP K is required for ERV silencing. RNAi-mediated knockdown of hnRNP K led to depletion of H3K9me3 at ERVs, concomitant with de-repression of proviral reporter constructs and specific ERV subfamilies, as well as a cohort of germline-specific genes directly targeted by SETDB1. While hnRNP K recruitment to ERVs is dependent upon KAP1, SETDB1 binding at these elements requires hnRNP K. Furthermore, an intact SUMO conjugation pathway is necessary for SETDB1 recruitment to proviral chromatin and depletion of hnRNP K resulted in reduced SUMOylation at ERVs. Taken together, these findings reveal a novel regulatory hierarchy governing SETDB1 recruitment and in turn, transcriptional silencing in mESCs.

An interactive analysis and exploration tool for epigenomic data

In this design study, we present an analysis and abstraction of the data and tasks related to the domain of epigenomics, and the design and implementation of an interactive tool to facilitate data analysis and visualization in this domain. Epigenomic data can be grouped into subsets either by k‐means clustering or by querying for combinations of presence or absence of signal (on/off) in different epigenomic experiments. These steps can easily be interleaved and the comparison of different workflows is explicitly supported. We took special care to contain the exponential expansion of possible on/off combinations by creating a novel querying interface. An interactive heat map facilitates the exploration and comparison of different clusters. We validated our iterative design by working closely with two groups of biologists on different biological problems. Both groups quickly found new insight into their data as well as claimed that our tool would save them several hours or days of work over using existing tools.

VisRseq: R-based visual framework for analysis of sequencing data

BackgroundSeveral tools have been developed to enable biologists to perform initial browsing and exploration of sequencing data. However the computational tool set for further analyses often requires significant computational expertise to use and many of the biologists with the knowledge needed to interpret these data must rely on programming experts.ResultsWe present VisRseq, a framework for analysis of sequencing datasets that provides a computationally rich and accessible framework for integrative and interactive analyses without requiring programming expertise. We achieve this aim by providing R apps, which offer a semi-auto generated and unified graphical user interface for computational packages in R and repositories such as Bioconductor. To address the interactivity limitation inherent in R libraries, our framework includes several native apps that provide exploration and brushing operations as well as an integrated genome browser. The apps can be chained together to create more powerful analysis workflows.ConclusionsTo validate the usability of VisRseq for analysis of sequencing data, we present two case studies performed by our collaborators and report their workflow and insights.

ALEA: a toolbox for allele-specific epigenomics analysis

The assessment of expression and epigenomic status using sequencing based methods provides an unprecedented opportunity to identify and correlate allelic differences with epigenomic status. We present ALEA, a computational toolbox for allele-specific epigenomics analysis, which incorporates allelic variation data within existing resources, allowing for the identification of significant associations between epigenetic modifications and specific allelic variants in human and mouse cells. ALEA provides a customizable pipeline of command line tools for allele-specific analysis of next-generation sequencing data (ChIP-seq, RNA-seq, etc.) that takes the raw sequencing data and produces separate allelic tracks ready to be viewed on genome browsers. The pipeline has been validated using human and hybrid mouse ChIP-seq and RNA-seq data. AVAILABILITY The package, test data and usage instructions are available online at http://www.bcgsc.ca/platform/bioinfo/software/alea CONTACT: : mkarimi1@interchange.ubc.ca or sjones@bcgsc.ca Supplementary information: Supplementary data are available at Bioinformatics online.