Search

Cell type identity is encoded by gene regulatory networks (GRN), in which transcription factors (TFs) bind to enhancers to regulate target gene expression. In the mammalian liver, lineage TFs have been characterized for the main cell types, including hepatocytes. Hepatocytes cover a relatively broad cellular state space, as they differ significantly in their metabolic state, and function, depending on their position with respect to the central or portal vein in a liver lobule. It is unclear whether this spatially defined cellular state space, called zonation, is also governed by a well-defined gene regulatory code. To address this challenge, we have mapped enhancer-GRNs across liver cell types at high resolution, using a combination of single cell multiomics, spatial omics, GRN inference, and deep learning. We found that cell state changes in transcription and chromatin accessibility in hepatocytes, liver sinusoidal endothelial cells and hepatic stellate cells depend on zonation. Enhancer-GRN mapping suggests that zonation states in hepatocytes are driven by the repressors Tcf7l1 and Tbx3, that modulate the core hepatocyte GRN, controlled by Hnf4a, Cebpa, Hnf1a, Onecut1 and Foxa1, among others. To investigate how these TFs cooperate with cell type TFs, we performed anin vivomassively parallel reporter assay on 12,000 hepatocyte enhancers and used these data to train a hierarchical deep learning model (called DeepLiver) that exploits both enhancer accessibility and activity. DeepLiver confirms Cebpa, Onecut, Foxa1, Hnf1a and Hnf4a as drivers of enhancer specificity in hepatocytes; Tcf7l1/2 and Tbx3 as regulators of the zonation state; and Hnf4a, Hnf1a, AP-1 and Ets as activators. Finally, taking advantage ofin silicomutagenesis predictions from DeepLiver and enhancer assays, we confirmed that the destruction of Tcf7l1/2 or Tbx3 motifs in zonated enhancers abrogates their zonation bias. Our study provides a multi-modal understanding of the regulatory code underlying hepatocyte identity and their zonation state, that can be exploited to engineer enhancers with specific activity levels and zonation patterns.

Joint profiling of chromatin accessibility and gene expression of individual cells provides an opportunity to decipher enhancer-driven gene regulatory networks (eGRN). Here we present a new method for the inference of eGRNs, called SCENIC+. SCENIC+ predicts genomic enhancers along with candidate upstream transcription factors (TF) and links these enhancers to candidate target genes. Specific TFs for each cell type or cell state are predicted based on the concordance of TF binding site accessibility, TF expression, and target gene expression. To improve both recall and precision of TF identification, we curated and clustered more than 40,000 position weight matrices that we could associate with 1,553 human TFs. We validated and benchmarked each of the SCENIC+ components on diverse data sets from different species, including human peripheral blood mononuclear cell types, ENCODE cell lines, human melanoma cell states, and Drosophila retinal development. Next, we exploit SCENIC+ predictions to study conserved TFs, enhancers, and GRNs between human and mouse cell types in the cerebral cortex. Finally, we provide new capabilities that exploit the inferred eGRNs to study the dynamics of gene regulation along differentiation trajectories; to map regulatory activities onto tissues using spatial omics data; and to predict the effect of TF perturbations on cell state. SCENIC+ provides critical insight into gene regulation, starting from multiome atlases of scATAC-seq and scRNA-seq. The SCENIC+ suite is available as a set of Python modules at https://scenicplus.readthedocs.io.

Wound response programs are often activated during neoplastic growth in tumors. In both wound repair and tumor growth, cells respond to acute stress and balance the activation of multiple programs including apoptosis, proliferation, and cell migration. Central to those responses are the activation of the JNK/MAPK and JAK/STAT signaling pathways. Yet, to what extent these signaling cascades interact at the cis-regulatory level, and how they orchestrate different regulatory and phenotypic responses is still unclear. Here, we aim to characterize the regulatory states that emerge and cooperate in the wound response, using the Drosophila melanogaster wing disc as a model system, and compare these with cancer cell states induced by rasV12scrib-/- in the eye disc. We used single-cell multiome profiling to derive enhancer Gene Regulatory Networks (eGRNs) by integrating chromatin accessibility and gene expression signals. We identify a “proliferative” eGRN, active in the majority of wounded cells and controlled by AP-1 and STAT. In a smaller, but distinct population of wound cells, a “senescent” eGRN is activated and driven by C/EBP-like transcription factors (Irbp18, Xrp1, Slow border, and Vrille) and Scalloped. In tumor cells on the other hand, both eGRN signatures are simultaneously active within the same cell. Our single-cell multiome and eGRNs resource offers an in-depth characterisation of the senescence markers, together with a new perspective on the shared gene regulatory programs acting during wound response and oncogenesis.

Melanoma cells can switch between a melanocytic and a mesenchymal-like state. Scattered evidence indicates that additional intermediate state(s) may exist. Here, to search for such states and decipher their underlying gene regulatory network (GRN), we studied 10 melanoma cultures using single-cell RNA sequencing (RNA-seq) as well as 26 additional cultures using bulk RNA-seq. Although each culture exhibited a unique transcriptome, we identified shared GRNs that underlie the extreme melanocytic and mesenchymal states and the intermediate state. This intermediate state is corroborated by a distinct chromatin landscape and is governed by the transcription factors SOX6, NFATC2, EGR3, ELF1 and ETV4. Single-cell migration assays confirmed the intermediate migratory phenotype of this state. Using time-series sampling of single cells after knockdown of SOX10, we unravelled the sequential and recurrent arrangement of GRNs during phenotype switching. Taken together, these analyses indicate that an intermediate state exists and is driven by a distinct and stable 'mixed' GRN rather than being a symbiotic heterogeneous mix of cells. ispartof: NATURE CELL BIOLOGY vol:22 issue:8 pages:986-+ ispartof: location:England status: published

The Drosophila brain is a frequently used model in neuroscience. Single-cell transcriptome analysis1-6, three-dimensional morphological classification7 and electron microscopy mapping of the connectome8,9 have revealed an immense diversity of neuronal and glial cell types that underlie an array of functional and behavioural traits in the fly. The identities of these cell types are controlled by gene regulatory networks (GRNs), involving combinations of transcription factors that bind to genomic enhancers to regulate their target genes. Here, to characterize GRNs at the cell-type level in the fly brain, we profiled the chromatin accessibility of 240,919 single cells spanning 9 developmental timepoints and integrated these data with single-cell transcriptomes. We identify more than 95,000 regulatory regions that are used in different neuronal cell types, of which 70,000 are linked to developmental trajectories involving neurogenesis, reprogramming and maturation. For 40 cell types, uniquely accessible regions were associated with their expressed transcription factors and downstream target genes through a combination of motif discovery, network inference and deep learning, creating enhancer GRNs. The enhancer architectures revealed by DeepFlyBrain lead to a better understanding of neuronal regulatory diversity and can be used to design genetic driver lines for cell types at specific timepoints, facilitating their characterization and manipulation. ispartof: NATURE vol:601 issue:7894 pages:630-+ ispartof: location:England status: published

Single-cell RNA-seq and single-cell assay for transposase-accessible chromatin (ATAC-seq) technologies are used extensively to create cell type atlases for a wide range of organisms, tissues, and disease processes. To increase the scale of these atlases, lower the cost and pave the way for more specialized multiome assays, custom droplet microfluidics may provide solutions complementary to commercial setups. We developed HyDrop, a flexible and open-source droplet microfluidic platform encompassing three protocols. The first protocol involves creating dissolvable hydrogel beads with custom oligos that can be released in the droplets. In the second protocol, we demonstrate the use of these beads for HyDrop-ATAC, a low-cost noncommercial scATAC-seq protocol in droplets. After validating HyDrop-ATAC, we applied it to flash-frozen mouse cortex and generated 7996 high-quality single-cell chromatin accessibility profiles in a single run. In the third protocol, we adapt both the reaction chemistry and the capture sequence of the barcoded hydrogel bead to capture mRNA, and demonstrate a significant improvement in throughput and sensitivity compared to previous open-source droplet-based scRNA-seq assays (Drop-seq and inDrop). Similarly, we applied HyDrop-RNA to flash-frozen mouse cortex and generated 9508 single-cell transcriptomes closely matching reference single-cell gene expression data. Finally, we leveraged HyDrop-RNA's high capture rate to analyze a small population of fluorescence-activated cell sorted neurons from theScientists are now able to determine the order of chemical blocks, or nucleic acids, that make up the genetic code. These sequencing tools can be used to identify which genes are active within a biological sample. They do this by extracting and analysing open chromatin (regions of DNA that are accessible to the cell’s machinery), or sequences of RNA (the molecular templates cells use to translate genes into working proteins). Initially, most sequencing tools could only provide an ‘averaged-out’ profile of the genes activated in bulk pieces of tissue which contain multiple types of cell. However, advances in technology have led to new methods that can extract and analyse open chromatin or RNA from individual cells. First, the cells are separated, via a technique called microfluidics, into tiny droplets of water along with a single bead that carries a unique barcode. The cell is then broken apart inside the droplet and the barcode within the bead gets released and attaches itself to the genetic material extracted from the cell. All the genetic material inside the droplets is then pooled together and sequenced. Researchers then use the barcode tags to identify which bits of RNA or DNA belong to each cell. Single-cell sequencing has many advantages, including being able to pinpoint precise genetic differences between healthy and abnormal cells, and to create cell atlases of whole organisms, tissues and microbial communities. But existing methods for extracting chromatin are very expensive, and there were no openly available tools for processing thousands of cells at speed. Furthermore, while several single-cell RNA sequencing tools are already freely available, they are not very sensitive or practical to use. Here, De Rop et al. have developed a new open-source platform called HyDrop that overcomes these barriers. The method entails a new type of barcoded bead and optimised elements of existing microfluidics protocols using open-source reagents. These changes created a more user-friendly workflow and increased sensitivity of sequencing at no additional cost. De Rop et al. used their new platform to screen the RNA and open chromatin of thousands of individuals cells from the brains of mice and flies. HyDrop outperformed other open-source methods when working in RNA-sequencing mode. It also provides the first open-source tool for sequencing open chromatin in single cells. Further improvements are expected as researchers tweak the platform, which for now provides an affordable alternative to existing methods.

Single-cell RNA-seq and single-cell ATAC-seq technologies are being used extensively to create cell type atlases for a wide range of organisms, tissues, and disease processes. To increase the scale of these atlases, lower the cost, and allow for more specialized multi-ome assays, custom droplet microfluidics may provide complementary solutions to commercial setups. We developed HyDrop, a flexible and generic droplet microfluidic platform encompassing three protocols. The first protocol involves creating dissolvable hydrogel beads with custom oligos that can be released in the droplets. In the second protocol, we demonstrate the use of these beads for HyDrop-ATAC, a low-cost non-commercial scATAC-seq protocol in droplets. After validating HyDrop-ATAC, we applied it to flash-frozen mouse cortex and generated 8,502 high-quality single-cell chromatin accessibility profiles in a single run. In the third protocol, we adapt both the reaction chemistry and the capture sequence of the barcoded hydrogel bead to capture mRNA, and demonstrate a significant improvement in throughput and sensitivity compared to previous open-source droplet-based scRNA-seq assays (Drop-seq and inDrop). Similarly, we applied HyDrop-RNA to flash-frozen mouse cortex and generated 9,508 single-cell transcriptomes closely matching reference single-cell gene expression data. Finally, we leveraged HyDrop-RNA’s high capture rate to analyse a small population of FAC-sorted neurons from the Drosophila brain, confirming the protocol’s applicability to low-input samples and small cells. HyDrop is currently capable of generating single-cell data in high throughput and at a reduced cost compared to commercial methods, and we envision that HyDrop can be further developed to be compatible with novel (multi-) omics protocols.

We present cisTopic, a probabilistic framework used to simultaneously discover coaccessible enhancers and stable cell states from sparse single-cell epigenomics data ( http://github.com/aertslab/cistopic ). Using a compendium of single-cell ATAC-seq datasets from differentiating hematopoietic cells, brain and transcription factor perturbations, we demonstrate that topic modeling can be exploited for robust identification of cell types, enhancers and relevant transcription factors. cisTopic provides insight into the mechanisms underlying regulatory heterogeneity in cell populations. As an unsupervised Bayesian framework, cisTopic classifies regions in scATAC-seq data into regulatory topics, which are used for clustering.

Although single-cell RNA-seq is revolutionizing biology, data interpretation remains a challenge. We present SCENIC for the simultaneous reconstruction of gene regulatory networks and identification of cell states. We apply SCENIC to a compendium of single-cell data from tumors and brain, and demonstrate that the genomic regulatory code can be exploited to guide the identification of transcription factors and cell states. SCENIC provides critical biological insights into the mechanisms driving cellular heterogeneity.

Melanoma is notorious for its cellular heterogeneity, which is at least partly due to its ability to transition between alternate cell states. Similarly to EMT, melanoma cells with a melanocytic phenotype can switch to a mesenchymal-like phenotype. However, scattered emerging evidence indicates that additional, intermediate state(s) may exist. In order to search for such new melanoma states and decipher their underlying gene regulatory network (GRN), we extensively studied ten patient-derived melanoma cultures by single-cell RNA-seq of >39,000 cells. Although each culture exhibited a unique transcriptome, we identified shared gene regulatory networks that underlie the extreme melanocytic and mesenchymal cell states, as well as one (stable) intermediate state. The intermediate state was corroborated by a distinct open chromatin landscape and governed by the transcription factors EGR3, NFATC2, and RXRG. Single-cell migration assays established that this “transition” state exhibits an intermediate migratory phenotype. Through a dense time-series sampling of single cells and dynamic GRN inference, we unraveled the sequential and recurrent arrangement of transcriptional programs at play during phenotype switching that ultimately lead to the mesenchymal cell state. We provide the scRNA-Seq data with 39,263 melanoma cells on our SCope platform and the ATAC-seq data on a UCSC hub to jointly serve as a resource for the melanoma field. Together, this exhaustive analysis of melanoma cell state diversity indicates that additional states exists between the two extreme melanocytic and mesenchymal-like states. The GRN we identified may serve as a new putative target to prevent the switch to mesenchymal cell state and thereby, acquisition of metastatic and drug resistant potential.

Summary The diversity of cell types and regulatory states in the brain, and how these change during aging, remains largely unknown. We present a single-cell transcriptome atlas of the entire adult Drosophila melanogaster brain sampled across its lifespan. Cell clustering identified 87 initial cell clusters that are further subclustered and validated by targeted cell-sorting. Our data show high granularity and identify a wide range of cell types. Gene network analyses using SCENIC revealed regulatory heterogeneity linked to energy consumption. During aging, RNA content declines exponentially without affecting neuronal identity in old brains. This single-cell brain atlas covers nearly all cells in the normal brain and provides the tools to study cellular diversity alongside other Drosophila and mammalian single-cell datasets in our unique single-cell analysis platform: SCope (http://scope.aertslab.org). These results, together with SCope, allow comprehensive exploration of all transcriptional states of an entire aging brain., Graphical Abstract, Highlights • A single-cell atlas of the adult fly brain during aging • Network inference reveals regulatory states related to oxidative phosphorylation • Cell identity is retained during aging despite exponential decline of gene expression • SCope: An online tool to explore and compare single-cell datasets across species, A single-cell atlas of adult fly brains identifies the ensemble of neuronal and glial cell types and their dynamic changes during aging.

Single-cell epigenomics provides new opportunities to decipher genomic regulatory programs from heterogeneous samples and dynamic processes. We present a probabilistic framework called cisTopic, to simultaneously discover “cis-regulatory topics” and stable cell states from sparse single-cell epigenomics data. After benchmarking cisTopic on single-cell ATAC-seq data, single-cell DNA methylation data, and semi-simulated single-cell ChIP-seq data, we use cisTopic to predict regulatory programs in the human brain and validate these by aligning them with co-expression networks derived from single-cell RNA-seq data. Next, we performed a time-series single-cell ATAC-seq experiment after SOX10 perturbations in melanoma cultures, where cisTopic revealed dynamic regulatory topics driven by SOX10 and AP-1. Finally, machine learning and enhancer modelling approaches allowed to predict cell type specific SOX10 and SOX9 binding sites based on topic specific co-regulatory motifs. cisTopic is available as an R/Bioconductor package at http://github.com/aertslab/cistopic.

Transcriptional enhancers function as docking platforms for combinations of transcription factors (TFs) to control gene expression. How enhancer sequences determine nucleosome occupancy, TF recruitment and transcriptional activation in vivo remains unclear. Using ATAC–seq across a panel of Drosophila inbred strains, we found that SNPs affecting binding sites of the TF Grainy head (Grh) causally determine the accessibility of epithelial enhancers. We show that deletion and ectopic expression of Grh cause loss and gain of DNA accessibility, respectively. However, although Grh binding is necessary for enhancer accessibility, it is insufficient to activate enhancers. Finally, we show that human Grh homologs—GRHL1, GRHL2 and GRHL3—function similarly. We conclude that Grh binding is necessary and sufficient for the opening of epithelial enhancers but not for their activation. Our data support a model positing that complex spatiotemporal expression patterns are controlled by regulatory hierarchies in which pioneer factors, such as Grh, establish tissue-specific accessible chromatin landscapes upon which other factors can act. ispartof: NATURE GENETICS vol:50 issue:7 pages:1011-1020 ispartof: location:United States status: published

Summary Inferring a Gene Regulatory Network (GRN) from gene expression data is a computationally expensive task, exacerbated by increasing data sizes due to advances in high-throughput gene profiling technology, such as single-cell RNA-seq. To equip researchers with a toolset to infer GRNs from large expression datasets, we propose GRNBoost2 and the Arboreto framework. GRNBoost2 is an efficient algorithm for regulatory network inference using gradient boosting, based on the GENIE3 architecture. Arboreto is a computational framework that scales up GRN inference algorithms complying with this architecture. Arboreto includes both GRNBoost2 and an improved implementation of GENIE3, as a user-friendly open source Python package. Availability and implementation Arboreto is available under the 3-Clause BSD license at http://arboreto.readthedocs.io. Supplementary information Supplementary data are available at Bioinformatics online.

SummaryThe diversity of cell types and regulatory states in the brain, and how these change during ageing, remains largely unknown. Here, we present a single-cell transcriptome catalogue of the entire adult Drosophila melanogaster brain sampled across its lifespan. Both neurons and glia age through a process of “regulatory erosion”, characterized by a strong decline of RNA content, and accompanied by increasing transcriptional and chromatin noise. We identify more than 50 cell types by specific transcription factors and their downstream gene regulatory networks. In addition to neurotransmitter types and neuroblast lineages, we find a novel neuronal cell state driven by datilografo and prospero. This state relates to neuronal birth order, the metabolic profile, and the activity of a neuron. Our single-cell brain catalogue reveals extensive regulatory heterogeneity linked to ageing and brain function and will serve as a reference for future studies of genetic variation and disease mutations.

Yeast and cancer cells share the unusual characteristic of favoring fermentation of sugar over respiration. We now reveal an evolutionary conserved mechanism linking fermentation to activation of Ras, a major regulator of cell proliferation in yeast and mammalian cells, and prime proto-oncogene product. A yeast mutant (tps1∆) with overactive influx of glucose into glycolysis and hyperaccumulation of Fru1,6bisP, shows hyperactivation of Ras, which causes its glucose growth defect by triggering apoptosis. Fru1,6bisP is a potent activator of Ras in permeabilized yeast cells, likely acting through Cdc25. As in yeast, glucose triggers activation of Ras and its downstream targets MEK and ERK in mammalian cells. Biolayer interferometry measurements show that physiological concentrations of Fru1,6bisP stimulate dissociation of the pure Sos1/H-Ras complex. Thermal shift assay confirms direct binding to Sos1, the mammalian ortholog of Cdc25. Our results suggest that the Warburg effect creates a vicious cycle through Fru1,6bisP activation of Ras, by which enhanced fermentation stimulates oncogenic potency., Yeast and cancer cells both favor sugar fermentation in aerobic conditions. Here the authors describe a conserved mechanism from yeast to mammals where the glycolysis intermediate fructose-1,6-bisphosphate binds Cdc25/Sos1 and couples increased glycolytic flux to increased Ras proto-oncoprotein activity.

Single-cell RNA-seq allows building cell atlases of any given tissue and infer the dynamics of cellular state transitions during developmental or disease trajectories. Both the maintenance and transitions of cell states are encoded by regulatory programs in the genome sequence. However, this regulatory code has not yet been exploited to guide the identification of cellular states from single-cell RNA-seq data. Here we describe a computational resource, called SCENIC (Single Cell rEgulatory Network Inference and Clustering), for the simultaneous reconstruction of gene regulatory networks (GRNs) and the identification of stable cell states, using single-cell RNA-seq data. SCENIC outperforms existing approaches at the level of cell clustering and transcription factor identification. Importantly, we show that cell state identification based on GRNs is robust towards batch-effects and technical-biases. We applied SCENIC to a compendium of single-cell data from the mouse and human brain and demonstrate that the proper combinations of transcription factors, target genes, enhancers, and cell types can be identified. Moreover, we used SCENIC to map the cell state landscape in melanoma and identified a gene regulatory network underlying a proliferative melanoma state driven by MITF and STAT and a contrasting network controlling an invasive state governed by NFATC2 and NFIB. We further validated these predictions by showing that two transcription factors are predominantly expressed in early metastatic sentinel lymph nodes. In summary, SCENIC is the first method to analyze scRNA-seq data using a network-centric, rather than cell-centric approach. SCENIC is generic, easy to use, and flexible, and allows for the simultaneous tracing of genomic regulatory programs and the mapping of cellular identities emerging from these programs. Availability: SCENIC is available as an R workflow based on three new R/Bioconductor packages: GENIE3, RcisTarget and AUCell. As scalable alternative to GENIE3, we also provide GRNboost, paving the way towards the network analysis across millions of single cells.

Single-cell technologies allow measuring chromatin accessibility and gene expression in each cell, but jointly utilizing both layers to map bona fide gene regulatory networks and enhancers remains challenging. Here, we generate independent single-cell RNA-seq and single-cell ATAC-seq atlases of the Drosophila eye-antennal disc and spatially integrate the data into a virtual latent space that mimics the organization of the 2D tissue using ScoMAP (Single-Cell Omics Mapping into spatial Axes using Pseudotime ordering). To validate spatially predicted enhancers, we use a large collection of enhancer-reporter lines and identify ~ 85% of enhancers in which chromatin accessibility and enhancer activity are coupled. Next, we infer enhancer-to-gene relationships in the virtual space, finding that genes are mostly regulated by multiple, often redundant, enhancers. Exploiting cell type-specific enhancers, we deconvolute cell type-specific effects of bulk-derived chromatin accessibility QTLs. Finally, we discover that Prospero drives neuronal differentiation through the binding of a GGG motif. In summary, we provide a comprehensive spatial characterization of gene regulation in a 2D tissue. ispartof: MOLECULAR SYSTEMS BIOLOGY vol:16 issue:5 ispartof: location:England status: published