Functional Genomics Identify a Regulatory Risk Variation rs4420550 in the 16p11.2 Schizophrenia-Associated Locus

Hong Chang,Xin Cai,Hui-Juan Li,Wei-Peng Liu,Li-Juan Zhao,Chu-Yi Zhang,Jun-Yang Wang,Jie-Wei Liu,Xiao-Lei Ma,Lu Wang,Yong-Gang Yao,Xiong-Jian Luo,Ming Li,Xiao Xiao
DOI: https://doi.org/10.1016/j.biopsych.2020.09.016
IF: 12.81
2021-02-01
Biological Psychiatry
Abstract:<h3 class="u-h4 u-margin-m-top u-margin-xs-bottom">Background</h3><p>Genome-wide association studies (GWAS) have reported hundreds of genomic loci associated with schizophrenia, yet identifying the functional risk variations is a key step in elucidating the underlying mechanisms.</p><h3 class="u-h4 u-margin-m-top u-margin-xs-bottom">Methods</h3><p>We applied multiple bioinformatics and molecular approaches, including expression quantitative trait loci (eQTL) analyses, epigenome signature identification, luciferase reporter assay, chromatin conformation capture (3C), homology-directed genome editing by CRISPR/Cas9, RNA-sequencing (RNA-Seq) and assay for transposase-accessible chromatin using sequencing (ATAC-Seq).</p><h3 class="u-h4 u-margin-m-top u-margin-xs-bottom">Results</h3><p>We found that the schizophrenia GWAS risk variations at 16p11.2 were significantly associated with mRNA levels of multiple genes in human brain, and one of the leading eQTL genes, <em>MAPK3</em>, located ∼200-kb away from these risk variations in the genome. Further analyses based on the epigenome marks in human brain and cell lines suggested that a noncoding single nucleotide polymorphism (SNP) rs4420550 (p=2.36×10<sup>–9</sup> in schizophrenia GWAS) was within a DNA enhancer region, which was validated via <em>in vitro</em> luciferase reporter assays. The 3C experiment showed that the rs4420550 region physically interacted with the <em>MAPK3</em> promoter and <em>TAOK2</em> promoter. Precise CRISPR/Cas9 editing of single base pair in cells followed by RNA-Seq further confirmed the regulatory effects of rs4420550 on the transcription of 16p11.2 genes, and ATAC-Seq demonstrated that rs4420550 affected chromatin accessibility at the 16p11.2 region. The rs4420550-[A/A] cells showed significantly higher proliferation rates compared with rs4420550-[G/G] cells.</p><h3 class="u-h4 u-margin-m-top u-margin-xs-bottom">Conclusions</h3><p>These results together suggest that rs4420550 is a functional risk variation, and this study illustrates an example of comprehensive functional characterization of schizophrenia GWAS risk loci.</p>
neurosciences,psychiatry
What problem does this paper attempt to address?