Student or educator? Sign up with your school email and get 6 months of Pro free.
CoChat
Sign Up
Cited by 73

The DOI matches its Crossref record; 34 of the 34 references it cites were found in a registry or an index.

Schizophrenia-associated NRXN1 deletions induce developmental-timing- and cell-type-specific vulnerabilities in human brain organoids

Sebastian, Rebecca, Jin, Kang, Pavon, Narciso, Bansal, Ruby, Potter, Andrew, Song, Yoonjae, Babu, Juliana, Gabriel, Rafael, Sun, Yubing, Aronow, Bruce, Pak, ChangHui·2023·Nature Communications
DOI 10.1038/s41467-023-39420-6PMC PMC10290702
Open this paper in CoChat — free account

Abstract

De novo mutations and copy number deletions in NRXN1 (2p16.3) pose a significant risk for schizophrenia (SCZ). It is unclear how NRXN1 deletions impact cortical development in a cell type-specific manner and disease background modulates these phenotypes. Here, we leveraged human pluripotent stem cell-derived forebrain organoid models carrying NRXN1 heterozygous deletions in isogenic and SCZ patient genetic backgrounds and conducted single-cell transcriptomic analysis over the course of brain organoid development from 3 weeks to 3.5 months. Intriguingly, while both deletions similarly impacted molecular pathways associated with ubiquitin-proteasome system, alternative splicing, and synaptic signaling in maturing glutamatergic and GABAergic neurons, SCZ- NRXN1 deletions specifically perturbed developmental trajectories of early neural progenitors and accumulated disease-specific transcriptomic signatures. Using calcium imaging, we found that both deletions led to long-lasting changes in spontaneous and synchronous neuronal networks, implicating synaptic dysfunction. Our study reveals developmental-timing- and cell-type-dependent actions of NRXN1 deletions in unique genetic contexts.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-023-39420-6

Cites

34 of 34 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Generation of forebrain organoids from genetically engineered NRXN1 cKO hESCs and donor derived iPSCs from healthy controls and SCZ patients. a Schematic of brain organoid generation protocol and the corresponding representative brightfield images over development. Scale bars – 250 μm for d-1, 0, 3; 100 μm for d35+. b Schematic showing genotypes used for scRNAseq and other experiments in the study. c Uniform Manifold Approximation and Projection (UMAP) showing distributions of cell classes of the integrated single-cell data. Abbreviations: NECs neural precursor cells, oRG outer radial glial cells, RG-like radial glial-like cells, IPC intermediate precursor cells, CN cortical excitatory neurons, and IN cortical GABAergic inhibitory neurons. d Heatmap of gene signatures across cell classes and cell groups. Student’s t tests were used to compute differentially expressed genes (DEGs) of each cell class in the integrated scRNA-seq data. Top 50 DEGs for each cell type are shown in a diagonal manner while normalized expression levels are scaled for each row. e UMAP showing distribution of pseudotime values of integrated single-cell data inferred by Monocle3. f NRXN1 mRNA expression (Seurat scaled expression values) across pseudotime in control donor-derived brain organoids. A higher pseudotime value indicates greater maturity as indicated by the number of various neuronal cell classes.
Fig. 1 Generation of forebrain organoids from genetically engineered NRXN1 cKO hESCs and donor derived iPSCs from healthy controls and SCZ patients. a Schematic of brain organoid generation protocol and the corresponding representative brightfield images over development. Scale bars – 250 μm for d-1, 0, 3; 100 μm for d35+. b Schematic showing genotypes used for scRNAseq and other experiments in the study. c Uniform Manifold Approximation and Projection (UMAP) showing distributions of cell classes of the integrated single-cell data. Abbreviations: NECs neural precursor cells, oRG outer radial glial cells, RG-like radial glial-like cells, IPC intermediate precursor cells, CN cortical excitatory neurons, and IN cortical GABAergic inhibitory neurons. d Heatmap of gene signatures across cell classes and cell groups. Student’s t tests were used to compute differentially expressed genes (DEGs) of each cell class in the integrated scRNA-seq data. Top 50 DEGs for each cell type are shown in a diagonal manner while normalized expression levels are scaled for each row. e UMAP showing distribution of pseudotime values of integrated single-cell data inferred by Monocle3. f NRXN1 mRNA expression (Seurat scaled expression values) across pseudotime in control donor-derived brain organoids. A higher pseudotime value indicates greater maturity as indicated by the number of various neuronal cell classes.CC BY · PubMed Central
Fig. 2 Perturbation effects of NRXN1 engineered deletions across development. a UMAPs showing distributions of cell classes (left), time points (top right), and genotypes (bottom right) of NRXN1 cKO engineered brain organoids. A total of 10 samples were processed for deletion (Cre) control (Flp) – n = 1 for 3 weeks (3 wk) Cre and Flp; n = 2 for 2 months (2 mo Cre and Flp; n = 2 for 3.5 months (3.5 mo) Cre and Flp. b Ridge plots showing the density of cell abundance across the dimension of pseudotime for 3 wk, 2 mo, and 3.5 mo engineered brain organoids. c, d Dot plots showing significant differential gene expression in deletion vs. control across time points and cell types. The size and color of each dot show the number and significance of overlapping DEGs for each comparison of two cell classes. The significance was measured by −log10(FDRadjustedpvalues) of hypergeometric tests (see Methods). e Representative enriched gene sets were shown for GSEA results of DEGs of multiple cell classes using ToppGene. Different categories of Gene Ontology are shown in different colors (blue: molecular function, orange: biological process, and red: pathway). Enrichment scores were defined as −log10(FDRadjustedpvalues) to represent the associations between DEG sets and Gene Ontology gene sets.
Fig. 2 Perturbation effects of NRXN1 engineered deletions across development. a UMAPs showing distributions of cell classes (left), time points (top right), and genotypes (bottom right) of NRXN1 cKO engineered brain organoids. A total of 10 samples were processed for deletion (Cre) control (Flp) – n = 1 for 3 weeks (3 wk) Cre and Flp; n = 2 for 2 months (2 mo Cre and Flp; n = 2 for 3.5 months (3.5 mo) Cre and Flp. b Ridge plots showing the density of cell abundance across the dimension of pseudotime for 3 wk, 2 mo, and 3.5 mo engineered brain organoids. c, d Dot plots showing significant differential gene expression in deletion vs. control across time points and cell types. The size and color of each dot show the number and significance of overlapping DEGs for each comparison of two cell classes. The significance was measured by −log10(FDRadjustedpvalues) of hypergeometric tests (see Methods). e Representative enriched gene sets were shown for GSEA results of DEGs of multiple cell classes using ToppGene. Different categories of Gene Ontology are shown in different colors (blue: molecular function, orange: biological process, and red: pathway). Enrichment scores were defined as −log10(FDRadjustedpvalues) to represent the associations between DEG sets and Gene Ontology gene sets.CC BY · PubMed Central
Fig. 3 Perturbation effects of SCZ associated NRXN1 deletions across development. a UMAPs showing distributions of cell classes (left), time points (top right), and genotypes (bottom right) of SCZ-NRXN1del donor derived brain organoids. A total of 16 samples were processed – 2 SCZ-NRXN1del donors and 2 Ctrl donors for 3 weeks (3 wk); 2 SCZ-NRXN1del donors and 2 Ctrl donors for 2 months (2 mo); 4 SCZ-NRXN1del donors and 4 Ctrl donors for 3.5 months (3.5 mo). b Ridge plots showing the density of cell abundance across the dimension of pseudotime for 3 wk, 2 mo, and 3.5 mo engineered brain organoids. c, d Dot plots showing significant differential gene expression in deletion vs. control across time points and cell types. The size and color of each dot show the number and significance of overlapping DEGs for each comparison of two cell classes. The significance was measured by −log10(FDRadjustedpvalues) of hypergeometric tests (see Methods). e Representative enriched gene sets were shown for GSEA results of DEGs of multiple cell classes using ToppGene. Different categories of Gene Ontology are shown in different colors (blue: molecular function, orange: biological process, and red: pathway). Enrichment scores were defined as −log10(FDRadjustedpvalues) to represent the associations between DEG sets and Gene Ontology gene sets.
Fig. 3 Perturbation effects of SCZ associated NRXN1 deletions across development. a UMAPs showing distributions of cell classes (left), time points (top right), and genotypes (bottom right) of SCZ-NRXN1del donor derived brain organoids. A total of 16 samples were processed – 2 SCZ-NRXN1del donors and 2 Ctrl donors for 3 weeks (3 wk); 2 SCZ-NRXN1del donors and 2 Ctrl donors for 2 months (2 mo); 4 SCZ-NRXN1del donors and 4 Ctrl donors for 3.5 months (3.5 mo). b Ridge plots showing the density of cell abundance across the dimension of pseudotime for 3 wk, 2 mo, and 3.5 mo engineered brain organoids. c, d Dot plots showing significant differential gene expression in deletion vs. control across time points and cell types. The size and color of each dot show the number and significance of overlapping DEGs for each comparison of two cell classes. The significance was measured by −log10(FDRadjustedpvalues) of hypergeometric tests (see Methods). e Representative enriched gene sets were shown for GSEA results of DEGs of multiple cell classes using ToppGene. Different categories of Gene Ontology are shown in different colors (blue: molecular function, orange: biological process, and red: pathway). Enrichment scores were defined as −log10(FDRadjustedpvalues) to represent the associations between DEG sets and Gene Ontology gene sets.CC BY · PubMed Central
Fig. 4 Dysregulated alternative splicing in NRXN1 mutant glutamatergic human neurons. a Differential gene expression (adjusted p values < 0.05) of neuronal splicing regulators in SCZ-NRXN1del1 compared to Ctrl.1 iN cells (n = 3 replicates). Log2 fold changes are plotted on the y-axis with example DEGs on the x-axis. Blue bars represent downregulation and red bars represent upregulation. b Venn diagram representing differential isoform abundance derived from RSEM analysis across iN cells from engineered NRXN1 deletions (n = 3 replicates from NRXN1 cKO iPSC) and donor-derived SCZ-NRXN1 deletions (3 healthy controls and 3 SCZ-NRXN1 deletion carriers; n = 3 replicates per line, 18 replicates total). Number of differentially expressed isoforms (184 for engineered iNs and 213 donor derived iNs) mapped to 160 and 148 unique gene IDs respectively. Among the unique genes, example SYNGO-mapped genes are highlighted. c, d Percentage of isoform ratio change (differential transcript usage, see methods) for NRXN1 and CHD4 in NRXN1 deletions from engineered (eng.; yellow) and donor (green) backgrounds. Error bars represent +/−S.E.M. e Detection and visualization of local differential splicing using Majiq and ggSashimi. Change in percent spliced in index (dPSI) is calculated for differential exon usage in the NRXN1 gene (exons 19–21) in the SCZ-NRXN1 deletion vs. NRXN1 cKO engineered deletion iNs. Number of detected reads is plotted on the y-axis with location of genomic coordinates on the x-axis.
Fig. 4 Dysregulated alternative splicing in NRXN1 mutant glutamatergic human neurons. a Differential gene expression (adjusted p values < 0.05) of neuronal splicing regulators in SCZ-NRXN1del1 compared to Ctrl.1 iN cells (n = 3 replicates). Log2 fold changes are plotted on the y-axis with example DEGs on the x-axis. Blue bars represent downregulation and red bars represent upregulation. b Venn diagram representing differential isoform abundance derived from RSEM analysis across iN cells from engineered NRXN1 deletions (n = 3 replicates from NRXN1 cKO iPSC) and donor-derived SCZ-NRXN1 deletions (3 healthy controls and 3 SCZ-NRXN1 deletion carriers; n = 3 replicates per line, 18 replicates total). Number of differentially expressed isoforms (184 for engineered iNs and 213 donor derived iNs) mapped to 160 and 148 unique gene IDs respectively. Among the unique genes, example SYNGO-mapped genes are highlighted. c, d Percentage of isoform ratio change (differential transcript usage, see methods) for NRXN1 and CHD4 in NRXN1 deletions from engineered (eng.; yellow) and donor (green) backgrounds. Error bars represent +/−S.E.M. e Detection and visualization of local differential splicing using Majiq and ggSashimi. Change in percent spliced in index (dPSI) is calculated for differential exon usage in the NRXN1 gene (exons 19–21) in the SCZ-NRXN1 deletion vs. NRXN1 cKO engineered deletion iNs. Number of detected reads is plotted on the y-axis with location of genomic coordinates on the x-axis.CC BY · PubMed Central
Fig. 5 Impaired neuronal network activities in brain organoids carrying NRXN1 deletion. a Representative fluorescent images of intact donor derived organoids (Ctrl.2 and NRXN1del2) expressing soma-GCaMP6f2 at ~120 days prior to live Ca2+ imaging. b Colored raw intensity traces are shown in the boxed graphs with averaged intensities plotted in bolded black. c Averaged data for synchronous firing rates (number of detected synchronous spikes/minute) representative of network activity, as well as amplitudes (dF/F0) and frequencies (total number of detected peaks/minute) of spontaneous spike activity, are shown in scatter plots. Each data point represents averaged data from a single field of view (FOV) consisting of 4-5 ROIs per FOV. At least 4–6 FOVs were taken from each organoid and 5 organoids per genotype were used for experiments. N = 29–30 for Ctrl; N = 23–25 for NRXN1del2. d Experimental flow of live Ca2+ imaging with glycine stimulation. e Synchronous firing rate is plotted in genotypes with or without stimulation. Each data point represents averaged data from a single field of view (FOV) consisting of 4–5 ROIs per FOV. At least 4–6 FOVs were taken from each organoid and 5 organoids per genotype were used for experiments. N = 24–25 for Ctrl; N = 16–19 for NRXN1del2. Error bars represent +/− S.E.M. Statistical significance is represented by asterisks: *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001. Student’s t test (two-sided) was performed for data shown for c and one-way ANOVA with multiple comparisons (Tukey’s test) was used for data shown on e. Source data and statistics are provided as a Source Data file.
Fig. 5 Impaired neuronal network activities in brain organoids carrying NRXN1 deletion. a Representative fluorescent images of intact donor derived organoids (Ctrl.2 and NRXN1del2) expressing soma-GCaMP6f2 at ~120 days prior to live Ca2+ imaging. b Colored raw intensity traces are shown in the boxed graphs with averaged intensities plotted in bolded black. c Averaged data for synchronous firing rates (number of detected synchronous spikes/minute) representative of network activity, as well as amplitudes (dF/F0) and frequencies (total number of detected peaks/minute) of spontaneous spike activity, are shown in scatter plots. Each data point represents averaged data from a single field of view (FOV) consisting of 4-5 ROIs per FOV. At least 4–6 FOVs were taken from each organoid and 5 organoids per genotype were used for experiments. N = 29–30 for Ctrl; N = 23–25 for NRXN1del2. d Experimental flow of live Ca2+ imaging with glycine stimulation. e Synchronous firing rate is plotted in genotypes with or without stimulation. Each data point represents averaged data from a single field of view (FOV) consisting of 4–5 ROIs per FOV. At least 4–6 FOVs were taken from each organoid and 5 organoids per genotype were used for experiments. N = 24–25 for Ctrl; N = 16–19 for NRXN1del2. Error bars represent +/− S.E.M. Statistical significance is represented by asterisks: *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001. Student’s t test (two-sided) was performed for data shown for c and one-way ANOVA with multiple comparisons (Tukey’s test) was used for data shown on e. Source data and statistics are provided as a Source Data file.CC BY · PubMed Central
Fig. 6 Differential enrichment of disease-associated transcriptomic signatures. a Heatmap showing results from a subset of gene enrichment analyses of DEGs performed on cells of both brain organoid types (engineered and donor derived brain organoids) across all time points using neurological disorder dysregulated gene sets in several categories (ASD autism spectrum disorders, BP bipolar disorder, MDD mood disorder, SCZ schizophrenia) (columns). Data shown here are from 3.5 mo time point only. Body Mass Index (BMI) was used as control (−). Gene Set Enrichment Analysis (GSEA) was employed using the package GSEAPY, and significance scores were defined as −log10(FDRadjustedpvalues) to represent the associations between DEG sets and neurological disorders. Scores were trimmed to 0~10 (see Methods). Significance levels were represented by numbers of asterisks (*: adjusted p values < 0.05; **: adjusted p values < 0.01; ***: adjusted p values < 0.001). b Significance of differential expression of prioritized genes obtained from PGC wave 3 and SCHMEA consortium51,63. The size of each dot represents the level of DE significance of each gene in each cell class of 3.5 mo donor brain organoids (left; n = 4 for Ctrls and n = 4 for SCZ-NRXN1dels) and 3.5 mo engineered brain organoids (right; n = 2 for Flp and n = 2 for Cre). Two-sided Wilcoxon rank-sum tests were employed for the DE tests, with FDR-adjusted p-values utilized.
Fig. 6 Differential enrichment of disease-associated transcriptomic signatures. a Heatmap showing results from a subset of gene enrichment analyses of DEGs performed on cells of both brain organoid types (engineered and donor derived brain organoids) across all time points using neurological disorder dysregulated gene sets in several categories (ASD autism spectrum disorders, BP bipolar disorder, MDD mood disorder, SCZ schizophrenia) (columns). Data shown here are from 3.5 mo time point only. Body Mass Index (BMI) was used as control (−). Gene Set Enrichment Analysis (GSEA) was employed using the package GSEAPY, and significance scores were defined as −log10(FDRadjustedpvalues) to represent the associations between DEG sets and neurological disorders. Scores were trimmed to 0~10 (see Methods). Significance levels were represented by numbers of asterisks (*: adjusted p values < 0.05; **: adjusted p values < 0.01; ***: adjusted p values < 0.001). b Significance of differential expression of prioritized genes obtained from PGC wave 3 and SCHMEA consortium51,63. The size of each dot represents the level of DE significance of each gene in each cell class of 3.5 mo donor brain organoids (left; n = 4 for Ctrls and n = 4 for SCZ-NRXN1dels) and 3.5 mo engineered brain organoids (right; n = 2 for Flp and n = 2 for Cre). Two-sided Wilcoxon rank-sum tests were employed for the DE tests, with FDR-adjusted p-values utilized.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 34 references

Research works printed but not linked: 63

  1. Lowther, C. et al. Molecular characterization of NRXN1 deletions from 19,263 clinical microarray cases identifies exons important for neurodevelopmental disease expression. Genet Med. 19, 53–61 (2017). · DOI
  2. Fuccillo, M. V. & Pak, C. Copy number variants in neurexin genes: phenotypes and mechanisms. Curr. Opin. Genet. Dev. 68, 64–70 (2021). · DOI
  3. Sudhof, T. C. Synaptic neurexin complexes: a molecular code for the logic of neural circuits. Cell 171, 745–769 (2017). · DOI · Paper page
  4. Dai, J., Aoto, J. & Südhof, T. C. Alternative splicing of presynaptic neurexins differentially controls postsynaptic NMDA and AMPA receptor responses. Neuron 102, 993–1008.e5 (2019). · DOI
  5. Davatolhagh, M. F. & Fuccillo, M. V. Neurexin1α differentially regulates synaptic efficacy within striatal circuits. Cell Rep. 34, 108773 (2021). · DOI
  6. Chen, L. Y., Jiang, M., Zhang, B., Gokce, O. & Südhof, T. C. Conditional deletion of all neurexins defines diversity of essential synaptic organizer functions for neurexins. Neuron 94, 611–625.e4 (2017). · DOI
  7. Ullrich, B., Ushkaryov, Y. A. & Südhof, T. C. Cartography of neurexins: more than 1000 isoforms generated by alternative splicing and expressed in distinct subsets of neurons. Neuron 14, 497–507 (1995). · DOI
  8. Schreiner, D. et al. Targeted combinatorial alternative splicing generates brain region-specific repertoires of neurexins. Neuron 84, 386–398 (2014). · DOI
  9. Jenkins, A. K. et al. Neurexin 1 (NRXN1) splice isoform expression during human neocortical development and aging. Mol. Psychiatry 21, 701–706 (2016). · DOI
  10. Zeng, L. et al. Functional impacts of NRXN1 knockdown on neurodevelopment in stem cell models. PloS One 8, e59685 (2013). · DOI
  11. Lam, M. et al. Single cell analysis of autism patient with bi-allelic NRXN1-alpha deletion reveals skewed fate choice in neural progenitors and impaired neuronal functionality. Exp. Cell Res. 383, 111469 (2019). · DOI
  12. Willsey, H. R. et al. Parallel in vivo analysis of large-effect autism genes implicates cortical neurogenesis and estrogen in risk and resilience. Neuron 109, 788–804.e8 (2021). · DOI
  13. Xiang, Y. et al. Fusion of regionally specified hPSC-derived organoids models human brain development and interneuron migration. Cell Stem Cell 21, 383–398.e7 (2017). · DOI
  14. Sebastian, R., Song, Y. & Pak, C. Probing the molecular and cellular pathological mechanisms of schizophrenia using human induced pluripotent stem cell models. Schizophr. Res S0920-9964, 00263–00268 (2022).
  15. Gordon, A. et al. Long-term maturation of human cortical organoids matches key early postnatal transitions. Nat. Neurosci. 24, 331–342 (2021). · DOI
  16. Kathuria, A. et al. Transcriptomic landscape and functional characterization of induced pluripotent stem cell-derived cerebral organoids in schizophrenia. JAMA Psychiatry 77, 745–754 (2020). · DOI
  17. Pak, C. et al. Human neuropsychiatric disease modeling using conditional deletion reveals synaptic transmission defects caused by heterozygous mutations in NRXN1. Cell Stem Cell 17, 316–328 (2015). · DOI · Paper page
  18. Pak, C. et al. Cross-platform validation of neurotransmitter release impairments in schizophrenia patient-derived NRXN1-mutant neurons. Proc. Natl Acad. Sci. USA 118, e2025598118 (2021). · DOI · Paper page
  19. Tanaka, Y., Cakir, B., Xiang, Y., Sullivan, G. J. & Park, I.-H. Synthetic analyses of single-cell transcriptomes from multiple brain organoids and fetal brain. Cell Rep. 30, 1682–1689.e3 (2020). · DOI
  20. Zhong, S. et al. A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature 555, 524–528 (2018). · DOI · Paper page
  21. Bhaduri, A. et al. An atlas of cortical arealization identifies dynamic molecular signatures. Nature 598, 200–204 (2021). · DOI
  22. Kanton, S. et al. Organoid single-cell genomic atlas uncovers human-specific features of brain development. Nature 574, 418–422 (2019). · DOI · Paper page
  23. Uzquiano, A. et al. Proper acquisition of cell class identity in organoids allows definition of fate specification programs of the human cerebral cortex. Cell. 185, 3770–3788.e27 (2022). · DOI
  24. Chen, J., Bardes, E. E., Aronow, B. J. & Jegga, A. G. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 37, W305–W311 (2009). · DOI · Paper page
  25. Charizanis, K. et al. Muscleblind-like 2-mediated alternative splicing in the developing brain and dysregulation in myotonic dystrophy. Neuron 75, 437–450 (2012). · DOI
  26. Li, Q. et al. The splicing regulator PTBP2 controls a program of embryonic splicing required for neuronal maturation. Elife 3, e01201 (2014). · DOI · Paper page
  27. Lukacsovich, D. et al. Single-cell RNA-seq reveals developmental origins and ontogenetic stability of neurexin alternative splicing profiles. Cell Rep. 27, 3752–3759.e4 (2019). · DOI
  28. Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014). · DOI · Paper page
  29. Koopmans, F. et al. SynGO: an evidence-based, expert-curated knowledge base for the synapse. Neuron 103, 217–234.e4 (2019). · DOI · Paper page
  30. Vaquero-Garcia, J. et al. A new view of transcriptome complexity and regulation through the lens of local splicing variations. eLife 5, e11752 (2016). · DOI · Paper page
  31. Shemesh, O. A. et al. Precision calcium imaging of dense neural populations via a cell-body-targeted calcium indicator. Neuron 107, 470–486.e11 (2020). · DOI
  32. Watanabe, A. et al. Fabp7 maps to a quantitative trait locus for a schizophrenia endophenotype. PLoS Biol. 5, e297 (2007). · DOI
  33. Maekawa, M. et al. Polymorphism screening of brain-expressed FABP7, 5 and 3 genes and association studies in autism and schizophrenia in Japanese subjects. J. Hum. Genet. 55, 127–130 (2010). · DOI
  34. Koga, M. et al. Plasma fatty acid-binding protein 7 concentration correlates with depression/anxiety, cognition, and positive symptom in patients with schizophrenia. J. Psychiatr. Res. 144, 304–311 (2021). · DOI
  35. Hu, C., Chen, W., Myers, S. J., Yuan, H. & Traynelis, S. F. Human GRIN2B variants in neurodevelopmental disorders. J. Pharmacol. Sci. 132, 115–121 (2016). · DOI
  36. Takasaki, Y. et al. Mutation screening of GRIN2B in schizophrenia and autism spectrum disorder in a Japanese population. Sci. Rep. 6, 33311 (2016). · DOI
  37. Coyle, J. T. NMDA receptor and schizophrenia: a brief history. Schizophr. Bull. 38, 920–926 (2012). · DOI
  38. Nakazawa, K. & Sapkota, K. The origin of NMDA receptor hypofunction in schizophrenia. Pharm. Ther. 205, 107426 (2020). · DOI
  39. Shimamoto, C. et al. Functional characterization of FABP3, 5 and 7 gene variants identified in schizophrenia and autism spectrum disorder and mouse behavioral studies. Hum. Mol. Genet. 23, 6495–6511 (2014). · DOI
  40. Lu, W. et al. Activation of synaptic NMDA receptors induces membrane insertion of new AMPA receptors and LTP in cultured hippocampal neurons. Neuron 29, 243–254 (2001). · DOI
  41. Wang, L. et al. Analyses of the autism-associated neuroligin-3 R451C mutation in human neurons reveal a gain-of-function synaptic mechanism. Mol. Psychiatry https://doi.org/10.1038/s41380-022-01834-x (2022). · DOI
  42. Yu, Y. et al. Rare loss of function mutations in N-methyl-D-aspartate glutamate receptors and their contributions to schizophrenia susceptibility. Transl. Psychiatry 8, 1–9 (2018). · DOI
  43. Yeo, G., Holste, D., Kreiman, G. & Burge, C. B. Variation in alternative splicing across human tissues. Genome Biol. 5, R74 (2004). · DOI · Paper page
  44. Raj, B. & Blencowe, B. J. Alternative splicing in the mammalian nervous system: recent insights into mechanisms and functional roles. Neuron 87, 14–27 (2015). · DOI · Paper page
  45. Tan, W. et al. Molecular cloning of a brain-specific, developmentally regulated neuregulin 1 (NRG1) isoform and identification of a functional promoter variant associated with schizophrenia. J. Biol. Chem. 282, 24343–24351 (2007). · DOI
  46. Law, A. J., Kleinman, J. E., Weinberger, D. R. & Weickert, C. S. Disease-associated intronic variants in the ErbB4 gene are related to altered ErbB4 splice-variant expression in the brain in schizophrenia. Hum. Mol. Genet. 16, 129–141 (2007). · DOI
  47. Sartorius, L. J. et al. Expression of a GRM3 splice variant is increased in the dorsolateral prefrontal cortex of individuals carrying a schizophrenia risk SNP. Neuropsychopharmacology 33, 2626–2634 (2008). · DOI
  48. Flaherty, E. et al. Neuronal impact of patient-specific aberrant NRXN1α splicing. Nat. Genet. 51, 1679–1690 (2019). · DOI · Paper page
  49. Luza, S. et al. The ubiquitin proteasome system and schizophrenia. Lancet Psychiatry 7, 528–537 (2020). · DOI
  50. Nishimura, A. et al. The carbohydrate deposits detected by histochemical methods in the molecular layer of the dentate gyrus in the hippocampal formation of patients with schizophrenia, Down’s syndrome and dementia, and aged person. Glycoconj. J. 17, 815–822 (2000). · DOI
  51. Altar, C. A. et al. Deficient hippocampal neuron expression of proteasome, ubiquitin, and mitochondrial genes in multiple schizophrenia cohorts. Biol. Psychiatry 58, 85–96 (2005). · DOI · Paper page
  52. Rubio, M. D., Wood, K., Haroutunian, V. & Meador-Woodruff, J. H. Dysfunction of the ubiquitin proteasome and ubiquitin-like systems in schizophrenia. Neuropsychopharmacology 38, 1910–1920 (2013). · DOI
  53. Bousman, C. A. et al. Elevated ubiquitinated proteins in brain and blood of individuals with schizophrenia. Sci. Rep. 9, 2307 (2019). · DOI
  54. Tai, H.-C., Besche, H., Goldberg, A. L. & Schuman, E. M. Characterization of the brain 26S proteasome and its interacting proteins. Front. Mol. Neurosci. 3, 12 (2010).
  55. Ferreira, J. S. et al. GluN2B-containing NMDA receptors regulate AMPA receptor traffic through anchoring of the synaptic proteasome. J. Neurosci. 35, 8462–8479 (2015). · DOI · Paper page
  56. Huo, Y. et al. The deubiquitinating enzyme USP46 regulates AMPA receptor ubiquitination and trafficking. J. Neurochem. 134, 1067–1080 (2015). · DOI
  57. Fu, A. K. Y. et al. APCCdh1 mediates EphA4-dependent downregulation of AMPA receptors in homeostatic plasticity. Nat. Neurosci. 14, 181–189 (2011). · DOI
  58. Weinberger, D. R. Future of days past: neurodevelopment and schizophrenia. Schizophr. Bull. 43, 1164–1168 (2017). · DOI
  59. Raedler, T. J., Knable, M. B. & Weinberger, D. R. Schizophrenia as a developmental disorder of the cerebral cortex. Curr. Opin. Neurobiol. 8, 157–161 (1998). · DOI
  60. Garrido-Martín, D., Palumbo, E., Guigó, R. & Breschi, A. ggsashimi: Sashimi plot revised for browser- and annotation-independent splicing visualization. PLOS Comput. Biol. 14, e1006360 (2018). · DOI
  61. Wolf, F. A., Angerer, P. & Theis, F. J. SCANPY: large-scale single-cell gene expression data analysis. Genome Biol. 19, 15 https://doi.org/10.1186/s13059-017-1382-0 (2018). · DOI · Paper page
  62. Jin, K. et al. An interactive single cell web portal identifies gene and cell networks in COVID-19 host responses. iScience 24, 103115 https://doi.org/10.1016/j.isci.2021.103115 (2021). · DOI
  63. Kaimal, V. et al. ToppCluster: a multiple gene list feature analyzer for comparative enrichment clustering and network-based dissection of biological systems. Nucleic Acids Res. 38, W96–102 https://doi.org/10.1093/nar/gkq418 (2010). · DOI

Cited by

Metadata

Authors
Sebastian, R., Jin, K., Pavon, N., Bansal, R., Potter, A., Song, Y., Babu, J., Gabriel, R., Sun, Y., Aronow, B., & Pak, C.
Venue
Nature Communications
Publisher
Springer Science and Business Media LLC
Published
2023-06-24
Volume
14
Issue
1
ISSN
2041-1723
DOI
10.1038/s41467-023-39420-6
PMC
PMC10290702
Cited by
73
References
34
Language
en
Licence
Funders
  • U.S. Department of Health & Human Services | NIH | National Institute of Mental Health (MH122519, MH130843)
  • Tourette Association of America
  • U.S. Department of Health & Human Services | NIH | National Institute of General Medical Sciences (GM135096)

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review

Data sources

Bibliographic data: Semantic Scholar, OpenAlex, Crossref.

Contains information from the Semantic Scholar Academic Graph Datasets, made available under the ODC Attribution License (ODC-BY 1.0).

Semantic Scholar’s data is described in: Kinney, R., et al. (2023). The Semantic Scholar Open Data Platform. arXiv:2301.10140.