Transcriptional Regulatory Elements Research Articles

Study of hub nodes of transcription factor-target gene regulatory network and immune mechanism for type 2 diabetes based on chip analysis of GEO database

Identification of novel therapeutic targets for type 2 diabetes is a key area of contemporary research. In this study, we screened differentially expressed genes in type 2 diabetes through the GEO database and sought to identify the key virulence factors for type 2 diabetes through a transcription factor regulatory network. Our findings may help identify new therapeutic targets for type 2 diabetes. Data pertaining to the humoral (whole blood) gene expression profile of diabetic patients were obtained from the NCBI’s GEO Datasets database and gene sets with differential expression were identified. Subsequently, the TRED transcriptional regulatory element database was integrated to build a gene regulatory network for type 2 diabetes. Functional analysis (GO-Analysis) and Pathway-analysis of differentially expressed genes were performed using the DAVID database and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Finally, gene-disease correlation analysis was performed using the DAVID online annotation tool. A total of 236 pathogenic genes, four transcription factors related to the pathogenic genes, and 261 corresponding target genes were identified. A transcription factor-target gene regulatory network for type 2 diabetes was constructed. Most of the key factors of the transcription factor-target gene regulatory network for type 2 diabetes were found closely related to the immune metabolic system and the functions of cell proliferation and transformation.

ChIP-Atlas 3.0: a data-mining suite to explore chromosome architecture together with large-scale regulome data.

ChIP-Atlas (https://chip-atlas.org/) presents a suite of data-mining tools for analyzing epigenomic landscapes, powered by the comprehensive integration of over 376000 public ChIP-seq, ATAC-seq, DNase-seqand Bisulfite-seq experiments from six representative model organisms. To unravel the intricacies of chromatin architecture that mediates the regulome-initiated generation of transcriptional and phenotypic diversity within cells, we report ChIP-Atlas 3.0 that enhances clarity by incorporating additional tracks for genomic and epigenomic features within a newly consolidated 'annotation track' section. The tracks include chromosomal conformation (Hi-C and eQTL datasets), transcriptional regulatory elements (ChromHMM and FANTOM5 enhancers), and genomic variants associated with diseases and phenotypes (GWAS SNPs and ClinVar variants). These annotation tracks are easily accessible alongside other experimental tracks, facilitating better elucidation of chromatin architecture underlying the diversification of transcriptional and phenotypic traits. Furthermore, 'Diff Analysis,' a new online tool, compares the query epigenome data to identify differentially bound, accessible, and methylated regions using ChIP-seq, ATAC-seq and DNase-seq, and Bisulfite-seq datasets, respectively. The integration of annotation tracks and the Diff Analysis tool, coupled with continuous data expansion, renders ChIP-Atlas 3.0 a robust resource for mining the landscape of transcriptional regulatory mechanisms, thereby offering valuable perspectives, particularly for genetic disease research and drug discovery.

Retrotransposons in embryogenesis and neurodevelopment.

Retrotransposable elements (RTEs) are genetic elements that can replicate and insert new copies into different genomic locations. RTEs have long been identified as 'parasitic genes', as their mobilization can cause mutations, DNA damage, and inflammation. Interestingly, high levels of retrotransposon activation are observed in early embryogenesis and neurodevelopment, suggesting that RTEs may possess functional roles during these stages of development. Recent studies demonstrate that RTEs can function as transcriptional regulatory elements through mechanisms such as chromatin organization and noncoding RNAs. It is clear, however, that RTE expression and activity must be restrained at some level during development, since overactivation of RTEs during neurodevelopment is associated with several developmental disorders. Further investigation is needed to understand the importance of RTE expression and activity during neurodevelopment and the balance between RTE-regulated development and RTE-mediated pathogenesis.

The choreography of chromatin in RNA polymerase III regulation.

Regulation of eukaryotic gene expression involves a dynamic interplay between the core transcriptional machinery, transcription factors, and chromatin organization and modification. While this applies to transcription by all RNA polymerase complexes, RNA polymerase III (RNAPIII) seems to be atypical with respect to its mechanisms of regulation. One distinctive feature of most RNAPIII transcribed genes is that they are devoid of nucleosomes, which relates to the high levels of transcription. Moreover, most of the regulatory sequences are not outside but within the transcribed open chromatin regions. Yet, several lines of evidence suggest that chromatin factors affect RNAPIII dynamics and activity and that gene sequence alone does not explain the observed regulation of RNAPIII. Here we discuss the role of chromatin modification and organization of RNAPIII transcribed genes and how they interact with the core transcriptional RNAPIII machinery and regulatory DNA elements in and around the transcribed genes.

Avian leukosis virus usurps the cellular SERBP1 protein to enhance its transcription and promote productive infections in avian cells

Avian leukosis virus subgroup K (ALV-K) is composed of newly emerging isolates, which cluster separately from the well-characterized subgroups A, B, C, D, E, and J in sequence analysis, and exhibits a specific host range and a unique pattern of superinfection interference. Avian leukosis virus subgroup K replicate more slowly in avian cells than other ALV strains, leading to escaped detection during ALV eradication, but the underlying mechanism are largely unknown. In our previous study, we have reported that JS11C1 and most of other suspected ALV-K strains possessed unique mutations in the U3 region. Here, we selected 5 mutations in some important transcriptional regulation elements to explore the possible factor contributing for the lower activity of LTR, including CA-TG mutation in the CAAT box, 21 nt deletion in the CAAT box, A-G and A-T mutations in the CArG boxes, 11 nt insertion in the PRE boxes, and C-T mutation in the TATA box. On the basis of infectious clone of JS11C1, we demonstrated that the 11 nt fragment in the PRE boxes was associated with the transcription activity of LTR, the enhancer ability of U3, and the replication capacity of the virus. Notably, we determined the differential U3-protein interaction profile of ALVs and found that the 11 nt fragment specifically binds to cellular SERPINE1 mRNA binding protein 1 (SERBP1) to increase the LTR activity and enhance virus replication. Collectively, these findings reveal that a 11 nt fragment in the U3 gene contributed to its binding ability to the cellular SERBP1 to enhance its transcription and the infectious virus productions in avian cells. This study highlighted the vital role of host factor in retrovirus replication and thus provides a new perspective to elucidate the interaction between retrovirus and its host and a molecular basis to develop efficient strategies against retroviruses.

GEO dataset mining analysis reveals novel Staphylococcus aureus virulence gene regulatory networks and diagnostic targets in mice.

Staphylococcus (S.) aureus infection is a serious, worldwide health concern, particularly in many communities and hospitals. Understanding the S. aureus pathogenetic regulatory network will provide significant insights into diagnostic target screening to improve clinical treatment of diseases caused by S. aureus. We screened differentially expressed genes between normal mice and S. aureus-infected mice. We used the Gene Expression Omnibus (GEO) DataSets database for functional analysis (GO-analysis) and the DAVID and KEGG databases for signaling pathway analyses. We next integrated the gene and pathway analyses with Transcriptional Regulatory Element Database (TRED) to build an antimicrobial resistance gene regulatory network of S. aureus. We performed association analysis of network genes and diseases using DAVID online annotation tools. We identified a total of 437 virulence genes and 15 transcription factors (TFs), as well as 444 corresponding target genes, in the S. aureus TF regulatory network. We screened seven key network nodes (Met, Mmp13, Il12b, Il4, Tnf, Ptgs2, and Ctsl), four key transcription factors (Jun, C3, Spil, and Il6) and an important signaling pathway (TNF). We hypothesized that the cytokine activity and growth factor activity of S. aureus are combinatorically cross-regulated by Met, Mmp13, Il12b, Il4, Tnf, Ptgs2, and Ctsl genes, the TFs Jun, C3, Spi1, and Il6, as well as the immune response, cellular response to lipopolysaccharide, and inflammatory response. Our study provides information and reference values for the molecular understanding of the S. aureus pathogenetic gene regulatory network.

Enhancement of PRMT6 binding to a novel germline GATA1 mutation associated with congenital anemia.

Mutations in the master hematopoietic transcription factor GATA1 are often associated with functional defects in erythropoiesis and megakaryopoiesis. In this study, we identified a novel GATA1 germline mutation (c.1162delGG, p.Leu387Leufs*62) in a patient with congenital anemia and occasional thrombocytopenia. The C-terminal GATA1, a rarely studied mutational region, undergoes frameshifting translation as a consequence of this double-base deletion mutation. To investigate the specific function and pathogenic mechanism of this mutant, in vitro mutant models of stable re-expression cells were generated. The mutation was subsequently validated to cause diminished transcriptional activity of GATA1 and defective differentiation of erythroid and megakaryocytes. Using proximity labeling and mass spectrometry, we identified selective alterations in the proximal protein networks of the mutant, revealing decreased binding to a set of normal GATA1-interaction proteins, including the essential co-factor FOG1. Notably, our findings further demonstrated enhanced recruitment of the protein arginine methyltransferase PRMT6, which mediates histone modification at H3R2me2a and represses transcription activity. We also found an enhanced binding of this mutant GATA1/PRMT6 complex to the transcriptional regulatory elements of GATA1's target genes. Moreover, treatment of the PRMT6 inhibitor MS023 could partially rescue the inhibited transcriptional and impaired erythroid differentiation caused by the GATA1 mutation. Taken together, our results provide molecular insights into erythropoiesis in which mutation leads to partial loss of GATA1 function and the broader role of PRMT6 and its inhibitor MS023 in congenital anemia, highlighting PRMT6 binding as a negative factor of GATA1 transcriptional activity in aberrant hematopoiesis.

Abstract B087: Mutant-p53 amplifies Cxcl1 expression from distal enhancers blunting immune checkpoint inhibition efficacy in pancreatic cancer

Abstract Pancreatic ductal Adenocarcinoma (PDAC) is an aggressive malignancy complicated by poor early diagnosis and a lack of response to traditional treatments. It is characterized by a desmoplastic stroma, a lack of infiltration and activation of T cells, and a low mutational burden. The genetic landscape of PDAC is defined by activating KRAS mutations (~90%) and p53 alterations (~70%), but the molecular switches perturbed by these genetic aberrations remain unclear. p53 missense mutations, unlike mutations resulting in the loss of p53, are considered to acquire tumor-supporting functions. But these novel functions remain uncharacterized. The ambiguity in the molecular mechanism of mutant-p53 is further exacerbated by the existence of various types of p53 mutations. The majority of p53 mutations are missense mutations in the DNA binding domain. The repertoire of transcription factors (TFs) it can interact with and the vast regulatory landscape of each TF—composed of gene promoters and distal enhancers—present obstacles in understanding the molecular mechanisms promoting PDAC. In this study, we examined how a common p53 missense mutation in PDAC plays a role in weakening the Immune checkpoint Inhibitors (ICIs) efficacy. Using cells derived from a genetically engineered mouse model (GEMM) of PDAC with activating KRAS mutation (KrasG12D/+) and a p53 missense mutation (p53R172H/-), we found that the PDAC tumorigenesis and resistance to ICIs are dependent on the mutant-p53. We used isogenic p53-null PDAC cells and the restoration of p53R172H in p53-null cells to demonstrate the role of p53R172H in controlling the expression of immunosuppressive chemokine genes such as Cxcl1. p53R172H deletion attenuated PDAC tumor growth, increased the influx of cytotoxic T-cells, and sensitized the tumor to ICIs. The p53R172H-mediated TME reprogramming was replicated by the deletion of the Cxcl1 gene, suggesting the anti-tumorigenic effect of p53R172H was mediated by the Cxcl1 gene. We probed the mechanism of Cxcl1 expression dependence on p53R172H. We found that in conjunction with NF-kB, p53R172H occupies the distal transcription regulatory elements (dTREs) of the Cxcl1 gene harboring NF-kB binding sites. Strikingly, deletion of the Cxcl1 dTREs in PDAC cells recapitulates the phenotypes of p53R172H deletion and Cxcl1 deletion in terms of tumor size, immune landscape of the TME, and ICI responsiveness. Furthermore, we examined the interplay between p53R172H and NF-kB and found that the p53R172H physically interacts with the NF-kB subunit RelA and facilitates its nuclear translocation. Overall, we characterize how a common p53 mutation in PDAC co-opts non-coding regulatory DNA to augment the expression of selective chemokine genes and establishes an immunosuppressive TME to shield the therapeutic benefits of ICIs. Citation Format: Dig B. Mahat, Heena Kumra, Emily Metcalf, Sarah Castro, Kim Nguyen1, Arundeep Singh, William W. Ho, Ivy Chen, Brandon Sullivan, Leon Yim, Enrico Moiso, Vikash Chauhan, Hernandez Moura Silva, Stefani Spranger, Rakesh Jain, Phillip A. Sharp. Mutant-p53 amplifies Cxcl1 expression from distal enhancers blunting immune checkpoint inhibition efficacy in pancreatic cancer [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Pancreatic Cancer; 2023 Sep 27-30; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(2 Suppl):Abstract nr B087.

Abstract PR009: SMARCAL1 is a dual regulator of innate immune signaling and PD-L1 expression that promotes tumor immune evasion

Abstract Genomic instability can trigger cancer-intrinsic innate immune responses that promote tumor rejection. However, cancer cells often evade these responses by overexpressing immune checkpoint regulators, such as PD-L1. Here, we identify the SNF2-family DNA translocase SMARCAL1 as a factor that favors tumor immune evasion by a dual mechanism involving both the suppression of innate immune signaling and the induction of PD-L1-mediated immune checkpoint responses. Mechanistically, SMARCAL1 relieves endogenous DNA damage and suppresses cGAS-STING-dependent immune signaling during cancer cell growth. Simultaneously, it cooperates with the AP-1 family member JUN to maintain chromatin accessibility at a transcriptional regulatory element in the PD-L1 gene, thereby promoting PD-L1 expression in cancer cells. Loss of SMARCAL1 enhances anti-tumor immune responses and sensitizes tumors to immune checkpoint blockade in a mouse melanoma model. Collectively, these studies uncover SMARCAL1 as a valuable target for cancer immunotherapy. Citation Format: Giuseppe Leuzzi, Alessandro Vasciaveo, Angelo Taglialatela, Xiao Chen, Tessa M. Firestone, Allison R. Hickman, Wendy Mao, Tanay Thakar, Alina Vaitsiankova, Jen-Wei Huang, Raquel Cuella-Martin, Samuel B. Hayward, Jordan S. Kesner, Ali Ghasemzadeh, Tarun S. Nambiar, Patricia Ho, Alexander Rialdi, Maxime Hebrard, Yinglu Li, Jinmei Gao, Saarang Gopinath, Oluwatobi A. Adeleke, Bryan Venters, Charles G. Drake, Richard Baer, Benjamin Izar, Ernesto Guccione, Michael-Christopher Keogh, Raphael Guerois, Lu Sun, Chao Lu, Andrea Califano, Alberto Ciccia. SMARCAL1 is a dual regulator of innate immune signaling and PD-L1 expression that promotes tumor immune evasion [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Damage Repair: From Basic Science to Future Clinical Application; 2024 Jan 9-11; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2024;84(1 Suppl):Abstract nr PR009.

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm.

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

LPS and type I and II interferons have opposing effects on epigenetic regulation of LAIR1 expression in mouse and human macrophages.

Inhibitory immune receptors are important for maintaining immune homeostasis. We identified epigenetic alterations in 2 members of this group, LAIR1 and LAIR2, in lymphoma patients with inflammatory tissue damage and susceptibility to infection. We predicted that the expression of LAIR genes is controlled by immune mediators acting on transcriptional regulatory elements. Using flow cytometry, quantitative reverse-transcription polymerase chain reaction, and RNA sequencing, we measured LAIR1 and LAIR2 in human and murine immune cell subsets at baseline and posttreatment with immune mediators, including type I and II interferons, tumor necrosis factor α, and lipopolysaccharide (LPS). We identified candidate regulatory elements using epigenome profiling and measured their regulatory activity using luciferase reporters. LAIR1 expression substantially increases during monocyte differentiation to macrophages in both species. In contrast, murine and human macrophages exhibited opposite changes in LAIR1 in response to immune stimuli: human LAIR1 increased with LPS while mouse LAIR1 increased with interferon γ. LAIR genes had distinct patterns of enhancer activity with variable responses to immune stimuli. To identify relevant transcription factors (TFs), we developed integrative bioinformatic techniques applied to TF chromatin immunoprecipitation sequencing, RNA sequencing, and luciferase activity, revealing distinct sets of TFs for each LAIR gene. Most strikingly, LAIR1 TFs include nuclear factor kappa B factors RELA and RELB, while Lair1 and LAIR2 instead include STAT3 and/or STAT5. Regulation by nuclear factor kappa B factors may therefore explain the LPS-induced increase in LAIR1 expression, in contrast to Lair1 decrease. Our findings reveal new insights into transcriptional mechanisms that control distinct expression patterns of LAIR genes in response to inflammatory stimuli in human and murine myeloid and lymphoid cells.

TE-TSS: an integrated data resource of human and mouse transposable element (TE)-derived transcription start site (TSS).

Transposable elements (TEs) are abundant in the genome and serve as crucial regulatory elements. Some TEs function as epigenetically regulated promoters, and these TE-derived transcription start sites (TSSs) play a crucial role in regulating genes associated with specific functions, such as cancer and embryogenesis. However, the lack of an accessible database that systematically gathers TE-derived TSS data is a current research gap. To address this, we established TE-TSS, an integrated data resource of human and mouse TE-derived TSSs (http://xozhanglab.com/TETSS). TE-TSS has compiled 2681 RNA sequencing datasets, spanning various tissues, cell linesand developmental stages. From these, we identified 5768 human TE-derived TSSs and 2797 mouse TE-derived TSSs, with 47% and 38% being experimentally validated, respectively. TE-TSS enables comprehensive exploration of TSS usage in diverse samples, providing insights into tissue-specific gene expression patterns and transcriptional regulatory elements. Furthermore, TE-TSS compares TE-derived TSS regions across 15 mammalian species, enhancing our understanding of their evolutionary and functional aspects. The establishment of TE-TSS facilitates further investigations into the roles of TEs in shaping the transcriptomic landscape and offers valuable resources for comprehending their involvement in diverse biological processes.

ERNA co-expression network uncovers TF dependency and convergent cooperativity

Enhancer RNAs (eRNAs) are non-coding RNAs produced by transcriptional enhancers that are highly correlated with their activity. Using a capped nascent RNA sequencing (PRO-cap) dataset in human lymphoblastoid cell lines across 67 individuals, we identified inter-individual variation in the expression of over 80 thousand transcribed transcriptional regulatory elements (tTREs), in both enhancers and promoters. Co-expression analysis of eRNAs from tTREs across individuals revealed how enhancers are associated with each other and with promoters. Mid- to long-range co-expression showed a distance-dependent decay that was modified by TF occupancy. In particular, we found a class of “bivalent” TFs, including Cohesin, that both facilitate and isolate the interaction between enhancers and/or promoters, depending on their topology. At short distances, we observed strand-specific correlations between nearby eRNAs in both convergent and divergent orientations. Our results support a cooperative model of convergent eRNAs, consistent with eRNAs facilitating adjacent enhancers rather than interfering with each other. Therefore, our approach to infer functional interactions from co-expression analyses provided novel insights into the principles of enhancer interactions as a function of distance, orientation, and binding landscapes of TFs.

Genome-wide profiling of transcription factor activity in primary liver cancer using single-cell ATAC sequencing

Primary liver cancer (PLC) consists of two main histological subtypes; hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA). The role of transcription factors (TFs) in malignant hepatobiliary lineage commitment between HCC and iCCA remains underexplored. Here, we present genome-wide profiling of transcription regulatory elements of 16 PLC patients using single-cell assay for transposase accessible chromatin sequencing. Single-cell open chromatin profiles reflect the compositional diversity of liver cancer, identifying both malignant and microenvironment component cells. TF motif enrichment levels of 31 TFs strongly discriminate HCC from iCCA tumors. These TFs are members of the nuclear/retinoid receptor, POU, or ETS motif families. POU factors are associated with prognostic features in iCCA. Overall, nuclear receptors, ETS and POU TF motif families delineate transcription regulation between HCC and iCCA tumors, which may be relevant to development and selection of PLC subtype-specific therapeutics.

Genome-wide DNase I-hypersensitive site assay reveals distinct genomic distributions and functional features of open chromatin in autopolyploid sugarcane.

The characterization of cis-regulatory DNA elements (CREs) is essential for deciphering the regulation of gene expression in eukaryotes. Although there have been endeavors to identify CREs in plants, the properties of CREs in polyploid genomes are still largely unknown. Here, we conducted the genome-wide identification of DNase I-hypersensitive sites (DHSs) in leaf and stem tissues of the auto-octoploid species Saccharum officinarum. We revealed that DHSs showed highly similar distributions in the genomes of these two S. officinarum tissues. Notably, we observed that approximately 74% of DHSs were located in distal intergenic regions, suggesting considerable differences in the abundance of distal CREs between S. officinarum and other plants. Leaf- and stem-dependent transcriptional regulatory networks were also developed by mining the binding motifs of transcription factors (TFs) from tissue-specific DHSs. Four TEOSINTE BRANCHED 1, CYCLOIDEA, and PCF1 (TCP) TFs (TCP2, TCP4, TCP7, and TCP14) and two ethylene-responsive factors (ERFs) (ERF109 and ERF03) showed strong causal connections with short binding distances from each other, pointing to their possible roles in the regulatory networks of leaf and stem development. Through functional validation in transiently transgenic protoplasts, we isolate a set of tissue-specific promoters. Overall, the DHS maps presented here offer a global view of the potential transcriptional regulatory elements in polyploid sugarcane and can be expected to serve as a valuable resource for both transcriptional network elucidation and genome editing in sugarcane breeding.

Human-specific evolutionary markers linked to foetal neurodevelopment modulate brain surface area in schizophrenia

Schizophrenia may represent a trade-off in the evolution of human-specific ontogenetic mechanisms that guide neurodevelopment. Human Accelerated Regions (HARs) are evolutionary markers functioning as neurodevelopmental transcription enhancers that have been associated with brain configuration, neural information processing, and schizophrenia risk. Here, we have investigated the influence of HARs’ polygenic load on neuroanatomical measures through a case-control approach (128 patients with schizophrenia and 115 controls). To this end, we have calculated the global schizophrenia Polygenic Risk Score (Global PRSSZ) and that specific to HARs (HARs PRSSZ). We have also estimated the polygenic burden restricted to the HARs linked to transcriptional regulatory elements active in the foetal brain (FB-HARs PRSSZ) and the adult brain (AB-HARs PRSSZ). We have explored the main effects of the PRSs and the PRSs x diagnosis interactions on brain regional cortical thickness (CT) and surface area (SA). The results indicate that a higher FB-HARs PRSSZ is associated with patients’ lower SA in the lateral orbitofrontal cortex, the superior temporal cortex, the pars triangularis and the paracentral lobule. While noHARs-derived PRSs show an effect on the risk, our neuroanatomical findings suggest that the human-specific transcriptional regulation during the prenatal period underlies SA variability, highlighting the role of these evolutionary markers in the schizophrenia genomic architecture.

Epigenomic analysis of formalin-fixed paraffin-embedded samples by CUT&Tag

For more than a century, formalin-fixed paraffin-embedded (FFPE) sample preparation has been the preferred method for long-term preservation of biological material. However, the use of FFPE samples for epigenomic studies has been difficult because of chromatin damage from long exposure to high concentrations of formaldehyde. Previously, we introduced Cleavage Under Targeted Accessible Chromatin (CUTAC), an antibody-targeted chromatin accessibility mapping protocol based on CUT&Tag. Here we show that simple modifications of our CUTAC protocol either in single tubes or directly on slides produce high-resolution maps of paused RNA Polymerase II at enhancers and promoters using FFPE samples. We find that transcriptional regulatory element differences produced by FFPE-CUTAC distinguish between mouse brain tumors and identify and map regulatory element markers with high confidence and precision, including microRNAs not detectable by RNA-seq. Our simple workflows make possible affordable epigenomic profiling of archived biological samples for biomarker identification, clinical applications and retrospective studies.

TMEM127 suppresses tumor development by promoting RET ubiquitination, positioning, and degradation

The TMEM127 gene encodes a transmembrane protein of poorly known function that is mutated in pheochromocytomas, neural crest-derived tumors of adrenomedullary cells. Here, we report that, at single-nucleus resolution, TMEM127-mutant tumors share precursor cells and transcription regulatory elements with pheochromocytomas carrying mutations of the tyrosine kinase receptor RET. Additionally, TMEM127-mutant pheochromocytomas, human cells, and mouse knockout models of TMEM127 accumulate RET and increase its signaling. TMEM127 contributes to RET cellular positioning, trafficking, and lysosome-mediated degradation. Mechanistically, TMEM127 binds to RET and recruits the NEDD4 E3 ubiquitin ligase for RET ubiquitination and degradation via TMEM127 C-terminal PxxY motifs. Lastly, increased cell proliferation and tumor burden after TMEM127 loss can be reversed by selective RET inhibitors invitro and invivo. Our results define TMEM127 as a component of the ubiquitin system and identify aberrant RET stabilization as a likely mechanism through which TMEM127 loss-of-function mutations cause pheochromocytoma.

Widespread priming of transcriptional regulatory elements by incipient accessibility or RNA polymerase II pause in early embryos of the sea urchin Strongylocentrotus purpuratus.

Transcriptional regulatory elements (TREs) are the primary nodes that control developmental gene regulatory networks. In embryo stages, larvae, and adult differentiated red spherule cells of the sea urchin Strongylocentrotus purpuratus, transcriptionally engaged TREs are detected by Precision Run-On Sequencing (PRO-seq), which maps genome-wide at base pair resolution the location of paused or elongating RNA polymerase II (Pol II). In parallel, TRE accessibility is estimated by the Assay for Transposase-Accessible Chromatin using Sequencing (ATAC-seq). Our analysis identifies surprisingly early and widespread TRE accessibility in 4-cell cleavage embryos that is not necessarily followed by concurrent or subsequent transcription. TRE transcriptional differences identified by PRO-seq provide more contrast among embryonic stages than ATAC-seq accessibility differences, in agreement with the apparent excess of accessible but inactive TREs during embryogenesis. Global TRE accessibility reaches a maximum around the 20-hour late blastula stage, which coincides with the consolidation of major embryo regionalizations and peak histone variant H2A.Z expression. A transcriptional potency model based on labile nucleosome TRE occupancy driven by DNA sequences and the prevalence of histone variants is proposed in order to explain the basal accessibility of transcriptionally inactive TREs during embryogenesis. However, our results would not reconcile well with labile nucleosome models based on simple A/T sequence enrichment. In addition, a large number of distal TREs become transcriptionally disengaged during developmental progression, in support of an early Pol II paused model for developmental gene regulation that eventually resolves in transcriptional activation or silencing. Thus, developmental potency in early embryos may be facilitated by incipient accessibility and transcriptional pause at TREs.

Chromatin accessibility memory of donor cells disrupts bovine somatic cell nuclear transfer blastocysts development.

The post-transfer developmental capacity of bovine somatic cell nuclear transfer (SCNT) blastocysts is reduced, implying that abnormalities in gene expression regulation are present at blastocyst stage. Chromatin accessibility, as an indicator for transcriptional regulatory elements mediating gene transcription activity, has heretofore been largely unexplored in SCNT embryos, especially at blastocyst stage. In the present study, single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) of in vivo and SCNT blastocysts were conducted to segregate lineages and demonstrate the aberrant chromatin accessibility of transcription factors (TFs) related to inner cell mass (ICM) development in SCNT blastocysts. Pseudotime analysis of lineage segregation further reflected dysregulated chromatin accessibility dynamics of TFs in the ICM of SCNT blastocysts compared to their in vivo counterparts. ATAC- and ChIP-seq results of SCNT donor cells revealed that the aberrant chromatin accessibility in the ICM of SCNT blastocysts was due to the persistence of chromatin accessibility memory at corresponding loci in the donor cells, with strong enrichment of trimethylation of histone H3 at lysine 4 (H3K4me3) at these loci. Correction of the aberrant chromatin accessibility through demethylation of H3K4me3 by KDM5B diminished the expression of related genes (e.g., BCL11B) and significantly improved the ICM proliferation in SCNT blastocysts. This effect was confirmed by knocking down BCL11B in SCNT embryos to down-regulate p21 and alleviate the inhibition of ICM proliferation. These findings expand our understanding of the chromatin accessibility abnormalities in SCNT blastocysts and BCL11B may be a potential target to improve SCNT efficiency.