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.

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.

Single-cell transcriptional regulation and genetic evolution of neuroendocrine prostate cancer.

Neuroendocrine prostate cancer (NEPC) is a lethal subtype of prostate cancer, with a 10% five-year survival rate. However, little is known about its origin and the mechanisms governing its emergence. Our study characterized ADPC and NEPC in prostate tumors from 7 patients using scRNA-seq. First, we identified two NEPC gene expression signatures representing different phases of trans-differentiation. New marker genes we identified may be used for clinical diagnosis. Second, integrative analyses combining expression and subclonal architecture revealed different paths by which NEPC diverges from the original ADPC, either directly from treatment-naïve tumor cells or from specific intermediate states of treatment-resistance. Third, we inferred a hierarchical transcription factor (TF) network underlying the progression, which involves constitutive regulation by ASCL1, FOXA2, and selective regulation by NKX2-2, POU3F2, and SOX2. Together, these results defined the complex expression profiles and advanced our understanding of the genetic and transcriptomic mechanisms leading to NEPC differentiation.

Enhancer dependence of cell-type-specific gene expression increases with developmental age.

How overall principles of cell-type-specific gene regulation (the "logic") may change during ontogeny is largely unexplored. We compared transcriptomic, epigenomic, and three-dimensional (3D) genomic profiles in embryonic (EryP) and adult (EryD) erythroblasts. Despite reduced chromatin accessibility compared to EryP, distal chromatin of EryD is enriched in H3K27ac, Gata1, and Myb occupancy. EryP-/EryD-shared enhancers are highly correlated with red blood cell identity genes, whereas cell-type-specific regulation employs different cis elements in EryP and EryD cells. In contrast to EryP-specific genes, which exhibit promoter-centric regulation through Gata1, EryD-specific genes rely more on distal enhancers for regulation involving Myb-mediated enhancer activation. Gata1 HiChIP demonstrated an overall increased enhancer-promoter interactions at EryD-specific genes, whereas genome editing in selected loci confirmed distal enhancers are required for gene expression in EryD but not in EryP. Applying a metric for enhancer dependence of transcription, we observed a progressive reliance on cell-specific enhancers with increasing ontogenetic age among diverse tissues of mouse and human origin. Our findings highlight fundamental and conserved differences at distinct developmental stages, characterized by simpler promoter-centric regulation of cell-type-specific genes in embryonic cells and increased combinatorial enhancer-driven control in adult cells.