Identifying novel gene dysregulation associated with opioid overdose death: A meta-analysis of differential gene expression in human prefrontal cortex.

Only recently have human postmortem brain studies of differential gene expression (DGE) associated with opioid overdose death (OOD) been published; sample sizes from these studies have been modest (N = 40-153). To increase statistical power to identify OOD-associated genes, we leveraged human prefrontal cortex RNAseq data from four independent OOD studies and conducted a transcriptome-wide DGE meta-analysis (N = 285). Using a unified gene expression data processing and analysis framework across studies, we meta-analyzed 20 098 genes and found 335 significant differentially expressed genes (DEGs) by OOD status (false discovery rate < 0.05). Of these, 66 DEGs were among the list of 303 genes reported as OOD-associated in prior prefrontal cortex molecular studies, including genes/gene families (e.g., OPRK1, NPAS4, DUSP, EGR). The remaining 269 DEGs were not previously reported (e.g., NR4A2, SYT1, HCRTR2, BDNF). There was little evidence of genetic drivers for the observed differences in gene expression between opioid addiction cases and controls. Enrichment analyses for the DEGs across molecular pathway and biological process databases highlight an interconnected set of genes and pathways from orexin and tyrosine kinase receptors through MEK/ERK/MAPK signaling to affect neuronal plasticity.

Convergence of case-specific epigenetic alterations identify a confluence of genetic vulnerabilities tied to opioid overdose.

Opioid use disorder is a highly heterogeneous disease driven by a variety of genetic and environmental risk factors which have yet to be fully elucidated. Opioid overdose, the most severe outcome of opioid use disorder, remains the leading cause of accidental death in the United States. We interrogated the effects of opioid overdose on the brain using ChIP-seq to quantify patterns of H3K27 acetylation in dorsolateral prefrontal cortical neurons isolated from 51 opioid-overdose cases and 51 accidental death controls. Among opioid cases, we observed global hypoacetylation and identified 388 putative enhancers consistently depleted for H3K27ac. Machine learning on H3K27ac patterns predicted case-control status with high accuracy. We focused on case-specific regulatory alterations, revealing 81,399 hypoacetylation events, uncovering vast inter-patient heterogeneity. We developed a strategy to decode this heterogeneity based on convergence analysis, which leveraged promoter-capture Hi-C to identify five genes over-burdened by alterations in their regulatory network or "plexus": ASTN2, KCNMA1, DUSP4, GABBR2, ENOX1. These convergent loci are enriched for opioid use disorder risk genes and heritability for generalized anxiety, number of sexual partners, and years of education. Overall, our multi-pronged approach uncovers neurobiological aspects of opioid use disorder and captures genetic and environmental factors perpetuating the opioid epidemic.

Cell Type-Specific Intralocus Interactions Reveal Oligodendrocyte Mechanisms in MS.

Multiple sclerosis (MS) is an autoimmune disease characterized by attack on oligodendrocytes within the central nervous system (CNS). Despite widespread use of immunomodulatory therapies, patients may still face progressive disability because of failure of myelin regeneration and loss of neurons, suggesting additional cellular pathologies. Here, we describe a general approach for identifying specific cell types in which a disease allele exerts a pathogenic effect. Applying this approach to MS risk loci, we pinpoint likely pathogenic cell types for 70%. In addition to T cell loci, we unexpectedly identified myeloid- and CNS-specific risk loci, including two sites that dysregulate transcriptional pause release in oligodendrocytes. Functional studies demonstrated inhibition of transcriptional elongation is a dominant pathway blocking oligodendrocyte maturation. Furthermore, pause release factors are frequently dysregulated in MS brain tissue. These data implicate cell-intrinsic aberrations outside of the immune system and suggest new avenues for therapeutic development. VIDEO ABSTRACT.

Rapid functional genetics of the oligodendrocyte lineage using pluripotent stem cells.

Oligodendrocyte dysfunction underlies many neurological disorders, but rapid assessment of mutation-specific effects in these cells has been impractical. To enable functional genetics in oligodendrocytes, here we report a highly efficient method for generating oligodendrocytes and their progenitors from mouse embryonic and induced pluripotent stem cells, independent of mouse strain or mutational status. We demonstrate that this approach, when combined with genome engineering, provides a powerful platform for the expeditious study of genotype-phenotype relationships in oligodendrocytes.

Publisher Correction: ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters.

This corrects the article DOI: 10.1038/ncomms7186.

Hotspots of aberrant enhancer activity punctuate the colorectal cancer epigenome.

In addition to mutations in genes, aberrant enhancer element activity at non-coding regions of the genome is a key driver of tumorigenesis. Here, we perform epigenomic enhancer profiling of a cohort of more than forty genetically diverse human colorectal cancer (CRC) specimens. Using normal colonic crypt epithelium as a comparator, we identify enhancers with recurrently gained or lost activity across CRC specimens. Of the enhancers highly recurrently activated in CRC, most are constituents of super enhancers, are occupied by AP-1 and cohesin complex members, and originate from primed chromatin. Many activate known oncogenes, and CRC growth can be mitigated through pharmacologic inhibition or genome editing of these loci. Nearly half of all GWAS CRC risk loci co-localize to recurrently activated enhancers. These findings indicate that the CRC epigenome is defined by highly recurrent epigenetic alterations at enhancers which activate a common, aberrant transcriptional programme critical for CRC growth and survival.

Modeling disease risk through analysis of physical interactions between genetic variants within chromatin regulatory circuitry.

SNPs associated with disease susceptibility often reside in enhancer clusters, or super-enhancers. Constituents of these enhancer clusters cooperate to regulate target genes and often extend beyond the linkage disequilibrium (LD) blocks containing risk SNPs identified in genome-wide association studies (GWAS). We identified 'outside variants', defined as SNPs in weak LD with GWAS risk SNPs that physically interact with risk SNPs as part of a target gene's regulatory circuitry. These outside variants further explain variation in target gene expression beyond that explained by GWAS-associated SNPs. Additionally, the clinical risk associated with GWAS SNPs is considerably modified by the genotype of outside variants. Collectively, these findings suggest a potential model in which outside variants and GWAS SNPs that physically interact in 3D chromatin collude to influence target transcript levels as well as clinical risk. This model offers an additional hypothesis for the source of missing heritability for complex traits.

ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters.

Chromatin interactions connect distal regulatory elements to target gene promoters guiding stimulus- and lineage-specific transcription. Few factors securing chromatin interactions have so far been identified. Here, by integrating chromatin interaction maps with the large collection of transcription factor-binding profiles provided by the ENCODE project, we demonstrate that the zinc-finger protein ZNF143 preferentially occupies anchors of chromatin interactions connecting promoters with distal regulatory elements. It binds directly to promoters and associates with lineage-specific chromatin interactions and gene expression. Silencing ZNF143 or modulating its DNA-binding affinity using single-nucleotide polymorphisms (SNPs) as a surrogate of site-directed mutagenesis reveals the sequence dependency of chromatin interactions at gene promoters. We also find that chromatin interactions alone do not regulate gene expression. Together, our results identify ZNF143 as a novel chromatin-looping factor that contributes to the architectural foundation of the genome by providing sequence specificity at promoters connected with distal regulatory elements.

Long ncRNA expression associates with tissue-specific enhancers.

Long non-coding RNAs (ncRNA) have recently been demonstrated to be expressed from a subset of enhancers and to be required for the distant regulation of gene expression. Several approaches to predict enhancers have been developed based on various chromatin marks and occupancy of enhancer-binding proteins. Despite the rapid advances in the field, no consensus how to define tissue specific enhancers yet exists. Here, we identify 2,695 long ncRNAs annotated by ENCODE (corresponding to 28% of all ENCODE annotated long ncRNAs) that overlap tissue-specific enhancers. We use a recently developed algorithm to predict tissue-specific enhancers, PreSTIGE, that is based on the H3K4me1 mark and tissue specific expression of mRNAs. The expression of the long ncRNAs overlapping enhancers is significantly higher when the enhancer is predicted as active in a specific cell line, suggesting a general interdependency of active enhancers and expression of long ncRNAs. This dependency is not identified using previous enhancer prediction algorithms that do not account for expression of their downstream targets. The predicted enhancers that overlap annotated long ncRNAs generally have a lower ratio of H3K4me1 to H3K4me3, suggesting that enhancers expressing long ncRNAs might be associated with specific epigenetic marks. In conclusion, we demonstrate the tissue-specific predictive power of PreSTIGE and provide evidence for thousands of long ncRNAs that are expressed from active tissue-specific enhancers, suggesting a particularly important functional relationship between long ncRNAs and enhancer activity in determining tissue-specific gene expression.

Enhancer variants: evaluating functions in common disease.

Gene enhancer elements are noncoding segments of DNA that play a central role in regulating transcriptional programs that control development, cell identity, and evolutionary processes. Recent studies have shown that noncoding single nucleotide polymorphisms (SNPs) that have been associated with risk for numerous common diseases through genome-wide association studies frequently lie in cell-type-specific enhancer elements. These enhancer variants probably influence transcriptional output, thereby offering a mechanistic basis to explain their association with risk for many common diseases. This review focuses on the identification and interpretation of disease-susceptibility variants that influence enhancer function. We discuss strategies for prioritizing the study of functional enhancer SNPs over those likely to be benign, review experimental and computational approaches to identifying the gene targets of enhancer variants, and highlight efforts to quantify the impact of enhancer variants on target transcript levels and cellular phenotypes. These studies are beginning to provide insights into the mechanistic basis of many common diseases, as well as into how we might translate this knowledge for improved disease diagnosis, prevention and treatments. Finally, we highlight five major challenges often associated with interpreting enhancer variants, and discuss recent technical advances that may help to surmount these challenges.

Epigenomic comparison reveals activation of "seed" enhancers during transition from naive to primed pluripotency.

Naive mouse embryonic stem cells (mESCs) and primed epiblast stem cells (mEpiSCs) represent successive snapshots of pluripotency during embryogenesis. Using transcriptomic and epigenomic mapping we show that a small fraction of transcripts are differentially expressed between mESCs and mEpiSCs and that these genes show expected changes in chromatin at their promoters and enhancers. Unexpectedly, the cis-regulatory circuitry of genes that are expressed at identical levels between these cell states also differs dramatically. In mESCs, these genes are associated with dominant proximal enhancers and dormant distal enhancers, which we term seed enhancers. In mEpiSCs, the naive-dominant enhancers are lost, and the seed enhancers take up primary transcriptional control. Seed enhancers have increased sequence conservation and show preferential usage in downstream somatic tissues, often expanding into super enhancers. We propose that seed enhancers ensure proper enhancer utilization and transcriptional fidelity as mammalian cells transition from naive pluripotency to a somatic regulatory program.

Combinatorial effects of multiple enhancer variants in linkage disequilibrium dictate levels of gene expression to confer susceptibility to common traits.

DNA variants (SNPs) that predispose to common traits often localize within noncoding regulatory elements such as enhancers. Moreover, loci identified by genome-wide association studies (GWAS) often contain multiple SNPs in linkage disequilibrium (LD), any of which may be causal. Thus, determining the effect of these multiple variant SNPs on target transcript levels has been a major challenge. Here, we provide evidence that for six common autoimmune disorders (rheumatoid arthritis, Crohn's disease, celiac disease, multiple sclerosis, lupus, and ulcerative colitis), the GWAS association arises from multiple polymorphisms in LD that map to clusters of enhancer elements active in the same cell type. This finding suggests a "multiple enhancer variant" hypothesis for common traits, where several variants in LD impact multiple enhancers and cooperatively affect gene expression. Using a novel method to delineate enhancer-gene interactions, we show that multiple enhancer variants within a given locus typically target the same gene. Using available data from HapMap and B lymphoblasts as a model system, we provide evidence at numerous loci that multiple enhancer variants cooperatively contribute to altered expression of their gene targets. The effects on target transcript levels tend to be modest and can be either gain- or loss-of-function. Additionally, the genes associated with multiple enhancer variants encode proteins that are often functionally related and enriched in common pathways. Overall, the multiple enhancer variant hypothesis offers a new paradigm by which noncoding variants can confer susceptibility to common traits.

Epigenomic enhancer profiling defines a signature of colon cancer.

Cancer is characterized by gene expression aberrations. Studies have largely focused on coding sequences and promoters, even though distal regulatory elements play a central role in controlling transcription patterns. We used the histone mark H3K4me1 to analyze gain and loss of enhancer activity genome-wide in primary colon cancer lines relative to normal colon crypts. We identified thousands of variant enhancer loci (VELs) that comprise a signature that is robustly predictive of the in vivo colon cancer transcriptome. Furthermore, VELs are enriched in haplotype blocks containing colon cancer genetic risk variants, implicating these genomic regions in colon cancer pathogenesis. We propose that reproducible changes in the epigenome at enhancer elements drive a specific transcriptional program to promote colon carcinogenesis.