Cis-regulatory Regions Research Articles

Mechanistic insights into genetic susceptibility to prostate cancer.

Prostate cancer (PCa) is the second most common cancer in men and is a highly heritable disease that affects millions of individuals worldwide. Genome-wide association studies have to date discovered nearly 270 genetic loci harboring hundreds of single nucleotide polymorphisms (SNPs) that are associated with PCa susceptibility. In contrast, the functional characterization of the mechanisms underlying PCa risk association is still growing. Given that PCa risk-associated SNPs are highly enriched in noncoding cis-regulatory genomic regions, accumulating evidence suggests a widespread modulation of transcription factor chromatin binding and allelic enhancer activity by these noncoding SNPs, thereby dysregulating gene expression. Emerging studies have shown that a proportion of noncoding variants can modulate the formation of transcription factor complexes at enhancers and CTCF-mediated 3D genome architecture. Interestingly, DNA methylation-regulated CTCF binding could orchestrate a long-range chromatin interaction between PCa risk enhancer and causative genes. Additionally, one-causal-variant-two-risk genes or multiple-risk-variant-multiple-genes are prevalent in some PCa risk-associated loci. In this review, we will discuss the current understanding of the general principles of SNP-mediated gene regulation, experimental advances, and functional evidence supporting the mechanistic roles of several PCa genetic loci with potential clinical impact on disease prevention and treatment.

Interspecies transcriptomics identify genes that underlie disproportionate foot growth in jerboas

Despite the great diversity of vertebrate limb proportion and our deep understanding of the genetic mechanisms that drive skeletal elongation, little is known about how individual bones reach different lengths in any species. Here, we directly compare the transcriptomes of homologous growth cartilages of the mouse (Mus musculus) and bipedal jerboa (Jaculus jaculus), the latter of which has "mouse-like" arms but extremely long metatarsals of the feet. Intersecting gene-expression differences in metatarsals and forearms of the two species revealed that about 10% of orthologous genes are associated with the disproportionately rapid elongation of neonatal jerboa feet. These include genes and enriched pathways not previously associated with endochondral elongation as well as those that might diversify skeletal proportion in addition to their known requirements for bone growth throughout the skeleton. We also identified transcription regulators that might act as "nodes" for sweeping differences in genome expression between species. Among these, Shox2, which is necessary for proximal limb elongation, has gained expression in jerboa metatarsals where it has not been detected in other vertebrates. We show that Shox2 is sufficient to increase mouse distal limb length, and a nearby putative cis-regulatory region is preferentially accessible in jerboa metatarsals. In addition to mechanisms that might directly promote growth, we found evidence that jerboa foot elongation may occur in part by de-repressing latent growth potential. The genes and pathways that we identified here provide a framework to understand the modular genetic control of skeletal growth and the remarkable malleability of vertebrate limb proportion.

Mapping cis-regulatory elements in the midgestation mouse placenta

The placenta is a temporary organ that provides the developing fetus with nutrients, oxygen, and protection in utero. Defects in its development, which may be caused by misregulated gene expression, can lead to devastating outcomes for the mother and fetus. In mouse, placental defects during midgestation commonly lead to embryonic lethality. However, the regulatory mechanisms controlling expression of genes during this period have not been thoroughly investigated. Therefore, we generated and analyzed ChIP-seq data for multiple histone modifications known to mark cis-regulatory regions. We annotated active and poised promoters and enhancers, as well as regions generally associated with repressed gene expression. We found that poised promoters were associated with neuronal development genes, while active promoters were largely associated with housekeeping genes. Active and poised enhancers were associated with placental development genes, though only active enhancers were associated with genes that have placenta-specific expression. Motif analysis within active enhancers identified a large network of transcription factors, including those that have not been previously studied in the placenta and are candidates for future studies. The data generated and genomic regions annotated provide researchers with a foundation for future studies, aimed at understanding how specific genes in the midgestation mouse placenta are regulated.

EPCO-05. GENOME-WIDE ANALYSIS OF TEAD1 OCCUPANCY IN BIOLOGICALLY DISTINCT GLIOBLASTOMA SAMPLES

Abstract The diffusely infiltrative nature of glioblastoma (GBM) cells is a major contributor to the disease’s aggressive behavior, including its rapid progression and therapeutic resistance. Moreover, current treatment options do not target the invasive nature of GBM. Recent chromatin accessibility studies prioritized enrichment of the TEAD1 transcription factor motif in glioblastoma stem cell biology and subsequent knockout and overexpression studies confirmed a critical role for TEAD1 in GBM migration, in vitro and in vivo. However, the downstream mechanisms through which TEAD1 regulates GBM cell migration remain poorly understood. In this study, we performed chromatin immunoprecipitation (ChIP-seq) using TEAD1-specific antibody and IgG as non-specific binding control, to characterize TEAD1 occupancy across GBM samples with unique genomic alterations. ChIP-seq peaks were called using MACS2, filtered for duplicates and blacklisted regions, and normalized per sample to their respective genomic input. Initial functional enrichment analyses were performed on three GBM samples with the highest number of TEAD1 occupancy peaks using CistromeGO, which ranked genes based on their TEAD1-specific regulatory potential (RP) score, as a function of peak number and distance from the transcription start site. Analyses of the top 1000 genes with highest TEAD1 RP scores identified 132 common targets across all samples, including known TEAD target genes ETV1 and Cyr61, which related to angiogenesis, cadherin and integrin signaling, cell adhesion, and chromatin regulation gene ontology terms, among others. Interestingly, KEGG pathway analysis also revealed Hippo pathway enrichment across all GBM samples, suggesting a possible TEAD1 regulatory feedback loop in GBM. Analysis of TEAD1 ChIP-seq peaks in non-GBM negative control tissue did not show functional enrichment for any of the terms seen in the GBM samples. Ongoing analyses are focused on characterizing TEAD1 occupancy at active cis-regulatory regions using parallel H3K27ac ChIP-seq data, in order to prioritize the most salient TEAD1-regulatory targets in GBM.

Shaping modern human skull through epigenetic, transcriptional and post-transcriptional regulation of the RUNX2 master bone gene

RUNX2 encodes the master bone transcription factor driving skeletal development in vertebrates, and playing a specific role in craniofacial and skull morphogenesis. The anatomically modern human (AMH) features sequence changes in the RUNX2 locus compared with archaic hominins’ species. We aimed to understand how these changes may have contributed to human skull globularization occurred in recent evolution. We compared in silico AMH and archaic hominins’ genomes, and used mesenchymal stromal cells isolated from skull sutures of craniosynostosis patients for in vitro functional assays. We detected 459 and 470 nucleotide changes in noncoding regions of the AMH RUNX2 locus, compared with the Neandertal and Denisovan genomes, respectively. Three nucleotide changes in the proximal promoter were predicted to alter the binding of the zinc finger protein Znf263 and long-distance interactions with other cis-regulatory regions. By surface plasmon resonance, we selected nucleotide substitutions in the 3’UTRs able to affect miRNA binding affinity. Specifically, miR-3150a-3p and miR-6785-5p expression inversely correlated with RUNX2 expression during in vitro osteogenic differentiation. The expression of two long non-coding RNAs, AL096865.1 and RUNX2-AS1, within the same locus, was modulated during in vitro osteogenic differentiation and correlated with the expression of specific RUNX2 isoforms. Our data suggest that RUNX2 may have undergone adaptive phenotypic evolution caused by epigenetic and post-transcriptional regulatory mechanisms, which may explain the delayed suture fusion leading to the present-day globular skull shape.

DNA Methylation Signature in Mononuclear Cells and Proinflammatory Cytokines May Define Molecular Subtypes in Sporadic Meniere Disease.

Meniere Disease (MD) is a multifactorial disorder of the inner ear characterized by vertigo attacks associated with sensorineural hearing loss and tinnitus with a significant heritability. Although MD has been associated with several genes, no epigenetic studies have been performed on MD. Here we performed whole-genome bisulfite sequencing in 14 MD patients and six healthy controls, with the aim of identifying an MD methylation signature and potential disease mechanisms. We observed a high number of differentially methylated CpGs (DMC) when comparing MD patients to controls (n= 9545), several of them in hearing loss genes, such as PCDH15, ADGRV1 and CDH23. Bioinformatic analyses of DMCs and cis-regulatory regions predicted phenotypes related to abnormal excitatory postsynaptic currents, abnormal NMDA-mediated receptor currents and abnormal glutamate-mediated receptor currents when comparing MD to controls. Moreover, we identified various DMCs in genes previously associated with cochleovestibular phenotypes in mice. We have also found 12 undermethylated regions (UMR) that were exclusive to MD, including two UMR in an inter CpG island in the PHB gene. We suggest that the DNA methylation signature allows distinguishing between MD patients and controls. The enrichment analysis confirms previous findings of a chronic inflammatory process underlying MD.

Transcription factor MEF2D is required for the maintenance of MLL-rearranged acute myeloid leukemia

AbstractAcute myeloid leukemia (AML) with MLL-rearrangement (MLL-r) comprises ∼10% of all AML cases and portends poor outcomes. Much remains uncovered on how MLL-r AML drives leukemia development while preventing cells from normal myeloid differentiation. Here, we identified that transcription factor MEF2D is a super-enhancer-associated, highly expressed gene in MLL-r AML. Knockout of MEF2D profoundly impaired leukemia growth, induced myeloid differentiation, and delayed oncogenic progression in vivo. Mechanistically, MEF2D loss led to robust activation of a CEBPE-centered myeloid differentiation program in AML cells. Chromatin profiling revealed that MEF2D binds to and suppresses the chromatin accessibility of CEBPE cis-regulatory regions. In human acute leukemia samples, MEF2D expression showed a strong negative correlation with the expression of CEBPE. Depletion of CEBPE partially rescued the cell growth defect and myeloid cell differentiation induced by the loss of MEF2D. Lastly, we show that MEF2D is positively regulated by HOXA9, and downregulation of MEF2D is an important mechanism for DOT1L inhibitor-induced antileukemia effects. Collectively, our findings suggest that MEF2D plays a critical role in human MLL-r AML and uncover the MEF2D-CEBPE axis as a crucial transcriptional mechanism regulating leukemia cell self-renewal and differentiation block.

Cis-regulatory elements of the cholinergic gene locus in the silkworm Bombyx mori.

Genes of choline acetyltransferase (ChAT) and vesicular acetylcholine transporter are encoded in the same gene locus, called the cholinergic gene locus. They are essential in cholinergic neurons to maintain their functional phenotype. The genomic structure of the cholinergic gene locus is conserved among invertebrates to mammals. However, the cholinergic gene expression in a specific subset of neurons is unknown in insects except for Drosophila melanogaster. In this study, we analysed the upstream sequence of cholinergic gene locus in the silkworm Bombyx mori to identify specific cis-regulatory regions. We found multiple enhancer regions that are localized within 1 kb upstream of the cholinergic gene locus. The combination of promoter assays using small deletions and bioinformatic analysis among insect species illuminates two conserved sequences in the cis-regulatory region: TGACGTA and CCAAT, which are known as the cAMP response element and CAAT box, respectively. We found that dibutyryl-cAMP, an analogue of cAMP, influences the expression of ChAT in B. mori. Tissue-specific expression analysis of transcriptional factors identified potential candidates that control the cholinergic gene locus expression. Our investigation provides new insight into the regulation mechanism of cholinergic neuron-specific gene machinery in this lepidopteran insect.

Differential regulation of mammalian and avian ATOH1 by E2F1 and its implication for hair cell regeneration in the inner ear

The mammalian inner ear has a limited capacity to regenerate its mechanosensory hair cells. This lack of regenerative capacity underlies the high incidence of age-related hearing loss in humans. In contrast, non-mammalian vertebrates can form new hair cells when damage occurs, a mechanism that depends on re-activation of expression of the pro-hair cell transcription factor Atoh1. Here, we show that members of the E2F transcription factor family, known to play a key role in cell cycle progression, regulate the expression of Atoh1. E2F1 activates chicken Atoh1 by directly interacting with a cis-regulatory region distal to the avian Atoh1 gene. E2F does not activate mouse Atoh1 gene expression, since this regulatory element is absent in mammals. We also show that E2F1 expression changes dynamically in the chicken auditory epithelium during ototoxic damage and hair cell regeneration. Therefore, we propose a model in which the mitotic regeneration of non-mammalian hair cells is due to E2F1-mediated activation of Atoh1 expression, a mechanism which has been lost in mammals.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

Transgenic reporter analysis of ChIP-Seq-defined enhancers identifies novel target genes for the terminal selector UNC-3/Collier/Ebf.

Terminal selector-type transcription factors are key regulators of neuronal identity and function (Hobert and Kratsios, 2019; Kratsios and Hobert, 2018). Mechanistically, terminal selectors are thought to act directly through binding at the cis-regulatory region of genes (termed “terminal identity genes”) that encode, among others, neurotransmitter [NT] synthesis proteins, ion channels, neuropeptides, and cell adhesion molecules (Hobert and Kratsios, 2019; Kratsios and Hobert, 2018). Although dozens of terminal selectors have been described thus far for individual neuron types of the nematode C. elegans (Hobert, 2016), the identification of their target genes has primarily relied on candidate approaches and availability of markers for neuronal terminal identity. Hence, unbiased methods are needed to identify the full spectrum of terminal selector target genes in individual neuron types. This study focuses on the phylogenetically conserved terminal selector UNC-3/Ebf (member of the Collier/Olf/Ebf family), which controls cholinergic motor neuron (MN) identity in the ventral nerve cord of the nematode C. elegans. To identify novel UNC-3 target genes, we took advantage of the genome-wide binding map of UNC-3 from our previous Chromatin Immunoprecipitation followed by Sequencing (ChIP-Seq) analysis (Li et al., 2020). We generated transgenic reporter lines for ten putative terminal identity genes (pxd-1, cal-2, lgc-4, ldb-1, nep-21, D2007.2, dmsr-2, ncs-2, npr-29, drn-1), whose expression patterns were largely unknown in C. elegans. Six of these reporter lines showed expression in ventral nerve cord MNs (nep-21, D2007.2, dmsr-2, ncs-2, npr-29, drn-1), whereas the remaining four (pxd-1, cal-2, lgc-4, ldb-1) showed expression in head and tail neurons, as well as some non-neuronal cells. Importantly, the number of ventral nerve cord MNs showing expression of the nep-21, D2007.2, and dmsr-2 reporters was significantly reduced in unc-3 null mutant animals, thereby expanding the repertoire of known UNC-3 target genes in these cells. Altogether, this study demonstrates that transgenic reporter analysis guided by ChIP-Seq results is a relatively efficient approach for the identification and validation of transcription factor target genes.

The \u03b22Tubulin, Rad50-ATPase and enolase cis-regulatory regions mediate male germline expression in Tribolium castaneum

Genetics-based pest management processes, including the sterile insect technique, are an effective method for the control of some pest insects. However, current SIT methods are not directly transferable to many important pest insect species due to the lack of genetic sexing strains. Genome editing is revolutionizing the way we conduct genetics in insects, including in Tribolium castaneum, an important genetic model and agricultural pest. We identified orthologues of β2Tubulin, Rad50-ATPase and enolase in T. castaneum. Using RT-PCR, we confirmed that these genes are predominantly expressed in the testis. PiggyBac-based transformation of T. castaneum cis-regulatory regions derived from Tc-β2t, Tc-rad50 or Tc-eno resulted in EGFP expression specifically in the T. castaneum testis. Additionally, we determined that each of these regulatory regions regulates EGFP expression in different cell types of the male gonad. Cis-regulatory regions from Tc-β2t produced EGFP expression throughout spermatogenesis and also in mature sperms; Tc-rad50 resulted in expression only in the haploid spermatid, while Tc-eno expressed EGFP in late spermatogenesis. In summary, the regulatory cis-regions characterized in this study are not only suited to study male gonadal function but could be used for development of transgenic sexing strains that produce one sex in pest control strategies.

Analysis of the proximal promoter of the human colon-specific B4GALNT2 (Sda synthase) gene: B4GALNT2 is transcriptionally regulated by ETS1

Background: The Sda antigen and corresponding biosynthetic enzyme B4GALNT2 are primarily expressed in normal colonic mucosa and are down-regulated to a variable degree in colon cancer tissues. Although their expression profile is well studied, little is known about the underlying regulatory mechanisms. Methods: To clarify the molecular basis of Sda expression in the human gastrointestinal tract, we investigated the transcriptional regulation of the human B4GALNT2 gene. The proximal promoter region was delineated using luciferase assays and essential trans-acting factors were identified through transient overexpression and silencing of several transcription factors. Results: A short cis-regulatory region restricted to the −72 to +12 area upstream of the B4GALNT2 short-type transcript variant contained the essential promoter activity that drives the expression of the human B4GALNT2 regardless of the cell type. We further showed that B4GALNT2 transcriptional activation mostly requires ETS1 and to a lesser extent SP1. Conclusions: Results presented herein are expected to provide clues to better understand B4GALNT2 regulatory mechanisms.

Abstract MP37: Epigenetic Changes Following Maternal Gut Microbiota Improvement With Dietary Fibre Lead To Cardio-protection In The Offspring

Dietary fibre is fermented by the gut microbiota and protects against the development of cardiovascular disease (CVD) through the production of gut microbial metabolites. We hypothesised dietary fibre intake during pregnancy may prevent the development of CVD in the offspring via in utero epigenetic mechanisms. To investigate this, we fed C57BL/6J female mice diets high or low in resistant starches (‘high-fibre’ and ‘low-fibre’, respectively) during gestation. At 6-weeks of age, we performed single-cell RNA-sequencing in the offspring (n=8/group) or they were challenged with saline (sham) or angiotensin II (Ang II, 0.25mg/kg/day, n=18-23/group). Maternal diet resulted in a distinct gut microbial composition ( P =0.001). This was still evident in the adult offspring, with high-fibre offspring having a different gut microbial colonisation ( P =0.001), irrespective of sham/Ang II treatment. Maternal fibre intake significantly changed the cardiac cellular and molecular landscape and promoted differential gene signatures in the offspring. This included upregulation of genes associated with extracellular matrix production in the offspring from low-fibre mothers. When challenged with Ang II, low-fibre offspring developed increased cardiac hypertrophy ( P =0.034) and fibrosis ( P =0.01) compared to high-fibre offspring. This was accompanied with decreased ejection fraction ( P =0.001) and increased left ventricular posterior wall thickness ( P= 0.017). These changes were independent of blood pressure. High-fibre offspring had decreased expression of natriuretic peptides ( Nppa , P =0.03, Nppb , P =0.002) compared to low-fibre offspring. Chromatin-immunoprecipitation assay revealed decrease in H3-acetylation at the cis-regulatory region of Nppa gene in Ang II-treated high-fibre offspring (P=0.002), suggesting that maternal fibre intake influences the epigenetic changes of the Nppa gene in the offspring’s heart. Together, these data reveal maternal high-fibre intake leads to foetal epigenetic reprogramming, likely through maternal to foetal transfer of gut microbial-derived metabolites.

Integration of transcription coregulator complexes with sequence-specific DNA-binding factor interactomes

The domain of transcription regulation has been notoriously difficult to annotate in the Gene Ontology, partly because of the intricacies of gene regulation which involve molecular interactions with DNA as well as amongst protein complexes. The molecular function ‘transcription coregulator activity’ is a part of the biological process ‘regulation of transcription, DNA-templated’ that occurs in the cellular component ‘chromatin’. It can mechanistically link sequence-specific DNA-binding transcription factor (dbTF) regulatory DNA target sites to coactivator and corepressor target sites through the molecular function ‘cis-regulatory region sequence-specific DNA binding’. Many questions arise about transcription coregulators (coTF). Here, we asked how many unannotated, putative coregulators can be identified in protein complexes? Therefore, we mined the CORUM and hu.MAP protein complex databases with known and strongly presumed human transcription coregulators. In addition, we trawled the BioGRID and IntAct molecular interaction databases for interactors of the known 1457 human dbTFs annotated by the GREEKC and GO consortia. This yielded 1093 putative transcription factor coregulator complex subunits, of which 954 interact directly with a dbTF. This substantially expands the set of coTFs that could be annotated to ‘transcription coregulator activity’ and sets the stage for renewed annotation and wet-lab research efforts. To this end, we devised a prioritisation score based on existing GO annotations of already curated transcription coregulators as well as interactome representation. Since all the proteins that we mined are parts of protein complexes, we propose to concomitantly engage in annotation of the putative transcription coregulator-containing complexes in the Complex Portal database.

K-mer Content Changes with Node Degree in Promoter-Enhancer Network of Mouse ES Cells.

Maps of Hi-C contacts between promoters and enhancers can be analyzed as networks, with cis-regulatory regions as nodes and their interactions as edges. We checked if in the published promoter–enhancer network of mouse embryonic stem (ES) cells the differences in the node type (promoter or enhancer) and the node degree (number of regions interacting with a given promoter or enhancer) are reflected by sequence composition or sequence similarity of the interacting nodes. We used counts of all k-mers (k = 4) to analyze the sequence composition and the Euclidean distance between the k-mer count vectors (k-mer distance) as the measure of sequence (dis)similarity. The results we obtained with 4-mers are interpretable in terms of dinucleotides. Promoters are GC-rich as compared to enhancers, which is known. Enhancers are enriched in scaffold/matrix attachment regions (S/MARs) patterns and depleted of CpGs. Furthermore, we show that promoters are more similar to their interacting enhancers than vice-versa. Most notably, in both promoters and enhancers, the GC content and the CpG count increase with the node degree. As a consequence, enhancers of higher node degree become more similar to promoters, whereas higher degree promoters become less similar to enhancers. We confirmed the key results also for human keratinocytes.

ChIP-GSM: Inferring active transcription factor modules to predict functional regulatory elements.

Transcription factors (TFs) often function as a module including both master factors and mediators binding at cis-regulatory regions to modulate nearby gene transcription. ChIP-seq profiling of multiple TFs makes it feasible to infer functional TF modules. However, when inferring TF modules based on co-localization of ChIP-seq peaks, often many weak binding events are missed, especially for mediators, resulting in incomplete identification of modules. To address this problem, we develop a ChIP-seq data-driven Gibbs Sampler to infer Modules (ChIP-GSM) using a Bayesian framework that integrates ChIP-seq profiles of multiple TFs. ChIP-GSM samples read counts of module TFs iteratively to estimate the binding potential of a module to each region and, across all regions, estimates the module abundance. Using inferred module-region probabilistic bindings as feature units, ChIP-GSM then employs logistic regression to predict active regulatory elements. Validation of ChIP-GSM predicted regulatory regions on multiple independent datasets sharing the same context confirms the advantage of using TF modules for predicting regulatory activity. In a case study of K562 cells, we demonstrate that the ChIP-GSM inferred modules form as groups, activate gene expression at different time points, and mediate diverse functional cellular processes. Hence, ChIP-GSM infers biologically meaningful TF modules and improves the prediction accuracy of regulatory region activities.

Cortex cis-regulatory switches establish scale colour identity and pattern diversity in Heliconius.

In Heliconius butterflies, wing colour pattern diversity and scale types are controlled by a few genes of large effect that regulate colour pattern switches between morphs and species across a large mimetic radiation. One of these genes, cortex, has been repeatedly associated with colour pattern evolution in butterflies. Here we carried out CRISPR knockouts in multiple Heliconius species and show that cortex is a major determinant of scale cell identity. Chromatin accessibility profiling and introgression scans identified cis-regulatory regions associated with discrete phenotypic switches. CRISPR perturbation of these regions in black hindwing genotypes recreated a yellow bar, revealing their spatially limited activity. In the H. melpomene/timareta lineage, the candidate CRE from yellow-barred phenotype morphs is interrupted by a transposable element, suggesting that cis-regulatory structural variation underlies these mimetic adaptations. Our work shows that cortex functionally controls scale colour fate and that its cis-regulatory regions control a phenotypic switch in a modular and pattern-specific fashion.

The Novel Non-coding Transcriptional Regulator Gm18840 Drives Cardiomyocyte Apoptosis in Myocardial Infarction Post Ischemia/Reperfusion.

BackgroundIschemia/reperfusion-mediated myocardial infarction (MIRI) is a major pathological factor implicated in the progression of ischemic heart disease, but the key factors dysregulated during MIRI have not been fully elucidated, especially those essential non-coding factors required for cardiovascular development.MethodsA murine MIRI model and RNA sequencing (RNA-seq) were used to identify key lncRNAs after myocardial infarction. qRT-PCR was used to validate expression in cardiac muscle tissues and myocardial cells. The role of Gm18840 in HL-1 cell growth was determined by flow cytometry experiments in vitro. Full-length Gm18840 was identified by using a rapid amplification of cDNA ends (RACE) assay. The subcellular distribution of Gm18840 was examined by nuclear/cytoplasmic RNA fractionation and qRT-PCR. RNA pulldown and RNA immunoprecipitation (RIP)-qPCR assays were performed to identify Gm18840-interacting proteins. Chromatin isolation by RNA purification (ChIRP)-seq (chromatin isolation by RNA purification) was used to identify the genome-wide binding of Gm18840 to chromatin. The regulatory activity of Gm18840 in transcriptional regulation was examined by a luciferase reporter assay and qRT-PCR.ResultsGm18840 was upregulated after myocardial infarction in both in vivo and in vitro MIRI models. Gm18840 was 1,471 nt in length and localized in both the cytoplasm and the nucleus of HL-1 cells. Functional studies showed that the knockdown of Gm18840 promoted the apoptosis of HL-1 cells. Gm18840 directly interacts with histones, including H2B, highlighting a potential function in transcriptional regulation. Further ChIRP-seq and RNA-seq analyses showed that Gm18840 is directly bound to the cis-regulatory regions of genes involved in developmental processes, such as Junb, Rras2, and Bcl3.ConclusionGm18840, a novel transcriptional regulator, promoted the apoptosis of myocardial cells via direct transcriptional regulation of essential genes and might serve as a novel therapeutic target for MIRI.

UniBind: maps of high-confidence direct TF-DNA interactions across nine species

BackgroundTranscription factors (TFs) bind specifically to TF binding sites (TFBSs) at cis-regulatory regions to control transcription. It is critical to locate these TF-DNA interactions to understand transcriptional regulation. Efforts to predict bona fide TFBSs benefit from the availability of experimental data mapping DNA binding regions of TFs (chromatin immunoprecipitation followed by sequencing - ChIP-seq).ResultsIn this study, we processed ~ 10,000 public ChIP-seq datasets from nine species to provide high-quality TFBS predictions. After quality control, it culminated with the prediction of ~ 56 million TFBSs with experimental and computational support for direct TF-DNA interactions for 644 TFs in > 1000 cell lines and tissues. These TFBSs were used to predict > 197,000 cis-regulatory modules representing clusters of binding events in the corresponding genomes. The high-quality of the TFBSs was reinforced by their evolutionary conservation, enrichment at active cis-regulatory regions, and capacity to predict combinatorial binding of TFs. Further, we confirmed that the cell type and tissue specificity of enhancer activity was correlated with the number of TFs with binding sites predicted in these regions. All the data is provided to the community through the UniBind database that can be accessed through its web-interface (https://unibind.uio.no/), a dedicated RESTful API, and as genomic tracks. Finally, we provide an enrichment tool, available as a web-service and an R package, for users to find TFs with enriched TFBSs in a set of provided genomic regions.ConclusionsUniBind is the first resource of its kind, providing the largest collection of high-confidence direct TF-DNA interactions in nine species.