Neurotransmission and bipolar disorder: a systematic family-based association study

Neurotransmission pathways/systems have been proposed to be involved in the pathophysiology and treatment of bipolar disorder for over 40 years. In order to test the hypothesis that common variants of genes in one or more of five neurotransmission systems confer risk for bipolar disorder, we analyzed 1,005 tag single nucleotide polymorphisms in 90 genes from dopaminergic, serotonergic, noradrenergic, GABAergic, and glutamatergic neurotransmitter systems in 101 trios and 203 quads from Caucasian bipolar families. Our sample has 80% power to detect ORs >or= 1.82 and >or=1.57 for minor allele frequencies of 0.1 and 0.5, respectively. Nominally significant allelic and haplotypic associations were found for genes from each neurotransmission system, with several reaching gene-wide significance (allelic: GRIA1, GRIN2D, and QDPR; haplotypic: GRIN2C, QDPR, and SLC6A3). However, none of these associations survived correction for multiple testing in an individual system, or in all systems considered together. Significant single nucleotide polymorphism associations were not found with sub-phenotypes (alcoholism, psychosis, substance abuse, and suicide attempts) or significant gene-gene interactions. These results suggest that, within the detectable odds ratios of this study, common variants of the selected genes in the five neurotransmission systems do not play major roles in influencing the risk for bipolar disorder or comorbid sub-phenotypes.

Distinct Phenotypes of Inflammation Associated Macrophages and Microglia in the Prefrontal Cortex Schizophrenia Compared to Controls

Approximately 40% of people with schizophrenia are classified as having "high inflammation." This subgroup has worse neuropathology than patients with "low inflammation." Thus, one would expect the resident microglia and possibly monocyte-derived macrophages infiltrating from the periphery to be "activated" in those with schizophrenia with elevated neuroinflammation. To test whether microglia and/or macrophages are associated with increased inflammatory signaling in schizophrenia, we measured microglia- and macrophage-associated transcripts in the postmortem dorsolateral prefrontal cortex of 69 controls and 72 people with schizophrenia. Both groups were stratified by neuroinflammatory status based on cortical mRNA levels of cytokines and SERPINA3. We found microglial mRNAs levels were either unchanged (IBA1 and Hexb, p > 0.20) or decreased (CD11c, <62% p < 0.001) in high inflammation schizophrenia compared to controls. Conversely, macrophage CD163 mRNA levels were increased in patients, substantially so in the high inflammation schizophrenia subgroup compared to low inflammation subgroup (>250%, p < 0.0001). In contrast, high inflammation controls did not have elevated CD163 mRNA compared to low inflammation controls (p > 0.05). The pro-inflammatory macrophage marker (CD64 mRNA) was elevated (>160%, all p < 0.05) and more related to CD163 mRNA in the high inflammation schizophrenia subgroup compared to high inflammation controls, while anti-inflammatory macrophage and cytokine markers (CD206 and IL-10 mRNAs) were either unchanged or decreased in schizophrenia. Finally, macrophage recruitment chemokine CCL2 mRNA was increased in schizophrenia (>200%, p < 0.0001) and CCL2 mRNA levels positively correlated with CD163 mRNA (r = 0.46, p < 0.0001). Collectively, our findings support the co-existence of quiescent microglia and increased pro-inflammatory macrophages in the cortex of people with schizophrenia.

Neuroimmune transcriptome changes in patient brains of psychiatric and neurological disorders

Neuroinflammation has been implicated in multiple brain disorders but the extent and the magnitude of change in immune-related genes (IRGs) across distinct brain disorders has not been directly compared. In this study, 1275 IRGs were curated and their expression changes investigated in 2467 postmortem brains of controls and patients with six major brain disorders, including schizophrenia (SCZ), bipolar disorder (BD), autism spectrum disorder (ASD), major depressive disorder (MDD), Alzheimer's disease (AD), and Parkinson's disease (PD). There were 865 IRGs present across all microarray and RNA-seq datasets. More than 60% of the IRGs had significantly altered expression in at least one of the six disorders. The differentially expressed immune-related genes (dIRGs) shared across disorders were mainly related to innate immunity. Moreover, sex, tissue, and putative cell type were systematically evaluated for immune alterations in different neuropsychiatric disorders. Co-expression networks revealed that transcripts of the neuroimmune systems interacted with neuronal-systems, both of which contribute to the pathology of brain disorders. However, only a few genes with expression changes were also identified as containing risk variants in genome-wide association studies. The transcriptome alterations at gene and network levels may clarify the immune-related pathophysiology and help to better define neuropsychiatric and neurological disorders.

Genetic tests of biologic systems in affective disorders

To liberate candidate gene analyses from criticisms of inexhaustiveness of examination of specific candidate genes, or incompleteness in the choice of candidate genes to study for specific neurobiological pathways, study of sizeable sets of genes pertinent to each putative pathophysiological pathway is required. For many years, genes have been tested in a 'one by one' manner for association with major affective disorders, primarily bipolar illness. However, it is conceivable that not individual genes but abnormalities in several genes within a system or in several neuronal, neural, or hormonal systems are implicated in the functional hypotheses for etiology of affective disorders. Compilation of candidate genes for entire pathways is a challenge, but can reasonably be carried out for the major affective disorders as discussed here. We present here five groupings of genes implicated by neuropharmacological and other evidence, which suggest 252 candidate genes worth examining. Inexhaustiveness of gene interrogation would apply to many studies in which only one polymorphism per gene is analyzed. In contrast to whole-genome association studies, a study of a limited number of candidate genes can readily exploit information on genomic sequence variations obtained from databases and/or resequencing, and has an advantage of not having the complication of an extremely stringent statistical criterion for association.

A comparative study of the genetic components of three subcategories of autism spectrum disorder

The fifth edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-5) controversially combined previously distinct subcategories of autism spectrum disorder (ASD) into a single diagnostic category. However, genetic convergences and divergences between different ASD subcategories are unclear. By retrieving 1725 exonic de novo mutations (DNMs) from 1628 subjects with autistic disorder (AD), 1873 from 1564 subjects with pervasive developmental disorder not otherwise specified (PDD-NOS), 276 from 247 subjects with Asperger's syndrome (AS), and 2077 from 2299 controls, we found that rates of putative functional DNMs (loss-of-function, predicted deleterious missense, and frameshift) in all three subcategories were significantly higher than those in control. We then investigated the convergences and divergences of the three ASD subcategories based on four genetic aspects: whether any two ASD subcategories (1) shared significantly more genes with functional DNMs, (2) exhibited similar spatio-temporal expression patterns, (3) shared significantly more candidate genes, and (4) shared some ASD-associated functional pathways. It is revealed that AD and PDD-NOS were broadly convergent in terms of all four genetic aspects, suggesting these two ASD subcategories may be genetically combined. AS was divergent to AD and PDD-NOS for aspects of functional DNMs and expression patterns, whereas AS and AD/PDD-NOS were convergent for aspects of candidate genes and functional pathways. Our results indicated that the three ASD subcategories present more genetic convergences than divergences, favouring DSM-5's new classification. This study suggests that specifically defined genotypes and their corresponding phenotypes should be integrated analyzed for precise diagnosis of complex disorders, such as ASD.

Positional effects revealed in Illumina methylation array and the impact on analysis

AIM

We aimed to prove the existence of positional effects in the Illumina methylation beadchip data and to find an optimal correction method.

MATERIALS & METHODS

Three HumanMethylation450, three HumanMethylation27 datasets and two EPIC datasets were analyzed. ComBat, linear regression, functional normalization and single-sample Noob were used for minimizing positional effects. The corrected results were evaluated by four methods.

RESULTS

We detected 52,988 CpG loci significantly associated with sample positions, 112 remained after ComBat correction in the primary dataset. The pre- and postcorrection comparisons indicate the positional effects could alter the measured methylation values and downstream analysis results.

CONCLUSION

Positional effects exist in the Illumina methylation array and may bias the analyses. Using ComBat to correct positional effects is recommended.

Brain expression quantitative trait locus mapping informs genetic studies of psychiatric diseases

Genome-wide association study (GWAS) can be used to identify genes that increase the risk of psychiatric diseases. However, much of the disease heritability is still unexplained, suggesting that there are genes to be discovered. Functional annotation of the genetic variants may increase the power of GWAS to identify disease genes, by providing prior information that can be used in Bayesian analysis or in reducing the number of tests. Expression quantitative trait loci (eQTLs) are genomic loci that regulate gene expression. Genetic mapping of eQTLs can help reveal novel functional effects of thousands of single nucleotide polymorphisms (SNPs). The present review mainly focused on the current knowledge on brain eQTL mapping, and discussed some major methodological issues and their possible solutions. The frequently ignored problems of batch effects, covariates, and multiple testing were emphasized, since they can lead to false positives and false negatives. The future application of eQTL data in GWAS analysis was also discussed.

Genome-wide association study of irritable vs. elated mania suggests genetic differences between clinical subtypes of bipolar disorder

The use of clinical features to define subtypes of a disorder may aid in gene identification for complex diseases. In particular, clinical subtypes of mania may distinguish phenotypic subgroups of bipolar subjects that may also differ genetically. To assess this possibility, we performed a genome-wide association study using genotype data from the Bipolar Genome Study (BiGS) and subjects that were categorized as having either irritable or elated mania during their most severe episode. A bipolar case-only analysis in the GAIN bipolar sample identified several genomic regions that differed between irritable and elated subjects, the most significant of which was for 33 SNPs on chromosome 13q31 (peak p = 2×10(-7)). This broad peak is in a relative gene desert over an unknown EST and between the SLITRK1 and SLITRK6 genes. Evidence for association to this region came predominantly from subjects in the sample that were originally collected as part of a family-based bipolar linkage study, rather than those collected as bipolar singletons. We then genotyped an additional sample of bipolar singleton cases and controls, and the analysis of irritable vs. elated mania in this new sample did not replicate our previous findings. However, this lack of replication is likely due to the presence of significant differences in terms of clinical co-morbity that were identified between these singleton bipolar cases and those that were selected from families segregating the disorder. Despite these clinical differences, analysis of the combined sample provided continued support for 13q31 and other regions from our initial analysis. Though genome-wide significance was not achieved, our results suggest that irritable mania results from a distinct set of genes, including a region on chromosome 13q31.

Expression of the G72/G30 gene in transgenic mice induces behavioral changes

The G72/G30 gene complex is a candidate gene for schizophrenia and bipolar disorder. However, G72 and G30 mRNAs are expressed at very low levels in human brain, with only rare splicing forms observed. We report here G72/G30 expression profiles and behavioral changes in a G72/G30 transgenic mouse model. A human BAC clone containing the G72/G30 genomic region was used to establish the transgenic mouse model, on which gene expression studies, western blot and behavioral tests were performed. Relative to their minimal expression in humans, G72 and G30 mRNAs were highly expressed in the transgenic mice, and had a more complex splicing pattern. The highest G72 transcript levels were found in testis, followed by cerebral cortex, with very low or undetectable levels in other tissues. No LG72 (the long putative isoform of G72) protein was detected in the transgenic mice. Whole-genome expression profiling identified 361 genes differentially expressed in transgenic mice compared with wild-type, including genes previously implicated in neurological and psychological disorders. Relative to wild-type mice, the transgenic mice exhibited fewer stereotypic movements in the open field test, higher baseline startle responses in the course of the prepulse inhibition test, and lower hedonic responses in the sucrose preference test. The transcriptome profile changes and multiple mouse behavioral effects suggest that the G72 gene may play a role in modulating behaviors relevant to psychiatric disorders.

Human forebrain organoids reveal connections between valproic acid exposure and autism risk

Valproic acid (VPA) exposure as an environmental factor that confers risk of autism spectrum disorder (ASD), its functional mechanisms in the human brain remain unclear since relevant studies are currently restricted to two-dimensional cell cultures and animal models. To identify mechanisms by which VPA contribute to ASD risk in human, here we used human forebrain organoids (hFOs), in vitro derived three-dimensional cell cultures that recapitulate key human brain developmental features. We identified that VPA exposure in hFOs affected the expression of genes enriched in neural development, synaptic transmission, oxytocin signaling, calcium, and potassium signaling pathways, which have been implicated in ASD. Genes (e.g., CAMK4, CLCN4, DPP10, GABRB3, KCNB1, PRKCB, SCN1A, and SLC24A2) that affected by VPA were significantly overlapped with those dysregulated in brains or organoids derived from ASD patients, and known ASD risk genes, as well as genes in ASD risk-associated gene coexpression modules. Single-cell RNA sequencing analysis showed that VPA exposure affected the expression of genes in choroid plexus, excitatory neuron, immature neuron, and medial ganglionic eminence cells annotated in hFOs. Microelectrode array further identified that VPA exposure in hFOs disrupted synaptic transmission. Taken together, this study connects VPA exposure to ASD pathogenesis using hFOs, which is valuable for illuminating the etiology of ASD and screening for potential therapeutic targets.

An evaluation of the assembly of an approximately 15-Mb region on human chromosome 13q32-q33 linked to bipolar disorder and schizophrenia

The human 13q32-q33 region has been linked to both bipolar disorder and schizophrenia. Before completion of the draft sequences, we developed an approximately 15-Mb comprehensive map for the region extending from D13S1300 to ATA35H12. This map was assembled using publicly available mapping data and sequence-tagged site (STS)-based PCR confirmation. We then compared this map with the NCBI, Celera Genomics, and UCSC Golden Path data in February, June, and September 2001. All data sets showed gaps, misassignment of STSs, and errors in orientation and marker order. Surprisingly, the completed sequences of many bacterial artificial chromosomes (BACs) had been truncated. Of 21 gaps that were detected, 4 were present in both the NCBI and Celera databases. All gaps could be filled using 1-2 BAC clones. A total of 39 loci mapped to additional sites within the human genome, providing evidence of segmental duplications. Additionally, 61 unique cDNA clones were sequenced to increase available transcribed sequence, and 11,353 reference single-nucleotide polymorphisms (SNPs) with an average density of 1 SNP/3720 bases were identified. Overall, integration of the data from multiple sources is still needed for complete assembly of the 13q32-q33 region. (c)

GWAS significance thresholds for deep phenotyping studies can depend upon minor allele frequencies and sample size

An important issue affecting genome-wide association studies with deep phenotyping (multiple correlated phenotypes) is determining the suitable family-wise significance threshold. Straightforward family-wise correction (Bonferroni) of p < 0.05 for 4.3 million genotypes and 335 phenotypes would give a threshold of p < 3.46E-11. This would be too conservative because it assumes all tests are independent. The effective number of tests, both phenotypic and genotypic, must be adjusted for the correlations between them. Spectral decomposition of the phenotype matrix and LD-based correction of the number of tested SNPs are currently used to determine an effective number of tests. In this paper, we compare these calculated estimates with permutation-determined family-wise significance thresholds. Permutations are performed by shuffling individual IDs of the genotype vector for this dataset, to preserve correlation of phenotypes. Our results demonstrate that the permutation threshold is influenced by minor allele frequency (MAF) of the SNPs, and by the number of individuals tested. For the more common SNPs (MAF > 0.1), the permutation family-wise threshold was in close agreement with spectral decomposition methods. However, for less common SNPs (0.05 < MAF ≤ 0.1), the permutation threshold calculated over all SNPs was off by orders of magnitude. This applies to the number of individuals studied (here 777) but not to very much larger numbers. Based on these findings, we propose that the threshold to find a particular level of family-wise significance may need to be established using separate permutations of the actual data for several MAF bins.

Integrative analyses prioritize GNL3 as a risk gene for bipolar disorder

Genome-wide association studies (GWASs) have identified numerous single nucleotide polymorphisms (SNPs) associated with bipolar disorder (BD), but what the causal variants are and how they contribute to BD is largely unknown. In this study, we used FUMA, a GWAS annotation tool, to pinpoint potential causal variants and genes from the latest BD GWAS findings, and performed integrative analyses, including brain expression quantitative trait loci (eQTL), gene coexpression network, differential gene expression, protein-protein interaction, and brain intermediate phenotype association analysis to identify the functions of a prioritized gene and its connection to BD. Convergent lines of evidence prioritized protein-coding gene G Protein Nucleolar 3 (GNL3) as a BD risk gene, with integrative analyses revealing GNL3's roles in cell proliferation, neuronal functions, and brain phenotypes. We experimentally revealed that BD-related eQTL SNPs rs10865973, rs12635140, and rs4687644 regulate GNL3 expression using dual luciferase reporter assay and CRISPR interference experiment in human neural progenitor cells. We further identified that GNL3 knockdown and overexpression led to aberrant neuronal proliferation and differentiation, using two-dimensional human neural cell cultures and three-dimensional forebrain organoid model. This study gathers evidence that BD-related genetic variants regulate GNL3 expression which subsequently affects neuronal proliferation and differentiation.

NRXN1 is associated with enlargement of the temporal horns of the lateral ventricles in psychosis

Schizophrenia, Schizoaffective, and Bipolar disorders share behavioral and phenomenological traits, intermediate phenotypes, and some associated genetic loci with pleiotropic effects. Volumetric abnormalities in brain structures are among the intermediate phenotypes consistently reported associated with these disorders. In order to examine the genetic underpinnings of these structural brain modifications, we performed genome-wide association analyses (GWAS) on 60 quantitative structural brain MRI phenotypes in a sample of 777 subjects (483 cases and 294 controls pooled together). Genotyping was performed with the Illumina PsychChip microarray, followed by imputation to the 1000 genomes multiethnic reference panel. Enlargement of the Temporal Horns of Lateral Ventricles (THLV) is associated with an intronic SNP of the gene NRXN1 (rs12467877, P = 6.76E-10), which accounts for 4.5% of the variance in size. Enlarged THLV is associated with psychosis in this sample, and with reduction of the hippocampus and enlargement of the choroid plexus and caudate. Eight other suggestively significant associations (P < 5.5E-8) were identified with THLV and 5 other brain structures. Although rare deletions of NRXN1 have been previously associated with psychosis, this is the first report of a common SNP variant of NRXN1 associated with enlargement of the THLV in psychosis.

Review of multi-omics data resources and integrative analysis for human brain disorders

In the last decade, massive omics datasets have been generated for human brain research. It is evolving so fast that a timely update is urgently needed. In this review, we summarize the main multi-omics data resources for the human brains of both healthy controls and neuropsychiatric disorders, including schizophrenia, autism, bipolar disorder, Alzheimer's disease, Parkinson's disease, progressive supranuclear palsy, etc. We also review the recent development of single-cell omics in brain research, such as single-nucleus RNA-seq, single-cell ATAC-seq and spatial transcriptomics. We further investigate the integrative multi-omics analysis methods for both tissue and single-cell data. Finally, we discuss the limitations and future directions of the multi-omics study of human brain disorders.

Three dysconnectivity patterns in treatment-resistant schizophrenia patients and their unaffected siblings

Among individuals diagnosed with schizophrenia, approximately 20%–33% are recognized as treatment-resistant schizophrenia (TRS) patients. These TRS patients suffer more severely from the disease but struggle to benefit from existing antipsychotic treatments. A few recent studies suggested that schizophrenia may be caused by impaired synaptic plasticity that manifests as functional dysconnectivity in the brain, however, few of those studies focused on the functional connectivity changes in the brains of TRS groups. In this study, we compared the whole brain connectivity variations in TRS patients, their unaffected siblings, and healthy controls. Connectivity network features between and within the 116 automated anatomical labeling (AAL) brain regions were calculated and compared using maps created with three contrasts: patient vs. control, patient vs. sibling, and sibling vs. control. To evaluate the predictive power of the selected features, we performed a multivariate classification approach. We also evaluated the influence of six important clinical measures (e.g. age, education level) on the connectivity features. This study identified abnormal significant connectivity changes of three patterns in TRS patients and their unaffected siblings: 1) 69 patient-specific connectivity (PCN); 2) 102 shared connectivity (SCN); and 3) 457 unshared connectivity (UCN). While the first two patterns were widely reported by previous non-TRS specific studies, we were among the first to report widespread significant connectivity differences between TRS patient groups and their healthy sibling groups. Observations of this study may provide new insights for the understanding of the neurophysiological mechanisms of TRS.

•

We first compared global functional connectivity in treatment-resistant schizophrenia and their unaffected siblings.

•

Widespread unshared significant functional connectivity in unaffected siblings of treatment-resistant schizophrenia

•

We studied the association of brain connectivity to clinical measures.

Polygenic risk for anxiety influences anxiety comorbidity and suicidal behavior in bipolar disorder

Bipolar disorder is often comorbid with anxiety, which is itself associated with poorer clinical outcomes, including suicide. A better etiologic understanding of this comorbidity could inform diagnosis and treatment. The present study aims to test whether comorbid anxiety in bipolar disorder reflects shared genetic risk factors. We also sought to assess the contribution of genetic risk for anxiety to suicide attempts in bipolar disorder. Polygenic risk scores (PRS) were calculated from published genome-wide association studies of samples of controls and cases with anxiety (n = 83,566) or bipolar disorder (n = 51,710), then scored in independent target samples (total n = 3369) of individuals with bipolar disorder who reported or denied lifetime anxiety disorders or suicidal attempts in research interviews. Participants were recruited from clinical and nonclinical settings and genotyped for common genetic variants. The results show that polygenic risk for anxiety was associated with comorbid anxiety disorders and suicide attempts in bipolar disorder, while polygenic risk for bipolar disorder was not associated with any of these variables. Our findings point out that comorbid anxiety disorders in bipolar disorder reflect a dual burden of bipolar and anxiety-related genes; the latter may also contribute to suicide attempts. Clinical care that recognizes and addresses this dual burden may help improve outcomes in people living with comorbid bipolar and anxiety disorders.

Applying polygenic risk scoring for psychiatric disorders to a large family with bipolar disorder and major depressive disorder

Psychiatric disorders are thought to have a complex genetic pathology consisting of interplay of common and rare variation. Traditionally, pedigrees are used to shed light on the latter only, while here we discuss the application of polygenic risk scores to also highlight patterns of common genetic risk. We analyze polygenic risk scores for psychiatric disorders in a large pedigree (n ~ 260) in which 30% of family members suffer from major depressive disorder or bipolar disorder. Studying patterns of assortative mating and anticipation, it appears increased polygenic risk is contributed by affected individuals who married into the family, resulting in an increasing genetic risk over generations. This may explain the observation of anticipation in mood disorders, whereby onset is earlier and the severity increases over the generations of a family. Joint analyses of rare and common variation may be a powerful way to understand the familial genetics of psychiatric disorders.

Integration of multi-omics data for integrative gene regulatory network inference

Gene regulatory networks provide comprehensive insights and indepth understanding of complex biological processes. The molecular interactions of gene regulatory networks are inferred from a single type of genomic data, e.g., gene expression data in most research. However, gene expression is a product of sequential interactions of multiple biological processes, such as DNA sequence variations, copy number variations, histone modifications, transcription factors, and DNA methylations. The recent rapid advances of high-throughput omics technologies enable one to measure multiple types of omics data, called 'multi-omics data', that represent the various biological processes. In this paper, we propose an Integrative Gene Regulatory Network inference method (iGRN) that incorporates multi-omics data and their interactions in gene regulatory networks. In addition to gene expressions, copy number variations and DNA methylations were considered for multi-omics data in this paper. The intensive experiments were carried out with simulation data, where iGRN's capability that infers the integrative gene regulatory network is assessed. Through the experiments, iGRN shows its better performance on model representation and interpretation than other integrative methods in gene regulatory network inference. iGRN was also applied to a human brain dataset of psychiatric disorders, and the biological network of psychiatric disorders was analysed.

Preliminary genetic imaging study of the association between estrogen receptor-α gene polymorphisms and harsh human maternal parenting

A failure of neural changes initiated by the estrogen surge in late pregnancy to reverse the valence of infant stimuli from aversive to rewarding is associated with dysfunctional maternal behavior in nonhuman mammals. Estrogen receptor-α plays the crucial role in mediating these neural effects of estrogen priming. This preliminary study examines associations between estrogen receptor-α gene polymorphisms and human maternal behavior. Two polymorphisms were associated with human negative maternal parenting. Furthermore, hemodynamic responses in functional magnetic resonance imaging to child stimuli in neural regions associated with social cognition fully mediated the association between genetic variation and negative parenting. This suggests testable hypotheses regarding a biological pathway between genetic variants and dysfunctional human maternal parenting.

Transcription factor POU3F2 regulates TRIM8 expression contributing to cellular functions implicated in schizophrenia

Schizophrenia (SCZ) is a neuropsychiatric disorder with aberrant expression of multiple genes. However, identifying its exact causal genes remains a considerable challenge. The brain-specific transcription factor POU3F2 (POU domain, class 3, transcription factor 2) has been recognized as a risk factor for SCZ, but our understanding of its target genes and pathogenic mechanisms are still limited. Here we report that POU3F2 regulates 42 SCZ-related genes in knockdown and RNA-sequencing experiments of human neural progenitor cells (NPCs). Among those SCZ-related genes, TRIM8 (Tripartite motif containing 8) is located in SCZ-associated genetic locus and is aberrantly expressed in patients with SCZ. Luciferase reporter and electrophoretic mobility shift assays (EMSA) showed that POU3F2 induces TRIM8 expression by binding to the SCZ-associated SNP (single nucleotide polymorphism) rs5011218, which affects POU3F2-binding efficiency at the promoter region of TRIM8. We investigated the cellular functions of POU3F2 and TRIM8 as they co-regulate several pathways related to neural development and synaptic function. Knocking down either POU3F2 or TRIM8 promoted the proliferation of NPCs, inhibited their neuronal differentiation, and impaired the excitatory synaptic transmission of NPC-derived neurons. These results indicate that POU3F2 regulates TRIM8 expression through the SCZ-associated SNP rs5011218, and both genes may be involved in the etiology of SCZ by regulating neural development and synaptic function.

Analyses of the autism-associated neuroligin-3 R451C mutation in human neurons reveal a gain-of-function synaptic mechanism

Mutations in many synaptic genes are associated with autism spectrum disorders (ASD), suggesting that synaptic dysfunction is a key driver of ASD pathogenesis. Among these mutations, the R451C substitution in the NLGN3 gene that encodes the postsynaptic adhesion molecule Neuroligin-3 is noteworthy because it was the first specific mutation linked to ASDs. In mice, the corresponding Nlgn3 R451C-knockin mutation recapitulates social interaction deficits of ASD patients and produces synaptic abnormalities, but the impact of the NLGN3 R451C mutation on human neurons has not been investigated. Here, we generated human knockin neurons with the NLGN3 R451C and NLGN3 null mutations. Strikingly, analyses of NLGN3 R451C-mutant neurons revealed that the R451C mutation decreased NLGN3 protein levels but enhanced the strength of excitatory synapses without affecting inhibitory synapses; meanwhile NLGN3 knockout neurons showed reduction in excitatory synaptic strengths. Moreover, overexpression of NLGN3 R451C recapitulated the synaptic enhancement in human neurons. Notably, the augmentation of excitatory transmission was confirmed in vivo with human neurons transplanted into mouse forebrain. Using single-cell RNA-seq experiments with co-cultured excitatory and inhibitory NLGN3 R451C-mutant neurons, we identified differentially expressed genes in relatively mature human neurons corresponding to synaptic gene expression networks. Moreover, gene ontology and enrichment analyses revealed convergent gene networks associated with ASDs and other mental disorders. Our findings suggest that the NLGN3 R451C mutation induces a gain-of-function enhancement in excitatory synaptic transmission that may contribute to the pathophysiology of ASD.