The landscape of histone modifications across 1% of the human genome in five human cell lines

Epigenetic biomarkers indicate islet cell death in xenotransplantation

BACKGROUND

Xenotransplantation of porcine islets has emerged in recent decades as a potential treatment for type 1 diabetes (T1D). Current methods of detection, indicative of successful engraftment, occur downstream of actual islet death. Epigenetic biomarkers can be detected in circulating cell-free DNA (cfDNA) to provide an earlier indication of graft dysfunction.

AIMS

The present study identified a biomarker of islet death using differential methylation of the insulin gene, INS, originating from β-cells in porcine islets.

MATERIALS & METHODS

Pyrosequencing primers specific for porcine INS were designed to quantify hypomethylation along 12 cysteine-guanine dinucleotide (CpG) sites, including three sites in the cyclic adenosine monophosphate (cAMP) response element (CRE) binding protein 2 (CRE2) binding region of the 5' untranslated region (UTR) and nine sites within intron 2.

RESULTS

PCR amplification of bisulfite-converted DNA combined with pyrosequencing data support the conclusion that hypomethylated porcine INS is specific to islet origin.

CONCLUSION

Moreover, the results of this study indicate a highly specific epigenetic biomarker, capable of detecting a single islet, supporting the measurement of cfDNA as a biomarker for transplanted islet death. Defining the epigenetic characteristics of porcine-derived islets within cfDNA will be crucial to develop a better understanding of graft survival immunology for transplantation.

A comprehensive review of computational prediction of genome-wide features

There are significant correlations among different types of genetic, genomic and epigenomic features within the genome. These correlations make the in silico feature prediction possible through statistical or machine learning models. With the accumulation of a vast amount of high-throughput data, feature prediction has gained significant interest lately, and a plethora of papers have been published in the past few years. Here we provide a comprehensive review on these published works, categorized by the prediction targets, including protein binding site, enhancer, DNA methylation, chromatin structure and gene expression. We also provide discussions on some important points and possible future directions.

Comparative Proteomic Analysis of Histone Modifications upon Acridone Derivative -Induced CCRF-CEM Cells by Data Independent Acquisition

The lead compound acridone derivative 8a showed potent antiproliferative activity by inducing DNA damage through direct stacking with DNA bases and triggering ROS in CCRF-CEM cells. To define the chromatin alterations during DNA damage sensing and repair, a detailed quantitative map of single and coexisting histone post-translational modifications (PTMs) in CCRF-CEM cells affected by 8a was performed by the Data Independent Acquisition (DIA) method on QE-plus. A total of 79 distinct and 164 coexisting histone PTMs were quantified, of which 16 distinct histone PTMs were significantly altered when comparing 8a-treated cells with vehicle control cells. The changes in histone PTMs were confirmed by Western blotting analysis for three H3 and one H4 histone markers. The up-regulated dimethylation on H3K9, H3K36, and H4K20 implied that CCRF-CEM cells might accelerate DNA damage repair to counteract the DNA lesion induced by 8a, which was verified by an increment in the 53BP1 foci localization at the damaged DNA. Most of the significantly altered PTMs were involved in transcriptional regulation, including down-regulated acetylation on H3K18, H3K27, and H3K122, and up-regulated di- and trimethylation on H3K9 and H3K27. This transcription-silencing phenomenon was associated with G2/M cell cycle arrest after 8a treatment by flow cytometry. This study shows that the DIA proteomics strategy provides a sensitive and accurate way to characterize the coexisting histone PTMs changes and their cross-talk in CCRF-CEM cells after 8a treatment. Specifically, histone PTMs rearrange transcription-silencing, and cell cycle arrest DNA damage repair may contribute to the mechanism of epigenetic response affected by 8a.

Regulation of human trophoblast syncytialization by histone demethylase LSD1

A successful pregnancy is critically dependent upon proper placental development and function. During human placentation, villous cytotrophoblast (CTB) progenitors differentiate to form syncytiotrophoblasts (SynTBs), which provide the exchange surface between the mother and fetus and secrete hormones to ensure proper progression of pregnancy. However, epigenetic mechanisms that regulate SynTB differentiation from CTB progenitors are incompletely understood. Here, we show that lysine-specific demethylase 1 (LSD1; also known as KDM1A), a histone demethylase, is essential to this process. LSD1 is expressed both in CTB progenitors and differentiated SynTBs in first-trimester placental villi; accordingly, expression in SynTBs is maintained throughout gestation. Impairment of LSD1 function in trophoblast progenitors inhibits induction of endogenous retrovirally encoded genes SYNCYTIN1/endogenous retrovirus group W member 1, envelope (ERVW1) and SYNCYTIN2/endogenous retrovirus group FRD member 1, envelope (ERVFRD1), encoding fusogenic proteins critical to human trophoblast syncytialization. Loss of LSD1 also impairs induction of chorionic gonadotropin α (CGA) and chorionic gonadotropin β (CGB) genes, which encode α and β subunits of human chorionic gonadotrophin (hCG), a hormone essential to modulate maternal physiology during pregnancy. Mechanistic analyses at the endogenous ERVW1, CGA, and CGB loci revealed a regulatory axis in which LSD1 induces demethylation of repressive histone H3 lysine 9 dimethylation (H3K9Me2) and interacts with transcription factor GATA2 to promote RNA polymerase II (RNA-POL-II) recruitment and activate gene transcription. Our study reveals a novel LSD1-GATA2 axis, which regulates human trophoblast syncytialization.

Nuclear Pore Proteins in Regulation of Chromatin State

Nuclear pore complexes (NPCs) are canonically known to regulate nucleocytoplasmic transport. However, research efforts over the last decade have demonstrated that NPCs and their constituent nucleoporins (Nups) also interact with the genome and perform important roles in regulation of gene expression. It has become increasingly clear that many Nups execute these roles specifically through regulation of chromatin state, whether through interactions with histone modifiers and downstream changes in post-translational histone modifications, or through relationships with chromatin-remodeling proteins that can result in physical changes in nucleosome occupancy and chromatin compaction. This review focuses on these findings, highlighting the functional connection between NPCs/Nups and regulation of chromatin structure, and how this connection can manifest in regulation of transcription.

DNA methylation, environmental exposures and early embryo development

The first crucial step in the developmental program occurs during pre-implantation, the time after the oocyte has been fertilized and before the embryo implants in the uterus. This period represents a vulnerable window as the epigenome undergoes dynamic changes in DNA methylation profiles. Alterations in the early embryonic reprogramming wave can impair DNA methylation patterns and induce permanent changes to the developmental program, leading to the onset of adverse health outcomes in offspring. Although there is an increasing body of evidence indicating that harmful exposures during pre-implantation embryo development can trigger lasting epigenetic alterations in offspring, the mechanisms are still not fully understood. Since physiological or pathological changes in DNA methylation can occur as a response to environmental cues, proper environmental milieu plays a critical role in the success of embryonic development. In this review, we depict the mechanisms behind the embryonic epigenetic reprogramming of DNA methylation and highlight how maternal environmental stressors (e.g., alcohol, heat stress, nutrient availability) during pre-implantation and assisted reproductive technology procedures affect development and DNA methylation marks.

Functions and mechanisms of epigenetic inheritance in animals

The idea that epigenetic determinants such as DNA methylation, histone modifications or RNA can be passed to the next generation through meiotic products (gametes) is long standing. Such meiotic epigenetic inheritance (MEI) is fairly common in yeast, plants and nematodes, but its extent in mammals has been much debated. Advances in genomics techniques are now driving the profiling of germline and zygotic epigenomes, thereby improving our understanding of MEI in diverse species. Whereas the role of DNA methylation in MEI remains unclear, insights from genome-wide studies suggest that a previously underappreciated fraction of mammalian genomes bypass epigenetic reprogramming during development. Notably, intergenerational inheritance of histone modifications, tRNA fragments and microRNAs can affect gene regulation in the offspring. It is important to note that MEI in mammals rarely constitutes transgenerational epigenetic inheritance (TEI), which spans multiple generations. In this Review, we discuss the examples of MEI in mammals, including mammalian epigenome reprogramming, and the molecular mechanisms of MEI in vertebrates in general. We also discuss the implications of the inheritance of histone modifications and small RNA for embryogenesis in metazoans, with a particular focus on insights gained from genome-wide studies.

Transcriptional control by enhancers and enhancer RNAs

The regulation of gene expression is a fundamental cellular process and its misregulation is a key component of disease. Enhancers are one of the most salient regulatory elements in the genome and help orchestrate proper spatiotemporal gene expression during development, in homeostasis, and in response to signaling. Notably, molecular aberrations at enhancers, such as translocations and single nucleotide polymorphisms, are emerging as an important source of human variation and susceptibility to disease. Herein we discuss emerging paradigms addressing how genes are regulated by enhancers, common features of active enhancers, and how non-coding enhancer RNAs (eRNAs) can direct gene expression programs that underlie cellular phenotypes. We survey the current evidence, which suggests that eRNAs can bind to transcription factors, mediate enhancer-promoter interactions, influence RNA Pol II elongation, and act as decoys for repressive cofactors. Furthermore, we discuss current methodologies for the identification of eRNAs and novel approaches to elucidate their functions.

Effects of histone deacetylase inhibitors on ATP-binding cassette transporters in lung cancer A549 and colorectal cancer HCT116 cells

Histone deacetylase (HDAC) inhibitors and DNA alkylators are effective components used in combination chemotherapy. In the present study, the effects of HDAC inhibitors on the expression of ATP-binding cassette (ABC) transporters were investigated. It was observed that HDAC inhibitors induced the expression of multidrug-resistant ABC transporters differently in lung cancer A549 cells than in colorectal cancer HCT116 cells. In these two cell lines, the HDAC inhibitors suberoylanilide hydroxamic acid (SAHA) and trichostatin A (TSA) significantly increased ABCB1 expression at the mRNA and protein levels, whereas they had no evident effect on ABCG2 protein expression. SAHA and TSA decreased ABCG2 mRNA expression in A549 cells and had no evident effect on ABCG2 mRNA expression in HCT116 cells. Notably, SAHA and TSA increased the mRNA expression levels of ABCC5, ABCC6, ABCC10, ABCC11 and ABCC12, as well as the protein expression levels of ABCC2, ABCC10 and ABCC12. By contrast, these inhibitors decreased the mRNA expression levels of ABCC1, ABCC2, ABCC3 and ABCC4, as well as the expression of ABCC1 and ABCC3 proteins. Furthermore, SAHA and TSA were found to downregulate HDAC3 and HDAC4, but not HDAC1 and HDAC2. Taken together, the results suggested that HDAC inhibitors work synergistically with DNA alkylators, in part, due to the inhibitory effect of these inhibitors on ABCC1 expression, which translocates these alkylators from inside to outside of cancer cells. These results further suggested the possibility of antagonism when HDAC inhibitors are combined with anthracyclines and other ABCB1 drug ligands in chemotherapy.

Antiproliferative effect of the histone demethylase inhibitor GSK-J4 in chondrosarcomas

Chondrosarcoma (CS) is the second most common malignant bone sarcoma. Its treatment remains an issue, because this tumor is radio- and chemo-resistant. In the present study, we investigated the antitumoral potential of GSK-J4, a small molecule described as an inhibitor of histone demethylases UTX and JMJD3 (KDM6A and KDM6B), alone or in combination with cisplatin in CSs. Human CS-derived cell lines were treated with GSK-J4 in the presence or not of cisplatin. Survival curves were established and cell proliferation and cycle were evaluated by flow cytometry using dividing cell tracking technique utilizing carboxyfluorescein succinimidyl ester labeling, or DNA staining by propidium iodide. Apoptosis and senescence were also investigated. GSK-J4 decreased proliferation of CS cells. Additionally, it induced apoptosis in CH2879 and JJ012 cells, but not in SW1353 CSs. In addition, its association with cisplatin decreased cell proliferation more than drugs alone, whereas it did not increase apoptosis compared to cisplatin alone. Interestingly, GSK-J4 alone as well as in association with cisplatin did not affect chondrocyte survival or proliferation. In conclusion, this study suggests that demethylase inhibitors may be useful in improving therapy for CS in reducing its proliferation.

Genetic Contexts Characterize Dynamic Histone Modification Patterns Among Cell Types

Histone modifications play critical roles in mammalian development, regulating chromatin structure and gene expression. Dynamic histone modifications among cell types have been shown to associate with changes in mammalian development. However, how to quantitatively measure the histone modification alterations and how histone modifications vary across cell types under different genetic contexts remain largely unexplored and whether these changes are related to the primary DNA sequence remains limited. Here, we employed an entropy-based method to measure histone modification alterations in six definite genomic regions across five cell types and identified lineage-specific histone modification genes. We observed that histone modification alterations prefer to enrich in 5'-UTR exons, and also in 3'-UTR exons and its downstream. Then we built a model to predict the histone modification patterns from the primary DNA sequence. We found that the frequencies of k-mer sequence compositions are predictive of histone modification patterns, suggesting that the primary DNA sequence correlated with the histone modification alterations among cell types. Additionally, the lineage-specific histone modification genes display a higher conservation and lower GC-content. Together, we performed a systematic analysis for histone modification alterations and demonstrated how to identify genomic region-specific elements of epigenetic and genetic regulation and histone modification patterns across different cell types.

DeepTACT: predicting 3D chromatin contacts via bootstrapping deep learning

Interactions between regulatory elements are of crucial importance for the understanding of transcriptional regulation and the interpretation of disease mechanisms. Hi-C technique has been developed for genome-wide detection of chromatin contacts. However, unless extremely deep sequencing is performed on a very large number of input cells, which is technically limited and expensive, current Hi-C experiments do not have high enough resolution to resolve contacts between regulatory elements. Here, we develop DeepTACT, a bootstrapping deep learning model, to integrate genome sequences and chromatin accessibility data for the prediction of chromatin contacts between regulatory elements. DeepTACT can infer not only promoter-enhancer interactions, but also promoter-promoter interactions. In tests based on promoter capture Hi-C data, DeepTACT shows better performance over existing methods. DeepTACT analysis also identifies a class of hub promoters, which are correlated with transcriptional activation across cell lines, enriched in housekeeping genes, functionally related to fundamental biological processes, and capable of reflecting cell similarity. Finally, the utility of chromatin contacts in the study of human diseases is illustrated by the association of IFNA2 to coronary artery disease via an integrative analysis of GWAS data and interactions predicted by DeepTACT.

Potential epigenomic co-management in rare diseases and epigenetic therapy

The purpose of this review is to highlight the impact of the alternative splicing process on human disease. Epigenetic regulation determines not only what parts of the genome are expressed but also how they are spliced. The recent progress in the field of epigenetics has important implications for the study of rare diseases. The role of epigenetics in rare diseases is a key issue in molecular physiology and medicine because not only rare diseases can benefit from epigenetic research, but can also provide useful principles for other common and complex disorders such as cancer, cardiovascular, type 2 diabetes, obesity, and neurological diseases. Predominantly, epigenetic modifications include DNA methylation, histone modification, and RNA-associated silencing. These modifications in the genome regulate numerous cellular activities. Disruption of epigenetic regulation process can contribute to the etiology of numerous diseases during both prenatal and postnatal life. Here, I discuss current knowledge about this matter including some current epigenetic therapies and future directions in the field by emphasizing on the RNA-based therapy via antisense oligonucleotides to correct splicing defects.

Transcriptional repression by androgen receptor: roles in castration-resistant prostate cancer

Androgen receptor (AR), a hormonal transcription factor, plays important roles during prostate cancer progression and is a key target for therapeutic interventions. While androgen-deprivation therapies are initially successful in regressing prostate tumors, the disease ultimately comes back as castration-resistant prostate cancer (CRPC) or at the late stage as neuroendocrine prostate cancer (NEPC). CRPC remains largely dependent on hyperactive AR signaling in the milieu of low androgen, while NEPC is negative of AR expression but positive of many AR-repressed genes. Recent technological advances in genome-wide analysis of transcription factor binding sites have revealed an unprecedented set of AR target genes. In addition to its well-known function in activating gene expression, AR is increasingly known to also act as a transcriptional repressor. Here, we review the molecular mechanisms by which AR represses gene expression. We also summarize AR-repressed genes that are aberrantly upregulated in CRPC and NEPC and represent promising targets for therapeutic intervention.

A new sequence logo plot to highlight enrichment and depletion

BACKGROUND

Sequence logo plots have become a standard graphical tool for visualizing sequence motifs in DNA, RNA or protein sequences. However standard logo plots primarily highlight enrichment of symbols, and may fail to highlight interesting depletions. Current alternatives that try to highlight depletion often produce visually cluttered logos.

RESULTS

We introduce a new sequence logo plot, the EDLogo plot, that highlights both enrichment and depletion, while minimizing visual clutter. We provide an easy-to-use and highly customizable R package Logolas to produce a range of logo plots, including EDLogo plots. This software also allows elements in the logo plot to be strings of characters, rather than a single character, extending the range of applications beyond the usual DNA, RNA or protein sequences. And the software includes new Empirical Bayes methods to stabilize estimates of enrichment and depletion, and thus better highlight the most significant patterns in data. We illustrate our methods and software on applications to transcription factor binding site motifs, protein sequence alignments and cancer mutation signature profiles.

CONCLUSIONS

Our new EDLogo plots and flexible software implementation can help data analysts visualize both enrichment and depletion of characters (DNA sequence bases, amino acids, etc.) across a wide range of applications.

The conservation landscape of the human ribosomal RNA gene repeats

Ribosomal RNA gene repeats (rDNA) encode ribosomal RNA, a major component of ribosomes. Ribosome biogenesis is central to cellular metabolic regulation, and several diseases are associated with rDNA dysfunction, notably cancer, However, its highly repetitive nature has severely limited characterization of the elements responsible for rDNA function. Here we make use of phylogenetic footprinting to provide a comprehensive list of novel, potentially functional elements in the human rDNA. Complete rDNA sequences for six non-human primate species were constructed using de novo whole genome assemblies. These new sequences were used to determine the conservation profile of the human rDNA, revealing 49 conserved regions in the rDNA intergenic spacer (IGS). To provide insights into the potential roles of these conserved regions, the conservation profile was integrated with functional genomics datasets. We find two major zones that contain conserved elements characterised by enrichment of transcription-associated chromatin factors, and transcription. Conservation of some IGS transcripts in the apes underpins the potential functional significance of these transcripts and the elements controlling their expression. Our results characterize the conservation landscape of the human IGS and suggest that noncoding transcription and chromatin elements are conserved and important features of this unique genomic region.

The Histone H3 Lysine 4 Presenter WDR5 as an Oncogenic Protein and Novel Epigenetic Target in Cancer

The histone H3 lysine 4 (H3K4) presenter WDR5 forms protein complexes with H3K4 methyltransferases MLL1-MLL4 and binding partner proteins including RBBP5, ASH2L, and DPY30, and plays a key role in histone H3K4 trimethylation, chromatin remodeling, transcriptional activation of target genes, normal biology, and diseases such as MLL-rearranged leukemia. By forming protein complexes with other proteins such as Myc, WDR5 induces transcriptional activation of key oncogenes, tumor cell cycle progression, DNA replication, cell proliferation, survival, tumor initiation, progression, invasion, and metastasis of cancer of a variety of organ origins. Several small molecule MLL/WDR5 protein-protein interaction inhibitors, such as MM-401, MM-589, WDR5-0103, Piribedil, and OICR-9429, have been confirmed to reduce H3K4 trimethylation, oncogenic gene expression, cell cycle progression, cancer cell proliferation, survival and resistance to chemotherapy without general toxicity to normal cells. Derivatives of the MLL/WDR5 interaction inhibitors with improved pharmacokinetic properties and in vivo bioavailability are expected to have the potential to be trialed in cancer patients.

Loss of SETDB1 decompacts the inactive X chromosome in part through reactivation of an enhancer in the IL1RAPL1 gene

Background

The product of dosage compensation in female mammals is the inactive X chromosome (Xi). Xi facultative heterochromatin is organized into two different types, one of which is defined by histone H3 trimethylated at lysine 9 (H3K9me3). The rationale for this study was to assess SET domain bifurcated 1 (SETDB1) as a candidate for maintaining this repressive modification at the human Xi.

Results

Here, we show that loss of SETDB1 does not result in large-scale H3K9me3 changes at the Xi, but unexpectedly we observed striking decompaction of the Xi territory. Close examination revealed a 0.5 Mb region of the Xi that transitioned from H3K9me3 heterochromatin to euchromatin within the 3′ end of the IL1RAPL1 gene that is part of a common chromosome fragile site that is frequently deleted or rearranged in patients afflicted with intellectual disability and other neurological ailments. Centrally located within this interval is a powerful enhancer adjacent to an ERVL-MaLR element. In the absence of SETDB1, the enhancer is reactivated on the Xi coupled with bidirectional transcription from the ERVL-MaLR element. Xa deletion of the enhancer/ERVL-MaLR resulted in loss of full-length IL1RAPL1 transcript in cis, coupled with trans decompaction of the Xi chromosome territory, whereas Xi deletion increased detection of full-length IL1RAPL1 transcript in trans, but did not impact Xi compaction.

Conclusions

These data support a critical role for SETDB1 in maintaining the ERVL-MaLR element and adjacent enhancer in the 3′ end of the IL1RAPL1 gene in a silent state to facilitate Xi compaction.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0218-9) contains supplementary material, which is available to authorized users.

Glucocorticoid receptor recruits to enhancers and drives activation by motif-directed binding

Glucocorticoids are potent steroid hormones that regulate immunity and metabolism by activating the transcription factor (TF) activity of glucocorticoid receptor (GR). Previous models have proposed that DNA binding motifs and sites of chromatin accessibility predetermine GR binding and activity. However, there are vast excesses of both features relative to the number of GR binding sites. Thus, these features alone are unlikely to account for the specificity of GR binding and activity. To identify genomic and epigenetic contributions to GR binding specificity and the downstream changes resultant from GR binding, we performed hundreds of genome-wide measurements of TF binding, epigenetic state, and gene expression across a 12-h time course of glucocorticoid exposure. We found that glucocorticoid treatment induces GR to bind to nearly all pre-established enhancers within minutes. However, GR binds to only a small fraction of the set of accessible sites that lack enhancer marks. Once GR is bound to enhancers, a combination of enhancer motif composition and interactions between enhancers then determines the strength and persistence of GR binding, which consequently correlates with dramatic shifts in enhancer activation. Over the course of several hours, highly coordinated changes in TF binding and histone modification occupancy occur specifically within enhancers, and these changes correlate with changes in the expression of nearby genes. Following GR binding, changes in the binding of other TFs precede changes in chromatin accessibility, suggesting that other TFs are also sensitive to genomic features beyond that of accessibility.

Modulation of matrix metabolism by ATP-citrate lyase in articular chondrocytes

Certain dysregulated chondrocyte metabolic adaptive responses such as decreased activity of the master regulator of energy metabolism AMP-activated protein kinase (AMPK) promote osteoarthritis (OA). Metabolism intersects with epigenetic and transcriptional responses. Hence, we studied chondrocyte ATP-citrate lyase (ACLY), which generates acetyl-CoA from mitochondrial-derived citrate, and modulates acetylation of histones and transcription factors. We assessed ACLY in normal and OA human knee chondrocytes and cartilages by Western blotting and immunohistochemistry, and quantified acetyl-CoA fluorometrically. We examined histone and transcription factor lysine acetylation by Western blotting, and assessed histone H3K9 and H3K27 occupancy of iNOS, MMP3, and MMP13 promoters by chromatin immunoprecipitation (ChIP) and quantitative PCR (qPCR). We analyzed iNOS, MMP3, MMP13, aggrecan (ACAN), and Col2a1 gene expression by RT-qPCR. Glucose availability regulated ACLY expression and function, nucleocytosolic acetyl-CoA, and histone acetylation. Human knee OA chondrocytes exhibited increased ACLY activation (assessed by Ser-455 phosphorylation), associated with increased H3K9 and H3K27 acetylation. Inhibition of ACLY attenuated IL-1β-induced transcription of iNOS, MMP3, and MMP13 by suppressing acetylation of p65 NF-κB, H3K9, and H3K27, blunted release of NO, MMP3, and MMP13, and also reduced SOX9 acetylation that promoted SOX9 nuclear translocation, leading to increased aggrecan and Col2a1 mRNA expression. ACLY is a novel player involved in regulation of cartilage matrix metabolism. Increased ACLY activity in OA chondrocytes increased nucleocytosolic acetyl-CoA, leading to increased matrix catabolism via dysregulated histone and transcription factor acetylation. Pharmacologic ACLY inhibition in OA chondrocytes globally reverses these changes and stimulates matrix gene expression and AMPK activation, supporting translational investigation in OA.

Differential effects of histone deacetylase inhibitors on cellular drug transporters and their implications for using epigenetic modifiers in combination chemotherapy

HDAC inhibitors, DNA alkylators and nucleoside analogs are effective components of combination chemotherapy. To determine a possible mechanism of their synergism, we analyzed the effects of HDAC inhibitors on the expression of drug transporters which export DNA alkylators. Exposure of PEER lymphoma T-cells to 15 nM romidepsin (Rom) resulted in 40%-50% reduction in mRNA for the drug transporter MRP1 and up to ~500-fold increase in the MDR1 mRNA within 32-48 hrs. MRP1 protein levels concomitantly decreased while MDR1 increased. Other HDAC inhibitors - panobinostat, belinostat and suberoylanilide hydroxamic acid (SAHA) - had similar effects on these transporters. The protein level of MRP1 correlated with cellular resistance to busulfan and chlorambucil, and Rom exposure sensitized cells to these DNA alkylators. The decrease in MRP1 correlated with decreased cellular drug export activity, and increased level of MDR1 correlated with increased export of daunorubicin. A similar decrease in the level of MRP1 protein, and increase in MDR1, were observed when mononuclear cells derived from patients with T-cell malignancies were exposed to Rom. Decreased MRP1 and increased MDR1 expressions were also observed in blood mononuclear cells from lymphoma patients who received SAHA-containing chemotherapy in a clinical trial. This inhibitory effect of HDAC inhibitors on the expression of MRP1 suggests that their synergism with DNA alkylating agents is partly due to decreased efflux of these alkylators. Our results further imply the possibility of antagonistic effects when HDAC inhibitors are combined with anthracyclines and other MDR1 drug ligands in chemotherapy.

BCL11B Drives Human Mammary Stem Cell Self-Renewal In Vitro by Inhibiting Basal Differentiation

The epithelial compartment of the mammary gland contains basal and luminal cell lineages, as well as stem and progenitor cells that reside upstream in the differentiation hierarchy. Stem and progenitor cell differentiation is regulated to maintain adult tissue and mediate expansion during pregnancy and lactation. The genetic factors that regulate the transition of cells between differentiation states remain incompletely understood. Here, we present a genome-scale method to discover genes driving cell-state specification. Applying this method, we identify a transcription factor, BCL11B, which drives stem cell self-renewal in vitro, by inhibiting differentiation into the basal lineage. To validate BCL11B's functional role, we use two-dimensional colony-forming and three-dimensional tissue differentiation assays to assess the lineage differentiation potential and functional abilities of primary human mammary cells. These findings show that BCL11B regulates mammary cell differentiation and demonstrate the utility of our proposed genome-scale strategy for identifying lineage regulators in mammalian tissues.

•

Novel strategy to identify lineage-specifying genes

•

BCL11B drives primary human mammary stem cell self-renewal in vitro

•

BCL11B inhibits basal differentiation of mammary stem cells in vitro

Factors Affecting the Immunity to Respiratory Syncytial Virus: From Epigenetics to Microbiome

Respiratory syncytial virus (RSV) is a common pathogen that infects virtually all children by 2 years of age and is the leading cause of hospitalization of infants worldwide. While most children experience mild symptoms, some children progress to severe lower respiratory tract infection. Those children with severe disease have a much higher risk of developing childhood wheezing later in life. Many risk factors are known to result in exacerbated disease, including premature birth and early age of RSV infection, when the immune system is relatively immature. The development of the immune system before and after birth may be altered by several extrinsic and intrinsic factors that could lead to severe disease predisposition in children who do not exhibit any currently known risk factors. Recently, the role of the microbiome and the resulting metabolite profile has been an area of intense study in the development of lung disease, including viral infection and asthma. This review explores both known risk factors that can lead to severe RSV-induced disease as well as emerging topics in the development of immunity to RSV and the long-term consequences of severe infection.

Coordinated activation of a cluster of MMP genes in response to UVB radiation

Ultraviolet (UV) B radiation is a dangerous environmental stressor, which can lead to photoaging, inflammation, immune suppression and tumour formation. A recent report has shown the transcriptional activation of several skin-specific genes including matrix metalloproteases (MMPs) in response to UV irradiation. Here, we use a novel human keratinocyte model, HKerE6SFM, to demonstrate that UVB activates the transcription of most members of the 11q22.3 MMP gene cluster including MMP13, MMP12, MMP3, MMP1 and MMP10. Curiously, the expression of the well-characterized UVB-inducible MMP9, which is located outside of the cluster, remains unchanged. In accordance with the increased expression of the MMP gene cluster upon UVB irradiation, RNA polymerase II showed increased occupancy at their promoters following UVB irradiation. The results also demonstrate increased acetylated histone H3K9 levels at the promoters of the MMP13, MMP12, MMP3, MMP1 and MMP10 genes. These findings suggest a coordinated transcriptional activation of genes in the MMP cluster at 11q22.3 and that acetylation of histone H3 at lysine 9 has an important role in the UVB-dependent enhancement of transcription of MMP genes in this region.

Identification and validation of regulatory SNPs that modulate transcription factor chromatin binding and gene expression in prostate cancer

Prostate cancer (PCa) is the second most common solid tumor for cancer related deaths in American men. Genome wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with the increased risk of PCa. Because most of the susceptibility SNPs are located in noncoding regions, little is known about their functional mechanisms. We hypothesize that functional SNPs reside in cell type-specific regulatory elements that mediate the binding of critical transcription factors (TFs), which in turn result in changes in target gene expression. Using PCa-specific functional genomics data, here we identify 38 regulatory candidate SNPs and their target genes in PCa. Through risk analysis by incorporating gene expression and clinical data, we identify 6 target genes (ZG16B, ANKRD5, RERE, FAM96B, NAALADL2 and GTPBP10) as significant predictors of PCa biochemical recurrence. In addition, 5 SNPs (rs2659051, rs10936845, rs9925556, rs6057110 and rs2742624) are selected for experimental validation using Chromatin immunoprecipitation (ChIP), dual-luciferase reporter assay in LNCaP cells, showing allele-specific enhancer activity. Furthermore, we delete the rs2742624-containing region using CRISPR/Cas9 genome editing and observe the drastic downregulation of its target gene UPK3A. Taken together, our results illustrate that this new methodology can be applied to identify regulatory SNPs and their target genes that likely impact PCa risk. We suggest that similar studies can be performed to characterize regulatory variants in other diseases.

Treponema denticola increases MMP-2 expression and activation in the periodontium via reversible DNA and histone modifications

Host-derived matrix metalloproteinases (MMPs) and bacterial proteases mediate destruction of extracellular matrices and supporting alveolar bone in periodontitis. The Treponema denticola dentilisin protease induces MMP-2 expression and activation in periodontal ligament (PDL) cells, and dentilisin-mediated activation of pro-MMP-2 is required for cellular fibronectin degradation. Here, we report that T. denticola regulates MMP-2 expression through epigenetic modifications in the periodontium. PDL cells were treated with epigenetic enzyme inhibitors before or after T. denticola challenge. Fibronectin fragmentation, MMP-2 expression, and activation were assessed by immunoblot, zymography, and qRT-PCR, respectively. Chromatin modification enzyme expression in T. denticola-challenged PDL cells and periodontal tissues were evaluated using gene arrays. Several classes of epigenetic enzymes showed significant alterations in transcription in diseased tissue and T. denticola-challenged PDL cells. T. denticola-mediated MMP-2 expression and activation were significantly reduced in PDL cells treated with inhibitors of aurora kinases and histone deacetylases. In contrast, DNA methyltransferase inhibitors had little effect, and inhibitors of histone acetyltransferases, methyltransferases, and demethylases exacerbated T. denticola-mediated MMP-2 expression and activation. Chronic epigenetic changes in periodontal tissues mediated by T. denticola or other oral microbes may contribute to the limited success of conventional treatment of chronic periodontitis and may be amenable to therapeutic reversal.

Structure and Chromosomal Organization of Yeast Genes Regulated by Topoisomerase II

Cellular DNA topoisomerases (topo I and topo II) are highly conserved enzymes that regulate the topology of DNA during normal genome transactions, such as DNA transcription and replication. In budding yeast, topo I is dispensable whereas topo II is essential, suggesting fundamental and exclusive roles for topo II, which might include the functions of the topo IIa and topo IIb isoforms found in mammalian cells. In this review, we discuss major findings of the structure and chromosomal organization of genes regulated by topo II in budding yeast. Experimental data was derived from short (10 min) and long term (120 min) responses to topo II inactivation in top-2 ts mutants. First, we discuss how short term responses reveal a subset of yeast genes that are regulated by topo II depending on their promoter architecture. These short term responses also uncovered topo II regulation of transcription across multi-gene clusters, plausibly by common DNA topology management. Finally, we examine the effects of deactivated topo II on the elongation of RNA transcripts. Each study provides an insight into the particular chromatin structure that interacts with the activity of topo II. These findings are of notable clinical interest as numerous anti-cancer therapies interfere with topo II activity.

Chromatin Immunoprecipitation of Skeletal Muscle Tissue

Chromatin immunoprecipitation (ChIP) is an invaluable method for studying interactions between histone proteins and genomic DNA regions and transcriptional regulation using antibodies to enrich genomic regions associated with these epitopes. Either to monitor the presence of histones with post-translational modifications at specific genomic locations or to measure transcription factor interactions with a candidate target gene, protein-DNA complexes are most commonly crosslinked using formaldehyde, which stabilizes these transient interactions. Chromatin is then fragmented to allow separation of genomic fragments bound by the histone or transcription factor of interest away from those that are unbound. Following immunoprecipitation, formaldehyde crosslinks are reversed and enriched DNA fragments are purified. While some investigators have successfully performed ChIP experiments from crosslinked skeletal muscle in cell culture, the process is relatively inefficient compared to whole tissue. This chapter provides protocols specifically designed for the crosslinking and immunoprecipitation of human skeletal muscle biopsy samples in preparation for chromatin immunoprecipitation-sequencing (ChIP-seq).

Predicting enhancers with deep convolutional neural networks

Background

With the rapid development of deep sequencing techniques in the recent years, enhancers have been systematically identified in such projects as FANTOM and ENCODE, forming genome-wide landscapes in a series of human cell lines. Nevertheless, experimental approaches are still costly and time consuming for large scale identification of enhancers across a variety of tissues under different disease status, making computational identification of enhancers indispensable.

Results

To facilitate the identification of enhancers, we propose a computational framework, named DeepEnhancer, to distinguish enhancers from background genomic sequences. Our method purely relies on DNA sequences to predict enhancers in an end-to-end manner by using a deep convolutional neural network (CNN). We train our deep learning model on permissive enhancers and then adopt a transfer learning strategy to fine-tune the model on enhancers specific to a cell line. Results demonstrate the effectiveness and efficiency of our method in the classification of enhancers against random sequences, exhibiting advantages of deep learning over traditional sequence-based classifiers. We then construct a variety of neural networks with different architectures and show the usefulness of such techniques as max-pooling and batch normalization in our method. To gain the interpretability of our approach, we further visualize convolutional kernels as sequence logos and successfully identify similar motifs in the JASPAR database.

Conclusions

DeepEnhancer enables the identification of novel enhancers using only DNA sequences via a highly accurate deep learning model. The proposed computational framework can also be applied to similar problems, thereby prompting the use of machine learning methods in life sciences.

Functional annotation of structural ncRNAs within enhancer RNAs in the human genome: implications for human disease

Enhancer RNAs (eRNAs) are a novel class of non-coding RNA (ncRNA) molecules transcribed from the DNA sequences of enhancer regions. Despite extensive efforts devoted to revealing the potential functions and underlying mechanisms of eRNAs, it remains an open question whether eRNAs are mere transcriptional noise or relevant biologically functional species. Here, we identified a catalogue of eRNAs in a broad range of human cell/tissue types and extended our understanding of eRNAs by demonstrating their multi-omic signatures. Gene Ontology (GO) analysis revealed that eRNAs play key roles in human cell identity. Furthermore, we detected numerous known and novel functional RNA structures within eRNA regions. To better characterize the cis-regulatory effects of non-coding variation in these structural ncRNAs, we performed a comprehensive analysis of the genetic variants of structural ncRNAs in eRNA regions that are associated with inflammatory autoimmune diseases. Disease-associated variants of the structural ncRNAs were disproportionately enriched in immune-specific cell types. We also identified riboSNitches in lymphoid eRNAs and investigated the potential pathogenic mechanisms by which eRNAs might function in autoimmune diseases. Collectively, our findings offer valuable insights into the function of eRNAs and suggest that eRNAs might be effective diagnostic and therapeutic targets for human diseases.

Epigenetic modifications of gene expression by lifestyle and environment

Epigenetics oftenly described as the heritable changes in gene expression independent of changes in DNA sequence. Various environmental factors such as nutrition-dietary components, lifestyle, exercise, physical activity, toxins, and other contributing factors remodel the genome either in a constructive or detrimental way. Since epigenetic changes are reversible and nutrition is one of the many epigenetic regulators that modify gene expression without changing the DNA sequence, dietary nutrients and bioactive food components contribute to epigenetic phenomena either by directly suppressing DNA methylation or histone catalyzing enzymes or by changing the availability of substrates required for enzymatic reactions. Diets that contain catechol-dominant polyphenols are reported to suppress enzyme activity and activate epigenetically silenced genes. Furthermore, several dietary nutrients play a crucial role in one-carbon metabolism including folate, cobalamin, riboflavin, pyridoxine, and methionine by directly affecting S-adenosyl-L-methionine. Soy polyphenols block DNA methyltransferases and histone deacetylases to reverse aberrant CpG island methylation. Organosulfur rich compounds such as the sulforaphane found in broccoli appear to normalize DNA methylation and activate miR-140 expression, which represses SOX9 and ALDH1 and decreases tumor growth. The purpose of this short communication is to overview the epigenetic regulatory mechanisms of diet and other environmental factors. We discuss the epigenetic contributions of dietary components with a particular focus on nutritional polyphenols and flavonoids as epigenetic mediators that modify epigenetic tags and control gene expression. These mechanisms provide new insights to better understand the influence of dietary nutrients on epigenetic modifications and gene expression.

Identification of active miRNA promoters from nuclear run-on RNA sequencing

The genome-wide identification of microRNA transcription start sites (miRNA TSSs) is essential for understanding how miRNAs are regulated in development and disease. In this study, we developed mirSTP (mirna transcription Start sites Tracking Program), a probabilistic model for identifying active miRNA TSSs from nascent transcriptomes generated by global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq). MirSTP takes advantage of characteristic bidirectional transcription signatures at active TSSs in GRO/PRO-seq data, and provides accurate TSS prediction for human intergenic miRNAs at a high resolution. MirSTP performed better than existing generalized and experiment specific methods, in terms of the enrichment of various promoter-associated marks. MirSTP analysis of 27 human cell lines in 183 GRO-seq and 28 PRO-seq experiments identified TSSs for 480 intergenic miRNAs, indicating a wide usage of alternative TSSs. By integrating predicted miRNA TSSs with matched ENCODE transcription factor (TF) ChIP-seq data, we connected miRNAs into the transcriptional circuitry, which provides a valuable source for understanding the complex interplay between TF and miRNA. With mirSTP, we not only predicted TSSs for 72 miRNAs, but also identified 12 primary miRNAs with significant RNA polymerase pausing alterations after JQ1 treatment; each miRNA was further validated through BRD4 binding to its predicted promoter. MirSTP is available at http://bioinfo.vanderbilt.edu/mirSTP/.

Linkage disequilibrium-dependent architecture of human complex traits shows action of negative selection

Recent work has hinted at the linkage disequilibrium (LD)-dependent architecture of human complex traits, where SNPs with low levels of LD (LLD) have larger per-SNP heritability. Here we analyzed summary statistics from 56 complex traits (average N = 101,401) by extending stratified LD score regression to continuous annotations. We determined that SNPs with low LLD have significantly larger per-SNP heritability and that roughly half of this effect can be explained by functional annotations negatively correlated with LLD, such as DNase I hypersensitivity sites (DHSs). The remaining signal is largely driven by our finding that more recent common variants tend to have lower LLD and to explain more heritability (P = 2.38 × 10^-104); the youngest 20% of common SNPs explain 3.9 times more heritability than the oldest 20%, consistent with the action of negative selection. We also inferred jointly significant effects of other LD-related annotations and confirmed via forward simulations that they jointly predict deleterious effects.

TMPRSS2-ERG fusion co-opts master transcription factors and activates NOTCH signaling in primary prostate cancer

TMPRSS2-ERG (T2E) structural rearrangements typify ∼50% of prostate tumors and result in overexpression of the ERG transcription factor. Using chromatin, genomic and expression data, we show distinct cis-regulatory landscapes between T2E-positive and non-T2E primary prostate tumors, which include clusters of regulatory elements (COREs). This difference is mediated by ERG co-option of HOXB13 and FOXA1, implementing a T2E-specific transcriptional profile. We also report a T2E-specific CORE on the structurally rearranged ERG locus arising from spreading of the TMPRSS2 locus pre-existing CORE, assisting in its overexpression. Finally, we show that the T2E-specific cis-regulatory landscape underlies a vulnerability against the NOTCH pathway. Indeed, NOTCH pathway inhibition antagonizes the growth and invasion of T2E-positive prostate cancer cells. Taken together, our work shows that overexpressed ERG co-opts master transcription factors to deploy a unique cis-regulatory landscape, inducing a druggable dependency on NOTCH signaling in T2E-positive prostate tumors.

Analysis of Histone Antibody Specificity with Peptide Microarrays

Post-translational modifications (PTMs) on histone proteins are widely studied for their roles in regulating chromatin structure and gene expression. The mass production and distribution of antibodies specific to histone PTMs has greatly facilitated research on these marks. As histone PTM antibodies are key reagents for many chromatin biochemistry applications, rigorous analysis of antibody specificity is necessary for accurate data interpretation and continued progress in the field. This protocol describes an integrated pipeline for the design, fabrication and use of peptide microarrays for profiling the specificity of histone antibodies. The design and analysis aspects of this procedure are facilitated by ArrayNinja, an open-source and interactive software package we recently developed to streamline the customization of microarray print formats. This pipeline has been used to screen a large number of commercially available and widely used histone PTM antibodies, and data generated from these experiments are freely available through an online and expanding Histone Antibody Specificity Database. Beyond histones, the general methodology described herein can be applied broadly to the analysis of PTM-specific antibodies.

p39 Is Responsible for Increasing Cdk5 Activity during Postnatal Neuron Differentiation and Governs Neuronal Network Formation and Epileptic Responses

Two distinct protein cofactors, p35 and p39, independently activate Cyclin-dependent kinase 5 (Cdk5), which plays diverse roles in normal brain function and the pathogenesis of many neurological diseases. The initial discovery that loss of p35 impairs neuronal migration in the embryonic brain prompted intensive research exploring the function of p35-dependent Cdk5 activity. In contrast, p39 expression is restricted to the postnatal brain and its function remains poorly understood. Despite the robustly increased Cdk5 activity during neuronal differentiation, which activator is responsible for enhancing Cdk5 activation and how the two distinct activators direct Cdk5 signaling to govern neuronal network formation and function still remains elusive. Here we report that p39, but not p35, is selectively upregulated by histone acetylation-mediated transcription, which underlies the robust increase of Cdk5 activity during rat and mouse neuronal differentiation. The loss of p39 attenuates overall Cdk5 activity in neurons and preferentially affects phosphorylation of specific Cdk5 targets, leading to aberrant axonal growth and impaired dendritic spine and synapse formation. In adult mouse brains, p39 deficiency results in dysregulation of p35 and Cdk5 targets in synapses. Moreover, in contrast to the proepileptic phenotype caused by the lack of p35, p39 loss leads to deficits in maintaining seizure activity and induction of immediate early genes that control hippocampal excitability. Together, our studies demonstrate essential roles of p39 in neuronal network development and function. Furthermore, our data support a model in which Cdk5 activators play nonoverlapping and even opposing roles to govern balanced Cdk5 signaling in the postnatal brain.

SIGNIFICANCE STATEMENT

Neuronal network development requires tightly regulated activation of Cyclin-dependent kinase 5 (Cdk5) by two distinct cofactors, p35 and p39. Despite the well-known p35-dependent Cdk5 function, why postnatal neurons express abundant p39 in addition to p35 remained unknown for decades. In this study, we discovered that selective upregulation of p39 is the underlying mechanism that accommodates the increased functional requirement of Cdk5 activation during neuronal differentiation. In addition, we demonstrated that p39 selectively directs Cdk5 to phosphorylate protein substrates essential for axonal development, dendritic spine formation, and synaptogenesis. Moreover, our studies suggest opposing roles of p39 and p35 in synaptic Cdk5 function and epileptic responses, arguing that cooperation between Cdk5 activators maintains balanced Cdk5 signing, which is crucial for postnatal brain function.

Sequence based predictor for discrimination of enhancer and their types by applying general form of Chou's trinucleotide composition

BACKGROUND AND OBJECTIVES

Enhancers are pivotal DNA elements, which are widely used in eukaryotes for activation of transcription genes. On the basis of enhancer strength, they are further classified into two groups; strong enhancers and weak enhancers. Due to high availability of huge amount of DNA sequences, it is needed to develop fast, reliable and robust intelligent computational method, which not only identify enhancers but also determines their strength. Considerable progress has been achieved in this regard; however, timely and precisely identification of enhancers is still a challenging task.

METHODS

Two-level intelligent computational model for identification of enhancers and their subgroups is proposed. Two different feature extraction techniques including di-nucleotide composition and tri-nucleotide composition were adopted for extraction of numerical descriptors. Four classification methods including probabilistic neural network, support vector machine, k-nearest neighbor and random forest were utilized for classification.

RESULTS

The proposed method yielded 77.25% of accuracy for dataset S₁ contains enhancers and non-enhancers, whereas 64.70% of accuracy for dataset S₂ comprises of strong enhancer and weak enhancer sequences using jackknife cross-validation test.

CONCLUSION

The predictive results validated that the proposed method is better than that of existing approaches so far reported in the literature. It is thus highly observed that the developed method will be useful and expedient for basic research and academia.

Epigenetic Modifications to H3K9 in Renal Tubulointerstitial Cells after Unilateral Ureteric Obstruction and TGF-β1 Stimulation

Introduction: Epigenetic regulation of fibrogenesis through post-translational histone modifications (marks) may be a key determinant of progression in renal disease. In this study, we examined the distribution and acquisition of histone 3 Lysine 9 (H3K9) marks after injury and stimulation with the pro-fibrotic cytokine TGF-β1. Our focus was on their presence in activated fibroblasts (myofibroblasts) and epithelial cells (epithelial-mesenchymal transition).

Methods and Results: Immunofluorescent microscopy was used to examine global H3K9 acetylation (H3K9Ac) and tri-methylation (H3K9Me3) after unilateral ureteric obstruction (UUO) in mice. Confocal, super resolution microscopy and flow cytometry were used to determine the in vitro effect of TGF-β1 on structural arrangement of these marks, and their relationship with α-smooth muscle actin (αSMA) expression, a marker of myofibroblasts and early EMT. The number of individual histone marks was increased 10 days after UUO (p < 0.05 vs. control), with both marks clearly seen in various cell types including proximal tubules and myofibroblasts. Sub-nuclear microscopy in primary rat renal fibroblasts and a proximal tubule cell line (NRK-52e) showed that H3K9Ac was co-localized with phosphorylated-Ser2 RNA polymerase II (pRNAPol II), while H3K9Me3 was not, consistent with permissive and repressive effects on gene expression respectively. In both cell types H3K9Ac was diffusely distributed throughout the nucleus, while H3K9Me3 was found in compartments resembling the nucleolus, and in the case of the fibroblast, also juxtapositioned with the nuclear membrane. TGF-β1 had no effect on H3K9Ac marks in either cell, but resulted in a redistribution of H3K9Me3 within the fibroblast nucleus. This was unrelated to any change in mitogenesis, but was associated with increased αSMA expression.

Conclusion: These findings highlight why it is important to consider the epigenetics of each cell individually, because whilst no overall enrichment occurred, renal myofibroblast differentiation was accompanied by distinct changes in histone mark arrangements.

DNA methylation at the mu-1 opioid receptor gene (OPRM1) promoter predicts preoperative, acute, and chronic postsurgical pain after spine fusion

Introduction

The perioperative pain experience shows great interindividual variability and is difficult to predict. The mu-1 opioid receptor gene (OPRM1) is known to play an important role in opioid-pain pathways. Since deoxyribonucleic acid (DNA) methylation is a potent repressor of gene expression, DNA methylation was evaluated at the OPRM1 promoter, as a predictor of preoperative, acute, and chronic postsurgical pain (CPSP).

Methods

A prospective observational cohort study was conducted in 133 adolescents with idiopathic scoliosis undergoing spine fusion under standard protocols. Data regarding pain, opioid consumption, anxiety, and catastrophizing (using validated questionnaires) were collected before and 2–3 months postsurgery. Outcomes evaluated were preoperative pain, acute postoperative pain (area under curve [AUC] for pain scores over 48 hours), and CPSP (numerical rating scale >3/10 at 2–3 months postsurgery). Blood samples collected preoperatively were analyzed for DNA methylation by pyrosequencing of 22 CpG sites at the OPRM1 gene promoter. The association of each pain outcome with the methylation percentage of each CpG site was assessed using multivariable regression, adjusting for significant (P<0.05) nongenetic variables.

Results

Majority (83%) of the patients reported no pain preoperatively, while CPSP occurred in 36% of the subjects (44/121). Regression on dichotomized preoperative pain outcome showed association with methylation at six CpG sites (1, 3, 4, 9, 11, and 17) (P<0.05). Methylation at CpG sites 4, 17, and 18 was associated with higher AUC after adjusting for opioid consumption and preoperative pain score (P<0.05). After adjusting for postoperative opioid consumption and preoperative pain score, methylation at CpG sites 13 and 22 was associated with CPSP (P<0.05).

Discussion

Novel CPSP biomarkers were identified in an active regulatory region of the OPRM1 gene that binds multiple transcription factors. Inhibition of binding by DNA methylation potentially decreases the OPRM1 gene expression, leading to a decreased response to endogenous and exogenous opioids, and an increased pain experience.

IDH1 or -2 mutations do not predict outcome and do not cause loss of 5-hydroxymethylcytosine or altered histone modifications in central chondrosarcomas

BACKGROUND

Mutations in isocitrate dehydrogenase (IDH)1 or -2 are found in ~50% of conventional central chondrosarcomas and in up to 87% of their assumed benign precursors enchondromas. The mutant enzyme acquires the activity to convert α-ketoglutarate into the oncometabolite d-2-hydroxyglutarate (d-2-HG), which competitively inhibits α-ketoglutarate dependent enzymes such as histone- and DNA demethylases.

METHODS

We therefore evaluated the effect of IDH1 or -2 mutations on histone modifications (H3K4me3, H3K9me3 and H3K27me3), chromatin remodeler ATRX expression, DNA modifications (5-hmC and 5-mC), and TET1 subcellular localization in a genotyped cohort (IDH, succinate dehydrogenase (SDH) and fumarate hydratase (FH)) of enchondromas and central chondrosarcomas (n = 101) using immunohistochemistry.

RESULTS

IDH1 or -2 mutations were found in 60.8% of the central cartilaginous tumours, while mutations in FH and SDH were absent. The mutation status did not correlate with outcome. Chondrosarcomas are strongly positive for the histone modifications H3K4me3, H3K9me3 and H3K27me3, which was independent of the IDH1 or -2 mutation status. Two out of 36 chondrosarcomas (5.6%) show complete loss of ATRX. Levels of 5-hmC and 5-mC are highly variable in central cartilaginous tumours and are not associated with mutation status. In tumours with loss of 5-hmC, expression of TET1 was more prominent in the cytoplasm than the nucleus (p = 0.0001).

CONCLUSIONS

In summary, in central chondrosarcoma IDH1 or -2 mutations do not affect immunohistochemical levels of 5-hmC, 5mC, trimethylation of H3K4, -K9 and K27 and outcome, as compared to wildtype.

Effect of irinotecan on HMGB1, MMP9 expression, cell cycle, and cell growth in breast cancer (MCF-7) cells

Irinotecan is a natural alkaloid agent widely used in cancer therapy. High-mobility group protein B1 as a non-histone chromosomal protein plays a fundamental role in gene expression and inflammation. In this study, the effect of irinotecan on high-mobility group protein B1 and MMP9 content, gene expression, cell cycle, and cell growth in human breast cancer cells (MCF-7) was investigated. The cells were exposed to various concentrations of irinotecan and the viability determined by trypan blue exclusion and 3-(4,5-dimethylthiazal-2-yl)-2,5-diphenyltetrazolium bromide assays. High-mobility group B proteins were extracted from the control and drug-treated cells and analyzed by immunoblot. High-mobility group protein B1 and MMP9 messenger RNA expression was studied by reverse transcription polymerase chain reaction. The results demonstrated reduction of cell viability upon increasing irinotecan concentration, up-regulated high-mobility group protein B1 gene expression, and down-regulated MMP9 mRNA. Although the content of high-mobility group protein B1 was decreased in chromatin extract upon drug action, no high-mobility group protein B1 release to extracellular space was detected by immunoblot analysis. Irinotecan decreased H3K9 acetylation and increased poly ADP-ribose polymerase fragmentation to 89 kDa and anion superoxide production suggesting induction of apoptosis in these cells. Propidium iodide staining of the cells 24 h after the drug treatment revealed arrest of the cells in S-phase. From the results, it is concluded that overexpression of high-mobility group protein B1 in the presence of irinotecan precedes breast cancer cells into apoptosis and in this response the binding of irinotecan to chromatin or high-mobility group protein B1 may condense/aggregate chromatin, preventing high-mobility group protein B1 release from chromatin.

Characterization of Complete Histone Tail Proteoforms Using Differential Ion Mobility Spectrometry

Histone proteins are subject to dynamic post-translational modifications (PTMs) that cooperatively modulate the chromatin structure and function. Nearly all functional PTMs are found on the N-terminal histone domains (tails) of ∼50 residues protruding from the nucleosome core. Using high-definition differential ion mobility spectrometry (FAIMS) with electron transfer dissociation, we demonstrate rapid baseline gas-phase separation and identification of tails involving monomethylation, trimethylation, acetylation, or phosphorylation in biologically relevant positions. These are by far the largest variant peptides resolved by any method, some with PTM contributing just 0.25% to the mass. This opens the door to similar separations for intact proteins and in top-down proteomics.

l-Arginine-Dependent Epigenetic Regulation of Interleukin-10, but Not Transforming Growth Factor-β, Production by Neonatal Regulatory T Lymphocytes

A growing number of diseases in humans, including trauma, certain cancers, and infection, are known to be associated with l-arginine deficiency. In addition, l-arginine must be supplemented by diet during pregnancy to aid fetal development. In conditions of l-arginine depletion, T cell proliferation is impaired. We have previously shown that neonatal blood has lower l-arginine levels than adult blood, which is associated with poor neonatal lymphocyte proliferation, and that l-arginine enhances neonatal lymphocyte proliferation through an interleukin (IL)-2-independent pathway. In this study, we have further investigated how exogenous l-arginine enhances neonatal regulatory T-cells (Tregs) function in relation to IL-10 production under epigenetic regulation. Results showed that cord blood mononuclear cells (CBMCs) produced higher levels of IL-10 than adult peripheral blood mononuclear cells (PBMCs) by phytohemagglutinin stimulation but not by anti-CD3/anti-CD28 stimulation. Addition of exogenous l-arginine had no effect on transforming growth factor-β production by PBMCs or CBMCs, but enhanced IL-10 production by neonatal CD4⁺CD25⁺FoxP3⁺ Tregs. Further studies showed that IL-10 promoter DNA hypomethylation, rather than histone modification, corresponded to the l-arginine-induced increase in IL-10 production by neonatal CD4⁺ T cells. These results suggest that l-arginine modulates neonatal Tregs through the regulation of IL-10 promoter DNA methylation. l-arginine supplementation may correct the Treg function in newborns with l-arginine deficiency.

Deregulation of α-synuclein in Parkinson's disease: Insight from epigenetic structure and transcriptional regulation of SNCA

Understanding regulation of α-synuclein has long been a central focus for Parkinson's disease (PD) researchers. Accumulation of this protein in the Lewy body or neurites, mutations in the coding region of the gene and strong association of α-synuclein encoding gene multiplication (duplication/triplication) with familial form of PD have indicated the importance of this molecule in pathogenesis of the disease. Several years of research identified many potential faulty pathways associated with accumulation of α-synuclein inside dopaminergic neurons and its transmission to neighboring ones. Concurrently, an appreciable body of research is growing to understand the epigenetic and genetic deregulation of α-synuclein that might contribute to the disease pathology. Completion of the ENCODE (Encyclopedia of DNA Elements) project and recent advancement made in the epigenetic and trans factor mediated regulation of each gene, has tremendously accelerated the need to carefully understand the epigenetic structure of the gene (SNCA) encoding α-synuclein protein in order to decipher the regulation and contribution of α-synuclein to the pathogenesis of PD. We have also analyzed the detailed epigenetic structure of this gene with knowledge from ENCODE database, which may open new avenues in α-synuclein research. Interestingly, we have found that the gene contains several transcriptionally activate histone modifications and associated potential transcription factor binding sites in the non-coding areas that strongly suggest alternative regulatory pathways. Altogether this review will provide interesting insight of α-synuclein gene regulation from epigenetic, genetic and post-transcriptional perspectives and their potential implication in the PD pathogenesis.

How does chromatin package DNA within nucleus and regulate gene expression?

The human body is made up of 60 trillion cells, each cell containing 2 millions of genomic DNA in its nucleus. How is this genomic deoxyribonucleic acid [DNA] organised into nuclei? Around 1880, W. Flemming discovered a nuclear substance that was clearly visible on staining under primitive light microscopes and named it 'chromatin'; this is now thought to be the basic unit of genomic DNA organization. Since long before DNA was known to carry genetic information, chromatin has fascinated biologists. DNA has a negatively charged phosphate backbone that produces electrostatic repulsion between adjacent DNA regions, making it difficult for DNA to fold upon itself. In this article, we will try to shed light on how does chromatin package DNA within nucleus and regulate gene expression?

ERK-Induced Activation of TCF Family of SRF Cofactors Initiates a Chromatin Modification Cascade Associated with Transcription

We investigated the relationship among ERK signaling, histone modifications, and transcription factor activity, focusing on the ERK-regulated ternary complex factor family of SRF partner proteins. In MEFs, activation of ERK by TPA stimulation induced a common pattern of H3K9acS10ph, H4K16ac, H3K27ac, H3K9acK14ac, and H3K4me3 at hundreds of transcription start site (TSS) regions and remote regulatory sites. The magnitude of the increase in histone modification correlated well with changes in transcription. H3K9acS10ph preceded the other modifications. Most induced changes were TCF dependent, but TCF-independent TSSs exhibited the same hierarchy, indicating that it reflects gene activation per se. Studies with TCF Elk-1 mutants showed that TCF-dependent ERK-induced histone modifications required Elk-1 to be phosphorylated and competent to activate transcription. Analysis of direct TCF-SRF target genes and chromatin modifiers confirmed this and showed that H3S10ph required only Elk-1 phosphorylation. Induction of histone modifications following ERK stimulation is thus directed by transcription factor activation and transcription.

Aorta macrophage inflammatory and epigenetic changes in a murine model of obstructive sleep apnea: Potential role of CD36

Obstructive sleep apnea (OSA) affects 8-10% of the population, is characterized by chronic intermittent hypoxia (CIH), and causally associates with cardiovascular morbidities. In CIH-exposed mice, closely mimicking the chronicity of human OSA, increased accumulation and proliferation of pro-inflammatory metabolic M1-like macrophages highly expressing CD36, emerged in aorta. Transcriptomic and MeDIP-seq approaches identified activation of pro-atherogenic pathways involving a complex interplay of histone modifications in functionally-relevant biological pathways, such as inflammation and oxidative stress in aorta macrophages. Discontinuation of CIH did not elicit significant improvements in aorta wall macrophage phenotype. However, CIH-induced aorta changes were absent in CD36 knockout mice, Our results provide mechanistic insights showing that CIH exposures during sleep in absence of concurrent pro-atherogenic settings (i.e., genetic propensity or dietary manipulation) lead to the recruitment of CD36(+)^high macrophages to the aortic wall and trigger atherogenesis. Furthermore, long-term CIH-induced changes may not be reversible with usual OSA treatment.

ChIPBase v2.0: decoding transcriptional regulatory networks of non-coding RNAs and protein-coding genes from ChIP-seq data

The abnormal transcriptional regulation of non-coding RNAs (ncRNAs) and protein-coding genes (PCGs) is contributed to various biological processes and linked with human diseases, but the underlying mechanisms remain elusive. In this study, we developed ChIPBase v2.0 (http://rna.sysu.edu.cn/chipbase/) to explore the transcriptional regulatory networks of ncRNAs and PCGs. ChIPBase v2.0 has been expanded with ∼10 200 curated ChIP-seq datasets, which represent about 20 times expansion when comparing to the previous released version. We identified thousands of binding motif matrices and their binding sites from ChIP-seq data of DNA-binding proteins and predicted millions of transcriptional regulatory relationships between transcription factors (TFs) and genes. We constructed 'Regulator' module to predict hundreds of TFs and histone modifications that were involved in or affected transcription of ncRNAs and PCGs. Moreover, we built a web-based tool, Co-Expression, to explore the co-expression patterns between DNA-binding proteins and various types of genes by integrating the gene expression profiles of ∼10 000 tumor samples and ∼9100 normal tissues and cell lines. ChIPBase also provides a ChIP-Function tool and a genome browser to predict functions of diverse genes and visualize various ChIP-seq data. This study will greatly expand our understanding of the transcriptional regulations of ncRNAs and PCGs.