Using a comprehensive atlas and predictive models to reveal the complexity and evolution of brain-active regulatory elements.

Most genetic variants associated with psychiatric disorders are located in noncoding regions of the genome. To investigate their functional implications, we integrate epigenetic data from the PsychENCODE Consortium and other published sources to construct a comprehensive atlas of candidate brain cis-regulatory elements. Using deep learning, we model these elements' sequence syntax and predict how binding sites for lineage-specific transcription factors contribute to cell type-specific gene regulation in various types of glia and neurons. The elements' evolutionary history suggests that new regulatory information in the brain emerges primarily via smaller sequence mutations within conserved mammalian elements rather than entirely new human- or primate-specific sequences. However, primate-specific candidate elements, particularly those active during fetal brain development and in excitatory neurons and astrocytes, are implicated in the heritability of brain-related human traits. Additionally, we introduce PsychSCREEN, a web-based platform offering interactive visualization of PsychENCODE-generated genetic and epigenetic data from diverse brain cell types in individuals with psychiatric disorders and healthy controls.

Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Multicenter integrated analysis of noncoding CRISPRi screens.

The ENCODE Consortium's efforts to annotate noncoding cis-regulatory elements (CREs) have advanced our understanding of gene regulatory landscapes. Pooled, noncoding CRISPR screens offer a systematic approach to investigate cis-regulatory mechanisms. The ENCODE4 Functional Characterization Centers conducted 108 screens in human cell lines, comprising >540,000 perturbations across 24.85 megabases of the genome. Using 332 functionally confirmed CRE-gene links in K562 cells, we established guidelines for screening endogenous noncoding elements with CRISPR interference (CRISPRi), including accurate detection of CREs that exhibit variable, often low, transcriptional effects. Benchmarking five screen analysis tools, we find that CASA produces the most conservative CRE calls and is robust to artifacts of low-specificity single guide RNAs. We uncover a subtle DNA strand bias for CRISPRi in transcribed regions with implications for screen design and analysis. Together, we provide an accessible data resource, predesigned single guide RNAs for targeting 3,275,697 ENCODE SCREEN candidate CREs with CRISPRi and screening guidelines to accelerate functional characterization of the noncoding genome.

An encyclopedia of enhancer-gene regulatory interactions in the human genome.

Identifying transcriptional enhancers and their target genes is essential for understanding gene regulation and the impact of human genetic variation on disease1-6. Here we create and evaluate a resource of >13 million enhancer-gene regulatory interactions across 352 cell types and tissues, by integrating predictive models, measurements of chromatin state and 3D contacts, and largescale genetic perturbations generated by the ENCODE Consortium7. We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,411 elementgene pairs measured in CRISPR perturbation experiments, >30,000 fine-mapped eQTLs, and 569 fine-mapped GWAS variants linked to a likely causal gene. Using this framework, we develop a new predictive model, ENCODE-rE2G, that achieves state-of-the-art performance across multiple prediction tasks, demonstrating a strategy involving iterative perturbations and supervised machine learning to build increasingly accurate predictive models of enhancer regulation. Using the ENCODE-rE2G model, we build an encyclopedia of enhancer-gene regulatory interactions in the human genome, which reveals global properties of enhancer networks, identifies differences in the functions of genes that have more or less complex regulatory landscapes, and improves analyses to link noncoding variants to target genes and cell types for common, complex diseases. By interpreting the model, we find evidence that, beyond enhancer activity and 3D enhancer-promoter contacts, additional features guide enhancerpromoter communication including promoter class and enhancer-enhancer synergy. Altogether, these genome-wide maps of enhancer-gene regulatory interactions, benchmarking software, predictive models, and insights about enhancer function provide a valuable resource for future studies of gene regulation and human genetics.

The ENCODE4 long-read RNA-seq collection reveals distinct classes of transcript structure diversity.

The majority of mammalian genes encode multiple transcript isoforms that result from differential promoter use, changes in exonic splicing, and alternative 3' end choice. Detecting and quantifying transcript isoforms across tissues, cell types, and species has been extremely challenging because transcripts are much longer than the short reads normally used for RNA-seq. By contrast, long-read RNA-seq (LR-RNA-seq) gives the complete structure of most transcripts. We sequenced 264 LR-RNA-seq PacBio libraries totaling over 1 billion circular consensus reads (CCS) for 81 unique human and mouse samples. We detect at least one full-length transcript from 87.7% of annotated human protein coding genes and a total of 200,000 full-length transcripts, 40% of which have novel exon junction chains. To capture and compute on the three sources of transcript structure diversity, we introduce a gene and transcript annotation framework that uses triplets representing the transcript start site, exon junction chain, and transcript end site of each transcript. Using triplets in a simplex representation demonstrates how promoter selection, splice pattern, and 3' processing are deployed across human tissues, with nearly half of multi-transcript protein coding genes showing a clear bias toward one of the three diversity mechanisms. Evaluated across samples, the predominantly expressed transcript changes for 74% of protein coding genes. In evolution, the human and mouse transcriptomes are globally similar in types of transcript structure diversity, yet among individual orthologous gene pairs, more than half (57.8%) show substantial differences in mechanism of diversification in matching tissues. This initial large-scale survey of human and mouse long-read transcriptomes provides a foundation for further analyses of alternative transcript usage, and is complemented by short-read and microRNA data on the same samples and by epigenome data elsewhere in the ENCODE4 collection.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites.

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

Downregulation of Hsp90 and the antimicrobial peptide Mtk suppresses poly(GR)-induced neurotoxicity in C9ORF72-ALS/FTD.

GGGGCC repeat expansion in the C9ORF72 gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Repeat RNAs can be translated into dipeptide repeat proteins, including poly(GR), whose mechanisms of action remain largely unknown. In an RNA-seq analysis of poly(GR) toxicity in Drosophila, we found that several antimicrobial peptide genes, such as metchnikowin (Mtk), and heat shock protein (Hsp) genes are activated. Mtk knockdown in the fly eye or in all neurons suppresses poly(GR) neurotoxicity. These findings suggest a cell-autonomous role of Mtk in neurodegeneration. Hsp90 knockdown partially rescues both poly(GR) toxicity in flies and neurodegeneration in C9ORF72 motor neurons derived from induced pluripotent stem cells (iPSCs). Topoisomerase II (TopoII) regulates poly(GR)-induced upregulation of Hsp90 and Mtk. TopoII knockdown also suppresses poly(GR) toxicity in Drosophila and improves survival of C9ORF72 iPSC-derived motor neurons. These results suggest potential novel therapeutic targets for C9ORF72-ALS/FTD.

Therapeutic Trial of anle138b in Mouse Models of Genetic Prion Disease.

Phenotypic screening has yielded small-molecule inhibitors of prion replication that are effective in vivo against certain prion strains but not others. Here, we sought to test the small molecule anle138b in multiple mouse models of prion disease. In mice inoculated with the RML strain of prions, anle138b doubled survival and durably suppressed astrogliosis measured by live-animal bioluminescence imaging. In knock-in mouse models of the D178N and E200K mutations that cause genetic prion disease, however, we were unable to identify a clear, quantifiable disease endpoint against which to measure therapeutic efficacy. Among untreated animals, the mutations did not impact overall survival, and bioluminescence remained low out to >20 months of age. Vacuolization and PrP deposition were observed in some brain regions in a subset of mutant animals but appeared to be unable to carry the weight of a primary endpoint in a therapeutic study. We conclude that not all animal models of prion disease are suited to well-powered therapeutic efficacy studies, and care should be taken in choosing the models that will support drug development programs. IMPORTANCE There is an urgent need to develop drugs for prion disease, a currently untreatable neurodegenerative disease. In this effort, there is a debate over which animal models can best support a drug development program. While the study of prion disease benefits from excellent animal models because prions naturally afflict many different mammals, different models have different capabilities and limitations. Here, we conducted a therapeutic efficacy study of the drug candidate anle138b in mouse models with two of the most common mutations that cause genetic prion disease. In a more typical model where prions are injected directly into the brain, we found anle138b to be effective. In the genetic models, however, the animals never reached a clear, measurable point of disease onset. We conclude that not all prion disease animal models are ideally suited to drug efficacy studies, and well-defined, quantitative disease metrics should be a priority.

FAVOR: functional annotation of variants online resource and annotator for variation across the human genome.

Large biobank-scale whole genome sequencing (WGS) studies are rapidly identifying a multitude of coding and non-coding variants. They provide an unprecedented resource for illuminating the genetic basis of human diseases. Variant functional annotations play a critical role in WGS analysis, result interpretation, and prioritization of disease- or trait-associated causal variants. Existing functional annotation databases have limited scope to perform online queries and functionally annotate the genotype data of large biobank-scale WGS studies. We develop the Functional Annotation of Variants Online Resources (FAVOR) to meet these pressing needs. FAVOR provides a comprehensive multi-faceted variant functional annotation online portal that summarizes and visualizes findings of all possible nine billion single nucleotide variants (SNVs) across the genome. It allows for rapid variant-, gene- and region-level queries of variant functional annotations. FAVOR integrates variant functional information from multiple sources to describe the functional characteristics of variants and facilitates prioritizing plausible causal variants influencing human phenotypes. Furthermore, we provide a scalable annotation tool, FAVORannotator, to functionally annotate large-scale WGS studies and efficiently store the genotype and their variant functional annotation data in a single file using the annotated Genomic Data Structure (aGDS) format, making downstream analysis more convenient. FAVOR and FAVORannotator are available at https://favor.genohub.org.

Building integrative functional maps of gene regulation.

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

Recombinant Aspergillus fumigatus antigens Asp f 3 and Asp f 9 in liposomal vaccine protect mice against invasive pulmonary aspergillosis.

Invasive pulmonary aspergillosis caused by the ubiquitous mold Aspergillus fumigatus is a major threat to immunocompromised patients, causing unacceptably high mortality despite standard of care treatment, and costing an estimated $1.2 billion annually. Treatment for this disease has been complicated by the emergence of azole resistant strains of A. fumigatus, rendering first-line antifungal therapy ineffective. The difficulties in treating infected patients using currently available drugs make immunotherapeutic vaccination an attractive option. Here, we demonstrate the efficacy of VesiVax® adjuvant liposomes, consisting of a combination of two individual liposome preparations, to which two recombinant A. fumigatus surface antigens, Asp f 3 and Asp f 9 (VesiVax® Af3/9), have been chemically conjugated. Using a murine model, we demonstrate that VesiVax® Af3/9 is protective against infection by azole resistant strains of A. fumigatus in both steroid-suppressed and neutropenic mice as quantified by improved survival and reduced fungal burden in the lungs. This protection correlates with upregulation of IL-4 produced by splenocytes, and the presence of Asp f 3 and Asp f 9 specific IgG2a antibodies in the serum of mice given VesiVax® Af3/9. Furthermore, mice given VesiVax® Af3/9 with a subsequent course of liposomal amphotericin B (AmBisome®) had improved survival over those given either treatment alone, indicating a benefit to VesiVax® Af3/9 vaccination even in the case of infections that require follow-up antifungal treatment. These data demonstrate that prophylactic vaccination with VesiVax® Af3/9 is a promising method of protection against invasive pulmonary aspergillosis even as the changing face of the disease renders current therapies ineffective.

Understanding Nonsuicidal Self-Injury.

ABSTRACT: Nonsuicidal self-injury has existed since biblical times. It is increasing and gaining acceptance particularly among adolescents. Spiritual discontent or doubt is often a component. Therapeutic interventions include cognitive and behavior practices and developing strong, sustained relationships with patients to encourage healthy transformation.

Regional variability and genotypic and pharmacodynamic effects on PrP concentration in the CNS.

Prion protein (PrP) concentration controls the kinetics of prion replication and is a genetically and pharmacologically validated therapeutic target for prion disease. In order to evaluate PrP concentration as a pharmacodynamic biomarker and assess its contribution to known prion disease risk factors, we developed and validated a plate-based immunoassay reactive for PrP across 6 species of interest and applicable to brain and cerebrospinal fluid (CSF). PrP concentration varied dramatically across different brain regions in mice, cynomolgus macaques, and humans. PrP expression did not appear to contribute to the known risk factors of age, sex, or common PRNP genetic variants. CSF PrP was lowered in the presence of rare pathogenic PRNP variants, with heterozygous carriers of P102L displaying 55%, and D178N just 31%, of the CSF PrP concentration of mutation-negative controls. In rodents, pharmacologic reduction of brain Prnp RNA was reflected in brain parenchyma PrP and, in turn in CSF PrP, validating CSF as a sampling compartment for the effect of PrP-lowering therapy. Our findings support the use of CSF PrP as a pharmacodynamic biomarker for PrP-lowering drugs and suggest that relative reduction from individual baseline CSF PrP concentration may be an appropriate marker for target engagement.

Factorbook: an updated catalog of transcription factor motifs and candidate regulatory motif sites.

The human genome contains ∼2000 transcriptional regulatory proteins, including ∼1600 DNA-binding transcription factors (TFs) recognizing characteristic sequence motifs to exert regulatory effects on gene expression. The binding specificities of these factors have been profiled both in vitro, using techniques such as HT-SELEX, and in vivo, using techniques including ChIP-seq. We previously developed Factorbook, a TF-centric database of annotations, motifs, and integrative analyses based on ChIP-seq data from Phase II of the ENCODE Project. Here we present an update to Factorbook which significantly expands the breadth of cell type and TF coverage. The update includes an expanded motif catalog derived from thousands of ENCODE Phase II and III ChIP-seq experiments and HT-SELEX experiments; this motif catalog is integrated with the ENCODE registry of candidate cis-regulatory elements to annotate a comprehensive collection of genome-wide candidate TF binding sites. The database also offers novel tools for applying the motif models within machine learning frameworks and using these models for integrative analysis, including annotation of variants and disease and trait heritability. Factorbook is publicly available at www.factorbook.org; we will continue to expand the resource as ENCODE Phase IV data are released.

Integration of high-resolution promoter profiling assays reveals novel, cell type-specific transcription start sites across 115 human cell and tissue types.

Accurate transcription start site (TSS) annotations are essential for understanding transcriptional regulation and its role in human disease. Gene collections such as GENCODE contain annotations for tens of thousands of TSSs, but not all of these annotations are experimentally validated nor do they contain information on cell type-specific usage. Therefore, we sought to generate a collection of experimentally validated TSSs by integrating RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression (RAMPAGE) data from 115 cell and tissue types, which resulted in a collection of approximately 50 thousand representative RAMPAGE peaks. These peaks are primarily proximal to GENCODE-annotated TSSs and are concordant with other transcription assays. Because RAMPAGE uses paired-end reads, we were then able to connect peaks to transcripts by analyzing the genomic positions of the 3' ends of read mates. Using this paired-end information, we classified the vast majority (37 thousand) of our RAMPAGE peaks as verified TSSs, updating TSS annotations for 20% of GENCODE genes. We also found that these updated TSS annotations are supported by epigenomic and other transcriptomic data sets. To show the utility of this RAMPAGE rPeak collection, we intersected it with the NHGRI/EBI genome-wide association study (GWAS) catalog and identified new candidate GWAS genes. Overall, our work shows the importance of integrating experimental data to further refine TSS annotations and provides a valuable resource for the biological community.

Genetic and epigenetic features of promoters with ubiquitous chromatin accessibility support ubiquitous transcription of cell-essential genes.

Gene expression is controlled by regulatory elements within accessible chromatin. Although most regulatory elements are cell type-specific, a subset is accessible in nearly all the 517 human and 94 mouse cell and tissue types assayed by the ENCODE consortium. We systematically analyzed 9000 human and 8000 mouse ubiquitously-accessible candidate cis-regulatory elements (cCREs) with promoter-like signatures (PLSs) from ENCODE, which we denote ubi-PLSs. These are more CpG-rich than non-ubi-PLSs and correspond to genes with ubiquitously high transcription, including a majority of cell-essential genes. ubi-PLSs are enriched with motifs of ubiquitously-expressed transcription factors and preferentially bound by transcriptional cofactors regulating ubiquitously-expressed genes. They are highly conserved between human and mouse at the synteny level but exhibit frequent turnover of motif sites; accordingly, ubi-PLSs show increased variation at their centers compared with flanking regions among the ∼186 thousand human genomes sequenced by the TOPMed project. Finally, ubi-PLSs are enriched in genes implicated in Mendelian diseases, especially diseases broadly impacting most cell types, such as deficiencies in mitochondrial functions. Thus, a set of roughly 9000 mammalian promoters are actively maintained in an accessible state across cell types by a distinct set of transcription factors and cofactors to ensure the transcriptional programs of cell-essential genes.