A single-cell atlas of chromatin accessibility in the human genome.

Current catalogs of regulatory sequences in the human genome are still incomplete and lack cell type resolution. To profile the activity of gene regulatory elements in diverse cell types and tissues in the human body, we applied single-cell chromatin accessibility assays to 30 adult human tissue types from multiple donors. We integrated these datasets with previous single-cell chromatin accessibility data from 15 fetal tissue types to reveal the status of open chromatin for ∼1.2 million candidate cis-regulatory elements (cCREs) in 222 distinct cell types comprised of >1.3 million nuclei. We used these chromatin accessibility maps to delineate cell-type-specificity of fetal and adult human cCREs and to systematically interpret the noncoding variants associated with complex human traits and diseases. This rich resource provides a foundation for the analysis of gene regulatory programs in human cell types across tissues, life stages, and organ systems.

Systematic analysis of binding of transcription factors to noncoding variants.

Many sequence variants have been linked to complex human traits and diseases1, but deciphering their biological functions remains challenging, as most of them reside in noncoding DNA. Here we have systematically assessed the binding of 270 human transcription factors to 95,886 noncoding variants in the human genome using an ultra-high-throughput multiplex protein-DNA binding assay, termed single-nucleotide polymorphism evaluation by systematic evolution of ligands by exponential enrichment (SNP-SELEX). The resulting 828 million measurements of transcription factor-DNA interactions enable estimation of the relative affinity of these transcription factors to each variant in vitro and evaluation of the current methods to predict the effects of noncoding variants on transcription factor binding. We show that the position weight matrices of most transcription factors lack sufficient predictive power, whereas the support vector machine combined with the gapped k-mer representation show much improved performance, when assessed on results from independent SNP-SELEX experiments involving a new set of 61,020 sequence variants. We report highly predictive models for 94 human transcription factors and demonstrate their utility in genome-wide association studies and understanding of the molecular pathways involved in diverse human traits and diseases.

Mediator complex proximal Tail subunit MED30 is critical for Mediator core stability and cardiomyocyte transcriptional network.

Dysregulation of cardiac transcription programs has been identified in patients and families with heart failure, as well as those with morphological and functional forms of congenital heart defects. Mediator is a multi-subunit complex that plays a central role in transcription initiation by integrating regulatory signals from gene-specific transcriptional activators to RNA polymerase II (Pol II). Recently, Mediator subunit 30 (MED30), a metazoan specific Mediator subunit, has been associated with Langer-Giedion syndrome (LGS) Type II and Cornelia de Lange syndrome-4 (CDLS4), characterized by several abnormalities including congenital heart defects. A point mutation in MED30 has been identified in mouse and is associated with mitochondrial cardiomyopathy. Very recent structural analyses of Mediator revealed that MED30 localizes to the proximal Tail, anchoring Head and Tail modules, thus potentially influencing stability of the Mediator core. However, in vivo cellular and physiological roles of MED30 in maintaining Mediator core integrity remain to be tested. Here, we report that deletion of MED30 in embryonic or adult cardiomyocytes caused rapid development of cardiac defects and lethality. Importantly, cardiomyocyte specific ablation of MED30 destabilized Mediator core subunits, while the kinase module was preserved, demonstrating an essential role of MED30 in stability of the overall Mediator complex. RNAseq analyses of constitutive cardiomyocyte specific Med30 knockout (cKO) embryonic hearts and inducible cardiomyocyte specific Med30 knockout (icKO) adult cardiomyocytes further revealed critical transcription networks in cardiomyocytes controlled by Mediator. Taken together, our results demonstrated that MED30 is essential for Mediator stability and transcriptional networks in both developing and adult cardiomyocytes. Our results affirm the key role of proximal Tail modular subunits in maintaining core Mediator stability in vivo.

CDP-DAG synthase 1 and 2 regulate lipid droplet growth through distinct mechanisms.

Lipid droplets (LDs) are evolutionarily conserved organelles that play critical roles in mammalian lipid storage and metabolism. However, the molecular mechanisms governing the biogenesis and growth of LDs remain poorly understood. Phosphatidic acid (PA) is a precursor of phospholipids and triacylglycerols and substrate of CDP-diacylglycerol (CDP-DAG) synthase 1 (CDS1) and CDS2, which catalyze the formation of CDP-DAG. Here, using siRNA-based gene knockdowns and CRISPR/Cas9-mediated gene knockouts, along with immunological, molecular, and fluorescence microscopy approaches, we examined the role of CDS1 and CDS2 in LD biogenesis and growth. Knockdown of either CDS1 or CDS2 expression resulted in the formation of giant or supersized LDs in cultured mammalian cells. Interestingly, down-regulation of cell death-inducing DFF45-like effector C (CIDEC), encoding a prominent regulator of LD growth in adipocytes, restored LD size in CDS1- but not in CDS2-deficient cells. On the other hand, reducing expression of two enzymes responsible for triacylglycerol synthesis, diacylglycerol O-acyltransferase 2 (DGAT2) and glycerol-3-phosphate acyltransferase 4 (GPAT4), rescued the LD phenotype in CDS2-deficient, but not CDS1-deficient, cells. Moreover, CDS2 deficiency, but not CDS1 deficiency, promoted the LD association of DGAT2 and GPAT4 and impaired initial LD maturation. Finally, although both CDS1 and CDS2 appeared to regulate PA levels on the LD surface, CDS2 had a stronger effect. We conclude that CDS1 and CDS2 regulate LD dynamics through distinct mechanisms.

Phosphatidylserine regulates plasma membrane repair through tetraspanin-enriched macrodomains.

The integrity of the plasma membrane is critical to cell function and survival. Cells have developed multiple mechanisms to repair damaged plasma membranes. A key process during plasma membrane repair is to limit the size of the damage, which is facilitated by the presence of tetraspanin-enriched rings surrounding damage sites. Here, we identify phosphatidylserine-enriched rings surrounding damaged sites of the plasma membrane, resembling tetraspanin-enriched rings. Importantly, the formation of both the phosphatidylserine- and tetraspanin-enriched rings requires phosphatidylserine and its transfer proteins ORP5 and ORP9. Interestingly, ORP9, but not ORP5, is recruited to the damage sites, suggesting cells acquire phosphatidylserine from multiple sources upon plasma membrane damage. We further demonstrate that ORP9 contributes to efficient plasma membrane repair. Our results thus unveil a role for phosphatidylserine and its transfer proteins in facilitating the formation of tetraspanin-enriched macrodomains and plasma membrane repair.

Multi-scale signaling and tumor evolution in high-grade gliomas.

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

Single-cell DNA methylation and 3D genome architecture in the human brain.

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

TMEM41B and VMP1 are phospholipid scramblases.

TMEM41B and VMP1, two endoplasmic reticulum (ER)-resident transmembrane proteins, play important roles in regulating the formation of lipid droplets (LDs), autophagy initiation, and viral infection. However, the biochemical functions of TMEM41B and VMP1 are unclear. A lipids distribution screen suggested TMEM41B and VMP1 are critical to the normal distribution of cholesterol and phosphatidylserine. Biochemical analyses unveiled that TMEM41B and VMP1 have scramblase activity. These findings shed light on the mechanism by which TMEM41B and VMP1 regulate LD formation, lipids distribution, macroautophagy, and viral infection.

Cardiac cell type-specific gene regulatory programs and disease risk association.

Misregulated gene expression in human hearts can result in cardiovascular diseases that are leading causes of mortality worldwide. However, the limited information on the genomic location of candidate cis-regulatory elements (cCREs) such as enhancers and promoters in distinct cardiac cell types has restricted the understanding of these diseases. Here, we defined >287,000 cCREs in the four chambers of the human heart at single-cell resolution, which revealed cCREs and candidate transcription factors associated with cardiac cell types in a region-dependent manner and during heart failure. We further found cardiovascular disease-associated genetic variants enriched within these cCREs including 38 candidate causal atrial fibrillation variants localized to cardiomyocyte cCREs. Additional functional studies revealed that two of these variants affect a cCRE controlling KCNH2/HERG expression and action potential repolarization. Overall, this atlas of human cardiac cCREs provides the foundation for illuminating cell type-specific gene regulation in human hearts during health and disease.