Next Article in Journal
The Preservation of PPARγ Genome Duplicates in Some Teleost Lineages: Insights into Lipid Metabolism and Xenobiotic Exploitation
Next Article in Special Issue
Toxoplasma gondii Seropositivity Interacts with Catechol-O-methyltransferase Val105/158Met Variation Increasing the Risk of Schizophrenia
Previous Article in Journal
Comparative Transcriptomic Analyses of Different Jujube Cultivars Reveal the Co-Regulation of Multiple Pathways during Fruit Cracking
Previous Article in Special Issue
Exploring the Contribution to ADHD of Genes Involved in Mendelian Disorders Presenting with Hyperactivity and/or Inattention
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 1
10.3390/genes13010106
genes-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

Adolescent Verbal Memory as a Psychosis Endophenotype: A Genome-Wide Association Study in an Ancestrally Diverse Sample

by
Baihan Wang
1,*,
Olga Giannakopoulou
1,2,
Isabelle Austin-Zimmerman
1,3,
Haritz Irizar
1,4,
Jasmine Harju-Seppänen
1,
Eirini Zartaloudi
1,
Anjali Bhat
1,
Andrew McQuillin
1,
Karoline Kuchenbäcker
1,2 and
Elvira Bramon
1,*
1
Department of Mental Health Neuroscience, Division of Psychiatry, University College London, London W1T 7BN, UK
2
UCL Genetics Institute, Division of Biosciences, University College London, London WC1E 6BT, UK
3
Social, Genetic and Developmental Psychiatry Centre, Institute of Psychiatry, Psychology and Neuroscience, King’s College London, London SE5 8AF, UK
4
Department of Genetics & Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA
*
Authors to whom correspondence should be addressed.
Genes 2022, 13(1), 106; https://doi.org/10.3390/genes13010106
Submission received: 26 November 2021 / Revised: 16 December 2021 / Accepted: 29 December 2021 / Published: 3 January 2022
(This article belongs to the Special Issue Advances in Genetics of Psychiatric Disorders)
Browse Figure
Review Reports Versions Notes

Abstract

:
Verbal memory impairment is one of the most prominent cognitive deficits in psychosis. However, few studies have investigated the genetic basis of verbal memory in a neurodevelopmental context, and most genome-wide association studies (GWASs) have been conducted in European-ancestry populations. We conducted a GWAS on verbal memory in a maximum of 11,017 participants aged 8.9 to 11.1 years in the Adolescent Brain Cognitive Development Study®, recruited from a diverse population in the United States. Verbal memory was assessed by the Rey Auditory Verbal Learning Test, which included three measures of verbal memory: immediate recall, short-delay recall, and long-delay recall. We adopted a mixed-model approach to perform a joint GWAS of all participants, adjusting for ancestral background and familial relatedness. The inclusion of participants from all ancestries increased the power of the GWAS. Two novel genome-wide significant associations were found for short-delay and long-delay recall verbal memory. In particular, one locus (rs9896243) associated with long-delay recall was mapped to the NSF (N-Ethylmaleimide Sensitive Factor, Vesicle Fusing ATPase) gene, indicating the role of membrane fusion in adolescent verbal memory. Based on the GWAS in the European subset, we estimated the SNP-heritability to be 15% to 29% for the three verbal memory traits. We found that verbal memory was genetically correlated with schizophrenia, providing further evidence supporting verbal memory as an endophenotype for psychosis.

1. Introduction

Psychotic disorders, including schizophrenia and bipolar disorder, are highly heritable (heritability estimated to be approximately 80%) [1,2]. Although the latest genome-wide association studies (GWASs) have discovered 270 genomic loci associated with schizophrenia and 64 associated with bipolar disorder [3,4], the mechanisms through which genetic variation in those loci contributes to disease risk are unknown. One potential solution to this issue is by investigating the endophenotypes of psychosis, which are biomarkers of the disease that lie closer to the genetic effects and are thought to be part of the mechanistic path from genetic variation to clinical manifestation [5,6]. Endophenotypes can offer insights into the biological mechanisms that are most pertinent to psychosis, thus informing our understanding of its complex aetiology.
Cognitive deficits of different severity are common in psychosis, and often manifest before symptom onset [7,8]. Therefore, psychosis has been suggested to be an outcome of neurodevelopmental abnormalities that occur in early life [9,10,11]. Among all cognitive domains, verbal memory impairment is one of the best characterised in psychosis. Previous meta-analyses have found prominent deficits in verbal memory in patients with schizophrenia (Cohen’s d = −1.20 to −0.85 [12]) and bipolar disorder (Cohen’s d = −0.56 to −0.50 [13]). Milder deficits have also been reported by meta-analyses on the relatives of patients with schizophrenia (Cohen’s d = −0.54 to −0.44 [14,15]) and bipolar disorder (Cohen’s d = −0.42 to −0.33 [16,17]). Those findings are consistent with verbal memory being associated with the genetic risk of psychosis, making it a potential endophenotype. Additionally, there is evidence that such deficits in verbal memory is already present in childhood and adolescence for those with increased genetic risk for psychosis [18,19]. Therefore, verbal memory impairment may reflect underlying neurodevelopmental abnormalities that convey the genetic risk of psychosis before its symptom onset or formal diagnosis.
Previous GWASs have found over 200 loci associated with general cognitive ability, as well as its negative genetic correlation with schizophrenia [20,21,22]. Chromosomal aberrations, including copy number variants, can also influence cognition [23,24,25]. Nevertheless, the population correlation between verbal memory and general cognitive ability at the phenotypic level was estimated to be 0.24 to 0.39, indicating that they are overlapping but different constructs [26]. Investigating the genetic basis of specific cognitive domains, such as verbal memory, could pinpoint the biological mechanisms that are most relevant to the aetiology of psychosis. However, there have been only a few GWASs conducted on verbal memory. The earliest findings came from Papassotiropoulos et al., who conducted GWASs in two samples from two Swiss cohorts (351 and 1073 participants) and found that delayed verbal episodic memory was associated with alleles in the KIBRA and CTNNBL1 genes [27,28]. Debette et al. conducted a GWAS meta-analysis in 29,076 older adults, and found that verbal declarative memory was associated with APOE, which has also been identified in GWASs on Alzheimer’s disease [29]. The association of verbal memory with APOE variants was replicated in a later GWAS by Arpawong et al. in 7486 older adults, with two additional significant SNPs identified in TOMM40 [30]. With regards to schizophrenia, in a sample of 127 patients and 136 controls, Nakahara et al. found that verbal memory was associated with a SNP near NDUFS4 [31]. Greenwood et al. conducted a GWAS on 11 endophenotypes for schizophrenia in a sample of 1533 participants, but found no genome-wide significant associations for verbal memory [32]. However, since the two largest GWASs on verbal memory were both conducted in older populations, their findings might be more relevant to neurodegenerative disorders such as Alzheimer’s disease instead of psychosis. Little is known about the biological basis of verbal memory in a neurodevelopmental context, and no GWASs on verbal memory have been conducted in children or young people so far.
The lack of diversity has long been an issue in GWASs, as 78% of published GWASs until 2019 were conducted only in European individuals [33]. Unsurprisingly, most previous GWASs on verbal memory have also been conducted in homogeneous European samples. This is mainly due to the lack of diversity in existing databases, as well as the statistical challenges of analysing samples from multiple ancestral backgrounds [34,35]. Such lack of diversity means we have a limited understanding of disease aetiology in underrepresented populations. By recruiting and including participants of all ancestries, we can maximise the power of GWASs, and may identify important genetic variants that are too rare to study in European populations [36]. One solution to these statistical difficulties is a meta-analysis of results of GWASs in multiple ancestries. However, this approach requires each of the primary GWASs in the meta-analysis to be conducted in a relatively homogeneous sample, and individuals of mixed ancestry are usually excluded from the analysis. Moreover, the sample size of each homogeneous ancestral group might be too small to estimate the small effects of individual variants. For heterogeneous multi-ancestry samples, joint mixed-model methods might be the optimal approach [34]. In the mixed models, modelling a genetic relationship matrix as a random effect allows the inclusion of individuals of mixed ancestry. This approach is not limited by the sample size of each homogeneous ancestral group and maximises statistical power [34,37,38]. One example of this is the linear mixed models with orthogonally partitioned structure included in the GENetic EStimation and Inference in Structured samples (GENESIS) R/Bioconductor package [37,38]. The GENESIS package partitions genetic structure into distant ancestral background and recent familial relatedness, thus accounting for ancestry better than other traditional mixed-model methods [37,38].
All in all, previous GWASs on verbal memory have mostly been conducted in European older adults, making them less relevant to psychosis as a neurodevelopmental condition. Therefore, the current study aims to investigate the genetic basis of verbal memory in adolescents as a neurodevelopmental endophenotype for psychosis. We conducted a GWAS on three verbal memory traits measured by the Adolescent Brain Cognitive Development (ABCD) Study® in over 10,000 adolescents of diverse ancestries. We adopted a novel mixed-model approach using the GENESIS package to account for the diverse ancestries in the sample. We also estimated the heritability of the three verbal memory traits and their genetic correlations with other psychiatric traits, in order to investigate the shared genetic basis between verbal memory and psychosis.

2. Materials and Methods

2.1. Participants

The ABCD Study® is a longitudinal study that recruited 11,878 participants aged about 9–11 years at baseline from 21 sites across the United States [39,40,41]. The 21 sites were nationally distributed and participants were recruited from schools in nearby catchment areas. Some schools were oversampled in order to obtain a sample representative of the national population, in terms of gender, ethnicity, socioeconomic status, and urbanicity [39]. Interested families were invited to the research sites and all children underwent a comprehensive set of demographic, physical, cognitive, mental health, and neuroimaging assessments, along with their blood and saliva samples collected for genotyping. Participants were excluded from the ABCD Study® if they were not fluent in English, had a history of traumatic brain injury, or had a current diagnosis of schizophrenia, moderate/severe autism spectrum disorder, intellectual disability, or alcohol/substance use disorders [42,43]. For this GWAS, we used the phenotypic and genetic data at baseline from ABCD Data Release 3.0, which have now been anonymised and released to bona fide researchers after a registration process. Details on the recruitment procedures of the ABCD Study® can be found in a previous publication [39].

2.2. Verbal Memory Assessment

The Rey Auditory Verbal Learning Test (RAVLT) is a verbal memory test widely used in research and clinical practice [44]. In the ABCD Study®, participants were firstly verbally presented with a list of 15 words (List A) and immediately asked to recall as many words from the list as possible. The same procedure was repeated four times for the same list (List A) after the first trial. The total number of words correctly recalled in the five trials was recorded as the RAVLT immediate recall score. After learning List A, participants were presented with a different list (List B, the interference list) and asked to recall words from List B. Right after recalling List B, participants were asked to recall List A again, and the correct number of words recalled in List A was recorded as the RAVLT short-delay recall score. Participants were assessed by other cognitive tests in the following 30 min, and were asked to recall list A at the end of the assessment. The final number of words correctly recalled from List A was recorded as the RAVLT long-delay recall score.

2.3. Genotyping, Imputation, and Quality Control

Participants were genotyped by the Rutgers University Cell and DNA Repository (RUCDR) using the Affymetrix NIDA SmokeScreen Array [45]. All genotyped data underwent standard quality control following the RICOPILI pipeline [46]. They were checked against the Haplotype Reference Consortium and the 1000 Genomes Project for consistencies [47,48], and then uploaded to the TOPMed server for imputation using mixed ancestry and Eagle v2.4 phasing [49,50]. We downloaded the imputed data from the ABCD Data Repository and used bcftools to annotate them with rs IDs based on dbSNP153 [51]. Imputed SNPs with dosage levels were converted to best-guess genotype format using PLINK v2.0 with a hard call threshold of 0.1 [52]. We performed further post-imputation quality control by removing SNPs with imputation quality score (r2) < 0.3 or minor allele frequency (MAF) < 0.01. The final sample included 11,017 participants and 11,229,083 variants after quality control. Details of genotyping, imputation, and quality control can be found in the Supplementary Methods.

2.4. Relationship Inference and Principal Component Analysis

Since the ABCD Study® included an ancestrally diverse sample with a combination of singletons, siblings, and twins, we used the GENESIS package in R-4.0.3 to conduct a linear mixed-model GWAS [38,53]. Based on the genotyped data after quality control, we firstly used KING-robust 2.2.5 to infer the unadjusted kinship coefficient for each pair of participants in the dataset [54]. We then used the PC-AiR package [38] and the output from KING-robust to conduct a principal component analysis of ancestral background that accounts for familial relatedness. This was performed on the genotyped data after pruning using the SNPRelate package in R [55]. Finally, to get a more accurate estimation of familial relatedness, we performed a PC-Relate analysis to generate a new kinship matrix with kinship coefficients adjusted for ancestral background [56]. To perform an ancestry-specific analysis in comparison to the multi-ancestry analysis, we also identified a subset of participants who were ancestrally European based on the principal component analysis. The same procedure was repeated to estimate the principal components and familial relatedness within the European subset. However, we were unable to analyse participants of other ancestries due to their small sample sizes. Details of relationship inference and principal component analysis can be found in the Supplementary Methods.

2.5. Genome-Wide Association Testing

A mixed-model genome-wide association test was conducted using the GENESIS package for the overall sample and the European subset [37]. The RAVLT immediate recall score, short-delay recall score, or long-delay recall score was added as an outcome. Age, sex, the first eight (for the overall sample) or four (for the European subset) ancestrally representative principal components were included as fixed effects, and the ancestry-adjusted kinship matrix was included as a random effect. The number of principal components was chosen based on visual inspection of the principal component analysis plots. Principal components were included in the model if they showed variability across ancestries. The genome-wide significant p-value threshold was set to 5 × 10−8. Additionally, we checked if the lead SNPs in loci that reached genome-wide significance for one verbal memory trait also reached the nominal significance level for the other two traits.

2.6. Locus Definition and Gene Mapping

Summary statistics of the GWAS results of the three verbal memory traits in the whole sample were uploaded to the platform of Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA) [57]. A reference panel from all populations in the 1000 Genomes Project (Phase 3) was used to obtain r2 [47]. In FUMA, independent significant SNPs were firstly defined as SNPs with p < 5 × 10−8 independent from each other at r2 < 0.6. Lead SNPs of genomic loci were then defined as a subset of independent significant SNPs at r2 < 0.1. Lead SNPs with a distance closer than 250 kb were merged to one locus. SNPs with r2≥ 0.6 were considered to be in linkage disequilibrium (LD) with the lead SNPs and thus their good proxies.
We used three strategies to map the lead SNPs and their good proxies to protein-coding genes. (1) Positional mapping: SNPs were mapped to genes if the SNPs are located within a 10 kb window around the gene boundaries. (2) Expression quantitative trait loci (eQTL) mapping: SNPs were mapped to genes if the SNPs have a significant impact on the gene expression. eQTLs were selected from GTEx v8 brain tissues and only significant SNP-gene pairs in GTEx (FDR < 0.05) were used [58]. Mapping was based on cis-eQTLs with a 1Mb window around each gene. (3) Chromatin interaction mapping: SNPs were mapped to genes if there are significant chromatin interactions between the SNPs and nearby or distant genes. Chromatin interaction information of four tissues was obtained from Hi-C, including adult cortex, fetal cortex, dorsolateral prefrontal cortex, and hippocampus [59,60]. Only significant interactions were used based on an FDR cutoff of 10−6. Promoter regions were defined as 250 bp upstream and 500 bp downstream from the transcriptional start site. SNPs were filtered by enhancers and genes were filtered by promoters based on annotations from the Roadmap Epigenomics Project brain tissue data [61].

2.7. Heritability and Genetic Correlations

Based on the GWAS results of the European subset in the ABCD Study®, we used linkage disequilibrium score regression (LDSC) to estimate the SNP heritability (h2; phenotypic variance explained by common SNPs) of the three verbal memory traits. We also used LDSC to test their genetic correlations (rg) with each other, as well as with three psychiatric traits (schizophrenia, bipolar disorder, and major depressive disorder) and educational attainment. LD scores calculated from the European participants in the 1000 Genomes Project (Phase 3) were used [47]. We used the latest GWAS summary statistics of schizophrenia, bipolar disorder, and major depressive disorder from the Psychiatric Genomics Consortium (PGC; https://www.med.unc.edu/pgc/; accessed on 2 November 2021) [3,4,62]. The GWAS summary statistics of major depressive disorder excluded the 23andMe sample due to the data transfer agreement between the PGC and 23andMe. The GWAS summary statistics of educational attainment were retrieved from the Social Science Genetic Association Consortium (https://www.thessgac.org/; accessed on 5 October 2021) [63]. All GWAS summary statistics included only participants of European ancestries.

3. Results

3.1. Sample Characteristics

As shown in Table 1, among the 11,017 participants who passed genetic quality control in the ABCD Study®, there were 5176 females (47%) and 5805 males (53%). The mean age of the participants was 9.9 years (SD = 0.6 years, range = 8.9 to 11.1 years). Based on broadly defined ethnicity reported by parents, 5872 (53%) participants were reported as White, 2099 (19%) as Hispanic, 1684 (14%) as Black, and 1% as Asian. The remaining 11% included children of American Indian/Alaska Native, Native Hawaiian and other Pacific Islander, and mixed ethnicity. Participants scored on average 44.2 for RAVLT immediate recall (SD = 9.9), 9.7 for RAVLT short-delay recall (SD = 3.0), and 9.2 for RAVLT long-delay recall (SD = 3.2). The three verbal memory traits were highly correlated at the phenotype level: r = 0.72 between RAVLT immediate recall and short-delay recall (p < 0.001); r = 0.73 between RAVLT immediate recall and long-delay recall (p < 0.001); r = 0.79 between RAVLT short-delay recall and long-delay recall (p < 0.001). All three verbal memory traits were approximately normally distributed in the sample (Figure S1).

3.2. Principal Component Analysis and Relatedness Estimation

32 principal components were returned by PC-AiR for the whole dataset. The distribution of the samples along the first eight principal components coloured by reported ethnicity can be seen in Figure S2. Based on the plots of the principal components and participants’ ethnicity reported by their parents, we identified a subset of 5763 participants of European ancestry (Supplementary Methods). In the whole sample, PC-Relate revealed 416 pairs of monozygotic twins (kinship coefficient > 0.354) and 1356 pairs of dizygotic twins or siblings (0.176 < kinship coefficient ≤ 0.354; Figure S3).

3.3. Genome-Wide Association Testing

The GWAS in the European subset included 5635 participants for RAVLT immediate recall, 5667 participants for RAVLT short-delay recall, and 5640 participants for RAVLT long-delay recall. No genome-wide significant associations were found for any of the three verbal memory traits (Figure S4). In contrast, the inclusion of non-European participants increased the sample sizes to 10,726 participants for RAVLT immediate recall, 10,801 participants for RAVLT short-delay recall, and 10,757 participants for RAVLT long-delay recall. The GWAS results in the whole sample are shown in Figure 1. No genome-wide significant SNPs were identified for RAVLT immediate recall. For RAVLT short-delay recall, one SNP (rs73984566) in an intergenic region of chromosome 2 reached genome-wide significance (beta = −0.847, standard error (SE) = 0.151, p = 1.86 × 10−8). For RAVLT long-delay recall, we found one genome-wide significant SNP (rs9896243) located in the intron of the NSF gene in chromosome 17 (beta = −0.309, SE = 0.055, p = 2.22 × 10−8). Table 2 shows the detailed information about each genome-wide significant locus. Both genome-wide significant SNPs associated with one verbal memory trait also reached the nominal significance threshold for the other two verbal memory traits, with an effect size in the same direction (Table S1).

3.4. Gene Mapping

For the gene mapping we considered the lead SNPs in the two novel loci, as well as their good proxies. In the locus associated with RAVLT short-delay recall (lead SNP: rs73984566), there was strong evidence for two genes, NABP1 and TMEFF2, based on chromatin interaction data in fetal cortex. The lead SNP (rs9896243) in the locus associated with RAVLT long-delay recall is located in the intron of the NSF gene, which was supported by positional mapping of the locus in FUMA. Additionally, a good proxy of the lead SNP (rs199533; r2 = 0.61) results in a synonymous amino acid change in this gene according to HaploReg v4.1 [64]. This SNP (rs199533) was genotyped in the study and showed an association with RAVLT long-delay recall at the nominal significance level (beta = −0.285, p = 1.76 × 10−6). Other variants in the locus have also been shown to affect the expression of the same gene in brain tissues. Additionally, the locus was mapped to seven other genes in the locus (17q21.31) based on brain eQTLs. Details of all mapped genes are summarised in Table S2.

3.5. Heritability and Genetic Correlations

Based on the GWAS in the European subset, LDSC estimated the SNP heritability (h2) to be 0.15 (SE = 0.07) for RAVLT immediate recall, 0.29 (SE = 0.07) for RAVLT short-delay recall, and 0.21 (SE = 0.07) for RAVLT long-delay recall. The genetic correlations among the three verbal memory traits were very high: rg = 0.99 (SE = 0.12) between RAVLT immediate recall and short-delay recall (p < 0.001); rg = 0.91 (SE = 0.12) between RAVLT immediate recall and long-delay recall (p < 0.001); rg = 0.93 (SE = 0.06) between RAVLT short-delay recall and long-delay recall (p < 0.001). The results of their genetic correlations with other psychiatric traits are shown in Table 3. All three verbal memory traits were significantly and positively correlated with educational attainment (rg = 0.31–0.53, p ≤ 0.001). Schizophrenia was found to be significantly and negatively correlated with both RAVLT immediate recall (rg = −0.29, SE = 0.10, p = 0.003) and RAVLT short-delay recall (rg = −0.15, SE = 0.06, p = 0.017), as well as RAVLT long-delay recall at a trend level (rg = −0.15, SE = 0.08, p = 0.063). We also found a significant negative genetic correlation between RAVLT short-delay recall and major depressive disorder (rg = −0.14, SE = 0.07, p = 0.050).

4. Discussion

This is one of the first studies investigating the genetic basis of verbal memory in adolescents. It benefited greatly from the diversity within the ABCD Study®, which allowed our GWAS to identify two genome-wide significant associations for RAVLT short-delay recall and long-delay recall. Genetic correlations estimated by LDSC based on the European GWAS revealed negative genetic correlations between schizophrenia and the three verbal memory traits, providing further evidence that supports verbal memory as a psychosis endophenotype.
Our study included all participants in the ABCD Study®, which increased the representation of diverse populations and the power to detect genome-wide significant associations [65]. An alternative way to perform GWASs in diverse populations is to stratify participants into different ancestral groups, conduct the GWASs separately, and meta-analyse the results together. However, because of the large proportion of admixed individuals in the ABCD Study®, such stratification could cause substantial sample loss. Instead, we employed a mixed-model approach using the GENESIS package to analyse the whole sample altogether while adjusting for ancestral background and familial relatedness. We successfully distinguished individual relatedness due to family structure from their ancestral backgrounds. This allowed for better adjustment for ancestry, especially for SNPs with relatively small or large differences in allele frequency across ancestries, which are usually not fully controlled for by other mixed-model methods [37]. The inclusion of all participants also boosted the power of the GWAS and it is only with this larger sample we identified two genome-wide significant associations. The associations were not identified in the European subset, most probably due to having just half the sample size and the low allele frequency of rs73984566 in European populations (MAF < 0.001).
We found one locus located in chromosome 2 (lead SNP: rs73984566) associated with RAVLT short-delay recall. Chromatin interaction data provided evidence for two genes associated with this locus, NABP1 and TMEFF2. NABP1 (Nucleic Acid Binding Protein 1) is involved in many ubiquitous DNA metabolic processes [66], but its role in brain development remains unclear and needs further research. More information is available for the TMEFF2 (transmembrane protein with EGF-like and two follistatin-like domains 2) gene, which is widely expressed in the brain and has been proposed to be a survival factor for neurons in the hippocampus and midbrain [67]. The neuroprotective effect of TMEFF2 has also been found in Alzheimer’s disease as it binds Amyloid-β oligomer and Amyloid-β protein precursor [68]. One study found that methylation of the TMEFF2 promoter region was associated with poor outcomes in patients with glioma, a type of brain tumour [69]. However, few studies have investigated the role of TMEFF2 in a neurodevelopmental context, and more research is needed to clarify its role in adolescent verbal memory.
We also found RAVLT short-delay recall to be associated with the 17q21.31 locus. The lead SNP, rs9896243, is located in the intronic region of the NSF (N-ethylmaleimide sensitive factor, vesicle fusing ATPase) gene. eQTL data provided further evidence for the role of this gene. NSF is an important factor in synaptic neurotransmission that facilitates the recycling of SNARE (soluble N-ethylmaleimide sensitive factor attachment protein receptors) proteins, which mediate membrane fusion [70]. Previous GWASs have reported many associations between the NSF gene and brain volume or structure measures [71,72], Parkinson’s disease [73], and neuroticism [74]. Notably, the NSF gene has also been identified by previous GWASs on general cognitive ability [20,21], which is in line with the phenotypic correlation observed between general cognitive ability and verbal memory [26]. This suggests that the heritability of general cognitive ability could be further partitioned into specific cognitive domains, although unique genetic influences on specific cognitive domains may also exist. Moreover, the same SNP (rs9896243) was found to be associated with worry measurement at the genome-wide significance level in a GWAS on neuroticism [74], which is a risk factor for various psychiatric conditions including psychosis. Indeed, increased SNARE protein–protein interactions have also been found in postmortem brain samples of patients with schizophrenia, supporting the hypothesis of general synaptic dysfunction in schizophrenia [75]. Since adolescent verbal memory is an endophenotype for psychosis [14,15], our finding suggests that the NSF might play a role in the synaptic dysfunction in psychosis by modulating SNARE complex activity during neurodevelopment. The 17q21.31 locus is also associated with Koolen-de Vries syndrome [76,77,78], a disorder characterised by developmental delay and intellectual disability, thus providing further evidence for the involvement of this locus in neurodevelopment.
Based on the results of the GWAS in the European subset, we estimated the SNP heritability to be 15%, 29%, and 21% for RAVLT immediate recall, short-delay recall, and long-delay recall. This supports the influence of common variants on verbal memory in adolescents. A previous family study using a similar measure (the California Verbal Learning Test) estimated the heritability to be 35%, 50%, and 73% for the three verbal memory traits, respectively, [79]. The “missing heritability” estimated from the current GWAS is expected, given that more SNPs with small effects are yet to be found, as well as other factors that GWASs could not account for (e.g., rare variants) [80]. Moreover, heritability estimates could also be influenced by age, as the sample in Bertisch et al. included adults while the ABCD Study® was only composed of children and adolescents [79]. There is evidence that the heritability of general cognitive ability increases during adolescence due to children interacting with their environments based on genetic propensities [14,15,16,17], although it is unclear if this pattern also exists for specific cognitive domains. Nevertheless, investigating the genetic basis of verbal memory during adolescence is still important, as it can help to identify relevant neurodevelopmental processes that manifest prior to the onset of many psychiatric conditions.
We found significant and positive genetic correlations between all three verbal memory traits and educational attainment, which is expected given the close relationship between the two phenotypes. We also found negative genetic correlations between all three verbal memory traits and schizophrenia, although the association was trend-level for RAVLT long-delay recall. However, as the genetic correlations within the three traits were very high and only European participants were included in this analysis, we believe the trend-level association will reach statistical significance with bigger sample sizes. This finding is in line with previous studies that found deficits in verbal memory amongst the relatives of patients with psychosis, supporting the role of verbal memory as a psychosis endophenotype [14,15,16,17]. Moreover, the negative genetic correlation between RAVLT short-delay recall and depression indicates that verbal memory deficits might not be specific to psychosis. Given the shared heritability within psychiatric disorders [81], verbal memory deficits in adolescence may represent general neurodevelopmental alterations that can manifest into different forms of psychiatric conditions later in life. It is worth noting that after correction by the number of external traits tested (new significance threshold: 0.05/4 = 0.0125), only the association between RAVLT immediate recall and schizophrenia remained significant. This could indicate that immediate verbal memory is a better psychosis endophenotype than the other two verbal memory traits. Thus, our findings should be viewed with caution and need to be replicated in the future.
The current study has limitations. Although benefiting from the large adolescent sample of the ABCD Study®, our GWAS is still underpowered to detect SNPs with small effect sizes. Due to the lack of public data available for adolescent verbal memory in diverse populations, we were also not able to conduct replication analysis in an independent sample. Although promising, our findings are preliminary, and future studies with diverse adolescent samples are still needed to replicate our findings. Moreover, as the current methods for estimating SNP heritability and genetic correlations requires LD references from a certain population, we were only able to perform those analyses in the European subset due to the limited sample sizes of other ancestries in the ABCD Study®. Methodological advances are needed in the future to allow for the inclusion of more diverse samples in such analyses. Furthermore, potential sampling bias might exist in the current GWAS. It is possible that children and adolescents with better general cognitive ability were more likely to participate in the ABCD Study®, causing the sample to be biased towards this demographic. Nonetheless, we believe that the study sample was representative of the general population at least in RAVLT performance, as the mean scores were comparable to previous results obtained from children in similar age groups [82,83]. Finally, we did not analyse chromosomal aberrations (e.g., copy number variants) in our study as this requires a different technique. As chromosomal aberrations have been associated with cognitive impairment [24,25], future studies should continue to investigate their influence on specific cognitive domains, such as verbal memory.
In summary, this GWAS, one of the first in children, identified two compelling loci associated with verbal memory that now warrant mechanistic research. We benefited from the diversity within the ABCD Study® and conducted a joint analysis of participants from all ancestries, which maximised the statistical power enabling genome-wide significant findings. Notably, the NSF gene was identified by both SNP-based association testing and gene mapping, indicating that synaptic neurotransmission and membrane fusion might play an important role in verbal memory during neurodevelopment. We also found negative genetic correlations between verbal memory and schizophrenia, suggesting their shared underlying mechanisms. Future research employing novel approaches, including multi-trait analysis of GWAS summary statistics and pathway-specific polygenic risk scores [84,85], will be beneficial in revealing the shared biological pathways between verbal memory, psychosis, and other neuropsychiatric traits.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes1010000/s1. Supplementary methods. Figure S1: Histograms of the three verbal memory traits. Figure S2: The first eight ancestrally representative principal components coloured by ethnicity reported by parents. Figure S3: Relatedness estimation adjusted for ancestral background. Figure S4: Manhattan plots of the genome-wide association analysis of the three verbal memory traits in the European subset. Table S1: Results of the genome-wide significant SNPs across the three verbal memory traits. Table S2: FUMA-mapped genes based on the identified genomic loci.

Author Contributions

Conceptualization: B.W., O.G., H.I. and E.B.; Methodology, B.W., O.G., I.A.-Z. and E.B.; Formal Analysis, B.W.; Investigation, B.W., I.A.-Z., J.H.-S., E.Z., A.B. and E.B.; Data Curation, B.W.; Writing—Original Draft Preparation, B.W.; Writing—Review and Editing, all; Visualization, B.W.; Supervision, E.B., A.M. and K.K.; Project Administration, B.W. and E.B.; Funding Acquisition, B.W. and E.B. All authors have read and agreed to the published version of the manuscript.

Funding

Baihan Wang was supported by the China Scholarship Council-University College London Joint Research Scholarship. Isabelle Austin-Zimmerman, Anjali Bhat, and Jasmine Harju-Seppänen were supported by the Medical Research Council doctoral studentships. Haritz Irizar has received funding from the European Union’s Horizon 220 research and innovation programme under the Marie Sklodowska-Curie grant agreement no 747429 and is currently supported by a grant from the National Institute of Allergy and Infectious Diseases, National Institutes of Health. Elvira Bramon has received the following funding that supported this work: National Institute of Health Research UK (NIHR200756); Mental Health Research UK John Grace QC Scholarship 2018; Economic Social Research Council UK (ESRC) co-funded doctoral award. BMA Margaret Temple Fellowship 2016; Medical Research Council New Investigator and Centenary Awards (G0901310, G1100583), Medical Research Council (G1100583); Wellcome Trust awards (085475/B/08/Z, 085475/Z/08/Z) and NIHR Biomedical Research Centre at University College London Hospitals NHS Foundation Trust and University College London (UCLH BRC—Mental Health Theme).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki. It was exempted from ethical approval at University College London as the research only involved anonymised records and data sets that exist in the public domain. However, the ABCD Study® was approved by the Institutional Review Board of University of California San Diego (IRB# 160091).

Informed Consent Statement

Informed consent was obtained from all participants and their parents in the ABCD Study®.

Data Availability Statement

The data presented in this study are openly available in the ABCD Study® at doi:10.15154/1524269. Summary statistics can be provided by the authors upon request.

Acknowledgments

Data used in the preparation of this article were obtained from the Adolescent Brain Cognitive DevelopmentSM (ABCD) Study (https://abcdstudy.org; accessed on 26 November 2021), held in the NIMH Data Archive (NDA). This is a multisite, longitudinal study designed to recruit more than 10,000 children age 9–10 and follow them over 10 years into early adulthood. The ABCD Study® is supported by the National Institutes of Health and additional federal partners under award numbers U01DA041048, U01DA050989, U01DA051016, U01DA041022, U01DA051018, U01DA051037, U01DA050987, U01DA041174, U01DA041106, U01DA041117, U01DA041028, U01DA041134, U01DA050988, U01DA051039, U01DA041156, U01DA041025, U01DA041120, U01DA051038, U01DA041148, U01DA041093, U01DA041089, U24DA041123, U24DA041147. A full list of supporters is available at https://abcdstudy.org/federal-partners.html (accessed on 26 November 2021). A listing of participating sites and a complete listing of the study investigators can be found at https://abcdstudy.org/consortium_members/ (accessed on 26 November 2021). ABCD consortium investigators designed and implemented the study and/or provided data but did not necessarily participate in the analysis or writing of this report. This manuscript reflects the views of the authors and may not reflect the opinions or views of the NIH or ABCD consortium investigators. The ABCD data repository grows and changes over time. The ABCD data used in this report came from doi:10.15154/1524269. DOIs can be found at https://nda.nih.gov/abcd (accessed on 26 November 2021).

Conflicts of Interest

O.G. became a full-time employee of UCB prior to the completion of this project.

References

  1. Sullivan, P.F.; Kendler, K.S.; Neale, M.C. Schizophrenia as a complex trait: Evidence from a meta-analysis of twin studies. Arch. Gen. Psychiatry 2003, 60, 1187–1192. [Google Scholar] [CrossRef] [Green Version]
  2. Barnett, J.H.; Smoller, J.W. The genetics of bipolar disorder. Neuroscience 2009, 164, 331–343. [Google Scholar] [CrossRef] [Green Version]
  3. 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. [Google Scholar] [CrossRef]
  4. Schizophrenia Working Group of the Psychiatric Genomics Consortium; Ripke, S.; Walters, J.T.; O’Donovan, M.C. Mapping genomic loci prioritises genes and implicates synaptic biology in schizophrenia. medRxiv 2020. [Google Scholar] [CrossRef]
  5. Gottesman, I.I.; Gould, T.D. The endophenotype concept in psychiatry: Etymology and strategic intentions. Am. J. Psychiatry 2003, 160, 636–645. [Google Scholar] [CrossRef]
  6. Gould, T.D.; Gottesman, I.I. Psychiatric endophenotypes and the development of valid animal models. Genes Brain Behav. 2006, 5, 113–119. [Google Scholar] [CrossRef] [PubMed]
  7. Bora, E.; Murray, R.M. Meta-analysis of cognitive deficits in ultra-high risk to psychosis and first-episode psychosis: Do the cognitive deficits progress over, or after, the onset of psychosis? Schizophr. Bull. 2014, 40, 744–755. [Google Scholar] [CrossRef] [Green Version]
  8. Kahn, R.S.; Keefe, R.S. Schizophrenia is a cognitive illness: Time for a change in focus. JAMA Psychiatry 2013, 70, 1107–1112. [Google Scholar] [CrossRef] [PubMed]
  9. Lewis, D.A.; Levitt, P. Schizophrenia as a disorder of neurodevelopment. Annu. Rev. Neurosci. 2002, 25, 409–432. [Google Scholar] [CrossRef]
  10. Rapoport, J.L.; Addington, A.M.; Frangou, S.; Psych, M.R. The neurodevelopmental model of schizophrenia: Update 2005. Mol. Psychiatry 2005, 10, 434–449. [Google Scholar] [CrossRef]
  11. Rapoport, J.L.; Giedd, J.N.; Gogtay, N. Neurodevelopmental model of schizophrenia: Update 2012. Mol. Psychiatry 2012, 17, 1228–1238. [Google Scholar] [CrossRef]
  12. Mesholam-Gately, R.I.; Giuliano, A.J.; Goff, K.P.; Faraone, S.V.; Seidman, L.J. Neurocognition in First-Episode Schizophrenia: A Meta-Analytic Review. Neuropsychology 2009, 23, 315–336. [Google Scholar] [CrossRef] [Green Version]
  13. Bora, E.; Pantelis, C. Meta-analysis of cognitive impairment in first-episode bipolar disorder: Comparison with first-episode schizophrenia and healthy controls. Schizophr. Bull. 2015, 41, 1095–1104. [Google Scholar] [CrossRef] [PubMed]
  14. Sitskoorn, M.M.; Aleman, A.; Ebisch, S.J.; Appels, M.C.; Kahn, R.S. Cognitive deficits in relatives of patients with schizophrenia: A meta-analysis. Schizophr. Res. 2004, 71, 285–295. [Google Scholar] [CrossRef]
  15. Snitz, B.E.; MacDonald, A.W.; Carter, C.S. Cognitive deficits in unaffected first-degree relatives of schizophrenia patients: A meta-analytic review of putative endophenotypes. Schizophr. Bull. 2006, 32, 179–194. [Google Scholar] [CrossRef] [Green Version]
  16. Arts, B.; Jabben, N.; Krabbendam, L.; van Os, J. Meta-analyses of cognitive functioning in euthymic bipolar patients and their first-degree relatives. Psychol. Med. 2008, 38, 771–785. [Google Scholar] [CrossRef]
  17. Bora, E.; Yucel, M.; Pantelis, C. Cognitive endophenotypes of bipolar disorder: A meta-analysis of neuropsychological deficits in euthymic patients and their first-degree relatives. J. Affect. Disord. 2009, 113, 1–20. [Google Scholar] [CrossRef] [PubMed]
  18. Erlenmeyer-Kimling, L.; Rock, D.; Roberts, S.A.; Janal, M.; Kestenbaum, C.; Cornblatt, B.; Adamo, U.H.; Gottesman, I.I. Attention, memory, and motor skills as childhood predictors of schizophrenia-related psychoses: The New York high-risk project. Am. J. Psychiatry 2000, 157, 1416–1422. [Google Scholar] [CrossRef]
  19. Dickson, H.; Cullen, A.E.; Reichenberg, A.; Hodgins, S.; Campbell, D.D.; Morris, R.G.; Laurens, K.R. Cognitive impairment among children at-risk for schizophrenia. J. Psychiatr. Res. 2014, 50, 92–99. [Google Scholar] [CrossRef] [PubMed]
  20. Davies, G.; Lam, M.; Harris, S.E.; Trampush, J.W.; Luciano, M.; Hill, W.D.; Hagenaars, S.P.; Ritchie, S.J.; Marioni, R.E.; Fawns-Ritchie, C.; et al. Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function. Nat. Commun. 2018, 9, 1–16. [Google Scholar] [CrossRef] [Green Version]
  21. Savage, J.E.; Jansen, P.R.; Stringer, S.; Watanabe, K.; Bryois, J.; de Leeuw, C.A.; Nagel, M.; Awasthi, S.; Barr, P.B.; Coleman, J.R.I.; et al. Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence. Nat. Genet. 2018, 50, 912–919. [Google Scholar] [CrossRef] [Green Version]
  22. Trampush, J.W.; Yang, M.L.Z.; Yu, J.; Knowles, E.; Davies, G.; Liewald, D.C.; Starr, J.M.; Djurovic, S.; Melle, I.; Sundet, K.; et al. GWAS meta-analysis reveals novel loci and genetic correlates for general cognitive function: A report from the COGENT consortium. Mol. Psychiatry 2017, 22, 336–345. [Google Scholar] [CrossRef]
  23. Huguet, G.; Schramm, C.; Douard, E.; Jiang, L.; Labbe, A.; Tihy, F.; Mathonnet, G.; Nizard, S.; Lemyre, E.; Mathieu, A.; et al. Measuring and Estimating the Effect Sizes of Copy Number Variants on General Intelligence in Community-Based Samples. JAMA Psychiatry 2018, 75, 447–457. [Google Scholar] [CrossRef]
  24. Ropers, H.H. Genetics of Early Onset Cognitive Impairment. Annu. Rev. Genom. Hum. Genet. 2010, 11, 161–187. [Google Scholar] [CrossRef] [Green Version]
  25. Figura, M.; Geremek, M.; Milanowski, L.M.; Meisner-Kramarz, I.; Duszynska-Was, K.; Szlufik, S.; Rózanski, D.; Smyk, M.; Koziorowski, D. Movement disorders associated with chromosomal aberrations diagnosed in adult patients. Neurol. I Neurochir. Pol. 2021, 55, 300–305. [Google Scholar] [CrossRef] [PubMed]
  26. Ackerman, P.L.; Beier, M.E.; Boyle, M.O. Working memory and intelligence: The same or different constructs? Psychol. Bull. 2005, 131, 30–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Papassotiropoulos, A.; Stephan, D.A.; Huentelman, M.J.; Hoerndli, F.J.; Craig, D.W.; Pearson, J.V.; Huynh, K.D.; Brunner, F.; Corneveaux, J.; Osborne, D.; et al. Common Kibra alleles are associated with human memory performance. Science 2006, 314, 475–478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Papassotiropoulos, A.; Stefanova, E.; Vogler, C.; Gschwind, L.; Ackermann, S.; Spalek, K.; Rasch, B.; Heck, A.; Aerni, A.; Hanser, E.; et al. A genome-wide survey and functional brain imaging study identify CTNNBL1 as a memory-related gene. Mol. Psychiatry 2013, 18, 255–263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Debette, S.; Ibrahim Verbaas, C.A.; Bressler, J.; Schuur, M.; Smith, A.; Bis, J.C.; Davies, G.; Wolf, C.; Gudnason, V.; Chibnik, L.B.; et al. Genome-wide studies of verbal declarative memory in nondemented older people: The Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium. Biol. Psychiatry 2015, 77, 749–763. [Google Scholar] [CrossRef] [Green Version]
  30. Arpawong, T.E.; Pendleton, N.; Mekli, K.; McArdle, J.J.; Gatz, M.; Armoskus, C.; Knowles, J.A.; Prescott, C.A. Genetic variants specific to aging-related verbal memory: Insights from GWASs in a population-based cohort. PLoS ONE 2017, 12, e0182448. [Google Scholar] [CrossRef]
  31. Nakahara, S.; Medland, S.; Turner, J.A.; Calhoun, V.D.; Lim, K.O.; Mueller, B.A.; Bustillo, J.R.; O’Leary, D.S.; Vaidya, J.G.; McEwen, S.; et al. Polygenic risk score, genome-wide association, and gene set analyses of cognitive domain deficits in schizophrenia. Schizophr. Res. 2018, 201, 393–399. [Google Scholar] [CrossRef]
  32. Greenwood, T.A.; Lazzeroni, L.C.; Maihofer, A.X.; Swerdlow, N.R.; Calkins, M.E.; Freedman, R.; Green, M.F.; Light, G.A.; Nievergelt, C.M.; Nuechterlein, K.H.; et al. Genome-wide Association of Endophenotypes for Schizophrenia From the Consortium on the Genetics of Schizophrenia (COGS) Study. JAMA Psychiatry 2019, 170, 521–532. [Google Scholar] [CrossRef] [PubMed]
  33. Sirugo, G.; Williams, S.M.; Tishkoff, S.A. The missing diversity in human genetic studies. Cell 2019, 177, 26–31. [Google Scholar] [CrossRef] [Green Version]
  34. Peterson, R.E.; Kuchenbaecker, K.; Walters, R.K.; Chen, C.Y.; Popejoy, A.B.; Periyasamy, S.; Lam, M.; Iyegbe, C.; Strawbridge, R.J.; Brick, L.; et al. Genome-wide association studies in ancestrally diverse populations: Opportunities, methods, pitfalls, and recommendations. Cell 2019, 179, 589–603. [Google Scholar] [CrossRef] [PubMed]
  35. Popejoy, A.B.; Fullerton, S.M. Genomics is failing on diversity. Nature 2016, 538, 161–164. [Google Scholar] [CrossRef] [Green Version]
  36. Dennison, C.A.; Legge, S.E.; Pardiñas, A.F.; Walters, J.T. Genome-wide association studies in schizophrenia: Recent advances, challenges and future perspective. Schizophr. Res. 2020, 217, 4–12. [Google Scholar] [CrossRef] [PubMed]
  37. Conomos, M.P.; Reiner, A.P.; McPeek, M.S.; Thornton, T.A. Genome-wide control of population structure and relatedness in genetic association studies via linear mixed models with orthogonally partitioned structure. bioRxiv 2018, 409953. [Google Scholar] [CrossRef]
  38. Gogarten, S.M.; Sofer, T.; Chen, H.; Yu, C.; Brody, J.A.; Thornton, T.A.; Rice, K.M.; Conomos, M.P. Genetic association testing using the GENESIS R/Bioconductor package. Bioinformatics 2019, 35, 5346–5348. [Google Scholar] [CrossRef]
  39. Garavan, H.; Bartsch, H.; Conway, K.; Decastro, A.; Goldstein, R.Z.; Heeringa, S.; Jernigan, T.; Potter, A.; Thompson, W.; Zahs, D. Recruiting the ABCD sample: Design considerations and procedures. Dev. Cogn. Neurosci. 2018, 32, 16–22. [Google Scholar] [CrossRef]
  40. Barch, D.M.; Albaugh, M.D.; Avenevoli, S.; Chang, L.; Clark, D.B.; Glantz, M.D.; Hudziak, J.J.; Jernigan, T.L.; Tapert, S.F.; Yurgelun-Todd, D.; et al. Demographic, physical and mental health assessments in the adolescent brain and cognitive development study: Rationale and description. Dev. Cogn. Neurosci. 2018, 32, 55–66. [Google Scholar] [CrossRef]
  41. Uban, K.A.; Horton, M.K.; Jacobus, J.; Heyser, C.; Thompson, W.K.; Tapert, S.F.; Madden, P.A.; Sowell, E.R. Biospecimens and the ABCD study: Rationale, methods of collection, measurement and early data. Dev. Cogn. Neurosci. 2018, 32, 97–106. [Google Scholar] [CrossRef]
  42. Thompson, W.K.; Barch, D.M.; Bjork, J.M.; Gonzalez, R.; Nagel, B.J.; Nixon, S.J.; Luciana, M. The structure of cognition in 9 and 10 year-old children and associations with problem behaviors: Findings from the ABCD study’s baseline neurocognitive battery. Dev. Cogn. Neurosci. 2019, 36, 100606. [Google Scholar] [CrossRef]
  43. Karcher, N.R.; O’Brien, K.J.; Kandala, S.; Barch, D.M. Resting-State Functional Connectivity and Psychotic-like Experiences in Childhood: Results From the Adolescent Brain Cognitive Development Study. Biol. Psychiatry 2019, 86, 7–15. [Google Scholar] [CrossRef]
  44. Schmidt, M. Rey Auditory Verbal Learning Test: A Handbook; Western Psychological Services: Torrance, CA, USA, 1996. [Google Scholar]
  45. Baurley, J.W.; Edlund, C.K.; Pardamean, C.I.; Conti, D.V.; Bergen, A.W. Smokescreen: A targeted genotyping array for addiction research. BMC Genom. 2016, 17, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Lam, M.; Awasthi, S.; Watson, H.J.; Goldstein, J.; Panagiotaropoulou, G.; Trubetskoy, V.; Karlsson, R.; Frei, O.; Fan, C.C.; De Witte, W.; et al. RICOPILI: Rapid Imputation for COnsortias PIpeLIne. Bioinformatics 2020, 36, 930–933. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Auton, A.; Abecasis, G.R.; Altshuler, D.M.; Durbin, R.M.; Bentley, D.R.; Chakravarti, A.; Clark, A.G.; Donnelly, P.; Eichler, E.E.; Flicek, P.; et al. A global reference for human genetic variation. Nature 2015, 526, 68–74. [Google Scholar] [CrossRef] [Green Version]
  48. McCarthy, S.; Das, S.; Kretzschmar, W.; Delaneau, O.; Wood, A.R.; Teumer, A.; Kang, H.M.; Fuchsberger, C.; Danecek, P.; Sharp, K.; et al. A reference panel of 64,976 haplotypes for genotype imputation. Nat. Genet. 2016, 48, 1279–1283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Loh, P.R.; Palamara, P.F.; Price, A.L. Fast and accurate long-range phasing in a UK Biobank cohort. Nat. Genet. 2016, 48, 811–816. [Google Scholar] [CrossRef] [Green Version]
  50. Das, S.; Forer, L.; Schönherr, S.; Sidore, C.; Locke, A.E.; Kwong, A.; Vrieze, S.I.; Chew, E.Y.; Levy, S.; McGue, M.; et al. Next-generation genotype imputation service and methods. Nat. Genet. 2016, 48, 1284–1287. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 2011, 27, 2987–2993. [Google Scholar] [CrossRef] [Green Version]
  52. Purcell, S.M.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2020. [Google Scholar]
  54. Manichaikul, A.; Mychaleckyj, J.C.; Rich, S.S.; Daly, K.; Sale, M.; Chen, W.M. Robust relationship inference in genome-wide association studies. Bioinformatics 2010, 26, 2867–2873. [Google Scholar] [CrossRef] [Green Version]
  55. Zheng, X.; Levine, D.; Shen, J.; Gogarten, S.M.; Laurie, C.; Weir, B.S. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics 2012, 28, 3326–3328. [Google Scholar] [CrossRef] [Green Version]
  56. Conomos, M.P.; Reiner, A.P.; Weir, B.S.; Thornton, T.A. Model-free Estimation of Recent Genetic Relatedness. Am. J. Hum. Genet. 2016, 98, 127–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Watanabe, K.; Taskesen, E.; Van Bochoven, A.; Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 2017, 8, 1–11. [Google Scholar] [CrossRef] [Green Version]
  58. Lonsdale, J.; Thomas, J.; Salvatore, M.; Phillips, R.; Lo, E.; Shad, S.; Hasz, R.; Walters, G.; Garcia, F.; Young, N.; et al. The Genotype-Tissue Expression (GTEx) project. Nat. Genet. 2013, 45, 580–585. [Google Scholar] [CrossRef]
  59. Schmitt, A.D.; Hu, M.; Jung, I.; Xu, Z.; Qiu, Y.; Tan, C.L.; Li, Y.; Lin, S.; Lin, Y.; Barr, C.L.; et al. A Compendium of Chromatin Contact Maps Reveals Spatially Active Regions in the Human Genome. Cell Rep. 2016, 17, 2042–2059. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Giusti-Rodríguez, P.; Lu, L.; Yang, Y.; Crowley, C.A.; Liu, X.; Juric, I.; Martin, J.S.; Abnousi, A.; Allred, S.C.; Ancalade, N.; et al. Using three-dimensional regulatory chromatin interactions from adult and fetal cortex to interpret genetic results for psychiatric disorders and cognitive traits. bioRxiv 2019, 18, 406330. [Google Scholar] [CrossRef] [Green Version]
  61. Kundaje, A.; Meuleman, W.; Ernst, J.; Bilenky, M.; Yen, A.; Heravi-Moussavi, A.; Kheradpour, P.; Zhang, Z.; Wang, J.; Ziller, M.J.; et al. Integrative analysis of 111 reference human epigenomes. Nature 2015, 518, 317–330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Howard, D.M.; Adams, M.J.; Clarke, T.K.; Hafferty, J.D.; Gibson, J.; Shirali, M.; Coleman, J.R.; Hagenaars, S.P.; Ward, J.; Wigmore, E.M.; et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat. Neurosci. 2019, 22, 343–352. [Google Scholar] [CrossRef] [Green Version]
  63. Lee, J.J.; Wedow, R.; Okbay, A.; Kong, E.; Maghzian, O.; Zacher, M.; Nguyen-Viet, T.A.; Bowers, P.; Sidorenko, J.; Karlsson Linnér, R.; et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat. Genet. 2018, 50, 1112–1121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Ward, L.D.; Kellis, M. HaploReg v4: Systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 2016, 44, D877–D881. [Google Scholar] [CrossRef] [PubMed]
  65. Ben-Eghan, C.; Sun, R.; Hleap, J.S.; Diaz-Papkovich, A.; Munter, H.M.; Grant, A.V.; Dupras, C.; Gravel, S. Don’t ignore genetic data from minority populations. Nature 2020, 585, 184–186. [Google Scholar] [CrossRef] [PubMed]
  66. Richard, D.J.; Bolderson, E.; Cubeddu, L.; Wadsworth, R.I.; Savage, K.; Sharma, G.G.; Nicolette, M.L.; Tsvetanov, S.; McIlwraith, M.J.; Pandita, R.K.; et al. Single-stranded DNA-binding protein hSSB1 is critical for genomic stability. Nature 2008, 453, 677–681. [Google Scholar] [CrossRef] [PubMed]
  67. Horie, M.; Mitsumoto, Y.; Kyushiki, H.; Kanemoto, N.; Watanabe, A.; Taniguchi, Y.; Nishino, N.; Okamoto, T.; Kondo, M.; Mori, T.; et al. Identification and characterization of TMEFF2, a novel survival factor for hippocampal and mesencephalic neurons. Genomics 2000, 67, 146–152. [Google Scholar] [CrossRef] [PubMed]
  68. Hong, H.S.; Maezawa, I.; Petrlova, J.; Zhao, X.Y.; Voss, J.C.; Jin, L.W. Tomoregulin (TMEFF2) Binds Alzheimer’s Disease Amyloid-β (Aβ) Oligomer and AβPP and Protects Neurons from Aβ-Induced Toxicity. J. Alzheimer’s Dis. 2015, 48, 731–743. [Google Scholar] [CrossRef] [Green Version]
  69. Xie, S.; Zhang, Y.; Peng, T.; Guo, J.; Cao, Y.; Guo, J.; Shi, X.; Li, Y.; Liu, Y.; Qi, S.; et al. TMEFF2 promoter hypermethylation is an unfavorable prognostic marker in gliomas. Cancer Cell Int. 2021, 21, 1–11. [Google Scholar] [CrossRef]
  70. Hay, J.C.; Scheller, R.H. SNARES and NSF in targeted membrane fusion. Curr. Opin. Cell Biol. 1997, 9, 505–512. [Google Scholar] [CrossRef]
  71. Zhao, B.; Luo, T.; Li, T.; Li, Y.; Zhang, J.; Shan, Y.; Wang, X.; Yang, L.; Zhou, F.; Zhu, Z.; et al. Genome-wide association analysis of 19,629 individuals identifies variants influencing regional brain volumes and refines their genetic co-architecture with cognitive and mental health traits. Nat. Genet. 2019, 51, 1637–1644. [Google Scholar] [CrossRef]
  72. Hofer, E.; Roshchupkin, G.V.; Adams, H.H.; Knol, M.J.; Lin, H.; Li, S.; Zare, H.; Ahmad, S.; Armstrong, N.J.; Satizabal, C.L.; et al. Genetic correlations and genome-wide associations of cortical structure in general population samples of 22,824 adults. Nat. Commun. 2020, 11, 1–16. [Google Scholar] [CrossRef] [PubMed]
  73. Nalls, M.A.; Blauwendraat, C.; Vallerga, C.L.; Heilbron, K.; Bandres-Ciga, S.; Chang, D.; Tan, M.; Kia, D.A.; Noyce, A.J.; Xue, A.; et al. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: A meta-analysis of genome-wide association studies. Lancet Neurol. 2019, 18, 1091–1102. [Google Scholar] [CrossRef]
  74. Nagel, M.; Jansen, P.R.; Stringer, S.; Watanabe, K.; De Leeuw, C.A.; Bryois, J.; Savage, J.E.; Hammerschlag, A.R.; Skene, N.G.; Muñoz-Manchado, A.B.; et al. Meta-analysis of genome-wide association studies for neuroticism in 449,484 individuals identifies novel genetic loci and pathways. Nat. Genet. 2018, 50, 920–927. [Google Scholar] [CrossRef] [Green Version]
  75. Ramos-Miguel, A.; Beasley, C.L.; Dwork, A.J.; Mann, J.J.; Rosoklija, G.; Barr, A.M.; Honer, W.G. Increased SNARE Protein-Protein Interactions in Orbitofrontal and Anterior Cingulate Cortices in Schizophrenia. Biol. Psychiatry 2015, 78, 361–373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Sharp, A.J.; Hansen, S.; Selzer, R.R.; Cheng, Z.; Regan, R.; Hurst, J.A.; Stewart, H.; Price, S.M.; Blair, E.; Hennekam, R.C.; et al. Discovery of previously unidentified genomic disorders from the duplication architecture of the human genome. Nat. Genet. 2006, 38, 1038–1042. [Google Scholar] [CrossRef]
  77. Shaw-Smith, C.; Pittman, A.M.; Willatt, L.; Martin, H.; Rickman, L.; Gribble, S.; Curley, R.; Cumming, S.; Dunn, C.; Kalaitzopoulos, D.; et al. Microdeletion encompassing MAPT at chromosome 17q21.3 is associated with developmental delay and learning disability. Nat. Genet. 2006, 38, 1032–1037. [Google Scholar] [CrossRef]
  78. Koolen, D.A.; Vissers, L.E.; Pfundt, R.; De Leeuw, N.; Knight, S.J.; Regan, R.; Kooy, R.F.; Reyniers, E.; Romano, C.; Fichera, M.; et al. A new chromosome 17q21.31 microdeletion syndrome associated with a common inversion polymorphism. Nat. Genet. 2006, 38, 999–1001. [Google Scholar] [CrossRef] [PubMed]
  79. Bertisch, H.; Li, D.; Hoptman, M.J.; DeLisi, L.E. Heritability Estimates for Cognitive Factors and Brain White Matter Integrity as Markers of Schizophrenia. Am. J. Med. Genet. 2010, 153B, 885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Manolio, T.A.; Collins, F.S.; Cox, N.J.; Goldstein, D.B.; Hindorff, L.A.; Hunter, D.J.; McCarthy, M.I.; Ramos, E.M.; Cardon, L.R.; Chakravarti, A.; et al. Finding the missing heritability of complex diseases. Nature 2009, 461, 747–753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Anttila, V.; Bulik-Sullivan, B.; Finucane, H.K.; Walters, R.K.; Bras, J.; Duncan, L.; Escott-Price, V.; Falcone, G.J.; Gormley, P.; Malik, R.; et al. Analysis of shared heritability in common disorders of the brain. Science 2018, 360, 1–40. [Google Scholar]
  82. van den Burg, W.; Kingma, A. Performance of 225 Dutch School Children on Rey’s Auditory Verbal Learning Test (AVLT): Parallel Test-Retest Reliabilities with an Interval of 3 Months and Normative Data. Arch. Clin. Neuropsychol. 1999, 14, 545–559. [Google Scholar] [CrossRef]
  83. Forrester, G.; Geffen, G. Performance measures of 7- to 15-year-old children on the Auditory Verbal Learning Test. Clin. Neuropsychol. 1991, 5, 345–359. [Google Scholar] [CrossRef]
  84. Turley, P.; Walters, R.K.; Maghzian, O.; Okbay, A.; Lee, J.J.; Fontana, M.A.; Nguyen-Viet, T.A.; Wedow, R.; Zacher, M.; Furlotte, N.A.; et al. Multi-trait analysis of genome-wide association summary statistics using MTAG. Nat. Genet. 2018, 50, 229–237. [Google Scholar] [CrossRef] [PubMed]
  85. Darst, B.F.; Koscik, R.L.; Racine, A.M.; Oh, J.M.; Krause, R.A.; Carlsson, C.M.; Zetterberg, H.; Blennow, K.; Christian, B.T.; Bendlin, B.B.; et al. Pathway-Specific Polygenic Risk Scores as Predictors of Amyloid-β Deposition and Cognitive Function in a Sample at Increased Risk for Alzheimer’s Disease. J. Alzheimer’s Dis. 2017, 55, 473–484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 13 00106 g001 550
Figure 1. Manhattan plots of the genome-wide association analysis of the three verbal memory traits in the whole sample. The dotted line indicates the genome-wide significance threshold of p = 5 × 10−8. RAVLT, Rey Auditory Verbal Learning Test.
Figure 1. Manhattan plots of the genome-wide association analysis of the three verbal memory traits in the whole sample. The dotted line indicates the genome-wide significance threshold of p = 5 × 10−8. RAVLT, Rey Auditory Verbal Learning Test.
Genes 13 00106 g001
Table
Table 1. Demographic and cognitive characteristics of the Adolescent Brain Cognitive Development sample.
Table 1. Demographic and cognitive characteristics of the Adolescent Brain Cognitive Development sample.
Variablen (11,017 in Total)%
Sex
Female517647
Male580553
Mean (SD) Age (years)9.90.6
Ethnicity
White587253
Hispanic209919
Black168415
Asian1551
AIAN/NHPI and mixed116911
Mean (SD) RAVLT score
Immediate recall44.29.9
Short-delay recall9.73
Long-delay recall9.23.2
Note. There were 36 missing data for sex and age, 38 for ethnicity, 291 for RAVLT immediate recall, 216 for RAVLT short-delay recall, and 260 for RAVLT long-delay recall. SD, standard deviation. AIAN, American Indian/Alaska Native. NHPI, Native Hawaiian and other Pacific Islander. RAVLT, Rey Auditory Verbal Learning Test.
Table
Table 2. Genomic loci associated with verbal memory.
Table 2. Genomic loci associated with verbal memory.
PhenotypeLead SNPChromosome (Position)Reference AlleleEffect Allele (Frequency)Beta (SE)pr2
RAVLT short-delay recallrs739845662 (191566282)GA (0.019)−0.847 (0.151)1.86 × 10−80.91
RAVLT long-delay recallrs989624317 (46748690)CG (0.185)−0.309 (0.055)2.22 × 10−80.96
Note. Genomic positions are based on GRCh38. Frequency represents allele frequency in the sample. RAVLT, Rey Auditory Verbal Learning Test. SE, standard error. r2, imputation quality score.
Table
Table 3. Genetic correlations between three verbal memory traits and relevant psychiatric traits or educational attainment.
Table 3. Genetic correlations between three verbal memory traits and relevant psychiatric traits or educational attainment.
Verbal Memory TraitsSchizophreniaBipolar DisorderMajor Depressive DisorderEducational Attainment
rg (SE)prg (SE)prg (SE)prg (SE)p
RAVLT immediate recall−0.29 (0.10)0.003−0.18 (0.10)0.075−0.14 (0.09)0.0920.53 (0.14)0.001
RAVLT short-delay recall−0.15 (0.06)0.017−0.11 (0.06)0.082−0.14 (0.07)0.0500.31 (0.06)<0.001
RAVLT long-delay recall−0.15 (0.08)0.063−0.13 (0.08)0.083−0.07 (0.08)0.3700.40 (0.08)<0.001
Note. RAVLT, Rey Auditory Verbal Learning Test. rg, genetic correlation. SE, standard error.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, B.; Giannakopoulou, O.; Austin-Zimmerman, I.; Irizar, H.; Harju-Seppänen, J.; Zartaloudi, E.; Bhat, A.; McQuillin, A.; Kuchenbäcker, K.; Bramon, E. Adolescent Verbal Memory as a Psychosis Endophenotype: A Genome-Wide Association Study in an Ancestrally Diverse Sample. Genes 2022, 13, 106. https://doi.org/10.3390/genes13010106

AMA Style

Wang B, Giannakopoulou O, Austin-Zimmerman I, Irizar H, Harju-Seppänen J, Zartaloudi E, Bhat A, McQuillin A, Kuchenbäcker K, Bramon E. Adolescent Verbal Memory as a Psychosis Endophenotype: A Genome-Wide Association Study in an Ancestrally Diverse Sample. Genes. 2022; 13(1):106. https://doi.org/10.3390/genes13010106

Chicago/Turabian Style

Wang, Baihan, Olga Giannakopoulou, Isabelle Austin-Zimmerman, Haritz Irizar, Jasmine Harju-Seppänen, Eirini Zartaloudi, Anjali Bhat, Andrew McQuillin, Karoline Kuchenbäcker, and Elvira Bramon. 2022. "Adolescent Verbal Memory as a Psychosis Endophenotype: A Genome-Wide Association Study in an Ancestrally Diverse Sample" Genes 13, no. 1: 106. https://doi.org/10.3390/genes13010106

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