Epigenetics of Genes Preferentially Expressed in Dissimilar Cell Populations: Myoblasts and Cerebellum.

While studying myoblast methylomes and transcriptomes, we found that CDH15 had a remarkable preference for expression in both myoblasts and cerebellum. To understand how widespread such a relationship was and its epigenetic and biological correlates, we systematically looked for genes with similar transcription profiles and analyzed their DNA methylation and chromatin state and accessibility profiles in many different cell populations. Twenty genes were expressed preferentially in myoblasts and cerebellum (Myob/Cbl genes). Some shared DNA hypo- or hypermethylated regions in myoblasts and cerebellum. Particularly striking was ZNF556, whose promoter is hypomethylated in expressing cells but highly methylated in the many cell populations that do not express the gene. In reporter gene assays, we demonstrated that its promoter's activity is methylation sensitive. The atypical epigenetics of ZNF556 may have originated from its promoter's hypomethylation and selective activation in sperm progenitors and oocytes. Five of the Myob/Cbl genes (KCNJ12, ST8SIA5, ZIC1, VAX2, and EN2) have much higher RNA levels in cerebellum than in myoblasts and displayed myoblast-specific hypermethylation upstream and/or downstream of their promoters that may downmodulate expression. Differential DNA methylation was associated with alternative promoter usage for Myob/Cbl genes MCF2L, DOK7, CNPY1, and ANK1. Myob/Cbl genes PAX3, LBX1, ZNF556, ZIC1, EN2, and VAX2 encode sequence-specific transcription factors, which likely help drive the myoblast and cerebellum specificity of other Myob/Cbl genes. This study extends our understanding of epigenetic/transcription associations related to differentiation and may help elucidate relationships between epigenetic signatures and muscular dystrophies or cerebellar-linked neuropathologies.

SRCAP mutations drive clonal hematopoiesis through epigenetic and DNA repair dysregulation.

Somatic mutations accumulate in all cells with age and can confer a selective advantage, leading to clonal expansion over time. In hematopoietic cells, mutations in a subset of genes regulating DNA repair or epigenetics frequently lead to clonal hematopoiesis (CH). Here, we describe the context and mechanisms that lead to enrichment of hematopoietic stem cells (HSCs) with mutations in SRCAP, which encodes a chromatin remodeler that also influences DNA repair. We show that SRCAP mutations confer a selective advantage in human cells and in mice upon treatment with the anthracycline-class chemotherapeutic doxorubicin and bone marrow transplantation. Furthermore, Srcap mutations lead to a lymphoid-biased expansion, driven by loss of SRCAP-regulated H2A.Z deposition and increased DNA repair. Altogether, we demonstrate that SRCAP operates at the intersection of multiple pathways in stem and progenitor cells, offering a new perspective on the functional impact of genetic variants that promote stem cell competition in the hematopoietic system.

Universal NicE-Seq: A Simple and Quick Method for Accessible Chromatin Detection in Fixed Cells.

Genome-wide accessible chromatin sequencing and identification has enabled deciphering the epigenetic information encoded in chromatin, revealing accessible promoters, enhancers, nucleosome positioning, transcription factor occupancy, and other chromosomal protein binding. The starting biological materials are often fixed using formaldehyde crosslinking. Here, we describe accessible chromatin library preparation from low numbers of formaldehyde-crosslinked cells using a modified nick translation method, where a nicking enzyme nicks one strand of DNA and DNA polymerase incorporates biotin-conjugated dATP, dCTP, and methyl-dCTP. Once the DNA is labeled, it can be isolated for NGS library preparation. We termed this method as universal NicE-seq (nicking enzyme-assisted sequencing). We also demonstrate a single tube method that enables direct NGS library preparation from low cell numbers without DNA purification. Furthermore, we demonstrated universal NicE-seq on FFPE tissue section sample.

Measuring Inaccessible Chromatin Genome-Wide Using Protect-seq.

Chromatin accessibility has been an immensely powerful metric for identifying and understanding regulatory elements in the genome. Many important regulatory elements, such as enhancers and transcriptional start sites, are characterized by "open" or nucleosome-free regions. Understanding the areas of the genome that are not considered open chromatin has been more difficult. Protect-seq is a genomics technique that aims to identify inaccessible chromatin associated with the nuclear periphery. These regions are enriched for histone modifications associated with transcriptional repression and correlate with loci identified by other techniques measuring heterochromatin and peripheral localization. Here, we discuss the protocol and best practices to perform Protect-seq.

NicE-viewSeq: An Integrative Visualization and Genomics Method to Detect Accessible Chromatin in Fixed Cells.

A novel genome-wide accessible chromatin visualization, quantitation, and sequencing method is described, which allows in situ fluorescence visualization and sequencing of the accessible chromatin in the mammalian cell. The cells are fixed by formaldehyde crosslinking, and processed using a modified nick translation method, where a nicking enzyme nicks one strand of DNA, and DNA polymerase incorporates biotin-conjugated dCTP, 5-methyl-dCTP, Fluorescein-12-dATP or Texas Red-5-dATP, dGTP, and dTTP. This allows accessible chromatin DNA to be labeled for visualization and on bead NGS library preparation. This technology allows cellular level chromatin accessibility quantification and genomic analysis of the epigenetic information in the chromatin, particularly accessible promoter, enhancers, nucleosome positioning, transcription factor occupancy, and other chromosomal protein binding.

Diverse silent chromatin states modulate genome compartmentalization and loop extrusion barriers.

The relationships between chromosomal compartmentalization, chromatin state and function are poorly understood. Here by profiling long-range contact frequencies in HCT116 colon cancer cells, we distinguish three silent chromatin states, comprising two types of heterochromatin and a state enriched for H3K9me2 and H2A.Z that exhibits neutral three-dimensional interaction preferences and which, to our knowledge, has not previously been characterized. We find that heterochromatin marked by H3K9me3, HP1α and HP1ß correlates with strong compartmentalization. We demonstrate that disruption of DNA methyltransferase activity greatly remodels genome compartmentalization whereby domains lose H3K9me3-HP1α/ß binding and acquire the neutrally interacting state while retaining late replication timing. Furthermore, we show that H3K9me3-HP1α/ß heterochromatin is permissive to loop extrusion by cohesin but refractory to CTCF binding. Together, our work reveals a dynamic structural and organizational diversity of the silent portion of the genome and establishes connections between the regulation of chromatin state and chromosome organization, including an interplay between DNA methylation, compartmentalization and loop extrusion.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage.

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Poly ADP-ribosylation of SET8 leads to aberrant H4K20 methylation in mammalian nuclear genome.

In mammalian cells, SET8 mediated Histone H4 Lys 20 monomethylation (H4K20me1) has been implicated in regulating mitotic condensation, DNA replication, DNA damage response, and gene expression. Here we show SET8, the only known enzyme for H4K20me1 is post-translationally poly ADP-ribosylated by PARP1 on lysine residues. PARP1 interacts with SET8 in a cell cycle-dependent manner. Poly ADP-ribosylation on SET8 renders it catalytically compromised, and degradation via ubiquitylation pathway. Knockdown of PARP1 led to an increase of SET8 protein levels, leading to aberrant H4K20me1 and H4K20me3 domains in the genome. H4K20me1 is associated with higher gene transcription levels while the increase of H4K20me3 levels was predominant in DNA repeat elements. Hence, SET8 mediated chromatin remodeling in mammalian cells are modulated by poly ADP-ribosylation by PARP1.

Polyglutamylation: biology and analysis.

Polyglutamylation is a posttranslational modification (PTM) that adds several glutamates on glutamate residues in the form of conjugated peptide chains by a family of enzymes known as polyglutamylases. Polyglutamylation is well documented in microtubules. Polyglutamylated microtubules consist of different α- and ß-tubulin subunits with varied number of added glutamate residues. Kinetic control and catalytic rates of tubulin modification by polyglutamylases influence the polyglutamylation pattern of functional microtubules. The recent studies uncovered catalytic mechanisms of the glutamylation enzymes family, particularly tubulin tyrosine ligase-like (TTLL). Variable length polyglutamylation of primary sequence glutamyl residues have been mapped with a multitude of protein chemistry and proteomics approaches. Although polyglutamylation was initially considered a tubulin-specific modification, the recent studies have uncovered a calmodulin-dependent glutamylase, SidJ. Nano-electrospray ionization (ESI) proteomic approaches have identified quantifiable polyglutamylated sites in specific substrates. Indeed, conjugated glutamylated peptides were used in nano-liquid chromatography gradient delivery due to their relative hydrophobicity for their tandem mass spectrometry (MS/MS) characterization. The recent polyglutamylation characterization has revealed three major sites: E445 in α-tubulin, E435 in ß-tubulin, and E860 in SdeA. In this review, we have summarized the progress made using proteomic approaches for large-scale detection of polyglutamylated peptides, including biology and analysis.

One-pot universal NicE-seq: all enzymatic downstream processing of 4% formaldehyde crosslinked cells for chromatin accessibility genomics.

BACKGROUND: Accessible chromatin landscape allows binding of transcription factors, and remodeling of promoter and enhancer elements during development. Chromatin accessibility along with integrated multiomics approaches have been used for determining molecular subtypes of cancer in patient samples. RESULTS: One-pot Universal NicE-seq (One-pot UniNicE-seq) is an improved accessible chromatin profiling method that negate DNA purification and incorporate sonication free enzymatic fragmentation before library preparation and is suited to a variety of mammalian cells. One-pot UniNicE-seq is versatile, capable of profiling 4% formaldehyde fixed chromatin in as low as 25 fixed cells. Accessible chromatin profile is more efficient on formaldehyde-fixed cells using one-pot UniNicE-seq compared to Tn5 transposon mediated methods, demonstrating its versatility. CONCLUSION: One-pot UniNicE-seq allows the entire process of accessible chromatin labeling and enrichment in one pot at 4% formaldehyde cross-linking conditions. It doesn't require enzyme titration, compared to other technologies, since accessible chromatin is labelled with 5mC incorporation and deter degradation by nicking enzyme, thus opening the possibility for automation.

Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA.

Bisulfite sequencing detects 5mC and 5hmC at single-base resolution. However, bisulfite treatment damages DNA, which results in fragmentation, DNA loss, and biased sequencing data. To overcome these problems, enzymatic methyl-seq (EM-seq) was developed. This method detects 5mC and 5hmC using two sets of enzymatic reactions. In the first reaction, TET2 and T4-BGT convert 5mC and 5hmC into products that cannot be deaminated by APOBEC3A. In the second reaction, APOBEC3A deaminates unmodified cytosines by converting them to uracils. Therefore, these three enzymes enable the identification of 5mC and 5hmC. EM-seq libraries were compared with bisulfite-converted DNA, and each library type was ligated to Illumina adaptors before conversion. Libraries were made using NA12878 genomic DNA, cell-free DNA, and FFPE DNA over a range of DNA inputs. The 5mC and 5hmC detected in EM-seq libraries were similar to those of bisulfite libraries. However, libraries made using EM-seq outperformed bisulfite-converted libraries in all specific measures examined (coverage, duplication, sensitivity, etc.). EM-seq libraries displayed even GC distribution, better correlations across DNA inputs, increased numbers of CpGs within genomic features, and accuracy of cytosine methylation calls. EM-seq was effective using as little as 100 pg of DNA, and these libraries maintained the described advantages over bisulfite sequencing. EM-seq library construction, using challenging samples and lower DNA inputs, opens new avenues for research and clinical applications.

Aba-Seq: High-Resolution Enzymatic Mapping of Genomic 5-Hydroxymethylcytosine.

Aba-Seq (DNA modification-dependent restriction endonuclease AbaSI coupled with sequencing) provides a cost-effective and non-chemical based method for the high-resolution mapping of genomic 5-hydroxymethylcytosine (5hmC). The high specificity of the AbaSI enzyme allows sensitive detection of 5hmC in the genome. Here, we describe the Aba-Seq technique that was used for the high-resolution mapping of 5hmC in the genome of mouse embryonic stem cells (E14).

Nondestructive enzymatic deamination enables single-molecule long-read amplicon sequencing for the determination of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution.

The predominant methodology for DNA methylation analysis relies on the chemical deamination by sodium bisulfite of unmodified cytosine to uracil to permit the differential readout of methylated cytosines. Bisulfite treatment damages the DNA, leading to fragmentation and loss of long-range methylation information. To overcome this limitation of bisulfite-treated DNA, we applied a new enzymatic deamination approach, termed enzymatic methyl-seq (EM-seq), to long-range sequencing technologies. Our methodology, named long-read enzymatic modification sequencing (LR-EM-seq), preserves the integrity of DNA, allowing long-range methylation profiling of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) over multikilobase length of genomic DNA. When applied to known differentially methylated regions (DMRs), LR-EM-seq achieves phasing of >5 kb, resulting in broader and better defined DMRs compared with that previously reported. This result showed the importance of phasing methylation for biologically relevant questions and the applicability of LR-EM-seq for long-range epigenetic analysis at single-molecule and single-nucleotide resolution.

Mapping of polyglutamylation in tubulins using nanoLC-ESI-MS/MS.

Tubulin polyglutamylation is a polymeric modification that extends from the carboxyl-terminus of tubulins. Molecular description of amino acids and their branching polyglutamyls is a hallmark of tubulin in microtubules. There are different chemical approaches for detecting these polymeric structures, mostly reported prior to development of nESI peptide analysis. Here we demonstrate a novel and simple approach to detect shared regions of amino acid ions from tubulin polyglutamylated peptides in nanoLC-MS/MS. This involves two parallel in gel digestions with trypsin and subtilisin followed by mapping of di- and triglutamyl modifications of α- and ß-tubulins using a routine proteomics assay. We present three levels of information: i) identification of proteomics MS/MS data, ii) description of internal fragment ion series common across digests, and iii) extracted ion chromatograms mapped relative to retention time standards for confirmation of relative hydrophobicity values. Our nanoLC assay positive ion ESI detects up to 3 conjugated glutamates in tubulins. We implemented an analytical column only bottom up approach that characterizes molecular features of polyglutamylated tubulins.

Universal NicE-seq for high-resolution accessible chromatin profiling for formaldehyde-fixed and FFPE tissues.

Accessible chromatin plays a central role in gene expression and chromatin architecture. Current accessible chromatin approaches depend on limited digestion/cutting and pasting adaptors at the accessible DNA, thus requiring additional materials and time for optimization. Universal NicE-seq (UniNicE-seq) is an improved accessible chromatin profiling method that negates the optimization step and is suited to a variety of mammalian cells and tissues. Addition of 5-methyldeoxycytidine triphosphate during accessible chromatin labeling and an on-bead library making step substantially improved the signal to noise ratio while protecting the accessible regions from repeated nicking in cell lines, mouse T cells, mouse kidney, and human frozen tissue sections. We also demonstrate one tube UniNicE-seq for the FFPE tissue section for direct NGS library preparation without sonication and DNA purification steps. These refinements allowed reliable mapping of accessible chromatin for high-resolution genomic feature studies.

Visualization and Sequencing of Accessible Chromatin Reveals Cell Cycle and Post-HDAC inhibitor Treatment Dynamics.

Chromatin accessibility is a predictor of gene expression, cell division, and cell type specificity. NicE-viewSeq (Nicking Enzyme-assisted viewing and Sequencing) allows accessible chromatin visualization and sequencing with overall lower mitochondrial DNA and duplicated sequences interference relative to ATAC-see. Using NicE-viewSeq, we interrogated the accessibility of chromatin in a cell cycle (G1, S, and G2/M)-specific manner using mammalian cells. Despite DNA replication and subsequent condensation of chromatin to chromosomes, chromatin accessibility remained generally preserved with minimal subtle alterations. Genome-wide alteration of chromatin accessibility within TSS and enhancer elements gradually decreased as cells progressed from G1 to G2M, with distinct differential accessibility near consensus transcription factors sites. Inhibition of histone deacetylases promoted accessible chromatin within gene bodies, correlating with apoptotic gene expression. In addition, reduced chromatin accessibility for the MYC oncogene pathway correlated with downregulation of pertinent genes. Surprisingly, repetitive RNA loci expression remained unaltered following histone acetylation-mediated increased accessibility. Therefore, we suggest that subtle changes in chromatin accessibility are a prerequisite during the cell cycle and histone deacetylase inhibitor-mediated therapeutics.

Heterochromatin establishment during early mammalian development is regulated by pericentromeric RNA and characterized by non-repressive H3K9me3.

Following fertilization in mammals, the gametes are reprogrammed to create a totipotent zygote, a process that involves de novo establishment of chromatin domains. A major feature occurring during preimplantation development is the dramatic remodelling of constitutive heterochromatin, although the functional relevance of this is unknown. Here, we show that heterochromatin establishment relies on the stepwise expression and regulated activity of SUV39H enzymes. Enforcing precocious acquisition of constitutive heterochromatin results in compromised development and epigenetic reprogramming, which demonstrates that heterochromatin remodelling is essential for natural reprogramming at fertilization. We find that de novo H3K9 trimethylation (H3K9me3) in the paternal pronucleus after fertilization is catalysed by SUV39H2 and that pericentromeric RNAs inhibit SUV39H2 activity and reduce H3K9me3. De novo H3K9me3 is initially non-repressive for gene expression, but instead bookmarks promoters for compaction. Overall, we uncover the functional importance for the restricted transmission of constitutive heterochromatin during reprogramming and a non-repressive role for H3K9me3.

The microtubule-associated histone methyltransferase SET8, facilitated by transcription factor LSF, methylates α-tubulin.

Microtubules are cytoskeletal structures critical for mitosis, cell motility, and protein and organelle transport and are a validated target for anticancer drugs. However, how tubulins are regulated and recruited to support these distinct cellular processes is incompletely understood. Posttranslational modifications of tubulins are proposed to regulate microtubule function and dynamics. Although many of these modifications have been investigated, only one prior study reports tubulin methylation and an enzyme responsible for this methylation. Here we used in vitro radiolabeling, MS, and immunoblotting approaches to monitor protein methylation and immunoprecipitation, immunofluorescence, and pulldown approaches to measure protein-protein interactions. We demonstrate that N-lysine methyltransferase 5A (KMT5A or SET8/PR-Set7), which methylates lysine 20 in histone H4, bound α-tubulin and methylated it at a specific lysine residue, Lys311 Furthermore, late SV40 factor (LSF)/CP2, a known transcription factor, bound both α-tubulin and SET8 and enhanced SET8-mediated α-tubulin methylation in vitro In addition, we found that the ability of LSF to facilitate this methylation is countered by factor quinolinone inhibitor 1 (FQI1), a specific small-molecule inhibitor of LSF. These findings suggest the general model that microtubule-associated proteins, including transcription factors, recruit or stimulate protein-modifying enzymes to target tubulins. Moreover, our results point to dual functions for SET8 and LSF not only in chromatin regulation but also in cytoskeletal modification.

Protect-seq: genome-wide profiling of nuclease inaccessible domains reveals physical properties of chromatin.

In metazoan cell nuclei, heterochromatin constitutes large chromatin domains that are in close contact with the nuclear lamina. These heterochromatin/lamina-associated domains (LADs) domains are difficult to profile and warrants a simpler and direct method. Here we report a new method, Protect-seq, aimed at identifying regions of heterochromatin via resistance to nuclease degradation followed by next-generation sequencing (NGS). We performed Protect-seq on the human colon cancer cell line HCT-116 and observed overlap with previously curated LADs. We provide evidence that these protected regions are enriched for and can distinguish between the repressive histone modification H3K9me3, H3K9me2 and H3K27me3. Moreover, in human cells the loss of H3K9me3 leads to an increase in chromatin accessibility and loss of Protect-seq signal. For further validation, we performed Protect-seq in the fibrosarcoma cell line HT1080 and found a similar correlation with previously curated LADs and repressive histone modifications. In sum, Protect-seq is an efficient technique that allows rapid identification of nuclease resistant chromatin, which correlate with heterochromatin and radial positioning.