KLF4 inhibits early neural differentiation of ESCs by coordinating specific 3D chromatin structure.

Neural differentiation of embryonic stem cells (ESCs) requires precisely orchestrated gene regulation, a process governed in part by changes in 3D chromatin structure. How these changes regulate gene expression in this context remains unclear. In this study, we observed enrichment of the transcription factor KLF4 at some poised or closed enhancers at TSS-linked regions of genes associated with neural differentiation. Combination analysis of ChIP, HiChIP and RNA-seq data indicated that KLF4 loss in ESCs induced changes in 3D chromatin structure, including increased chromatin interaction loops between neural differentiation-associated genes and active enhancers, leading to upregulated expression of neural differentiation-associated genes and therefore early neural differentiation. This study suggests KLF4 inhibits early neural differentiation by regulation of 3D chromatin structure, which is a new mechanism of early neural differentiation.

MMP-9 facilitates selective proteolysis of the histone H3 tail at genes necessary for proficient osteoclastogenesis.

Although limited proteolysis of the histone H3 N-terminal tail (H3NT) is frequently observed during mammalian differentiation, the specific genomic sites targeted for H3NT proteolysis and the functional significance of H3NT cleavage remain largely unknown. Here we report the first method to identify and examine H3NT-cleaved regions in mammals, called chromatin immunoprecipitation (ChIP) of acetylated chromatin (ChIPac). By applying ChIPac combined with deep sequencing (ChIPac-seq) to an established cell model of osteoclast differentiation, we discovered that H3NT proteolysis is selectively targeted near transcription start sites of a small group of genes and that most H3NT-cleaved genes displayed significant expression changes during osteoclastogenesis. We also discovered that the principal H3NT protease of osteoclastogenesis is matrix metalloproteinase 9 (MMP-9). In contrast to other known H3NT proteases, MMP-9 primarily cleaved H3K18-Q19 in vitro and in cells. Furthermore, our results support CBP/p300-mediated acetylation of H3K18 as a central regulator of MMP-9 H3NT protease activity both in vitro and at H3NT cleavage sites during osteoclastogenesis. Importantly, we found that abrogation of H3NT proteolysis impaired osteoclastogenic gene activation concomitant with defective osteoclast differentiation. Our collective results support the necessity of MMP-9-dependent H3NT proteolysis in regulating gene pathways required for proficient osteoclastogenesis.

A HOTAIR regulatory element modulates glioma cell sensitivity to temozolomide through long-range regulation of multiple target genes.

Temozolomide (TMZ) is a frequently used chemotherapy for glioma; however, chemoresistance is a major problem limiting its effectiveness. Thus, knowledge of mechanisms underlying this outcome could improve patient prognosis. Here, we report that deletion of a regulatory element in the HOTAIR locus increases glioma cell sensitivity to TMZ and alters transcription of multiple genes. Analysis of a combination of RNA-seq, Capture Hi-C, and patient survival data suggests that CALCOCO1 and ZC3H10 are target genes repressed by the HOTAIR regulatory element and that both function in regulating glioma cell sensitivity to TMZ. Rescue experiments and 3C data confirmed this hypothesis. We propose a new regulatory mechanism governing glioma cell TMZ sensitivity.

Enhancer architecture-dependent multilayered transcriptional regulation orchestrates RA signaling-induced early lineage differentiation of ESCs.

Signaling pathway-driven target gene transcription is critical for fate determination of embryonic stem cells (ESCs), but enhancer-dependent transcriptional regulation in these processes remains poorly understood. Here, we report enhancer architecture-dependent multilayered transcriptional regulation at the Halr1-Hoxa1 locus that orchestrates retinoic acid (RA) signaling-induced early lineage differentiation of ESCs. We show that both homeobox A1 (Hoxa1) and Hoxa adjacent long non-coding RNA 1 (Halr1) are identified as direct downstream targets of RA signaling and regulated by RARA/RXRA via RA response elements (RAREs). Chromosome conformation capture-based screens indicate that RA signaling promotes enhancer interactions essential for Hoxa1 and Halr1 expression and mesendoderm differentiation of ESCs. Furthermore, the results also show that HOXA1 promotes expression of Halr1 through binding to enhancer; conversely, loss of Halr1 enhances interaction between Hoxa1 chromatin and four distal enhancers but weakens interaction with chromatin inside the HoxA cluster, leading to RA signaling-induced Hoxa1 overactivation and enhanced endoderm differentiation. These findings reveal complex transcriptional regulation involving synergistic regulation by enhancers, transcription factors and lncRNA. This work provides new insight into intrinsic molecular mechanisms underlying ESC fate determination during RA signaling-induced early differentiation.

Three-Dimensional Gene Regulation Network in Glioblastoma Ferroptosis.

Ferroptosis is an iron-dependent form of cell death, which is reported to be associated with glioma progression and drug sensitivity. Targeting ferroptosis is a potential therapeutic approach for glioma. However, the molecular mechanism of glioma cell ferroptosis is not clear. In this study, we profile the change of 3D chromatin structure in glioblastoma ferroptosis by using HiChIP and study the 3D gene regulation network in glioblastoma ferroptosis. A combination of an analysis of HiChIP and RNA-seq data suggests that change of chromatin loops mediated by 3D chromatin structure regulates gene expressions in glioblastoma ferroptosis. Genes that are regulated by 3D chromatin structures include genes that were reported to function in ferroptosis, like HDM2 and TXNRD1. We propose a new regulatory mechanism governing glioblastoma cell ferroptosis by 3D chromatin structure.

CTCF-binding element regulates ESC differentiation via orchestrating long-range chromatin interaction between enhancers and HoxA.

Proper expression of Homeobox A cluster genes (HoxA) is essential for embryonic stem cell (ESC) differentiation and individual development. However, mechanisms controlling precise spatiotemporal expression of HoxA during early ESC differentiation remain poorly understood. Herein, we identified a functional CTCF-binding element (CBE+47) closest to the 3'-end of HoxA within the same topologically associated domain (TAD) in ESC. CRISPR-Cas9-mediated deletion of CBE+47 significantly upregulated HoxA expression and enhanced early ESC differentiation induced by retinoic acid (RA) relative to wild-type cells. Mechanistic analysis by chromosome conformation capture assay (Capture-C) revealed that CBE+47 deletion decreased interactions between adjacent enhancers, enabling formation of a relatively loose enhancer-enhancer interaction complex (EEIC), which overall increased interactions between that EEIC and central regions of HoxA chromatin. These findings indicate that CBE+47 organizes chromatin interactions between its adjacent enhancers and HoxA. Furthermore, deletion of those adjacent enhancers synergistically inhibited HoxA activation, suggesting that these enhancers serve as an EEIC required for RA-induced HoxA activation. Collectively, these results provide new insight into RA-induced HoxA expression during early ESC differentiation, also highlight precise regulatory roles of the CTCF-binding element in orchestrating high-order chromatin structure.

A distal enhancer maintaining Hoxa1 expression orchestrates retinoic acid-induced early ESCs differentiation.

Retinoic acid (RA) induces rapid differentiation of embryonic stem cells (ESCs), partly by activating expression of the transcription factor Hoxa1, which regulates downstream target genes that promote ESCs differentiation. However, mechanisms of RA-induced Hoxa1 expression and ESCs early differentiation remain largely unknown. Here, we identify a distal enhancer interacting with the Hoxa1 locus through a long-range chromatin loop. Enhancer deletion significantly inhibited expression of RA-induced Hoxa1 and endoderm master control genes such as Gata4 and Gata6. Transcriptome analysis revealed that RA-induced early ESCs differentiation was blocked in Hoxa1 enhancer knockout cells, suggesting a requirement for the enhancer. Restoration of Hoxa1 expression partly rescued expression levels of ∼40% of genes whose expression changed following enhancer deletion, and ∼18% of promoters of those rescued genes were directly bound by Hoxa1. Our data show that a distal enhancer maintains Hoxa1 expression through long-range chromatin loop and that Hoxa1 directly regulates downstream target genes expression and then orchestrates RA-induced early differentiation of ESCs. This discovery reveals mechanisms of a novel enhancer regulating RA-induced Hoxa genes expression and early ESCs differentiation.

Smek promotes corticogenesis through regulating Mbd3's stability and Mbd3/NuRD complex recruitment to genes associated with neurogenesis.

The fate of neural progenitor cells (NPCs) during corticogenesis is determined by a complex interplay of genetic or epigenetic components, but the underlying mechanism is incompletely understood. Here, we demonstrate that Suppressor of Mek null (Smek) interact with methyl-CpG-binding domain 3 (Mbd3) and the complex plays a critical role in self-renewal and neuronal differentiation of NPCs. We found that Smek promotes Mbd3 polyubiquitylation and degradation, blocking recruitment of the repressive Mbd3/nucleosome remodeling and deacetylase (NuRD) complex at the neurogenesis-associated gene loci, and, as a consequence, increasing acetyl histone H3 activity and cortical neurogenesis. Furthermore, overexpression of Mbd3 significantly blocked neuronal differentiation of NPCs, and Mbd3 depletion rescued neurogenesis defects seen in Smek1/2 knockout mice. These results reveal a novel molecular mechanism underlying Smek/Mbd3/NuRD axis-mediated control of NPCs' self-renewal and neuronal differentiation during mammalian corticogenesis.

HOTAIRM1, an enhancer lncRNA, promotes glioma proliferation by regulating long-range chromatin interactions within HOXA cluster genes.

The long noncoding RNA HOTAIRM1 reportedly plays important roles in acute myeloid leukemia, gastric cancer and colorectal cancer. Here, we analyzed potential function of HOTAIRM1 in glioma and asked whether it participates in long-range chromatin interactions. We monitored expression of HOTAIRM1 in glioma tissues and correlated levels with patient survival using the TCGA dataset. HOTAIRM1 was highly expressed in glioma tissue, with high levels associated with shortened patient survival time. We then suppressed HOTAIRM1 activity in the human glioblastoma U251 line using CRISPR-cas9 to knock in a truncating polyA fragment. Reporter analysis of these and control cells confirmed that the HOTAIRM1 locus serves as an active enhancer. We then performed Capture-C analysis to identify target genes of that locus and applied RNA antisense purification to assess chromatin interactions between the HOTAIRM1 locus and HOXA cluster genes. HOTAIRM1 knockdown in glioma cells decreased proliferation and reduced expression of HOXA cluster genes. HOTAIRM1 regulates long-range interactions between the HOTAIRM1 locus and HOXA genes. Our work suggests a new mechanism by which HOTAIRM1 regulates glioma progression by regulating high-order chromatin structure and could suggest novel therapeutic targets to treat an intractable cancer.

Smek1/2 is a nuclear chaperone and cofactor for cleaved Wnt receptor Ryk, regulating cortical neurogenesis.

The receptor-like tyrosine kinase (Ryk), a Wnt receptor, is important for cell fate determination during corticogenesis. During neuronal differentiation, the Ryk intracellular domain (ICD) is cleaved. Cleavage of Ryk and nuclear translocation of Ryk-ICD are required for neuronal differentiation. However, the mechanism of translocation and how it regulates neuronal differentiation remain unclear. Here, we identified Smek1 and Smek2 as Ryk-ICD partners that regulate its nuclear localization and function together with Ryk-ICD in the nucleus through chromatin recruitment and gene transcription regulation. Smek1/2 double knockout mice displayed pronounced defects in the production of cortical neurons, especially interneurons, while the neural stem cell population increased. In addition, both Smek and Ryk-ICD bound to the Dlx1/2 intergenic regulator element and were involved in its transcriptional regulation. These findings demonstrate a mechanism of the Ryk signaling pathway in which Smek1/2 and Ryk-ICD work together to mediate neural cell fate during corticogenesis.

Valproic acid promotes the neuronal differentiation of spiral ganglion neural stem cells with robust axonal growth.

Hearing loss occurs with the loss of hair cells of the cochlea and subsequent degeneration of spiral ganglion neurons (SGNs). Regeneration of SGNs is a potentially promising therapeutic approach to hearing loss in addition to the use of a cochlear implant (CI), because this device stimulates SGNs directly to restore hearing bypassing the missing hair cells. The presence of SGN-neural stem cells (NSCs) has been reported in adult human and mice. These cells have the potential to become SGNs and thus represent a cellular foundation for regeneration therapies for hearing loss. Valproic acid (VPA) has been shown to influence the neural differentiation of NSCs through multiple signaling pathways involving glycogen synthase kinase3ß (GSK3ß). Our present study therefore aimed to modulate the neural differentiation potential of SGN-NSCs by treatment with VPA. We here report that a clinically relevant concentration of 1 mM VPA induced the differentiation of basic fibroblast growth factor (bFGF)-treated P1- and P14-SGN-NSCs into neuronal and glial cells, confirmed by neuronal marker (Tuj1 and MAP2) and glial cell marker (GFAP and S100ß) detection. VPA-treated cells also promoted much longer neurite outgrowth compared to differentiated cells cultured without bFGF. The effects of VPA on the regulation of differentiation may be related to the activation of the Wnt/ß-catenin signaling pathway, but not the inhibition of histone deacetylases (HDACs). We propose that VPA has the potential to convert SGN-NSCs into SGNs and thereby restore hearing when combined with a CI.

Circular RNA profile in gliomas revealed by identification tool UROBORUS.

Recent evidence suggests that many endogenous circular RNAs (circRNAs) may play roles in biological processes. However, the expression patterns and functions of circRNAs in human diseases are not well understood. Computationally identifying circRNAs from total RNA-seq data is a primary step in studying their expression pattern and biological roles. In this work, we have developed a computational pipeline named UROBORUS to detect circRNAs in total RNA-seq data. By applying UROBORUS to RNA-seq data from 46 gliomas and normal brain samples, we detected thousands of circRNAs supported by at least two read counts, followed by successful experimental validation on 24 circRNAs from the randomly selected 27 circRNAs. UROBORUS is an efficient tool that can detect circRNAs with low expression levels in total RNA-seq without RNase R treatment. The circRNAs expression profiling revealed more than 476 circular RNAs differentially expressed in control brain tissues and gliomas. Together with parental gene expression, we found that circRNA and its parental gene have diversified expression patterns in gliomas and control brain tissues. This study establishes an efficient and sensitive approach for predicting circRNAs using total RNA-seq data. The UROBORUS pipeline can be accessed freely for non-commercial purposes at http://uroborus.openbioinformatics.org/.

w4CSeq: software and web application to analyze 4C-seq data.

Circularized Chromosome Conformation Capture followed by deep sequencing (4C-Seq) is a powerful technique to identify genome-wide partners interacting with a pre-specified genomic locus. Here, we present a computational and statistical approach to analyze 4C-Seq data generated from both enzyme digestion and sonication fragmentation-based methods. We implemented a command line software tool and a web interface called w4CSeq, which takes in the raw 4C sequencing data (FASTQ files) as input, performs automated statistical analysis and presents results in a user-friendly manner. Besides providing users with the list of candidate interacting sites/regions, w4CSeq generates figures showing genome-wide distribution of interacting regions, and sketches the enrichment of key features such as TSSs, TTSs, CpG sites and DNA replication timing around 4C sites. AVAILABILITY AND IMPLEMENTATION: Users can establish their own web server by downloading source codes at https://github.com/WGLab/w4CSeq Additionally, a demo web server is available at http://w4cseq.wglab.org CONTACT: kaiwang@usc.edu or wangelu@usc.eduSupplementary information: Supplementary data are available at Bioinformatics online.

The Wnt receptor Ryk reduces neuronal and cell survival capacity by repressing FOXO activity during the early phases of mutant huntingtin pathogenicity.

The Wnt receptor Ryk is an evolutionary-conserved protein important during neuronal differentiation through several mechanisms, including γ-secretase cleavage and nuclear translocation of its intracellular domain (Ryk-ICD). Although the Wnt pathway may be neuroprotective, the role of Ryk in neurodegenerative disease remains unknown. We found that Ryk is up-regulated in neurons expressing mutant huntingtin (HTT) in several models of Huntington's disease (HD). Further investigation in Caenorhabditis elegans and mouse striatal cell models of HD provided a model in which the early-stage increase of Ryk promotes neuronal dysfunction by repressing the neuroprotective activity of the longevity-promoting factor FOXO through a noncanonical mechanism that implicates the Ryk-ICD fragment and its binding to the FOXO co-factor ß-catenin. The Ryk-ICD fragment suppressed neuroprotection by lin-18/Ryk loss-of-function in expanded-polyQ nematodes, repressed FOXO transcriptional activity, and abolished ß-catenin protection of mutant htt striatal cells against cell death vulnerability. Additionally, Ryk-ICD was increased in the nucleus of mutant htt cells, and reducing γ-secretase PS1 levels compensated for the cytotoxicity of full-length Ryk in these cells. These findings reveal that the Ryk-ICD pathway may impair FOXO protective activity in mutant polyglutamine neurons, suggesting that neurons are unable to efficiently maintain function and resist disease from the earliest phases of the pathogenic process in HD.

Comprehensive functional annotation of 77 prostate cancer risk loci.

Genome-wide association studies (GWAS) have revolutionized the field of cancer genetics, but the causal links between increased genetic risk and onset/progression of disease processes remain to be identified. Here we report the first step in such an endeavor for prostate cancer. We provide a comprehensive annotation of the 77 known risk loci, based upon highly correlated variants in biologically relevant chromatin annotations--we identified 727 such potentially functional SNPs. We also provide a detailed account of possible protein disruption, microRNA target sequence disruption and regulatory response element disruption of all correlated SNPs at r(2) ≥ 0.88%. 88% of the 727 SNPs fall within putative enhancers, and many alter critical residues in the response elements of transcription factors known to be involved in prostate biology. We define as risk enhancers those regions with enhancer chromatin biofeatures in prostate-derived cell lines with prostate-cancer correlated SNPs. To aid the identification of these enhancers, we performed genomewide ChIP-seq for H3K27-acetylation, a mark of actively engaged enhancers, as well as the transcription factor TCF7L2. We analyzed in depth three variants in risk enhancers, two of which show significantly altered androgen sensitivity in LNCaP cells. This includes rs4907792, that is in linkage disequilibrium (r(2) = 0.91) with an eQTL for NUDT11 (on the X chromosome) in prostate tissue, and rs10486567, the index SNP in intron 3 of the JAZF1 gene on chromosome 7. Rs4907792 is within a critical residue of a strong consensus androgen response element that is interrupted in the protective allele, resulting in a 56% decrease in its androgen sensitivity, whereas rs10486567 affects both NKX3-1 and FOXA-AR motifs where the risk allele results in a 39% increase in basal activity and a 28% fold-increase in androgen stimulated enhancer activity. Identification of such enhancer variants and their potential target genes represents a preliminary step in connecting risk to disease process.

Defective entry into mitosis 1 (Dim1) negatively regulates osteoclastogenesis by inhibiting the expression of nuclear factor of activated T-cells, cytoplasmic, calcineurin-dependent 1 (NFATc1).

Bone remodeling is a continuous process of osteoblastic bone formation and osteoclastic bone resorption to maintain normal bone mass. NFATc1 is the master regulator of osteoclastogensis and transcriptionally activated by c-Fos and NF-κB in response to receptor activator of NF-κB ligand (RANKL) treatment. Defective entry into mitosis 1 (Dim1) is a nuclear protein that is implicated in pre-mRNA splicing and cell cycle progression, but the possible role of Dim1 in regulating other cellular processes remains unknown. Here, we demonstrate that Dim1 attenuates RANKL-induced osteoclastogenesis by targeting NFATc1 signaling pathway. Expression levels of Dim1 and NFATc1 are significantly increased during the formation of multinucleated osteoclasts. RNAi-mediated knockdown of Dim1 markedly enhances the expression of NFATc1 and its target genes, leading to the increase of RANKL-induced osteoclastogenesis in bone marrow-derived macrophages. Conversely, ectopic expression of Dim1 decreases RANKL-induced osteoclast differentiation by silencing NFATc1 and its target genes, further linking Dim1 to the dynamic regulation of osteoclastogenesis. Consistent with this notion, ChIP and interaction studies show that Dim1 directly associates with c-Fos and prevents c-Fos from binding to the NFATc1 promoter, resulting in targeted inactivation of the NFATc1 gene. Therefore, our studies reveal an unrecognized role for Dim1 as a master modulator of osteoclast differentiation, as well as the molecular mechanism underlying its repressive action toward osteoclastogensis.

Regulation of induced pluripotent stem (iPS) cell induction by Wnt/ß-catenin signaling.

Wnt signaling has been implicated in promoting somatic cell reprogramming. However, its molecular mechanisms remain unknown. Here we report that Wnt/ß-catenin enhances iPSCs induction at the early stage of reprogramming. The augmented reprogramming induced by ß-catenin is not due to increased total cell population or activation of c-Myc. In addition, ß-catenin interacts with reprogramming factors Klf4, Oct4, and Sox2, further enhancing expression of pluripotency circuitry genes. These studies reveal novel mechanisms underlying the regulation of reprogramming somatic cells to pluripotency by Wnt/ß-catenin signaling.

Iron regulatory protein-1 protects against mitoferrin-1-deficient porphyria.

Mitochondrial iron is essential for the biosynthesis of heme and iron-sulfur ([Fe-S]) clusters in mammalian cells. In developing erythrocytes, iron is imported into the mitochondria by MFRN1 (mitoferrin-1, SLC25A37). Although loss of MFRN1 in zebrafish and mice leads to profound anemia, mutant animals showed no overt signs of porphyria, suggesting that mitochondrial iron deficiency does not result in an accumulation of protoporphyrins. Here, we developed a gene trap model to provide in vitro and in vivo evidence that iron regulatory protein-1 (IRP1) inhibits protoporphyrin accumulation. Mfrn1(+/gt);Irp1(-/-) erythroid cells exhibit a significant increase in protoporphyrin levels. IRP1 attenuates protoporphyrin biosynthesis by binding to the 5'-iron response element (IRE) of alas2 mRNA, inhibiting its translation. Ectopic expression of alas2 harboring a mutant IRE, preventing IRP1 binding, in Mfrn1(gt/gt) cells mimics Irp1 deficiency. Together, our data support a model whereby impaired mitochondrial [Fe-S] cluster biogenesis in Mfrn1(gt/gt) cells results in elevated IRP1 RNA-binding that attenuates ALAS2 mRNA translation and protoporphyrin accumulation.

The functional genetic link of NLGN4X knockdown and neurodevelopment in neural stem cells.

Genetic mutations in NLGN4X (neuroligin 4), including point mutations and copy number variants (CNVs), have been associated with susceptibility to autism spectrum disorders (ASDs). However, it is unclear how mutations in NLGN4X result in neurodevelopmental defects. Here, we used neural stem cells (NSCs) as in vitro models to explore the impacts of NLGN4X knockdown on neurodevelopment. Using two shRNAmir-based vectors targeting NLGN4X and one control shRNAmir vector, we modulated NLGN4X expression and differentiated these NSCs into mature neurons. We monitored the neurodevelopmental process at Weeks 0, 0.5, 1, 2, 4 and 6, based on morphological analysis and whole-genome gene expression profiling. At the cellular level, in NSCs with NLGN4X knockdown, we observed increasingly delayed neuronal development and compromised neurite formation, starting from Week 2 through Week 6 post differentiation. At the molecular level, we identified multiple pathways, such as neurogenesis, neuron differentiation and muscle development, which are increasingly disturbed in cells with NLGN4X knockdown. Notably, several postsynaptic genes, including DLG4, NLGN1 and NLGN3, also have decreased expression. Based on in vitro models, NLGN4X knockdown directly impacts neurodevelopmental process during the formation of neurons and their connections. Our functional genomics study highlights the utility of NSCs models in understanding the functional roles of CNVs in affecting neurodevelopment and conferring susceptibility to neurodevelopmental diseases.