Activity-Induced DNA Breaks Govern the Expression of Neuronal Early-Response Genes.

Neuronal activity causes the rapid expression of immediate early genes that are crucial for experience-driven changes to synapses, learning, and memory. Here, using both molecular and genome-wide next-generation sequencing methods, we report that neuronal activity stimulation triggers the formation of DNA double strand breaks (DSBs) in the promoters of a subset of early-response genes, including Fos, Npas4, and Egr1. Generation of targeted DNA DSBs within Fos and Npas4 promoters is sufficient to induce their expression even in the absence of an external stimulus. Activity-dependent DSB formation is likely mediated by the type II topoisomerase, Topoisomerase IIß (Topo IIß), and knockdown of Topo IIß attenuates both DSB formation and early-response gene expression following neuronal stimulation. Our results suggest that DSB formation is a physiological event that rapidly resolves topological constraints to early-response gene expression in neurons.

Psychiatric risk gene Transcription Factor 4 (TCF4) regulates the density and connectivity of distinct inhibitory interneuron subtypes.

Transcription factor 4 (TCF4) is a basic helix-loop-helix transcription factor that is implicated in a variety of psychiatric disorders including autism spectrum disorder (ASD), major depression, and schizophrenia. Autosomal dominant mutations in TCF4 are causal for a specific ASD called Pitt-Hopkins Syndrome (PTHS). However, our understanding of etiological and pathophysiological mechanisms downstream of TCF4 mutations is incomplete. Single cell sequencing indicates TCF4 is highly expressed in GABAergic interneurons (INs). Here, we performed cell-type specific expression analysis (CSEA) and cellular deconvolution (CD) on bulk RNA sequencing data from 5 different PTHS mouse models. Using CSEA we observed differentially expressed genes (DEGs) were enriched in parvalbumin expressing (PV+) INs and CD predicted a reduction in the PV+ INs population. Therefore, we investigated the role of TCF4 in regulating the development and function of INs in the Tcf4+/tr mouse model of PTHS. In Tcf4+/tr mice, immunohistochemical (IHC) analysis of subtype-specific IN markers and reporter mice identified reductions in PV+, vasoactive intestinal peptide (VIP+), and cortistatin (CST+) expressing INs in the cortex and cholinergic (ChAT+) INs in the striatum, with the somatostatin (SST+) IN population being spared. The reduction of these specific IN populations led to cell-type specific alterations in the balance of excitatory and inhibitory inputs onto PV+ and VIP+ INs and excitatory pyramidal neurons within the cortex. These data indicate TCF4 is a critical regulator of the development of specific subsets of INs and highlight the inhibitory network as an important source of pathophysiology in PTHS.

Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates.

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of COVID-19. The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of ACE2 sequences from 410 vertebrate species, including 252 mammals, to study the conservation of ACE2 and its potential to be used as a receptor by SARS-CoV-2. We designed a five-category binding score based on the conservation properties of 25 amino acids important for the binding between ACE2 and the SARS-CoV-2 spike protein. Only mammals fell into the medium to very high categories and only catarrhine primates into the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 spike protein binding and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (frequency <0.001) variants in 10/25 binding sites. In addition, we found significant signals of selection and accelerated evolution in the ACE2 coding sequence across all mammals, and specific to the bat lineage. Our results, if confirmed by additional experimental data, may lead to the identification of intermediate host species for SARS-CoV-2, guide the selection of animal models of COVID-19, and assist the conservation of animals both in native habitats and in human care.

Chromatin accessibility and microRNA expression in nephron progenitor cells during kidney development.

Mammalian nephrons originate from a population of nephron progenitor cells, and changes in these cells' transcriptomes contribute to the cessation of nephrogenesis, an important determinant of nephron number. To characterize microRNA (miRNA) expression and identify putative cis-regulatory regions, we collected nephron progenitor cells from mouse kidneys at embryonic day 14.5 and postnatal day zero and assayed small RNA expression and transposase-accessible chromatin. We detect expression of 1104 miRNA (114 with expression changes), and 46,374 chromatin accessible regions (2103 with changes in accessibility). Genome-wide, our data highlight processes like cellular differentiation, cell migration, extracellular matrix interactions, and developmental signaling pathways. Furthermore, they identify new candidate cis-regulatory elements for Eya1 and Pax8, both genes with a role in nephron progenitor cell differentiation. Finally, we associate expression-changing miRNAs, including let-7-5p, miR-125b-5p, miR-181a-2-3p, and miR-9-3p, with candidate cis-regulatory elements and target genes. These analyses highlight new putative cis-regulatory loci for miRNA in nephron progenitors.

Addiction-Associated Genetic Variants Implicate Brain Cell Type- and Region-Specific Cis-Regulatory Elements in Addiction Neurobiology.

Recent large genome-wide association studies have identified multiple confident risk loci linked to addiction-associated behavioral traits. Most genetic variants linked to addiction-associated traits lie in noncoding regions of the genome, likely disrupting cis-regulatory element (CRE) function. CREs tend to be highly cell type-specific and may contribute to the functional development of the neural circuits underlying addiction. Yet, a systematic approach for predicting the impact of risk variants on the CREs of specific cell populations is lacking. To dissect the cell types and brain regions underlying addiction-associated traits, we applied stratified linkage disequilibrium score regression to compare genome-wide association studies to genomic regions collected from human and mouse assays for open chromatin, which is associated with CRE activity. We found enrichment of addiction-associated variants in putative CREs marked by open chromatin in neuronal (NeuN+) nuclei collected from multiple prefrontal cortical areas and striatal regions known to play major roles in reward and addiction. To further dissect the cell type-specific basis of addiction-associated traits, we also identified enrichments in human orthologs of open chromatin regions of female and male mouse neuronal subtypes: cortical excitatory, D1, D2, and PV. Last, we developed machine learning models to predict mouse cell type-specific open chromatin, enabling us to further categorize human NeuN+ open chromatin regions into cortical excitatory or striatal D1 and D2 neurons and predict the functional impact of addiction-associated genetic variants. Our results suggest that different neuronal subtypes within the reward system play distinct roles in the variety of traits that contribute to addiction.SIGNIFICANCE STATEMENT We combine statistical genetic and machine learning techniques to find that the predisposition to for nicotine, alcohol, and cannabis use behaviors can be partially explained by genetic variants in conserved regulatory elements within specific brain regions and neuronal subtypes of the reward system. Our computational framework can flexibly integrate open chromatin data across species to screen for putative causal variants in a cell type- and tissue-specific manner for numerous complex traits.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin.

BACKGROUND: Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high. RESULTS: We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin. CONCLUSION: The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

Cell Type-Specific Oxidative Stress Genomic Signatures in the Globus Pallidus of Dopamine-Depleted Mice.

Neuron subtype dysfunction is a key contributor to neurologic disease circuits, but identifying associated gene regulatory pathways is complicated by the molecular complexity of the brain. For example, parvalbumin-expressing (PV+) neurons in the external globus pallidus (GPe) are critically involved in the motor deficits of dopamine-depleted mouse models of Parkinson's disease, where cell type-specific optogenetic stimulation of PV+ neurons over other neuron populations rescues locomotion. Despite the distinct roles these cell types play in the neural circuit, the molecular correlates remain unknown because of the difficulty of isolating rare neuron subtypes. To address this issue, we developed a new viral affinity purification strategy, Cre-Specific Nuclear Anchored Independent Labeling, to isolate Cre recombinase-expressing (Cre+) nuclei from the adult mouse brain. Applying this technology, we performed targeted assessments of the cell type-specific transcriptomic and epigenetic effects of dopamine depletion on PV+ and PV- cells within three brain regions of male and female mice: GPe, striatum, and cortex. We found GPe PV+ neuron-specific gene expression changes that suggested increased hypoxia-inducible factor 2α signaling. Consistent with transcriptomic data, regions of open chromatin affected by dopamine depletion within GPe PV+ neurons were enriched for hypoxia-inducible factor family binding motifs. The gene expression and epigenomic experiments performed on PV+ neurons isolated by Cre-Specific Nuclear Anchored Independent Labeling identified a transcriptional regulatory network mediated by the neuroprotective factor Hif2a as underlying neural circuit differences in response to dopamine depletion.SIGNIFICANCE STATEMENT Cre-Specific Nuclear Anchored Independent Labeling is an enhanced, virus-based approach to isolate nuclei of a specific cell type for transcriptome and epigenome interrogation that decreases dependency on transgenic animals. Applying this technology to GPe parvalbumin-expressing neurons in a mouse model of Parkinson's disease, we discovered evidence for an upregulation of the oxygen homeostasis maintaining pathway involving Hypoxia-inducible factor 2α. These results provide new insight into how neuron subtypes outside the substantia nigra pars compacta may be compensating at a molecular level for differences in the motor production neural circuit during the progression of Parkinson's disease. Furthermore, they emphasize the utility of cell type-specific technologies, such as Cre-Specific Nuclear Anchored Independent Labeling, for isolated assessment of specific neuron subtypes in complex systems.

HALPER facilitates the identification of regulatory element orthologs across species.

SUMMARY: Diverse traits have evolved through cis-regulatory changes in genome sequence that influence the magnitude, timing and cell type-specificity of gene expression. Advances in high-throughput sequencing and regulatory genomics have led to the identification of regulatory elements in individual species, but these genomic regions remain difficult to align across taxonomic orders due to their lack of sequence conservation relative to protein coding genes. The groundwork for tracing the evolution of regulatory elements is provided by the recent assembly of hundreds of genomes, the generation of reference-free Cactus multiple sequence alignments of these genomes, and the development of the halLiftover tool for mapping regions across these alignments. We present halLiftover Post-processing for the Evolution of Regulatory Elements (HALPER), a tool for constructing contiguous regulatory element orthologs from the outputs of halLiftover. We anticipate that this tool will enable users to efficiently identify orthologs of regulatory elements across hundreds of species, providing novel insights into the evolution of traits that have evolved through gene expression. AVAILABILITY AND IMPLEMENTATION: HALPER is implemented in python and available on github: https://github.com/pfenninglab/halLiftover-postprocessing. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Conserved epigenomic signals in mice and humans reveal immune basis of Alzheimer's disease.

Alzheimer's disease (AD) is a severe age-related neurodegenerative disorder characterized by accumulation of amyloid-ß plaques and neurofibrillary tangles, synaptic and neuronal loss, and cognitive decline. Several genes have been implicated in AD, but chromatin state alterations during neurodegeneration remain uncharacterized. Here we profile transcriptional and chromatin state dynamics across early and late pathology in the hippocampus of an inducible mouse model of AD-like neurodegeneration. We find a coordinated downregulation of synaptic plasticity genes and regulatory regions, and upregulation of immune response genes and regulatory regions, which are targeted by factors that belong to the ETS family of transcriptional regulators, including PU.1. Human regions orthologous to increasing-level enhancers show immune-cell-specific enhancer signatures as well as immune cell expression quantitative trait loci, while decreasing-level enhancer orthologues show fetal-brain-specific enhancer activity. Notably, AD-associated genetic variants are specifically enriched in increasing-level enhancer orthologues, implicating immune processes in AD predisposition. Indeed, increasing enhancers overlap known AD loci lacking protein-altering variants, and implicate additional loci that do not reach genome-wide significance. Our results reveal new insights into the mechanisms of neurodegeneration and establish the mouse as a useful model for functional studies of AD regulatory regions.

Integrative analysis of 111 reference human epigenomes.

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

High-throughput functional comparison of promoter and enhancer activities.

Promoters initiate RNA synthesis, and enhancers stimulate promoter activity. Whether promoter and enhancer activities are encoded distinctly in DNA sequences is unknown. We measured the enhancer and promoter activities of thousands of DNA fragments transduced into mouse neurons. We focused on genomic loci bound by the neuronal activity-regulated coactivator CREBBP, and we measured enhancer and promoter activities both before and after neuronal activation. We find that the same sequences typically encode both enhancer and promoter activities. However, gene promoters generate more promoter activity than distal enhancers, despite generating similar enhancer activity. Surprisingly, the greater promoter activity of gene promoters is not due to conventional core promoter elements or splicing signals. Instead, we find that particular transcription factor binding motifs are intrinsically biased toward the generation of promoter activity, whereas others are not. Although the specific biases we observe may be dependent on experimental or cellular context, our results suggest that gene promoters are distinguished from distal enhancers by specific complements of transcriptional activators.

The transcription factor calcium-response factor limits NMDA receptor-dependent transcription in the developing brain.

Neuronal activity sculpts brain development by inducing the transcription of genes such as brain-derived neurotrophic factor (Bdnf) that modulate the function of synapses. Sensory experience is transduced into changes in gene transcription via the activation of calcium signaling pathways downstream of both L-type voltage-gated calcium channels (L-VGCCs) and NMDA-type glutamate receptors (NMDARs). These signaling pathways converge on the regulation of transcription factors including calcium-response factor (CaRF). Although CaRF is dispensable for the transcriptional induction of Bdnf following the activation of L-VGCCs, here we show that the loss of CaRF leads to enhanced NMDAR-dependent transcription of Bdnf as well as Arc. We identify the NMDAR subunit-encoding gene Grin3a as a regulatory target of CaRF, and we show that expression of both Carf and Grin3a is depressed by the elevation of intracellular calcium, linking the function of this transcriptional regulatory pathway to neuronal activity. We find that light-dependent activation of Bdnf and Arc transcription is enhanced in the visual cortex of young CaRF knockout mice, suggesting a role for CaRF-dependent dampening of NMDAR-dependent transcription in the developing brain. Finally, we demonstrate that enhanced Bdnf expression in CaRF-lacking neurons increases inhibitory synapse formation. Taken together, these data reveal a novel role for CaRF as an upstream regulator of NMDAR-dependent gene transcription and synapse formation in the developing brain. NMDARs promote brain development by inducing the transcription of genes, including brain-derived neurotrophic factor (BDNF). We show that the transcription factor calcium-response factor (CaRF) limits NMDAR-dependent BDNF induction by regulating expression of the NMDAR subunit GluN3A. Loss of CaRF leads to enhanced BDNF-dependent GABAergic synapse formation indicating the importance of this process for brain development. Our observation that both CaRF and GluN3A are down-regulated by intracellular calcium suggests that this may be a mechanism for experience-dependent modulation of synapse formation.

The genome of a songbird.

The zebra finch is an important model organism in several fields with unique relevance to human neuroscience. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken-the only bird with a sequenced genome until now. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.

Relation of addiction genes to hypothalamic gene changes subserving genesis and gratification of a classic instinct, sodium appetite.

Sodium appetite is an instinct that involves avid specific intention. It is elicited by sodium deficiency, stress-evoked adrenocorticotropic hormone (ACTH), and reproduction. Genome-wide microarrays in sodium-deficient mice or after ACTH infusion showed up-regulation of hypothalamic genes, including dopamine- and cAMP-regulated neuronal phosphoprotein 32 kDa (DARPP-32), dopamine receptors-1 and -2, α-2C- adrenoceptor, and striatally enriched protein tyrosine phosphatase (STEP). Both DARPP-32 and neural plasticity regulator activity-regulated cytoskeleton associated protein (ARC) were up-regulated in lateral hypothalamic orexinergic neurons by sodium deficiency. Administration of dopamine D1 (SCH23390) and D2 receptor (raclopride) antagonists reduced gratification of sodium appetite triggered by sodium deficiency. SCH23390 was specific, having no effect on osmotic-induced water drinking, whereas raclopride also reduced water intake. D1 receptor KO mice had normal sodium appetite, indicating compensatory regulation. Appetite was insensitive to SCH23390, confirming the absence of off-target effects. Bilateral microinjection of SCH23390 (100 nM in 200 nL) into rats' lateral hypothalamus greatly reduced sodium appetite. Gene set enrichment analysis in hypothalami of mice with sodium appetite showed significant enrichment of gene sets previously linked to addiction (opiates and cocaine). This finding of concerted gene regulation was attenuated on gratification with perplexingly rapid kinetics of only 10 min, anteceding significant absorption of salt from the gut. Salt appetite and hedonic liking of salt taste have evolved over >100 million y (e.g., being present in Metatheria). Drugs causing pleasure and addiction are comparatively recent and likely reflect usurping of evolutionary ancient systems with high survival value by the gratification of contemporary hedonic indulgences. Our findings outline a molecular logic for instinctive behavior encoded by the brain with possible important translational-medical implications.

Single nuclei transcriptomics in human and non-human primate striatum in opioid use disorder.

In brain, the striatum is a heterogenous region involved in reward and goal-directed behaviors. Striatal dysfunction is linked to psychiatric disorders, including opioid use disorder (OUD). Striatal subregions are divided based on neuroanatomy, each with unique roles in OUD. In OUD, the dorsal striatum is involved in altered reward processing, formation of habits, and development of negative affect during withdrawal. Using single nuclei RNA-sequencing, we identified both canonical (e.g., dopamine receptor subtype) and less abundant cell populations (e.g., interneurons) in human dorsal striatum. Pathways related to neurodegeneration, interferon response, and DNA damage were significantly enriched in striatal neurons of individuals with OUD. DNA damage markers were also elevated in striatal neurons of opioid-exposed rhesus macaques. Sex-specific molecular differences in glial cell subtypes associated with chronic stress were found in OUD, particularly female individuals. Together, we describe different cell types in human dorsal striatum and identify cell type-specific alterations in OUD.

Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements.

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Integrative multi-dimensional characterization of striatal projection neuron heterogeneity in adult brain.

Striatal projection neurons (SPNs) are traditionally segregated into two subpopulations expressing dopamine (DA) D1-like or D2-like receptors. However, this dichotomy is challenged by recent evidence. Functional and expression studies raise important questions: do SPNs co-express different DA receptors, and do these differences reflect unique striatal spatial distributions and expression profiles? Using RNAscope in mouse striatum, we report heterogenous SPN subpopulations distributed across dorsal-ventral and rostral-caudal axes. SPN subpopulations co-express multiple DA receptors, including D1 and D2 (D1/2R) and D1 and D3. Our integrative approach using single-nuclei multi-omics analyses provides a simple consensus to describe SPNs across diverse datasets, connecting it to complementary spatial mapping. Combining RNAscope and multi-omics shows D1/2R SPNs further separate into distinct subtypes according to spatial organization and conserved marker genes. Each SPN cell type contributes uniquely to genetic risk for neuropsychiatric diseases. Our results bridge anatomy and transcriptomics to offer new understandings of striatal neuron heterogeneity.

Relating enhancer genetic variation across mammals to complex phenotypes using machine learning.

Protein-coding differences between species often fail to explain phenotypic diversity, suggesting the involvement of genomic elements that regulate gene expression such as enhancers. Identifying associations between enhancers and phenotypes is challenging because enhancer activity can be tissue-dependent and functionally conserved despite low sequence conservation. We developed the Tissue-Aware Conservation Inference Toolkit (TACIT) to associate candidate enhancers with species' phenotypes using predictions from machine learning models trained on specific tissues. Applying TACIT to associate motor cortex and parvalbumin-positive interneuron enhancers with neurological phenotypes revealed dozens of enhancer-phenotype associations, including brain size-associated enhancers that interact with genes implicated in microcephaly or macrocephaly. TACIT provides a foundation for identifying enhancers associated with the evolution of any convergently evolved phenotype in any large group of species with aligned genomes.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease.

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.