Developmental and evolutionary dynamics of cis-regulatory elements in mouse cerebellar cells

Single-nucleus multi-omics identifies shared and distinct pathways in Pick's and Alzheimer's disease

The study of neurodegenerative diseases, particularly tauopathies like Pick's disease (PiD) and Alzheimer's disease (AD), offers insights into the underlying regulatory mechanisms. By investigating epigenomic variations in these conditions, we identified critical regulatory changes driving disease progression, revealing potential therapeutic targets. Our comparative analyses uncovered disease-enriched non-coding regions and genome-wide transcription factor (TF) binding differences, linking them to target genes. Notably, we identified a distal human-gained enhancer (HGE) associated with E3 ubiquitin ligase (UBE3A), highlighting disease-specific regulatory alterations. Additionally, fine-mapping of AD risk genes uncovered loci enriched in microglial enhancers and accessible in other cell types. Shared and distinct TF binding patterns were observed in neurons and glial cells across PiD and AD. We validated our findings using CRISPR to excise a predicted enhancer region in UBE3A and developed an interactive database (http://swaruplab.bio.uci.edu/scROAD) to visualize predicted single-cell TF occupancy and regulatory networks.

Flexible use of conserved motif vocabularies constrains genome access in cell type evolution

Cell types evolve into a hierarchy with related types grouped into families. How cell type diversification is constrained by the stable separation between families over vast evolutionary times remains unknown. Here, integrating single-nucleus multiomic sequencing and deep learning, we show that hundreds of sequence features (motifs) divide into distinct sets associated with accessible genomes of specific cell type families. This division is conserved across highly divergent, early-branching animals including flatworms and cnidarians. While specific interactions between motifs delineate cell type relationships within families, surprisingly, these interactions are not conserved between species. Consistently, while deep learning models trained on one species can predict accessibility of other species' sequences, their predictions frequently rely on distinct, but synonymous, motif combinations. We propose that long-term stability of cell type families is maintained through genome access specified by conserved motif sets, or 'vocabularies', whereas cell types diversify through flexible use of motifs within each set.

Dissecting the regulatory logic of specification and differentiation during vertebrate embryogenesis

The interplay between transcription factors and chromatin accessibility regulates cell type diversification during vertebrate embryogenesis. To systematically decipher the gene regulatory logic guiding this process, we generated a single-cell multi-omics atlas of RNA expression and chromatin accessibility during early zebrafish embryogenesis. We developed a deep learning model to predict chromatin accessibility based on DNA sequence and found that a small number of transcription factors underlie cell-type-specific chromatin landscapes. While Nanog is well-established in promoting pluripotency, we discovered a new function in priming the enhancer accessibility of mesendodermal genes. In addition to the classical stepwise mode of differentiation, we describe instant differentiation, where pluripotent cells skip intermediate fate transitions and terminally differentiate. Reconstruction of gene regulatory interactions reveals that this process is driven by a shallow network in which maternally deposited regulators activate a small set of transcription factors that co-regulate hundreds of differentiation genes. Notably, misexpression of these transcription factors in pluripotent cells is sufficient to ectopically activate their targets. This study provides a rich resource for analyzing embryonic gene regulation and reveals the regulatory logic of instant differentiation.

A single-cell atlas of chromatin accessibility in mouse organogenesis

Organogenesis is a highly complex and precisely regulated process. Here we profiled the chromatin accessibility in >350,000 cells derived from 13 mouse embryos at four developmental stages from embryonic day (E) 10.5 to E13.5 by SPATAC-seq in a single experiment. The resulting atlas revealed the status of 830,873 candidate cis-regulatory elements in 43 major cell types. By integrating the chromatin accessibility atlas with the previous transcriptomic dataset, we characterized cis-regulatory sequences and transcription factors associated with cell fate commitment, such as Nr5a2 in the development of gastrointestinal tract, which was preliminarily supported by the in vivo experiment in zebrafish. Finally, we integrated this atlas with the previous single-cell chromatin accessibility dataset from 13 adult mouse tissues to delineate the developmental stage-specific gene regulatory programmes within and across different cell types and identify potential molecular switches throughout lineage development. This comprehensive dataset provides a foundation for exploring transcriptional regulation in organogenesis.

Single-nucleus transcriptional and chromatin accessible profiles reveal critical cell types and molecular architecture underlying chicken sex determination

INTRODUCTION

Understanding the sex determination mechanisms in birds has great significance for the biological sciences and production in the poultry industry. Sex determination in chickens is a complex process that involves fate decisions of supporting cells such as granulosa or Sertoli cells. However, a systematic understanding of the genetic regulation and cell commitment process underlying sex determination in chickens is still lacking.

OBJECTIVES

We aimed to dissect the molecular characteristics associated with sex determination in the gonads of chicken embryos.

METHODS

Single-nucleus RNA-seq (snRNA-seq) and ATAC-seq (snATAC-seq) analysis were conducted on the gonads of female and male chickens at embryonic day 3.5 (E3.5), E4.5, and E5.5.

RESULTS

Here, we provided a time-course transcriptional and chromatin accessible profiling of gonads during chicken sex determination at single-cell resolution. We uncovered differences in cell composition and developmental trajectories between female and male gonads and found that the divergence of transcription and accessibility in gonadal cells first emerged at E5.5. Furthermore, we revealed key cell-type-specific transcription factors (TFs) and regulatory networks that drive lineage commitment. Sex determination signaling pathways, dominated by BMP signaling, are preferentially activated in males during gonadal development. Further pseudotime analysis of the supporting cells indicated that granulosa cells were regulated mainly by the TEAD gene family and that Sertoli cells were driven by the DMRT1 regulons. Cross-species analysis suggested high conservation of both cell types and cell-lineage-specific TFs across the six vertebrates.

CONCLUSIONS

Overall, our study will contribute to accelerating the development of sex manipulation technology in the poultry industry and the application of chickens as a unique model for studying cell fate decisions.

The Hypothalamic Epigenetic Landscape in Dietary Obesity

The hypothalamus in the brain plays a pivotal role in controlling energy balance in vertebrates. Nutritional excess through high-fat diet (HFD) feeding can dysregulate hypothalamic signaling at multiple levels. Yet, it remains largely unknown in what magnitude HFD feeding may impact epigenetics in this brain region. Here, it is shown that HFD feeding can significantly alter hypothalamic epigenetic events, including posttranslational histone modifications, DNA methylation, and chromatin accessibility. The authors comprehensively analyze the chromatin immunoprecipitation-sequencing (ChIP-seq), methylated DNA immunoprecipitation-sequencing (MeDIP-seq), single nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), and RNA-seq data of the hypothalamus of C57 BL/6 mice fed with a chow or HFD for 1 to 6 months. The chromatins are categorized into 6 states using the obtained ChIP-seq data for H3K4me3, H3K27ac, H3K9me3, H3K27me3, and H3K36me3. A 1-month HFD feeding dysregulates histone modifications and DNA methylation more pronouncedly than that of 3- or 6-month. Besides, HFD feeding differentially impacts chromatin accessibility in hypothalamic cells. Thus, the epigenetic landscape is dysregulated in the hypothalamus of dietary obesity mice.

Single-cell multi-omics analysis of lineage development and spatial organization in the human fetal cerebellum

Human cerebellum encompasses numerous neurons, exhibiting a distinct developmental paradigm from cerebrum. Here we conducted scRNA-seq, scATAC-seq and spatial transcriptomic analyses of fetal samples from gestational week (GW) 13 to 18 to explore the emergence of cellular diversity and developmental programs in the developing human cerebellum. We identified transitory granule cell progenitors that are conserved across species. Special patterns in both granule cells and Purkinje cells were dissected multidimensionally. Species-specific gene expression patterns of cerebellar lobes were characterized and we found that PARM1 exhibited inconsistent distribution in human and mouse granule cells. A novel cluster of potential neuroepithelium at the rhombic lip was identified. We also resolved various subtypes of Purkinje cells and unipolar brush cells and revealed gene regulatory networks controlling their diversification. Therefore, our study offers a valuable multi-omics landscape of human fetal cerebellum and advances our understanding of development and spatial organization of human cerebellum.

scATAC-Ref: a reference of scATAC-seq with known cell labels in multiple species

Chromatin accessibility profiles at single cell resolution can reveal cell type-specific regulatory programs, help dissect highly specialized cell functions and trace cell origin and evolution. Accurate cell type assignment is critical for effectively gaining biological and pathological insights, but is difficult in scATAC-seq. Hence, by extensively reviewing the literature, we designed scATAC-Ref (https://bio.liclab.net/scATAC-Ref/), a manually curated scATAC-seq database aimed at providing a comprehensive, high-quality source of chromatin accessibility profiles with known cell labels across broad cell types. Currently, scATAC-Ref comprises 1 694 372 cells with known cell labels, across various biological conditions, >400 cell/tissue types and five species. We used uniform system environment and software parameters to perform comprehensive downstream analysis on these chromatin accessibility profiles with known labels, including gene activity score, TF enrichment score, differential chromatin accessibility regions, pathway/GO term enrichment analysis and co-accessibility interactions. The scATAC-Ref also provided a user-friendly interface to query, browse and visualize cell types of interest, thereby providing a valuable resource for exploring epigenetic regulation in different tissues and cell types.

Cellular development and evolution of the mammalian cerebellum

The expansion of the neocortex, a hallmark of mammalian evolution1,2, was accompanied by an increase in cerebellar neuron numbers3. However, little is known about the evolution of the cellular programmes underlying the development of the cerebellum in mammals. In this study we generated single-nucleus RNA-sequencing data for around 400,000 cells to trace the development of the cerebellum from early neurogenesis to adulthood in human, mouse and the marsupial opossum. We established a consensus classification of the cellular diversity in the developing mammalian cerebellum and validated it by spatial mapping in the fetal human cerebellum. Our cross-species analyses revealed largely conserved developmental dynamics of cell-type generation, except for Purkinje cells, for which we observed an expansion of early-born subtypes in the human lineage. Global transcriptome profiles, conserved cell-state markers and gene-expression trajectories across neuronal differentiation show that cerebellar cell-type-defining programmes have been overall preserved for at least 160 million years. However, we also identified many orthologous genes that gained or lost expression in cerebellar neural cell types in one of the species or evolved new expression trajectories during neuronal differentiation, indicating widespread gene repurposing at the cell-type level. In sum, our study unveils shared and lineage-specific gene-expression programmes governing the development of cerebellar cells and expands our understanding of mammalian brain evolution.

Epigenetic switch reshapes epithelial progenitor cell signatures and drives inflammatory pathogenesis in hidradenitis suppurativa

Hidradenitis suppurativa (HS) is a complex inflammatory skin disease with undefined mechanistic underpinnings. Here, we investigated HS epithelial cells and demonstrated that HS basal progenitors modulate their lineage restriction and give rise to pathogenic keratinocyte clones, resulting in epidermal hyperproliferation and dysregulated inflammation in HS. When comparing to healthy epithelial stem/progenitor cells, in HS, we identified changes in gene signatures that revolve around the mitotic cell cycle, DNA damage response and repair, as well as cell-cell adhesion and chromatin remodeling. By reconstructing cell differentiation trajectory and CellChat modeling, we identified a keratinocyte population specific to HS. This population is marked by S100A7/8/9 and KRT6 family members, triggering IL1, IL10, and complement inflammatory cascades. These signals, along with HS-specific proinflammatory cytokines and chemokines, contribute to the recruitment of certain immune cells during the disease progression. Furthermore, we revealed a previously uncharacterized role of S100A8 in regulating the local chromatin environment of target loci in HS keratinocytes. Through the integration of genomic and epigenomic datasets, we identified genome-wide chromatin rewiring alongside the switch of transcription factors (TFs), which mediated HS transcriptional profiles. Importantly, we identified numerous clinically relevant inflammatory enhancers and their coordinated TFs in HS basal CD49fhigh cells. The disruption of the S100A enhancer using the CRISPR/Cas9-mediated approach or the pharmacological inhibition of the interferon regulatory transcription factor 3 (IRF3) efficiently reduced the production of HS-associated inflammatory regulators. Our study not only uncovers the plasticity of epidermal progenitor cells in HS but also elucidates the epigenetic mechanisms underlying HS pathogenesis.

Single-cell analysis of chromatin accessibility in the adult mouse brain

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Single Cell Multimodal Analyses Reveal Epigenomic and Transcriptomic Basis for Birth Defects in Maternal Diabetes

Maternal diabetes mellitus is among the most frequent environmental contributors to congenital birth defects, including heart defects and craniofacial anomalies, yet the cell types affected and mechanisms of disruption are largely unknown. Using multi-modal single cell analyses, here we show that maternal diabetes affects the epigenomic landscape of specific subsets of cardiac and craniofacial progenitors during embryogenesis. A previously unrecognized cardiac progenitor subpopulation expressing the homeodomain-containing protein ALX3 showed prominent chromatin accessibility changes and acquired a more posterior identity. Similarly, a subpopulation of neural crest-derived cells in the second pharyngeal arch, which contributes to craniofacial structures, displayed abnormalities in the epigenetic landscape and axial patterning defects. Chromatin accessibility changes in both populations were associated with increased retinoic acid signaling, known to establish anterior-posterior identity. This work highlights how an environmental insult can have highly selective epigenomic consequences on discrete cell types leading to developmental patterning defects.

Single-cell epigenomics and spatiotemporal transcriptomics reveal human cerebellar development

Human cerebellar development is orchestrated by molecular regulatory networks to achieve cytoarchitecture and coordinate motor and cognitive functions. Here, we combined single-cell transcriptomics, spatial transcriptomics and single cell chromatin accessibility states to systematically depict an integrative spatiotemporal landscape of human fetal cerebellar development. We revealed that combinations of transcription factors and cis-regulatory elements (CREs) play roles in governing progenitor differentiation and cell fate determination along trajectories in a hierarchical manner, providing a gene expression regulatory map of cell fate and spatial information for these cells. We also illustrated that granule cells located in different regions of the cerebellar cortex showed distinct molecular signatures regulated by different signals during development. Finally, we mapped single-nucleotide polymorphisms (SNPs) of disorders related to cerebellar dysfunction and discovered that several disorder-associated genes showed spatiotemporal and cell type-specific expression patterns only in humans, indicating the cellular basis and possible mechanisms of the pathogenesis of neuropsychiatric disorders.

Lifelong restructuring of 3D genome architecture in cerebellar granule cells

The cerebellum contains most of the neurons in the human brain and exhibits distinctive modes of development and aging. In this work, by developing our single-cell three-dimensional (3D) genome assay-diploid chromosome conformation capture, or Dip-C-into population-scale (Pop-C) and virus-enriched (vDip-C) modes, we resolved the first 3D genome structures of single cerebellar cells, created life-spanning 3D genome atlases for both humans and mice, and jointly measured transcriptome and chromatin accessibility during development. We found that although the transcriptome and chromatin accessibility of cerebellar granule neurons mature in early postnatal life, 3D genome architecture gradually remodels throughout life, establishing ultra-long-range intrachromosomal contacts and specific interchromosomal contacts that are rarely seen in neurons. These results reveal unexpected evolutionarily conserved molecular processes that underlie distinctive features of neural development and aging across the mammalian life span.

Robust enhancer-gene regulation identified by single-cell transcriptomes and epigenomes

Single-cell sequencing could help to solve the fundamental challenge of linking millions of cell-type-specific enhancers with their target genes. However, this task is confounded by patterns of gene co-expression in much the same way that genetic correlation due to linkage disequilibrium confounds fine-mapping in genome-wide association studies (GWAS). We developed a non-parametric permutation-based procedure to establish stringent statistical criteria to control the risk of false-positive associations in enhancer-gene association studies (EGAS). We applied our procedure to large-scale transcriptome and epigenome data from multiple tissues and species, including the mouse and human brain, to predict enhancer-gene associations genome wide. We tested the functional validity of our predictions by comparing them with chromatin conformation data and causal enhancer perturbation experiments. Our study shows how controlling for gene co-expression enables robust enhancer-gene linkage using single-cell sequencing data.

Identification and characterization of transcribed enhancers during cerebellar development through enhancer RNA analysis

BACKGROUND

The development of the brain requires precise coordination of molecular processes across many cell-types. Underpinning these events are gene expression programs which require intricate regulation by non-coding regulatory sequences known as enhancers. In the context of the developing brain, transcribed enhancers (TEs) regulate temporally-specific expression of genes critical for cell identity and differentiation. Transcription of non-coding RNAs at active enhancer sequences, known as enhancer RNAs (eRNAs), is tightly associated with enhancer activity and has been correlated with target gene expression. TEs have been characterized in a multitude of developing tissues, however their regulatory role has yet to be described in the context of embryonic and early postnatal brain development. In this study, eRNA transcription was analyzed to identify TEs active during cerebellar development, as a proxy for the developing brain. Cap Analysis of Gene Expression followed by sequencing (CAGE-seq) was conducted at 12 stages throughout embryonic and early postnatal cerebellar development.

RESULTS

Temporal analysis of eRNA transcription identified clusters of TEs that peak in activity during either embryonic or postnatal times, highlighting their importance for temporally specific developmental events. Functional analysis of putative target genes identified molecular mechanisms under TE regulation revealing that TEs regulate genes involved in biological processes specific to neurons. We validate enhancer activity using in situ hybridization of eRNA expression from TEs predicted to regulate Nfib, a gene critical for cerebellar granule cell differentiation.

CONCLUSIONS

The results of this analysis provide a valuable dataset for the identification of cerebellar enhancers and provide insight into the molecular mechanisms critical for brain development under TE regulation. This dataset is shared with the community through an online resource ( https://goldowitzlab.shinyapps.io/trans-enh-app/ ).

Regulation of early cerebellar development

The study of cerebellar development has been at the forefront of neuroscience since the pioneering work of Wilhelm His Sr., Santiago Ramón y Cajal and many others since the 19th century. They laid the foundation to identify the circuitry of the cerebellum, already revealing its stereotypic three-layered cortex and discerning several of its neuronal components. Their work was fundamental in the acceptance of the neuron doctrine, which acknowledges the key role of individual neurons in forming the basic units of the nervous system. Increasing evidence shows that the cerebellum performs a variety of homeostatic and higher order neuronal functions beyond the mere control of motor behaviour. Over the last three decades, many studies have revealed the molecular machinery that regulates distinct aspects of cerebellar development, from the establishment of a cerebellar anlage in the posterior brain to the identification of cerebellar neuron diversity at the single cell level. In this review, we focus on summarizing our current knowledge on early cerebellar development with a particular emphasis on the molecular determinants that secure neuron specification and contribute to the diversity of cerebellar neurons.

A Mesp1-dependent developmental breakpoint in transcriptional and epigenomic specification of early cardiac precursors

Transcriptional networks governing cardiac precursor cell (CPC) specification are incompletely understood owing, in part, to limitations in distinguishing CPCs from non-cardiac mesoderm in early gastrulation. We leveraged detection of early cardiac lineage transgenes within a granular single-cell transcriptomic time course of mouse embryos to identify emerging CPCs and describe their transcriptional profiles. Mesp1, a transiently expressed mesodermal transcription factor, is canonically described as an early regulator of cardiac specification. However, we observed perdurance of CPC transgene-expressing cells in Mesp1 mutants, albeit mislocalized, prompting us to investigate the scope of the role of Mesp1 in CPC emergence and differentiation. Mesp1 mutant CPCs failed to robustly activate markers of cardiomyocyte maturity and crucial cardiac transcription factors, yet they exhibited transcriptional profiles resembling cardiac mesoderm progressing towards cardiomyocyte fates. Single-cell chromatin accessibility analysis defined a Mesp1-dependent developmental breakpoint in cardiac lineage progression at a shift from mesendoderm transcriptional networks to those necessary for cardiac patterning and morphogenesis. These results reveal Mesp1-independent aspects of early CPC specification and underscore a Mesp1-dependent regulatory landscape required for progression through cardiogenesis.

Comparative single cell epigenomic analysis of gene regulatory programs in the rodent and primate neocortex

Sequence divergence of cis- regulatory elements drives species-specific traits, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains to be elucidated. We investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset, and mouse with single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome, and chromosomal conformation profiles from a total of over 180,000 cells. For each modality, we determined species-specific, divergent, and conserved gene expression and epigenetic features at multiple levels. We find that cell type-specific gene expression evolves more rapidly than broadly expressed genes and that epigenetic status at distal candidate cis -regulatory elements (cCREs) evolves faster than promoters. Strikingly, transposable elements (TEs) contribute to nearly 80% of the human-specific cCREs in cortical cells. Through machine learning, we develop sequence-based predictors of cCREs in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Lastly, we show that epigenetic conservation combined with sequence similarity helps uncover functional cis -regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Cerebellar Granule Cells Develop Non-neuronal 3D Genome Architecture over the Lifespan

The cerebellum contains most of the neurons in the human brain, and exhibits unique modes of development, malformation, and aging. For example, granule cells-the most abundant neuron type-develop unusually late and exhibit unique nuclear morphology. Here, by developing our high-resolution single-cell 3D genome assay Dip-C into population-scale (Pop-C) and virus-enriched (vDip-C) modes, we were able to resolve the first 3D genome structures of single cerebellar cells, create life-spanning 3D genome atlases for both human and mouse, and jointly measure transcriptome and chromatin accessibility during development. We found that while the transcriptome and chromatin accessibility of human granule cells exhibit a characteristic maturation pattern within the first year of postnatal life, 3D genome architecture gradually remodels throughout life into a non-neuronal state with ultra-long-range intra-chromosomal contacts and specific inter-chromosomal contacts. This 3D genome remodeling is conserved in mice, and robust to heterozygous deletion of chromatin remodeling disease-associated genes (Chd8 or Arid1b). Together these results reveal unexpected and evolutionarily-conserved molecular processes underlying the unique development and aging of the mammalian cerebellum.

Complexity of enhancer networks predicts cell identity and disease genes revealed by single-cell multi-omics analysis

Many enhancers exist as clusters in the genome and control cell identity and disease genes; however, the underlying mechanism remains largely unknown. Here, we introduce an algorithm, eNet, to build enhancer networks by integrating single-cell chromatin accessibility and gene expression profiles. The complexity of enhancer networks is assessed by two metrics: the number of enhancers and the frequency of predicted enhancer interactions (PEIs) based on chromatin co-accessibility. We apply eNet algorithm to a human blood dataset and find cell identity and disease genes tend to be regulated by complex enhancer networks. The network hub enhancers (enhancers with frequent PEIs) are the most functionally important. Compared with super-enhancers, enhancer networks show better performance in predicting cell identity and disease genes. eNet is robust and widely applicable in various human or mouse tissues datasets. Thus, we propose a model of enhancer networks containing three modes: Simple, Multiple and Complex, which are distinguished by their complexity in regulating gene expression. Taken together, our work provides an unsupervised approach to simultaneously identify key cell identity and disease genes and explore the underlying regulatory relationships among enhancers in single cells.

Integration of ATAC-Seq and RNA-Seq identifies key genes and pathways involved in the neuroprotection of S-adenosylmethionine against perioperative neurocognitive disorder

The gene-regulatory landscape is highly dynamic in healthy and diseased brains. DNA methylation is a well-known epigenetic modification that regulates gene expression, and our previous study demonstrated that S-adenosylmethionine (SAM), a methylome modulator, was a neuroprotectant against perioperative neurocognitive disorder (PND). However, the underlying mechanism remains to be elucidated. Here, we integrated an assay for transposase-accessible chromatin by sequencing (ATAC-seq) and RNA sequencing (RNA-seq) to identify the key genes and pathways involved in the neuroprotection of SAM against PND. Our RNA-seq data demonstrated that genes involved in biological processes such as Wnt signaling, inflammatory response, transcription and long-term potentiation likely mediate the neuroprotection of SAM. Our ATAC-seq data provided comprehensive maps of chromatin accessibility changes induced by laparotomy and laparotomy + SAM treatment, and functional annotation of the regions with high variations in chromatin accessibility highlighted the role of the Wnt signaling pathway in PND pathogenesis and SAM treatment. Further motif analysis identified key transcription factors (e.g., CTCF, TFDP1, TCFL5, KLF15, ZBTB14, TFAP2E) that may participate in the neuroprotection of SAM. In conclusion, the current study provides an epigenomic perspective to understand the pathogenesis of PND and its treatment by SAM.

EKLF/Klf1 regulates erythroid transcription by its pioneering activity and selective control of RNA Pol II pause-release

EKLF/Klf1 is a zinc-finger transcription activator essential for erythroid lineage commitment and terminal differentiation. Using ChIP-seq, we investigate EKLF DNA binding and transcription activation mechanisms during mouse embryonic erythropoiesis. We utilize the Nan/+ mouse that expresses the EKLF-E339D (Nan) variant mutated in its conserved zinc-finger region and address the mechanism of hypomorphic and neomorphic changes in downstream gene expression. First, we show that Nan-EKLF limits normal EKLF binding to a subset of its sites. Second, we find that ectopic binding of Nan-EKLF occurs largely at enhancers and activates transcription through pioneering activity. Third, we find that for a subset of ectopic targets, gene activation is achieved in Nan/+ only by Nan-EKLF binding to distal enhancers, leading to RNA polymerase II pause-release. These results have general applicability to understanding how a DNA binding variant factor confers dominant disruptive effects on downstream gene expression even in the presence of its normal counterpart.

Single-cell chromatin accessibility identifies enhancer networks driving gene expression during spinal cord development in mouse

Spinal cord development is precisely orchestrated by spatiotemporal gene regulatory programs. However, the underlying epigenetic mechanisms remain largely elusive. Here, we profiled single-cell chromatin accessibility landscapes in mouse neural tubes spanning embryonic days 9.5-13.5. We identified neuronal-cell-cluster-specific cis-regulatory elements in neural progenitors and neurons. Furthermore, we applied a novel computational method, eNet, to build enhancer networks by integrating single-cell chromatin accessibility and gene expression data and identify the hub enhancers within enhancer networks. It was experimentally validated in vivo for Atoh1 that knockout of the hub enhancers, but not the non-hub enhancers, markedly decreased Atoh1 expression and reduced dp1/dI1 cells. Together, our work provides insights into the epigenetic regulation of spinal cord development and a proof-of-concept demonstration of enhancer networks as a general mechanism in transcriptional regulation.

Building integrative functional maps of gene regulation

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

Purkinje Cell Patterning-Insights from Single-Cell Sequencing

Despite their homogeneous appearance, Purkinje cells are remarkably diverse with respect to their molecular phenotypes, physiological properties, afferent and efferent connectivity, as well as their vulnerability to insults. Heterogeneity in Purkinje cells arises early in development, with molecularly distinct embryonic cell clusters present soon after Purkinje cell specification. Traditional methods have characterized cerebellar development and cell types, including Purkinje cell subtypes, based on knowledge of selected markers. However, recent single-cell RNA sequencing studies provide vastly increased resolution of the whole cerebellar transcriptome. Here we draw together the results of multiple single-cell transcriptomic studies in developing and adult cerebellum in both mouse and human. We describe how this detailed transcriptomic data has increased our understanding of the intricate development and function of Purkinje cells and provides first clues into features specific to human cerebellar development.

Cerebellum lineage allocation, morphogenesis and repair: impact of interplay amongst cells

The cerebellum has a simple cytoarchitecture consisting of a folded cortex with three cell layers that surrounds a nuclear structure housing the output neurons. The excitatory neurons are generated from a unique progenitor zone, the rhombic lip, whereas the inhibitory neurons and astrocytes are generated from the ventricular zone. The growth phase of the cerebellum is driven by lineage-restricted progenitor populations derived from each zone. Research during the past decade has uncovered the importance of cell-to-cell communication between the lineages through largely unknown signaling mechanisms for regulating the scaling of cell numbers and cell plasticity during mouse development and following injury in the neonatal (P0-P14) cerebellum. This Review focuses on how the interplay between cell types is key to morphogenesis, production of robust neural circuits and replenishment of cells after injury, and ends with a discussion of the implications of the greater complexity of the human cerebellar progenitor zones for development and disease.

Molecular diversity and evolution of neuron types in the amniote brain

The existence of evolutionarily conserved regions in the vertebrate brain is well established. The rules and constraints underlying the evolution of neuron types, however, remain poorly understood. To compare neuron types across brain regions and species, we generated a cell type atlas of the brain of a bearded dragon and compared it with mouse datasets. Conserved classes of neurons could be identified from the expression of hundreds of genes, including homeodomain-type transcription factors and genes involved in connectivity. Within these classes, however, there are both conserved and divergent neuron types, precluding a simple categorization of the brain into ancestral and novel areas. In the thalamus, neuronal diversification correlates with the evolution of the cortex, suggesting that developmental origin and circuit allocation are drivers of neuronal identity and evolution.

Temporal analysis of enhancers during mouse cerebellar development reveals dynamic and novel regulatory functions

We have identified active enhancers in the mouse cerebellum at embryonic and postnatal stages which provides a view of novel enhancers active during cerebellar development. The majority of cerebellar enhancers have dynamic activity between embryonic and postnatal development. Cerebellar enhancers were enriched for neural transcription factor binding sites with temporally specific expression. Putative gene targets displayed spatially restricted expression patterns, indicating cell-type specific expression regulation. Functional analysis of target genes indicated that enhancers regulate processes spanning several developmental epochs such as specification, differentiation and maturation. We use these analyses to discover one novel regulator and one novel marker of cerebellar development: Bhlhe22 and Pax3, respectively. We identified an enrichment of de novo mutations and variants associated with autism spectrum disorder in cerebellar enhancers. Furthermore, by comparing our data with relevant brain development ENCODE histone profiles and cerebellar single-cell datasets we have been able to generalize and expand on the presented analyses, respectively. We have made the results of our analyses available online in the Developing Mouse Cerebellum Enhancer Atlas (https://goldowitzlab.shinyapps.io/developing_mouse_cerebellum_enhancer_atlas/), where our dataset can be efficiently queried, curated and exported by the scientific community to facilitate future research efforts. Our study provides a valuable resource for studying the dynamics of gene expression regulation by enhancers in the developing cerebellum and delivers a rich dataset of novel gene-enhancer associations providing a basis for future in-depth studies in the cerebellum.

Histone Deacetylase 3 Inhibitor Alleviates Cerebellar Defects in Perinatal Hypothyroid Mice by Stimulating Histone Acetylation and Transcription at Thyroid Hormone-Responsive Gene Loci

Perinatal hypothyroidism impairs cerebellar organogenesis and results in motor coordination defects. The thyroid hormone receptor binds to corepressor complexes containing histone deacetylase (HDAC) 3 in the absence of ligands and acts as a transcriptional repressor. Although histone acetylation status is strongly correlated with transcriptional regulation, its role in cerebellar development remains largely unknown. We aimed to study whether the cerebellar developmental defects induced by perinatal hypothyroidism can be rescued by treatment with a specific HDAC3 inhibitor, RGFP966. Motor coordination was analyzed using three behavioral tests. The cerebella were subjected to RT-qPCR and chromatin immunoprecipitation assays for acetylated histone H3. The treatment with RGFP966 partially reversed the cerebellar morphological defects in perinatal hypothyroid mice. These findings were associated with the alleviation of motor coordination defects in these mice. In addition, the RGFP966 administration increased the mRNA levels of cerebellar thyroid hormone-responsive genes. These increases were accompanied by augmented histone acetylation status at these gene loci. These findings indicate that HDAC3 plays an important role in the cerebellar developmental defects induced by perinatal hypothyroidism. The HDAC3 inhibitor might serve as a novel therapeutic agent for hypothyroidism-induced cerebellar defects by acetylating histone tails and stimulating transcription at thyroid hormone-responsive gene loci.

A reference single-cell regulomic and transcriptomic map of cynomolgus monkeys

Non-human primates are attractive laboratory animal models that accurately reflect both developmental and pathological features of humans. Here we present a compendium of cell types across multiple organs in cynomolgus monkeys (Macaca fascicularis) using both single-cell chromatin accessibility and RNA sequencing data. The integrated cell map enables in-depth dissection and comparison of molecular dynamics, cell-type compositions and cellular heterogeneity across multiple tissues and organs. Using single-cell transcriptomic data, we infer pseudotime cell trajectories and cell-cell communications to uncover key molecular signatures underlying their cellular processes. Furthermore, we identify various cell-specific cis-regulatory elements and construct organ-specific gene regulatory networks at the single-cell level. Finally, we perform comparative analyses of single-cell landscapes among mouse, monkey and human. We show that cynomolgus monkey has strikingly higher degree of similarities in terms of immune-associated gene expression patterns and cellular communications to human than mouse. Taken together, our study provides a valuable resource for non-human primate cell biology.

Non-human primates are attractive laboratory animal models that can accurately reflect some developmental and pathological features of humans. Here the authors chart a reference cell map of cynomolgus monkeys using both scATAC-seq and scRNA-seq data across multiple organs, providing insights into the molecular dynamics and cellular heterogeneity of this organism.

Human Cerebellar Development and Transcriptomics: Implications for Neurodevelopmental Disorders

Developmental abnormalities of the cerebellum are among the most recognized structural brain malformations in human prenatal imaging. Yet reliable information regarding their cause in humans is sparse, and few outcome studies are available to inform prognosis. We know very little about human cerebellar development, in stark contrast to the wealth of knowledge from decades of research on cerebellar developmental biology of model organisms, especially mice. Recent studies show that multiple aspects of human cerebellar development significantly differ from mice and even rhesus macaques, a nonhuman primate. These discoveries challenge many current mouse-centric models of normal human cerebellar development and models regarding the pathogenesis of several neurodevelopmental phenotypes affecting the cerebellum, including Dandy-Walker malformation and medulloblastoma. Since we cannot model what we do not know, additional normative and pathological human developmental data are essential, and new models are needed.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

Integrated single-cell transcriptomic and epigenetic study of cell state transition and lineage commitment in embryonic mouse cerebellum

Recent studies using single-cell RNA-sequencing have revealed cellular heterogeneity in the developing mammalian cerebellum, yet the regulatory logic underlying this cellular diversity remains to be elucidated. Using integrated single-cell RNA and ATAC analyses, we resolved developmental trajectories of cerebellar progenitors and identified putative trans- and cis-elements that control cell state transition. We reverse engineered gene regulatory networks (GRNs) of each cerebellar cell type. Through in silico simulations and in vivo experiments, we validated the efficacy of GRN analyses and uncovered the molecular control of a posterior transitory zone (PTZ), a distinct progenitor zone residing immediately anterior to the morphologically defined rhombic lip (RL). We showed that perturbing cell fate specification in the PTZ and RL causes posterior cerebellar vermis hypoplasia, the most common cerebellar birth defect in humans. Our study provides a foundation for comprehensive studies of developmental programs of the mammalian cerebellum.

Molecular ontology of the parabrachial nucleus

This article has been removed because of a technical problem in the rendering of the PDF. 11 February 2022.

Gene regulatory networks controlling temporal patterning, neurogenesis, and cell-fate specification in mammalian retina

Gene regulatory networks (GRNs), consisting of transcription factors and their target sites, control neurogenesis and cell-fate specification in the developing central nervous system. In this study, we use integrated single-cell RNA and single-cell ATAC sequencing (scATAC-seq) analysis in developing mouse and human retina to identify multiple interconnected, evolutionarily conserved GRNs composed of cell-type-specific transcription factors that both activate genes within their own network and inhibit genes in other networks. These GRNs control temporal patterning in primary progenitors, regulate transition from primary to neurogenic progenitors, and drive specification of each major retinal cell type. We confirm that NFI transcription factors selectively activate expression of genes promoting late-stage temporal identity in primary retinal progenitors and identify other transcription factors that regulate rod photoreceptor specification in postnatal retina. This study inventories cis- and trans-acting factors that control retinal development and can guide cell-based therapies aimed at replacing retinal neurons lost to disease.