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.

Multimodal single-cell analysis reveals distinct radioresistant stem-like and progenitor cell populations in murine glioma.

Radiation therapy is part of the standard of care for gliomas and kills a subset of tumor cells, while also altering the tumor microenvironment. Tumor cells with stem-like properties preferentially survive radiation and give rise to glioma recurrence. Various techniques for enriching and quantifying cells with stem-like properties have been used, including the fluorescence activated cell sorting (FACS)-based side population (SP) assay, which is a functional assay that enriches for stem-like tumor cells. In these analyses, mouse models of glioma have been used to understand the biology of this disease and therapeutic responses, including the radiation response. We present combined SP analysis and single-cell RNA sequencing of genetically-engineered mouse models of glioma to show a time course of cellular response to radiation. We identify and characterize two distinct tumor cell populations that are inherently radioresistant and also distinct effects of radiation on immune cell populations within the tumor microenvironment.

Formation of chimeric genes by copy-number variation as a mutational mechanism in schizophrenia.

Chimeric genes can be caused by structural genomic rearrangements that fuse together portions of two different genes to create a novel gene. We hypothesize that brain-expressed chimeras may contribute to schizophrenia. Individuals with schizophrenia and control individuals were screened genome wide for copy-number variants (CNVs) that disrupted two genes on the same DNA strand. Candidate events were filtered for predicted brain expression and for frequency < 0.001 in an independent series of 20,000 controls. Four of 124 affected individuals and zero of 290 control individuals harbored such events (p = 0.002); a 47 kb duplication disrupted MATK and ZFR2, a 58 kb duplication disrupted PLEKHD1 and SLC39A9, a 121 kb duplication disrupted DNAJA2 and NETO2, and a 150 kb deletion disrupted MAP3K3 and DDX42. Each fusion produced a stable protein when exogenously expressed in cultured cells. We examined whether these chimeras differed from their parent genes in localization, regulation, or function. Subcellular localizations of DNAJA2-NETO2 and MAP3K3-DDX42 differed from their parent genes. On the basis of the expression profile of the MATK promoter, MATK-ZFR2 is likely to be far more highly expressed in the brain during development than the ZFR2 parent gene. MATK-ZFR2 includes a ZFR2-derived isoform that we demonstrate localizes preferentially to neuronal dendritic branch sites. These results suggest that the formation of chimeric genes is a mechanism by which CNVs contribute to schizophrenia and that, by interfering with parent gene function, chimeras may disrupt critical brain processes, including neurogenesis, neuronal differentiation, and dendritic arborization.

Germline mutations of inhibins in early-onset ovarian epithelial tumors.

To identify novel genetic bases of early-onset epithelial ovarian tumors, we used the trio exome sequencing strategy in a patient without familial history of cancer who presented metastatic serous ovarian adenocarcinomas at 21 years of age. We identified a single de novo mutation (c.1157A>G/p.Asn386Ser) within the INHBA gene encoding the ßA-subunit of inhibins/activins, which play a key role in ovarian development. In vitro, this mutation alters the ratio of secreted activins and inhibins. In a second patient with early-onset serous borderline papillary cystadenoma, we identified an unreported germline mutation (c.179G>T/p.Arg60Leu) of the INHA gene encoding the α-subunit, the partner of the ßA-subunit. This mutation also alters the secreted activin/inhibin ratio, by disrupting both inhibin A and inhibin B biosynthesis. In a cohort of 62 cases, we detected an additional unreported germline mutation of the INHBA gene (c.839G>A/p.Gly280Glu). Our results strongly suggest that inhibin mutations contribute to the genetic determinism of epithelial ovarian tumors.

Garrod's fourth inborn error of metabolism solved by the identification of mutations causing pentosuria.

Pentosuria is one of four conditions hypothesized by Archibald Garrod in 1908 to be inborn errors of metabolism. Mutations responsible for the other three conditions (albinism, alkaptonuria, and cystinuria) have been identified, but the mutations responsible for pentosuria remained unknown. Pentosuria, which affects almost exclusively individuals of Ashkenazi Jewish ancestry, is characterized by high levels of the pentose sugar L-xylulose in blood and urine and deficiency of the enzyme L-xylulose reductase. The condition is autosomal-recessive and completely clinically benign, but in the early and mid-20th century attracted attention because it was often confused with diabetes mellitus and inappropriately treated with insulin. Persons with pentosuria were identified from records of Margaret Lasker, who studied the condition in the 1930s to 1960s. In the DCXR gene encoding L-xylulose reductase, we identified two mutations, DCXR c.583ΔC and DCXR c.52(+1)G > A, each predicted to lead to loss of enzyme activity. Of nine unrelated living pentosuric subjects, six were homozygous for DCXR c.583ΔC, one was homozygous for DCXR c.52(+1)G > A, and two were compound heterozygous for the two mutant alleles. L-xylulose reductase was not detectable in protein lysates from subjects' cells and high levels of xylulose were detected in their sera, confirming the relationship between the DCXR genotypes and the pentosuric phenotype. The combined frequency of the two mutant DCXR alleles in 1,067 Ashkenazi Jewish controls was 0.0173, suggesting a pentosuria frequency of approximately one in 3,300 in this population. Haplotype analysis indicated that the DCXR c.52(+1)G > A mutation arose more recently than the DCXR c.583ΔC mutation.

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.

Context-specific synthetic T cell promoters from assembled transcriptional elements.

Genetic engineering of human lymphocytes for therapeutic applications is constrained by a lack of transgene transcriptional control, resulting in a compromised therapeutic index. Incomplete understanding of transcriptional logic limits the rational design of contextually responsive genetic modules1. Here, we juxtaposed rationally curated transcriptional response element (TRE) oligonucleotides by random concatemerization to generate a library from which we selected context-specific inducible synthetic promoters (iSynPros). Through functional selection, we screened an iSynPro library for "IF-THEN" logic-gated transcriptional responses in human CD8+ T cells expressing a 4-1BB second generation chimeric antigen receptor (CAR). iSynPros exhibiting stringent off-states in quiescent T cells and CAR activation-dependent transcriptional responsiveness were cloned and subjected to TRE composition and pattern analysis, as well as performance in regulating candidate antitumor potency enhancement modules. These data reveal synthetic TRE grammar can mediate logic-gated transgene transcription in human T cells that, when applied to CAR T cell engineering, enhance potency and improve therapeutic indices.

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.

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.

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.

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.

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.