Zic2 is an enhancer-binding factor required for embryonic stem cell specification

Epigenetic control of cell identities from epiblast to gastrulation

Epigenetic modifications of chromatin are essential for the establishment of cell identities during embryogenesis. Between embryonic days 3.5-7.5 of murine development, major cell lineage decisions are made that discriminate extraembryonic and embryonic tissues, and the embryonic primary germ layers are formed, thereby laying down the basic body plan. In this review, we cover the contribution of dynamic chromatin modifications by DNA methylation, changes of chromatin accessibility, and histone modifications, that in combination with transcription factors control gene expression programs of different cell types. We highlight the differences in regulation of enhancer and promoter marks and discuss their requirement in cell lineage specification. Importantly, in many cases, lineage-specific targeting of epigenetic modifiers is carried out by pioneer or master transcription factors, that in sum mediate the chromatin landscape and thereby control the transcription of cell-type-specific gene programs and thus, cell identities.

Disabling leading and lagging strand histone transmission results in parental histones loss and reduced cell plasticity and viability

In the process of DNA replication, the first steps in restoring the chromatin landscape involve parental histone recycling and new histone deposition. Disrupting histone recycling to either the leading or lagging strand induces asymmetric histone inheritance, affecting epigenome maintenance and cellular identity. However, the order and kinetics of these effects remain elusive. Here, we use inducible mutants to dissect the early and late consequences of impaired histone recycling. Simultaneous disruption of both leading (POLE4) and lagging strand (MCM2-2A) recycling pathways impairs the transmission of parental histones to newly synthesized DNA, releasing some parental histones to the soluble pool. Subsequently, H3K27me3 accumulates aberrantly during chromatin restoration in a manner preceding gene expression changes. Loss of histone inheritance and the ensuing chromatin restoration defects alter gene expression in embryonic stem cells and challenge differentiation programs and cell viability. Our findings demonstrate the importance of efficient transmission of histone-based information during DNA replication for maintaining chromatin landscapes, differentiation potential, and cellular viability.

The enhancer module of Integrator controls cell identity and early neural fate commitment

Lineage-specific transcription factors operate as master orchestrators of developmental processes by activating select cis-regulatory enhancers and proximal promoters. Direct DNA binding of transcription factors ultimately drives context-specific recruitment of the basal transcriptional machinery that comprises RNA polymerase II (RNAPII) and a host of polymerase-associated multiprotein complexes, including the metazoan-specific Integrator complex. Integrator is primarily known to modulate RNAPII processivity and to surveil RNA integrity across coding genes. Here we describe an enhancer module of Integrator that directs cell fate specification by promoting epigenetic changes and transcription factor binding at neural enhancers. Depletion of Integrator's INTS10 subunit upends neural traits and derails cells towards mesenchymal identity. Commissioning of neural enhancers relies on Integrator's enhancer module, which stabilizes SOX2 binding at chromatin upon exit from pluripotency. We propose that Integrator is a functional bridge between enhancers and promoters and a main driver of early development, providing new insight into a growing family of neurodevelopmental syndromes.

ZIC2 and ZIC3 promote SWI/SNF recruitment to safeguard progression towards human primed pluripotency

The primed epiblast acts as a transitional stage between the relatively homogeneous naïve epiblast and the gastrulating embryo. Its formation entails coordinated changes in regulatory circuits driven by transcription factors and epigenetic modifications. Using a multi-omic approach in human embryonic stem cell models across the spectrum of peri-implantation development, we demonstrate that the transcription factors ZIC2 and ZIC3 have overlapping but essential roles in opening primed-specific enhancers. Together, they are essential to facilitate progression to and maintain primed pluripotency. ZIC2/3 accomplish this by recruiting SWI/SNF to chromatin and loss of ZIC2/3 or degradation of SWI/SNF both prevent enhancer activation. Loss of ZIC2/3 also results in transcriptome changes consistent with perturbed Polycomb activity and a shift towards the expression of genes linked to differentiation towards the mesendoderm. Additionally, we find an intriguing dependency on the transcriptional machinery for sustained recruitment of ZIC2/3 over a subset of primed-hESC specific enhancers. Taken together, ZIC2 and ZIC3 regulate highly dynamic lineage-specific enhancers and collectively act as key regulators of human primed pluripotency.

Genome binding properties of Zic transcription factors underlie their changing functions during neuronal maturation

BACKGROUND

The Zic family of transcription factors (TFs) promote both proliferation and maturation of cerebellar granule neurons (CGNs), raising the question of how a single, constitutively expressed TF family can support distinct developmental processes. Here we use an integrative experimental and bioinformatic approach to discover the regulatory relationship between Zic TF binding and changing programs of gene transcription during postnatal CGN differentiation.

RESULTS

We first established a bioinformatic pipeline to integrate Zic ChIP-seq data from the developing mouse cerebellum with other genomic datasets from the same tissue. In newborn CGNs, Zic TF binding predominates at active enhancers that are co-bound by developmentally regulated TFs including Atoh1, whereas in mature CGNs, Zic TF binding consolidates toward promoters where it co-localizes with activity-regulated TFs. We then performed CUT&RUN-seq in differentiating CGNs to define both the time course of developmental shifts in Zic TF binding and their relationship to gene expression. Mapping Zic TF binding sites to genes using chromatin looping, we identified the set of Zic target genes that have altered expression in RNA-seq from Zic1 or Zic2 knockdown CGNs.

CONCLUSIONS

Our data show that Zic TFs are required for both induction and repression of distinct, developmentally regulated target genes through a mechanism that is largely independent of changes in Zic TF binding. We suggest that the differential collaboration of Zic TFs with other TF families underlies the shift in their biological functions across CGN development.

Dual-role transcription factors stabilize intermediate expression levels

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Alternative splicing of a chromatin modifier alters the transcriptional regulatory programs of stem cell maintenance and neuronal differentiation

Development of embryonic stem cells (ESCs) into neurons requires intricate regulation of transcription, splicing, and translation, but how these processes interconnect is not understood. We found that polypyrimidine tract binding protein 1 (PTBP1) controls splicing of DPF2, a subunit of BRG1/BRM-associated factor (BAF) chromatin remodeling complexes. Dpf2 exon 7 splicing is inhibited by PTBP1 to produce the DPF2-S isoform early in development. During neuronal differentiation, loss of PTBP1 allows exon 7 inclusion and DPF2-L expression. Different cellular phenotypes and gene expression programs were induced by these alternative DPF2 isoforms. We identified chromatin binding sites enriched for each DPF2 isoform, as well as sites bound by both. In ESC, DPF2-S preferential sites were bound by pluripotency factors. In neuronal progenitors, DPF2-S sites were bound by nuclear factor I (NFI), while DPF2-L sites were bound by CCCTC-binding factor (CTCF). DPF2-S sites exhibited enhancer modifications, while DPF2-L sites showed promoter modifications. Thus, alternative splicing redirects BAF complex targeting to impact chromatin organization during neuronal development.

Transcription factor ZIC2 regulates the tumorigenic phenotypes associated with both bulk and cancer stem cells in epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is the most lethal gynecologic malignancy in North America. Current therapeutic regimens are ineffective against advanced EOC. A better understanding of the molecular mechanisms that regulate the biology of EOC will be a critical step toward developing more efficacious therapies against EOC. Herein, we demonstrate that elevated expression of transcription factor ZIC2 was associated with lower survival of EOC patients. Knockout of endogenous ZIC2 in EOC cells attenuated the tumorigenic phenotypes associated with both bulk and cancer stem cells in vitro and in vivo, indicating a pro-tumorigenic role of ZIC2 in EOC. On the other hand, however, overexpression of ZIC2 in EOC cells that do not express endogenous ZIC2 promoted cell migration and sphere formation, but inhibited cell growth and colony formation in vitro and tumor growth in vivo, indicating that the role for ZIC2 in EOC is context dependent. Our transcriptomic analysis showed that ZIC2-regulated genes were involved in multiple biological processes and signaling pathways associated with tumor progression. In conclusion, our findings reveal a context-dependent role for ZIC2 in regulating tumorigenic phenotypes in EOC, providing evidence that ZIC2 can be a potential therapeutic target for EOCs that express a high level of ZIC2.

Joubert syndrome-derived induced pluripotent stem cells show altered neuronal differentiation in vitro

Joubert syndrome (JS) is a recessively inherited congenital ataxia characterized by hypotonia, psychomotor delay, abnormal ocular movements, intellectual disability, and a peculiar cerebellar and brainstem malformation, the "molar tooth sign." Over 40 causative genes have been reported, all encoding for proteins implicated in the structure or functioning of the primary cilium, a subcellular organelle widely present in embryonic and adult tissues. In this paper, we developed an in vitro neuronal differentiation model using patient-derived induced pluripotent stem cells (iPSCs), to evaluate possible neurodevelopmental defects in JS. To this end, iPSCs from four JS patients harboring mutations in distinct JS genes (AHI1, CPLANE1, TMEM67, and CC2D2A) were differentiated alongside healthy control cells to obtain mid-hindbrain precursors and cerebellar granule cells. Differentiation was monitored over 31 days through the detection of lineage-specific marker expression by qRT-PCR, immunofluorescence, and transcriptomics analysis. All JS patient-derived iPSCs, regardless of the mutant gene, showed a similar impairment to differentiate into mid-hindbrain and cerebellar granule cells when compared to healthy controls. In addition, analysis of primary cilium count and morphology showed notable ciliary defects in all differentiating JS patient-derived iPSCs compared to controls. These results confirm that patient-derived iPSCs are an accessible and relevant in vitro model to analyze cellular phenotypes connected to the presence of JS gene mutations in a neuronal context.

PRAMEL7 and CUL2 decrease NuRD stability to establish ground-state pluripotency

Pluripotency is established in E4.5 preimplantation epiblast. Embryonic stem cells (ESCs) represent the immortalization of pluripotency, however, their gene expression signature only partially resembles that of developmental ground-state. Induced PRAMEL7 expression, a protein highly expressed in the ICM but lowly expressed in ESCs, reprograms developmentally advanced ESC+serum into ground-state pluripotency by inducing a gene expression signature close to developmental ground-state. However, how PRAMEL7 reprograms gene expression remains elusive. Here we show that PRAMEL7 associates with Cullin2 (CUL2) and this interaction is required to establish ground-state gene expression. PRAMEL7 recruits CUL2 to chromatin and targets regulators of repressive chromatin, including the NuRD complex, for proteasomal degradation. PRAMEL7 antagonizes NuRD-mediated repression of genes implicated in pluripotency by decreasing NuRD stability and promoter association in a CUL2-dependent manner. Our data link proteasome degradation pathways to ground-state gene expression, offering insights to generate in vitro models to reproduce the in vivo ground-state pluripotency.

Genome binding properties of Zic transcription factors underlie their changing functions during neuronal maturation

Background

The Zic family of transcription factors (TFs) promote both proliferation and maturation of cerebellar granule neurons (CGNs), raising the question of how a single, constitutively expressed TF family can support distinct developmental processes. Here we use an integrative experimental and bioinformatic approach to discover the regulatory relationship between Zic TF binding and changing programs of gene transcription during CGN differentiation.

Results

We first established a bioinformatic pipeline to integrate Zic ChIP-seq data from the developing mouse cerebellum with other genomic datasets from the same tissue. In newborn CGNs, Zic TF binding predominates at active enhancers that are co-bound by developmentally-regulated TFs including Atoh1, whereas in mature CGNs, Zic TF binding consolidates toward promoters where it co-localizes with activity-regulated TFs. We then performed CUT&RUN-seq in differentiating CGNs to define both the time course of developmental shifts in Zic TF binding and their relationship to gene expression. Mapping Zic TF binding sites to genes using chromatin looping, we identified the set of Zic target genes that have altered expression in RNA-seq from Zic1 or Zic2 knockdown CGNs.

Conclusion

Our data show that Zic TFs are required for both induction and repression of distinct, developmentally regulated target genes through a mechanism that is largely independent of changes in Zic TF binding. We suggest that the differential collaboration of Zic TFs with other TF families underlies the shift in their biological functions across CGN development.

Ablation of OCT4 function in cattle embryos by double electroporation of CRISPR-Cas for DNA and RNA targeting (CRISPR-DART)

CRISPR-Cas ribonucleoproteins (RNPs) are important tools for gene editing in preimplantation embryos. However, the inefficient production of biallelic deletions in cattle zygotes has hindered mechanistic studies of gene function. In addition, the presence of maternal RNAs that support embryo development until embryonic genome activation may cause confounding phenotypes. Here, we aimed to improve the efficiency of biallelic deletions and deplete specific maternal RNAs in cattle zygotes using CRISPR-Cas editing technology. Two electroporation sessions with Cas9D10A RNPs targeting exon 1 and the promoter of OCT4 produced biallelic deletions in 91% of the embryos tested. In most cases, the deletions were longer than 1,000 nucleotides long. Electroporation of Cas13a RNPs prevents the production of the corresponding proteins. We electroporated Cas9D10A RNPs targeting exon 1, including the promoter region, of OCT4 in two sessions with inclusion of Cas13a RNPs targeting OCT4 mRNAs in the second session to ablate OCT4 function in cattle embryos. A lack of OCT4 resulted in embryos arresting development prior to blastocyst formation at a greater proportion (13%) than controls (31.6%, P < 0.001). The few embryos that developed past the morula stage did not form a normal inner cell mass. Transcriptome analysis of single blastocysts, confirmed to lack exon 1 and promoter region of OCT4, revealed a significant (False Discovery Rate, FDR < 0.1) reduction in transcript abundance of many genes functionally connected to stemness, including markers of pluripotency (CADHD1, DPPA4, GNL3, RRM2). The results confirm that OCT4 is a key regulator of genes that modulate pluripotency and is required to form a functional blastocyst in cattle.

Young LINE-1 transposon 5' UTRs marked by elongation factor ELL3 function as enhancers to regulate naïve pluripotency in embryonic stem cells

LINE-1s are the major clade of retrotransposons with autonomous retrotransposition activity. Despite the potential genotoxicity, LINE-1s are highly activated in early embryos. Here we show that a subset of young LINE-1s, L1Md_Ts, are marked by the RNA polymerase II elongation factor ELL3, and function as enhancers in mouse embryonic stem cells. ELL3 depletion dislodges the DNA hydroxymethylase TET1 and the co-repressor SIN3A from L1Md_Ts, but increases the enrichment of the Bromodomain protein BRD4, leading to loss of 5hmC, gain of H3K27ac, and upregulation of the L1Md_T nearby genes. Specifically, ELL3 occupies and represses the L1Md_T-based enhancer located within Akt3, which encodes a key regulator of AKT pathway. ELL3 is required for proper ERK activation and efficient shutdown of naïve pluripotency through inhibiting Akt3 during naïve-primed transition. Our study reveals that the enhancer function of a subset of young LINE-1s controlled by ELL3 in transcription regulation and mouse early embryo development.

The AT-hook is an evolutionarily conserved auto-regulatory domain of SWI/SNF required for cell lineage priming

The SWI/SNF ATP-dependent chromatin remodeler is a master regulator of the epigenome, controlling pluripotency and differentiation. Towards the C-terminus of the catalytic subunit of SWI/SNF is a motif called the AT-hook that is evolutionary conserved. The AT-hook is present in many chromatin modifiers and generally thought to help anchor them to DNA. We observe however that the AT-hook regulates the intrinsic DNA-stimulated ATPase activity aside from promoting SWI/SNF recruitment to DNA or nucleosomes by increasing the reaction velocity a factor of 13 with no accompanying change in substrate affinity (KM). The changes in ATP hydrolysis causes an equivalent change in nucleosome movement, confirming they are tightly coupled. The catalytic subunit's AT-hook is required in vivo for SWI/SNF remodeling activity in yeast and mouse embryonic stem cells. The AT-hook in SWI/SNF is required for transcription regulation and activation of stage-specific enhancers critical in cell lineage priming. Similarly, growth assays suggest the AT-hook is required in yeast SWI/SNF for activation of genes involved in amino acid biosynthesis and metabolizing ethanol. Our findings highlight the importance of studying SWI/SNF attenuation versus eliminating the catalytic subunit or completely shutting down its enzymatic activity.

Zfp296 knockout enhances chromatin accessibility and induces a unique state of pluripotency in embryonic stem cells

The Zfp296 gene encodes a zinc finger-type protein. Its expression is high in mouse embryonic stem cells (ESCs) but rapidly decreases following differentiation. Zfp296-knockout (KO) ESCs grew as flat colonies, which were reverted to rounded colonies by exogenous expression of Zfp296. KO ESCs could not form teratomas when transplanted into mice but could efficiently contribute to germline-competent chimeric mice following blastocyst injection. Transcriptome analysis revealed that Zfp296 deficiency up- and down-regulates a distinct group of genes, among which Dppa3, Otx2, and Pou3f1 were markedly downregulated. Chromatin immunoprecipitation sequencing demonstrated that ZFP296 binding is predominantly seen in the vicinity of the transcription start sites (TSSs) of a number of genes, and ZFP296 was suggested to negatively regulate transcription. Consistently, chromatin accessibility assay clearly showed that ZFP296 binding reduces the accessibility of the TSS regions of target genes. Zfp296-KO ESCs showed increased histone H3K9 di- and trimethylation. Co-immunoprecipitation analyses revealed interaction of ZFP296 with G9a and GLP. These results show that ZFP296 plays essential roles in maintaining the global epigenetic state of ESCs through multiple mechanisms including activation of Dppa3, attenuation of chromatin accessibility, and repression of H3K9 methylation, but that Zfp296-KO ESCs retain a unique state of pluripotency while lacking the teratoma-forming ability.

CHD-associated enhancers shape human cardiomyocyte lineage commitment

Enhancers orchestrate gene expression programs that drive multicellular development and lineage commitment. Thus, genetic variants at enhancers are thought to contribute to developmental diseases by altering cell fate commitment. However, while many variant-containing enhancers have been identified, studies to endogenously test the impact of these enhancers on lineage commitment have been lacking. We perform a single-cell CRISPRi screen to assess the endogenous roles of 25 enhancers and putative cardiac target genes implicated in genetic studies of congenital heart defects (CHDs). We identify 16 enhancers whose repression leads to deficient differentiation of human cardiomyocytes (CMs). A focused CRISPRi validation screen shows that repression of TBX5 enhancers delays the transcriptional switch from mid- to late-stage CM states. Endogenous genetic deletions of two TBX5 enhancers phenocopy epigenetic perturbations. Together, these results identify critical enhancers of cardiac development and suggest that misregulation of these enhancers could contribute to cardiac defects in human patients.

Butyrate Induces Modifications of the CTCF-Binding Landscape in Cattle Cells

Butyrate is produced in the rumen from microbial fermentation and is related to several functions, including cell differentiation and proliferation. Butyrate supplementation in calves can accelerate rumen development. DNA-protein interactions, such as the CCCTC-binding factor (CTCF), play essential roles in chromatin organization and gene expression regulation. Although CTCF-binding sites have been identified recently in cattle, a deeper characterization, including differentially CTCF-binding sites (DCBS), is vital for a better understanding of butyrate’s role in the chromatin landscape. This study aimed to identify CTCF-binding regions and DCBS under a butyrate-induced condition using ChIP-seq in bovine cells; 61,915 CTCF peaks were identified in the butyrate and 51,347 in the control. From these regions, 2265 DCBS were obtained for the butyrate vs. control comparison, comprising ~90% of induced sites. Most of the butyrate DCBS were in distal intergenic regions, showing a potential role as insulators. Gene ontology enrichment showed crucial terms for the induced DCBS, mainly related to cellular proliferation, cell adhesion, and growth regulation. Interestingly, the ECM-receptor interaction pathway was observed for the induced DCBS. Motif enrichment analysis further identified transcription factors, including CTCF, BORIS, TGIF2, and ZIC3. When DCBS was integrated with RNA-seq data, putative genes were identified for the repressed DCBS, including GATA4. Our study revealed promising candidate genes in bovine cells by a butyrate-induced condition that might be related to the regulation of rumen development, such as integrins, keratins, and collagens. These results provide a better understanding of the function of butyrate in cattle rumen development and chromatin landscape regulation.

Epigenetics as "conductor" in "orchestra" of pluripotent states

Pluripotent character is described as the potency of cells to differentiate into all three germ layers. The best example to reinstate the term lies in the context of embryonic stem cells (ESCs). Pluripotent ESC describes the in vitro status of those cells that originate during the complex process of embryogenesis. Pre-implantation to post-implantation development of embryo embrace cells with different levels of stemness. Currently, four states of pluripotency have been recognized, in the progressing order of "naïve," "poised," "formative," and "primed." Epigenetics act as the "conductor" in this "orchestra" of transition in pluripotent states. With a distinguishable gene expression profile, these four states associate with different epigenetic signatures, sometimes distinct while otherwise overlapping. The present review focuses on how epigenetic factors, including DNA methylation, bivalent chromatin, chromatin remodelers, chromatin/nuclear architecture, and microRNA, could dictate pluripotent states and their transition among themselves.

Dynamic Transcriptome Profiling Reveals LncRNA-Centred Regulatory Networks in the Modulation of Pluripotency

Long noncoding RNAs (lncRNAs) have emerged as vital regulators of gene expression during embryonic stem cell (ESC) self-renewal and differentiation. Here, we systemically analyzed the differentially regulated lncRNAs during ESC-derived cardiomyocyte (CM) differentiation. We established a perspicuous profile of lncRNA expression at four critical developmental stages and found that the differentially expressed lncRNAs were grouped into six distinct clusters. The cluster with specific expression in ESC enriches the largest number of lncRNAs. Investigation of lncRNA-protein interaction network revealed that they are not only controlled by classic key transcription factors, but also modulated by epigenetic and epitranscriptomic factors including N6-methyladenosine (m6A) effector machineries. A detailed inspection revealed that 28 out of 385 lncRNAs were modified by methylation as well as directly recruited by the nuclear m6A reader protein Ythdc1. Unlike other 27 non-coding transcripts, the ESC-specific lncRNA Gm2379, located in both nucleus and cytoplasm, becomes dramatically upregulated in response to the depletion of m6A or Ythdc1. Consistent with the role of m6A in cell fate regulation, depletion of Gm2379 results in dysregulated expressions of pluripotent genes and crucial genes required for the formation of three germ layers. Collectively, our study provides a foundation for understanding the dynamic regulation of lncRNA transcriptomes during ESC differentiation and identifies the interplay between epitranscriptomic modification and key lncRNAs in the regulation of cell fate decision.

The DNA dioxygenase Tet1 regulates H3K27 modification and embryonic stem cell biology independent of its catalytic activity

Tet enzymes (Tet1/2/3) oxidize 5-methylcytosine to promote DNA demethylation and partner with chromatin modifiers to regulate gene expression. Tet1 is highly expressed in embryonic stem cells (ESCs), but its enzymatic and non-enzymatic roles in gene regulation are not dissected. We have generated Tet1 catalytically inactive (Tet1m/m) and knockout (Tet1-/-) ESCs and mice to study these functions. Loss of Tet1, but not loss of its catalytic activity, caused aberrant upregulation of bivalent (H3K4me3+; H3K27me3+) developmental genes, leading to defects in differentiation. Wild-type and catalytic-mutant Tet1 occupied similar genomic loci which overlapped with H3K27 tri-methyltransferase PRC2 and the deacetylase complex Sin3a at promoters of bivalent genes and with the helicase Chd4 at active genes. Loss of Tet1, but not loss of its catalytic activity, impaired enrichment of PRC2 and Sin3a at bivalent promoters leading to reduced H3K27 trimethylation and deacetylation, respectively, in absence of any changes in DNA methylation. Tet1-/-, but not Tet1m/m, embryos expressed higher levels of Gata6 and were developmentally delayed. Thus, the critical functions of Tet1 in ESCs and early development are mediated through its non-catalytic roles in regulating H3K27 modifications to silence developmental genes, and are more important than its catalytic functions in DNA demethylation.

Transcriptional Control of Axon Guidance at Midline Structures

The development of the nervous system is a time-ordered and multi-stepped process that includes neurogenesis and neuronal specification, axonal navigation, and circuits assembly. During axonal navigation, the growth cone, a dynamic structure located at the tip of the axon, senses environmental signals that guide axons towards their final targets. The expression of a specific repertoire of receptors on the cell surface of the growth cone together with the activation of a set of intracellular transducing molecules, outlines the response of each axon to specific guidance cues. This collection of axon guidance molecules is defined by the transcriptome of the cell which, in turn, depends on transcriptional and epigenetic regulators that modify the structure and DNA accessibility to determine what genes will be expressed to elicit specific axonal behaviors. Studies focused on understanding how axons navigate intermediate targets, such as the floor plate of vertebrates or the mammalian optic chiasm, have largely contributed to our knowledge of how neurons wire together during development. In fact, investigations on axon navigation at these midline structures led to the identification of many of the currently known families of proteins that act as guidance cues and their corresponding receptors. Although the transcription factors and the regulatory mechanisms that control the expression of these molecules are not well understood, important advances have been made in recent years in this regard. Here we provide an updated overview on the current knowledge about the transcriptional control of axon guidance and the selection of trajectories at midline structures.

β-catenin perturbations control differentiation programs in mouse embryonic stem cells

The Wnt/β-catenin pathway is involved in development, cancer, and embryonic stem cell (ESC) maintenance; its dual role in stem cell self-renewal and differentiation is still controversial. Here, by applying an in vitro system enabling inducible gene expression control, we report that moderate induction of transcriptionally active exogenous β-catenin in β-catenin null mouse ESCs promotes epiblast-like cell (EpiLC) derivation in vitro. Instead, in wild-type cells, moderate chemical pre-activation of the Wnt/β-catenin pathway promotes EpiLC in vitro derivation. Finally, we suggest that moderate β-catenin levels in β-catenin null mouse ESCs favor early stem cell commitment toward mesoderm if the exogenous protein is induced only in the “ground state” of pluripotency condition, or endoderm if the induction is maintained during the differentiation. Overall, our results confirm previous findings about the role of β-catenin in pluripotency and differentiation, while indicating a role for its doses in promoting specific differentiation programs.

•

Moderate β-catenin levels promote EpiLCs derivation in vitro

•

Chemical pre-activation of the Wnt pathway enhances ESC-EpiLC transition

•

β-catenin overexpression tips the balance between mesoderm and endoderm

•

Cell fate is influenced by the extent of β-catenin induction

The Molecular Landscape of Medulloblastoma in Teenagers and Young Adults

Medulloblastoma (MB) is a childhood malignant brain tumour but also occurs in teenagers and young adults (TYA). Considering that MB is heterogeneous, this study aimed to define the molecular landscape of MBs in TYAs. We collated more than 2000 MB samples that included 287 TYA patients (13-24 years). We performed computational analyses consisting of genome-wide methylation and transcriptomic profiles and developed a prognostics model for the TYAs with MB. We identified that TYAs predominantly comprised of Group 4 (40%) and Sonic Hedgehog (SHH)-activated (33%) tumours, with Wingless-type (WNT, 17%) and Group 3 (10%) being less common. TYAs with SHH tumours displayed significantly more gene expression alterations, whereas no gene was detected in the Group 4 tumours. Across MB subgroups, we identified unique and shared sets of TYA-specific differentially methylated probes and DNA-binding motifs. Finally, a 22-gene signature stratified TYA patients into high- and low-risk groups, and the prognostic significance of these risk groups persisted in multivariable regression models (P = 0.001). This study is an important step toward delineating the molecular landscape of TYAs with MB. The emergence of novel genes and pathways may provide a basis for improved clinical management of TYA with MB.

ANANSE: an enhancer network-based computational approach for predicting key transcription factors in cell fate determination

Proper cell fate determination is largely orchestrated by complex gene regulatory networks centered around transcription factors. However, experimental elucidation of key transcription factors that drive cellular identity is currently often intractable. Here, we present ANANSE (ANalysis Algorithm for Networks Specified by Enhancers), a network-based method that exploits enhancer-encoded regulatory information to identify the key transcription factors in cell fate determination. As cell type-specific transcription factors predominantly bind to enhancers, we use regulatory networks based on enhancer properties to prioritize transcription factors. First, we predict genome-wide binding profiles of transcription factors in various cell types using enhancer activity and transcription factor binding motifs. Subsequently, applying these inferred binding profiles, we construct cell type-specific gene regulatory networks, and then predict key transcription factors controlling cell fate transitions using differential networks between cell types. This method outperforms existing approaches in correctly predicting major transcription factors previously identified to be sufficient for trans-differentiation. Finally, we apply ANANSE to define an atlas of key transcription factors in 18 normal human tissues. In conclusion, we present a ready-to-implement computational tool for efficient prediction of transcription factors in cell fate determination and to study transcription factor-mediated regulatory mechanisms. ANANSE is freely available at https://github.com/vanheeringen-lab/ANANSE.

Loss of KMT2C reprograms the epigenomic landscape in hPSCs resulting in NODAL overexpression and a failure of hemogenic endothelium specification

Germline or somatic variation in the family of KMT2 lysine methyltransferases have been associated with a variety of congenital disorders and cancers. Notably, KMT2A-fusions are prevalent in 70% of infant leukaemias but fail to phenocopy short latency leukaemogenesis in mammalian models, suggesting additional factors are necessary for transformation. Given the lack of additional somatic mutation, the role of epigenetic regulation in cell specification, and our prior results of germline KMT2C variation in infant leukaemia patients, we hypothesized that germline dysfunction of KMT2C altered haematopoietic specification. In isogenic KMT2C KO hPSCs, we found genome-wide differences in histone modifications at active and poised enhancers, leading to gene expression profiles akin to mesendoderm rather than mesoderm highlighted by a significant increase in NODAL expression and WNT inhibition, ultimately resulting in a lack of in vitro hemogenic endothelium specification. These unbiased multi-omic results provide new evidence for germline mechanisms increasing risk of early leukaemogenesis.

Downregulation of LINC00665 suppresses the progression of lung adenocarcinoma via regulating miR-181c-5p/ZIC2 axis

Long non-coding RNA (lncRNA) LINC00665 was demonstrated to be upregulated in lung adenocarcinoma (LUAD) and target miR-181c-5p. ZIC2, which is upregulated in LUAD, serves as a putative target of miR-181c-5p. In this study, we aimed to reveal whether LINC00665 regulates miR-181c-5p/ZIC2 axis to promote LUAD progression. The results showed that LINC00665, HOXA1, ZIC2, and HOXA11 levels were increased in LUAD tissues, while miR-181c-5p level was decreased when compared to the adjacent normal tissues. High expression levels of LINC00665, ZIC2, HOXA1 and HOXA11, and low expression of miR-181c-5p were closely linked to poor prognosis of LUAD patients. Knockdown of LINC00665 induced obvious inhibitions in cell viability, clone formation, invasion and tumorigenesis in LUAD cells, whereas miR-181c-5p downregulation significantly neutralized these effects. In addition, downregulation of ZIC2 obviously reversed the enhancements of cell viability, clone formation, invasion and tumorigenesis induced by miR-181c-5p knockdown. In summary, the present study reveals that silencing of LINC00665 suppresses LUAD progression through targeting miR-181c-5p/ZIC2 axis.

Accessible chromatin reveals regulatory mechanisms underlying cell fate decisions during early embryogenesis

This study was conducted to investigate epigenetic landscape across multiple species and identify transcription factors (TFs) and their roles in controlling cell fate decision events during early embryogenesis. We made a comprehensively joint-research of chromatin accessibility of five species during embryogenesis by integration of ATAC-seq and RNA-seq datasets. Regulatory roles of candidate early embryonic TFs were investigated. Widespread accessible chromatin in early embryos overlapped with putative cis-regulatory sequences. Sets of cell-fate-determining TFs were identified. YOX1, a key cell cycle regulator, were found to homologous to clusters of TFs that are involved in neuron and epidermal cell-fate determination. Our research provides an intriguing insight into evolution of cell-fate decision during early embryogenesis among organisms.

Integrative genomic analysis of early neurogenesis reveals a temporal genetic program for differentiation and specification of preplate and Cajal-Retzius neurons

Neurogenesis in the developing neocortex begins with the generation of the preplate, which consists of early-born neurons including Cajal-Retzius (CR) cells and subplate neurons. Here, utilizing the Ebf2-EGFP transgenic mouse in which EGFP initially labels the preplate neurons then persists in CR cells, we reveal the dynamic transcriptome profiles of early neurogenesis and CR cell differentiation. Genome-wide RNA-seq and ChIP-seq analyses at multiple early neurogenic stages have revealed the temporal gene expression dynamics of early neurogenesis and distinct histone modification patterns in early differentiating neurons. We have identified a new set of coding genes and lncRNAs involved in early neuronal differentiation and validated with functional assays in vitro and in vivo. In addition, at E15.5 when Ebf2-EGFP+ cells are mostly CR neurons, single-cell sequencing analysis of purified Ebf2-EGFP+ cells uncovers molecular heterogeneities in CR neurons, but without apparent clustering of cells with distinct regional origins. Along a pseudotemporal trajectory these cells are classified into three different developing states, revealing genetic cascades from early generic neuronal differentiation to late fate specification during the establishment of CR neuron identity and function. Our findings shed light on the molecular mechanisms governing the early differentiation steps during cortical development, especially CR neuron differentiation.

Neural stem cells and progenitor cells in the embryonic brain give rise to neurons following a precise temporal order after initial expansion. Early-born neurons including Cajal-Retzius (CR) cells and subplate neurons form the preplate in the developing cerebral cortex, then CR neurons occupy the layer 1, playing an important role in cortical histogenesis. The molecular mechanisms governing the early neuronal differentiation processes remain to be explored. Here, by genome-wide approaches including bulk RNA-seq, single-cell RNA-seq and ChIP-seq, we comprehensively characterized the temporal dynamic gene expression profile and epigenetic status at different stages during early cortical development and uncovered molecularly heterogeneous subpopulations within the CR cells. We revealed CR neuron signatures and cell type-specific histone modification patterns along early neuron specification. Using in vitro and in vivo assays, we identified novel lncRNAs as potential functional regulators in preplate differentiation and CR neuron identity establishment. Our study provides a comprehensive analysis of the genetic and epigenetic programs during neuronal differentiation and would help bring new insights into the early cortical neurogenesis process, particularly the differentiation of CR neurons.

Epigenomic landscape of enhancer elements during Hydra head organizer formation

BACKGROUND

Axis patterning during development is accompanied by large-scale gene expression changes. These are brought about by changes in the histone modifications leading to dynamic alterations in chromatin architecture. The cis regulatory DNA elements also play an important role towards modulating gene expression in a context-dependent manner. Hydra belongs to the phylum Cnidaria where the first asymmetry in the body plan was observed and the oral-aboral axis originated. Wnt signaling has been shown to determine the head organizer function in the basal metazoan Hydra.

RESULTS

To gain insights into the evolution of cis regulatory elements and associated chromatin signatures, we ectopically activated the Wnt signaling pathway in Hydra and monitored the genome-wide alterations in key histone modifications. Motif analysis of putative intergenic enhancer elements from Hydra revealed the conservation of bilaterian cis regulatory elements that play critical roles in development. Differentially regulated enhancer elements were identified upon ectopic activation of Wnt signaling and found to regulate many head organizer specific genes. Enhancer activity of many of the identified cis regulatory elements was confirmed by luciferase reporter assay. Quantitative chromatin immunoprecipitation analysis upon activation of Wnt signaling further confirmed the enrichment of H3K27ac on the enhancer elements of Hv_Wnt5a, Hv_Wnt11 and head organizer genes Hv_Bra1, CnGsc and Hv_Pitx1. Additionally, perturbation of the putative H3K27me3 eraser activity using a specific inhibitor affected the ectopic activation of Wnt signaling indicating the importance of the dynamic changes in the H3K27 modifications towards regulation of the genes involved in the head organizer activity.

CONCLUSIONS

The activation-associated histone marks H3K4me3, H3K27ac and H3K9ac mark chromatin in a similar manner as seen in bilaterians. We identified intergenic cis regulatory elements which harbor sites for key transcription factors involved in developmental processes. Differentially regulated enhancers exhibited motifs for many zinc-finger, T-box and ETS related TFs whose homologs have a head specific expression in Hydra and could be a part of the pioneer TF network in the patterning of the head. The ability to differentially modify the H3K27 residue is critical for the patterning of Hydra axis revealing a dynamic acetylation/methylation switch to regulate gene expression and chromatin architecture.

ZIC2 promotes cancer stem cell traits via up-regulating OCT4 expression in lung adenocarcinoma cells

Background: Accumulating evidence has revealed the importance of cancer stem cells (CSCs) in self-renewal and chemoresistance. Previous studies reported high expression of ZIC2 was closely associated with tumorigenesis and CSC traits. However, the role of ZIC2 as a crucial factor for regulating CSC properties in lung adenocarcinoma (LAC) remains elusive. Methods: RT-PCR and WB assay were employed to assess ZIC2 expression in 20 LAC tumor tissues and the matched non-cancerous tissues. The role of ZIC2 in LAC CSC were analyzed by evaluation of CSC-related markers expression and spheroid formation in vitro. Cisplatin and paclitaxel resistance capacities were evaluated by CCK8 assay, colony formation assay, and flow cytometry analysis. Subcutaneous NOD/SCID mice models were generated to assess in vivo CSC features. Results: High expression of ZIC2 was found in LAC tumor tissues and indicated a poor overall survival in LAC patients. ZIC2 upregulated an array of CSCs-related genes, including EpCAM, OCT4, SOX2, NANOG, C-Myc and Bmi-1. Knockdown of ZIC2 inhibited sphere-forming capacity and decreased cisplatin and paclitaxel resistance. However, overexpression of ZIC2 achieved opposite effects. Mechanically, ZIC2 acts upstream of OCT4 to promote its expression, resulting in enhancement of CSC traits in LAC. Conclusion: Our results demonstrated that ZIC2 was crucial for promoting CSC traits in LAC cells, and served as a potential biomarker for predicting prognosis. The ZIC2-OCT4 network will facilitate the evaluation of the potential therapeutic efficacy of chemotherapy and predict patient sensitivity to treatment.

Systematized reporter assays reveal ZIC protein regulatory abilities are Subclass-specific and dependent upon transcription factor binding site context

The ZIC proteins are a family of transcription regulators with a well-defined zinc finger DNA-binding domain and there is evidence that they elicit functional DNA binding at a ZIC DNA binding site. Little is known, however, regarding domains within ZIC proteins that confer trans-activation or -repression. To address this question, a new cell-based trans-activation assay system suitable for ZIC proteins in HEK293T cells was constructed. This identified two previously unannotated evolutionarily conserved regions of ZIC3 that are necessary for trans-activation. These domains are found in all Subclass A ZIC proteins, but not in the Subclass B proteins. Additionally, the Subclass B proteins fail to elicit functional binding at a multimerised ZIC DNA binding site. All ZIC proteins, however, exhibit functional binding when the ZIC DNA binding site is embedded in a multiple transcription factor locus derived from ZIC target genes in the mouse genome. This ability is due to several domains, some of which are found in all ZIC proteins, that exhibit context dependent trans-activation or -repression activity. This knowledge is valuable for assessing the likely pathogenicity of variant ZIC proteins associated with human disorders and for determining factors that influence functional transcription factor binding.

ZIC3 Controls the Transition from Naive to Primed Pluripotency

Embryonic stem cells (ESCs) must transition through a series of intermediate cell states before becoming terminally differentiated. Here, we investigated the early events in this transition by determining the changes in the open chromatin landscape as naive mouse ESCs transition to epiblast-like cells (EpiLCs). Motif enrichment analysis of the newly opening regions coupled with expression analysis identified ZIC3 as a potential regulator of this cell fate transition. Chromatin binding and genome-wide transcriptional profiling following Zic3 depletion confirmed ZIC3 as an important regulatory transcription factor, and among its targets are genes encoding a number of transcription factors. Among these is GRHL2, which acts through enhancer switching to maintain the expression of a subset of genes from the ESC state. Our data therefore place ZIC3 upstream of a set of pro-differentiation transcriptional regulators and provide an important advance in our understanding of the regulatory factors governing the early steps in ESC differentiation.

•

Transcription factor ZIC3 regulates gene expression during the ESC to EpiLC transition

•

Extensive changes occur in the open chromatin landscape as ESCs progress to EpiLCs

•

ZIC3 activates the expression of a network of transcription factors

•

ZIC3-activated genes in EpiLCs are upregulated in the post-implantation epiblast

MEDEA: analysis of transcription factor binding motifs in accessible chromatin

Deciphering the interplay between chromatin accessibility and transcription factor (TF) binding is fundamental to understanding transcriptional regulation, control of cellular states, and the establishment of new phenotypes. Recent genome-wide chromatin accessibility profiling studies have provided catalogs of putative open regions, where TFs can recognize their motifs and regulate gene expression programs. Here, we present motif enrichment in differential elements of accessibility (MEDEA), a computational tool that analyzes high-throughput chromatin accessibility genomic data to identify cell-type-specific accessible regions and lineage-specific motifs associated with TF binding therein. To benchmark MEDEA, we used a panel of reference cell lines profiled by ENCODE and curated by the ENCODE Project Consortium for the ENCODE-DREAM Challenge. By comparing results with RNA-seq data, ChIP-seq peaks, and DNase-seq footprints, we show that MEDEA improves the detection of motifs associated with known lineage specifiers. We then applied MEDEA to 610 ENCODE DNase-seq data sets, where it revealed significant motifs even when absolute enrichment was low and where it identified novel regulators, such as NRF1 in kidney development. Finally, we show that MEDEA performs well on both bulk and single-cell ATAC-seq data. MEDEA is publicly available as part of our Glossary-GENRE suite for motif enrichment analysis.

Zygotic pioneer factor activity of Odd-paired/Zic is necessary for late function of the Drosophila segmentation network

Because chromatin determines whether information encoded in DNA is accessible to transcription factors, dynamic chromatin states in development may constrain how gene regulatory networks impart embryonic pattern. To determine the interplay between chromatin states and regulatory network function, we performed ATAC-seq on Drosophila embryos during the establishment of the segmentation network, comparing wild-type and mutant embryos in which all graded maternal patterning inputs are eliminated. While during the period between zygotic genome activation and gastrulation many regions maintain stable accessibility, cis-regulatory modules (CRMs) within the network undergo extensive patterning-dependent changes in accessibility. A component of the network, Odd-paired (opa), is necessary for pioneering accessibility of late segmentation network CRMs. opa-driven changes in accessibility are accompanied by equivalent changes in gene expression. Interfering with the timing of opa activity impacts the proper patterning of expression. These results indicate that dynamic systems for chromatin regulation directly impact the reading of embryonic patterning information.

Coordination of germ layer lineage choice by TET1 during primed pluripotency

Here, Luo et al. show how the DNA dioxygenase Tet1 plays a pivotal role upstream of germ layer lineage bifurcation. A permissive role for Tet1 in neural fate induction is identified, and involves Zic2-dependent engagement at neural target genes at lineage priming, is dependent on the signaling environment during gastrulation, and impacts neural tube closure after gastrulation.

Gastrulation in the early postimplantation stage mammalian embryo begins when epiblast cells ingress to form the primitive streak or develop as the embryonic ectoderm. The DNA dioxygenase Tet1 is highly expressed in the epiblast and yet continues to regulate lineage specification during gastrulation when its expression is diminished. Here, we show how Tet1 plays a pivotal role upstream of germ layer lineage bifurcation. During the transition from naive pluripotency to lineage priming, a global reconfiguration redistributes Tet1 from Oct4-cobound promoters to distal regulatory elements at lineage differentiation genes, which are distinct from high-affinity sites engaged by Oct4. An altered chromatin landscape in Tet1-deficient primed epiblast-like cells is associated with enhanced Oct4 expression and binding to Nodal and Wnt target genes, resulting in collaborative signals that enhance mesendodermal and inhibit neuroectodermal gene expression during lineage segregation. A permissive role for Tet1 in neural fate induction involves Zic2-dependent engagement at neural target genes at lineage priming, is dependent on the signaling environment during gastrulation, and impacts neural tube closure after gastrulation. Our findings provide mechanistic information for epigenetic integration of pluripotency and signal-induced differentiation cues.

Establishment of chromatin accessibility by the conserved transcription factor Grainy head is developmentally regulated

The dramatic changes in gene expression required for development necessitate the establishment of cis-regulatory modules defined by regions of accessible chromatin. Pioneer transcription factors have the unique property of binding closed chromatin and facilitating the establishment of these accessible regions. Nonetheless, much of how pioneer transcription factors coordinate changes in chromatin accessibility during development remains unknown. To determine whether pioneer-factor function is intrinsic to the protein or whether pioneering activity is developmentally modulated, we studied the highly conserved, essential transcription factor Grainy head (Grh). Prior work established that Grh is expressed throughout Drosophila development and is a pioneer factor in the larva. We demonstrated that Grh remains bound to mitotic chromosomes, a property shared with other pioneer factors. By assaying chromatin accessibility in embryos lacking maternal and/or zygotic Grh at three stages of development, we discovered that Grh is not required for chromatin accessibility in early embryogenesis, in contrast to its essential functions later in development. Our data reveal that the pioneering activity of Grh is temporally regulated and likely influenced by additional factors expressed at a given developmental stage.

Genome-wide CRISPR screen identifies ZIC2 as an essential gene that controls the cell fate of early mesodermal precursors to human heart progenitors

Cardiac progenitor formation is one of the earliest committed steps of human cardiogenesis and requires the cooperation of multiple gene sets governed by developmental signaling cascades. To determine the key regulators for cardiac progenitor formation, we have developed a two‐stage genome‐wide CRISPR‐knockout screen. We mimicked the progenitor formation process by differentiating human pluripotent stem cells (hPSCs) into cardiomyocytes, monitored by two distinct stage markers of early cardiac mesodermal formation and commitment to a multipotent heart progenitor cell fate: MESP1 and ISL1, respectively. From the screen output, we compiled a list of 15 candidate genes. After validating seven of them, we identified ZIC2 as an essential gene for cardiac progenitor formation. ZIC2 is known as a master regulator of neurogenesis. hPSCs with ZIC2 mutated still express pluripotency markers. However, their ability to differentiate into cardiomyocytes was greatly attenuated. RNA‐Seq profiling of the ZIC2‐mutant cells revealed that the mutants switched their cell fate alternatively to the noncardiac cell lineage. Further, single cell RNA‐seq analysis showed the ZIC2 mutants affected the apelin receptor‐related signaling pathway during mesoderm formation. Our results provide a new link between ZIC2 and human cardiogenesis and document the potential power of a genome‐wide unbiased CRISPR‐knockout screen to identify the key steps in human mesoderm precursor cell‐ and heart progenitor cell‐fate determination during in vitro hPSC cardiogenesis.

Zfp281 orchestrates interconversion of pluripotent states by engaging Ehmt1 and Zic2

Developmental cell fate specification is a unidirectional process that can be reverted in response to injury or experimental reprogramming. Whether differentiation and de-differentiation trajectories intersect mechanistically is unclear. Here, we performed comparative screening in lineage-related mouse naïve embryonic stem cells (ESCs) and primed epiblast stem cells (EpiSCs), and identified the constitutively expressed zinc finger transcription factor (TF) Zfp281 as a bidirectional regulator of cell state interconversion. We showed that subtle chromatin binding changes in differentiated cells translate into activation of the histone H3 lysine 9 (H3K9) methyltransferase Ehmt1 and stabilization of the zinc finger TF Zic2 at enhancers and promoters. Genetic gain-of-function and loss-of-function experiments confirmed a critical role of Ehmt1 and Zic2 downstream of Zfp281 both in driving exit from the ESC state and in restricting reprogramming of EpiSCs. Our study reveals that cell type-invariant chromatin association of Zfp281 provides an interaction platform for remodeling the cis-regulatory network underlying cellular plasticity.

Drosophila Zic family member odd-paired is needed for adult post-ecdysis maturation

Specific neuropeptides regulate in arthropods the shedding of the old cuticle (ecdysis) followed by maturation of the new cuticle. In Drosophila melanogaster, the last ecdysis occurs at eclosion from the pupal case, with a post-eclosion behavioural sequence that leads to wing extension, cuticle stretching and tanning. These events are highly stereotyped and are controlled by a subset of crustacean cardioactive peptide (CCAP) neurons through the expression of the neuropeptide Bursicon (Burs). We have studied the role of the transcription factor Odd-paired (Opa) during the post-eclosion period. We report that opa is expressed in the CCAP neurons of the central nervous system during various steps of the ecdysis process and in peripheral CCAP neurons innerving the larval muscles involved in adult ecdysis. We show that its downregulation alters Burs expression in the CCAP neurons. Ectopic expression of Opa, or the vertebrate homologue Zic2, in the CCAP neurons also affects Burs expression, indicating an evolutionary functional conservation. Finally, our results show that, independently of its role in Burs regulation, Opa prevents death of CCAP neurons during larval development.

Inactivation of the Wnt/β-catenin signaling pathway underlies inhibitory role of microRNA-129-5p in epithelial-mesenchymal transition and angiogenesis of prostate cancer by targeting ZIC2

Background

Prostate cancer (PCa) is a common disease that often occurs among older men and a frequent cause of malignancy associated death in this group. microRNA (miR)-129-5p has been identified as an essential regulator with a significant role in the prognosis of PC. Therefore, this study aimed to investigate roles of miR-129-5p in PCa.

Methods

Microarray analysis was conducted to identify PCa-related genes. The expression of miR-129-5p and ZIC2 in PCa tissues was investigated. To understand the role of miR-129-5p and ZIC2 in PCa, DU145 cells were transfected with mimic or inhibitor of miR-129-5p, or si-ZIC2 and the expression of Wnt, β-catenin, E-cadherin, vimentin, N-cadherin, vascular endothelial growth factor (VEGF), and CD31, as well as the extent of β-catenin phosphorylation was determined. In addition, cell proliferation, migration, invasion, angiogenesis, apoptosis and tumorigenesis were detected.

Results

miR-129-5p was poorly expressed and ZIC2 was highly expressed in PCa tissues. Down-regulation of ZIC2 or overexpression of miR-129-5p reduced the expression of ZIC2, Wnt, β-catenin, N-cadherin, vimentin, and β-catenin phosphorylation but increased the expression of E-cadherin. Importantly, miR-129-5p overexpression significantly reduced cell migration, invasion, angiogenesis and tumorigenesis while increasing cell apoptosis.

Conclusions

The findings of the present study indicated that overexpression of miR-129-5p or silencing of ZIC2 could inhibit epithelial–mesenchymal transition (EMT) and angiogenesis in PCa through blockage of the Wnt/β-catenin signaling pathway.

A sparse differential clustering algorithm for tracing cell type changes via single-cell RNA-sequencing data

Cell types in cell populations change as the condition changes: some cell types die out, new cell types may emerge and surviving cell types evolve to adapt to the new condition. Using single-cell RNA-sequencing data that measure the gene expression of cells before and after the condition change, we propose an algorithm, SparseDC, which identifies cell types, traces their changes across conditions and identifies genes which are marker genes for these changes. By solving a unified optimization problem, SparseDC completes all three tasks simultaneously. SparseDC is highly computationally efficient and demonstrates its accuracy on both simulated and real data.

β-Catenin/Tcf7l2-dependent transcriptional regulation of GLUT1 gene expression by Zic family proteins in colon cancer

The zinc finger of the cerebellum (ZIC) proteins has been implicated to function in normal tissue development. Recent studies have described the critical functions of Zic proteins in cancers and the potential tumor-suppressive functions in colon cancer development and progression. To elucidate the functional roles of Zic proteins in colorectal cancer, we knocked out the Zic5 gene and analyzed the chromatin localization pattern and transcriptional regulation of target gene expression. We found that Zic5 regulates glucose metabolism, and Zic5 knockout is accompanied by an increased glycolytic state and tolerance to a low-glucose condition. Furthermore, loss of β-catenin or TCF7l2 diminishes the chromatin binding of Zic5 globally. Our studies suggest that the Wnt/β-catenin signaling pathway has a strong influence on the function of Zic proteins and glucose metabolism in colorectal cancers through GLUT1. Interfering Wnt/-catenin-Zic5 axis-regulated aerobic glycolysis represents a potentially effective strategy to selectively target colon cancer cells.

CG dinucleotides enhance promoter activity independent of DNA methylation

Most mammalian RNA polymerase II initiation events occur at CpG islands, which are rich in CpGs and devoid of DNA methylation. Despite their relevance for gene regulation, it is unknown to what extent the CpG dinucleotide itself actually contributes to promoter activity. To address this question, we determined the transcriptional activity of a large number of chromosomally integrated promoter constructs and monitored binding of transcription factors assumed to play a role in CpG island activity. This revealed that CpG density significantly improves motif-based prediction of transcription factor binding. Our experiments also show that high CpG density alone is insufficient for transcriptional activity, yet results in increased transcriptional output when combined with particular transcription factor motifs. However, this CpG contribution to promoter activity is independent of DNA methyltransferase activity. Together, this refines our understanding of mammalian promoter regulation as it shows that high CpG density within CpG islands directly contributes to an environment permissive for full transcriptional activity.

The miR-96 and RARγ signaling axis governs androgen signaling and prostate cancer progression

Expression levels of retinoic acid receptor gamma (NR1B3/RARG, encodes RARγ) are commonly reduced in prostate cancer (PCa). Therefore, we sought to establish the cellular and gene regulatory consequences of reduced RARγ expression, and determine RARγ regulatory mechanisms. RARG shRNA approaches in non-malignant (RWPE-1 and HPr1-AR) and malignant (LNCaP) prostate models revealed that reducing RARγ levels, rather than adding exogenous retinoid ligand, had the greatest impact on prostate cell viability and gene expression. ChIP-Seq defined the RARγ cistrome, which was significantly enriched at active enhancers associated with AR binding sites. Reflecting a significant genomic role for RARγ to regulate androgen signaling, RARγ knockdown in HPr1-AR cells significantly regulated the magnitude of the AR transcriptome. RARγ downregulation was explained by increased miR-96 in PCa cell and mouse models, and TCGA PCa cohorts. Biochemical approaches confirmed that miR-96 directly regulated RARγ expression and function. Capture of the miR-96 targetome by biotin-miR-96 identified that RARγ and a number of RARγ interacting co-factors including TACC1 were all targeted by miR-96, and expression of these genes were prominently altered, positively and negatively, in the TCGA-PRAD cohort. Differential gene expression analyses between tumors in the TCGA-PRAD cohort with lower quartile expression levels of RARG and TACC1 and upper quartile miR-96, compared to the reverse, identified a gene network including several RARγ target genes (e.g., SOX15) that significantly associated with worse disease-free survival (hazard ratio 2.23, 95% CI 1.58 to 2.88, p = 0.015). In summary, miR-96 targets a RARγ network to govern AR signaling, PCa progression and disease outcome.

Interactome Rewiring Following Pharmacological Targeting of BET Bromodomains

Targeting bromodomains (BRDs) of the bromo-and-extra-terminal (BET) family offers opportunities for therapeutic intervention in cancer and other diseases. Here, we profile the interactomes of BRD2, BRD3, BRD4, and BRDT following treatment with the pan-BET BRD inhibitor JQ1, revealing broad rewiring of the interaction landscape, with three distinct classes of behavior for the 603 unique interactors identified. A group of proteins associate in a JQ1-sensitive manner with BET BRDs through canonical and new binding modes, while two classes of extra-terminal (ET)-domain binding motifs mediate acetylation-independent interactions. Last, we identify an unexpected increase in several interactions following JQ1 treatment that define negative functions for BRD3 in the regulation of rRNA synthesis and potentially RNAPII-dependent gene expression that result in decreased cell proliferation. Together, our data highlight the contributions of BET protein modules to their interactomes allowing for a better understanding of pharmacological rewiring in response to JQ1.

Zinc finger protein 32 promotes breast cancer stem cell-like properties through directly promoting GPER transcription

Breast cancer is one of the leading causes of death in women. Due to the existence of a small fraction of stem cell-like subpopulations, some breast cancer subtypes exhibit very high malignancy and resistance to multiple therapies. The underlying mechanisms of how these subtypes acquire stem cell-like properties and progress more aggressively remain largely unknown. Zinc finger protein 32 (ZNF32), a newly discovered transcription factor, has been reported to be associated with breast cancer progression. However, many questions remain about its target genes and its exact mechanisms in regulating stem cell-like properties and drug resistance. In the present study, we examined the relationship between ZNF32 and GPER, a membrane-associated estrogen receptor, and we addressed their roles in stemness regulation in human breast cancer cell lines. Our results showed that ZNF32 could induce expansion of stem cell-like subpopulations and increase drug resistance by upregulating GPER expression, in which ERK activation was also implicated. We also illustrated that ZNF32 induced GPER expression via a ZNF32 binding sequence located within the GPER promoter region. A correlation between ZNF32/GPER expression and increased tumor incidence and burden was observed in xenograft mouse models. We conclude that ZNF32 can engage GPER/ERK signalling and confer breast cancer stem cell-like properties, which may indicate poor prognosis of breast cancer patients. ZNF32 and GPER targeted therapies might provide new solutions for breast cancer treatment.

NuRD-interacting protein ZFP296 regulates genome-wide NuRD localization and differentiation of mouse embryonic stem cells

The nucleosome remodeling and deacetylase (NuRD) complex plays an important role in gene expression regulation, stem cell self-renewal, and lineage commitment. However, little is known about the dynamics of NuRD during cellular differentiation. Here, we study these dynamics using genome-wide profiling and quantitative interaction proteomics in mouse embryonic stem cells (ESCs) and neural progenitor cells (NPCs). We find that the genomic targets of NuRD are highly dynamic during differentiation, with most binding occurring at cell-type specific promoters and enhancers. We identify ZFP296 as an ESC-specific NuRD interactor that also interacts with the SIN3A complex. ChIP-sequencing in Zfp296 knockout (KO) ESCs reveals decreased NuRD binding both genome-wide and at ZFP296 binding sites, although this has little effect on the transcriptome. Nevertheless, Zfp296 KO ESCs exhibit delayed induction of lineage-specific markers upon differentiation to embryoid bodies. In summary, we identify an ESC-specific NuRD-interacting protein which regulates genome-wide NuRD binding and cellular differentiation.

The NuRD complex plays an important role in regulating lineage commitment and cell fate during early embryonic development. Here the authors present an integrative analysis of MBD3/NuRD composition and binding in mouse embryonic stem cells and neural progenitor cells, providing a molecular basis for genome-wide NuRD localization

A Requirement for Zic2 in the Regulation of Nodal Expression Underlies the Establishment of Left-Sided Identity

ZIC2 mutation is known to cause holoprosencephaly (HPE). A subset of ZIC2 HPE probands harbour cardiovascular and visceral anomalies suggestive of laterality defects. 3D-imaging of novel mouse Zic2 mutants uncovers, in addition to HPE, laterality defects in lungs, heart, vasculature and viscera. A strong bias towards right isomerism indicates a failure to establish left identity in the lateral plate mesoderm (LPM), a phenotype that cannot be explained simply by the defective ciliogenesis previously noted in Zic2 mutants. Gene expression analysis showed that the left-determining NODAL-dependent signalling cascade fails to be activated in the LPM, and that the expression of Nodal at the node, which normally triggers this event, is itself defective in these embryos. Analysis of ChiP-seq data, in vitro transcriptional assays and mutagenesis reveals a requirement for a low-affinity ZIC2 binding site for the activation of the Nodal enhancer HBE, which is normally active in node precursor cells. These data show that ZIC2 is required for correct Nodal expression at the node and suggest a model in which ZIC2 acts at different levels to establish LR asymmetry, promoting both the production of the signal that induces left side identity and the morphogenesis of the cilia that bias its distribution.

Genomic positional conservation identifies topological anchor point RNAs linked to developmental loci

BACKGROUND

The mammalian genome is transcribed into large numbers of long noncoding RNAs (lncRNAs), but the definition of functional lncRNA groups has proven difficult, partly due to their low sequence conservation and lack of identified shared properties. Here we consider promoter conservation and positional conservation as indicators of functional commonality.

RESULTS

We identify 665 conserved lncRNA promoters in mouse and human that are preserved in genomic position relative to orthologous coding genes. These positionally conserved lncRNA genes are primarily associated with developmental transcription factor loci with which they are coexpressed in a tissue-specific manner. Over half of positionally conserved RNAs in this set are linked to chromatin organization structures, overlapping binding sites for the CTCF chromatin organiser and located at chromatin loop anchor points and borders of topologically associating domains (TADs). We define these RNAs as topological anchor point RNAs (tapRNAs). Characterization of these noncoding RNAs and their associated coding genes shows that they are functionally connected: they regulate each other's expression and influence the metastatic phenotype of cancer cells in vitro in a similar fashion. Furthermore, we find that tapRNAs contain conserved sequence domains that are enriched in motifs for zinc finger domain-containing RNA-binding proteins and transcription factors, whose binding sites are found mutated in cancers.

CONCLUSIONS

This work leverages positional conservation to identify lncRNAs with potential importance in genome organization, development and disease. The evidence that many developmental transcription factors are physically and functionally connected to lncRNAs represents an exciting stepping-stone to further our understanding of genome regulation.