Impact of IQ on the diagnostic yield of chromosomal microarray in a community sample of adults with schizophrenia

Background Schizophrenia is a severe psychiatric disorder associated with IQ deficits. Rare copy number variations (CNVs) have been established to play an important role in the etiology of schizophrenia. Several of the large rare CNVs associated with schizophrenia have been shown to negatively affect IQ in population-based controls where no major neuropsychiatric disorder is reported. The aim of this study was to examine the diagnostic yield of microarray testing and the functional impact of genome-wide rare CNVs in a community ascertained cohort of adults with schizophrenia and low (< 85) or average (≥ 85) IQ. Methods We recruited 546 adults of European ancestry with schizophrenia from six community psychiatric clinics in Canada. Each individual was assigned to the low or average IQ group based on standardized tests and/or educational attainment. We used rigorous methods to detect genome-wide rare CNVs from high-resolution microarray data. We compared the burden of rare CNVs classified as pathogenic or as a variant of unknown significance (VUS) between each of the IQ groups and the genome-wide burden and functional impact of rare CNVs after excluding individuals with a pathogenic CNV. Results There were 39/546 (7.1%; 95% confidence interval [CI] = 5.2–9.7%) schizophrenia participants with at least one pathogenic CNV detected, significantly more of whom were from the low IQ group (odds ratio [OR] = 5.01 [2.28–11.03], p = 0.0001). Secondary analyses revealed that individuals with schizophrenia and average IQ had the lowest yield of pathogenic CNVs (n = 9/325; 2.8%), followed by those with borderline intellectual functioning (n = 9/130; 6.9%), non-verbal learning disability (n = 6/29; 20.7%), and co-morbid intellectual disability (n = 15/62; 24.2%). There was no significant difference in the burden of rare CNVs classified as a VUS between any of the IQ subgroups. There was a significantly (p=0.002) increased burden of rare genic duplications in individuals with schizophrenia and low IQ that persisted after excluding individuals with a pathogenic CNV. Conclusions Using high-resolution microarrays we were able to demonstrate for the first time that the burden of pathogenic CNVs in schizophrenia differs significantly between IQ subgroups. The results of this study have implications for clinical practice and may help inform future rare variant studies of schizophrenia using next-generation sequencing technologies. Electronic supplementary material The online version of this article (doi:10.1186/s13073-017-0488-z) contains supplementary material, which is available to authorized users.


Background
Schizophrenia is a severe psychiatric disorder associated with significant impairments in cognitive functioning [1]. On average, full scale IQ (FSIQ) is 7-8 points lower in cohorts with schizophrenia compared to general population norms [2] and the risk for schizophrenia has been shown to increase by 3.8% per 1-point decrease in FSIQ [3,4]. However, this risk appears to be greatest for individuals with FSIQ < 85, and for those with a significantly lower performance IQ (PIQ) than verbal IQ (VIQ) (i.e.~7-point difference or greater in the two main components of FSIQ) [4][5][6]. More extreme VIQ > PIQ discrepancies (i.e. ≥ 15 points) are clinically relevant and represent a neuropsychological hallmark of non-verbal learning disability (NVLD), a condition characterized by deficits in visualspatial perception, complex psychomotor skills, nonverbal problem solving, arithmetic, and social judgment [7,8]. The prevalence of schizophrenia in individuals with intellectual disability (ID; generally, IQ < 70) is threefold to fivefold higher than the general population prevalence of 1% [3,9]. Taken together, these data suggest that the underlying genetic mechanisms that predispose individuals to schizophrenia may be stronger in those with low FSIQ, particularly low PIQ, than in those with higher IQ. Given that the IQ deficits in schizophrenia are associated with functional outcome [1], further study of genetic risk variants for schizophrenia in the context of the intellectual profile appears warranted.
Rare copy number variations (CNVs) have been identified to play an important role in the etiology of schizophrenia and developmental disability and/or ID (DD/ ID) [10,11]. Several large rare CNVs, including deletions at 2p16.3 overlapping NRXN1, 15q13.3 (BP4-BP5) deletions, and 16p11.2 deletions/duplications, have been identified in schizophrenia and DD/ID [12][13][14]. Additionally, CNVs associated with schizophrenia have been shown to negatively affect IQ in population-based controls without any major neuropsychiatric disorder [15]. The widespread use of clinical microarray testing in DD/ ID has established the yield of pathogenic CNVs to be 15-20% [16]. In contrast, there have been significantly fewer diagnostic yield studies in schizophrenia [10,17], possibly due to the lack of guidelines endorsing routine clinical microarray testing in this complex adult-onset condition [18]. Since most rare CNV studies of schizophrenia do not report IQ and/or have excluded participants with co-morbid ID [13,19], the yield of pathogenic CNVs and the underlying genetic architecture of schizophrenia in the context of low IQ (schizophrenia-LIQ) remains unknown. Further, there have been no studies examining the genome-wide burden and/or functional impact of rare CNVs on schizophrenia while taking into account IQ, and after removing those CNVs that are deemed pathogenic.
Identifying sub-populations of individuals with schizophrenia who may be at an increased risk for a clinically reportable CNV, classified as pathogenic or a variant of unknown significance, would be useful for clinical practice. The primary aims of this study were twofold: (1) to compare the genome-wide burden of clinically reportable CNVs between individuals with schizophrenia-LIQ and schizophrenia-average IQ; and (2) to compare the genome-wide burden and functional impact of rare CNVs, beyond those that are currently deemed pathogenic, between individuals with schizophrenia-LIQ and average to higher IQ. Secondary analyses were aimed at identifying the yield of clinically reportable CNVs in schizophrenia across a wider range of IQ groups, including for those with a NVLD.

Schizophrenia sample collection and ascertainment
We recruited 688 adults who met the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, diagnostic criteria for schizophrenia or schizoaffective disorder. Our detailed ascertainment strategy is described elsewhere [10]; however, it should be noted that the majority of the individuals recruited were chronically ill and therefore unlikely to include individuals in the first onset of illness whose diagnosis may change over time. There were 644 participants ascertained from six community mental health clinics across Central and Eastern Canada. In order to increase the number of individuals with schizophrenia at the lower end of the IQ spectrum we recruited an additional 44 participants with schizophrenia and ID from two outpatient mental health clinics that specialize in treating adults with a dual diagnosis (ID and a psychiatric disorder). However, of these 44 individuals, only 19 (43.2%) were included in the final cohort of 546 unrelated participants of European ancestry with adequate IQ data. The CNV data for a subset of the individuals with schizophrenia (n = 459; 66.7%) were previously published [10], although without the associated IQ data. Consent was obtained from all participants and surrogate consent was provided by an individual with power of attorney or equivalent for health decisions for individuals deemed incapable of providing informed consent. This study was approved by local institutional research ethics boards at the Centre for Addiction and Mental Health, Saint John Horizon Health Network, Humber River Hospital, Queen Elizabeth Hospital, Hamilton Health Services, and Bethesda Services.

Clinical assessment of IQ level in individuals with schizophrenia
Similar to previous studies [20], we used a combination of previous IQ testing and educational attainment data to assign individuals with schizophrenia to an IQ subgroup.
We also performed a comprehensive screening interview with each individual and/or his or her relative(s) to obtain medical, developmental, educational, and psychiatric history in addition to detailed demographic information. We retrospectively reviewed the available lifetime medical and psychiatric records for all 688 participants, blind to CNV status, and recorded results from all previous IQ and clinical genetic testing. These previous genetic and IQ results were not known at the time of recruitment. There were 212 of 546 (38.8%) individuals in the final sample with IQ scores (n = 136; 19.8%) and/or descriptive IQ ranges (n = 76; 11.0%) available (collectively referred to as IQ scores in the remaining text), 202 (36.9%) of whom had age at testing and schizophrenia age at onset both available. The majority of these IQ scores (n = 164/202; 81.2%) were obtained during the five years preceding the first onset of psychotic illness or within the 15 years after onset. Eighteen (8.9%) individuals had IQ testing completed more than five years before the first onset of psychotic illness and 20 (9.9%) had testing completed more than 15 years after onset. Individuals with IQ data had to be stable enough (e.g. with respect to psychotic symptoms) to be able to complete standardized IQ testing. There were no data on antipsychotic treatment at testing, but such treatment is unlikely to have affected IQ results [21].
We assigned individuals to the schizophrenia-LIQ or the schizophrenia-average IQ group if they had an IQ score of < 85 or ≥ 85, or an estimated IQ of borderline/ ID or average range, respectively. For secondary analyses, individuals in the schizophrenia-LIQ group were divided into borderline intellectual functioning (IQ 71-85) or ID (IQ ≤ 70) groups. Given that the risk for schizophrenia may be higher for individuals with a significant discrepancy between their PIQ and VIQ scores we assigned individuals meeting criteria for a NVLD (PIQ ≥ 15 points lower than VIQ; Additional file 1: Figure  S1) to a separate schizophrenia-NVLD category [6,7]. We also used educational attainment to assign participants to intellectual functioning groups. However, IQ scores were deemed the more accurate measure of intellectual ability when years of education appeared out of keeping with expectations and functioning. Examples included individuals with IQ < 70 yet 12 years of education in a modified curriculum (assigned to schizophrenia-ID group) and individuals with IQ of 90 who left school to work after just eight years of education (assigned to schizophrenia-average IQ group).
In the absence of IQ scores we used educational attainment, which has a 0.6-0.7 correlation with FSIQ in the general population [22] and/or additional clinical data to assign individuals to each group as follows: the schizophrenia-LIQ group comprised individuals with a history of special education and/or had ID noted repeatedly throughout the medical records (estimated mild/moderate ID) and individuals who had 8-11 years of formal education with reported difficulties in school (e.g. repeated grades, enrolled in general courses in high school; estimated borderline intellectual functioning) [22,23]. Years of education are not informative for individuals with schizophrenia-ID given that the majority of individuals are enrolled in special education and/or had modified academic curriculums. Individuals who had completed ≥ 12 years of education (graduated high school), had no reported difficulties in school, and had not repeated any grades were assigned to the schizophrenia-average IQ group [22,23]. However, there were a number of scenarios in which our detailed clinical data led us to believe that an individual's formal educational attainment did not reflect their true cognitive abilities. For example, we assigned individuals to the schizophrenia-average IQ group if they left school early due to incarceration, vocational and/or family responsibilities, or early onset of psychotic symptoms if they were reported to have done well academically up until that point. All assessments of IQ and educational attainment were performed blind to CNV status.

CNV detection and annotation
High-quality genomic DNA was available for 540/546 (98.9%) participants and was submitted to the Centre for Applied Genomics in Toronto, Canada for genotyping using either the Affymetrix® Genome-Wide Human SNP array 6.0 or the CytoScan HD array. All samples met the Affymetrix quality control cut-offs. Similar to previous studies [10,24], we only included CNVs that were > 10 kb, identified by at least two CNV calling algorithms (two of ChAS, iPattern, or Genotyping Console for the CytoScan HD array and two of iPattern, Birdsuite or Genotyping Console for the Affymetrix 6.0 array), spanning ten consecutive array probes, and overlapping < 75% of segmental duplications. Over 90% of CNVs called using these criteria validate using a second laboratory method [24]. The CytoScan HD array has a higher resolution than the Affymetrix 6.0 array; however, 90.0% of deletions ≥ 25 kb and spanning 25 consecutive array probes and duplications ≥ 50 kb and spanning 50 consecutive array probes are concordant between the two microarrays [25]. There was no significant difference in the proportion of individuals from the schizophrenia-LIQ or the schizophrenia-average IQ group analyzed on the Affymetrix 6.0 and CytoScanHD array (χ 2 = 1.50, df = 1, p = 0.219). There were six (1.1% of 546) participants with 22q11.2 deletions included in the cohort who did not have Affymetrix 6.0 or CytoScan HD microarray data available and were therefore only included in the analyses comparing the burden of pathogenic CNVs.
We used 10,113 population-based controls (Additional file 1: Table S1) to adjudicate CNV rarity in the schizophrenia-LIQ and schizophrenia-average IQ groups. As before [10,24,26], we used a conservative definition of "rare," defined as CNVs found in < 0.1% of these 10,113 independent controls using a 50% overlap criterion. Further quality control methods included removing locus-specific batch effects (i.e. CNVs with identical coordinates and copy number state that were present in > 1% of the sample) and manually joining large CNVs that appeared to be fragmented [13]. All CNV coordinates are given using the Genome Reference Consortium February 2009 build of the human genome (GRCh37/hg19).

Assessment of ancestry and relatedness
We genotyped the 549,374 SNPs that are common to both the Affymetrix 6.0 and CytoScan HD arrays for participants using Birdseed v2 or Chromosomal Analysis Suite 3.1, respectively. Genotype data from 293,511 unlinked SNPs were used to estimate ancestry for the individuals with schizophrenia using PLINK [27]. Genotype data from 778 HapMap participants were used as a known reference for ancestry. Of the 688 individuals with schizophrenia in the original sample, 617 (89.6%) were identified to be of European descent. Pair-wise identity by descent analyses for individuals with high-resolution microarray data revealed that none of these participants were related to one another (all PI_HAT values were < 0.1). The unrelated individuals of European descent with schizophrenia who had sufficient IQ/educational data available to be categorized by intellect (n=546; 88.5% of 617) comprised the sample for this study.

Clinical adjudication of rare CNVs in schizophrenia participants
All rare (< 0.1%) exonic CNVs > 100 kb and all noncoding CNVs > 500 kb were assessed for clinical relevance by a trained cytogeneticist following the American College of Medical Genetics (ACMG) guidelines for CNV interpretation [28]. CNVs were classified according to the five standard ACMG categories: (1) pathogenic; (2) variant of unknown significance (VUS) likely pathogenic; (3) VUS; (4) VUS likely benign; and (5) benign. We considered CNVs classified as pathogenic or VUS-likely pathogenic to be pathogenic. CNVs defined as clinically reportable included those classified as pathogenic, VUSlikely pathogenic, and VUS. The yield of pathogenic CNVs, VUS, and clinically reportable CNVs (pathogenic and VUS combined) were calculated based on the proportion of individuals in the schizophrenia-LIQ vs the schizophrenia-average IQ group with at least one of these CNV types, regardless of size or chromosomal location.

Genome-wide CNV burden and statistical analyses
In our primary analyses, we tested the hypothesis that the genome-wide burden of clinically reportable CNVs was greater for participants with schizophrenia-LIQ than for those in the schizophrenia-average IQ group. In addition, after excluding individuals who were identified to have a pathogenic CNV (Table 1), which tend to be large and overlap many genes, we performed a logistic regression analysis to compare the total number, total length, and genic content of rare autosomal CNVs (all, deletions and duplications separately) > 10 kb between the schizophrenia-LIQ and schizophrenia-average IQ groups. Sex and genotyping platform were included as covariates. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated using R 3.3.1 software. All tests were twosided with p < 0.05 defined for statistical significance, and uncorrected given limited multiple testing.

Gene-set enrichment analysis
We performed a gene-set enrichment analysis in order to determine if the functional impact of rare autosomal CNVs differed between the schizophrenia-LIQ and schizophrenia-average IQ groups. We tested 17 genesets that were postulated to play a role in the pathogenesis of schizophrenia and/or DD/ID. These included 15 sets that were significantly enriched for deletions (n = 15) or duplications (n = 1) in a recent large-scale CNV study of schizophrenia [13]. Briefly, these included two sets containing genes that are predicted to be targets of FMR1 [29,30], three sets containing genes coding for members of N-methyl-D-aspartate receptors (NMDAR), neuronal activity-regulated cytoskeleton-associated protein, and components of the postsynaptic density (PSD) [31], and ten sets associated with neuronal function, synaptic components, and/or neurological/neurodevelopmental phenotypes in humans (n = 7) or mice (n = 3) [13]. We also included two sets that comprised genes that were overlapped significantly more often by deletions (n = 1) or duplications (n = 1) in a clinically ascertained cohort with DD/ID compared to controls [12]. Detailed descriptions of how these 17 gene-sets were compiled are outlined in Additional file 2.
The gene-set enrichment analysis used a logistic regression deviance test [31] [R/Bioconductor package cnvGSA: Gene Set Analysis of (Rare) Copy Number Variants (version 1.18.0)] (https://www.bioconductor.org/packages/release/bioc/html/cnvGSA.html) to evaluate if the number of genes overlapped by rare exonic deletions or duplications in each individual for each of the gene-sets (i.e. gene-set specific exonic burden) was predictive of the participant being a member of the schizophrenia-LIQ or the schizophrenia-average IQ group. We included sex, genotyping platform, and the total number of genes overlapped by rare CNVs per individual as covariates. Multiple-testing

Clinical features of the cohort
Of the 546 unrelated participants with schizophrenia of European descent, there were 325 (59.5%) assigned to the schizophrenia-average IQ group, 192 (35.2%) assigned to the schizophrenia-LIQ group, 130 (67.7%) with borderline intellectual functioning and 62 with mild (n = 57) or moderate (n = 5) ID, and 29 (5.3%) assigned to the schizophrenia-NVLD group. Total years of education were significantly lower (Mann-Whitney U = 6453.5, p < 0.0001) in the schizophrenia-borderline intellectual functioning group (median = 10; range = 5-16 years) compared to the schizophrenia-average IQ group (median = 12; range = 5-19 years) and not significantly different between the schizophrenia-average IQ and the schizophrenia-NVLD (median = 12; range = 7-18 years) groups (p = 0.385). Before involvement in this study, only seven (1.3%) individuals from the entire cohort had previously received clinical genetic testing, all of whom had been recruited from a specialized dual diagnosis clinic. These included six (9.6%) individuals from the schizophrenia-ID group and one (0.8%) individual from the schizophrenia-borderline intellectual functioning group. Further demographic and clinical data for the cohort are provided in Additional file 1: Table S2.

Total burden of clinically reportable CNVs
There were 78 CNVs classified as VUS in 70 (12.8%; 95% CI 10.2-16.0%) schizophrenia participants (Additional file 3), five of whom also had a pathogenic CNV. In contrast to the pathogenic CNV results, there was no significant difference (p = 0.243) in the prevalence of individuals with one or more VUS between the schizophrenia-LIQ group (n = 26/192; 23.5%) and schizophrenia-average IQ group (n = 33/325; 10.2%). Secondary analyses revealed that there was also no significant difference in the prevalence of participants with a VUS between any of the IQ subgroups (Fig. 1). Of the 78 CNVs classified as VUS (median size 723 kb; range 115 kb to 4.3 Mb), there were slightly more duplications (n = 51; 65.3%) than deletions (n = 27; 34.7%), but this difference was non-significant (p = 0.057). Taken together, there were 99 (18.1%; 95% CI 15.0-21.7%) schizophrenia individuals with a clinically reportable CNV (pathogenic and/or VUS). Together there were 14 (2.6%) participants with two or more clinically reportable CNVs, significantly more of whom were in the schizophrenia-ID (n = 5/62; 8.1%) group compared to the schizophrenia-average IQ (n = 4/325; 1.2%) group (OR = 8.06 [2.21-32.93], p = 0.0018). There was no significant Fig. 1 Yield of clinically reportable CNVs in schizophrenia by IQ group. The figure depicts the percentage of individuals with schizophrenia for each of the IQ groups with one or more pathogenic (defined as pathogenic or VUS-likely pathogenic) CNV (a) or one or more CNV classified as a VUS (b), determined using the ACMG guidelines for CNV interpretation [28]. Individuals with more than one clinically reportable CNV were only counted once. Schizophrenia participants were assigned to each of the IQ sub-groups using the methods described in the manuscript. Average average, IQ intelligence quotient, BL borderline intellectual functioning, ID intellectual disability, NVLD non-verbal learning disability, VUS variant of unknown significance. Asterisks above horizontal brackets represent the significance of the between-group comparisons: *p < 0.05, **p < 0.01, ***p < 0.001, NS = p > 0.1. All other comparisons were non-significant difference (p = 0.135) between the schizophrenia borderline intellectual functioning group (n = 5/130; 3.8%) and the schizophrenia-average IQ group.

Burden of genome-wide rare CNVs
We also attempted to determine if the genome-wide burden of rare CNV was higher in an expanded schizophrenia-LIQ group compared to the schizophrenia-average IQ group, after excluding the 39 individuals with a pathogenic CNV. Given that the prevalence of pathogenic CNVs was similar for participants with ID and a NVLD, we added the 23 individuals with a NVLD (and no pathogenic CNV) to the original schizophrenia-LIQ group for the remaining analyses. After controlling for sex and genotyping platform there was no significant difference in the genome-wide burden, total genomic length, or total number of genes overlapped by rare autosomal CNVs between the two groups (Additional file 1: Table S3). However, there were significantly more genic CNVs in the expanded schizophrenia-LIQ group compared to the schizophrenia-average IQ group (OR = 1.19 [1.01-1.41], p = 0.042), primarily driven by an increased burden of genic duplications (OR = 1.42 [1.14-1.81], p = 0.002); findings for genic deletions did not reach significance (p = 0.129) (Additional file 1: Table S3).

Gene-set enrichment analysis
After multiple-test correction (BH-FDR < 10% and p value < 0.05), we detected no gene-sets that were significantly enriched for rare autosomal deletions in the expanded schizophrenia-LIQ group (comprising individuals with schizophrenia and ID, borderline intellectual functioning, or NVLD) compared to the schizophrenia-average IQ group (data not shown). There was one gene-set, GO nervous system development, that was significantly (p = 0.013) enriched for rare duplications in the schizophrenia-LIQ group that had a BH-FDR of 0.22 (Table 2). To see if we could improve the FDR, we reduced the 17 gene-sets to only those six reported to have a FDR < 30% for rare genic duplications in the recent Psychiatric Genomics Consortium schizophrenia case control CNV study [13]. This resulted in an improved FDR (from 0.22 to 0.07) for the GO nervous system development gene-set. The GO nervous system development gene-set became nonsignificant (p = 0.074, FDR = 0.37) after excluding n = 39 participants with an pathogenic CNV ( Table 2).
There were 44 rare duplications in 35 individuals in the expanded schizophrenia-LIQ group and 29 rare duplications in 28 individuals in the schizophrenia-average IQ group contributing to the GO nervous system development gene-set result (Additional file 3). The duplications not currently classified as pathogenic or a VUS in the schizophrenia-LIQ individuals overlapped several interesting neuropsychiatric candidate genes, such as CNTN4, NDUFV2, and RCAN1 [47][48][49]. There were also duplications in two participants from the expanded schizophrenia-LIQ group that overlapped two genes (ARSA and EIF2B1) associated with leukodystrophy, a progressive disease that causes abnormal development of and/or destruction of the myelin sheath and can present in adulthood with symptoms similar to that of schizophrenia [50,51].

Discussion
This is the first study to examine the burden of clinically reportable CNVs in schizophrenia by IQ group. Our results revealed that 7.1% of schizophrenia individuals ascertained from a community outpatient setting may have a pathogenic CNV detected by genome-wide microarray. However, this diagnostic yield was not uniformly distributed across the cohort as there was a The CNVs comprised four pathogenic 16p11.2 duplications in four individuals, for which the contributing gene for this gene-set result was MAPK3 LIQ low IQ, GO gene ontology, CNV copy number variation, p statistical result when all rare autosomal CNVs are included; BH-FDR Benjamini-Hochberg false discovery rate, Inf infinity, NMDAR N-methyl-D-aspartate receptor components significant increase in the yield of pathogenic CNVs as IQ decreased (Fig. 1). We also demonstrated for the first time that the prevalence of pathogenic CNVs may be similar for individuals with schizophrenia-ID and schizophrenia-NVLD. Further, we identified an increased burden of rare genic autosomal duplications in the schizophrenia-LIQ compared to the schizophreniaaverage IQ group, a finding that was not attributable to large rare pathogenic CNVs.
The importance of clinical microarray testing in the dual diagnosis adult population In the current study, the highest yield of pathogenic CNVs (24.1%) was identified in individuals with schizophrenia and co-morbid ID. This yield was higher than that that reported for epilepsy (~5-10%) [52] or ASD (~10-15%) [53] alone, and comparable to that for DD/ ID (~15-20%) [16]. There have been few studies that have examined the burden of pathogenic CNVs in adults with a dual diagnosis (ID plus one or more additional neurodevelopmental and/or neuropsychiatric conditions), and even fewer that have focused specifically on schizophrenia-ID. A recent study identified a pathogenic CNV in nine of 72 (12.5%) individuals with ID and psychosis [54], about half of the yield reported in the current study. Such discrepancies in the diagnostic yield may be due to differences in ID severity between these two cohorts (data unavailable for the Wolfe et al. cohort) [54] and/or differences in the application of guidelines for CNV interpretation that are vulnerable to subjective evaluations of the evidence supporting the classification of a variant [55]. Further studies examining the diagnostic yield of microarray in the schizophrenia-ID population will be needed to clarify this. With respect to other neuropsychiatric conditions, the rate of rare de novo CNVs has been shown to be increase with decreasing IQ individuals with ASD [56] but no studies have yet formally examined the diagnostic yield of clinical microarray in the ASD-ID population. Conversely, a recent study investigating adults with ID and pediatric-onset epilepsy revealed that 16.0% (n = 23/ 143) of the cohort had a pathogenic CNV [57], a higher prevalence than that reported for just epilepsy [52]. Taken together, these data suggest that adults with a dual diagnosis should be prioritized for clinical microarray testing. However, it is important to note that before inclusion in our study fewer than 10% of the dual diagnosis schizophrenia-ID participants in our cohort had received any type of clinical genetic testing despite meeting criteria for routine CMA testing based on the Miller et al. 2010 clinical recommendations [16]. This suggests that additional efforts to increase widespread CMA testing in the adult DD/ID population, particularly for those with a dual diagnosis, are needed.

Increased burden of clinically relevant CNVs in individuals with schizophrenia-NVLD
To our knowledge, this is the first study to report on the burden of clinically reportable CNVs in individuals with schizophrenia-NVLD. Individuals with a NVLD demonstrate significant deficits in visual-spatial organization, motor coordination, and social perception and interaction, but retain relatively well-developed verbal skills [7]. Determining a diagnosis of NVLD relies heavily on formal IQ testing, allowing for the detection of significant VIQ > PIQ discrepancies that are difficult to ascertain clinically. Perhaps unsurprisingly, given the emphasis on verbal skills in formal education, there was no significant difference in the total years of education between the schizophrenia-NVLD and schizophreniaaverage IQ groups (Additional file 1: Table S2). Yet, the schizophrenia-NVLD participants had a significantly higher burden of pathogenic CNVs compared to the schizophrenia-average IQ group. Indeed, the burden of pathogenic CNVs in the schizophrenia-NVLD group was more comparable to that of the schizophrenia-ID group (20.7% vs 24.1%). Interestingly, all of the individuals with schizophrenia-NVLD and a pathogenic CNV had a PIQ < 85, yet VIQ for all but one individual was in the average to above average range (Additional file 1: Figure S1). This finding has important clinical relevance because individuals with schizophrenia-NVLD would not generally be considered for clinical microarray testing [16].
NVLD has been previously described as being associated with several structural variants, including those underlying Turner syndrome [58], Williams syndrome [59], and the 22q11.2 deletion syndrome [60]. This study extends these findings to include other rare recurrent CNVs, including 15q11.2-q13.1 and 16p11.3 duplications. We also identified a novel 280-kb deletion at 3p26.1 overlapping genes SUMF1 and ITPR1 in a participant with schizophrenia-NVLD that has not been previously reported in the literature. ITPR1 encodes the inositol 1,4,5-triphosphate receptor that plays an important role in releasing Ca 2+ from the endoplasmic reticulum [61]. Interestingly, a recent whole-exome sequencing (WES) study identified 11 individuals with schizophrenia who had ultra-rare disrupting/damaging variants in ITPR1 [46]. There is also a 60-kb deletion overlapping the first four exons of ITPR1 in a case (nsv996226) with autistic behavior reported in the Clinical Genome Resource database (https://www.clinicalgenome.org). Heterozygous deletions and missense mutations in this gene have also been associated with adult onset spinocerebellar ataxia-15 (MIM 606658) and childhood onset spinocerebellar ataxia-29 (MIM 117360) [62,63]. Additionally, homozygous/compound heterozygous truncating mutations and heterozygous deletions and missense mutations in ITPR1 have been associated with Gillespie syndrome (MIM 206700), a disorder characterized by hypotonia, progressive hypoplasia, ataxia, and variable cognitive impairment with onset occurring within the first year of life [61]. Notably, many genes associated with spinocerebellar ataxia are also reported to play a possible role in schizophrenia, such as ATAXN1, ATAXN2, and ATAXN10 [64][65][66].
Delineating the genetic architecture of schizophrenia Data from this study demonstrate that the genetic architecture of schizophrenia-LIQ, and probably schizophrenia-NVLD, differs significantly from schizophrenia-average IQ, even after excluding large rare pathogenic CNVs that have a well-documented impact on cognition in the general population [15]. Despite having a relatively small sample size, we were able to identify an increased burden of rare exonic duplications in the schizophrenia-LIQ group that may overlap more genes involved in nervous system development ( Table 2). Larger sample sizes could help provide improved statistical support for these findings as well as potentially identify additional pathways relevant to schizophrenia-LIQ. Data from a recent WES study in schizophrenia identified a significantly increased burden of rare damaging variants in loss-of-function (LoF) intolerant genes associated with developmental disorders in individuals with schizophrenia-ID compared to those with schizophrenia and average IQ [20]. Interestingly, the burden of these damaging variants in LoF intolerant genes was also increased in schizophrenia-average IQ compared to controls, suggesting that they contribute to risk for developing schizophrenia but to a lesser degree than that for schizophrenia-ID [20].

Advantages and limitations
Significant strengths of our study included the robust CNV detection methods used and systematic application of established guidelines for CNV interpretation [28]. Also, the community-based ascertainment strategy more closely reflects the general schizophrenia population than a strictly hospital-based recruitment strategy, thus allowing for more generalizable diagnostic yield estimates. The comprehensive phenotyping protocol facilitated our ability to stratify the schizophrenia cohort by IQ. Yet, despite substantial efforts, the numbers of individuals with schizophrenia-ID remained relatively small. We therefore targeted specialized dual diagnosis clinics in order to increase recruitment at the lower end of the IQ range. This appeared to increase the yield of pathogenic CNVs detected in the borderline intellectual functioning group, possibly reflecting an ascertainment bias in which more severely affected and/or psychiatrically complex individuals may be referred to such specialized clinics. Efforts to identify more individuals from community settings are needed.
Another limitation of our study was that although there were formal IQ data for a substantial number of individuals, educational data were used to group the remainder of the cohort. While the correlation between FSIQ and educational attainment is high (0.6-0.7) [22], it is possible that we have misclassified some participants, particularly those with a learning disability (such as NVLD) that cannot be delineated by educational attainment alone. Further, given that there are some data to suggest that there is an IQ decline associated with schizophrenia illness onset [67] and that the majority of the IQ scores used in this study were obtained just before or sometime after the onset of illness, it is possible that some individuals were incorrectly assigned to an IQ group lower than their "premorbid" IQ score would have placed them. However, we based IQ group placement on multiple pieces of evidence, including IQ scores, educational attainment, functioning, and personal circumstances. This process mirrors the multiple factors taken into consideration when making a clinical diagnosis of intellectual functioning. Also, each IQ group covers a fairly wide range of IQ scores. We therefore believe that IQ group misclassification is likely to be low for this study.

Conclusions
Using high-resolution microarrays, the results indicate that the burden of pathogenic CNVs is significantly greater for individuals with schizophrenia and low IQ compared to those with normal to superior IQ. These data have important clinical and research implications, including demonstrating that participants with schizophrenia and low IQ should be prioritized for clinical microarray testing and highlighting the importance of taking IQ into account for the interpretation of future rare variant studies of schizophrenia. The data also suggest that individuals with schizophrenia-NVLD, comprising 5.3% of the sample, may also have an increased burden of pathogenic CNVs. Data from next-generation sequencing will allow for the detection of sequence variants and smaller structural variants that may help shed light on more specific mechanisms related to schizophrenia-LIQ. While we identified several potential candidate genes for schizophrenia-LIQ, larger samples will be required to provide sufficient statistical support for any given loci.

Additional files
Additional file 1: A word document that contains one figure that depicts the verbal and performance IQ scores for 29 individuals with schizophrenia and a NVLD ( Figure S1), and three tables, including: (1) a list of 10,113 population-based controls used to adjudicate CNV rarity in schizophrenia participants (Table S1); (2) the demographic and clinical information for 546 probands with schizophrenia of European ancestry (Table S2); and (3) the genome-wide burden of all rare autosomal CNVs > 10 kb between the expanded schizophrenia-LIQ and schizophrenia-average IQ groups (Table S3) Availability of data and materials All relevant microarray data have been deposited into the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) database and can be obtained using the accession number GSE106818. The rare (< 0.1% in population controls) CNV calls from this cohort are also available in Additional file 3.
Author's contributions CL and ASB designed the study. CL and GC recruited participants and collected the clinical data and samples. JW and KB provided the 44 additional individuals with schizophrenia-ID. CL and ASB collected the cognitive and educational data and performed all statistical analyses. AL, JW, CRM, and SWS performed the genomic analyses. AN, DJS, and MS adjudicated the CNV data for clinical relevance. CL and DM designed and performed the gene-set enrichment analysis. CL wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate This research project complies with the Declaration of Helsinki. Consent was obtained from all participants and surrogate consent was provided by an individual with power of attorney or equivalent for health decisions for individuals deemed incapable of providing informed consent. This study was approved by local institutional research ethics boards at the Centre for Addiction and Mental Health, Saint John Horizon Health Network, Humber River Hospital, Queen Elizabeth Hospital, Hamilton Health Services, and Bethesda Services.

Consent for publication Not applicable
Competing interests DM is an employee of Deep Genomics, Inc. The remaining authors declare that they have no competing interests.

Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.