Genetic copy number variants, cognition and psychosis: a meta analysis and a family study

Molecular Psychiatry (2021) 26:5307–5319 https://doi.org/10.1038/s41380-020-0820-7 ARTICLE Genetic copy number variants, cognition and psychosis: a metaanalysis and a family study Johan H. Thygesen 1●Amelia Presman1●Jasmine Harju-Seppänen 1●Haritz Irizar1●Rebecca Jones1● Karoline Kuchenbaecker 1,2 ●Kuang Lin3,4 ●Behrooz Z. Alizadeh5,6 ●Isabelle Austin-Zimmerman1● Agna Bartels-Velthuis5●Anjali Bhat 1●Richard Bruggeman5,7 ●Wiepke Cahn8,9 ●Stella Calafato1● Benedicto Crespo-Facorro10,11,12 ●Liewe de Haan13,14 ●Sonja M. C. de Zwarte8●Marta Di Forti3● Álvaro Díez-Revuelta1,15 ●Jeremy Hall 16 ●Mei-Hua Hall17 ●Conrad Iyegbe3●Assen Jablensky18 ●Rene Kahn 8,19 ● Luba Kalaydjieva20 ●Eugenia Kravariti3●Stephen Lawrie 21 ●Jurjen J. Luykx 8,22,23 ●Igancio Mata10,24 ● Colm McDonald 25 ●Andrew M. McIntosh 21,26 ●Andrew McQuillin 1●Rebecca Muir1●Roel Ophoff27,28,29 ● Marco Picchioni3●Diana P. Prata3,30,31 ●Siri Ranlund1,3 ●Dan Rujescu32,33 ●Bart P. F. Rutten 34,35 ●Katja Schulze3,36 ● Madiha Shaikh37,38 ●Frederike Schirmbeck13,14 ●Claudia J. P. Simons34,39 ●Timothea Toulopoulou3,40 ● Therese van Amelsvoort34 ●Neeltje van Haren8,41 ●Jim van Os3,34,42 ●Ruud van Winkel43 ●Evangelos Vassos 3● Muriel Walshe1●Matthias Weisbrod44,45 ●Eirini Zartaloudi1●Vaughan Bell1●John Powell3●Cathryn M. Lewis 3● Robin M. Murray3,36 ●Elvira Bramon 1,3,46 Received: 7 August 2018 / Revised: 11 June 2020 / Accepted: 11 June 2020 / Published online: 27 July 2020 © The Author(s) 2020. This article is published with open access Abstract The burden of large and rare copy number genetic variants (CNVs) as well as certain specific CNVs increase the risk of developing schizophrenia. Several cognitive measures are purported schizophrenia endophenotypes and may represent an intermediate point between genetics and the illness. This paper investigates the influence of CNVs on cognition. We conducted a systematic review and meta-analysis of the literature exploring the effect of CNV burden on general intelligence. We included ten primary studies with a total of 18,847 participants and found no evidence of association. In a new psychosis family study, we investigated the effects of CNVs on specific cognitive abilities. We examined the burden of large and rare CNVs (>200 kb, <1% MAF) as well as known schizophrenia-associated CNVs in patients with psychotic disorders, their unaffected relatives and controls (N=3428) from the Psychosis Endophenotypes International Consortium (PEIC). The carriers of specific schizophrenia-associated CNVs showed poorer performance than non-carriers in immediate (P=0.0036) and delayed (P=0.0115) verbal recall. We found suggestive evidence that carriers of schizophrenia-associated CNVs had poorer block design performance (P=0.0307). We do not find any association between CNV burden and cognition. Our findings show that the known high-risk CNVs are not only associated with schizophrenia and other neurodevelopmental disorders, but are also a contributing factor to impairment in cognitive domains such as memory and perceptual reasoning, and act as intermediate biomarkers of disease risk. Introduction Copy number variants (CNVs) occur if sections of DNA sequence become deleted or duplicated [1–3]. Although many CNVs are benign and contribute to natural human variation [4], larger and rarer variants are more likely to be pathogenic and under negative selection pressure [5,6]. The phenotypic effects of CNVs are not fully understood, but they influence neurodevelopment, cognitive abilities and the risk of several common brain disorders [6]. These authors contributed equally: Johan H. Thygesen, Amelia Presman, Jasmine Harju-Seppänen *Elvira Bramon [email protected] Extended author information available on the last page of the article Supplementary information The online version of this article (https:// doi.org/10.1038/s41380-020-0820-7) contains supplementary material, which is available to authorized users. 1234567890();,: 1234567890();,: Specific CNV loci are associated with increased risk of developing schizophrenia [7–11]. A recent large CNV meta-analysis by the Psychiatric Genomics Consortium showed robust genome-wide significant associations for eight loci as well as suggestive support for an additional nine [12]. Schizophrenia-associated CNVs have incomplete penetrance [6,13], and are rare, hence most people with schizophrenia are not carriers. However, schizophreniaassociated CNVs have odds ratios ranging from 2 to 30 [12,14] and thus constitute some of the strongest known risk factors for the illness. An increased burden of large and rare CNVs has also been associated with schizophrenia [15,16]. Studies have shown that, compared with healthy controls, individuals with schizophrenia carry a greater number of rare (<1% frequency) CNVs of over 20 kilobases (kb) [12], 100 kb [15–17], 200 kb [16,17], 500 kb [16,17] and 1 Mb [18]. The largest study to date, by the Psychiatric Genomics Consortium, further shows that the burden is enriched for genes associated with synaptic function and that deletions assert greater effect than duplications [12]. Despite strong evidence that CNVs are risk factors for schizophrenia and other developmental disorders, the mechanisms by which CNVs lead to disease onset remain unclear. Endophenotypes are biomarkers that characterise illnesses and indicate genetic liability, as intermediate steps on the pathway from genes to disease [19,20]. Cognitive function is one such endophenotype for schizophrenia and extensive literature shows that individuals with schizophrenia display reduced performance across a range of tests of specific and general cognition [21,22]. This is not simply due to the effects of antipsychotic medication [23], and nor is it just an epiphenomenon of the symptoms of schizophrenia; cognition is impaired before illness onset [24,25] as well as amongst the unaffected relatives of patients with schizophrenia [26–28]. A recent genome-wide association study with more than 269,000 samples shows a bidirectional effect with intelligence having a strong protective effect towards schizophrenia risk, and a smaller reverse effect, with schizophrenia predisposing to impaired cognitive functioning [29]. IQ and general cognitive ability are heritable [30–33]; however despite the identification of 205 loci affecting over 1000 genes associated with intelligence [29], they only explain ~5% of the inter-individual variability in intelligence. Part of the unexplained heritability of intelligence could be attributable to copy number variants. Many CNVs affect genes involved in neurodevelopment [15,34,35], providing a mechanism by which specific CNVs and CNV burden could affect cognition. There is evidence linking several specific CNVs with schizophrenia, other neurodevelopmental disorders, educational attainment [36,37] and with impaired cognition [38–41]. Furthermore, Stefansson et al. [42] showed that healthy carriers of any of 26 neuropsychiatric CNVs collectively performed at an intermediate level between healthy non-carriers and schizophrenia patients in several cognitive tests. This indicates that, while the risk CNVs may not have full penetrance for disease, most carriers will exhibit some degree of phenotypic change such as impaired cognition. A large study on the UK Biobank further supports this effect of neuropsychiatric CNVs impairing cognition in healthy carriers [37]. While the detrimental effects of specific schizophreniaassociated CNVs on cognition are well characterised, the influence of CNV burden on cognition is less clear. Some evidence, both in clinical samples and healthy populations, suggests that increased CNV burden is associated with lower IQ [30,31,43,44], while other studies have failed to find this association [35,45–47]. Until now, few studies have reported the effects of schizophreniaassociated CNVs or CNV burden on specific cognitive abilities [37,42,48]. Firstly, we conducted a systematic review and metaanalysis of the literature examining the relationship between CNV burden and general cognitive ability. We then present data from a new family study from the Psychosis Endophenotypes International Consortium (PEIC) [49] investigating the influence of CNVs (both burden and loci) on cognitive endophenotypes for schizophrenia [27,50]. Methods Meta-analysis of published association studies of CNVs and general IQ We conducted a literature search using the databases Pubmed, Medline, and PsychINFO using the following search terms: “(CNV* OR copy number OR copy-number) AND (IQ OR intelligence OR cogniti*)”. The time window included any papers published before 1st April 2019. The reference and citations lists of relevant papers were examined to identify other relevant papers. We imposed no restriction on participant age, geographical location, or article language. In addition to IQ, we included papers examining other measures of general intelligence as they are thought to be closely linked [51]. Papers investigating both patient and healthy populations were included. Where different papers used the same sample of participants, only the study with the most relevant phenotype was included. If multiple measures of intelligence were included (for example see references [47,52]) the one deemed closest to the other studies was used for the meta-analysis. Similarly, if measures for both common and rare CNV-burden were reported [52] we included the latter for the meta-analysis. 5308 J. H. Thygesen et al. Titles and abstracts of all relevant papers were screened to assess whether they met the inclusion criteria. Where necessary, we contacted the authors to request additional information needed to include the study in the metaanalysis. Supplementary Table 1 shows the data extracted from papers. Meta-analysis of the literature A random effects meta-analysis was conducted using StatsDirect version 3.0 [53] to calculate an overall estimate of the correlation between CNV burden and IQ for the included studies. Where primary studies reported Spearman’s correlations or standardised coefficients they were converted to Pearson’s correlation coefficients to ensure the studies were as comparable as possible [54,55]. A random effects meta-analysis was chosen due to the variability in methods of the included studies (including different participant samples and inclusion criteria for CNVs). Statistical heterogeneity was measured using Cochran’s Q statistic. CNV analysis of the Psychosis Endophenotypes International Consortium sample The initial dataset (prior to quality control) included 5597 participants from the PEIC family study [49], including people with schizophrenia, bipolar disorder with psychotic symptoms and other forms of psychosis, their unaffected relatives and unrelated controls. Participants were of European ancestry and assessments were conducted at nine centres: Amsterdam, Edinburgh, Groningen, London, Maastricht, Munich, Pamplona, Perth and Utrecht (see Supplementary Table 2 for further detail). Participants were recruited through clinical teams, voluntary organisations and press advertisements, and contributed both genetic data and cognitive performance measures [49]. All participants provided written informed consent and the study was approved by the respective ethical committees at each of the participating centres. Details of diagnostic classifications can be found in the Supplementary Methods. Genotyping and quality control DNA was obtained from blood for all participants and sent to the Wellcome Trust Sanger Institute (Cambridge, United Kingdom). Samples were genotyped with the Human SNP Array 6.0 at the Affymetrix Services Laboratory (www. affymetrix.com). We applied standard quality control procedures as described in the Supplementary Methods and in Bramon et al. [49]. CNVs were identified using PennCNV [56] and Affymetrix Power Tools, following the PennCNV- Affy protocol to calculate log R ratio (LRR) and B-allele frequency (BAF). Standard PennCNV settings were used and data were adjusted for genomic waves [57] using Affymetrix 6.0 GC-model file. Individual-based quality control for the CNVs was performed using statistics calculated with PennCNV: quality control thresholds were determined based on inspection of the frequency distributions of the BAF-drift, LRR-standard deviation, and number of CNVs per participant, respectively. Individuals with either BAF-drift of >0.01, LRR- standard deviation of >0.5 or more than 300 CNV calls were removed. CNV-level quality control was performed by excluding CNVs with ten or fewer SNPs and by iteratively merging adjacent calls together if the length between calls was <20% of the combined length. Calls made by PennCNV in the pseudo-autosomal regions of the X- chromosome (10,000–2,781,479 bp and 155,701,382–156,030,895 bp, hg18) were excluded. All CNVs included in the analysis were visually inspected by two researchers blind to clinical data using an in-house script to visualise BAF and LRR patterns. A consensus decision on inclusion was made between two researchers, both blind to clinical data, based on the comparison of the observed LRR and BAF of the affected region with the expected for a CNV with the given copy state. PennCNV frequently made CNV predictions that did not fit with the expected allelic and/or intensity pattern of the given copy state and thus 72% of CNV calls were discarded. CNV burden analysis Only rare (<1% frequency in the sample) and large (>200 kb) CNVs were included in the CNV burden. Frequency of CNVs were determined by identifying common CNV loci, through independent mapping of start and stop positions of all CNVs. CNVs whose start positions were within 300,000 bp of each other, where the stop position was also within a 300,000 bp bin were considered to be the same loci (see Supplementary Fig. 1 for details). In addition to calculations using total length of CNVs, we also measured burden as numbers of genes affected, as described in Marshall et al. [12]. We did this by annotating CNVs with RefSeq genes (hg18), including genes where at least one base pair of an exon overlapped with the CNV and adding up all unique genes affected in each individual. Total burden, and deletion and duplication burdens were analysed separately. Analysis of schizophrenia-associated CNVs We searched for carriers of 27 CNVs with good evidence of an association with schizophrenia as described by Marshall et al. [12], Kirov et al. [6] and Stefansson et al. [42] (see Supplementary Table 3). For the analysis of schizophreniaassociated CNVs we considered all samples prior to any Genetic copy number variants, cognition and psychosis: a meta-analysis and a family study 5309 sample/CNV-level quality control and performed visual inspection of all samples with CNV calls that overlapped with >10% of a schizophrenia-associated locus. Samples identified with true schizophrenia-associated CNVs by consensus of two researchers blind to clinical data, were included in the analysis even if they failed sample/CNV- level quality control. For 2p16 deletions, all CNV calls overlapping the region were inspected, and participants with validated CNVs affecting exons of the causative gene NRXN1 [58] were identified as carriers. Cognitive measures Cognitive measures collected from participants included block design [59,60] (a test of perceptual reasoning), the combined digit span (measuring attention and working memory), and the Rey Auditory Verbal Learning Task (RAVLT) immediate and delayed recall (measuring short and long-term verbal memory, respectively). As different versions of these tests were used across centres, participants’raw scores were converted into percentages by dividing each participant’s score by the maximum achievable score and multiplying by 100. Supplementary Table 2 details the number of participants for each cognitive measure. Kinship matrix The kinship coefficient is a probabilistic estimate that a random gene from a subject is “identical by descent”to a gene in the same locus from another subject. For “n”subjects, these probabilities can be assembled in an n×n “kinship matrix”, which can be used to model the covariance (or “relatedness”) between individuals and the population structure in a dataset. A kinship matrix based on a LD-pruned set of SNPs (102,112 SNPs selected with pruning parameters: r2=0.2; window =1000 kb) was generated using LDAK [61] and added as a random effect to the linear and logistic mixed model regressions. Statistical analysis Firstly, the CNV burden data were not normally distributed showing clear zero inflation as expected. Therefore we performed a Wilcoxon rank sum test to compare the CNV burden of individuals with and without cognitive data available. Secondly, the association between the CNV measures and clinical group was investigated. For this analysis, we performed mixed effects logistic regressions with disease status as outcome (either patients versus controls or relatives versus controls), age, gender and study centre as fixed effects and the kinship matrix as a random effect. Thirdly, the relationships between known schizophreniaassociated CNVs and quantitative cognitive measures were examined using linear mixed models. The outcome variable was the cognitive measure and the predictor was the participants’carrier status of schizophrenia-associated CNVs (carriers versus non-carriers). Given that schizophrenia-associated CNVs are very rare, with frequencies ranging from 0.01 to 0.3% [6,42], only a combined analysis including several such CNVs was feasible. The analysis for a particular cognitive measure was only performed if data from at least ten carriers of schizophrenia-associated CNVs with cognitive measures were available. Finally, we also used linear mixed models to investigate the association between cognitive measures and CNV burden. Age, gender, clinical group (patient, relative and control), study centre and a kinship matrix were included as covariates in all linear mixed models. Analyses were performed using R version 3.5.0 [62]. All mixed model regressions including the kinship matrix as a random effect were performed using the lme4qtl R package [63]. As is standard practice in genetic association studies, we adjusted the significance threshold for multiple testing. We used two different analysis approaches depending on the whether the outcome variable was quantitative or categorical. Firstly, for the linear mixed models, we tested the correlations between all the cognitive outcomes and divided the significance threshold (0.05) by the effective number of traits, as per standard method [64,65]. We investigated four cognitive measures: digit span, block design, RAVLT immediate, RAVLT delayed. These measures were strongly correlated [20], particularly, the last two (0.79), which was reflected by a calculation of the eigenvalues for the correlation matrix, which identified the number of effective traits as three. Secondly, for the mixed effects logistic regression analyses investigating categorical clinical group as outcomes, the number of effective traits were two, since we performed separate logistic regressions comparing cases versus controls and relatives versus controls. We present all uncorrected pvalues throughout the paper. However, interpretation of what constitute significant findings was based exclusively on the multiple-testing adjusted pvalue thresholds of 0.017 (0.05/3 effective traits) for cognition and of 0.025 (0.05/2 effective traits) for clinical group outcomes. All tests performed were two sided and the R-code used is available upon request from the corresponding author. Results Findings from the systematic review and metaanalysis of the literature The literature search returned 905 results. Screening of titles and abstracts revealed 13 papers that were assessed 5310 J. H. Thygesen et al. for eligibility. See PRISMA diagram and details on the literature search in Supplementary Fig. 2. Eleven primary papers of similar quality were found to examine the association between CNV burden and intelligence [30,31,35,43–47,66,67] (Supplementary Tables 1 and 4). Ten studies, with a total of 18,847 participants, provided the required data to conduct a meta-analysis and were included in the random effects meta-analysis [30,35,43–45,66,67]. Forest plots for analyses of length of deletions (N=18,658) and length of duplications (N=18,580) are displayed in Fig. 1, additional forest plots can be found in Supplementary Fig. 3. None of the meta-analyses showed evidence for an association between their measure of CNV burden and IQ. The pooled correlations between IQ and length of deletions or length of duplications were −0.04 (CI =−0.07, −0.01) and −0.002 (CI =−0.02, 0.02) respectively. Cochran’sQ- statistic revealed evidence for between-study heterogeneity in the correlation between length of deletions and IQ (χ2= 21.56, P=0.0058) and number of deletions and IQ (χ2= 33.77, P< 0.0001), number of duplications (χ2=11.06, P=0.0114). There was no evidence for study heterogeneity for the other measures. One study included in the systematic review did not provide data suitable for the meta-analysis. However, its findings are consistent with the meta-analysis since Van Scheltinga et al. [46] found no association between their measures of CNV burden and IQ. Results from the Psychosis Endophenotypes International Consortium sample The full sample consisted of 5597 participants. One thousand, three hundred three participants had CNV data failing quality control due to one or more of the following reasons: 1107 participants had more than 300 CNVs, 478 had BAF drift of >0.01 and 400 had Log R ratio standard deviation of >0.5. In our study, 77% of samples passed stringent quality control criteria for CNV calling. This call rate is comparable to other CNV studies as exemplified by the latest Psychiatric Genomics Consortium CNV large mega-analysis reporting an overall 72% call rate across 43 primary studies [12]. The challenges with CNV calling from SNP microarrays, are well known [68], and is to a large extent technical in origin. SNP microarrays were not originally designed for CNV detection, and although it is possible as demonstrated by numerus publications, CNV detection from SNP arrays is sensitive to quality issues especially with regards to the intensity measures captured by the probes on the assay [68]. The 23% of samples excluded on quality control grounds did not differ significantly from those included in the study on key parameters including clinical group distribution and sex. The only significant difference was in age, where the excluded samples on average were 5 years younger (see full details in Supplementary Table 5). Age is included as a covariate in all analysis. Of the 4294 participants who passed quality control, 3426 were included in the CNV burden analysis as they had at least one cognitive measure and full information on the included covariates. There were no significant differences in CNV burden measures between samples with and without cognitive data available (see Supplementary Table 6). An additional two participants that failed the CNV burden quality control were identified from two independent blind visual inspections as carriers of schizophrenia-associated CNVs, and included in that analysis (N=3428). This sample included 769 patients with psychotic disorders (576 with schizophrenia (74.9%), 89 with bipolar disorder Fig. 1 Forest plots for the meta-analyses investigating length of deletions (N=18,658) and length of duplications (N=18,580). For additional plots see Supplementary Fig. 3. Genetic copy number variants, cognition and psychosis: a meta-analysis and a family study 5311 (11.6%) and 104 with other psychoses (13.5%), 646 unaffected relatives and 2013 healthy controls (see Table 1). In our sample, we identified 29 participants who carried one schizophrenia-associated CNV each (see loci in Supplementary Table 3). Table 2shows the analyses of schizophrenia-associated CNVs and cognition, adjusted for age, gender, clinical group, centre and genetic relatedness. We found evidence of an association between schizophrenia-associated CNVs and RAVLT idiate (regression coefficient =−8.0, 95% CI =−13.3 to −2.6, P=0.0036) and delayed (regression coefficient =−3.3, 95% CI =−5.8, −0.7, P=0.0115) recall. This indicates that participants with a schizophrenia-associated CNV had a mean RAVLT immediate recall score that was 8.0% lower than non-carriers, as well as a mean RAVLT delayed recall score that was 3.3% lower. We also found suggestive evidence that carriers of a schizophrenia-associated CNV had poorer scores for block design than non-carriers (mean difference =−10.0, 95% CI =−19.2 to −0.9, P=0.031) but only at the uncorrected level of significance. As a sensitivity analysis we performed the same associations using only patients with a schizophrenia diagnosis, their relatives and healthy controls, see Supplementary Table 7. In that analysis the association between CNV carrier status and RAVLT immediate recall remained significant (P= 0.004), and we observed a weaker association with RAVLT delayed recall score at the uncorrected significance level (P=0.025). Table 1 Sample description. Characteristic CNV carrier (N) Noncarrier (N) Total (N) CNV carrier (%) Clinical group Controls 16 1997 2013 0.8% Relatives 2 644 646 0.3% Patients 11 758 769 1.5% Centre Munich 3 949 952 0.3% Perth 6 548 554 1.1% London 2 478 480 0.5% Maastricht 4 395 399 1.0% Amsterdam 4 329 333 1.2% Utrecht 5 304 309 1.6% Groningen 5 310 315 1.6% Pamplona 0 44 44 0.0% Edinburgh 0 42 42 0.0% Sex Males 16 1743 1759 0.9% Females 13 1656 1669 0.8% Characteristic CNV carrier (N) CNV carrier (mean/SD) Noncarrier (N) Non-carrier (mean/SD) Total (N) Total (mean/SD) Age 29 39.9 (15.6) 3399 43.7 (15.9) 3428 43.7 (15.9) Cognitive performance Block design 21 50.6 (24.9) 2726 57.4 (23.6) 2747 57.4 (42.0) Digit span 5 60.5 (4.5) 1265 50.7 (15.1) 1270 50.7 (9.4) RAVLT immediate Recall 24 46.5 (16.6) 1882 54.7 (14.5) 1906 54.6 (14.9) RAVLT delayed recall 24 14.6 (7.8) 1865 17.7 (6.9) 1889 17.7 (6.6) Demographic information for participants with data on schizophrenia-associated CNVs. Two participants with schizophrenia-associated CNV were identified in the samples failing QC for CNV burden, these were included only in the analysis of schizophrenia-associated CNVs; thus giving a sample of 3428 participants for the analysis of schizophrenia-associated CNVs and 3426 for the analysis of CNV burden. CNV carriers refer to individuals who were identified as carrying a known schizophrenia-associated CNV (see Supplementary Table 3). RAVLT Rey Auditory Verbal Learning Test, SD standard deviation. 5312 J. H. Thygesen et al. Our CNV burden analysis showed that the mean total length of DNA affected by deletions and duplications was 117.8 and 219.1 kb, respectively and on average 1.8 genes per participant were affected by these CNVs. We did not find evidence for an increased total CNV burden measured as length of DNA affected by CNVs, or as number of genes affected by CNVs, in people with psychosis compared to healthy controls. Since Marshall et al. [12] found that the association between CNV burden and schizophrenia was more significant when indexed as number of genes affected, we focused our subsequent analyses on this measure. For analyses using CNV length see Supplementary Tables 8 and 9. We found no evidence for an association between the four cognitive measures and any of the three CNV burden measures, see Table 2.Stratified analyses were also conducted by group (Supplementary Table 10), which showed an association between schizophrenia-associated CNVs and RAVLT immediate recall (regression coefficient =−11.8, 95% CI = −20.2 to −3.4, P=0.006), and a weaker association with the RAVLT delayed recall (P=0.02) in the patient group but only at the uncorrected level of significance. There was no other evidence for an association between the various cognitive measures and burden when examining the groups separately. The mixed effects logistic regression suggests that there was no evidence in the difference of having a schizophrenia-associated CNV amongst patients, relatives and controls. Similarly for CNV burden, the number of genes affected by large CNVs did not differ between the three clinical groups. See Supplementary Table 11 for details of CNV comparisons between clinical groups. Cognitive performance was impaired in patients compared with controls for all variables examined, as expected, and was worse in relatives than controls for block design and digit span. See Supplementary Table 12 for adjusted analyses. As a follow-up analysis we examined five loci found to protect carriers from developing schizophrenia [12,69]. We identified 41 carriers (22q11.21.dup (N=1), 7q11.21.del (N=21), 7q11.21.dup (N=7), 13q12.11.dup (N=7) and Xq28.dup (N=5)), but found no significant association Table 2 Associations between schizophrenia-associated CNVs and CNV burden with cognition. Predictor Cognitive measure Participants (CNV Carrier) Regression coefficient 95% CI Sig. Schizophreniaassociated CNVs Block design 2747 (21) −10.1 −19.2, −0.9 0.031 Digit span –––– RAVLT immediate Recall 1906 (24) −8.0 −13.3, −2.6 0.0036 RAVLT delayed recall 1889 (24) −3.3 −5.8, −0.7 0.0115 Genes affected by all CNVs (burden) Block design 2747 −0.1 −0.3, 0.1 0.343 Digit span 1270 0.03 −0.2, 0.2 0.775 RAVLT immediate recall 1906 −0.03 −0.2, 0.1 0.634 RAVLT delayed recall 1889 −0.01 −0.1, 0.1 0.854 Genes affected by deletions (burden) Block design 2747 −0.4 −0.8, 0.01 0.056 Digit span 1270 −0.1 −0.6, 0.3 0.544 RAVLT immediate recall 1906 −0.2 −0.5, −0.1 0.0119 RAVLT delayed recall 1889 −0.1 −0.2, −0.01 0.0661 Genes affected by duplications (burden) Block design 2747 −0.001 −0.2, 0.2 0.992 Digit span 1270 0.06 −0.1, 0.3 0.564 RAVLT immediate recall 1906 0.02 −0.1, 0.2 0.774 RAVLT delayed recall 1889 0.03 −0.04, 0.1 0.412 Figures in bold are below the multiple testing adjusted p-value threshold of 0.017. Associations between known schizophrenia-associated CNVs and CNV burden with cognitive performance. For the schizophrenia-associated CNV analysis digit span was not examined as fewer than ten CNV carriers had available data. CNV burden was measured as number of genes affected by CNVs larger than >200 kb, with <1% frequency. All analyses are adjusted for the covariates age, sex, clinical group, centre and genetic relatedness (kinship matrix). Genetic copy number variants, cognition and psychosis: a meta-analysis and a family study 5313 between performance in cognitive tests and carrier status, or in number of carriers between patients, relatives or controls. Discussion This study aimed to investigate: (1) the influence of CNV burden on general cognitive ability (IQ) based on a systematic review and meta-analysis of the literature and (2) the influence of schizophrenia-associated CNVs and CNV burden on specific cognitive skills in our family study from the PEIC. The meta-analysis of published studies found no associations between any CNV burden measures and overall IQ. The PEIC sample revealed that carriers of specific schizophrenia-associated CNVs had clear impairments in immediate and delayed verbal recall. Verbal memory performance has been found to index cortical thinning in medial temporal and prefrontal regions in schizophrenia [70,71] and has been found to be a cognitive predictor of outcome in schizophrenia and first episode psychosis [72,73] supporting its role as a plausible endophenotype in psychosis. We also found suggestive evidence that carriers of schizophrenia-associated CNVs perform worse in block design although further investigation is needed to verify the link between schizophrenia-associated CNVs and this and similar measures of perceptual reasoning. To the best of our knowledge, this is the first and only meta-analysis investigating the associations between CNV burden and IQ. We found ten relevant studies with a total of 18,847 participants. Evidence suggests that larger and rarer CNVs are more likely to be pathogenic [17]. Therefore, we hypothesised that the larger and rarer CNVs would have greater effects on intelligence. However, the primary studies included in our systematic review and meta-analysis did not always follow this pattern. Indeed the four studies, out of nine, that did find association with deletion length and intelligence included CNVs down to 100 kb or smaller. Two of these studies reported on fewer than 80 participants each. The two remaining studies [52,67] were among the largest available and performed a rigorous quality control. They both found evidence that all of their rare deletion burden measures (CNV length and number of genes affected) resulted in lower IQ. We believe that at least part of the reason why none of the metaanalyses we performed found association between any of the measures of CNV burden and intelligence, is due to the heterogeneity in the methodology and CNV criteria of the studies conducted in this area, and it seems the field is still in need of a more consistent and stringent way to analyse CNV burden. Along with Gialluisi et al. [48] our study is one of the first to report on associations between measures of large and rare CNV burden and performance on specific cognitive tests, and to our knowledge the first to examine this in both patients with psychosis and their unaffected relatives. The main limitation of our family study was the modest sample size, which limited our power to detect very small effects. CNV detection is known to be prone to false positive calls [74,75], therefore a strength is the thoroughness taken in calling CNVs, as all calls included in our burden analysis were visually inspected by two researchers blind to all clinical data in order to ensure their accuracy. We discarded 72% of the CNVs predicted by PennCNV after visual inspection, which shows the importance of such checks. In our sample, we found no evidence of selection bias in relation to quality control or to availability of cognitive data. Thus, the percentage of individuals passing CNV quality control was similar between patients, relatives and controls and comparable to that of Marshall et al. [12]. Furthermore, the CNV burden did not differ between individuals with and without cognitive data available. Our findings suggest that CNVs have greater effects on specific aspects of cognition rather than on general intelligence. Our family study investigated the burden of large (>200 kb), rare (<1% frequency) CNVs, which in general are expected to have greater phenotypic effects [5,12,76] than the smaller and more frequent CNVs included in many of the studies in the meta-analysis. Furthermore, the genetic contribution to IQ increases with age [77,78], and in the meta-analysis there was substantial diversity in ages and clinical group examined across the studies. Finally, the relationship between CNV burden and intelligence may vary substantially with clinical status. Indeed, the majority of studies of the meta-analysis that included patients with neuropsychiatric conditions [30,31,44,66] found associations between at least one measure of CNV burden and intelligence, whereas the majority of studies examining healthy participants [35,45,47] did not. However, our family study included 22.5% of participants with psychosis, compared to 0.8% participants in the meta-analysis, and still did not find an association with any burden measures. Current available studies, including ours, may still be underpowered to detect a small effect of CNV burden on cognition. Furthermore, it is also important to consider the crudity of the current CNV burden measures definedbysize and frequency criteria. As our understanding of the phenotypic effects of CNV improves, more specific burden measures targeting neurodevelopment and brain diseases will emerge. Limited power in our PEIC family study is likely to explain why we did not replicate the association between poorer digitspan and schizophrenia-associated CNVs, as reported by Kendall et al. [37] in the much larger UK Biobank study. The RAVLT immediate recall and delayed recall, which measures working memory and long-term memory respectively, are robust endophenotypes of schizophrenia [27]. Thus, as hypothesised, we found they were impaired amongst the carriers of known schizophrenia-associated CNVs in the PEIC family based sample. The existing literature [37,67,79] as well as our data provide consistent 5314 J. H. Thygesen et al. evidence that the carriers of specific CNVs that increase schizophrenia risk have cognitive impairments. It is widely agreed that a better understanding of the genetics of psychosis is essential for developing new diagnostic and therapeutic interventions. Animal and cellular models will provide essential evidence to understand the mechanisms of the implicated genetic loci, but are only available for a few CNVs. Studying endophenotypes in the human in vivo is non-invasive and one of the best tools available to elucidate the role and mechanisms of genetic variants that increase the risk of developing neuropsychiatric disorders. Acknowledgements We would like to thank all the patients, relatives and controls who took part in this research, as well as the clinical staff who facilitated their involvement. This work was supported by the Medical Research Council (G0901310) and the Wellcome Trust (grants 085475/B/08/Z, 085475/Z/08/Z). We thank the UCL Computer Science Cluster team for their excellent IT support. This study was supported by the NIHR Biomedical Research Centre at University College London Hospitals NHS Foundation Trust and University College London and by the NIHR Biomedical Research Centre for Mental Health at the South London and Maudsley NHS Foundation Trust at King’s College London. Further support to EB: Mental Health Research UK’s John Grace QC award, BMA Margaret Temple grants 2016 and 2006, MRC—Korean Health Industry Development Institute Partnering Award (MC_PC_16014), MRC New Investigator Award and a MRC Centenary Award (G0901310), National Institute of Health Research UK post-doctoral fellowship, the Psychiatry Research Trust, the Schizophrenia Research Fund, the Brain and Behaviour Research foundation’s NARSAD Young Investigator Awards 2005, 2008, Wellcome Trust Research Training Fellowship, the NIHR Biomedical Research Centre at UCLH, and the NIHR Biomedical Research Centre for Mental Health at the South London and Maudsley NHS Foundation Trust and Institute of Psychiatry King’s College London. Further support to co-authors: The Brain and Behaviour Research foundation’s (NARSAD’s) Young Investigator Award (Grant 22604, awarded to CI). The BMA Margaret Temple grant 2016 to JT. A 2014 European Research Council Marie Curie award to A Díez-Revuelta. HI has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 747429. A Medical Research Council doctoral studentship to JH-S, IA-Z and AB. A Mental Health Research UK studentship to RM. VB is supported by a Wellcome Trust Seed Award in Science (200589/Z/16/Z). FWO Senior Clinical Fellowship to RvW. The infrastructure for the GROUP consortium is funded through the Geestkracht programme of the Dutch Health Research Council (ZON-MW, grant number 10-000-1001), and matching funds from participating pharmaceutical companies (Lundbeck, AstraZeneca, Eli Lilly, Janssen Cilag) and universities and mental health care organisations (Amsterdam: Academic Psychiatric Centre of the Academic Medical Centre and the mental health institutions: GGZ Ingeest, Arkin, Dijk en Duin, GGZ Rivierduinen, Erasmus Medical Centre, GGZ Noord Holland Noord. Groningen: University Medical Centre Groningen and the mental health institutions: Lentis, GGZ Friesland, GGZ Drenthe, Dimence, Mediant, GGNet Warnsveld, Yulius Dordrecht and Parnassia psycho-medical centre The Hague. Maastricht: Maastricht University Medical Centre and the mental health institutions: GGZ Eindhoven en De Kempen, GGZ Breburg, GGZ Oost-Brabant, Vincent van Gogh voor Geestelijke Gezondheid, Mondriaan, Virenze riagg, Zuyderland GGZ, MET ggz, Universitair Centrum Sint-Jozef Kortenberg, CAPRI University of Antwerp, PC Ziekeren Sint-Truiden, PZ Sancta Maria Sint-Truiden, GGZ Overpelt, OPZ Rekem. Utrecht: University Medical Centre Utrecht and the mental health institutions Altrecht, GGZ Centraal and Delta). The Santander cohort was supported by Instituto de Salud Carlos III (PI020499, PI050427, PI060507), SENY Fundació (CI 2005-0308007), Fundacion Ramón Areces and Fundacion Marqués de Valdecilla (API07/011, API10/13). We thank Valdecilla Biobank for providing the biological PAFIP samples and associated data included in this study and for its help in the technical execution of this work; we also thank IDIVAL Neuroimaging Unit for its help in the acquisition and processing of imaging PAFIP data. Membership of Wellcome Trust Case Control Consortium 2 (WTCCC2) Management Committee Peter Donnelly (Chair)1,2, Ines Barroso (Deputy Chair)3, Jenefer M Blackwell4,5, Elvira Bramon6, Matthew A Brown7, Juan P Casas8, Aiden Corvin9, Panos Deloukas3, Audrey Duncanson10, Janusz Jankowski11, Hugh S Markus12, Christopher G Mathew13, Colin NA Palmer14, Robert Plomin15, Anna Rautanen1, Stephen J Sawcer16, Richard C Trembath13, Ananth C Viswanathan17, Nicholas W Wood18 Data and Analysis Group Chris C A Spencer1, Gavin Band1, Céline Bellenguez1, Colin Freeman1, Garrett Hellenthal1, Eleni Giannoulatou1, Matti Pirinen1, Richard Pearson1, Amy Strange1, Zhan Su1, Damjan Vukcevic1, Peter Donnelly1,2 DNA, Genotyping, Data QC and Informatics Group Cordelia Langford3, Sarah E Hunt3, Sarah Edkins3, Rhian Gwilliam3, Hannah Blackburn3, Suzannah J Bumpstead3, Serge Dronov3, Matthew Gillman3, Emma Gray3, Naomi Hammond3, Alagurevathi Jayakumar3, Owen T McCann3, Jennifer Liddle3, Simon C Potter3, Radhi Ravindrarajah3, Michelle Ricketts3, Avazeh Tashakkori- Ghanbaria3, Matthew Waller3, Paul Weston3, Sara Widaa3, Pamela Whittaker3, Ines Barroso3, Panos Deloukas3. Publications Committee Christopher G Mathew (Chair)13, Jenefer M Blackwell4,5, Matthew A Brown7, Aiden Corvin9, Mark I McCarthy19, Chris C A Spencer1 Affiliations of all WTCCC2 members: 1Wellcome Trust Centre for Human Genetics, Roosevelt Drive, Oxford OX3 7BN, UK; 2Department of Statistics, University of Oxford, Oxford OX1 3TG, UK; 3Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, UK; 4Telethon Institute for Child Health Research, Centre for Child Health Research, University of Western Australia, 100 Roberts Road, Subiaco, Western Australia 6008; 5Cambridge Institute for Medical Research, University of Cambridge School of Clinical Medicine, Cambridge CB2 0XY, UK; 6University College London, UCL Mental Health Sciences Unit and UCL Institute of Cognitive Neurosciences, London W1W 7EJ, UK; 7Diamantina Institute of Cancer, Immunology and Metabolic Medicine, Princess Alexandra Hospital, University of Queensland, Brisbane, Queensland, Australia; 8Department of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London WC1E 7HT and Department of Epidemiology and Public Health, University College London WC1E 6BT, UK; 9Neuropsychiatric Genetics Research Group, Institute of Molecular Medicine, Trinity College Dublin, Dublin 2, Ireland; 10Molecular and Physiological Sciences, The Wellcome Trust, London NW1 2BE; 11Centre for Digestive Diseases, Queen Mary University of London, London E1 2AD, UK and Digestive Diseases Centre, Leicester Royal Infirmary, Leicester LE7 7HH, UK and Department of Clinical Pharmacology, Old Road Campus, University of Oxford, Oxford OX3 7DQ, UK; 12Clinical Neurosciences, St George's University of London, London SW17 0RE; 13King’s College London Department of Medical and Molecular Genetics, School of Medicine, Guy’s Hospital, London SE1 9RT, UK; 14Biomedical Research Centre, Ninewells Hospital and Medical School, Dundee DD1 9SY, Genetic copy number variants, cognition and psychosis: a meta-analysis and a family study 5315