Next Article in Journal
Targeting Innate Immunity in Glioma Therapy
Previous Article in Journal
Chloroquine and Chemotherapeutic Compounds in Experimental Cancer Treatment
Previous Article in Special Issue
Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 2
10.3390/ijms25020946
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Association of Hippocampal Long-Term Potentiation-Induced Gene Expression with Genetic Risk for Psychosis

by
Natalie L. Wellard
1,
Nicholas E. Clifton
1,2,*,
Elliott Rees
1,
Kerrie L. Thomas
1 and
Jeremy Hall
1
1
Division of Psychological Medicine and Clinical Neurosciences, Cardiff University, Cardiff CF24 4HQ, UK
2
Medical School, Faculty of Health and Life Sciences, University of Exeter, Exeter EX4 4QF, UK
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(2), 946; https://doi.org/10.3390/ijms25020946
Submission received: 5 December 2023 / Revised: 8 January 2024 / Accepted: 11 January 2024 / Published: 12 January 2024
(This article belongs to the Special Issue Gene Regulation in Brain Development and Physiology)
Browse Figures
Versions Notes

Abstract

:
Genomic studies focusing on the contribution of common and rare genetic variants of schizophrenia and bipolar disorder support the view that substantial risk is conferred through molecular pathways involved in synaptic plasticity in the neurons of cortical and subcortical brain regions, including the hippocampus. Synaptic long-term potentiation (LTP) is central to associative learning and memory and depends on a pattern of gene expression in response to neuronal stimulation. Genes related to the induction of LTP have been associated with psychiatric genetic risk, but the specific cell types and timepoints responsible for the association are unknown. Using published genomic and transcriptomic datasets, we studied the relationship between temporally defined gene expression in hippocampal pyramidal neurons following LTP and enrichment for common genetic risk for schizophrenia and bipolar disorder, and for copy number variants (CNVs) and de novo coding variants associated with schizophrenia. We observed that upregulated genes in hippocampal pyramidal neurons at 60 and 120 min following LTP induction were enriched for common variant association with schizophrenia and bipolar disorder subtype I. At 60 min, LTP-induced genes were enriched in duplications from patients with schizophrenia, but this association was not specific to pyramidal neurons, perhaps reflecting the combined effects of CNVs in excitatory and inhibitory neuron subtypes. Gene expression following LTP was not related to enrichment for de novo coding variants from schizophrenia cases. Our findings refine our understanding of the role LTP-related gene sets play in conferring risk to conditions causing psychosis and provide a focus for future studies looking to dissect the molecular mechanisms associated with this risk.

1. Introduction

Genomic studies have provided critical insights into the molecular pathways contributing to risk for neuropsychiatric disorders. There is increasing evidence that many susceptibility genes for schizophrenia and related disorders impact synaptic plasticity mechanisms [1,2]; however, the nature of these associations, including the specific genes, physiological processes, cell types and timepoints implicated, remain to be further clarified.
In psychiatric genomics, genetic variants conferring risk are typically placed into one of two categories based on their prevalence: common variants with individually low effect sizes, typically identified through genome-wide association studies (GWASs), and rare variants, often with larger effect sizes, including copy number variants (CNVs) and rare coding variants. An extended analysis of the latest and largest schizophrenia GWAS meta-analysis from the Psychiatric Genomics Consortium (PGC3) identified common variants associated with schizophrenia in 287 distinct genomic loci [3]. Further genomic studies have identified rare CNVs at several loci that associate with schizophrenia, including deletions at 1q21.1, 3q29, 22q11.2 and NRXN1, and duplications at 16p11.2 [4,5,6,7,8]. Deleterious rare coding mutations associated with schizophrenia, including those arising de novo, have been identified through the exome sequencing of cases and controls [9,10,11].
Several of the genetic variants identified through common and rare variant analyses of schizophrenia impact synaptic processes [1]. In 2012, Kirov and colleagues combined CNV and proteomics datasets to conduct a systematic analysis of synaptic protein complexes [12]. They found an enrichment of post-synaptic density (PSD) genes within schizophrenia-associated CNVs, including a particular enrichment of genes involved in the N-methyl-D-aspartate (NMDA) receptor complex and activity-regulated cytoskeleton-associated protein (Arc) interactors. A further larger investigation confirmed the association between patient CNVs with NMDA receptor- and Arc-associated genes, and also found that CNVs were enriched more broadly for other regulators of synaptic function, including GABA receptors and calcium channels, and genes activated during synaptic plasticity and associative learning [13]. Recent evidence from a Schizophrenia Exome Meta-Analysis (SCHEMA) consortium study also implicates protein-damaging mutations affecting synaptic genes, including the NMDA receptor subunit GRIN2A, the AMPA receptor subunit GRIA3 and a synaptic voltage-gated calcium channel (CACNA1G) [10]. Large-scale common variant (GWAS) studies have also strongly implicated genes encoding synaptic proteins in risk for schizophrenia [3] and bipolar disorder [14,15]. Importantly, in some cases, converging evidence from both common and rare variant studies implicates the same synaptic genes, including GRIN2A in schizophrenia [1,3,10].
Through the use of curated databases such as Gene Ontology (GO), previous genomic studies have suggested an enrichment of long-term potentiation (LTP)-induced genes for genetic risk associated with schizophrenia mediated through both common variants [16] and rare copy number variants [13]. First described by Bliss and Lomo at excitatory synapses in the hippocampus [17], LTP is a key form of activity-dependent plasticity, giving rise to persistent increases in synaptic strength and connectivity that underly associative learning and memory [18]. While there are multiple forms of LTP discernible by their stimulation requirements for induction and initial cellular mechanisms for expression, a key distinction is that persistent LTP (late LTP or LTP3) depends on protein synthesis, whereas early LTP (or LTP/LTP2) not sustained beyond a few hours does not [19,20]. There is a general consensus that long-term memory requires de novo mRNA and protein synthesis [21,22]. Patterns of de novo gene expression persist for hours after LTP induction [23], but at which timepoints or in which cell types schizophrenia risk variants impact LTP-associated molecular pathways are unknown.
Several approaches have been used to investigate hippocampal gene expression following the induction of LTP. Early studies focused on candidate gene expression; however, as technology progressed, transcriptome-wide studies began to identify more extensive lists of genes expressed following LTP, including those involved in transmitter transport, growth factors, ion channels, Arc interactors and regulation of the cytoskeleton [24,25,26,27,28]. In particular, the advent of Translating Ribosome Affinity Purification (TRAP) has allowed for the profiling of translating RNA from select populations of cells [29,30]. TRAP by RiboTag [31] involves genetically tagging a ribosomal protein under the control of a tissue-specific promotor and isolating ribosome-associated mRNAs by immunoprecipitation for transcriptome profiling. Chen and colleagues (2017) utilized this technique to profile the genes expressed 30, 60 and 120 min after chemically induced LTP in hippocampal slices [23] (Supplementary Figure S1). The ribosome tag was driven by a αCaMKII promotor [32] to profile LTP gene expression specifically in CA1 excitatory pyramidal neurons, a cell type genetically associated with schizophrenia through their unique repertoire of gene expression [33]. They highlight a range of synaptic and non-synaptic pathways stimulated by LTP and present specific sets of genes expressed at each timepoint.
Here, we explored the relationship between common and rare genomic risk factors for schizophrenia and related disorders with genes expressed at different timepoints after LTP induction in mouse hippocampal CA1 excitatory neurons [23].

2. Results

2.1. Common Variant Association

Transcriptomic changes following LTP were obtained from published datasets derived from mouse hippocampal CA1 using total RNA-seq and cell-type-specific TRAP-seq [23]. To assess the enrichment in LTP-related genes for common genetic association with schizophrenia and related disorders, we used GWAS summary statistics to perform gene set enrichment analyses [34] on sets of genes captured at distinct timepoints following LTP (Supplementary Figure S1), conditioning on all tissue-expressed genes. There was a significant enrichment for association in genes identified through TRAP-seq at 60 and 120 min after LTP induction but not at 30 min (30 min: β = −0.032, Bonferroni-adjusted p (p.bonf) = 1.0; 60 min: β = 0.12, p.bonf = 0.0017; 120 min: β = 0.17, p.bonf = 2.1 × 10−5) (Figure 1). At both 60 min and 120 min, the association was driven by the upregulated subset (60 min: β = 0.18, p.bonf = 3.0 × 10−4; 120 min: β = 0.17, 2.9 × 10−4), whilst downregulated genes were not significantly enriched for association (60 min: β = 0.037, p.bonf = 1.0; 120 min: β = 0.17, p.bonf = 0.10). Among the differentially expressed genes (DEGs) derived from total RNA-seq, there was no evidence of enrichment for common schizophrenia risk at any time point (30 min: β = −0.23, p.bonf = 1.0; 60 min: β = 0.076, p.bonf = 1.0, 120 min: β = 0.067, p.bonf = 0.75).
Similarly, gene sets identified through TRAP-seq at 60 min and 120 min after LTP were also associated with bipolar disorder risk (TRAP-seq: 30 min: β = −0.17, p.bonf = 1.0, 60 min: β = 0.11, p.bonf = 7.8 × 10−4, 120 min: β = 0.095 p.bonf = 0.017). Whilst the association at 120 min was driven by upregulated genes (upregulated: β = 0.12, p.bonf = 0.0041; downregulated: β = −0.019, p.bonf = 1.0), the association at 60 min was prevalent in both up- and downregulated genes (upregulated: β = 0.10, p.bonf = 0.033; downregulated: β = 0.10, p.bonf = 0.041). When the two bipolar subtypes were analysed separately, we observed that the association signal from LTP-induced genes was restricted to bipolar subtype I (Table 1). In addition, there was no significant enrichment of any gene sets for common variant association with chronic kidney disease (TRAP-seq 30 min: β = −0.024, p.bonf = 1.0, 60 min: β = 0.0061, p.bonf = 1.0, 120 min: β = −0.0029, p.bonf = 1.0; Total RNA-seq 30 min: β = −0.11, p.bonf = 1.0, 60 min: β = −0.11, p.bonf = 1.0, 120 min: β = 0.033, p.bonf = 1.0).

2.2. Schizophrenia CNV Enrichment

To investigate the enrichment of LTP gene sets in schizophrenia-associated CNVs, we performed CNV enrichment analysis using deletions and duplications from cases and controls. Similar to the common variant enrichment, genes expressed 60 min post-LTP in CA1 excitatory neurons captured by TRAP were found to be enriched in CNVs from patients (β = 0.12, p.unadjusted = 8.75 × 10−6, p.adjusted = 0.025) (Figure 2), again driven by the upregulated subset (β = 0.14, p.unadjusted = 1.4 × 10−4, p.adjusted = 0.046). Of genes deriving from total RNAseq, those differentially expressed at 30 min and 60 min were also enriched in case CNVs (30 min: β = 1.12, p.unadjusted = 0.0063, p.adjusted = 0.046; 60 min: β = 0.55, p.unadjusted = 5.1 × 10−4, p.adjusted = 0.029) (Figure 2).
Deletion and duplication CNVs were then analysed separately to investigate whether our results were driven by a particular subset of CNVs. Interestingly, only the 60 min downregulated gene set from TRAP-seq was significantly enriched for deletions (β = 0.25, p.unadjusted = 2.0 × 10−5, p.adjusted = 0.045), whilst the 60 min upregulated genes from TRAP-seq or total RNA-seq were enriched for duplications (TRAP-seq: β = 0.21, p.unadjusted = 1.3 × 10−6, p.adjusted = 0.013; total RNA-seq: β = 0.96, p.unadjusted = 6.8 × 10−7, p.adjusted = 0.016).

2.3. Schizophrenia De Novo Rare Coding Variant Enrichment

Lastly, the enrichment of de novo rare coding variants from schizophrenia cases in the LTP gene sets was examined, focusing on loss-of-function protein truncating variants (PTVs) and missense mutations with an MPC score of greater than 2. TRAP-acquired genes induced 120 min after LTP were enriched for case variants of these classes with nominal significance, but there was no significant enrichment of case variants in the TRAP- or total RNA-seq conditions when corrected for multiple comparisons (Table 2).

3. Discussion

In the present study, we tested sets of LTP-related genes for genetic association with schizophrenia and bipolar disorder. The analysis of transcriptomes from total tissue and CA1 pyramidal neurons specifically, and at multiple experimental timepoints, allowed for an in-depth investigation into the manner in which psychiatric risk variants impact on LTP-induced genes. Our primary finding was an enrichment of LTP-induced genes from 60 and 120 min timepoints in common genetic variants for both schizophrenia and bipolar disorder type I and also an enrichment of LTP-induced genes for CNVs in schizophrenia.
Previous studies of psychiatric genomics found an association between schizophrenia risk and genetic pathways implicating synaptic function, including LTP-related gene sets [3,5,12,13,16,35]. Here, we expand on these findings by demonstrating the temporal, spatial and directional specificity of these associations. We report that there was an enrichment of LTP-associated genes for association with schizophrenia, as conferred through common variation, at 60 and 120 min post-LTP-induction. This association was restricted to upregulated genes and was only identified using data from tissue-specific TRAP-seq and not that from total hippocampal samples, suggesting that LTP-regulated transcripts implicated in schizophrenia through common variant association are specifically enriched in CA1 pyramidal neurons. These results are consistent with previous studies, showing that CA1 pyramidal cells are impacted by schizophrenia risk variants [33]. Furthermore, the enrichment of schizophrenia risk variants, particularly those with synaptic annotations, in mouse and human CA1 [3] and in LTP-regulated transcripts in CA1 (here) highlight the functional impact of the risk variants on plasticity mechanisms underlying specific learning and memory functions in the mature brain. The selective enrichment highlights the benefit of using the tissue-specific TRAP method to interrogate LTP-related transcriptomics and identify a localised association in a functional setting.
The transition to persistent forms of LTP and from short-term to long-term memory involves the transcription of genes and proteins that are associated with synaptic growth [36,37]. The distinctive, temporal progression of gene expression measured by either total RNA seq or TRAP-Seq by Chen and colleagues (2017) brings to light the co-ordinated programme of gene expression triggered by LTP induction. These data are consistent with the early regulation of genes such as transcription factors that have effector roles in the subsequent expression of genes, supporting the synaptic structural changes associated with long-term synaptic plasticity and memory storage [36]. The enrichment of LTP-associated genes for association with schizophrenia at 60 and 120 min but not the earliest time point investigated concurs with the observation that synaptic proteins, consequential for maintaining plasticity and memory, are also enriched for association with schizophrenia [3]. It is possible that the impact of schizophrenia risk variants on adult plasticity may preferentially affect the magnitude and persistence of synaptic changes after induction or following learning.
Gene expression captured in CA1 pyramidal neurons 60 and 120 min following LTP induction was similarly associated with risk for bipolar disorder. Furthermore, when the GWAS data were analysed separately for each bipolar disorder subtype, we observed that the association was restricted to bipolar I cases. Clinically and genomically, bipolar disorder type I is more strongly related to schizophrenia, consistent with their shared association with LTP-induced genes, whereas bipolar disorder II is more aligned with the clinical and genomic characteristics of major depressive disorder [14,15,38]. Common variants associated with bipolar disorder were recently found to be enriched in genes expressed in the hippocampus, including hippocampal pyramidal neurons, and in GO terms relating to synaptic signalling [15]. Thus, the current finding adds to the growing literature implicating dynamic synaptic processes in the aetiology of both schizophrenia and bipolar disorder, although we note that genetic discovery sample sizes were larger for bipolar disorder type I than bipolar disorder type II, potentially limiting power to identify associations with bipolar disorder type II. By testing for associations in a chronic kidney disease GWAS of comparable sample size to the schizophrenia and bipolar disorder GWAS, we provide further confidence that our primary results were not inflated by methodological artefacts of using large GWAS samples.
To further investigate the relationship between LTP and schizophrenia genetic risk, we examined the enrichment of LTP-induced genes in case CNVs. Previously, CNVs in people with schizophrenia were shown to be enriched for genes involved in synaptic plasticity, including those encoding Arc interactors and NMDA receptor complexes, and genes associated with abnormal synaptic transmission and LTP [12,13]. In combined analyses of both deletion and duplication CNVs, an enrichment of genes upregulated 60 min after LTP was observed. The enrichment was mirrored in gene sets derived from both total and TRAP RNA-seq. This could reflect a lack of specificity of the risk conferred by CNVs to CA1 pyramidal neurons, instead impacting a range of hippocampal cell types involved in LTP induction. Alternatively, the enrichment in case CNVs may be driven by a small number of genes with large expression changes in pyramidal neurons, detectable in total RNA-seq. By separating analyses of CNV subtypes, we found that the enrichment of LTP-induced genes in case CNVs was particularly driven by duplications.
Taken together, these results suggest that genetic risk for schizophrenia, as conferred through CNVs, impacts a different subset of LTP-related genes compared to common variants. The difference is firstly highlighted by the temporal pattern of enrichment, whereby CNVs (duplications) were enriched for genes expressed at 60 min but not 120 min. Whilst 60 min was also characterised by the largest gene set, permutation analyses adjusted for this. Secondly, the specificity of the genetic association to TRAP-derived gene sets observed in analyses of GWAS data was not seen in tests of CNV enrichment. Considering previous findings demonstrating a contribution of genes involved in inhibitory synaptic complexes to the risk associated with CNVs [13], our results may reflect differences in the contribution of inhibitory neurons to the association signals derived from case CNVs vs. common variants. However, some of the association attributable to CNVs may be missed through gaps in probe coverage and filtering of CNVs by size, which could remove the signal from some genes also impacted by common variants.
Further inconsistencies across types of genetic variant emerged when examining the association of LTP-related genes sets with schizophrenia through de novo coding variants. We found only nominal enrichment for de novo coding variants from schizophrenia probands in LTP-induced genes at 120 min from TRAP but no overall association regardless of timepoint or transcriptomic approach. Previous research reported that de novo nonsynonymous mutations from patients with schizophrenia are enriched in synaptic gene sets, the postsynaptic density and associated protein complexes [39]. That we observe no such enrichment in LTP sets suggests that these variants confer risk through the disruption of synaptic pathways other than those related to potentiation. However, the current study of de novo variation might be underpowered to detect these effects. Increased sample sizes from exome sequencing will contribute to increased confidence in understanding how de novo coding variants play a role in modulating plasticity mechanisms in schizophrenia.
Our interpretations of these data assume that post-LTP transcriptomic changes in humans sufficiently mirror those observed in mice. Whilst there may be small differences in the temporal patterns or precise transcripts expressed following LTP, the composition and function of neuronal synapses are strongly conserved between the species [40,41]. Since our conclusions are based on genetic associations of groups of co-expressed genes and not individual genes, minor effects of species are not expected to alter our main findings.
Our findings build on our understanding of genetic associations in LTP-related gene sets in schizophrenia and bipolar disorder by refining these associations according to their temporal expression patterns, cell types and classes of genomic risk variant. In the case of bipolar disorder, we provide additional specificity to the disorder subtype in which the associations manifest, supporting an association specifically with bipolar disorder type I. The data highlight a reduced set of LTP-induced genes for the focus of downstream studies into plasticity mechanisms relevant to the development of psychiatric disorders.

4. Materials and Methods

4.1. LTP-Induced Gene Sets

RNA-sequencing datasets were generated by Chen and colleagues [23]. Briefly, RiboTag mice were generated with a floxed HA-tagged ribosomal protein L22 (The Jackson Laboratory: 011029) under CaMKIIα-cre mediated expression (The Jackson Laboratory: 005359) to direct expression to the CA1 pyramidal cell layer of the hippocampus. At 8–12 weeks old, a chemical LTP induction protocol was applied to ex vivo hippocampal slices with the dentate gyrus (DG) removed, using forskolin [42,43]. Chen and colleagues sequenced ribosome-bound mRNA from CA1 excitatory neurons and total mRNA from bulk hippocampal tissue at baseline and at three timepoints after LTP induction (30 min, 60 min and 120 min).
Processed RNA sequencing datasets were obtained from the NCBI Gene Expression Omnibus (GEO), accession number GSE79790. These datasets contained the log fold change (logFC) and false-discovery rate (FDR) for each gene, for each contrast of interest. These contrasts were: “TRAP-LTP vs. TRAP-basal” (TRAP-seq) and “Total-LTP vs. Total-basal” (total RNA-seq), for each of the three time points. These datasets were further filtered to include only those genes with an FDR < 0.01 to form the gene sets. Note that this significance threshold is more stringent than that used in the original study [23]. To facilitate comparisons with human data, mouse Entrez IDs were converted to their human homologs. Genes that did not have a unique human homolog were excluded from further analyses.
After filtering of total RNA-seq data, there were 21 differentially expressed genes (DEGs) at 30 min, 97 DEGs at 60 min and 474 DEGs at 120 min following LTP. From TRAP-seq data, there were 67 DEGs at 30 min, 1499 DEGs at 60 min and 1082 DEGs at 120 min (Supplementary Figure S1).

4.2. Genotype Data

SNP genotype data were obtained from published case and control samples. Schizophrenia summary statistics were acquired from the primary meta-analysis from wave 3 of the psychiatric genomics consortium, using genotype data from 74,776 case and 101,023 control subjects of European, East Asian, African American and Latino ancestry [3,44]. Summary statistics for bipolar disorder, and each subtype, were acquired from a recent study of 41,917 cases and 371,549 controls of European ancestry [15]. A bipolar disorder type I meta-analysis was performed using 25,060 cases and 449,978 controls, and a bipolar disorder type II meta-analysis using 6781 cases and 364,075 controls [15]. Lastly, GWAS summary statistics from a study of chronic kidney disease, used for control purposes, derived from a meta-analysis of 41,395 cases and 439,303 controls of European ancestry [45].
Ultra-rare de novo coding variants from 3444 schizophrenia proband parent trios were obtained from published exome sequencing studies [46,47,48,49,50,51,52,53,54]. Protein truncating variants (PTVs), defined as frameshift variants, stop-gain variants and donor/acceptor splice site variants, as well as missense variants with a “Missense badness, Polyphen-2, constraint” (MPC) score of 2 or greater, were taken forward for analysis [55].
Copy number variants, including deletions and duplications, from schizophrenia cases and controls were collated from CLOZUK and Cardiff cognition in Schizophrenia samples [6,8], International Schizophrenia Consortium [56] and Molecular Genetics of Schizophrenia [57] samples, totalling 17,565 cases and 24,830 control subjects.

4.3. Gene Set Enrichment Analysis

Gene set enrichment analysis for common variation was undertaken using Multi-marker Analysis of GenoMic Annotation (MAGMA) [34]. GWAS summary statistic SNP files were filtered to remove SNP with an info score of less than 0.8 or an allele frequency of less than 1%. SNPs were mapped to genes using the NCBI build 37 gene location files as a reference, and gene-wide analysis performed with linkage disequilibrium calculated using the 1000 genomes reference file (European ancestry). The resulting gene-wide p values were used in gene set enrichment analyses, with all expressed genes included as a covariate within the model to control for any background association from brain-expressed genes. Expressed genes were defined independently for TRAP and total RNA-seq datasets and were defined as having a normalised expression count greater than 10 in at least 3 basal samples. Gene set p-values were subjected to Bonferroni correction where appropriate.
De novo coding variant gene set enrichment was determined using a two sample Poisson rate ratio test, which compared the observed rate of variants in the gene set with the expected rate under the null (https://github.com/reeseg/SZ_NDD_pleiotropy_analysis (accessed on 8 March 2023)), while correcting for the background enrichment of these variants from brain-expressed genes. The expected number of de novo coding variants was estimated from per-gene mutation rates [58].
CNV datasets were filtered to only include those larger than 100 kb in size and covered by at least 15 probes. Protein coding genes overlapping each CNV were identified using the NCBI build appropriate to the dataset (ISC: build 35, MGS: build 36, CLOZUK: build 37), and the number of genes overlapping each CNV was counted. Logistic regression analysis was used to investigate the relative enrichment of a set of genes in case CNVs compared to controls using a general linear model:
Case–control status ~ number of genes in CNV + CNV size + number of
probes + CNV study + number of gene set hits
To adjust for p-value inflation with larger gene set sizes, permutation correction was applied. The background gene set was permutated 2000 times to generate a null distribution of p-values from random size-matched sets of brain-expressed genes. This procedure yielded an empirical p-value for each LTP gene set.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25020946/s1.

Author Contributions

Conceptualization, K.L.T. and J.H.; methodology, N.E.C. and E.R.; software, E.R.; validation, N.L.W. and E.R.; formal analysis, N.L.W. and N.E.C.; investigation, N.L.W. and N.E.C.; resources, N.E.C. and E.R.; data curation, N.L.W., N.E.C. and E.R.; writing—original draft preparation, N.L.W.; writing—review and editing, N.E.C., K.L.T. and J.H.; visualization, N.L.W. and N.E.C.; supervision, K.L.T., J.H. and N.E.C.; project administration, N.E.C. and J.H.; funding acquisition, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by a Wellcome Trust PhD studentship to N.L.W. E.R. was funded by a UKRI Future Leaders Fellowship Grant (MR/T018712/1). N.E.C. was supported by a UKRI Career Development Award (MR/W017156/1).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and cited works.

Acknowledgments

We thank the Psychiatric Genomics Consortium for making available genomics datasets used in this study. Analyses of copy number variation described in this paper used datasets from the Molecular Genetics of Schizophrenia (MGS; dbGAP phs000021.v3.p2 and phs000167.v1.p1) and the International Schizophrenia Consortium (ISC). The CLOZUK datasets contain data obtained from outside sources: dbGaP phs000404.v1.p1, phs000187.v1. p1, phs000303.v1.p1, phs000179.v3.p2, phs000421.v1.p, phs000395.v1. p1, phs000519.v1.p1 and the Wellcome Trust Case Control Consortium 2 study.

Conflicts of Interest

E.R. reports receiving a grant from Akrivia Health outside the submitted work. Akrivia played no part in the conception, design, implementation, or interpretation of this study. J.H. reports receiving a grant from Takeda Pharmaceutical Company outside the submitted work. Takeda played no part in the conception, design, implementation, or interpretation of this study.

References

  1. Hall, J.; Bray, N.J. Schizophrenia genomics: Convergence on synaptic development, adult synaptic plasticity, or both? Biol. Psychiatry 2022, 91, 709–717. [Google Scholar] [CrossRef] [PubMed]
  2. Mould, A.W.; Hall, N.A.; Milosevic, I.; Tunbridge, E.M. Targeting synaptic plasticity in schizophrenia: Insights from genomic studies. Trends Mol. Med. 2021, 27, 1022–1032. [Google Scholar] [CrossRef] [PubMed]
  3. Trubetskoy, V.; Pardiñas, A.F.; Qi, T.; Panagiotaropoulou, G.; Awasthi, S.; Bigdeli, T.B.; Bryois, J.; Chen, C.-Y.; Dennison, C.A.; Hall, L.S.; et al. Mapping genomic loci implicates genes and synaptic biology in schizophrenia. Nature 2022, 604, 502–508. [Google Scholar] [CrossRef] [PubMed]
  4. Green, T.; Gothelf, D.; Glaser, B.; Debbane, M.; Frisch, A.; Kotler, M.; Weizman, A.; Eliez, S. Psychiatric Disorders and Intellectual Functioning Throughout Development in Velocardiofacial (22q11.2 Deletion) Syndrome. J. Am. Acad. Child Adolesc. Psychiatry 2009, 48, 1060–1068. [Google Scholar] [CrossRef] [PubMed]
  5. Marshall, C.R.; Howrigan, D.P.; Merico, D.; Thiruvahindrapuram, B.; Wu, W.; Greer, D.S.; Antaki, D.; Shetty, A.; Holmans, P.A.; Pinto, D.; et al. Contribution of copy number variants to schizophrenia from a genome-wide study of 41,321 subjects. Nat. Genet. 2017, 49, 27–35. [Google Scholar] [CrossRef]
  6. Rees, E.; Walters, J.T.R.; Georgieva, L.; Isles, A.R.; Chambert, K.D.; Richards, A.L.; Mahoney-Davies, G.; Legge, S.E.; Moran, J.L.; McCarroll, S.A.; et al. Analysis of copy number variations at 15 schizophrenia-associated loci. Br. J. Psychiatry 2014, 204, 108–114. [Google Scholar] [CrossRef] [PubMed]
  7. Stefansson, H.; Rujescu, D.; Cichon, S.; Pietiläinen, O.P.H.; Ingason, A.; Steinberg, S.; Fossdal, R.; Sigurdsson, E.; Sigmundsson, T.; Buiz er-Voskamp, J.E.; et al. Large recurrent microdeletions associated with schizophrenia. Nature 2008, 455, 232–236. [Google Scholar] [CrossRef]
  8. Rees, E.; Kendall, K.; Pardiñas, A.F.; Legge, S.E.; Pocklington, A.; Escott-Price, V.; MacCabe, J.H.; Collier, D.A.; Holmans, P.; O’donovan, M.C.; et al. Analysis of Intellectual Disability Copy Number Variants for Association With Schizophrenia. JAMA Psychiatry 2016, 73, 963–969. [Google Scholar] [CrossRef]
  9. Rees, E.; Carrera, N.; Morgan, J.; Hambridge, K.; Escott-Price, V.; Pocklington, A.J.; Richards, A.L.; Pardiñas, A.F.; McDonald, C.; Donohoe, G.; et al. Targeted Sequencing of 10,198 Samples Confirms Abnormalities in Neuronal Activity and Implicates Voltage-Gated Sodium Channels in Schizophrenia Pathogenesis. Biol. Psychiatry 2019, 85, 554–562. [Google Scholar] [CrossRef]
  10. Singh, T.; Poterba, T.; Curtis, D.; Akil, H.; Al Eissa, M.; Barchas, J.D.; Bass, N.; Bigdeli, T.B.; Breen, G.; Bromet, E.J.; et al. Rare coding variants in ten genes confer substantial risk for schizophrenia. Nature 2022, 604, 509–516. [Google Scholar] [CrossRef]
  11. Rees, E.; Creeth, H.D.J.; Hwu, H.-G.; Chen, W.J.; Tsuang, M.; Glatt, S.J.; Rey, R.; Kirov, G.; Walters, J.T.R.; Holmans, P.; et al. Schizophrenia, autism spectrum disorders and developmental disorders share specific disruptive coding mutations. Nat. Commun. 2021, 12, 5353. [Google Scholar] [CrossRef]
  12. Kirov, G.; Pocklington, A.J.; Holmans, P.; Ivanov, D.; Ikeda, M.; Ruderfer, D.; Moran, J.; Chambert, K.; Toncheva, D.; Georgieva, L.; et al. De novo CNV analysis implicates specific abnormalities of postsynaptic signalling complexes in the pathogenesis of schizophrenia. Mol. Psychiatry 2012, 17, 142–153. [Google Scholar] [CrossRef] [PubMed]
  13. Pocklington, A.J.; Rees, E.; Walters, J.T.R.; Han, J.; Kavanagh, D.H.; Chambert, K.D.; Holmans, P.; Moran, J.L.; McCarroll, S.A.; Kirov, G.; et al. Novel Findings from CNVs Implicate Inhibitory and Excitatory Signaling Complexes in Schizophrenia. Neuron 2015, 86, 1203–1214. [Google Scholar] [CrossRef] [PubMed]
  14. Stahl, E.A.; Breen, G.; Forstner, A.J.; McQuillin, A.; Ripke, S.; Trubetskoy, V.; Mattheisen, M.; Wang, Y.; Coleman, J.R.I.; Gaspar, H.A.; et al. Genome-wide association study identifies 30 loci associated with bipolar disorder. Nat. Genet. 2019, 51, 793–803. [Google Scholar] [CrossRef]
  15. Mullins, N.; Forstner, A.J.; O’connell, K.S.; Coombes, B.; Coleman, J.R.I.; Qiao, Z.; Als, T.D.; Bigdeli, T.B.; Børte, S.; Bryois, J.; et al. Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology. Nat. Genet. 2021, 53, 817–829. [Google Scholar] [CrossRef] [PubMed]
  16. Pardiñas, A.F.; Holmans, P.; Pocklington, A.J.; Escott-Price, V.; Ripke, S.; Carrera, N.; Legge, S.E.; Bishop, S.; Cameron, D.; Hamshere, M.L.; et al. Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. Nat. Genet. 2018, 50, 381–389. [Google Scholar] [CrossRef]
  17. Bliss, T.V.P.; Lømo, T. Long-lasting potentiation of synaptic transmission in the dentate area of the anaesthetized rabbit following stimulation of the perforant path. J. Physiol. 1973, 232, 331–356. [Google Scholar] [CrossRef]
  18. Whitlock, J.R.; Heynen, A.J.; Shuler, M.G.; Bear, M.F. Learning Induces Long-Term Potentiation in the Hippocampus. Science 2006, 313, 1093–1097. [Google Scholar] [CrossRef]
  19. Frey, U.; Krug, M.; Reymann, K.G.; Matthies, H. Anisomycin, an inhibitor of protein synthesis, blocks late phases of LTP phenomena in the hippocampal CA1 region in vitro. Brain Res. 1988, 452, 57–65. [Google Scholar] [CrossRef]
  20. Bradshaw, K.D.; Emptage, N.J.; Bliss, T.V.P. A role for dendritic protein synthesis in hippocampal late LTP. Eur. J. Neurosci. 2003, 18, 3150–3152. [Google Scholar] [CrossRef]
  21. Redondo, R.L.; Morris, R.G.M. Making memories last: The synaptic tagging and capture hypothesis. Nat. Rev. Neurosci. 2011, 12, 17–30. [Google Scholar] [CrossRef] [PubMed]
  22. Asok, A.; Leroy, F.; Rayman, J.B.; Kandel, E.R. Molecular Mechanisms of the Memory Trace. Trends Neurosci. 2019, 42, 14–22. [Google Scholar] [CrossRef]
  23. Chen, P.B.; Kawaguchi, R.; Blum, C.; Achiro, J.M.; Coppola, G.; O’Dell, T.J.; Martin, K.C. Mapping Gene Expression in Excitatory Neurons during Hippocampal Late-Phase Long-Term Potentiation. Front. Mol. Neurosci. 2017, 10, 39. [Google Scholar] [CrossRef] [PubMed]
  24. Thompson, K.J.; Orfila, J.E.; Achanta, P.; Martinez, J.L. Gene expression associated with in vivo induction of early phase-long-term potentiation (LTP) in the hippocampal mossy fiber-Cornus Ammonis (CA)3 pathway. Cell. Mol. Biol. 2003, 49, 1281–1287. [Google Scholar]
  25. Wibrand, K.; Messaoudi, E.; Håvik, B.; Steenslid, V.; Løvlie, R.; Steen, V.M.; Bramham, C.R. Identification of genes co-upregulated with Arc during BDNF-induced long-term potentiation in adult rat dentate gyrus in vivo. Eur. J. Neurosci. 2006, 23, 1501–1511. [Google Scholar] [CrossRef]
  26. Valor, L.M.; Barco, A. Hippocampal gene profiling: Toward a systems biology of the hippocampus. Hippocampus 2012, 22, 929–941. [Google Scholar] [CrossRef]
  27. Ryan, M.M.; Ryan, B.; Kyrke-Smith, M.; Logan, B.; Tate, W.P.; Abraham, W.C.; Williams, J.M. Temporal Profiling of Gene Networks Associated with the Late Phase of Long-Term Potentiation In Vivo. PLoS ONE 2012, 7, e40538. [Google Scholar] [CrossRef] [PubMed]
  28. Park, C.S.; Gong, R.; Stuart, J.; Tang, S.-J. Molecular Network and Chromosomal Clustering of Genes Involved in Synaptic Plasticity in the Hippocampus. J. Biol. Chem. 2006, 281, 30195–30211. [Google Scholar] [CrossRef]
  29. Doyle, J.P.; Dougherty, J.D.; Heiman, M.; Schmidt, E.F.; Stevens, T.R.; Ma, G.; Bupp, S.; Shrestha, P.; Shah, R.D.; Doughty, M.L.; et al. Application of a Translational Profiling Approach for the Comparative Analysis of CNS Cell Types. Cell 2008, 135, 749–762. [Google Scholar] [CrossRef]
  30. Haimon, Z.; Volaski, A.; Orthgiess, J.; Boura-Halfon, S.; Varol, D.; Shemer, A.; Yona, S.; Zuckerman, B.; David, E.; Chappell-Maor, L.; et al. Re-evaluating microglia expression profiles using RiboTag and cell isolation strategies. Nat. Immunol. 2018, 19, 636–644. [Google Scholar] [CrossRef]
  31. Sanz, E.; Bean, J.C.; Carey, D.P.; Quintana, A.; McKnight, G.S. RiboTag: Ribosomal Tagging Strategy to Analyze Cell-Type-Specific mRNA Expression In Vivo. Curr. Protoc. Neurosci. 2019, 88, e77. [Google Scholar] [CrossRef] [PubMed]
  32. Tsien, J.Z.; Chen, D.F.; Gerber, D.; Tom, C.; Mercer, E.H.; Anderson, D.J.; Mayford, M.; Kandel, E.R.; Tonegawa, S. Subregion- and Cell Type–Restricted Gene Knockout in Mouse Brain. Cell 1996, 87, 1317–1326. [Google Scholar] [CrossRef] [PubMed]
  33. Skene, N.G.; Bryois, J.; Bakken, T.E.; Breen, G.; Crowley, J.J.; Gaspar, H.A.; Giusti-Rodriguez, P.; Hodge, R.D.; Miller, J.A.; Muñoz-Manchado, A.B.; et al. Genetic identification of brain cell types underlying schizophrenia. Nat. Genet. 2018, 50, 825–833. [Google Scholar] [CrossRef] [PubMed]
  34. de Leeuw, C.A.; Mooij, J.M.; Heskes, T.; Posthuma, D. MAGMA: Generalized Gene-Set Analysis of GWAS Data. PLoS Comput. Biol. 2015, 11, e1004219. [Google Scholar] [CrossRef] [PubMed]
  35. Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature 2014, 511, 421–427. [Google Scholar] [CrossRef]
  36. Alberini, C.M.; Kandel, E.R. The Regulation of Transcription in Memory Consolidation. Cold Spring Harb. Perspect. Biol. 2014, 7, a021741. [Google Scholar] [CrossRef]
  37. Yang, Y.; Liu, J.-J. Structural LTP: Signal transduction, actin cytoskeleton reorganization, and membrane remodeling of dendritic spines. Curr. Opin. Neurobiol. 2022, 74, 102534. [Google Scholar] [CrossRef]
  38. Lee, P.H.; Feng, Y.-C.A.; Smoller, J.W. Pleiotropy and Cross-Disorder Genetics Among Psychiatric Disorders. Biol. Psychiatry 2021, 89, 20–31. [Google Scholar] [CrossRef]
  39. Fromer, M.; Pocklington, A.J.; Kavanagh, D.H.; Williams, H.J.; Dwyer, S.; Gormley, P.; Georgieva, L.; Rees, E.; Palta, P.; Ruderfer, D.M.; et al. De novo mutations in schizophrenia implicate synaptic networks. Nature 2014, 506, 179–184. [Google Scholar] [CrossRef]
  40. Pembroke, W.G.; Hartl, C.L.; Geschwind, D.H. Evolutionary conservation and divergence of the human brain transcriptome. Genome Biol. 2021, 22, 52. [Google Scholar] [CrossRef]
  41. Emes, R.D.; Grant, S.G.N. Evolution of Synapse Complexity and Diversity. Annu. Rev. Neurosci. 2012, 35, 111–131. [Google Scholar] [CrossRef] [PubMed]
  42. Gobert, D.; Topolnik, L.; Azzi, M.; Huang, L.; Badeaux, F.; Desgroseillers, L.; Sossin, W.S.; Lacaille, J.-C. Forskolin induction of late-LTP and up-regulation of 5′ TOP mRNAs translation via mTOR, ERK, and PI3K in hippocampal pyramidal cells. J. Neurochem. 2008, 106, 1160–1174. [Google Scholar] [CrossRef] [PubMed]
  43. Otmakhov, N.; Khibnik, L.; Otmakhova, N.; Carpenter, S.; Riahi, S.; Asrican, B.; Lisman, J. Forskolin-Induced LTP in the CA1 Hippocampal Region Is NMDA Receptor Dependent. J. Neurophysiol. 2004, 91, 1955–1962. [Google Scholar] [CrossRef]
  44. Bigdeli, T.B.; Genovese, G.; Georgakopoulos, P.; Meyers, J.L.; Peterson, R.E.; Iyegbe, C.O.; Medeiros, H.; Valderrama, J.; Achtyes, E.D.; Kotov, R.; et al. Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry. Mol. Psychiatry 2020, 25, 2455–2467. [Google Scholar] [CrossRef]
  45. Wuttke, M.; Li, Y.; Li, M.; Sieber, K.B.; Feitosa, M.F.; Gorski, M.; Tin, A.; Wang, L.; Chu, A.Y.; Hoppmann, A.; et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat. Genet. 2020, 51, 957–972. [Google Scholar] [CrossRef]
  46. Rees, E.; Han, J.; Morgan, J.; Carrera, N.; Escott-Price, V.; Pocklington, A.J.; Duffield, M.; Hall, L.S.; Legge, S.E.; Pardiñas, F.; et al. De novo mutations identified by exome sequencing implicate rare missense variants in SLC6A1 in schizophrenia. Nat. Neurosci. 2020, 23, 179–184. [Google Scholar] [CrossRef] [PubMed]
  47. Ambalavanan, A.; Girard, S.L.; Ahn, K.; Zhou, S.; Dionne-Laporte, A.; Spiegelman, D.; Bourassa, C.V.; Gauthier, J.; Hamdan, F.F.; Xiong, L.; et al. De novo variants in sporadic cases of childhood onset schizophrenia. Eur. J. Hum. Genet. 2016, 24, 944–948. [Google Scholar] [CrossRef]
  48. Girard, S.L.; Gauthier, J.; Noreau, A.; Xiong, L.; Zhou, S.; Jouan, L.; Dionne-Laporte, A.; Spiegelman, D.; Henrion, E.; Diallo, O.; et al. Increased exonic de novo mutation rate in individuals with schizophrenia. Nat. Genet. 2011, 43, 860–863. [Google Scholar] [CrossRef]
  49. Guipponi, M.; Santoni, F.; Setola, V.; Gehrig, C.; Rotharmel, M.; Cuenca, M.; Guillin, O.; Dikeos, D.; Georgantopoulos, G.; Papadimitriou, G.; et al. Exome Sequencing in 53 Sporadic Cases of Schizophrenia Identifies 18 Putative Candidate Genes. PLoS ONE 2014, 9, e112745. [Google Scholar] [CrossRef]
  50. Gulsuner, S.; Walsh, T.; Watts, A.C.; Lee, M.K.; Thornton, A.M.; Casadei, S.; Rippey, C.; Shahin, H.; Nimgaonkar, V.L.; Go, R.C.; et al. Spatial and Temporal Mapping of De Novo Mutations in Schizophrenia to a Fetal Prefrontal Cortical Network. Cell 2013, 154, 518–529. [Google Scholar] [CrossRef]
  51. E McCarthy, S.; Gillis, J.; Kramer, M.; Lihm, J.; Yoon, S.; Berstein, Y.; Mistry, M.; Pavlidis, P.; Solomon, R.; Ghiban, E.; et al. De novo mutations in schizophrenia implicate chromatin remodeling and support a genetic overlap with autism and intellectual disability. Mol. Psychiatry 2014, 19, 652–658. [Google Scholar] [CrossRef]
  52. Takata, A.; Xu, B.; Ionita-Laza, I.; Roos, J.L.; Gogos, J.A.; Karayiorgou, M. Loss-of-Function Variants in Schizophrenia Risk and SETD1A as a Candidate Susceptibility Gene. Neuron 2014, 82, 773–780. [Google Scholar] [CrossRef]
  53. Wang, Q.; Li, M.; Yang, Z.; Hu, X.; Wu, H.-M.; Ni, P.; Ren, H.; Deng, W.; Li, M.; Ma, X.; et al. Increased co-expression of genes harboring the damaging de novo mutations in Chinese schizophrenic patients during prenatal development. Sci. Rep. 2015, 5, 18209. [Google Scholar] [CrossRef]
  54. Xu, B.; Ionita-Laza, I.; Roos, J.L.; Boone, B.; Woodrick, S.; Sun, Y.; Levy, S.; Gogos, J.A.; Karayiorgou, M. De novo gene mutations highlight patterns of genetic and neural complexity in schizophrenia. Nat. Genet. 2012, 44, 1365–1369. [Google Scholar] [CrossRef]
  55. Samocha, K.E.; Kosmicki, J.A.; Karczewski, K.J.; O’Donnell-Luria, A.H.; Pierce-Hoffman, E.; MacArthur, D.G.; Neale, B.M.; Daly, M.J. Regional missense constraint improves variant deleteriousness prediction. BioRxiv 2017. Preprint. [Google Scholar] [CrossRef]
  56. International Schizophrenia Consortium. Rare chromosomal deletions and duplications increase risk of schizophrenia. Nature 2008, 455, 237–241. [Google Scholar] [CrossRef]
  57. Levinson, D.F.; Duan, J.; Oh, S.; Wang, K.; Sanders, A.R.; Shi, J.; Zhang, N.; Mowry, B.J.; Olincy, A.; Amin, F.; et al. Copy Number Variants in Schizophrenia: Confirmation of Five Previous Findings and New Evidence for 3q29 Microdeletions and VIPR2 Duplications. Am. J. Psychiatry 2011, 168, 302–316. [Google Scholar] [CrossRef] [PubMed]
  58. Samocha, K.E.; Robinson, E.B.; Sanders, S.J.; Stevens, C.; Sabo, A.; McGrath, L.M.; Kosmicki, J.A.; Rehnström, K.; Mallick, S.; Kirby, A.; et al. A framework for the interpretation of de novo mutation in human disease. Nat. Genet. 2014, 46, 944–950. [Google Scholar] [CrossRef]
Ijms 25 00946 g001
Figure 1. Gene set enrichment analysis of LTP datasets for common variant association with schizophrenia, bipolar disorder, and chronic kidney disease (negative control). Gene sets represent differentially expressed genes at 30, 60 and 120 min following LTP compared to baseline expression in the hippocampus. p-values have been log10 transformed. Dotted lines represent the threshold for statistical significance after correction for multiple testing.
Figure 1. Gene set enrichment analysis of LTP datasets for common variant association with schizophrenia, bipolar disorder, and chronic kidney disease (negative control). Gene sets represent differentially expressed genes at 30, 60 and 120 min following LTP compared to baseline expression in the hippocampus. p-values have been log10 transformed. Dotted lines represent the threshold for statistical significance after correction for multiple testing.
Ijms 25 00946 g001
Ijms 25 00946 g002
Figure 2. Enrichment of LTP gene sets in schizophrenia-associated CNVs. Logistic regression analysis was used to test the relative enrichment of each gene set in case CNVs compared to control CNVs. Deletions and duplicates were analysed together (All CNVs) and separately. p-values resulting from logistic regression analysis have undergone permutation correction to derive empirical p-values and adjust for inflation and background enrichment. Dotted lines represent the threshold for significance following permutation correction.
Figure 2. Enrichment of LTP gene sets in schizophrenia-associated CNVs. Logistic regression analysis was used to test the relative enrichment of each gene set in case CNVs compared to control CNVs. Deletions and duplicates were analysed together (All CNVs) and separately. p-values resulting from logistic regression analysis have undergone permutation correction to derive empirical p-values and adjust for inflation and background enrichment. Dotted lines represent the threshold for significance following permutation correction.
Ijms 25 00946 g002
Table 1. Enrichment for common variant association with bipolar disorder, split by bipolar disorder subtypes (I and II). Results obtained from gene set enrichment analyses in MAGMA, conditional on all tissue-expressed genes. Bold text indicates significant enrichment of gene set. p-values were Bonferroni adjusted to account for testing of multiple gene sets. SE = Standard error.
Table 1. Enrichment for common variant association with bipolar disorder, split by bipolar disorder subtypes (I and II). Results obtained from gene set enrichment analyses in MAGMA, conditional on all tissue-expressed genes. Bold text indicates significant enrichment of gene set. p-values were Bonferroni adjusted to account for testing of multiple gene sets. SE = Standard error.
Bipolar Disorder Subtypeβ Value (SE)p-ValueAdjusted p-Value
30 min TRAP-seqBipolar I−0.090 (0.13)0.761.0
Bipolar II−0.090 (0.11)0.801.0
30 min total RNA-seqBipolar I−0.066 (0.22)0.621.0
Bipolar II−0.11 (0.18)0.721.0
60 min TRAP-seqBipolar I0.13 (0.029)3.1 × 10−61.9 × 10−5
Bipolar II0.0025 (0.025)0.461.0
60 min total RNA-seqBipolar I−0.027 (0.10)0.600.898
Bipolar II−0.096 (0.091)0.851.0
120 min TRAP-seqBipolar I0.094 (0.034)0.00270.016
Bipolar II0.0070 (0.029)0.401.0
120 min total RNA-seqBipolar I−0.028 (0.050)0.711.0
Bipolar II−0.028 (0.043)0.741.0
Table 2. LTP gene sets were tested for enrichment of de novo protein truncating variants (PTVs) or missense variants with an MPC score of greater than 2 from schizophrenia cases. Adjusted p-values have been corrected for background enrichment in tissue-expressed genes and multiple comparisons using the Bonferroni method.
Table 2. LTP gene sets were tested for enrichment of de novo protein truncating variants (PTVs) or missense variants with an MPC score of greater than 2 from schizophrenia cases. Adjusted p-values have been corrected for background enrichment in tissue-expressed genes and multiple comparisons using the Bonferroni method.
LTP Gene SetSchizophrenia Ultra-Rare Coding Variant Enrichment
Observed/Expectedp-ValueAdjusted p-ValueRate Ratio (95% Confidence Intervals)
30-min total2/0.4970.1160.7363.45 (0.416–12.5)
30-min TRAP5/3.450.6101.001.24 (0.401–2.91)
60-min total7/3.270.1140.6421.84 (0.737–3.82)
60-min TRAP80/57.40.1010.9311.23 (0.954–1.56)
120-min total17/14.20.9001.001.02 (0.592–1.66)
120-min TRAP80/51.40.008920.1811.39 (1.08–1.77)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wellard, N.L.; Clifton, N.E.; Rees, E.; Thomas, K.L.; Hall, J. The Association of Hippocampal Long-Term Potentiation-Induced Gene Expression with Genetic Risk for Psychosis. Int. J. Mol. Sci. 2024, 25, 946. https://doi.org/10.3390/ijms25020946

AMA Style

Wellard NL, Clifton NE, Rees E, Thomas KL, Hall J. The Association of Hippocampal Long-Term Potentiation-Induced Gene Expression with Genetic Risk for Psychosis. International Journal of Molecular Sciences. 2024; 25(2):946. https://doi.org/10.3390/ijms25020946

Chicago/Turabian Style

Wellard, Natalie L., Nicholas E. Clifton, Elliott Rees, Kerrie L. Thomas, and Jeremy Hall. 2024. "The Association of Hippocampal Long-Term Potentiation-Induced Gene Expression with Genetic Risk for Psychosis" International Journal of Molecular Sciences 25, no. 2: 946. https://doi.org/10.3390/ijms25020946

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop