Transcriptome outlier analysis implicates schizophrenia susceptibility genes and enriches putatively functional rare genetic variants.

We searched a gene expression dataset comprised of 634 schizophrenia (SZ) cases and 713 controls for expression outliers (i.e., extreme tails of the distribution of transcript expression values) with SZ cases overrepresented compared with controls. These outlier genes were enriched for brain expression and for genes known to be associated with neurodevelopmental disorders. SZ cases showed higher outlier burden (i.e., total outlier events per subject) than controls for genes within copy number variants (CNVs) associated with SZ or neurodevelopmental disorders. Outlier genes were enriched for CNVs and for rare putative regulatory variants, but this only explained a small proportion of the outlier subjects, highlighting the underlying presence of additional genetic and potentially, epigenetic mechanisms.

A rare functional noncoding variant at the GWAS-implicated MIR137/MIR2682 locus might confer risk to schizophrenia and bipolar disorder.

Schizophrenia (SZ) genome-wide association studies (GWASs) have identified common risk variants in >100 susceptibility loci; however, the contribution of rare variants at these loci remains largely unexplored. One of the strongly associated loci spans MIR137 (miR137) and MIR2682 (miR2682), two microRNA genes important for neuronal function. We sequenced ∼6.9 kb MIR137/MIR2682 and upstream regulatory sequences in 2,610 SZ cases and 2,611 controls of European ancestry. We identified 133 rare variants with minor allele frequency (MAF) <0.5%. The rare variant burden in promoters and enhancers, but not insulators, was associated with SZ (p = 0.021 for MAF < 0.5%, p = 0.003 for MAF < 0.1%). A rare enhancer SNP, 1:g.98515539A>T, presented exclusively in 11 SZ cases (nominal p = 4.8 × 10(-4)). We further identified its risk allele T in 2 of 2,434 additional SZ cases, 11 of 4,339 bipolar (BP) cases, and 3 of 3,572 SZ/BP study controls and 1,688 population controls; yielding combined p values of 0.0007, 0.0013, and 0.0001 for SZ, BP, and SZ/BP, respectively. The risk allele T of 1:g.98515539A>T reduced enhancer activity of its flanking sequence by >50% in human neuroblastoma cells, predicting lower expression of MIR137/MIR2682. Both empirical and computational analyses showed weaker transcription factor (YY1) binding by the risk allele. Chromatin conformation capture (3C) assay further indicated that 1:g.98515539A>T influenced MIR137/MIR2682, but not the nearby DPYD or LOC729987. Our results suggest that rare noncoding risk variants are associated with SZ and BP at MIR137/MIR2682 locus, with risk alleles decreasing MIR137/MIR2682 expression.

MicroRNA-9 and microRNA-326 regulate human dopamine D2 receptor expression, and the microRNA-mediated expression regulation is altered by a genetic variant.

The human dopamine receptor D2 (DRD2) has been implicated in the pathophysiology of schizophrenia and other neuropsychiatric disorders. Most antipsychotic drugs influence dopaminergic transmission through blocking dopamine receptors, primarily DRD2. We report here the post-transcriptional regulation of DRD2 expression by two brain-expressed microRNAs (miRs), miR-326 and miR-9, in an ex vivo mode, and show the relevance of miR-mediated DRD2 expression regulation in human dopaminergic neurons and in developing human brains. Both miRs targeted the 3'-UTR (untranslated region) of DRD2 in NT2 (neuron-committed teratocarcinoma, which endogenously expresses DRD2) and CHO (Chinese hamster ovary) cell lines, decreasing luciferase activity measured by a luciferase reporter gene assay. miR-326 overexpression reduced DRD2 mRNA and DRD2 receptor synthesis. Both antisense miR-326 and antisense miR-9 increased DRD2 protein abundance, suggesting an endogenous repression of DRD2 expression by both miRs. Furthermore, a genetic variant (rs1130354) within the DRD2 3'-UTR miR-targeting site interferes with miR-326-mediated repression of DRD2 expression. Finally, co-expression analysis identified an inverse correlation of DRD2 expression with both miR-326 and miR-9 in differentiating dopaminergic neurons derived from human induced pluripotent stem cells (iPSCs) and in developing human brain regions implicated in schizophrenia. Our study provides empirical evidence suggesting that miR-326 and miR-9 may regulate dopaminergic signaling, and miR-326 and miR-9 may be considered as potential drug targets for the treatment of disorders involving abnormal DRD2 function, such as schizophrenia.

Open chromatin dynamics reveals stage-specific transcriptional networks in hiPSC-based neurodevelopmental model.

Chromatin accessibility to transcription factors (TFs) strongly influences gene transcription and cell differentiation. However, a mechanistic understanding of the transcriptional control during the neuronal differentiation of human induced pluripotent stem cells (hiPSCs), a promising cellular model for mental disorders, remains elusive. Here, we carried out additional analyses on our recently published open chromatin regions (OCRs) profiling at different stages of hiPSC neuronal differentiation. We found that the dynamic changes of OCR during neuronal differentiation highlighted cell stage-specific gene networks, and the chromatin accessibility at the core promoter region of a gene correlates with the corresponding transcript abundance. Within the cell stage-specific OCRs, we identified the binding of cell stage-specific TFs and observed a lag of a neuronal TF binding behind the mRNA expression of the corresponding TF. Interestingly, binding footprints of NEUROD1 and NEUROG2, both of which induce high efficient conversion of hiPSCs to glutamatergic neurons, were among those most enriched in the relatively mature neurons. Furthermore, TF network analysis showed that both NEUROD1 and NEUROG2 were present in the same core TF network specific to more mature neurons, suggesting a pivotal mechanism of epigenetic control of neuronal differentiation and maturation. Our study provides novel insights into the epigenetic control of glutamatergic neurogenesis in the context of TF networks, which may be instrumental to improving hiPSC modeling of neuropsychiatric disorders.

Open Chromatin Profiling in hiPSC-Derived Neurons Prioritizes Functional Noncoding Psychiatric Risk Variants and Highlights Neurodevelopmental Loci.

Most disease variants lie within noncoding genomic regions, making their functional interpretation challenging. Because chromatin openness strongly influences transcriptional activity, we hypothesized that cell-type-specific open chromatin regions (OCRs) might highlight disease-relevant noncoding sequences. To investigate, we mapped global OCRs in neurons differentiating from hiPSCs, a cellular model for studying neurodevelopmental disorders such as schizophrenia (SZ). We found that the OCRs are highly dynamic and can stratify GWAS-implicated SZ risk variants. Of the more than 3,500 SZ-associated variants analyzed, we prioritized ∼100 putatively functional ones located in neuronal OCRs, including rs1198588, at a leading risk locus flanking MIR137. Excitatory neurons derived from hiPSCs with CRISPR/Cas9-edited rs1198588 or a rare proximally located SZ risk variant showed altered MIR137 expression, dendrite arborization, and synapse maturation. Our study shows that noncoding disease variants in OCRs can affect neurodevelopment, and that analysis of open chromatin regions can help prioritize functionally relevant noncoding variants identified by GWAS.