Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Single-cell census of human tooth development enables generation of human enamel

Tooth enamel secreted by ameloblasts (AMs) is the hardest material in the human body, acting as a shield to protect the teeth. However, the enamel is gradually damaged or partially lost in over 90% of adults and cannot be regenerated due to a lack of ameloblasts in erupted teeth. Here, we use single-cell combinatorial indexing RNA sequencing (sci-RNA-seq) to establish a spatiotemporal single-cell census for the developing human tooth and identify regulatory mechanisms controlling the differentiation process of human ameloblasts. We identify key signaling pathways involved between the support cells and ameloblasts during fetal development and recapitulate those findings in human ameloblast in vitro differentiation from induced pluripotent stem cells (iPSCs). We furthermore develop a disease model of amelogenesis imperfecta in a three-dimensional (3D) organoid system and show AM maturation to mineralized structure in vivo. These studies pave the way for future regenerative dentistry.

A single-cell multi-omic atlas spanning the adult rhesus macaque brain

Cataloging the diverse cellular architecture of the primate brain is crucial for understanding cognition, behavior, and disease in humans. Here, we generated a brain-wide single-cell multimodal molecular atlas of the rhesus macaque brain. Together, we profiled 2.58 M transcriptomes and 1.59 M epigenomes from single nuclei sampled from 30 regions across the adult brain. Cell composition differed extensively across the brain, revealing cellular signatures of region-specific functions. We also identified 1.19 M candidate regulatory elements, many previously unidentified, allowing us to explore the landscape of cis-regulatory grammar and neurological disease risk in a cell type-specific manner. Altogether, this multi-omic atlas provides an open resource for investigating the evolution of the human brain and identifying novel targets for disease interventions.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Multiregion transcriptomic profiling of the primate brain reveals signatures of aging and the social environment

Aging is accompanied by a host of social and biological changes that correlate with behavior, cognitive health and susceptibility to neurodegenerative disease. To understand trajectories of brain aging in a primate, we generated a multiregion bulk (N = 527 samples) and single-nucleus (N = 24 samples) brain transcriptional dataset encompassing 15 brain regions and both sexes in a unique population of free-ranging, behaviorally phenotyped rhesus macaques. We demonstrate that age-related changes in the level and variance of gene expression occur in genes associated with neural functions and neurological diseases, including Alzheimer's disease. Further, we show that higher social status in females is associated with younger relative transcriptional ages, providing a link between the social environment and aging in the brain. Our findings lend insight into biological mechanisms underlying brain aging in a nonhuman primate model of human behavior, cognition and health.

Genome-wide fetalization of enhancer architecture in heart disease

Heart disease is associated with re-expression of key transcription factors normally active only during prenatal development of the heart. However, the impact of this reactivation on the regulatory landscape in heart disease is unclear. Here, we use RNA-seq and ChIP-seq targeting a histone modification associated with active transcriptional enhancers to generate genome-wide enhancer maps from left ventricle tissue from up to 26 healthy controls, 18 individuals with idiopathic dilated cardiomyopathy (DCM), and five fetal hearts. Healthy individuals have a highly reproducible epigenomic landscape, consisting of more than 33,000 predicted heart enhancers. In contrast, we observe reproducible disease-associated changes in activity at 6,850 predicted heart enhancers. Combined analysis of adult and fetal samples reveals that the heart disease epigenome and transcriptome both acquire fetal-like characteristics, with 3,400 individual enhancers sharing fetal regulatory properties. We also provide a comprehensive data resource (http://heart.lbl.gov) for the mechanistic exploration of DCM etiology.

The Ties That Bind: Mapping the Dynamic Enhancer-Promoter Interactome

Coupling chromosome conformation capture to molecular enrichment for promoter-containing DNA fragments enables the systematic mapping of interactions between individual distal regulatory sequences and their target genes. In this Minireview, we describe recent progress in the application of this technique and related complementary approaches to gain insight into the lineage- and cell-type-specific dynamics of interactions between regulators and gene promoters.

Genome-wide compendium and functional assessment of in vivo heart enhancers

Whole-genome sequencing is identifying growing numbers of non-coding variants in human disease studies, but the lack of accurate functional annotations prevents their interpretation. We describe the genome-wide landscape of distant-acting enhancers active in the developing and adult human heart, an organ whose impairment is a predominant cause of mortality and morbidity. Using integrative analysis of >35 epigenomic data sets from mouse and human pre- and postnatal hearts we created a comprehensive reference of >80,000 putative human heart enhancers. To illustrate the importance of enhancers in the regulation of genes involved in heart disease, we deleted the mouse orthologs of two human enhancers near cardiac myosin genes. In both cases, we observe in vivo expression changes and cardiac phenotypes consistent with human heart disease. Our study provides a comprehensive catalogue of human heart enhancers for use in clinical whole-genome sequencing studies and highlights the importance of enhancers for cardiac function.

Identification of non-coding variants has outstripped our ability to annotate and interpret them. Dickel et al. present a compendium of over 80,000 putative human heart enhancers and demonstrate that two conserved enhancers are required for proper cardiac function in mice.

Genetic and functional analyses identify DISC1 as a novel callosal agenesis candidate gene

Agenesis of the corpus callosum (AgCC) is a congenital brain malformation that occurs in approximately 1:1,000-1:6,000 births. Several syndromes associated with AgCC have been traced to single gene mutations; however, the majority of AgCC causes remain unidentified. We investigated a mother and two children who all shared complete AgCC and a chromosomal deletion at 1q42. We fine mapped this deletion and show that it includes Disrupted-in-Schizophrenia 1 (DISC1), a gene implicated in schizophrenia and other psychiatric disorders. Furthermore, we report a de novo chromosomal deletion at 1q42.13 to q44, which includes DISC1, in another individual with AgCC. We resequenced DISC1 in a cohort of 144 well-characterized AgCC individuals and identified 20 sequence changes, of which 4 are rare potentially pathogenic variants. Two of these variants were undetected in 768 control chromosomes. One of these is a splice site mutation at the 5' boundary of exon 11 that dramatically reduces full-length mRNA expression of DISC1, but not of shorter forms. We investigated the developmental expression of mouse DISC1 and find that it is highly expressed in the embryonic corpus callosum at a critical time for callosal formation. Taken together our results suggest a significant role for DISC1 in corpus callosum development.