Integrative Analyses Followed by Functional Characterization Reveal TMEM180 as a Schizophrenia Risk Gene

Brain transcriptome-wide association study implicates novel risk genes underlying schizophrenia risk

BACKGROUND

To identify risk genes whose expression are regulated by the reported risk variants and to explore the potential regulatory mechanism in schizophrenia (SCZ).

METHODS

We systematically integrated three independent brain expression quantitative traits (eQTLs) (CommonMind, GTEx, and BrainSeq Phase 2, a total of 1039 individuals) and GWAS data (56 418 cases and 78 818 controls), with the use of transcriptome-wide association study (TWAS). Diffusion magnetic resonance imaging was utilized to quantify the integrity of white matter bundles and determine whether polygenic risk of novel genes linked to brain structure was present in patients with first-episode antipsychotic SCZ.

RESULTS

TWAS showed that eight risk genes (CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, PCDHA8, THOC7, and TYW5) reached transcriptome-wide significance (TWS) level. These findings were confirmed by an independent integrative approach (i.e. Sherlock). We further conducted conditional analyses and identified the potential risk genes that driven the TWAS association signal in each locus. Gene expression analysis showed that several TWS genes (including CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, THOC7 and TYW5) were dysregulated in the dorsolateral prefrontal cortex of SCZ cases compared with controls. TWS genes were mainly expressed on the surface of glutamatergic neurons, GABAergic neurons, and microglia. Finally, SCZ cases had a substantially greater TWS genes-based polygenic risk (PRS) compared to controls, and we showed that fractional anisotropy of the cingulum-hippocampus mediates the influence of TWS genes PRS on SCZ.

CONCLUSIONS

Our findings identified novel SCZ risk genes and highlighted the importance of the TWS genes in frontal-limbic dysfunctions in SCZ, indicating possible therapeutic targets.

Mendelian Randomization Study Using Dopaminergic Neuron-Specific eQTL Identifies Novel Risk Genes for Schizophrenia

Multiple integrative studies have been performed to identify the potential target genes of the non-coding schizophrenia (SCZ) risk variants. However, all the integrative studies used expression quantitative trait loci (eQTL) data from bulk tissues. Considering the cell type-specific regulatory effect of many genetic variants, it is important to conduct integrative studies using cell type-specific eQTL data. Here, we conduct a Mendelian randomization (MR) study by integrating genome-wide associations of SCZ (74,776 cases and 101,023 controls) and eQTL data (N = 215) from dopaminergic neurons, which were differentiated from human-induced pluripotent stem cell (iPSC) lines. For eQTL from young post-mitotic dopaminergic neurons (differentiation of iPSC for 30 days, D30), we identified 34 genes whose genetically regulated expression in dopaminergic neurons may have a causal role in SCZ. Among which, ARL3 showed the most significant associations with SCZ. For eQTL from more mature dopaminergic neurons (D52), we identified 37 potential SCZ causal genes, and ARL3 and GNL3 showed the most significant associations. Only 12 genes showed significant associations with SCZ in both D30 and D52 eQTL datasets, indicating the time point-specific genetic regulatory effects in young post-mitotic dopaminergic neurons and more mature dopaminergic neurons. Comparing the results from dopaminergic neurons with bulk brain tissues prioritized 2 high-confidence risk genes, including DDHD2 and GALNT10. Our study identifies multiple risk genes whose genetically regulated expression in dopaminergic neurons may have a causal role in SCZ. Further mechanistic investigation will provide pivotal insights into SCZ pathophysiology.

Identifying novel proteins underlying schizophrenia via integrating pQTLs of the plasma, CSF, and brain with GWAS summary data

BACKGROUND

Schizophrenia (SCZ) is a chronic and severe mental illness with no cure so far. Mendelian randomization (MR) is a genetic method widely used to explore etiologies of complex traits. In the current study, we aimed to identify novel proteins underlying SCZ with a systematic analytical approach.

METHODS

We integrated protein quantitative trait loci (pQTLs) of the brain, cerebrospinal fluid (CSF), and plasma with the latest and largest SCZ genome-wide association study (GWAS) via a systematic analytical framework, including two-sample MR analysis, Steiger filtering analysis, and Bayesian colocalization analysis.

RESULTS

The genetically determined protein level of C4A/C4B (OR = 0.70, p = 1.66E-07) in the brain and ACP5 (OR = 0.42, p = 3.73E-05), CNTN2 (OR = 0.62, p = 2.57E-04), and PLA2G7 (OR = 0.71, p = 1.48E-04) in the CSF was associated with a lower risk of SCZ, while the genetically determined protein level of TIE1 (OR = 3.46, p = 4.76E-05), BCL6 (OR = 3.63, p = 1.59E-07), and MICB (OR = 4.49, p = 2.31E-11) in the CSF were associated with an increased risk for SCZ. Pathway enrichment analysis indicated that genetically determined proteins suggestively associated with SCZ were enriched in the biological process of the immune response.

CONCLUSION

In conclusion, we identified one protein in the brain and six proteins in the CSF that showed supporting evidence of being potentially associated with SCZ, which could provide insights into future mechanistic studies to find new treatments for the disease. Our results also supported the important role of neuroinflammation in the pathogenesis of SCZ.

Functional variant rs2270363 on 16p13.3 confers schizophrenia risk by regulating NMRAL1

Recent genome-wide association studies have reported multiple schizophrenia risk loci, yet the functional variants and their roles in schizophrenia remain to be characterized. Here we identify a functional single nucleotide polymorphism (rs2270363: G>A) at the schizophrenia risk locus 16p13.3. rs2270363 lies in the E-box element of the promoter of NMRAL1 and disrupts binding of the basic helix-loop-helix leucine zipper family proteins, including USF1, MAX and MXI1. We validated the regulatory effects of rs2270363 using reporter gene assays and electrophoretic mobility shift assay. Besides, expression quantitative trait loci analysis showed that the risk allele (A) of rs2270363 was significantly associated with elevated NMRAL1 expression in the human brain. Transcription factors knockdown and CRISPR-Cas9-mediated editing further confirmed the regulatory effects of the genomic region containing rs2270363 on NMRAL1. Intriguingly, NMRAL1 was significantly downregulated in the brain of schizophrenia patients compared with healthy subjects, and knockdown of Nmral1 expression affected proliferation and differentiation of mouse neural stem cells, as well as genes and pathways associated with brain development and synaptic transmission. Of note, Nmral1 knockdown resulted in significant decrease of dendritic spine density, revealing the potential pathophysiological mechanisms of NMRAL1 in schizophrenia. Finally, we independently confirmed the association between rs2270363 and schizophrenia in the Chinese population and found that the risk allele of rs2270363 was the same in European and Chinese populations. These lines of evidence suggest that rs2270363 may confer schizophrenia risk by regulating NMRAL1, a gene whose expression dysregulation might be involved in the pathogenesis of schizophrenia by affecting neurodevelopment and synaptic plasticity.

Joint-Tissue Integrative Analysis Identified Hundreds of Schizophrenia Risk Genes

Genome-wide association studies (GWAS) have identified a large number of schizophrenia risk variants, and most of them are mapped to noncoding regions. By leveraging multiple joint-tissue gene expression data and GWAS data, we herein performed a transcriptome-wide association study (TWAS) and Mendelian randomization (MR) analysis and identified 144 genes whose mRNA levels were related to genetic risk of schizophrenia. Most of these genes exhibited diametrically opposite trends of expression in prenatal and postnatal brain tissues, despite that their expression levels in dorsolateral prefrontal cortex (DLPFC) tissues did not significantly differ between schizophrenics and healthy controls. We then found significant enrichment of these genes in dopamine-related pathways that were repeatedly implicated in schizophrenia pathogenesis and in the action of antipsychotic drugs. Gene expression analysis using single cell RNA-sequencing (scRNA-seq) data of mid-gestation fetal brains further revealed enrichment of these genes in glutamatergic excitatory neurons and cycling progenitors. These lines of evidence, in consistency with previous findings, confirmed the polygenic nature of schizophrenia and highlighted involvement of early neurodevelopment aberrations in this disorder. Further investigations using advanced algorithms in both bulk brain tissues and in single cells and at different developmental stages are necessary to characterize transcriptomic features of schizophrenia pathogenesis along brain development.