ChIP-seq Assays Research Articles

Functional characteristics and computational model of abundant hyperactive loci in the human genome.

Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

Functional characteristics and computational model of abundant hyperactive loci in the human genome

Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

Comparative genomics provides insights into adaptive evolution and demographics of bats.

Bats possess a range of distinctive characteristics, including flight, echolocation, impressive longevity, and the ability to harbor various zoonotic pathogens. Additionally, they account for the second-highest species diversity among mammalian orders, yet their phylogenetic relationships and demographic history remain underexplored. Here, we generated de novo assembled genomes for 17 bat species and two of their mammalian relatives (the Amur hedgehog and Chinese mole shrew), with 12 genomes reaching chromosome-level assembly. Comparative genomics and ChIP-seq assays identified newly gained genomic regions in bats potentially linked to the regulation of gene activity and expression. Notably, some antiviral infection related gene under positive selection exhibited the activity of suppressing cancer, evidencing the linkage between virus tolerance and cancer resistance in bats. By integrating published bat genome assemblies, phylogenetic reconstruction established the proximity of noctilionoid bats to vesper bats. Interestingly, we found two distinct patterns of ancient population dynamics in bats and population changes since the last-glacial maximum do not reflect species phylogenetic relationships. These findings enriched our understanding of adaptive mechanisms and demographic history of bats.

Functional characteristics and computational model of abundant hyperactive loci in the human genome.

Enhancers and promoters are classically considered to be bound by a small set of TFs in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1,003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected 5 distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

The circZEB1/miR-337–3p/OGT axis mediates angiogenesis and metastasis via O-GlcNAcylation and up-regulating YBX1 in breast cancer

BackgroundA growing corpus of research has revealed that circular RNAs (circRNAs) have become increasingly important for the growth of malignancies in recent years. CircRNAs as ideal candidates for breast cancer (BC) therapeutic targets is still absent. MethodsIn our study, the dysregulated circRNAs in BC progression were explored, we analysed the BC's circRNA expression profiles using publicly available datasets (GSE101124 and GSE101122). The expression of circZEB1 in BC and cell lines was investigated by qPCR. RNase and actinomycin D were used to examine the features of circZEB1. The function of circZEB1 was subsequently investigated through the utilisation of colony formation, tube formation, transwell assays, and xenograft animal models.RNA immunoprecipitation (RIP), luciferase reporter assays, immunoprecipitation (co-IP) test in conjunction with LC-MS, and ChIP-seq assay to investigate the molecular mechanism underlying the biological activity of circZEB1 in BC. ResultsAmong the circRNAs, we were particularly interested in hsa_circ_0000228, which is spliced from the oncogene ZEB1. In BC cell lines, CircZEB1 expression was upregulated. CircZEB1 knockdown prevented BC cells from migrating and invading, as well as HUVECs from forming tubes and developing. By sponging miR-337–3p, functional testing revealed that circZEB1 promoted O-GlcNAcylation, increased YBX1, and OGT expression. Moreover, circZEB1 overexpression is reversible, in contrast to YBX1 knockdown, which mostly results in the downregulation of multiple oncogenes. ConclusionOur study indicate that circZEB1 had oncogenic function in BC by focusing on circZEB1/miR-337–3p/OGT and YBX1. It might be inferred that circZEB1 could be a promising new target for BC treatment.

Microglial forkhead box O3a deficiency attenuates LPS-induced neuro-inflammation and depressive-like behaviour through regulating the expression of peroxisome proliferator-activated receptor-γ.

Depression is closely linked with microglial activation and neuro-inflammation. Peroxisome proliferator-activated receptor-γ (PPAR-γ) plays an important role in M2 activation of microglia. Forkhead box (FOX) O3a has been implicated in the regulation of mood-relevant behaviour. However, little is known about the inflammatory mechanisms of in the microglia of the brain. Here, we have investigated the role of microglial FOXO3a/PPAR-γ in the development of depression. The effect of FOXO3a on microglia inflammation was analysed in vitro and in lipopolysaccharide (LPS)-induced depression-like behaviours in vivo. ChIP-seq and Dual-luciferase reporter assays were used to confirm the interaction between FOXO3a and PPAR-γ. Behavioural changes were measured, while inflammatory cytokines, microglial phenotype and morphological properties were determined by ELISA, qRT-PCR, western blotting and immunostaining. Overexpression of FOXO3a significantly attenuated expression of PPAR-γ and enhanced the microglial polarization towards the M1 phenotype, while knockdown of FOXO3a had the opposite effect. FOXO3a binds to the promoters of PPAR-γ and decreases its transcription activity. Importantly, deacetylation and activation of FOXO3a regulate LPS-induced neuro-inflammation by inhibiting the expression of PPAR-γ in microglia cells, supporting the antidepressant potential of histone deacetylase inhibitors. Microglial FOXO3a deficiency in mice alleviated LPS-induced neuro-inflammation and depression-like behaviours but failed to reduce anxiety behaviour, whereas pharmacological inhibition of PPAR-γ by GW9662 restored LPS-induced microglial activation and depressive-like behaviours in microglial FOXO3a-deficient mice. FOXO3a/PPAR-γ axis plays an important role in microglial activation and depression, identifying a new therapeutic avenue for the treatment of major depression.

The transcriptional regulation of a putative hemicellulose gene, PtrPARVUS2 in poplar

The plant cell wall serves as a critical interface between the plant and its environment, offering protection against various stresses and contributing to biomass production. Hemicellulose is one of the major components of the cell wall, and understanding the transcriptional regulation of its production is essential to fully understanding cell wall formation. This study explores the regulatory mechanisms underlying one of the genes involved in hemicellulose biosynthesis, PtrPARVUS2. Six transcription factors (TFs) were identified from a xylem-biased library to negatively regulate PtrPARVUS2 expression. These TFs, belonging to diverse TF families, were confirmed to bind to specific cis-elements in the PtrPARVUS2 promoter region, as validated by Yeast One-Hybrid (Y1H) assays, transient expression analysis, and Chromatin Immunoprecipitation sequencing (ChIP-seq) assays. Furthermore, motif analysis identified putative cis-regulatory elements bound by these TFs, shedding light on the transcriptional regulation of SCW biosynthesis genes. Notably, several TFs targeted genes encoding uridine diphosphate glycosyltransferases (UGTs), crucial enzymes involved in hemicellulose glycosylation. Phylogenetic analysis of UGTs regulated by these TFs highlighted their diverse roles in modulating hemicellulose synthesis. Overall, this study identifies a set of TFs that regulate PARVUS2 in poplar, providing insights into the intricate coordination of TFs and PtrPARVUS2 in SCW formation. Understanding these regulatory mechanisms enhances our ability to engineer plant biomass for tailored applications, including biofuel production and bioproduct development.

Abstract PO1-24-05: EZH2 histone methyltransferase activity promotes spindle metaplastic breast carcinoma metastasis and is induced by CCN6 knockout/b-catenin/TCF axis

Abstract Background: Metaplastic breast carcinomas (mBrCAs) are a highly aggressive subtype of triple negative breast cancer with histological evidence of deregulated differentiation towards non-glandular components. Previous studies have demonstrated that human mBrCA often exhibit activation of differentiation pathways, including canonical Wnt/b-catenin and EZH2-mediated transcriptional repression. Our lab has identified CCN6 as a tumor suppressor in mBrCA. MMTV-Cre;Ccn6fl/fl (CCN6KO) mice develop spindle mBrCAs and CCN6 is reduced/lost in 68% of human mBrCAs. Recently, we have demonstrated that the tumor suppressor function of CCN6 in spindle mBrCA requires activation of Wnt/b-catenin signaling pathway. Here, we tested the hypothesis that CCN6 KO leads to the upregulation of EZH2 histone methyltransferase promoting mBrCA. Furthermore, we investigated the requirement for the activation of the Wnt/b-catenin pathway in this mechanism. Methods: To test the effect of CCN6/b-catenin on EZH2 gene and protein expression we performed IHC, IF, qRT-PCR, ChIP-Seq, RNA-seq, invasion and adhesion assays, EZH2 reporter assay and immunoblots in mBrCA cell lines and MMTV-Cre;Ccn6fl/fl tumors. To investigate the role of CCN6KO-induced b-catenin activation on EZH2 activity and neoplastic functions we employed three independent approaches: i) Expression of a dominant-negative Tcf4 (dnTcf4) rescued with EZH2-WT, dSET and dNLS domain mutants versus vector in MMTV-Cre;Ccn6fl/fl tumor-derived cells; ii) Expression of a constitutively active mutant (S33Y) b-catenin in concert with treatment with recombinant human CCN6 (rhCCN6; 500 ug/ml) versus control; iii) Syngeneic orthotopic mammary tumor transplants of MMTV-Cre;Ccn6fl/fl were used for in vivo rescue experiments with rhCCN6 or BSA. To assess the therapeutic benefit of inhibiting EZH2 methyltransferase activity, MMTV-Cre;Ccn6fl/fl tumor cells (CCN6KO cells) were implanted orthotopically or intracardially in FVB mice, followed by treatment with EPZ-6438 (a selective EZH2 methyltransferase inhibitor) or vehicle. We tested CCN6, b-catenin, and EZH2 expression by IHC in a cohort of 27 human mBrCA tumor samples. Results: CCN6KO-induced b-catenin/TCF activation mediates EZH2 transcriptional upregulation and the deposition of repressive H3K27me3 in spindle mBrCAs. We found that the invasion program triggered by CCN6KO-induced Wnt/b-catenin requires EZH2 catalytic activity as WT-EZH2 (but not dSET-EZH2) rescued the reduced invasion of CCN6KO-dnTcf4 cells. RNA-seq and 3H3K27 ChIP-seq of CCN6KO-dnTcf4 cells transduced with EZH2-WT or dSET-EZH2 identify specific CCN6KO-b-catenin/Tcf targets that require EZH2 transcriptional repressor function. In vivo, administration of CCN6 protein to MMTV-Cre;Ccn6fl/fl tumor transplants reduces tumor growth and nuclear b-catenin, EZH2 and 3H3K27 in the tumors. Pharmacologic inhibition of EZH2 reduces the growth and metastasis of CCN6KO mBrCA tumors and improves survival. We identify a subset of human spindle mBrCA (54%) that display a CCN6Low/nuclear b-cat/EZH2High phenotype. Conclusion: We found a critical role for EZH2 activation in CCN6-deficient mBrCA tumor phenotypes via b-catenin/TCF Wnt canonical signaling. We demonstrate the effectiveness of pharmacological inhibition of EZH2 methyltransferase activity in reducing primary tumor growth and distant metastasis in mouse models of spindle mBrCA. In clinical samples, low CCN6 is significantly associated with activated b-catenin and high EZH2 in spindle mBrCAs compared to other subtypes. These data reveal a novel tumor suppressor mechanism of CCN6 and provide compelling evidence supporting the potential therapeutic value of CCN6 restoration, b-catenin or EZH2 inhibition as promising approaches for the treatment of spindle mBrCAs. Citation Format: Maria Gonzalez, Ahmad Eido, GIUSEPPINA AUGIMERI, Celina Kleer. EZH2 histone methyltransferase activity promotes spindle metaplastic breast carcinoma metastasis and is induced by CCN6 knockout/b-catenin/TCF axis [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO1-24-05.

Abstract 2004: CDC7 inhibition prevents neuroendocrine transformation in the lung and prostate through MYC degradation

Abstract Neuroendocrine (NE) transformation occurs as a mechanism of resistance to targeted therapy in up to 14% and 30% of EGFR-mutant lung and AR-dependent prostate adenocarcinomas, respectively, leading to poor prognosis. Even if we know the tumor population at high risk of transformation (TP53/RB1-mutated), no therapies to prevent NE relapse are currently available. To identify therapeutic vulnerabilities for tumors undergoing NE transformation, we performed an in vitro CRISPR screen in a NE-transformed lung tumor model. This screen identified CDC7, involved in DNA replication and DNA damage response, as a potential therapeutic target in this setting. Proteogenomic analyses revealed CDC7 upregulation in lung and prostate clinical samples undergoing NE transformation, detected already in pre-transformation adenocarcinomas. These results indicated that CDC7 expression is induced already at early steps of transformation. Consistently, TP53/RB1-mutated lung and prostate adenocarcinomas exhibited higher CDC7 expression than their double wild-type counterparts. Combined ChIP-seq and promoter reporter assays indicated that CDC7 expression was directly regulated by TP53 and RB1 inactivation. Importantly, TP53/RB1-inactivation induced sensitivity to the CDC7 inhibitor simurosertib, unraveling a therapeutic vulnerability in tumors at high risk of NE transformation. Thus, we tested the combination of simurosertib with targeted therapy in vivo in different lung and prostate patient-derived models of NE transformation, namely TP53/RB1-knock out adenocarcinomas known to undergo NE transformation on targeted therapy. In these, simurosertib was able to suppress NE transformation and dramatically delay tumor relapse. Trajectory analysis on single-cell transcriptomic data for such models revealed a NE transformation transcriptional program occurring already in the untreated tumors before transformation. Remarkably, simurosertib treatment reverted this transcriptomic state and induced a reversion to the original adenocarcinoma transcriptomic profile. CDC7 inhibition led to increased proteasomal activity and degradation of MYC, a stemness transcription factor involved in NE transformation. Ectopic overexpression of MYCT58A, a proteasome degradation-resistant MYC isoform, rescued the NE phenotype in these transformation models, suggesting that CDC7 inhibition-induced MYC degradation is the mechanism by which NE transformation is prevented. In sum, CDC7 inhibition may suppress, or at least dramatically delay NE transformation in patients with lung and prostate adenocarcinomas at high risk of transformation, by inducing MYC proteasomal degradation. The clinical availability of CDC7 inhibitors, currently in phase II clinical trials after demonstrating tolerability and preliminary efficacy, will allow rapid translation of these results into the clinics. Citation Format: Alvaro Quintanal-Villalonga, Fathema Uddin, Kenta Kawasaki, Esther Redin, Vidushi Durani, Amin Sabet, Wouter Karthaus, Yingqian A. Zhan, Samir Zaidi, Moniquetta Shaffer, Harsha Sridhar, Juan Qiu, Parvathy Manoj, Elisa De Stanchina, Michael C. Haffner, Charles L. Sawyers, Charles M. Rudin. CDC7 inhibition prevents neuroendocrine transformation in the lung and prostate through MYC degradation [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2004.

WWP2 Regulates Renal Fibrosis and the Metabolic Reprogramming of Profibrotic Myofibroblasts.

Renal fibrosis is a common pathological endpoint in CKD that is challenging to reverse, and myofibroblasts are responsible for the accumulation of a fibrillar collagen-rich extracellular matrix (ECM). Recent studies have unveiled myofibroblasts diversity in terms of proliferative and fibrotic characteristics, which are linked to different metabolic states. We previously demonstrated the regulation of ECM genes and tissue fibrosis by WWP2, a multifunctional E3 ubiquitin-protein ligase. Here, we investigate WWP2 in renal fibrosis and in the metabolic reprograming of myofibroblasts in CKD. We used kidney samples from CKD patients and WWP2 -null kidney disease mice models, and leveraged single cell RNA-seq analysis to detail the cell-specific regulation of WWP2 in fibrotic kidneys. Experiments in primary cultured myofibroblasts by bulk-RNA seq, ChIP-seq, metabolomics and cellular metabolism assays, were used to study the metabolic regulation of WWP2 and its downstream signaling. The tubulointerstitial expression of WWP2 was associated with fibrotic progression in CKD patients and in murine kidney disease models. WWP2 deficiency promoted myofibroblast proliferation and halts pro-fibrotic activation, reducing the severity of kidney fibrosis in vivo . In renal myofibroblasts, WWP2 deficiency increased fatty acid oxidation and activated the pentose phosphate pathway, boosting mitochondrial respiration at the expense of glycolysis. WWP2 suppressed the transcription of PGC-1α, a metabolic mediator of fibrotic response, and pharmacological inhibition of PGC-1α partially abrogated the protective effects of WWP2 deficiency on myofibroblasts. WWP2 regulates the metabolic reprogramming of profibrotic myofibroblasts by a WWP2-PGC-1α axis, and WWP2 deficiency protects against kidney fibrosis in CKD.

Long non-coding RNA lung cancer-associated transcript-1 promotes glioblastoma progression by enhancing Hypoxia-inducible factor 1 alpha activity.

Hypoxia is associated with poor prognosis in many cancers including glioblastoma (GBM). Glioma stem-like cells (GSCs) often reside in hypoxic regions and serve as reservoirs for disease progression. Long non-coding RNAs (lncRNAs) have been implicated in GBM. However, the lncRNAs that modulate GSC adaptations to hypoxia are poorly understood. Identification of these lncRNAs may provide new therapeutic strategies to target GSCs under hypoxia. lncRNAs induced by hypoxia in GSCs were identified by RNA-seq. Lung cancer-associated transcript-1 (LUCAT1) expression was assessed by qPCR, RNA-seq, Northern blot, single molecule FISH in GSCs, and interrogated in IvyGAP, The Cancer Genome Atlas, and CGGA databases. LUCAT1 was depleted by shRNA, CRISPR/Cas9, and CRISPR/Cas13d. RNA-seq, Western blot, immunohistochemistry, co-IP, ChIP, ChIP-seq, RNA immunoprecipitation, and proximity ligation assay were performed to investigate mechanisms of action of LUCAT1. GSC viability, limiting dilution assay, and tumorigenic potential in orthotopic GBM xenograft models were performed to assess the functional consequences of depleting LUCAT1. A new isoform of Lucat1 is induced by Hypoxia inducible factor 1 alpha (HIF1α) and Nuclear factor erythroid 2-related factor 2 (NRF2) in GSCs under hypoxia. LUCAT1 is highly expressed in hypoxic regions in GBM. Mechanistically, LUCAT1 formed a complex with HIF1α and its co-activator CBP to regulate HIF1α target gene expression and GSC adaptation to hypoxia. Depletion of LUCAT1 impaired GSC self-renewal. Silencing LUCAT1 decreased tumor growth and prolonged mouse survival in GBM xenograft models. A HIF1α-LUCAT1 axis forms a positive feedback loop to amplify HIF1α signaling in GSCs under hypoxia. LUCAT1 promotes GSC self-renewal and GBM tumor growth. LUCAT1 is a potential therapeutic target in GBM.

Folate deficiency reduced aberrant level of DOT1L-mediated histone H3K79 methylation causes disruptive SHH gene expression involved in neural tube defects

BackgroundNeural tube defects (NTDs) are one of the most severe congenital abnormalities characterized by failures of the neural tube to close during early embryogenesis. Maternal folate deficiency could impact the occurrence of NTDs, however, the mechanisms involved in the cause of NTDs are poorly defined.ResultsHere, we report that histone H3 methyltransferase disruptor of telomeric silencing 1-like (DOT1L) expression was significantly downregulated, and low levels of H3K79me2 were found in the corresponding NTDs samples with their maternal serum folate under low levels. Using ChIP-seq assays, we found that a decrease of H3K79me2 downregulates the expression of Shh and Sufu in mouse embryonic stem cells (mESC) under folate deficiency. Interestingly, folate antagonist methotrexate treatment led to attenuation of H3K79me2 due to Dot1l, affecting Shh and Sufu genes regulation. Upon further analysis, we find that the genes Shh and Sufu are both downregulated in the brain tissues of mice and humans with NTDs. There was a positive correlation between the transcription levels of Shh, Sufu and the protein levels of DOT1L by Pearson correlation analysis.ConclusionOur results indicate that abnormal Shh and Sufu genes expression reduced by aberrant Dot1l-mediated H3K79me2 levels could be the cause of NTDs occurrence.

Non-Canonical Transcription Activator Activity of RNA-Binding Protein MBNL1 Uniquely Controls FLT3 expression in MLL-Rearranged Leukemias

MLL-rearranged (MLL-r) leukemias count for more than 80% of infant leukemia, ~5-10% of B-cell acute lymphoblastic leukemia (B-ALL), and ~10% of acute myeloid leukemia (AML) cases, where they confer a particularly poor outcome. Despite treatment with intensive multi-agent chemotherapy, most MLL-r patients ultimately relapsed after an initial remission. So far, the molecular mechanism by which MLL-r leukemia maintains progression and prevents differentiation remains largely unclear. Observations in the clinic support aberrant expression and concomitant activating driver mutations in the gene encoding the tyrosine kinase FLT3 that occur in leukemia, including the MLL-r subtype. To interrogate the novel MLL-r vulnerable genes regulating FLT3, we developed an algorithm so-called “context-dependent association analysis” (CDAA) to explore potential MLL-r dependent FLT3 regulators using the DepMap dataset. The RNA-binding protein MBNL1 was identified as the top candidate. MBNL1 has been recognized as an RNA-binding protein (RBP) involved in splicing, RNA export, and stability. However, our study suggested its non-canonical transcriptional activator function by uniquely controlling FLT3 expression in MLL-r leukemias. CRISPR/Cas9 disruption of the coding region of MBNL1 or CRISPR-interference (CRISPRi) against the MBNL1-bound site in the FLT3 distal enhancer notably decreased FLT3 mRNA expression in MLL-r leukemia cell lines. However, the FLT3 gene splicing defect was not observed upon MBNL1 loss, suggesting MBNL1 could be affecting FLT3 expression through a transcriptional regulation mechanism. Genome-wide ChIP-seq assays performed in MLL-r SEM and MOLM13 cells identified only two reproducible MBNL1-bound peaks at a genome-wide scale, one of which was located in the distal enhancer of FLT3, ~170kb away from the FLT3 promoter. CRISPR-interference (CRISPRi) against the MBNL1-bound site in the FLT3 distal enhancer notably decreased FLT3 mRNA expression in MLL-r leukemia cell lines without affecting enhancer-promoter looping. Using in vivo assay, SEM cells targeted with sgRNA against the MBNL1-bound site in the distal enhancer of FLT3 transplanted into NSG mice exhibited significant growth retardation. Furthermore, the impaired cell fitness crisis phenotype could be entirely rescued by overexpression of FLT3 cDNA. To identify the protein domains of MBNL1 required for its oncogenic function in MLL-r leukemia cells, an MBNL1 domain CRISPR screen was performed on the SEM cell line, stably expressing Cas9 and a pooled sgRNA library targeting the coding region of MBNL1. Essential domains observed from the screen included MBNL1's zinc finger domains. A biochemistry electrophoretic mobility shift assay (EMSA) was conducted using purified recombinant MBNL1 protein and a 207 bp DNA oligo synthesized based on the MBNL1 ChIP-seq peak at the FLT3 enhancer. The result clearly demonstrated that MBNL1 could directly bind to the enhancer DNA of FLT3 through its zinc fingers 1 and 2 domains, supporting its role as a novel transcriptional regulator of FLT3. Using state-of-art research tools, including genome editing, ChIP-seq, and EMSA assays, we have systematically interrogated an undocumented transcriptional regulation axis of MBNL1/FLT3 in MLL-r leukemias. This work will significantly advance our understanding of how transcription programs and novel factors reconstruct the regulatory network specifying MLL-r dependency, promoting the future development of alternative therapeutic targets for MLL-r leukemia in patients.

Polycom Group Protein Mel18 Inhibits Hematopoietic Stem Cell Self-Renewal through Repressing the Expression of Genes Important for Cellular Senescence and Proliferation

Polycomb group (PcG) proteins are epigenetic gene silencers that have been implicated in stem cell maintenance and cancer development. Genetic and biochemical studies indicate that Polycomb group proteins exist in at least two protein complexes, Polycomb repressive complex 2 (PRC2) and Polycomb repressive complex 1 (PRC1), that act in concert to initiate and maintain stable gene repression. Bmi1 and Mel18 are two major homologs of the PCGF subunit within the canonical PRC1 complex. While Bmi1 is a positive regulator of hematopoietic stem cell (HSC) and leukemia stem cell (LSC) self-renewal, the role of Mel18 in normal and malignant hematopoiesis is not fully understood. To further determine the role of Mel18 in hematopoiesis, we have generated Mel18 conditional knockout mice ( Mel18 f/f-Mx1Cre +). Acute deletion of Mel18 in the hematopoietic compartment did not affect the frequency and survival of hematopoietic stem and progenitor cells (HSPCs). To determine the impact of Mel18 deficiency on HSPC proliferation, we performed serial replating assays and found that loss of Mel18 increases the replating potential of HSPCs in vitro. To determine the role of Mel18 in HSC self-renewal, we first performed serial bone marrow (BM) transplantation assays and found that Mel18 null BM cells show increased engraftment in both primary and secondary transplantation assays. Limiting dilution transplantation showed that there were more functional HSCs in the BM of Mel18 knockout mice. We then performed HSC transplantation assays and found that loss Mel18 enhances HSC self-renewal. Thus, we demonstrated that Mel18 inhibits HSC self-renewal and decreases the number of functional HSCs in vivo. To understand the mechanism by which Mel18 inhibits HSC self-renewal and proliferation, we performed RNA-seq using Mel18 +/+ and Mel18 -/- HSCs. We then performed gene set enrichment analysis (GSEA) to identify potential Mel18 target genes important for HSC behavior. We found that HSC gene signatures are enriched in Mel18 -/- HSCs. We also found that leukemia stem cell gene signatures and CBF-mutant-AML genes are significantly enriched in Mel18 -/- HSCs, suggesting that loss of Mel18 may prime normal HSCs for leukemic transformation. Notably, our RNA-seq assay showed that senescence-related genes, including S100a9 and S100a8, were down-regulated in Mel18 -/- HSPCs. Importantly, loss of Mel18 significantly decreased the number of senescent HSPCs. To understand how Mel18 regulates gene transcription, we performed ATAC-seq assays in Mel18 +/+ and Mel18 -/- HSPCs. Mel18 null HSPCs showed decreased chromatin accessibility near both S100a8 and S100a9 genes. We confirmed that both S100a8 and S100a9 are downregulated in Mel18 -/- HSPCs compared to Mel18 +/+ HSPCs. Given that cell cycle arrest is a common feature of cellular senescence, we also examined the expression of cell cycle related genes in Mel18 -/- HSPCs. The expression of cell cycle regulators, including CDK4 and CCND2, is increased in Mel18 -/- HSPCs. Mel18 can bind to a specific DNA sequence, 5‘-GACTNGACT-3‘, thereby repressing the expression of target genes. Mel18 directly binds to the CCND2 gene locus as revealed by ChIP-seq assays in embryonic stem cells, suggesting that Mel18 may directly repress CCND2 expression in HSPCs. In summary, we demonstrate that Mel18 inhibits HSC self-renewal via repressing the expression of genes that are important for cellular senescence and proliferation.

Impaired Early Spliceosome Complex Assembly Underlies Gene Body Elongation Transcription Defect in SF3B1K700E

Introduction: Recurrent mutations in splicing factors (SFs) such as SF3B1, U2AF1 and SRSF2 are common in clonal myeloid disorders. These mutations are typically hemizygous, exclusive to each other, and non-synonymous, pointing to neomorphic functions. Molecular mechanisms connecting the mutations to disease pathogenesis remain unclear. Previous studies predominantly focused on the well-established role of splicing factors in pre-mRNA splicing, seeking to explain how mutations in these SFs leading to changes in the abundance of key regulatory genes through their effects on alternative splicing. These studies are limited by inconsistent and relatively low isoform changes across mutations and independent datasets. They also don't explain the mutual exclusivity of these mutations. In this study, we explored a novel mechanism that link SF mutations to disease biology: dysregulation of RNA Polymerase II (Pol2) transcriptional kinetics arising from impaired spliceosome assembly. Methods: We generated inducible isogenic K562 cell lines that express single mutant alleles (SF3B1 K700E and U2AF1 S34F,and corresponding wild-type controls). A novel combinatory approach of AAV-intron trap with CRISPR/Cas9 (Boddu et al, Comms Biol, 2021) was utilized for gene editing. RNAPII transcription kinetics changes were assessed, at an early post-induction time point of 4 days (Fig1A), using ChIP-seq, and nascent-transcriptome assays (including GRO-seq and transient transcriptome-time lapse sequencing or TT-TL seq). Spliceosome assembly kinetics were assayed with non-denaturing gels. Co-transcriptional splicing efficiency (CoSE) was assessed using nascent long-read sequencing (LRS) (Reimer, Mol Cell 2021). Key findings in K562 cell lines were validated in CD34+ progenitor populations from MDS-derived patient samples (and healthy volunteer controls). Results: ChIP-seq for Pol 2 (both total (Fig1B) and elongation-specific Ser2P (Fig1C) showed a redistribution of Pol2 into gene-bodies, suggesting a failure of Pol2 elongation in SF3B1 K700E. This Pol2 redistribution was dependent on gene length and intron number (long, multi-intron genes most affected). GRO-seq confirmed the ChIP-seq results (Fig1D). Additionally, spike-in normalized TT-seq showed significantly reduced nascent transcription in gene bodies in SF3B1 K700E (Fig 1E), consistent with a gene-body elongation rate defect. A similar redistribution of Pol2 into gene bodies was noted by ChIP-seq and TT-seq in U2AF1S34F, pointing to a convergent transcriptional response among distinct SF mutations. Genome-wide CoSE was noted to be lower in SF3B1 K700Eindicating an associated splicing assembly defect (Fig1F). Changes to transcription due to SF3B1 mutations were similarly observed, using low-input Pol2 CUT&RUN, in CD34+ progenitor populations isolated from human MDS patient samples (Fig1G). Of pathologic relevance, this transcriptional dysregulation led to generation of R-loops, DNA damage, and S-phase arrest. Given close coordination of splicing and transcription, and role of SF3B1 and U2AF1 in early spliceosome formation, we speculated if mutant SF3B1 leads to defective assembly of early spliceosomes. Compared to wild-type protein, SF3B1 K700E was noted to have significantly impaired interactions with U2AF1, U2AF2, HTATSF1 (Fig1H). We then determined how spliceosome assembly kinetics of SF3B1 K700E differed from those of SF3B1 WTusing an in vitro splicing assay. The transition from E- to A- complex intermediates was reduced in SF3B1 K700E(Fig1I). Furthermore, disrupting the early pre-spliceosome assembly (through loss of DDX46 or HTATSF1 using both siRNA and inducible shRNA (Fig1J) knockdown systems) resulted in a transcriptional response similar to SF3B1 K700E. Based on recent cryo-EM structures of Pol II-U1 snRNP interactions (Zhang et al, Science, 2021), we speculate that the impaired transition to A/B complex in SF3B1 K700E prevents the release of U1 snRNP and release of Pol2, leading to intron-predominant defective Pol II elongation (Fig1K). Conclusion: Here, we found that SF3B1 and U2AF1 mutations impair Pol2 transcription elongation along gene bodies, leading to pathologic R-loops and replication stress. This elongation results from disrupted pre-spliceosome assembly due to impaired protein-protein interactions of mutant SF3B1.

Lineage commitment of dermal fibroblast progenitors is controlled by Kdm6b-mediated chromatin demethylation.

Dermal Fibroblast Progenitors (DFPs) differentiate into distinct fibroblast lineages during skin development. However, the epigenetic mechanisms that regulate DFP differentiation are not known. Our objective was to use multimodal single-cell approaches, epigenetic assays, and allografting techniques to define a DFP state and the mechanism that governs its differentiation potential. Our initial results indicated that the overall transcription profile of DFPs is repressed by H3K27me3 and has inaccessible chromatin at lineage-specific genes. Surprisingly, the repressive chromatin profile of DFPs renders them unable to reform the skin in allograft assays despite their multipotent potential. We hypothesized that chromatin derepression was modulated by the H3K27me3 demethylase, Kdm6b/Jmjd3. Dermal fibroblast-specific deletion of Kdm6b/Jmjd3 in mice resulted in adipocyte compartment ablation and inhibition of mature dermal papilla functions, confirmed by additional single-cell RNA-seq, ChIP-seq, and allografting assays. We conclude that DFPs are functionally derepressed during murine skin development by Kdm6b/Jmjd3. Our studies therefore reveal a multimodal understanding of how DFPs differentiate into distinct fibroblast lineages and provide a novel publicly available multiomics search tool.

CALMODULIN6 negatively regulates cold tolerance by attenuating ICE1-dependent stress responses in tomato.

Chilling temperatures induce an increase in cytoplasmic calcium (Ca2+) ions to transmit cold signals, but the precise role of Calmodulins (CaMs), a type of Ca2+ sensor, in plant tolerance to cold stress remains elusive. In this study, we characterized a tomato (Solanum lycopersicum) CaM gene, CALMODULIN6 (CaM6), which responds to cold stimulus. Overexpressing CaM6 increased tomato sensitivity to cold stress whereas silencing CaM6 resulted in a cold-insensitive phenotype. We showed that CaM6 interacts with Inducer of CBF expression 1 (ICE1) in a Ca2+-independent process and ICE1 contributes to cold tolerance in tomato plants. By integrating RNA-sequencing (RNA-seq) and chromatin immunoprecipitation-sequencing (ChIP-seq) assays, we revealed that ICE1 directly altered the expression of 76 downstream cold-responsive (COR) genes that potentially confer cold tolerance to tomato plants. Moreover, the physical interaction of CaM6 with ICE1 attenuated ICE1 transcriptional activity during cold stress. These findings reveal that CaM6 attenuates the cold tolerance of tomato plants by suppressing ICE1-dependent COR gene expression. We propose a CaM6/ICE1 module in which ICE1 is epistatic to CaM6 under cold stress. Our study sheds light on the mechanism of plant response to cold stress and reveals CaM6 is involved in the regulation of ICE1.

GATA2 co-opts TGFβ1/SMAD4 oncogenic signaling and inherited variants at 6q22 to modulate prostate cancer progression

BackgroundAberrant somatic genomic alteration including copy number amplification is a hallmark of cancer genomes. We previously profiled genomic landscapes of prostate cancer (PCa), yet the underlying causal genes with prognostic potential has not been defined. It remains unclear how a somatic genomic event cooperates with inherited germline variants contribute to cancer predisposition and progression.MethodsWe applied integrated genomic and clinical data, experimental models and bioinformatic analysis to identify GATA2 as a highly prevalent metastasis-associated genomic amplification in PCa. Biological roles of GATA2 in PCa metastasis was determined in vitro and in vivo. Global chromatin co-occupancy and co-regulation of GATA2 and SMAD4 was investigated by coimmunoprecipitation, ChIP-seq and RNA-seq assays. Tumor cellular assays, qRT-PCR, western blot, ChIP, luciferase assays and CRISPR-Cas9 editing methods were performed to mechanistically understand the cooperation of GATA2 with SMAD4 in promoting TGFβ1 and AR signaling and mediating inherited PCa risk and progression.ResultsIn this study, by integrated genomics and experimental analysis, we identified GATA2 as a prevalent metastasis-associated genomic amplification to transcriptionally augment its own expression in PCa. Functional experiments demonstrated that GATA2 physically interacted and cooperated with SMAD4 for genome-wide chromatin co-occupancy and co-regulation of PCa genes and metastasis pathways like TGFβ signaling. Mechanistically, GATA2 was cooperative with SMAD4 to enhance TGFβ and AR signaling pathways, and activated the expression of TGFβ1 via directly binding to a distal enhancer of TGFβ1. Strinkingly, GATA2 and SMAD4 globally mediated inherited PCa risk and formed a transcriptional complex with HOXB13 at the PCa risk-associated rs339331/6q22 enhancer, leading to increased expression of the PCa susceptibility gene RFX6.ConclusionsOur study prioritizes causal genomic amplification genes with prognostic values in PCa and reveals the pivotal roles of GATA2 in transcriptionally activating the expression of its own and TGFβ1, thereby co-opting to TGFβ1/SMAD4 signaling and RFX6 at 6q22 to modulate PCa predisposition and progression.

Profibrotic role of the SOX9–MMP10–ECM biosynthesis axis in the tracheal fibrosis after injury and repair

Fibroblast activation and extracellular matrix (ECM) deposition play an important role in the tracheal abnormal repair process and fibrosis. As a transcription factor, SOX9 is involved in fibroblast activation and ECM deposition. However, the mechanism of how SOX9 regulates fibrosis after tracheal injury remains unclear. We investigated the role of SOX9 in TGF-β1-induced fibroblast activation and ECM deposition in rat tracheal fibroblast (RTF) cells. SOX9 overexpression adenovirus (Ad-SOX9) and siRNA were transfected into RTF cells. We found that SOX9 expression was up-regulated in RTF cells treated with TGF-β1. SOX9 overexpression activated fibroblasts and promoted ECM deposition. Silencing SOX9 inhibited cell proliferation, migration, and ECM deposition, induced G2 arrest, and increased apoptosis in RTF cells. RNA-seq and chromatin immunoprecipitation sequencing (ChIP-seq) assays identified MMP10, a matrix metalloproteinase involved in ECM deposition, as a direct target of SOX9, which promotes ECM degradation by increasing MMP10 expression through the Wnt/β-catenin signaling pathway. Furthermore, in vivo, SOX9 knockdown ameliorated granulation proliferation and tracheal fibrosis, as manifested by reduced tracheal stenosis. In conclusion, our findings indicate that SOX9 can drive fibroblast activation, cell proliferation, and apoptosis resistance in tracheal fibrosis via the Wnt/β-catenin signaling pathway. The SOX9–MMP10–ECM biosynthesis axis plays an important role in tracheal injury and repair. Targeting SOX9 and its downstream target MMP10 may represent a promising therapeutic approach for tracheal fibrosis.

Epigenomics Analysis of the Suppression Role of SIRT1 via H3K9 Deacetylation in Preadipocyte Differentiation.

Sirtuin 1 (SIRT1) overexpression significantly inhibits lipid deposition during yak intramuscular preadipocyte (YIMA) differentiation; however, the regulatory mechanism remains unknown. We elucidated the role of SIRT1 in YIMA differentiation using lentivirus-mediated downregulation technology and conducted mRNA-seq and ChIP-seq assays using H3K9ac antibodies after SIRT1 overexpression in order to reveal SIRT1 targets during YIMA adipogenesis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed in order to identify the functional annotation of common genes. In addition, a potential target of SIRT1 was selected to verify its effects on the differentiation and proliferation of YIMAs. SIRT1 interfered with lipid deposition and promoted YIMA differentiation. In total, 143,518 specific peaks were identified after SIRT1 overexpression, where genes associated with downregulation peaks were enriched in transcription, gene expression, lipid-related processes, and classical lipid-related pathways. The H3K9ac signal in the whole genome promoter region (2 kb upstream and downstream of the transcription start site (TSS)) was weakened, and the peaks were distributed across all gene functional regions. Genes that lost signals in their TSS region or gene body region were enriched in both biological processes and pathways associated with lipogenesis. The ChIP-seq results revealed 714 common differential genes in mRNA-seq, which were enriched in "MAPK signaling", "lipid and atherosclerosis", "mTOR signaling", and "FoxO signaling" pathways. A total of 445 genes were downregulated in both their H3K9ac signals and mRNA expression, and one of their most significantly enriched pathways was FoxO signaling. Nine genes (FBP2, FPGT, HSD17B11, KCNJ15, MAP3K20, SLC5A3, TRIM23, ZCCHC10, and ZMYM1) lost the H3K9ac signal in their TSS regions and had low mRNA expression, and three genes (KCNJ15, TGM3, and TRIM54) had low expression but lost their H3K9ac signal in the gene body region. The interference of TRIM23 significantly inhibited fat deposition during preadipocyte differentiation and promoted cell proliferation by increasing S-phase cell numbers. The present study provides new insights into the molecular mechanism of intramuscular fat content deposition and the epigenetic role of SIRT1 in adipocyte differentiation.