Jeremie Poschmann

The PsychENCODE project

Recent research on disparate psychiatric disorders has implicated rare variants in genes involved in global gene regulation and chromatin modification, as well as many common variants located primarily in regulatory regions of the genome. Understanding precisely how these variants contribute to disease will require a deeper appreciation for the mechanisms of gene regulation in the developing and adult human brain. The PsychENCODE project aims to produce a public resource of multidimensional genomic data using tissue- and cell type–specific samples from approximately 1,000 phenotypically well-characterized, high-quality healthy and disease-affected human post-mortem brains, as well as functionally characterize disease-associated regulatory elements and variants in model systems. We are beginning with a focus on autism spectrum disorder, bipolar disorder and schizophrenia, and expect that this knowledge will apply to a wide variety of psychiatric disorders. This paper outlines the motivation and design of PsychENCODE.

MiR‐135b is a direct PAX6 target and specifies human neuroectoderm by inhibiting TGF‐β/BMP signaling

Several transcription factors (TFs) have been implicated in neuroectoderm (NE) development, and recently, the TF PAX6 was shown to be critical for human NE specification. However, microRNA networks regulating human NE development have been poorly documented. We hypothesized that microRNAs activated by PAX6 should promote NE development. Using a genomics approach, we identified PAX6 binding sites and active enhancers genome‐wide in an in vitro model of human NE development that was based on neural differentiation of human embryonic stem cells (hESC). PAX6 binding to active enhancers was found in the proximity of several microRNAs, including hsa‐miR‐135b. MiR‐135b was activated during NE development, and ectopic expression of miR‐135b in hESC promoted differentiation toward NE. MiR‐135b promotes neural conversion by targeting components of the TGF‐β and BMP signaling pathways, thereby inhibiting differentiation into alternate developmental lineages. Our results demonstrate a novel TF‐miRNA module that is activated during human neuroectoderm development and promotes the irreversible fate specification of human pluripotent cells toward the neural lineage.

Sensitive detection of chromatin-altering polymorphisms reveals autoimmune disease mechanisms

Most disease associations detected by genome-wide association studies (GWAS) lie outside coding genes, but very few have been mapped to causal regulatory variants. Here, we present a method for detecting regulatory quantitative trait loci (QTLs) that does not require genotyping or whole-genome sequencing. The method combines deep, long-read chromatin immunoprecipitation–sequencing (ChIP-seq) with a statistical test that simultaneously scores peak height correlation and allelic imbalance: the genotype-independent signal correlation and imbalance (G-SCI) test. We performed histone acetylation ChIP-seq on 57 human lymphoblastoid cell lines and used the resulting reads to call 500,066 single-nucleotide polymorphisms de novo within regulatory elements. The G-SCI test annotated 8,764 of these as histone acetylation QTLs (haQTLs)—an order of magnitude larger than the set of candidates detected by expression QTL analysis. Lymphoblastoid haQTLs were highly predictive of autoimmune disease mechanisms. Thus, our method facilitates large-scale regulatory variant detection in any moderately sized cohort for which functional profiling data can be generated, thereby simplifying identification of causal variants within GWAS loci.

Pitfalls of haplotype phasing from amplicon-based long-read sequencing

The long-read sequencers from Pacific Bioscience (PacBio) and Oxford Nanopore Technologies (ONT) offer the opportunity to phase mutations multiple kilobases apart directly from sequencing reads. In this study, we used long-range PCR with ONT and PacBio sequencing to phase two variants 9 kb apart in the RET gene. We also re-analysed data from a recent paper which had apparently successfully used ONT to phase clinically important haplotypes at the CYP2D6 and HLA loci. From these analyses, we demonstrate PCR-chimera formation during PCR amplification and reference alignment bias are pitfalls that need to be considered when attempting to phase variants using amplicon-based long-read sequencing technologies. These methodological pitfalls need to be avoided if the opportunities provided by long-read sequencers are to be fully exploited.

ACSL1 Is Associated With Fetal Programming of Insulin Sensitivity and Cellular Lipid Content

Individuals who are born small for gestational age (SGA) have a risk to develop various metabolic diseases during their life course. The biological memory of the prenatal state of growth restricted individuals may be reflected in epigenetic alterations in stem cell populations. Mesenchymal stem cells (MSCs) from the Whartons jelly of umbilical cord tissue are multipotent, and we generated primary umbilical cord MSC isolates from SGA and normal neonates, which were subsequently differentiated into adipocytes. We established chromatin state maps for histone marks H3K27 acetylation and H3K27 trimethylation and tested whether enrichment of these marks was associated with gene expression changes. After validating gene expression levels for 10 significant chromatin immunoprecipitation sequencing candidate genes, we selected acyl-coenzyme A synthetase 1 (ACSL1) for further investigations due to its key roles in lipid metabolism. The ACSL1 gene was found to be highly associated with histone acetylation in adipocytes differentiated from MSCs with SGA background. In SGA-derived adipocytes, the ACSL1 expression level was also found to be associated with increased lipid loading as well as higher insulin sensitivity. ACSL1 depletion led to changes in expression of candidate genes such as proinflammatory chemokines and down-regulated both, the amount of cellular lipids and glucose uptake. Increased ACSL1, as well as modulated downstream candidate gene expression, may reflect the obese state, as detected in mice fed a high-fat diet. In summary, we believe that ACSL1 is a programmable mediator of insulin sensitivity and cellular lipid content and adipocytes differentiated from Whartons jelly MSCs recapitulate important physiological characteristics of SGA individuals.

The transcription factor SOX6 contributes to the developmental origins of obesity by promoting adipogenesis

An association between impaired fetal growth and the postnatal development of obesity has been established. Here, by comparing adipocytes differentiated from mesenchymal stem cells (MSCs) taken from the umbilical cord and derived from normal and growth-restricted neonates, we identified the transcription factor SOX6 as highly expressed only in growth-restricted individuals. We found that SOX6 regulates adipogenesis in vertebrate species by activating adipogenic regulators including PPARγ, C/EBPα and MEST. We further show that SOX6 interacts with β-catenin in adipocytes, suggesting an inhibition of WNT/β-catenin signaling, thereby promoting adipogenesis. The upstream regulatory region of the MEST gene in MSCs from growth-restricted subjects harbors hypomethylated CpGs next to SOX6 binding motifs, and we found that SOX6 binding is impaired by adjacent CpG methylation. In summary, we report that SOX6 is a novel regulator of adipogenesis synergizing with epigenetic mechanisms. Highlighted article: The developmental transcription factor SOX6 is a positive regulator of adipogenesis in vertebrates and is associated with the fetal origins of human obesity.

Comprehensive benchmarking reveals H2BK20 acetylation as a distinctive signature of cell-state-specific enhancers and promoters.

Although over 35 different histone acetylation marks have been described, the overwhelming majority of regulatory genomics studies focus exclusively on H3K27ac and H3K9ac. In order to identify novel epigenomic traits of regulatory elements, we constructed a benchmark set of validated enhancers by performing 140 enhancer assays in human T cells. We tested 40 chromatin signatures on this unbiased enhancer set and identified H2BK20ac, a little-studied histone modification, as the most predictive mark of active enhancers. Notably, we detected a novel class of functionally distinct enhancers enriched in H2BK20ac but lacking H3K27ac, which was present in all examined cell lines and also in embryonic forebrain tissue. H2BK20ac was also unique in highlighting cell-type-specific promoters. In contrast, other acetylation marks were present in all active promoters, regardless of cell-type specificity. In stimulated microglial cells, H2BK20ac was more correlated with cell-state-specific expression changes than H3K27ac, with TGF-beta signaling decoupling the two acetylation marks at a subset of regulatory elements. In summary, our study reveals a previously unknown connection between histone acetylation and cell-type-specific gene regulation and indicates that H2BK20ac profiling can be used to uncover new dimensions of gene regulation.

Genome-wide DNA methylation profiling identifies convergent molecular signatures associated with idiopathic and syndromic forms of autism in postmortem human brain tissue.

Autism spectrum disorder (ASD) encompasses a collection of complex neuropsychiatric disorders characterized by deficits in social functioning, communication and repetitive behavior. Building on recent studies supporting a role for developmentally moderated regulatory genomic variation in the molecular etiology of ASD, we quantified genome-wide patterns of DNA methylation in 233 post-mortem tissues samples isolated from three brain regions (prefrontal cortex, temporal cortex and cerebellum) dissected from 43 ASD patients and 38 non-psychiatric control donors. We identified widespread differences in DNA methylation associated with idiopathic ASD (iASD), with consistent signals in both cortical regions that were distinct to those observed in the cerebellum. Individuals carrying a duplication on chromosome 15q (dup15q), representing a genetically-defined subtype of ASD, were characterized by striking differences in DNA methylation across a discrete domain spanning an imprinted gene cluster within the duplicated region. In addition to the dramatic cis-effects on DNA methylation observed in dup15q carriers, we identified convergent methylomic signatures associated with both iASD and dup15q, reflecting the findings from previous studies of gene expression and H3K27ac. Cortical co-methylation network analysis identified a number of co-methylated modules significantly associated with ASD that are enriched for genomic regions annotated to genes involved in the immune system, synaptic signalling and neuronal regulation. Our study represents the first systematic analysis of DNA methylation associated with ASD across multiple brain regions, providing novel evidence for convergent molecular signatures associated with both idiopathic and syndromic autism.

A histone acetylome-wide association study of Alzheimer’s disease identifies disease-associated H3K27ac differences in the entorhinal cortex

We quantified genome-wide patterns of lysine H3K27 acetylation (H3K27ac) in entorhinal cortex samples from Alzheimer’s disease (AD) cases and matched controls using chromatin immunoprecipitation and highly parallel sequencing. We observed widespread acetylomic variation associated with AD neuropathology, identifying 4,162 differential peaks (false discovery rate < 0.05) between AD cases and controls. Differentially acetylated peaks were enriched in disease-related biological pathways and included regions annotated to genes involved in the progression of amyloid-β and tau pathology (for example, APP, PSEN1, PSEN2, and MAPT), as well as regions containing variants associated with sporadic late-onset AD. Partitioned heritability analysis highlighted a highly significant enrichment of AD risk variants in entorhinal cortex H3K27ac peak regions. AD-associated variable H3K27ac was associated with transcriptional variation at proximal genes including CR1, GPR22, KMO, PIM3, PSEN1, and RGCC. In addition to identifying molecular pathways associated with AD neuropathology, we present a framework for genome-wide studies of histone modifications in complex disease.Widespread differences in H3K27ac, a key histone modification, are associated with Alzheimer’s disease. H3K27ac differences were enriched in genomic regions containing loci involved in the progression of Aβ and tau pathology.