ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation

Integrated 3D genome, epigenome and transcriptome analyses reveal transcriptional coordination of circadian rhythm in rice

Photoperiods integrate with the circadian clock to coordinate gene expression rhythms and thus ensure plant fitness to the environment. Genome-wide characterization and comparison of rhythmic genes under different light conditions revealed delayed phase under constant darkness (DD) and reduced amplitude under constant light (LL) in rice. Interestingly, ChIP-seq and RNA-seq profiling of rhythmic genes exhibit synchronous circadian oscillation in H3K9ac modifications at their loci and long non-coding RNAs (lncRNAs) expression at proximal loci. To investigate how gene expression rhythm is regulated in rice, we profiled the open chromatin regions and transcription factor (TF) footprints by time-series ATAC-seq. Although open chromatin regions did not show circadian change, a significant number of TFs were identified to rhythmically associate with chromatin and drive gene expression in a time-dependent manner. Further transcriptional regulatory networks mapping uncovered significant correlation between core clock genes and transcription factors involved in light/temperature signaling. In situ Hi-C of ZT8-specific expressed genes displayed highly connected chromatin association at the same time, whereas this ZT8 chromatin connection network dissociates at ZT20, suggesting the circadian control of gene expression by dynamic spatial chromatin conformation. These findings together implicate the existence of a synchronization mechanism between circadian H3K9ac modifications, chromatin association of TF and gene expression, and provides insights into circadian dynamics of spatial chromatin conformation that associate with gene expression rhythms.

Evaluation of the determinants for improved pluripotency induction and maintenance by engineered SOX17

An engineered SOX17 variant with point mutations within its DNA binding domain termed SOX17FNV is a more potent pluripotency inducer than SOX2, yet the underlying mechanism remains unclear. Although wild-type SOX17 was incapable of inducing pluripotency, SOX17FNV outperformed SOX2 in mouse and human pluripotency reprogramming. In embryonic stem cells, SOX17FNV could replace SOX2 to maintain pluripotency despite considerable sequence differences and upregulated genes expressed in cleavage-stage embryos. Mechanistically, SOX17FNV co-bound OCT4 more cooperatively than SOX2 in the context of the canonical SoxOct DNA element. SOX2, SOX17, and SOX17FNV were all able to bind nucleosome core particles in vitro, which is a prerequisite for pioneer transcription factors. Experiments using purified proteins and in cellular contexts showed that SOX17 variants phase-separated more efficiently than SOX2, suggesting an enhanced ability to self-organise. Systematic deletion analyses showed that the N-terminus of SOX17FNV was dispensable for its reprogramming activity. However, the C-terminus encodes essential domains indicating multivalent interactions that drive transactivation and reprogramming. We defined a minimal SOX17FNV (miniSOX) that can support reprogramming with high activity, reducing the payload of reprogramming cassettes. This study uncovers the mechanisms behind SOX17FNV-induced pluripotency and establishes engineered SOX factors as powerful cell engineering tools.

A microwell platform for high-throughput longitudinal phenotyping and selective retrieval of organoids

Organoids are powerful experimental models for studying the ontogeny and progression of various diseases including cancer. Organoids are conventionally cultured in bulk using an extracellular matrix mimic. However, bulk-cultured organoids physically overlap, making it impossible to track the growth of individual organoids over time in high throughput. Moreover, local spatial variations in bulk matrix properties make it difficult to assess whether observed phenotypic heterogeneity between organoids results from intrinsic cell differences or differences in the microenvironment. Here, we developed a microwell-based method that enables high-throughput quantification of image-based parameters for organoids grown from single cells, which can further be retrieved from their microwells for molecular profiling. Coupled with a deep learning image-processing pipeline, we characterized phenotypic traits including growth rates, cellular movement, and apical-basal polarity in two CRISPR-engineered human gastric organoid models, identifying genomic changes associated with increased growth rate and changes in accessibility and expression correlated with apical-basal polarity. A record of this paper's transparent peer review process is included in the supplemental information.

Analysis of regulatory network modules in hundreds of human stem cell lines reveals complex epigenetic and genetic factors contribute to pluripotency state differences between subpopulations

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discovered 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which were highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlied the coordinated expression of genes in the GNMs. Epigenetic analyses revealed that regulatory networks underlying self-renewal and pluripotency have a surprising level of complexity. Genetic analyses identified thousands of regulatory variants that overlapped predicted transcription factor binding sites and were associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network were significantly enriched for regulatory variants with large effects, suggesting that they may play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work captures the coordinated activity of tens of thousands of regulatory elements in hiPSCs and bins these elements into discrete functionally characterized regulatory networks, shows that regulatory elements in pluripotency networks harbor variants with large effects, and provides a rich resource for future pluripotent stem cell research.

Common and divergent gene regulatory networks control injury-induced and developmental neurogenesis in zebrafish retina

Following acute retinal damage, zebrafish possess the ability to regenerate all neuronal subtypes. This regeneration requires Müller glia (MG) to reprogram and divide asymmetrically to produce a multipotent Müller glia-derived neuronal progenitor cell (MGPC). This raises three key questions. First, does loss of different retinal cell subtypes induce unique MG regeneration responses? Second, do MG reprogram to a developmental retinal progenitor cell state? And finally, to what extent does regeneration recapitulate retinal development? We examined these questions by performing single-nuclear and single-cell RNA-Seq and ATAC-Seq in both developing and regenerating retinas. While MG reprogram to a state similar to late-stage retinal progenitors in developing retinas, there are transcriptional differences between reprogrammed MG/MGPCs and late progenitors, as well as reprogrammed MG in outer and inner retinal damage models. Validation of candidate genes confirmed that loss of different subtypes induces differences in transcription factor gene expression and regeneration outcomes. This work identifies major differences between gene regulatory networks activated following the selective loss of different subtypes of retina neurons, as well as between retinal regeneration and development.

Fluid shear stress-modulated chromatin accessibility reveals the mechano-dependency of endothelial SMAD1/5-mediated gene transcription

Bone morphogenetic protein (BMP) signaling and fluid shear stress (FSS) mediate complementary functions in vascular homeostasis and disease development. It remains to be shown whether altered chromatin accessibility downstream of BMP and FSS offers a crosstalk level to explain changes in SMAD-dependent transcription. Here, we employed ATAC-seq to analyze arterial endothelial cells stimulated with BMP9 and/or FSS. We found that BMP9-sensitive regions harbor non-palindromic GC-rich SMAD-binding elements (GGCTCC) and 69.7% of these regions become BMP-insensitive in the presence of FSS. While GATA and KLF transcription factor (TF) motifs are unique to BMP9- and FSS-sensitive regions, respectively, SOX motifs are common to both. Finally, we show that both SOX(13/18) and GATA(2/3/6) family members are directly upregulated by SMAD1/5. These findings highlight the mechano-dependency of SMAD-signaling by a sequential mechanism of first elevated pioneer TF expression, allowing subsequent chromatin opening to eventually providing accessibility to novel SMAD binding sites.

Deciphering early human pancreas development at the single-cell level

Understanding pancreas development can provide clues for better treatments of pancreatic diseases. However, the molecular heterogeneity and developmental trajectory of the early human pancreas are poorly explored. Here, we performed large-scale single-cell RNA sequencing and single-cell assay for transposase accessible chromatin sequencing of human embryonic pancreas tissue obtained from first-trimester embryos. We unraveled the molecular heterogeneity, developmental trajectories and regulatory networks of the major cell types. The results reveal that dorsal pancreatic multipotent cells in humans exhibit different gene expression patterns than ventral multipotent cells. Pancreato-biliary progenitors that generate ventral multipotent cells in humans were identified. Notch and MAPK signals from mesenchymal cells regulate the differentiation of multipotent cells into trunk and duct cells. Notably, we identified endocrine progenitor subclusters with different differentiation potentials. Although the developmental trajectories are largely conserved between humans and mice, some distinct gene expression patterns have also been identified. Overall, we provide a comprehensive landscape of early human pancreas development to understand its lineage transitions and molecular complexity.

Single-cell multi-omics sequencing reveals the immunological disturbance underlying STAT3-V637M Hyper-IgE syndrome

Hyper-IgE syndrome (HIES) is a primary immunodeficiency characterized by, among others, the excessive production of IgE and repetitive bacterial/fungal infections. Mutations in STAT3, a transcription factor that orchestrates immune responses, may cause HIES, but the underlying mechanisms are not fully understood. Here, we used multi-omic approaches to comprehensively decipher the immune disturbance in a male HIES patient harboring STAT3-V637M. In his peripheral blood mononuclear cell (PBMC) we found significant clonal expansion of CD8 T cells (with increased CD8 subunits expression, potentially enhancing responsiveness to MHC I molecules), but not in his CD4 T cells and B cells. Although his B cells exhibited a higher potential in producing immunoglobulin, elevated SPIC binding might bias the products toward IgE isotype. Immune checkpoint inhibitors, including CTLA4, LAG3, were overexpressed in his PBMC-CD4 T cells, accompanied by reduced CD28 and IL6ST (gp130) expression. In his CD4 T cells, integrative analyses predicted upstream transcription factors (including ETV6, KLF13, and RORA) for LAG3, IL6ST, and CD28, respectively. The down-regulation of phagocytosis and nitric oxide synthesis-related genes in his PBMC-monocytes seem to be the culprit of his disseminated bacterial/fungal infection. Counterintuitively, in his PBMC we predicted increased STAT3 binding in both naïve and mature CD4 compartments, although this was not observed in most of his PBMC. In his bronchoalveolar lavage fluid (BALF), we found two macrophage subtypes with anti-bacterial properties, which were identified by CXCL8/S100A8/S100A9, or SOD2, respectively. Together, we described how the immune cell landscape was disturbed in STAT3-V637M HIES, providing a resource for further studies.

In vivo screening characterizes chromatin factor functions during normal and malignant hematopoiesis

Cellular differentiation requires extensive alterations in chromatin structure and function, which is elicited by the coordinated action of chromatin and transcription factors. By contrast with transcription factors, the roles of chromatin factors in differentiation have not been systematically characterized. Here, we combine bulk ex vivo and single-cell in vivo CRISPR screens to characterize the role of chromatin factor families in hematopoiesis. We uncover marked lineage specificities for 142 chromatin factors, revealing functional diversity among related chromatin factors (i.e. barrier-to-autointegration factor subcomplexes) as well as shared roles for unrelated repressive complexes that restrain excessive myeloid differentiation. Using epigenetic profiling, we identify functional interactions between lineage-determining transcription factors and several chromatin factors that explain their lineage dependencies. Studying chromatin factor functions in leukemia, we show that leukemia cells engage homeostatic chromatin factor functions to block differentiation, generating specific chromatin factor-transcription factor interactions that might be therapeutically targeted. Together, our work elucidates the lineage-determining properties of chromatin factors across normal and malignant hematopoiesis.

Low temperature-induced regulatory network rewiring via WRKY regulators during banana peel browning

Banana (Musa spp.) fruits, as typical tropical fruits, are cold sensitive, and lower temperatures can disrupt cellular compartmentalization and lead to severe browning. How tropical fruits respond to low temperature compared to the cold response mechanisms of model plants remains unknown. Here, we systematically characterized the changes in chromatin accessibility, histone modifications, distal cis-regulatory elements, transcription factor binding, and gene expression levels in banana peels in response to low temperature. Dynamic patterns of cold-induced transcripts were generally accompanied by concordant chromatin accessibility and histone modification changes. These upregulated genes were enriched for WRKY binding sites in their promoters and/or active enhancers. Compared to banana peel at room temperature, large amounts of banana WRKYs were specifically induced by cold and mediated enhancer-promoter interactions regulating critical browning pathways, including phospholipid degradation, oxidation, and cold tolerance. This hypothesis was supported by DNA affinity purification sequencing, luciferase reporter assays, and transient expression assay. Together, our findings highlight widespread transcriptional reprogramming via WRKYs during banana peel browning at low temperature and provide an extensive resource for studying gene regulation in tropical plants in response to cold stress, as well as potential targets for improving cold tolerance and shelf life of tropical fruits.

PRDM6 promotes medulloblastoma by repressing chromatin accessibility and altering gene expression

SNCAIP duplication may promote Group 4 medulloblastoma via induction of PRDM6, a poorly characterized member of the PRDF1 and RIZ1 homology domain-containing (PRDM) family of transcription factors. Here, we investigated the function of PRDM6 in human hindbrain neuroepithelial stem cells and tested PRDM6 as a driver of Group 4 medulloblastoma. We report that human PRDM6 localizes predominantly to the nucleus, where it causes widespread repression of chromatin accessibility and complex alterations of gene expression patterns. Genome-wide mapping of PRDM6 binding reveals that PRDM6 binds to chromatin regions marked by histone H3 lysine 27 trimethylation that are located within, or proximal to, genes. Moreover, we show that PRDM6 expression in neuroepithelial stem cells promotes medulloblastoma. Surprisingly, medulloblastomas derived from PRDM6-expressing neuroepithelial stem cells match human Group 3, but not Group 4, medulloblastoma. We conclude that PRDM6 expression has oncogenic potential but is insufficient to drive Group 4 medulloblastoma from neuroepithelial stem cells. We propose that both PRDM6 and additional factors, such as specific cell-of-origin features, are required for Group 4 medulloblastoma. Given the lack of PRDM6 expression in normal tissues and its oncogenic potential shown here, we suggest that PRDM6 inhibition may have therapeutic value in PRDM6-expressing medulloblastomas.

Unravelling human hematopoietic progenitor cell diversity through association with intrinsic regulatory factors

Hematopoietic stem and progenitor cell (HSPC) transplantation is an essential therapy for hematological conditions, but finer definitions of human HSPC subsets with associated function could enable better tuning of grafts and more routine, lower-risk application. To deeply phenotype HSPCs, following a screen of 328 antigens, we quantified 41 surface proteins and functional regulators on millions of CD34+ and CD34- cells, spanning four primary human hematopoietic tissues: bone marrow, mobilized peripheral blood, cord blood, and fetal liver. We propose more granular definitions of HSPC subsets and provide new, detailed differentiation trajectories of erythroid and myeloid lineages. These aspects of our revised human hematopoietic model were validated with corresponding epigenetic analysis and in vitro clonal differentiation assays. Overall, we demonstrate the utility of using molecular regulators as surrogates for cellular identity and functional potential, providing a framework for description, prospective isolation, and cross-tissue comparison of HSPCs in humans.

PTEN regulates hematopoietic lineage plasticity via PU.1-dependent chromatin accessibility

PTEN loss in fetal liver hematopoietic stem cells (HSCs) leads to alterations in myeloid, T-, and B-lineage potentials and T-lineage acute lymphoblastic leukemia (T-ALL) development. To explore the mechanism underlying PTEN-regulated hematopoietic lineage choices, we carry out integrated assay for transposase-accessible chromatin using sequencing (ATAC-seq), single-cell RNA-seq, and in vitro culture analyses using in vivo-isolated mouse pre-leukemic HSCs and progenitors. We find that PTEN loss alters chromatin accessibility of key lineage transcription factor (TF) binding sites at the prepro-B stage, corresponding to increased myeloid and T-lineage potentials and reduced B-lineage potential. Importantly, we find that PU.1 is an essential TF downstream of PTEN and that altering PU.1 levels can reprogram the chromatin accessibility landscape and myeloid, T-, and B-lineage potentials in Ptennull prepro-B cells. Our study discovers prepro-B as the key developmental stage underlying PTEN-regulated hematopoietic lineage choices and suggests a critical role of PU.1 in modulating the epigenetic state and lineage plasticity of prepro-B progenitors.

Common and divergent gene regulatory networks control injury-induced and developmental neurogenesis in zebrafish retina

Following acute retinal damage, zebrafish possess the ability to regenerate all neuronal subtypes. This regeneration requires Müller glia (MG) to reprogram and divide asymmetrically to produce a multipotent Müller glia-derived neuronal progenitor cell (MGPC). This raises three key questions. First, does loss of different retinal cell subtypes induce unique MG regeneration responses? Second, do MG reprogram to a developmental retinal progenitor cell state? And finally, to what extent does regeneration recapitulate retinal development? We examined these questions by performing single-nuclear and single-cell RNA-Seq and ATAC-Seq in both developing and regenerating retinas. While MG reprogram to a state similar to late-stage retinal progenitors in developing retinas, there are transcriptional differences between reprogrammed MG/MGPCs and late progenitors, as well as reprogrammed MG in outer and inner retinal damage models. Validation of candidate genes confirmed that loss of different subtypes induces differences in transcription factor gene expression and regeneration outcomes. This work identifies major differences between gene regulatory networks activated following the selective loss of different subtypes of retina neurons, as well as between retinal regeneration and 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.

Characterization of an eye field-like state during optic vesicle organoid development

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors. This crucial event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

FOXA2 controls the anti-oxidant response in FH-deficient cells

Hereditary leiomyomatosis and renal cell cancer (HLRCC) is a cancer syndrome caused by inactivating germline mutations in fumarate hydratase (FH) and subsequent accumulation of fumarate. Fumarate accumulation leads to profound epigenetic changes and the activation of an anti-oxidant response via nuclear translocation of the transcription factor NRF2. The extent to which chromatin remodeling shapes this anti-oxidant response is currently unknown. Here, we explored the effects of FH loss on the chromatin landscape to identify transcription factor networks involved in the remodeled chromatin landscape of FH-deficient cells. We identify FOXA2 as a key transcription factor that regulates anti-oxidant response genes and subsequent metabolic rewiring cooperating without direct interaction with the anti-oxidant regulator NRF2. The identification of FOXA2 as an anti-oxidant regulator provides additional insights into the molecular mechanisms behind cell responses to fumarate accumulation and potentially provides further avenues for therapeutic intervention for HLRCC.

Epigenomic mapping in B-cell acute lymphoblastic leukemia identifies transcriptional regulators and noncoding variants promoting distinct chromatin architectures

B-cell lineage acute lymphoblastic leukemia (B-ALL) is comprised of diverse molecular subtypes and while transcriptional and DNA methylation profiling of B-ALL subtypes has been extensively examined, the accompanying chromatin landscape is not well characterized for many subtypes. We therefore mapped chromatin accessibility using ATAC-seq for 10 B-ALL molecular subtypes in primary ALL cells from 154 patients. Comparisons with B-cell progenitors identified candidate B-ALL cell-of-origin and AP-1-associated cis-regulatory rewiring in B-ALL. Cis-regulatory rewiring promoted B-ALL-specific gene regulatory networks impacting oncogenic signaling pathways that perturb normal B-cell development. We also identified that over 20% of B-ALL accessible chromatin sites exhibit strong subtype enrichment, with transcription factor (TF) footprint profiling identifying candidate TFs that maintain subtype-specific chromatin architectures. Over 9000 inherited genetic variants were further uncovered that contribute to variability in chromatin accessibility among individual patient samples. Overall, our data suggest that distinct chromatin architectures are driven by diverse TFs and inherited genetic variants which promote unique gene regulatory networks that contribute to transcriptional differences among B-ALL subtypes.

Progesterone receptor mediates ovulatory transcription through RUNX transcription factor interactions and chromatin remodelling

Progesterone receptor (PGR) plays diverse roles in reproductive tissues and thus coordinates mammalian fertility. In the ovary, rapid acute induction of PGR is the key determinant of ovulation through transcriptional control of a unique set of genes that culminates in follicle rupture. However, the molecular mechanisms for this specialized PGR function in ovulation is poorly understood. We have assembled a detailed genomic profile of PGR action through combined ATAC-seq, RNA-seq and ChIP-seq analysis in wildtype and isoform-specific PGR null mice. We demonstrate that stimulating ovulation rapidly reprograms chromatin accessibility in two-thirds of sites, correlating with altered gene expression. An ovary-specific PGR action involving interaction with RUNX transcription factors was observed with 70% of PGR-bound regions also bound by RUNX1. These transcriptional complexes direct PGR binding to proximal promoter regions. Additionally, direct PGR binding to the canonical NR3C motif enable chromatin accessibility. Together these PGR actions mediate induction of essential ovulatory genes. Our findings highlight a novel PGR transcriptional mechanism specific to ovulation, providing new targets for infertility treatments or new contraceptives that block ovulation.

Single-cell chromatin accessibility of developing murine pancreas identifies cell state-specific gene regulatory programs

Numerous studies have characterized the existence of cell subtypes, along with their corresponding transcriptional profiles, within the developing mouse pancreas. The upstream mechanisms that initiate and maintain gene expression programs across cell states, however, remain largely unknown. Here, we generate single-nucleus ATAC-Sequencing data of developing murine pancreas and perform an integrated, multi-omic analysis of both chromatin accessibility and RNA expression to describe the chromatin landscape of the developing pancreas at both E14.5 and E17.5 at single-cell resolution. We identify candidate transcription factors regulating cell fate and construct gene regulatory networks of active transcription factor binding to regulatory regions of downstream target genes. This work serves as a valuable resource for the field of pancreatic biology in general and contributes to our understanding of lineage plasticity among endocrine cell types. In addition, these data identify which epigenetic states should be represented in the differentiation of stem cells to the pancreatic beta cell fate to best recapitulate in vitro the gene regulatory networks that are critical for progression along the beta cell lineage in vivo.

ETV2 primes hematoendothelial gene enhancers prior to hematoendothelial fate commitment

Mechanisms underlying distinct specification, commitment, and differentiation phases of cell fate determination remain undefined due to difficulties capturing these processes. Here, we interrogate the activity of ETV2, a transcription factor necessary and sufficient for hematoendothelial differentiation, within isolated fate intermediates. We observe transcriptional upregulation of Etv2 and opening of ETV2-binding sites, indicating new ETV2 binding, in a common cardiac-hematoendothelial progenitor population. Accessible ETV2-binding sites are active at the Etv2 locus but not at other hematoendothelial regulator genes. Hematoendothelial commitment coincides with the activation of a small repertoire of previously accessible ETV2-binding sites at hematoendothelial regulators. Hematoendothelial differentiation accompanies activation of a large repertoire of new ETV2-binding sites and upregulation of hematopoietic and endothelial gene regulatory networks. This work distinguishes specification, commitment, and sublineage differentiation phases of ETV2-dependent transcription and suggests that the shift from ETV2 binding to ETV2-bound enhancer activation, not ETV2 binding to target enhancers, drives hematoendothelial fate commitment.

Chi3l1 Is a Modulator of Glioma Stem Cell States and a Therapeutic Target in Glioblastoma

Chitinase 3-like 1 (Chi3l1) is a secreted protein that is highly expressed in glioblastoma. Here, we show that Chi3l1 alters the state of glioma stem cells (GSC) to support tumor growth. Exposure of patient-derived GSCs to Chi3l1 reduced the frequency of CD133+SOX2+ cells and increased the CD44+Chi3l1+ cells. Chi3l1 bound to CD44 and induced phosphorylation and nuclear translocation of β-catenin, Akt, and STAT3. Single-cell RNA sequencing and RNA velocity following incubation of GSCs with Chi3l1 showed significant changes in GSC state dynamics driving GSCs towards a mesenchymal expression profile and reducing transition probabilities towards terminal cellular states. ATAC-seq revealed that Chi3l1 increases accessibility of promoters containing a Myc-associated zinc finger protein (MAZ) transcription factor footprint. Inhibition of MAZ downregulated a set of genes with high expression in cellular clusters that exhibit significant cell state transitions after treatment with Chi3l1, and MAZ deficiency rescued the Chi3L-induced increase of GSC self-renewal. Finally, targeting Chi3l1 in vivo with a blocking antibody inhibited tumor growth and increased the probability of survival. Overall, this work suggests that Chi3l1 interacts with CD44 on the surface of GSCs to induce Akt/β-catenin signaling and MAZ transcriptional activity, which in turn upregulates CD44 expression in a pro-mesenchymal feed-forward loop. The role of Chi3l1 in regulating cellular plasticity confers a targetable vulnerability to glioblastoma.

SIGNIFICANCE

Chi3l1 is a modulator of glioma stem cell states that can be targeted to promote differentiation and suppress growth of glioblastoma.

Characterization and Optimization of Multiomic Single-Cell Epigenomic Profiling

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n = 18) and a solid tumor type, ovarian cancer (OC, n = 18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell type identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n = 6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.

Chromatin accessibility dynamics insight into crosstalk between regulatory landscapes in poplar responses to multiple treatments

Perennial trees develop and coordinate endogenous response signaling pathways, including their crosstalk and convergence, to cope with various environmental stresses which occur simultaneously in most cases. These processes are involved in gene transcriptional regulations that depend on dynamic interactions between regulatory proteins and corresponding chromatin regions, but the mechanisms remain poorly understood in trees. In this study, we detected chromatin regulatory landscapes of poplar under abscisic acid, methyl jasmonate, salicylic acid and sodium chloride (NaCl) treatment, through integrating ATAC-seq and RNA-seq data. Our results showed that the degree of chromatin accessibility for a given gene is closely related to its expression level. However, unlike the gene expression that shows treatment-specific response patterns, changes in chromatin accessibility exhibit high similarities under these treatments. We further proposed and experimentally validated that a homologous gene copy of RESPONSIVE TO DESICCATION 26 mediates the crosstalk between jasmonic acid and NaCl signaling pathways by directly regulating the stress-responsive genes and that circadian clock-related transcription factors like REVEILLE8 play a central role in response of poplar to these treatments. Overall, our study provides a chromatin insight into the molecular mechanism of transcription regulatory networks in response to different environmental stresses and raises the key roles of the circadian clock of poplar to adapt to adverse environments.

LIS1 RNA-binding orchestrates the mechanosensitive properties of embryonic stem cells in AGO2-dependent and independent ways

Lissencephaly-1 (LIS1) is associated with neurodevelopmental diseases and is known to regulate the molecular motor cytoplasmic dynein activity. Here we show that LIS1 is essential for the viability of mouse embryonic stem cells (mESCs), and it governs the physical properties of these cells. LIS1 dosage substantially affects gene expression, and we uncovered an unexpected interaction of LIS1 with RNA and RNA-binding proteins, most prominently the Argonaute complex. We demonstrate that LIS1 overexpression partially rescued the extracellular matrix (ECM) expression and mechanosensitive genes conferring stiffness to Argonaute null mESCs. Collectively, our data transforms the current perspective on the roles of LIS1 in post-transcriptional regulation underlying development and mechanosensitive processes.

Correction of transposase sequence bias in ATAC-seq data with rule ensemble modeling

Chromatin accessibility assays have revolutionized the field of transcription regulation by providing single-nucleotide resolution measurements of regulatory features such as promoters and transcription factor binding sites. ATAC-seq directly measures how well the Tn5 transposase accesses chromatinized DNA. Tn5 has a complex sequence bias that is not effectively scaled with traditional bias-correction methods. We model this complex bias using a rule ensemble machine learning approach that integrates information from many input k-mers proximal to the ATAC sequence reads. We effectively characterize and correct single-nucleotide sequence biases and regional sequence biases of the Tn5 enzyme. Correction of enzymatic sequence bias is an important step in interpreting chromatin accessibility assays that aim to infer transcription factor binding and regulatory activity of elements in the genome.

Epigenetic programming defines haematopoietic stem cell fate restriction

Haematopoietic stem cells (HSCs) are multipotent, but individual HSCs can show restricted lineage output in vivo. Currently, the molecular mechanisms and physiological role of HSC fate restriction remain unknown. Here we show that lymphoid fate is epigenetically but not transcriptionally primed in HSCs. In multi-lineage HSCs that produce lymphocytes, lymphoid-specific upstream regulatory elements (LymUREs) but not promoters are preferentially accessible compared with platelet-biased HSCs that do not produce lymphoid cell types, providing transcriptionally silent lymphoid lineage priming. Runx3 is preferentially expressed in multi-lineage HSCs, and reinstating Runx3 expression increases LymURE accessibility and lymphoid-primed multipotent progenitor 4 (MPP4) output in old, platelet-biased HSCs. In contrast, platelet-biased HSCs show elevated levels of epigenetic platelet-lineage priming and give rise to MPP2 progenitors with molecular platelet bias. These MPP2 progenitors generate platelets with faster kinetics and through a more direct cellular pathway compared with MPP2s derived from multi-lineage HSCs. Epigenetic programming therefore predicts both fate restriction and differentiation kinetics in HSCs.

Stem cells tightly regulate dead cell clearance to maintain tissue fitness

Macrophages and dendritic cells have long been appreciated for their ability to migrate to and engulf dying cells and debris, including some of the billions of cells that are naturally eliminated from our body daily. However, a substantial number of these dying cells are cleared by 'non-professional phagocytes', local epithelial cells that are critical to organismal fitness. How non-professional phagocytes sense and digest nearby apoptotic corpses while still performing their normal tissue functions is unclear. Here, we explore the molecular mechanisms underlying their multifunctionality. Exploiting the cyclical bouts of tissue regeneration and degeneration during the hair cycle, we show that stem cells can transiently become non-professional phagocytes when confronted with dying cells. Adoption of this phagocytic state requires both local lipids produced by apoptotic corpses to activate RXRα, and tissue-specific retinoids for RARγ activation. This dual factor dependency enables tight regulation of the genes requisite to activate phagocytic apoptotic clearance. The tunable phagocytic program we describe here offers an effective mechanism to offset phagocytic duties against the primary stem cell function of replenishing differentiated cells to preserve tissue integrity during homeostasis. Our findings have broad implications for other non-motile stem or progenitor cells which experience cell death in an immune-privileged niche.

A wound-induced differentiation trajectory for neurons

Animals capable of whole-body regeneration can replace any missing cell type and regenerate fully-functional new organs, de novo . The regeneration of a new brain requires the formation of diverse neuronal cell types and their assembly into an organized structure and correctly-wired circuits. Recent work in various regenerative animals has revealed transcriptional programs required for the differentiation of distinct neuronal subpopulations, however how these transcriptional programs are initiated upon amputation remains unknown. Here, we focused on the highly regenerative acoel worm, Hofstenia miamia , to study wound-induced transcriptional regulatory events that lead to the production of neurons. Footprinting analysis using chromatin accessibility data on an improved genome assembly revealed that binding sites for the NFY transcription factor complex were significantly bound during regeneration, showing a dynamic increase in binding within one hour upon amputation specifically in tail fragments, which will regenerate a new brain. Strikingly, NFY targets were highly enriched for genes with neuronal functional. Single-cell transcriptome analysis combined with functional studies identified sox4 + stem cells as the likely progenitor population for multiple neuronal subtypes. Further, we found that wound-induced sox4 expression is likely under direct transcriptional control by NFY, uncovering a mechanism for how early wound-induced binding of a transcriptional regulator results in the initiation of a neuronal differentiation pathway.

Highlights

A new chromosome-scale assembly for Hofstenia enables comprehensive analysis of transcription factor binding during regeneration NFY motifs become dynamically bound by 1hpa in regenerating tail fragments, particularly in the loci of neural genes A sox4 + neural-specialized stem cell is identified using scRNA-seq sox4 is wound-induced and required for differentiation of multiple neural cell types NFY regulates wound-induced expression of sox4 during regeneration.

Low-nutrient diet in Drosophila larvae stage causes enhancement in dopamine modulation in adult brain due epigenetic imprinting

Nutrient scarcity is a frequent adverse condition that organisms face during their development. This condition may lead to long-lasting effects on the metabolism and behaviour of adults due to developmental epigenetic modifications. Here, we show that reducing nutrient availability during larval development affects adult spontaneous activity and sleep behaviour, together with changes in gene expression and epigenetic marks in the mushroom bodies (MBs). We found that open chromatin regions map to 100 of 241 transcriptionally upregulated genes in the adult MBs, these new opening zones are preferentially located in regulatory zones such as promoter-TSS and introns. Importantly, opened chromatin at the Dopamine 1-like receptor 2 regulatory zones correlate with increased expression. In consequence, adult administration of a dopamine antagonist reverses increased spontaneous activity and diminished sleep time observed in response to early-life nutrient restriction. In comparison, reducing dop1R2 expression in MBs also ameliorates these effects, albeit to a lesser degree. These results lead to the conclusion that increased dopamine signalling in the MBs of flies reared in a poor nutritional environment underlies the behavioural changes observed due to this condition during development.

MORC proteins regulate transcription factor binding by mediating chromatin compaction in active chromatin regions

BACKGROUND

The microrchidia (MORC) proteins are a family of evolutionarily conserved GHKL-type ATPases involved in chromatin compaction and gene silencing. Arabidopsis MORC proteins act in the RNA-directed DNA methylation (RdDM) pathway, where they act as molecular tethers to ensure the efficient establishment of RdDM and de novo gene silencing. However, MORC proteins also have RdDM-independent functions although their underlying mechanisms are unknown.

RESULTS

In this study, we examine MORC binding regions where RdDM does not occur in order to shed light on the RdDM-independent functions of MORC proteins. We find that MORC proteins compact chromatin and reduce DNA accessibility to transcription factors, thereby repressing gene expression. We also find that MORC-mediated repression of gene expression is particularly important under conditions of stress. MORC-regulated transcription factors can in some cases regulate their own transcription, resulting in feedback loops.

CONCLUSIONS

Our findings provide insights into the molecular mechanisms of MORC-mediated chromatin compaction and transcription regulation.

High-throughput methods for the analysis of transcription factors and chromatin modifications: Low input, single cell and spatial genomic technologies

Genome-wide analysis of transcription factors and epigenomic features is instrumental to shed light on DNA-templated regulatory processes such as transcription, cellular differentiation or to monitor cellular responses to environmental cues. Two decades of technological developments have led to a rich set of approaches progressively pushing the limits of epigenetic profiling towards single cells. More recently, disruptive technologies using innovative biochemistry came into play. Assays such as CUT&RUN, CUT&Tag and variations thereof show considerable potential to survey multiple TFs or histone modifications in parallel from a single experiment and in native conditions. These are in the path to become the dominant assays for genome-wide analysis of TFs and chromatin modifications in bulk, single-cell, and spatial genomic applications. The principles together with pros and cons are discussed.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

The contributions of DNA accessibility and transcription factor occupancy to enhancer activity during cellular differentiation

The upregulation of gene expression by enhancers depends upon the interplay between the binding of sequence-specific transcription factors (TFs) and DNA accessibility. DNA accessibility is thought to limit the ability of TFs to bind to their sites, while TFs can increase accessibility to recruit additional factors that upregulate gene expression. Given this interplay, the causative regulatory events underlying the modulation of gene expression during cellular differentiation remain unknown for the vast majority of genes. We investigated the binding-site resolution dynamics of DNA accessibility and the expression dynamics of the enhancers of an important neutrophil gene, Cebpa, during macrophage-neutrophil differentiation. Reporter genes were integrated in a site-specific manner in PUER cells, which are progenitors that can be differentiated into neutrophils or macrophages in vitro by activating the pan-leukocyte TF PU.1. Time series data show that two enhancers upregulate reporter expression during the first 48 hours of neutrophil differentiation. Surprisingly, there is little or no increase in the total accessibility, measured by ATAC-Seq, of the enhancers during the same time period. Conversely, total accessibility peaks 96 hrs after PU.1 activation-consistent with its role as a pioneer-but the enhancers do not upregulate gene expression. Combining deeply sequenced ATAC-Seq data with a new bias-correction method allowed the profiling of accessibility at single-nucleotide resolution and revealed protected regions in the enhancers that match all previously characterized TF binding sites and ChIP-Seq data. Although the accessibility of most positions does not change during early differentiation, that of positions neighboring TF binding sites, an indicator of TF occupancy, did increase significantly. The localized accessibility changes are limited to nucleotides neighboring C/EBP-family TF binding sites, showing that the upregulation of enhancer activity during early differentiation is driven by C/EBP-family TF binding. These results show that increasing the total accessibility of enhancers is not sufficient for upregulating their activity and other events such as TF binding are necessary for upregulation. Also, TF binding can cause upregulation without a perceptible increase in total accessibility. Finally, this study demonstrates the feasibility of comprehensively mapping individual TF binding sites as footprints using high coverage ATAC-Seq and inferring the sequence of events in gene regulation by combining with time-series gene expression data.

ETV4 mediates dosage-dependent prostate tumor initiation and cooperates with p53 loss to generate prostate cancer

The mechanisms underlying ETS-driven prostate cancer initiation and progression remain poorly understood due to a lack of model systems that recapitulate this phenotype. We generated a genetically engineered mouse with prostate-specific expression of the ETS factor, ETV4, at lower and higher protein dosage through mutation of its degron. Lower-level expression of ETV4 caused mild luminal cell expansion without histologic abnormalities, and higher-level expression of stabilized ETV4 caused prostatic intraepithelial neoplasia (mPIN) with 100% penetrance within 1 week. Tumor progression was limited by p53-mediated senescence and Trp53 deletion cooperated with stabilized ETV4. The neoplastic cells expressed differentiation markers such as Nkx3.1 recapitulating luminal gene expression features of untreated human prostate cancer. Single-cell and bulk RNA sequencing showed that stabilized ETV4 induced a previously unidentified luminal-derived expression cluster with signatures of cell cycle, senescence, and epithelial-to-mesenchymal transition. These data suggest that ETS overexpression alone, at sufficient dosage, can initiate prostate neoplasia.

Three-dimensional chromatin re-organization during muscle stem cell aging

Age-related skeletal muscle atrophy or sarcopenia is a significant societal problem that is becoming amplified as the world's population continues to increase. The regeneration of damaged skeletal muscle is mediated by muscle stem cells, but in old age muscle stem cells become functionally attenuated. The molecular mechanisms that govern muscle stem cell aging encompass changes across multiple regulatory layers and are integrated by the three-dimensional organization of the genome. To quantitatively understand how hierarchical chromatin architecture changes during muscle stem cell aging, we generated 3D chromatin conformation maps (Hi-C) and integrated these datasets with multi-omic (chromatin accessibility and transcriptome) profiles from bulk populations and single cells. We observed that muscle stem cells display static behavior at global scales of chromatin organization during aging and extensive rewiring of local contacts at finer scales that were associated with variations in transcription factor binding and aberrant gene expression. These data provide insights into genome topology as a regulator of molecular function in stem cell aging.

Somatic chromosome pairing has a determinant impact on 3D chromatin organization

In the nucleus, chromatin is intricately structured into multiple layers of 3D organization important for genome activity. How distinct layers influence each other is not well understood. In particular, the contribution of chromosome pairing to 3D chromatin organization has been largely neglected. Here, we address this question in Drosophila, an organism that shows robust chromosome pairing in interphasic somatic cells. The extent of chromosome pairing depends on the balance between pairing and anti-pairing factors, with the anti-pairing activity of the CAP-H2 condensin II subunit being the best documented. Here, we identify the zinc-finger protein Z4 as a strong anti-pairer that interacts with and mediates the chromatin binding of CAP-H2. We also report that hyperosmotic cellular stress induces fast and reversible chromosome unpairing that depends on Z4/CAP-H2. And, most important, by combining Z4 depletion and osmostress, we show that chromosome pairing reinforces intrachromosomal 3D interactions. On the one hand, pairing facilitates RNAPII occupancy that correlates with enhanced intragenic gene-loop interactions. In addition, acting at a distance, pairing reinforces chromatin-loop interactions mediated by Polycomb (Pc). In contrast, chromosome pairing does not affect which genomic intervals segregate to active (A) and inactive (B) compartments, with only minimal effects on the strength of A-A compartmental interactions. Altogether, our results unveil the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions, unraveling the interwoven relationship between different layers of chromatin organization and the essential contribution of chromosome pairing.

Transcriptional Activation of Regenerative Hematopoiesis via Vascular Niche Sensing

Transition between activation and quiescence programs in hematopoietic stem and progenitor cells (HSC/HSPCs) is perceived to be governed intrinsically and by microenvironmental co-adaptation. However, HSC programs dictating both transition and adaptability, remain poorly defined. Single cell multiome analysis divulging differential transcriptional activity between distinct HSPC states, indicated for the exclusive absence of Fli-1 motif from quiescent HSCs. We reveal that Fli-1 activity is essential for HSCs during regenerative hematopoiesis. Fli-1 directs activation programs while manipulating cellular sensory and output machineries, enabling HSPCs co-adoptability with a stimulated vascular niche. During regenerative conditions, Fli-1 presets and enables propagation of niche-derived Notch1 signaling. Constitutively induced Notch1 signaling is sufficient to recuperate functional HSC impairments in the absence of Fli-1. Applying FLI-1 modified-mRNA transduction into lethargic adult human mobilized HSPCs, enables their vigorous niche-mediated expansion along with superior engraftment capacities. Thus, decryption of stem cell activation programs offers valuable insights for immune regenerative medicine.

Single-cell multi-scale footprinting reveals the modular organization of DNA regulatory elements

Cis-regulatory elements control gene expression and are dynamic in their structure, reflecting changes to the composition of diverse effector proteins over time1-3. Here we sought to connect the structural changes at cis-regulatory elements to alterations in cellular fate and function. To do this we developed PRINT, a computational method that uses deep learning to correct sequence bias in chromatin accessibility data and identifies multi-scale footprints of DNA-protein interactions. We find that multi-scale footprints enable more accurate inference of TF and nucleosome binding. Using PRINT with single-cell multi-omics, we discover wide-spread changes to the structure and function of candidate cis-regulatory elements (cCREs) across hematopoiesis, wherein nucleosomes slide, expose DNA for TF binding, and promote gene expression. Activity segmentation using the co-variance across cell states identifies "sub-cCREs" as modular cCRE subunits of regulatory DNA. We apply this single-cell and PRINT approach to characterize the age-associated alterations to cCREs within hematopoietic stem cells (HSCs). Remarkably, we find a spectrum of aging alterations among HSCs corresponding to a global gain of sub-cCRE activity while preserving cCRE accessibility. Collectively, we reveal the functional importance of cCRE structure across cell states, highlighting changes to gene regulation at single-cell and single-base-pair resolution.

Tcf21 marks visceral adipose mesenchymal progenitors and functions as a rate-limiting factor during visceral adipose tissue development

Distinct locations of different white adipose depots suggest anatomy-specific developmental regulation, a relatively understudied concept. Here, we report a population of Tcf21 lineage cells (Tcf21 LCs) present exclusively in visceral adipose tissue (VAT) that dynamically contributes to VAT development and expansion. During development, the Tcf21 lineage gives rise to adipocytes. In adult mice, Tcf21 LCs transform into a fibrotic or quiescent state. Multiomics analyses show consistent gene expression and chromatin accessibility changes in Tcf21 LC, based on which we constructed a gene-regulatory network governing Tcf21 LC activities. Furthermore, single-cell RNA sequencing (scRNA-seq) identifies the heterogeneity of Tcf21 LCs. Loss of Tcf21 promotes the adipogenesis and developmental progress of Tcf21 LCs, leading to improved metabolic health in the context of diet-induced obesity. Mechanistic studies show that the inhibitory effect of Tcf21 on adipogenesis is at least partially mediated via Dlk1 expression accentuation.

Integrative epigenetic analysis reveals AP-1 promotes activation of tumor-infiltrating regulatory T cells in HCC

Regulatory T (Treg) cells that infiltrate human tumors exhibit stronger immunosuppressive activity compared to peripheral blood Treg cells (PBTRs), thus hindering the induction of effective antitumor immunity. Previous transcriptome studies have identified a set of genes that are conserved in tumor-infiltrating Treg cells (TITRs). However, epigenetic profiles of TITRs have not yet been completely deciphered. Here, we employed ATAC-seq and CUT&Tag assays to integrate transcriptome profiles and identify functional regulatory elements in TITRs. We observed a global difference in chromatin accessibility and enhancer landscapes between TITRs and PBTRs. We identified two types of active enhancer formation in TITRs. The H3K4me1-predetermined enhancers are poised to be activated in response to tumor microenvironmental stimuli. We found that AP-1 family motifs are enriched at the enhancer regions of TITRs. Finally, we validated that c-Jun binds at regulatory regions to regulate signature genes of TITRs and AP-1 is required for Treg cells activation in vitro. High c-Jun expression is correlated with poor survival in human HCC. Overall, our results provide insights into the mechanism of AP-1-mediated activation of TITRs and can hopefully be used to develop new therapeutic strategies targeting TITRs in liver cancer treatment.

KDM6A Loss Triggers an Epigenetic Switch That Disrupts Urothelial Differentiation and Drives Cell Proliferation in Bladder Cancer

Disruption of KDM6A, a histone lysine demethylase, is one of the most common somatic alternations in bladder cancer. Insights into how KDM6A mutations affect the epigenetic landscape to promote carcinogenesis could help reveal potential new treatment approaches. Here, we demonstrated that KDM6A loss triggers an epigenetic switch that disrupts urothelial differentiation and induces a neoplastic state characterized by increased cell proliferation. In bladder cancer cells with intact KDM6A, FOXA1 interacted with KDM6A to activate genes instructing urothelial differentiation. KDM6A-deficient cells displayed simultaneous loss of FOXA1 target binding and genome-wide redistribution of the bZIP transcription factor ATF3, which in turn repressed FOXA1-target genes and activated cell-cycle progression genes. Importantly, ATF3 depletion reversed the cell proliferation phenotype induced by KDM6A deficiency. These data establish that KDM6A loss engenders an epigenetic state that drives tumor growth in an ATF3-dependent manner, creating a potentially targetable molecular vulnerability.

SIGNIFICANCE

A gain-of-function epigenetic switch that disrupts differentiation is triggered by inactivating KDM6A mutations in bladder cancer and can serve as a potential target for novel therapies.

Spontaneously evolved progenitor niches escape Yap oncogene addiction in advanced pancreatic ductal adenocarcinomas

Lineage plasticity has been proposed as a major source of intratumoral heterogeneity and therapeutic resistance. Here, by employing an inducible genetic engineered mouse model, we illustrate that lineage plasticity enables advanced Pancreatic Ductal Adenocarcinoma (PDAC) tumors to develop spontaneous relapse following elimination of the central oncogenic driver - Yap. Transcriptomic and immunohistochemistry analysis of a large panel of PDAC tumors reveals that within high-grade tumors, small niches of PDAC cells gradually evolve to re-activate pluripotent transcription factors (PTFs), which lessen their dependency on Yap. Comprehensive Cut&Tag analysis demonstrate that although acquisition of PTF expression is coupled with the process of epithelial-to-mesenchymal transition (EMT), PTFs form a core transcriptional regulatory circuitry (CRC) with Jun to overcome Yap dependency, which is distinct from the classic TGFb-induced EMT-TF network. A chemical-genetic screen and follow-up functional studies establish Brd4 as an epigenetic gatekeeper for the PTF-Jun CRC, and strong synergy between BET and Yap inhibitors in blocking PDAC growth.

Manipulation of the nucleoscaffold potentiates cellular reprogramming kinetics

Somatic cell fate is an outcome set by the activities of specific transcription factors and the chromatin landscape and is maintained by gene silencing of alternate cell fates through physical interactions with the nuclear scaffold. Here, we evaluate the role of the nuclear scaffold as a guardian of cell fate in human fibroblasts by comparing the effects of transient loss (knockdown) and mutation (progeria) of functional Lamin A/C, a core component of the nuclear scaffold. We observed that Lamin A/C deficiency or mutation disrupts nuclear morphology, heterochromatin levels, and increases access to DNA in lamina-associated domains. Changes in Lamin A/C were also found to impact the mechanical properties of the nucleus when measured by a microfluidic cellular squeezing device. We also show that transient loss of Lamin A/C accelerates the kinetics of cellular reprogramming to pluripotency through opening of previously silenced heterochromatin domains while genetic mutation of Lamin A/C into progerin induces a senescent phenotype that inhibits the induction of reprogramming genes. Our results highlight the physical role of the nuclear scaffold in safeguarding cellular fate.

RGT: a toolbox for the integrative analysis of high throughput regulatory genomics data

BACKGROUND

Massive amounts of data are produced by combining next-generation sequencing with complex biochemistry techniques to characterize regulatory genomics profiles, such as protein-DNA interaction and chromatin accessibility. Interpretation of such high-throughput data typically requires different computation methods. However, existing tools are usually developed for a specific task, which makes it challenging to analyze the data in an integrative manner.

RESULTS

We here describe the Regulatory Genomics Toolbox (RGT), a computational library for the integrative analysis of regulatory genomics data. RGT provides different functionalities to handle genomic signals and regions. Based on that, we developed several tools to perform distinct downstream analyses, including the prediction of transcription factor binding sites using ATAC-seq data, identification of differential peaks from ChIP-seq data, and detection of triple helix mediated RNA and DNA interactions, visualization, and finding an association between distinct regulatory factors.

CONCLUSION

We present here RGT; a framework to facilitate the customization of computational methods to analyze genomic data for specific regulatory genomics problems. RGT is a comprehensive and flexible Python package for analyzing high throughput regulatory genomics data and is available at: https://github.com/CostaLab/reg-gen . The documentation is available at: https://reg-gen.readthedocs.io.

The GRN concept as a guide for evolutionary developmental biology

Organismal phenotypes result largely from inherited developmental programs, usually executed during embryonic and juvenile life stages. These programs are not blank slates onto which natural selection can draw arbitrary forms. Rather, the mechanisms of development play an integral role in shaping phenotypic diversity and help determine the evolutionary trajectories of species. Modern evolutionary biology must, therefore, account for these mechanisms in both theory and in practice. The gene regulatory network (GRN) concept represents a potent tool for achieving this goal whose utility has grown in tandem with advances in "omic" technologies and experimental techniques. However, while the GRN concept is widely utilized, it is often less clear what practical implications it has for conducting research in evolutionary developmental biology. In this Perspective, we attempt to provide clarity by discussing how experiments and projects can be designed in light of the GRN concept. We first map familiar biological notions onto the more abstract components of GRN models. We then review how diverse functional genomic approaches can be directed toward the goal of constructing such models and discuss current methods for functionally testing evolutionary hypotheses that arise from them. Finally, we show how the major steps of GRN model construction and experimental validation suggest generalizable workflows that can serve as a scaffold for project design. Taken together, the practical implications that we draw from the GRN concept provide a set of guideposts for studies aiming at unraveling the molecular basis of phenotypic diversity.

Annelid functional genomics reveal the origins of bilaterian life cycles

Indirect development with an intermediate larva exists in all major animal lineages1, which makes larvae central to most scenarios of animal evolution2-11. Yet how larvae evolved remains disputed. Here we show that temporal shifts (that is, heterochronies) in trunk formation underpin the diversification of larvae and bilaterian life cycles. We performed chromosome-scale genome sequencing in the annelid Owenia fusiformis with transcriptomic and epigenomic profiling during the life cycles of this and two other annelids. We found that trunk development is deferred to pre-metamorphic stages in the feeding larva of O. fusiformis but starts after gastrulation in the non-feeding larva with gradual metamorphosis of Capitella teleta and the direct developing embryo of Dimorphilus gyrociliatus. Accordingly, the embryos of O. fusiformis develop first into an enlarged anterior domain that forms larval tissues and the adult head12. Notably, this also occurs in the so-called 'head larvae' of other bilaterians13-17, with which the O. fusiformis larva shows extensive transcriptomic similarities. Together, our findings suggest that the temporal decoupling of head and trunk formation, as maximally observed in head larvae, facilitated larval evolution in Bilateria. This diverges from prevailing scenarios that propose either co-option9,10 or innovation11 of gene regulatory programmes to explain larva and adult origins.

General transcription factor TAF4 antagonizes epigenetic silencing by Polycomb to maintain intestine stem cell functions

Taf4 (TATA-box binding protein-associated factor 4) is a subunit of the general transcription factor TFIID, a component of the RNA polymerase II pre-initiation complex that interacts with tissue-specific transcription factors to regulate gene expression. Properly regulated gene expression is particularly important in the intestinal epithelium that is constantly renewed from stem cells. Tissue-specific inactivation of Taf4 in murine intestinal epithelium during embryogenesis compromised gut morphogenesis and the emergence of adult-type stem cells. In adults, Taf4 loss impacted the stem cell compartment and associated Paneth cells in the stem cell niche, epithelial turnover and differentiation of mature cells, thus exacerbating the response to inflammatory challenge. Taf4 inactivation ex vivo in enteroids prevented budding formation and maintenance and caused broad chromatin remodeling and a strong reduction in the numbers of stem and progenitor cells with a concomitant increase in an undifferentiated cell population that displayed high activity of the Ezh2 and Suz12 components of Polycomb Repressive Complex 2 (PRC2). Treatment of Taf4-mutant enteroids with a specific Ezh2 inhibitor restored buddings, cell proliferation and the stem/progenitor compartment. Taf4 loss also led to increased PRC2 activity in cells of adult crypts associated with modification of the immune/inflammatory microenvironment that potentiated Apc-driven tumorigenesis. Our results reveal a novel function of Taf4 in antagonizing PRC2-mediated repression of the stem cell gene expression program to assure normal development, homeostasis, and immune-microenvironment of the intestinal epithelium.

Genomic profiling of HIV-1 integration in microglia cells links viral integration to the topologically associated domains

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary T cell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+ T cells, highlighting the importance of host genome organization in HIV-1 infection.