A missense variant in NDUFA6 confers schizophrenia risk by affecting YY1 binding and NAGA expression

Genome-wide association studies (GWASs) have revealed that genetic variants at the 22q13.2 risk locus were robustly associated with schizophrenia. However, the causal variants at this risk locus and their roles in schizophrenia remain elusive. Here we identify the risk missense variant rs1801311 (located in the 1st exon of NDUFA6 gene) as likely causal for schizophrenia at 22q13.2 by disrupting binding of YY1, TAF1, and POLR2A. We systematically elucidated the regulatory mechanisms of rs1801311 and validated the regulatory effect of this missense variant. Intriguingly, rs1801311 physically interacted with NAGA (encodes the alpha-N-acetylgalactosaminidase, which is mainly involved in regulating metabolisms of glycoproteins and glycolipids in lysosome) and showed the most significant association with NAGA expression in the human brain, with the risk allele (G) associated with higher NAGA expression. Consistent with eQTL analysis, expression analysis showed that NAGA was significantly upregulated in brains of schizophrenia cases compared with controls, further supporting that rs1801311 may confer schizophrenia risk by regulating NAGA expression. Of note, we found that NAGA regulates important neurodevelopmental processes, including proliferation and differentiation of neural stem cells. Transcriptome analysis corroborated that NAGA regulates pathways associated with neuronal differentiation. Finally, we independently confirmed the association between rs1801311 and schizophrenia in a large Chinese cohort. Our study elucidates the regulatory mechanisms of the missense schizophrenia risk variant rs1801311 and provides mechanistic links between risk variant and schizophrenia etiology. In addition, this study also revealed the novel role of coding variants in gene regulation and schizophrenia risk, i.e., genetic variant in coding region of a specific gene may confer disease risk through regulating distal genes (act as regulatory variant for distal genes).