Chromatin signatures of active enhancers

Histone H3K18 & H3K23 acetylation directs establishment of MLL-mediated H3K4 methylation

In an unmodified state, positively charged histone N-terminal tails engage nucleosomal DNA in a manner which restricts access to not only the underlying DNA but also key tail residues subject to binding and/or modification. Charge-neutralizing modifications, such as histone acetylation, serve to disrupt this DNA-tail interaction, facilitating access to such residues. We previously showed that a polyacetylation-mediated chromatin "switch" governs the read-write capability of H3K4me3 by the MLL1 methyltransferase complex. Here, we discern the relative contributions of site-specific acetylation states along the H3 tail and extend our interrogation to other chromatin modifiers. We show that the contributions of H3 tail acetylation to H3K4 methylation by MLL1 are highly variable, with H3K18 and H3K23 acetylation exhibiting robust stimulatory effects and that this extends to the related H3K4 methyltransferase complex, MLL4. We show that H3K4me1 and H3K4me3 are found preferentially co-enriched with H3 N-terminal tail proteoforms bearing dual H3K18 and H3K23 acetylation (H3{K18acK23ac}). We further show that this effect is specific to H3K4 methylation, while methyltransferases targeting other H3 tail residues (H3K9, H3K27, & H3K36), a methyltransferase targeting the nucleosome core (H3K79), and a kinase targeting a residue directly adjacent to H3K4 (H3T3) are insensitive to tail acetylation. Together, these findings indicate a unique and robust stimulation of H3K4 methylation by H3K18 and H3K23 acetylation and provide key insight into why H3K4 methylation is often associated with histone acetylation in the context of active gene expression.

Non-Coding RNA-Mediated Gene Regulation in Cardiovascular Disorders: Current Insights and Future Directions

Cardiovascular diseases (CVDs) pose a significant burden on global health. Developing effective diagnostic, therapeutic, and prognostic indicators for CVDs is critical. This narrative review explores the role of select non-coding RNAs (ncRNAs) and provides an in-depth exploration of the roles of miRNAs, lncRNAs, and circRNAs in different aspects of CVDs, offering insights into their mechanisms and potential clinical implications. The review also sheds light on the diverse functions of ncRNAs, including their modulation of gene expression, epigenetic modifications, and signaling pathways. It comprehensively analyzes the interplay between ncRNAs and cardiovascular health, paving the way for potential novel interventions. Finally, the review provides insights into the methodologies used to investigate ncRNA-mediated gene regulation in CVDs, as well as the implications and challenges associated with translating ncRNA research into clinical applications. Considering the broader implications, this research opens avenues for interdisciplinary collaborations, enhancing our understanding of CVDs across scientific disciplines.

CRISPR-Cas9-mediated deletion enhancer of MECOM play a tumor suppressor role in ovarian cancer

MDS1 and EVI1 complex locus (MECOM), a transcription factor encoding several variants, has been implicated in progression of ovarian cancer. The function of regulatory regions in regulating MECOM expression in ovarian cancer is not fully understood. In this study, MECOM expression was evaluated in ovarian cancer cell lines treated with bromodomain and extraterminal (BET) inhibitor JQ-1. Oncogenic phenotypes were assayed using assays of CCK-8, colony formation, wound-healing and transwell. Oncogenic phenotypes were estimated in stable sgRNA-transfected OVCAR3 cell lines. Xenograft mouse model was assayed via subcutaneous injection of enhancer-deleted OVCAR3 cell lines. The results displayed that expression of MECOM is downregulated in cell lines treated with JQ-1. Data from published ChIP-sequencing (H3K27Ac) in 3 ovarian cancer cell lines displayed a potential enhancer around the first exon. mRNA and protein expression were downregulated in OVCAR3 cells after deletion of the MECOM enhancer. Similarly, oncogenic phenotypes both in cells and in the xenograft mouse model were significantly attenuated. This study demonstrates that JQ-1 can inhibit the expression of MECOM and tumorigenesis. Deletion of the enhancer activity of MECOM has an indispensable role in inhibiting ovarian cancer progress, which sheds light on a promising opportunity for ovarian cancer treatment through the application of this non-coding DNA deletion.

Role of the histone variant H2A.Z.1 in memory, transcription, and alternative splicing is mediated by lysine modification

Creating long-lasting memories requires learning-induced changes in gene expression, which are impacted by epigenetic modifications of DNA and associated histone proteins. Post-translational modifications (PTMs) of histones are key regulators of transcription, with different PTMs producing unique effects on gene activity and behavior. Although recent studies implicate histone variants as novel regulators of memory, effects of PTMs on the function of histone variants are rarely considered. We previously showed that the histone variant H2A.Z suppresses memory, but it is unclear if this role is impacted by H2A.Z acetylation, a PTM that is typically associated with positive effects on transcription and memory. To answer this question, we used a mutation approach to manipulate acetylation on H2A.Z without impacting acetylation of other histone types. Specifically, we used adeno-associated virus (AAV) constructs to overexpress mutated H2A.Z.1 isoforms that either mimic acetylation (acetyl-mimic) by replacing lysines 4, 7 and 11 with glutamine (KQ), or H2A.Z.1 with impaired acetylation (acetyl-defective) by replacing the same lysines with alanine (KA). Expressing the H2A.Z.1 acetyl-mimic (H2A.Z.1KQ) improved memory under weak learning conditions, whereas expressing the acetyl-defective H2A.Z.1KA generally impaired memory, indicating that the effect of H2A.Z.1 on memory depends on its acetylation status. RNA sequencing showed that H2A.Z.1KQ and H2A.Z.1KA uniquely impact the expression of different classes of genes in both females and males. Specifically, H2A.Z.1KA preferentially impacts genes involved in synaptic function, suggesting that acetyl-defective H2A.Z.1 impairs memory by altering synaptic regulation. Finally, we describe, for the first time, that H2A.Z is also involved in alternative splicing of neuronal genes, whereby H2A.Z depletion, as well as expression of H2A.Z.1 lysine mutants influence transcription and splicing of different gene targets, suggesting that H2A.Z.1 can impact behavior through effects on both splicing and gene expression. This is the first study to demonstrate that direct manipulation of H2A.Z post-translational modifications regulates memory, whereby acetylation adds another regulatory layer by which histone variants can fine tune higher brain functions through effects on gene expression and splicing.

H3K18 & H3K23 acetylation directs establishment of MLL-mediated H3K4 methylation

In an unmodified state, positively charged histone N-terminal tails engage nucleosomal DNA in a manner which restricts access to not only the underlying DNA, but also key tail residues subject to binding and/or modification. Charge-neutralizing modifications, such as histone acetylation, serve to disrupt this DNA-tail interaction, facilitating access to such residues. We previously showed that a polyacetylation-mediated chromatin "switch" governs the read-write capability of H3K4me3 by the MLL1 methyltransferase complex. Here, we discern the relative contributions of site-specific acetylation states along the H3 tail and extend our interrogation to other chromatin modifiers. We show that the contributions of H3 tail acetylation to H3K4 methylation by MLL1 are highly variable, with H3K18 and H3K23 acetylation exhibiting robust stimulatory effects, and that this extends to the related H3K4 methyltransferase complex, MLL4. We show that H3K4me1 and H3K4me3 are found preferentially co-enriched with H3 N-terminal tail proteoforms bearing dual H3K18 and H3K23 acetylation (H3{K18acK23ac}). We further show that this effect is specific to H3K4 methylation, while methyltransferases targeting other H3 tail residues (H3K9, H3K27, & H3K36), a methyltransferase targeting the nucleosome core (H3K79), and a kinase targeting a residue directly adjacent to H3K4 (H3T3) are insensitive to tail acetylation. Together, these findings indicate a unique and robust stimulation of H3K4 methylation by H3K18 and H3K23 acetylation and provide key insight into why H3K4 methylation is often associated with histone acetylation in the context of active gene expression.

Prognostic analysis of mutated genes and insight into effects of DNA damage and repair on mutational strand asymmetries in gastric cancer

Gastric cancer (GACA) is a complex and multifaceted disease influenced by a variety of environmental and genetic factors. Somatic mutations play a major role in its development, and their characteristics, including the asymmetry between two DNA strands, are of great interest and appear as a signal of information and guidance, revealing mechanisms of DNA damage and repair. Here, we analyzed the impact of High-frequency mutated genes on patient prognosis and found that the effect of expression levels of tumor protein p53 (TP53) and lysine methyltransferase 2C (KMT2C) genes remained high throughout the development of GACA, with similar expression patterns. After investigating mutation asymmetry across mutagenic processes, we found that transcriptional asymmetry was dominated by T > G mutations under the influence of transcription couples repair and damage. The apolipoprotein B mRNA editing enzyme catalytic polypeptide like (APOBEC) enzyme that induces mutations during DNA replication has been identified here and we identified a replicative asymmetry, which was dominated by C > A mutations in left-replicating. Strand bias in different mutation classes at transcription factor binding sites and enhancer regions were also confirmed, which implies the important role of non-coding regulatory elements in the occurrence of mutations. This work systematically describes mutational strand asymmetries in specific genomic regions, shedding light on the DNA damage and repair mechanisms underlying somatic mutations in cohorts of GACA patients with gastric cancer.

Systematic identification of genotype-dependent enhancer variants in eosinophilic esophagitis

Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation, that develops in a small subset of people with food allergies. Genome-wide association studies (GWASs) have identified 9 independent EoE risk loci reaching genome-wide significance (p < 5 × 10-8) and 27 additional loci of suggestive significance (5 × 10-8 < p < 1 × 10-5). In the current study, we perform linkage disequilibrium (LD) expansion of these loci to nominate a set of 531 variants that are potentially causal. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their genomic sequence context cloned into a GFP reporter library. Analysis of reporter gene expression in TE-7, HaCaT, and Jurkat cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants, representing 6 genome-wide significant EoE loci and 7 suggestive EoE loci, that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in individuals with EoE. Transcription factor enrichment analyses reveal possible roles for cell-type-specific regulators, including GATA3. Our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing functional insights into the effects of genetic variation in this disease.

ChIP-DIP: A multiplexed method for mapping hundreds of proteins to DNA uncovers diverse regulatory elements controlling gene expression

Gene expression is controlled by the dynamic localization of thousands of distinct regulatory proteins to precise regions of DNA. Understanding this cell-type specific process has been a goal of molecular biology for decades yet remains challenging because most current DNA-protein mapping methods study one protein at a time. To overcome this, we developed ChIP-DIP (ChIP Done In Parallel), a split-pool based method that enables simultaneous, genome-wide mapping of hundreds of diverse regulatory proteins in a single experiment. We demonstrate that ChIP-DIP generates highly accurate maps for all classes of DNA-associated proteins, including histone modifications, chromatin regulators, transcription factors, and RNA Polymerases. Using these data, we explore quantitative combinations of protein localization on genomic DNA to define distinct classes of regulatory elements and their functional activity. Our data demonstrate that ChIP-DIP enables the generation of 'consortium level', context-specific protein localization maps within any molecular biology lab.

Cellular reprogramming is driven by widespread rewiring of promoter-enhancer interactions

BACKGROUND

Long-range interactions between promoters and cis-regulatory elements, such as enhancers, play critical roles in gene regulation. However, the role of three-dimensional (3D) chromatin structure in orchestrating changes in transcriptional regulation during direct cell reprogramming is not fully understood.

RESULTS

Here, we performed integrated analyses of chromosomal architecture, epigenetics, and gene expression using Hi-C, promoter Capture Hi-C (PCHi-C), ChIP-seq, and RNA-seq during trans-differentiation of Pre-B cells into macrophages with a β-estradiol inducible C/EBPαER transgene. Within 1h of β-estradiol induction, C/EBPα translocated from the cytoplasm to the nucleus, binding to thousands of promoters and putative regulatory elements, resulting in the downregulation of Pre-B cell-specific genes and induction of macrophage-specific genes. Hi-C results were remarkably consistent throughout trans-differentiation, revealing only a small number of TAD boundary location changes, and A/B compartment switches despite significant changes in the expression of thousands of genes. PCHi-C revealed widespread changes in promoter-anchored loops with decreased interactions in parallel with decreased gene expression, and new and increased promoter-anchored interactions in parallel with increased expression of macrophage-specific genes.

CONCLUSIONS

Overall, our data demonstrate that C/EBPα-induced trans-differentiation involves few changes in genome architecture at the level of TADs and A/B compartments, in contrast with widespread reorganization of thousands of promoter-anchored loops in association with changes in gene expression and cell identity.

Characterization of transcriptional enhancers in the chicken genome using CRISPR-mediated activation

DNA regulatory elements intricately control when, where, and how genes are activated. Therefore, understanding the function of these elements could unveil the complexity of the genetic regulation network. Genome-wide significant variants are predominantly found in non-coding regions of DNA, so comprehending the predicted functional regulatory elements is crucial for understanding the biological context of these genomic markers, which can be incorporated into breeding programs. The emergence of CRISPR technology has provided a powerful tool for studying non-coding regulatory elements in genomes. In this study, we leveraged epigenetic data from the Functional Annotation of Animal Genomes project to identify promoter and putative enhancer regions associated with three genes (HBBA, IRF7, and PPARG) in the chicken genome. To identify the enhancer regions, we designed guide RNAs targeting the promoter and candidate enhancer regions and utilized CRISPR activation (CRISPRa) with dCas9-p300 and dCas9-VPR as transcriptional activators in chicken DF-1 cells. By comparing the expression levels of target genes between the promoter activation and the co-activation of the promoter and putative enhancers, we were able to identify functional enhancers that exhibited augmented upregulation. In conclusion, our findings demonstrate the remarkable efficiency of CRISPRa in precisely manipulating the expression of endogenous genes by targeting regulatory elements in the chicken genome, highlighting its potential for functional validation of non-coding regions.

Rho enhancers play unexpectedly minor roles in Rhodopsin transcription and rod cell integrity

Enhancers function with a basal promoter to control the transcription of target genes. Enhancer regulatory activity is often studied using reporter-based transgene assays. However, unmatched results have been reported when selected enhancers are silenced in situ. In this study, using genomic deletion analysis in mice, we investigated the roles of two previously identified enhancers and the promoter of the Rho gene that codes for the visual pigment rhodopsin. The Rho gene is robustly expressed by rod photoreceptors of the retina, and essential for the subcellular structure and visual function of rod photoreceptors. Mutations in RHO cause severe vision loss in humans. We found that each Rho regulatory region can independently mediate local epigenomic changes, but only the promoter is absolutely required for establishing active Rho chromatin configuration and transcription and maintaining the cell integrity and function of rod photoreceptors. To our surprise, two Rho enhancers that enable strong promoter activation in reporter assays are largely dispensable for Rho expression in vivo. Only small and age-dependent impact is detectable when both enhancers are deleted. Our results demonstrate context-dependent roles of enhancers and highlight the importance of studying functions of cis-regulatory regions in the native genomic context.

Ovol1/2 loss-induced epidermal defects elicit skin immune activation and alter global metabolism

Skin epidermis constitutes the outer permeability barrier that protects the body from dehydration, heat loss, and myriad external assaults. Mechanisms that maintain barrier integrity in constantly challenged adult skin and how epidermal dysregulation shapes the local immune microenvironment and whole-body metabolism remain poorly understood. Here, we demonstrate that inducible and simultaneous ablation of transcription factor-encoding Ovol1 and Ovol2 in adult epidermis results in barrier dysregulation through impacting epithelial-mesenchymal plasticity and inflammatory gene expression. We find that aberrant skin immune activation then ensues, featuring Langerhans cell mobilization and T cell responses, and leading to elevated levels of secreted inflammatory factors in circulation. Finally, we identify failure to gain body weight and accumulate body fat as long-term consequences of epidermal-specific Ovol1/2 loss and show that these global metabolic changes along with the skin barrier/immune defects are partially rescued by immunosuppressant dexamethasone. Collectively, our study reveals key regulators of adult barrier maintenance and suggests a causal connection between epidermal dysregulation and whole-body metabolism that is in part mediated through aberrant immune activation.

A comprehensive revisit of the machine-learning tools developed for the identification of enhancers in the human genome

Enhancers are non-coding DNA elements that play a crucial role in enhancing the transcription rate of a specific gene in the genome. Experiments for identifying enhancers can be restricted by their conditions and involve complicated, time-consuming, laborious, and costly steps. To overcome these challenges, computational platforms have been developed to complement experimental methods that enable high-throughput identification of enhancers. Over the last few years, the development of various enhancer computational tools has resulted in significant progress in predicting putative enhancers. Thus, researchers are now able to use a variety of strategies to enhance and advance enhancer study. In this review, an overview of machine learning (ML)-based prediction methods for enhancer identification and related databases has been provided. The existing enhancer-prediction methods have also been reviewed regarding their algorithms, feature selection processes, validation techniques, and software utility. In addition, the advantages and drawbacks of these ML approaches and guidelines for developing bioinformatic tools have been highlighted for a more efficient enhancer prediction. This review will serve as a useful resource for experimentalists in selecting the appropriate ML tool for their study, and for bioinformaticians in developing more accurate and advanced ML-based predictors.

Chromatin accessibility differences between alpha, beta, and delta cells identifies common and cell type-specific enhancers

BACKGROUND

High throughput sequencing has enabled the interrogation of the transcriptomic landscape of glucagon-secreting alpha cells, insulin-secreting beta cells, and somatostatin-secreting delta cells. These approaches have furthered our understanding of expression patterns that define healthy or diseased islet cell types and helped explicate some of the intricacies between major islet cell crosstalk and glucose regulation. All three endocrine cell types derive from a common pancreatic progenitor, yet alpha and beta cells have partially opposing functions, and delta cells modulate and control insulin and glucagon release. While gene expression signatures that define and maintain cellular identity have been widely explored, the underlying epigenetic components are incompletely characterized and understood. However, chromatin accessibility and remodeling is a dynamic attribute that plays a critical role to determine and maintain cellular identity.

RESULTS

Here, we compare and contrast the chromatin landscape between mouse alpha, beta, and delta cells using ATAC-Seq to evaluate the significant differences in chromatin accessibility. The similarities and differences in chromatin accessibility between these related islet endocrine cells help define their fate in support of their distinct functional roles. We identify patterns that suggest that both alpha and delta cells are poised, but repressed, from becoming beta-like. We also identify patterns in differentially enriched chromatin that have transcription factor motifs preferentially associated with different regions of the genome. Finally, we not only confirm and visualize previously discovered common endocrine- and cell specific- enhancer regions across differentially enriched chromatin, but identify novel regions as well. We compiled our chromatin accessibility data in a freely accessible database of common endocrine- and cell specific-enhancer regions that can be navigated with minimal bioinformatics expertise.

CONCLUSIONS

Both alpha and delta cells appear poised, but repressed, from becoming beta cells in murine pancreatic islets. These data broadly support earlier findings on the plasticity in identity of non-beta cells under certain circumstances. Furthermore, differential chromatin accessibility shows preferentially enriched distal-intergenic regions in beta cells, when compared to either alpha or delta cells.

Аpplication of massive parallel reporter analysis in biotechnology and medicine

The development and functioning of an organism relies on tissue-specific gene programs. Genome regulatory elements play a key role in the regulation of such programs, and disruptions in their function can lead to the development of various pathologies, including cancers, malformations and autoimmune diseases. The emergence of high-throughput genomic studies has led to massively parallel reporter analysis (MPRA) methods, which allow the functional verification and identification of regulatory elements on a genome-wide scale. Initially MPRA was used as a tool to investigate fundamental aspects of epigenetics, but the approach also has great potential for clinical and practical biotechnology. Currently, MPRA is used for validation of clinically significant mutations, identification of tissue-specific regulatory elements, search for the most promising loci for transgene integration, and is an indispensable tool for creating highly efficient expression systems, the range of application of which extends from approaches for protein development and design of next-generation therapeutic antibody superproducers to gene therapy. In this review, the main principles and areas of practical application of high-throughput reporter assays will be discussed.

Accurate prediction of functional states of cis-regulatory modules reveals common epigenetic rules in humans and mice

BACKGROUND

Predicting cis-regulatory modules (CRMs) in a genome and their functional states in various cell/tissue types of the organism are two related challenging computational tasks. Most current methods attempt to simultaneously achieve both using data of multiple epigenetic marks in a cell/tissue type. Though conceptually attractive, they suffer high false discovery rates and limited applications. To fill the gaps, we proposed a two-step strategy to first predict a map of CRMs in the genome, and then predict functional states of all the CRMs in various cell/tissue types of the organism. We have recently developed an algorithm for the first step that was able to more accurately and completely predict CRMs in a genome than existing methods by integrating numerous transcription factor ChIP-seq datasets in the organism. Here, we presented machine-learning methods for the second step.

RESULTS

We showed that functional states in a cell/tissue type of all the CRMs in the genome could be accurately predicted using data of only 1~4 epigenetic marks by a variety of machine-learning classifiers. Our predictions are substantially more accurate than the best achieved so far. Interestingly, a model trained on a cell/tissue type in humans can accurately predict functional states of CRMs in different cell/tissue types of humans as well as of mice, and vice versa. Therefore, epigenetic code that defines functional states of CRMs in various cell/tissue types is universal at least in humans and mice. Moreover, we found that from tens to hundreds of thousands of CRMs were active in a human and mouse cell/tissue type, and up to 99.98% of them were reutilized in different cell/tissue types, while as small as 0.02% of them were unique to a cell/tissue type that might define the cell/tissue type.

CONCLUSIONS

Our two-step approach can accurately predict functional states in any cell/tissue type of all the CRMs in the genome using data of only 1~4 epigenetic marks. Our approach is also more cost-effective than existing methods that typically use data of more epigenetic marks. Our results suggest common epigenetic rules for defining functional states of CRMs in various cell/tissue types in humans and mice.

Identification of putative enhancer-like elements predicts regulatory networks active in planarian adult stem cells

Planarians have become an established model system to study regeneration and stem cells, but the regulatory elements in the genome remain almost entirely undescribed. Here, by integrating epigenetic and expression data we use multiple sources of evidence to predict enhancer elements active in the adult stem cell populations that drive regeneration. We have used ChIP-seq data to identify genomic regions with histone modifications consistent with enhancer activity, and ATAC-seq data to identify accessible chromatin. Overlapping these signals allowed for the identification of a set of high-confidence candidate enhancers predicted to be active in planarian adult stem cells. These enhancers are enriched for predicted transcription factor (TF) binding sites for TFs and TF families expressed in planarian adult stem cells. Footprinting analyses provided further evidence that these potential TF binding sites are likely to be occupied in adult stem cells. We integrated these analyses to build testable hypotheses for the regulatory function of TFs in stem cells, both with respect to how pluripotency might be regulated, and to how lineage differentiation programs are controlled. We found that our predicted GRNs were independently supported by existing TF RNAi/RNA-seq datasets, providing further evidence that our work predicts active enhancers that regulate adult stem cells and regenerative mechanisms.

M2-like tumor-associated macrophages transmit exosomal miR-27b-3p and maintain glioblastoma stem-like cell properties

There is growing evidence supporting the implications of exosomes-shuttled microRNAs (miRs) in the phenotypes of glioblastoma stem cells (GSCs), whilst the role of exosomal miR-27b-3p remains to be established. Herein, the aim of this study was to investigate the effect of M2 tumor-associated macrophage (TAM)-derived exosomal miR-27b-3p on the function of GSCs. Clinical glioblastoma (GBM) specimens were obtained and GSCs and M2-TAMs were isolated by fluorescence-activated cell sorting (FACS), and exosomes were separated from M2-TAMs. It was observed that M2-TAM-derived exosomes promoted the stem-like properties of GSCs. Gain- and loss- of function assays were then conducted to explore the effects of exosomal miR-27b-3p and the miR-27b-3p/MLL4/PRDM1 axis on GSC phenotypes. A xenograft tumor model of GBM was further established for in vivo substantiation. Inhibition of miR-27b-3p in M2-TAMs reduced exosomal miR-27b-3p transferred into GSCs and consequently diminished GSC viability in vitro and tumor-promoting effects of GSCs in vivo. The interaction among miR-27b-3p, mixed linked leukemia 4 (MLL4), positive regulatory domain I (PRDM1) was validated by dual-luciferase and ChIP assays. MLL4 positively regulated PRDM1 expression by inducing methylation in the PRDM1 enhancer region and ultimately reduced IL-33 expression. miR-27b-3p targeted MLL4/PRDM1 to activate IL-33 and maintain the stem-like function of GSCs. In conclusion, our study elucidated that M2-TAM-derived exosomal miR-27b-3p enhanced the tumorigenicity of GSCs through the MLL4/PRDM1/IL-33 axis.

Establishing gene regulatory networks from Parkinson's disease risk loci

The latest meta-analysis of genome-wide association studies identified 90 independent variants across 78 genomic regions associated with Parkinson's disease, yet the mechanisms by which these variants influence the development of the disease remains largely elusive. To establish the functional gene regulatory networks associated with Parkinson's disease risk variants, we utilized an approach combining spatial (chromosomal conformation capture) and functional (expression quantitative trait loci) data. We identified 518 genes subject to regulation by 76 Parkinson's variants across 49 tissues, whicih encompass 36 peripheral and 13 CNS tissues. Notably, one-third of these genes were regulated via trans-acting mechanisms (distal; risk locus-gene separated by >1 Mb, or on different chromosomes). Of particular interest is the identification of a novel trans-expression quantitative trait loci-gene connection between rs10847864 and SYNJ1 in the adult brain cortex, highlighting a convergence between familial studies and Parkinson's disease genome-wide association studies loci for SYNJ1 (PARK20) for the first time. Furthermore, we identified 16 neurodevelopment-specific expression quantitative trait loci-gene regulatory connections within the foetal cortex, consistent with hypotheses suggesting a neurodevelopmental involvement in the pathogenesis of Parkinson's disease. Through utilizing Louvain clustering we extracted nine significant and highly intraconnected clusters within the entire gene regulatory network. The nine clusters are enriched for specific biological processes and pathways, some of which have not previously been associated with Parkinson's disease. Together, our results not only contribute to an overall understanding of the mechanisms and impact of specific combinations of Parkinson's disease variants, but also highlight the potential impact gene regulatory networks may have when elucidating aetiological subtypes of Parkinson's disease.

Analysis of the landscape of human enhancer sequences in biological databases

The process of gene regulation extends as a network in which both genetic sequences and proteins are involved. The levels of regulation and the mechanisms involved are multiple. Transcription is the main control mechanism for most genes, being the downstream steps responsible for refining the transcription patterns. In turn, gene transcription is mainly controlled by regulatory events that occur at promoters and enhancers. Several studies are focused on analyzing the contribution of enhancers in the development of diseases and their possible use as therapeutic targets. The study of regulatory elements has advanced rapidly in recent years with the development and use of next generation sequencing techniques. All this information has generated a large volume of information that has been transferred to a growing number of public repositories that store this information. In this article, we analyze the content of those public repositories that contain information about human enhancers with the aim of detecting whether the knowledge generated by scientific research is contained in those databases in a way that could be computationally exploited. The analysis will be based on three main aspects identified in the literature: types of enhancers, type of evidence about the enhancers, and methods for detecting enhancer-promoter interactions. Our results show that no single database facilitates the optimal exploitation of enhancer data, most types of enhancers are not represented in the databases and there is need for a standardized model for enhancers. We have identified major gaps and challenges for the computational exploitation of enhancer data.

The Hippo effector YAP1/TEAD1 regulates EPHA3 expression to control cell contact and motility

The EPHA3 protein tyrosine kinase, a member of the ephrin receptor family, regulates cell fate, cell motility, and cell-cell interaction. These cellular events are critical for tissue development, immunological responses, and the processes of tumorigenesis. Earlier studies revealed that signaling via the STK4-encoded MST1 serine-threonine protein kinase, a core component of the Hippo pathway, attenuated EPHA3 expression. Here, we investigated the mechanism by which MST1 regulates EPHA3. Our findings have revealed that the transcriptional regulators YAP1 and TEAD1 are crucial activators of EPHA3 transcription. Silencing YAP1 and TEAD1 suppressed the EPHA3 protein and mRNA levels. In addition, we identified putative TEAD enhancers in the distal EPHA3 promoter, where YAP1 and TEAD1 bind and promote EPHA3 expression. Furthermore, EPHA3 knockout by CRISPR/Cas9 technology reduced cell-cell interaction and cell motility. These findings demonstrate that EPHA3 is transcriptionally regulated by YAP1/TEAD1 of the Hippo pathway, suggesting that it is sensitive to cell contact-dependent interactions.

Molecular Mechanisms of lncRNAs in the Dependent Regulation of Cancer and Their Potential Therapeutic Use

Deep whole genome and transcriptome sequencing have highlighted the importance of an emerging class of non-coding RNA longer than 200 nucleotides (i.e., long non-coding RNAs (lncRNAs)) that are involved in multiple cellular processes such as cell differentiation, embryonic development, and tissue homeostasis. Cancer is a prime example derived from a loss of homeostasis, primarily caused by genetic alterations both in the genomic and epigenetic landscape, which results in deregulation of the gene networks. Deregulation of the expression of many lncRNAs in samples, tissues or patients has been pointed out as a molecular regulator in carcinogenesis, with them acting as oncogenes or tumor suppressor genes. Herein, we summarize the distinct molecular regulatory mechanisms described in literature in which lncRNAs modulate carcinogenesis, emphasizing epigenetic and genetic alterations in particular. Furthermore, we also reviewed the current strategies used to block lncRNA oncogenic functions and their usefulness as potential therapeutic targets in several carcinomas.

Genome-wide identification of enhancers and transcription factors regulating the myogenic differentiation of bovine satellite cells

BACKGROUND

Satellite cells are the myogenic precursor cells in adult skeletal muscle. The objective of this study was to identify enhancers and transcription factors that regulate gene expression during the differentiation of bovine satellite cells into myotubes.

RESULTS

Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) was performed to identify genomic regions where lysine 27 of H3 histone is acetylated (H3K27ac), i.e., active enhancers, from bovine satellite cells before and during differentiation into myotubes. A total of 19,027 and 47,669 H3K27ac-marked enhancers were consistently identified from two biological replicates of before- and during-differentiation bovine satellite cells, respectively. Of these enhancers, 5882 were specific to before-differentiation, 35,723 to during-differentiation, and 13,199 common to before- and during-differentiation bovine satellite cells. Whereas most of the before- or during-differentiation-specific H3K27ac-marked enhancers were located distally to the transcription start site, the enhancers common to before- and during-differentiation were located both distally and proximally to the transcription start site. The three sets of H3K27ac-marked enhancers were associated with functionally different genes and enriched with different transcription factor binding sites. Specifically, many of the H3K27ac-marked enhancers specific to during-differentiation bovine satellite cells were associated with genes involved in muscle structure and development, and were enriched with binding sites for the MyoD, AP-1, KLF, TEAD, and MEF2 families of transcription factors. A positive role was validated for Fos and FosB, two AP-1 family transcription factors, in the differentiation of bovine satellite cells into myotubes by siRNA-mediated knockdown.

CONCLUSIONS

Tens of thousands of H3K27ac-marked active enhancers have been identified from bovine satellite cells before or during differentiation. These enhancers contain binding sites not only for transcription factors whose role in satellite cell differentiation is well known but also for transcription factors whose role in satellite cell differentiation is unknown. These enhancers and transcription factors are valuable resources for understanding the complex mechanism that mediates gene expression during satellite cell differentiation. Because satellite cell differentiation is a key step in skeletal muscle growth, the enhancers, the transcription factors, and their target genes identified in this study are also valuable resources for identifying and interpreting skeletal muscle trait-associated DNA variants in cattle.

Cell division- and DNA replication-free reprogramming of somatic nuclei for embryonic transcription

Nuclear transfer systems represent the efficient means to reprogram a cell and in theory provide a basis for investigating the development of endangered species. However, conventional nuclear transfer using oocytes of laboratory animals does not allow reprogramming of cross-species nuclei owing to defects in cell divisions and activation of embryonic genes. Here, we show that somatic nuclei transferred into mouse four-cell embryos arrested at the G2/M phase undergo reprogramming toward the embryonic state. Remarkably, genome-wide transcriptional reprogramming is induced within a day, and ZFP281 is important for this replication-free reprogramming. This system further enables transcriptional reprogramming of cells from Oryx dammah, now extinct in the wild. Thus, our findings indicate that arrested mouse embryos are competent to induce intra- and cross-species reprogramming. The direct induction of embryonic transcripts from diverse genomes paves a unique approach for identifying mechanisms of transcriptional reprogramming and genome activation from a diverse range of species.

A dominant-negative SOX18 mutant disrupts multiple regulatory layers essential to transcription factor activity

Few genetically dominant mutations involved in human disease have been fully explained at the molecular level. In cases where the mutant gene encodes a transcription factor, the dominant-negative mode of action of the mutant protein is particularly poorly understood. Here, we studied the genome-wide mechanism underlying a dominant-negative form of the SOX18 transcription factor (SOX18^RaOp) responsible for both the classical mouse mutant Ragged Opossum and the human genetic disorder Hypotrichosis-lymphedema-telangiectasia-renal defect syndrome. Combining three single-molecule imaging assays in living cells together with genomics and proteomics analysis, we found that SOX18^RaOp disrupts the system through an accumulation of molecular interferences which impair several functional properties of the wild-type SOX18 protein, including its target gene selection process. The dominant-negative effect is further amplified by poisoning the interactome of its wild-type counterpart, which perturbs regulatory nodes such as SOX7 and MEF2C. Our findings explain in unprecedented detail the multi-layered process that underpins the molecular aetiology of dominant-negative transcription factor function.

Enhancers in disease: molecular basis and emerging treatment strategies

Enhancers are genomic sequences that play a key role in regulating tissue-specific gene expression levels. An increasing number of diseases are linked to impaired enhancer function through chromosomal rearrangement, genetic variation within enhancers, or epigenetic modulation. Here, we review how these enhancer disruptions have recently been implicated in congenital disorders, cancers, and common complex diseases and address the implications for diagnosis and treatment. Although further fundamental research into enhancer function, target genes, and context is required, enhancer-targeting drugs and gene editing approaches show great therapeutic promise for a range of diseases.

The role of CTCF in the organization of the centromeric 11p15 imprinted domain interactome

DNA methylation, chromatin-binding proteins, and DNA looping are common components regulating genomic imprinting which leads to parent-specific monoallelic gene expression. Loss of methylation (LOM) at the human imprinting center 2 (IC2) on chromosome 11p15 is the most common cause of the imprinting overgrowth disorder Beckwith-Wiedemann Syndrome (BWS). Here, we report a familial transmission of a 7.6 kB deletion that ablates the core promoter of KCNQ1. This structural alteration leads to IC2 LOM and causes recurrent BWS. We find that occupancy of the chromatin organizer CTCF is disrupted proximal to the deletion, which causes chromatin architecture changes both in cis and in trans. We also profile the chromatin architecture of IC2 in patients with sporadic BWS caused by isolated LOM to identify conserved features of IC2 regulatory disruption. A strong interaction between CTCF sites around KCNQ1 and CDKN1C likely drive their expression on the maternal allele, while a weaker interaction involving the imprinting control region element may impede this connection and mediate gene silencing on the paternal allele. We present an imprinting model in which KCNQ1 transcription is necessary for appropriate CTCF binding and a novel chromatin conformation to drive allele-specific gene expression.

TempoMAGE: a deep learning framework that exploits the causal dependency between time-series data to predict histone marks in open chromatin regions at time-points with missing ChIP-seq datasets

MOTIVATION

Identifying histone tail modifications using ChIP-seq is commonly used in time-series experiments in development and disease. These assays, however, cover specific time-points leaving intermediate or early stages with missing information. Although several machine learning methods were developed to predict histone marks, none exploited the dependence that exists in time-series experiments between data generated at specific time-points to extrapolate these findings to time-points where data cannot be generated for lack or scarcity of materials (i.e., early developmental stages).

RESULTS

Here, we train a deep learning model named TempoMAGE, to predict the presence or absence of H3K27ac in open chromatin regions by integrating information from sequence, gene expression, chromatin accessibility and the estimated change in H3K27ac state from a reference time-point. We show that adding reference time-point information systematically improves the overall model's performance. Additionally, sequence signatures extracted from our method were exclusive to the training dataset indicating that our model learned data-specific features. As an application, TempoMAGE was able to predict the activity of enhancers from pre-validated in-vivo dataset highlighting its ability to be used for functional annotation of putative enhancers.

AVAILABILITY

TempoMAGE is freely available through GitHub at https://github.com/pkhoueiry/TempoMAGE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Preclinical evidence for the therapeutic value of TBX5 normalization in arrhythmia control

AIMS

Arrhythmias and sudden cardiac death (SCD) occur commonly in patients with heart failure. We found T-box 5 (TBX5) dysregulated in ventricular myocardium from heart failure patients and thus we hypothesized that TBX5 reduction contributes to arrhythmia development in these patients. To understand the underlying mechanisms, we aimed to reveal the ventricular TBX5-dependent transcriptional network and further test the therapeutic potential of TBX5 level normalization in mice with documented arrhythmias.

METHODS AND RESULTS

We used a mouse model of TBX5 conditional deletion in ventricular cardiomyocytes. Ventricular (v) TBX5 loss in mice resulted in mild cardiac dysfunction and arrhythmias and was associated with a high mortality rate (60%) due to SCD. Upon angiotensin stimulation, vTbx5KO mice showed exacerbated cardiac remodelling and dysfunction suggesting a cardioprotective role of TBX5. RNA-sequencing of a ventricular-specific TBX5KO mouse and TBX5 chromatin immunoprecipitation was used to dissect TBX5 transcriptional network in cardiac ventricular tissue. Overall, we identified 47 transcripts expressed under the control of TBX5, which may have contributed to the fatal arrhythmias in vTbx5KO mice. These included transcripts encoding for proteins implicated in cardiac conduction and contraction (Gja1, Kcnj5, Kcng2, Cacna1g, Chrm2), in cytoskeleton organization (Fstl4, Pdlim4, Emilin2, Cmya5), and cardiac protection upon stress (Fhl2, Gpr22, Fgf16). Interestingly, after TBX5 loss and arrhythmia development in vTbx5KO mice, TBX5 protein-level normalization by systemic adeno-associated-virus (AAV) 9 application, re-established TBX5-dependent transcriptome. Consequently, cardiac dysfunction was ameliorated and the propensity of arrhythmia occurrence was reduced.

CONCLUSIONS

This study uncovers a novel cardioprotective role of TBX5 in the adult heart and provides preclinical evidence for the therapeutic value of TBX5 protein normalization in the control of arrhythmia.

Putative Causal Variants Are Enriched in Annotated Functional Regions From Six Bovine Tissues

Genetic variants which affect complex traits (causal variants) are thought to be found in functional regions of the genome. Identifying causal variants would be useful for predicting complex trait phenotypes in dairy cows, however, functional regions are poorly annotated in the bovine genome. Functional regions can be identified on a genome-wide scale by assaying for post-translational modifications to histone proteins (histone modifications) and proteins interacting with the genome (e.g., transcription factors) using a method called Chromatin immunoprecipitation followed by sequencing (ChIP-seq). In this study ChIP-seq was performed to find functional regions in the bovine genome by assaying for four histone modifications (H3K4Me1, H3K4Me3, H3K27ac, and H3K27Me3) and one transcription factor (CTCF) in 6 tissues (heart, kidney, liver, lung, mammary and spleen) from 2 to 3 lactating dairy cows. Eighty-six ChIP-seq samples were generated in this study, identifying millions of functional regions in the bovine genome. Combinations of histone modifications and CTCF were found using ChromHMM and annotated by comparing with active and inactive genes across the genome. Functional marks differed between tissues highlighting areas which might be particularly important to tissue-specific regulation. Supporting the cis-regulatory role of functional regions, the read counts in some ChIP peaks correlated with nearby gene expression. The functional regions identified in this study were enriched for putative causal variants as seen in other species. Interestingly, regions which correlated with gene expression were particularly enriched for potential causal variants. This supports the hypothesis that complex traits are regulated by variants that alter gene expression. This study provides one of the largest ChIP-seq annotation resources in cattle including, for the first time, in the mammary gland of lactating cows. By linking regulatory regions to expression QTL and trait QTL we demonstrate a new strategy for identifying causal variants in cattle.

Novel lncRNAs Co-Expression Networks Identifies LINC00504 with Oncogenic Role in Luminal A Breast Cancer Cells

Long non-coding RNAs (lncRNAs) are functional transcripts with more than 200 nucleotides. These molecules exhibit great regulatory capacity and may act at different levels of gene expression regulation. Despite this regulatory versatility, the biology of these molecules is still poorly understood. Computational approaches are being increasingly used to elucidate biological mechanisms in which these lncRNAs may be involved. Co-expression networks can serve as great allies in elucidating the possible regulatory contexts in which these molecules are involved. Herein, we propose the use of the pipeline deposited in the RTN package to build lncRNAs co-expression networks using TCGA breast cancer (BC) cohort data. Worldwide, BC is the most common cancer in women and has great molecular heterogeneity. We identified an enriched co-expression network for the validation of relevant cell processes in the context of BC, including LINC00504. This lncRNA has increased expression in luminal subtype A samples, and is associated with prognosis in basal-like subtype. Silencing this lncRNA in luminal A cell lines resulted in decreased cell viability and colony formation. These results highlight the relevance of the proposed method for the identification of lncRNAs in specific biological contexts.

Predicting Genome Architecture: Challenges and Solutions

Genome architecture plays a pivotal role in gene regulation. The use of high-throughput methods for chromatin profiling and 3-D interaction mapping provide rich experimental data sets describing genome organization and dynamics. These data challenge development of new models and algorithms connecting genome architecture with epigenetic marks. In this review, we describe how chromatin architecture could be reconstructed from epigenetic data using biophysical or statistical approaches. We discuss the applicability and limitations of these methods for understanding the mechanisms of chromatin organization. We also highlight the emergence of new predictive approaches for scoring effects of structural variations in human cells.

Genome-wide analysis of chromatin structure changes upon MyoD binding in proliferative myoblasts during the cell cycle

MyoD, a myogenic differentiation protein, has been studied for its critical role in skeletal muscle differentiation. MyoD-expressing myoblasts have a potency to be differentiated with proliferation of ectopic cells. However, little is known about the effect on chromatin structure of MyoD binding in proliferative myoblasts. In this study, we evaluated the chromatin structure around MyoD-bound genome regions during the cell cycle by chromatin immunoprecipitation sequencing. Genome-wide analysis of histone modifications was performed in proliferative mouse C2C12 myoblasts during three phases (G1, S, G2/M) of the cell cycle. We found that MyoD-bound genome regions had elevated levels of active histone modifications, such as H3K4me1/2/3, and H3K27ac, compared with MyoD-unbound genome regions during the cell cycle. We also demonstrated that the elevated H3K4me2/3 modification level was maintained during the cell cycle, whereas the H3K27ac and H3K4me1 modification levels decreased to the same level as MyoD-unbound genome regions during the later phases. Immunoblot analysis revealed that MyoD abundance was high in the G1 phase then decreased in the S and G2/M phases. Our results suggest that MyoD binding formed selective epigenetic memories with H3K4me2/3 during the cell cycle in addition to myogenic gene induction via active chromatin formation coupled with transcription.

Bivalent Genes Targeting of Glioma Heterogeneity and Plasticity

Gliomas account for most primary Central Nervous System (CNS) neoplasms, characterized by high aggressiveness and low survival rates. Despite the immense research efforts, there is a small improvement in glioma survival rates, mostly attributed to their heterogeneity and complex pathophysiology. Recent data indicate the delicate interplay of genetic and epigenetic mechanisms in regulating gene expression and cell differentiation, pointing towards the pivotal role of bivalent genes. Bivalency refers to a property of chromatin to acquire more than one histone marks during the cell cycle and rapidly transition gene expression from an active to a suppressed transcriptional state. Although first identified in embryonal stem cells, bivalent genes have now been associated with tumorigenesis and cancer progression. Emerging evidence indicates the implication of bivalent gene regulation in glioma heterogeneity and plasticity, mainly involving Homeobox genes, Wingless-Type MMTV Integration Site Family Members, Hedgehog protein, and Solute Carrier Family members. These genes control a wide variety of cellular functions, including cellular differentiation during early organism development, regulation of cell growth, invasion, migration, angiogenesis, therapy resistance, and apoptosis. In this review, we discuss the implication of bivalent genes in glioma pathogenesis and their potential therapeutic targeting options.

NGS-Integrator: An efficient tool for combining multiple NGS data tracks using minimum Bayes' factors

BACKGROUND

Next-generation sequencing (NGS) is widely used for genome-wide identification and quantification of DNA elements involved in the regulation of gene transcription. Studies that generate multiple high-throughput NGS datasets require data integration methods for two general tasks: 1) generation of genome-wide data tracks representing an aggregate of multiple replicates of the same experiment; and 2) combination of tracks from different experimental types that provide complementary information regarding the location of genomic features such as enhancers.

RESULTS

NGS-Integrator is a Java-based command line application, facilitating efficient integration of multiple genome-wide NGS datasets. NGS-Integrator first transforms all input data tracks using the complement of the minimum Bayes' factor so that all values are expressed in the range [0,1] representing the probability of a true signal given the background noise. Then, NGS-Integrator calculates the joint probability for every genomic position to create an integrated track. We provide examples using real NGS data generated in our laboratory and from the mouse ENCODE database.

CONCLUSIONS

Our results show that NGS-Integrator is both time- and memory-efficient. Our examples show that NGS-Integrator can integrate information to facilitate downstream analyses that identify functional regulatory domains along the genome.

H. pylori infection confers resistance to apoptosis via Brd4-dependent BIRC3 eRNA synthesis

H. pylori infection is one of the leading causes of gastric cancer and the pathogenicity of H. pylori infection is associated with its ability to induce chronic inflammation and apoptosis resistance. While H. pylori infection-induced expression of pro-inflammatory cytokines for chronic inflammation is well studied, the molecular mechanism underlying the apoptosis resistance in infected cells is not well understood. In this study, we demonstrated that H. pylori infection-induced apoptosis resistance in gastric epithelial cells triggered by Raptinal, a drug that directly activates caspase-3. This resistance resulted from the induction of cIAP2 (encoded by BIRC3) since depletion of BIRC3 by siRNA or inhibition of cIAP2 via BV6 reversed H. pylori-suppressed caspase-3 activation. The induction of cIAP2 was regulated by H. pylori-induced BIRC3 eRNA synthesis. Depletion of BIRC3 eRNA decreased H. pylori-induced cIAP2 and reversed H. pylori-suppressed caspase-3 activation. Mechanistically, H. pylori stimulated the recruitment of bromodomain-containing factor Brd4 to the enhancer of BIRC3 and promoted BIRC3 eRNA and mRNA synthesis. Inhibition of Brd4 diminished the expression of BIRC3 eRNA and the anti-apoptotic response to H. pylori infection. Importantly, H. pylori isogenic cagA-deficient mutant failed to activate the synthesis of BIRC3 eRNA and the associated apoptosis resistance. Finally, in primary human gastric epithelial cells, H. pylori also induced resistance to Raptinal-triggered caspase-3 activation by activating the Brd4-dependent BIRC3 eRNA synthesis in a CagA-dependent manner. These results identify a novel function of Brd4 in H. pylori-mediated apoptosis resistance via activating BIRC3 eRNA synthesis, suggesting that Brd4 could be a potential therapeutic target for H. pylori-induced gastric cancer.

Promoter anchored interaction landscape of THP-1 macrophages captures early immune response processes

Macrophages are highly plastic immune cells with temporally distinct transcriptome changes upon lipopolysaccride (LPS) activation. However, to what extent transcriptome reprogramming is mediated via spatial chromatin looping is not well studied. We generated high resolution chromatin interaction maps for LPS-stimulated THP-1 macrophages (0 and 2 h) using capture Hi-C. Success of LPS stimulation was validated with transcriptome sequencing. Circa 2900 genes changed their interaction profile upon LPS stimulation and those gaining interactions were enriched for LPS response relevant processes, suggesting a substantial role for distal regulation. Immune and cardiovascular risk variants were enriched within the interacting regions, thereby providing insights into macrophage biology.

Embryonic Program Activated during Blast Crisis of Chronic Myelogenous Leukemia (CML) Implicates a TCF7L2 and MYC Cooperative Chromatin Binding

Chronic myeloid leukemia (CML) is characterized by an inherent genetic instability, which contributes to the progression of the disease towards an accelerated phase (AP) and blast crisis (BC). Several cytogenetic and genomic alterations have been reported in the progression towards BC, but the precise molecular mechanisms of this event are undetermined. Transcription Factor 7 like 2 (TFC7L2) is a member of the TCF family of proteins that are known to activate WNT target genes such as Cyclin D1. TCF7L2 has been shown to be overexpressed in acute myeloid leukemia (AML) and represents a druggable target. We report here that TCF7L2 transcription factor expression was found to be correlated to blast cell numbers during the progression of the disease. In these cells, TCF7L2 CHIP-sequencing highlighted distal cis active enhancer, such as elements in SMAD3, ATF5, and PRMT1 genomic regions and a proximal active transcriptional program of 144 genes. The analysis of CHIP-sequencing of MYC revealed a significant overlapping of TCF7L2 epigenetic program with MYC. The β-catenin activator lithium chloride and the MYC-MAX dimerization inhibitor 10058-F4 significantly modified the expression of three epigenetic targets in the BC cell line K562. These results suggest for the first time the cooperative role of TCF7L2 and MYC during CML-BC and they strengthen previous data showing a possible involvement of embryonic genes in this process.

Keeping RNA polymerase II on the run: Functions of MLL fusion partners in transcriptional regulation

Since the identification of key MLL fusion partners as transcription elongation factors regulating expression of HOX cluster genes during hematopoiesis, extensive work from the last decade has resulted in significant progress in our overall mechanistic understanding of role of MLL fusion partner proteins in transcriptional regulation of diverse set of genes beyond just the HOX cluster. In this review, we are going to detail overall understanding of role of MLL fusion partner proteins in transcriptional regulation and thus provide mechanistic insights into possible MLL fusion protein-mediated transcriptional misregulation leading to aberrant hematopoiesis and leukemogenesis.

Profiling of histone 3 lysine 27 acetylation reveals its role in a chronic DSS-induced colitis mouse model

Inflammatory bowel disease (IBD) is a chronic inflammatory condition of the gastrointestinal tract. In current dogma, pathogenesis of IBD is attributed to the dysregulated mucosal immune response to gut flora in genetically susceptible individuals, but the genetics evidence from GWAS studies so far is insufficient to explain the observed heritability in IBD. For this discordance, epigenetics has emerged to be one of the important causes. Recent studies have reported that histone acetylation is correlated with the development of IBD, whereas its role and underlying molecular mechanism in the disease still remain elusive. Here, we established a dextran sulfate sodium (DSS)-induced chronic colitis model and performed RNA-sequencing (RNA-seq) and Chromatin Immunoprecipitation followed by NGS sequencing (ChIP-seq) for H3K27ac in the mice colon tissues to investigate whether H3K27ac is involved in the development of intestinal inflammation. We found that the global H3K27ac level and distribution in colon tissue had no significant difference after DSS treatment, while H3K27ac signals were significantly enriched in the typical-enhancers of the DSS group compared with the control. By combining with RNA-seq data (fold change >2), we identified 56 candidate genes as potential target genes for H3K27ac change upon DSS treatment. We further predicted transcription factors (TFs) involved in DSS-induced colitis according to the enhancers with increased H3K27ac. H3K27ac increase in special typical-enhancers in the DSS group is possibly related to the development of intestinal inflammation by up-regulating adjacent gene expression and shifting TF networks, which will provide new insight into the pathogenesis and therapy of IBD.

Chromatin Dynamics and Genetic Variation Combine to Regulate Innate Immune Memory

Recent work by Ciernia et al. (2020) identified how genetic and epigenetic mechanisms interact to regulate innate immune memory in bone marrow derived macrophages. The authors examined the BTBR strain, a naturally occurring mouse model of Autism Spectrum Disorder (ASD) that captures the complex genetics, behavioral and immune dysregulation found in the human disorder. Immune cell cultures from the BTBR strain compared to the standard C57 showed hyper-responsive immune gene expression that was linked to altered chromatin accessibility at sites with genetic differences between the strains. Together, findings from this work demonstrated that multiple levels of gene regulation likely dictate the formation of innate immune memory and are likely disrupted in immune cells in ASD. Future work will be needed to extend these findings to immune gene regulation in the brain and how changes in immune function are related to abnormal behaviors in brain disorders.

Neurobiological functions of transcriptional enhancers

Transcriptional enhancers are regulatory DNA elements that underlie the specificity and dynamic patterns of gene expression. Over the past decade, large-scale functional genomics projects have driven transformative progress in our understanding of enhancers. These data have relevance for identifying mechanisms of gene regulation in the CNS, elucidating the function of non-coding regulatory sequences in neurobiology and linking sequence variation within enhancers to genetic risk for neurological and psychiatric disorders. However, the sheer volume and complexity of genomic data presents a challenge to interpreting enhancer function in normal and pathogenic neurobiological processes. Here, to advance the application of genome-scale enhancer data, we offer a primer on current models of enhancer function in the CNS, we review how enhancers regulate gene expression across the neuronal lifespan, and we suggest how emerging findings regarding the role of non-coding sequence variation offer opportunities for understanding brain disorders and developing new technologies for neuroscience.

Association between genetic polymorphisms of NRF2, KEAP1, MAFF, MAFK and anti-tuberculosis drug-induced liver injury: a nested case-control study

Reactive metabolites of anti-tuberculosis (anti-TB) drugs can result in excessive reactive oxygen species (ROS), which are responsible for drug-induced liver injury. The nuclear factor erythroid 2-related factor 2 (Nrf2) - antioxidant response elements (ARE) (Nrf2-ARE) signaling pathway plays a crucial role in protecting liver cells from ROS, inducing enzymes such as phase II metabolizing enzymes and antioxidant enzymes. Based on a Chinese anti-TB treatment cohort, a nested case-control study was performed to explore the association between 13 tag single-nucleotide polymorphisms (tagSNPs) in the NRF2, KEAP1, MAFF, MAFK genes in Nrf2-ARE signaling pathway and the risk of anti-TB drug-induced liver injury (ATLI) in 314 cases and 628 controls. Conditional logistic regression models were used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) after adjusting weight and usage of hepatoprotectant. Patients carrying the TC genotype at rs4243387 or haplotype C-C (rs2001350-rs6726395) in NRF2 were at an increased risk of ATLI (adjusted OR = 1.362, 95% CI: 1.017–1.824, P = 0.038; adjusted OR = 2.503, 95% CI: 1.273–4.921, P = 0.008, respectively), whereas patients carrying TC genotype at rs2267373 or haplotype C-G-C (rs2267373-rs4444637-rs4821767) in MAFF were at a reduced risk of ATLI (adjusted OR = 0.712, 95% CI: 0.532–0.953, P = 0.022; adjusted OR = 0.753, 95% CI: 0.587–0.965, P = 0.025, respectively). Subgroup analysis also detected a significant association between multiple tagSNPs (rs4821767 and rs4444637 in MAFF, rs4720833 in MAFK) and specific clinical patterns of liver injury under different genetic models. This study shows that genetic polymorphisms of NRF2, MAFF and MAFK may contribute to the susceptibility to ATLI in the Chinese anti-TB treatment population.

Dissecting Epistatic QTL for Blood Pressure in Rats: Congenic Strains versus Heterogeneous Stocks, a Reality Check

Advances in molecular genetics have provided well-defined physical genetic maps and large numbers of genetic markers for both model organisms and humans. It is now possible to gain a fundamental understanding of the genetic architecture underlying quantitative traits, of which blood pressure (BP) is an important example. This review emphasizes analytical techniques and results obtained using the Dahl salt-sensitive (S) rat as a model of hypertension by presenting results in detail for three specific chromosomal regions harboring genetic elements of increasing complexity controlling BP. These results highlight the critical importance of genetic interactions (epistasis) on BP at all levels of structure, intragenic, intergenic, intrachromosomal, interchromosomal, and across whole genomes. In two of the three examples presented, specific DNA structural variations leading to biochemical, physiological, and pathological mechanisms are well defined. This proves the usefulness of the techniques involving interval mapping followed by substitution mapping using congenic strains. These classic techniques are compared to newer approaches using sophisticated statistical analysis on various segregating or outbred model-organism populations, which in some cases are uniquely useful in demonstrating the existence of higher-order interactions. It is speculated that hypertension as an outlier quantitative phenotype is dependent on higher-order genetic interactions. The obstacle to the identification of genetic elements and the biochemical/physiological mechanisms involved in higher-order interactions is not theoretical or technical but the lack of future resources to finish the job of identifying the individual genetic elements underlying the quantitative trait loci for BP and ascertaining their molecular functions. © 2019 American Physiological Society. Compr Physiol 9:1305-1337, 2019.

Variational infinite heterogeneous mixture model for semi-supervised clustering of heart enhancers

Motivation

Mammalian genomes can contain thousands of enhancers but only a subset are actively driving gene expression in a given cellular context. Integrated genomic datasets can be harnessed to predict active enhancers. One challenge in integration of large genomic datasets is the increasing heterogeneity: continuous, binary and discrete features may all be relevant. Coupled with the typically small numbers of training examples, semi-supervised approaches for heterogeneous data are needed; however, current enhancer prediction methods are not designed to handle heterogeneous data in the semi-supervised paradigm.

Results

We implemented a Dirichlet Process Heterogeneous Mixture model that infers Gaussian, Bernoulli and Poisson distributions over features. We derived a novel variational inference algorithm to handle semi-supervised learning tasks where certain observations are forced to cluster together. We applied this model to enhancer candidates in mouse heart tissues based on heterogeneous features. We constrained a small number of known active enhancers to appear in the same cluster, and 47 additional regions clustered with them. Many of these are located near heart-specific genes. The model also predicted 1176 active promoters, suggesting that it can discover new enhancers and promoters.

Availability and implementation

We created the ‘dphmix’ Python package: https://pypi.org/project/dphmix/.

Supplementary information

Supplementary data are available at Bioinformatics online.

The basic helix-loop-helix transcription factor SHARP1 is an oncogenic driver in MLL-AF6 acute myelogenous leukemia

Acute Myeloid Leukemia (AML) with MLL gene rearrangements demonstrate unique gene expression profiles driven by MLL-fusion proteins. Here, we identify the circadian clock transcription factor SHARP1 as a novel oncogenic target in MLL-AF6 AML, which has the worst prognosis among all subtypes of MLL-rearranged AMLs. SHARP1 is expressed solely in MLL-AF6 AML, and its expression is regulated directly by MLL-AF6/DOT1L. Suppression of SHARP1 induces robust apoptosis of human MLL-AF6 AML cells. Genetic deletion in mice delays the development of leukemia and attenuated leukemia-initiating potential, while sparing normal hematopoiesis. Mechanistically, SHARP1 binds to transcriptionally active chromatin across the genome and activates genes critical for cell survival as well as key oncogenic targets of MLL-AF6. Our findings demonstrate the unique oncogenic role for SHARP1 in MLL-AF6 AML.

Gene fusions involving MLL and different partner genes define unique subgroups of acute myelogenous leukemia, but the mechanisms underlying specific subgroups are not fully clear. Here the authors elucidate the mechanisms of MLL-AF6 induced transformation, providing a distinct pathway that involves SHARP1 as a critical target.

Membrane and Nuclear Estrogen Receptor Alpha Actions: From Tissue Specificity to Medical Implications

Estrogen receptor alpha (ERα) has been recognized now for several decades as playing a key role in reproduction and exerting functions in numerous nonreproductive tissues. In this review, we attempt to summarize the in vitro studies that are the basis of our current understanding of the mechanisms of action of ERα as a nuclear receptor and the key roles played by its two activation functions (AFs) in its transcriptional activities. We then depict the consequences of the selective inactivation of these AFs in mouse models, focusing on the prominent roles played by ERα in the reproductive tract and in the vascular system. Evidence has accumulated over the two last decades that ERα is also associated with the plasma membrane and activates non-nuclear signaling from this site. These rapid/nongenomic/membrane-initiated steroid signals (MISS) have been characterized in a variety of cell lines, and in particular in endothelial cells. The development of selective pharmacological tools that specifically activate MISS and the generation of mice expressing an ERα protein impeded for membrane localization have begun to unravel the physiological role of MISS in vivo. Finally, we discuss novel perspectives for the design of tissue-selective ER modulators based on the integration of the physiological and pathophysiological roles of MISS actions of estrogens.

Changes in Gene Expression and Estrogen Receptor Cistrome in Mouse Liver Upon Acute E2 Treatment

Transcriptional regulation by the estrogen receptor-α (ER) has been investigated mainly in breast cancer cell lines, but estrogens such as 17β-estradiol (E2) exert numerous extrareproductive effects, particularly in the liver, where E2 exhibits both protective metabolic and deleterious thrombotic actions. To analyze the direct and early transcriptional effects of estrogens in the liver, we determined the E2-sensitive transcriptome and ER cistrome in mice after acute administration of E2 or placebo. These analyses revealed the early induction of genes involved in lipid metabolism, which fits with the crucial role of ER in the prevention of liver steatosis. Characterization of the chromatin state of ER binding sites (BSs) in mice expressing or not ER demonstrated that ER is not required per se for the establishment and/or maintenance of chromatin modifications at the majority of its BSs. This is presumably a consequence of a strong overlap between ER and hepatocyte nuclear factor 4α BSs. In contrast, 40% of the BSs of the pioneer factor forkhead box protein a (Foxa2) were dependent upon ER expression, and ER expression also affected the distribution of nucleosomes harboring dimethylated lysine 4 of Histone H3 around Foxa2 BSs. We finally show that, in addition to a network of liver-specific transcription factors including CCAAT/enhancer-binding protein and hepatocyte nuclear factor 4α, ER might be required for proper Foxa2 function in this tissue.

Transcription Factor Hepatocyte Nuclear Factor-1β Regulates Renal Cholesterol Metabolism

HNF-1β is a tissue-specific transcription factor that is expressed in the kidney and other epithelial organs. Humans with mutations in HNF-1β develop kidney cysts, and HNF-1β regulates the transcription of several cystic disease genes. However, the complete spectrum of HNF-1β-regulated genes and pathways is not known. Here, using chromatin immunoprecipitation/next generation sequencing and gene expression profiling, we identified 1545 protein-coding genes that are directly regulated by HNF-1β in murine kidney epithelial cells. Pathway analysis predicted that HNF-1β regulates cholesterol metabolism. Expression of dominant negative mutant HNF-1β or kidney-specific inactivation of HNF-1β decreased the expression of genes that are essential for cholesterol synthesis, including sterol regulatory element binding factor 2 (Srebf2) and 3-hydroxy-3-methylglutaryl-CoA reductase (Hmgcr). HNF-1β mutant cells also expressed lower levels of cholesterol biosynthetic intermediates and had a lower rate of cholesterol synthesis than control cells. Additionally, depletion of cholesterol in the culture medium mitigated the inhibitory effects of mutant HNF-1β on the proteins encoded by Srebf2 and Hmgcr, and HNF-1β directly controlled the renal epithelial expression of proprotein convertase subtilisin-like kexin type 9, a key regulator of cholesterol uptake. These findings reveal a novel role of HNF-1β in a transcriptional network that regulates intrarenal cholesterol metabolism.