An epigenetic barrier sets the timing of human neuronal maturation.

The pace of human brain development is highly protracted compared with most other species1-7. The maturation of cortical neurons is particularly slow, taking months to years to develop adult functions3-5. Remarkably, such protracted timing is retained in cortical neurons derived from human pluripotent stem cells (hPSCs) during in vitro differentiation or upon transplantation into the mouse brain4,8,9. Those findings suggest the presence of a cell-intrinsic clock setting the pace of neuronal maturation, although the molecular nature of this clock remains unknown. Here we identify an epigenetic developmental programme that sets the timing of human neuronal maturation. First, we developed a hPSC-based approach to synchronize the birth of cortical neurons in vitro which enabled us to define an atlas of morphological, functional and molecular maturation. We observed a slow unfolding of maturation programmes, limited by the retention of specific epigenetic factors. Loss of function of several of those factors in cortical neurons enables precocious maturation. Transient inhibition of EZH2, EHMT1 and EHMT2 or DOT1L, at progenitor stage primes newly born neurons to rapidly acquire mature properties upon differentiation. Thus our findings reveal that the rate at which human neurons mature is set well before neurogenesis through the establishment of an epigenetic barrier in progenitor cells. Mechanistically, this barrier holds transcriptional maturation programmes in a poised state that is gradually released to ensure the prolonged timeline of human cortical neuron maturation.

Scalable and unbiased sequence-informed embedding of single-cell ATAC-seq data with CellSpace.

Standard scATAC sequencing (scATAC-seq) analysis pipelines represent cells as sparse numeric vectors relative to an atlas of peaks or genomic tiles and consequently ignore genomic sequence information at accessible loci. Here we present CellSpace, an efficient and scalable sequence-informed embedding algorithm for scATAC-seq that learns a mapping of DNA k-mers and cells to the same space, to address this limitation. We show that CellSpace captures meaningful latent structure in scATAC-seq datasets, including cell subpopulations and developmental hierarchies, and can score transcription factor activities in single cells based on proximity to binding motifs embedded in the same space. Importantly, CellSpace implicitly mitigates batch effects arising from multiple samples, donors or assays, even when individual datasets are processed relative to different peak atlases. Thus, CellSpace provides a powerful tool for integrating and interpreting large-scale scATAC-seq compendia.

Multicenter integrated analysis of noncoding CRISPRi screens.

The ENCODE Consortium's efforts to annotate noncoding cis-regulatory elements (CREs) have advanced our understanding of gene regulatory landscapes. Pooled, noncoding CRISPR screens offer a systematic approach to investigate cis-regulatory mechanisms. The ENCODE4 Functional Characterization Centers conducted 108 screens in human cell lines, comprising >540,000 perturbations across 24.85 megabases of the genome. Using 332 functionally confirmed CRE-gene links in K562 cells, we established guidelines for screening endogenous noncoding elements with CRISPR interference (CRISPRi), including accurate detection of CREs that exhibit variable, often low, transcriptional effects. Benchmarking five screen analysis tools, we find that CASA produces the most conservative CRE calls and is robust to artifacts of low-specificity single guide RNAs. We uncover a subtle DNA strand bias for CRISPRi in transcribed regions with implications for screen design and analysis. Together, we provide an accessible data resource, predesigned single guide RNAs for targeting 3,275,697 ENCODE SCREEN candidate CREs with CRISPRi and screening guidelines to accelerate functional characterization of the noncoding genome.

Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis.

We present a gene-level regulatory model, single-cell ATAC + RNA linking (SCARlink), which predicts single-cell gene expression and links enhancers to target genes using multi-ome (scRNA-seq and scATAC-seq co-assay) sequencing data. The approach uses regularized Poisson regression on tile-level accessibility data to jointly model all regulatory effects at a gene locus, avoiding the limitations of pairwise gene-peak correlations and dependence on peak calling. SCARlink outperformed existing gene scoring methods for imputing gene expression from chromatin accessibility across high-coverage multi-ome datasets while giving comparable to improved performance on low-coverage datasets. Shapley value analysis on trained models identified cell-type-specific gene enhancers that are validated by promoter capture Hi-C and are 11× to 15× and 5× to 12× enriched in fine-mapped eQTLs and fine-mapped genome-wide association study (GWAS) variants, respectively. We further show that SCARlink-predicted and observed gene expression vectors provide a robust way to compute a chromatin potential vector field to enable developmental trajectory analysis.

ERG activates a stem-like proliferation-differentiation program in prostate epithelial cells with mixed basal-luminal identity.

To gain insight into how ERG translocations cause prostate cancer, we performed single cell transcriptional profiling of an autochthonous mouse model at an early stage of disease initiation. Despite broad expression of ERG in all prostate epithelial cells, proliferation was enriched in a small, stem-like population with mixed-luminal basal identity (called intermediate cells). Through a series of lineage tracing and primary prostate tissue transplantation experiments, we find that tumor initiating activity resides in a subpopulation of basal cells that co-express the luminal genes Tmprss2 and Nkx3.1 (called BasalLum) but not in the larger population of classical Krt8+ luminal cells. Upon ERG activation, BasalLum cells give rise to the highly proliferative intermediate state, which subsequently transitions to the larger population of Krt8+ luminal cells characteristic of ERG-positive human cancers. Furthermore, this proliferative population is characterized by an ERG-specific chromatin state enriched for NFkB, AP-1, STAT and NFAT binding, with implications for TF cooperativity. The fact that the proliferative potential of ERG is enriched in a small stem-like population implicates the chromatin context of these cells as a critical variable for unmasking its oncogenic activity.

CRISPR Screening Uncovers a Long-Range Enhancer for ONECUT1 in Pancreatic Differentiation and Links a Diabetes Risk Variant.

Functional enhancer annotation is a valuable first step for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants for investigation. However, unbiased enhancer discovery in physiologically relevant contexts remains a major challenge. To discover regulatory elements pertinent to diabetes, we conducted a CRISPR interference screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers uncovered, we focused on a long-range enhancer ∼664 kb from the ONECUT1 promoter, since coding mutations in ONECUT1 cause pancreatic hypoplasia and neonatal diabetes. Homozygous enhancer deletion in hPSCs was associated with a near-complete loss of ONECUT1 gene expression and compromised pancreatic differentiation. This enhancer contains a confidently fine-mapped type 2 diabetes associated variant (rs528350911) which disrupts a GATA motif. Introduction of the risk variant into hPSCs revealed substantially reduced binding of key pancreatic transcription factors (GATA4, GATA6 and FOXA2) on the edited allele, accompanied by a slight reduction of ONECUT1 transcription, supporting a causal role for this risk variant in metabolic disease. This work expands our knowledge about transcriptional regulation in pancreatic development through the characterization of a long-range enhancer and highlights the utility of enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

Thetis cells induce food-specific Treg cell differentiation and oral tolerance.

The intestinal immune system must establish tolerance to food antigens to prevent onset of allergic and inflammatory diseases. Peripherally generated regulatory T (pTreg) cells play an essential role in suppressing inflammatory responses to allergens; however, the antigen-presenting cell (APC) that instructs food-specific pTreg cells is not known. Here, we show that antigen presentation and TGF-ß activation by a subset of RORγt + antigen-presenting cells (APC), Thetis cells IV (TC IV), is required for food-induced pTreg cell differentiation and oral tolerance. By contrast, antigen presentation by dendritic cells (DCs) was dispensable for pTreg induction but required for T H 1 effector responses, highlighting a division of labor between tolerogenic TCs and pro-inflammatory DCs. While antigen presentation by TCs was required for food-specific pTreg generation both in early life and adulthood, the increased abundance of TCs in the peri-weaning period was associated with a window of opportunity for enhanced pTreg differentiation. These findings establish a critical role for TCs in oral tolerance and suggest that these cells may represent a key therapeutic target for the treatment of food-associated allergic and inflammatory diseases.

CRISPR screening uncovers a long-range enhancer for ONECUT1 in pancreatic differentiation and links a diabetes risk variant.

Functional enhancer annotation is critical for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants. However, unbiased enhancer discovery in disease-relevant contexts remains challenging. To identify enhancers pertinent to diabetes, we conducted a CRISPR interference (CRISPRi) screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers identified, we focused on an enhancer we named ONECUT1e-664kb, ∼664 kb from the ONECUT1 promoter. Previous studies have linked ONECUT1 coding mutations to pancreatic hypoplasia and neonatal diabetes. We found that homozygous deletion of ONECUT1e-664kb in hPSCs leads to a near-complete loss of ONECUT1 expression and impaired pancreatic differentiation. ONECUT1e-664kb contains a type 2 diabetes-associated variant (rs528350911) disrupting a GATA motif. Introducing the risk variant into hPSCs reduced binding of key pancreatic transcription factors (GATA4, GATA6, and FOXA2), supporting its causal role in diabetes. This work highlights the utility of unbiased enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer.

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

An epigenetically distinct HSC subset supports thymic reconstitution.

Hematopoietic stem cells (HSCs) with multilineage potential are critical for effective T cell reconstitution and restoration of the adaptive immune system after allogeneic Hematopoietic Cell Transplantation (allo-HCT). The Kit lo subset of HSCs is enriched for multipotential precursors, 1, 2 but their T-cell lineage potential has not been well-characterized. We therefore studied the thymic reconstituting and T-cell potential of Kit lo HSCs. Using a preclinical allo-HCT model, we demonstrate that Kit lo HSCs support better thymic recovery, and T-cell reconstitution resulting in improved T cell responses to infection post-HCT. Furthermore, Kit lo HSCs with augmented BM lymphopoiesis mitigate age-associated thymic alterations, thus enhancing T-cell recovery in middle-aged hosts. We find the frequency of the Kit lo subset declines with age, providing one explanation for the reduced frequency of T-competent HSCs and reduced T-lymphopoietic potential in BM precursors of aged mice. 3, 4, 5 Chromatin profiling revealed that Kit lo HSCs exhibit higher activity of lymphoid-specifying transcription factors (TFs), including Zbtb1 . Deletion of Zbtb1 in Kit lo HSCs diminished their T-cell potential, while reinstating Zbtb1 in megakaryocytic-biased Kit hi HSCs rescued T-cell potential, in vitro and in vivo . Finally, we discover an analogous Kit lo HSC subset with enhanced lymphoid potential in human bone marrow. Our results demonstrate that Kit lo HSCs with enhanced lymphoid potential have a distinct underlying epigenetic program.