Key PointsQuestion
What is the population-based prevalence and psychiatric risk associated with recurrent copy number variants (rCNVs) and how do they differ across outcomes and rCNVs?
Findings
In a genetic association study of the Danish population born in 1981-2008 and followed up until 2015, most rCNV risk profiles were similar across autism, attention-deficit/hyperactivity disorder, and schizophrenia. The magnitude of risk increase correlated with locus size and gene constraint but did not differ between deletions and duplications.
Meaning
Several rCNVs were more prevalent and conferred less risk of psychiatric disorders than estimated previously, which highlights the importance of population-based estimates for implementing genetics-based prediction in health care.
Importance
Recurrent copy number variants (rCNVs) have been associated with increased risk of psychiatric disorders in case-control studies, but their population-level impact is unknown.
Objective
To provide unbiased population-based estimates of prevalence and risk associated with psychiatric disorders for rCNVs and to compare risks across outcomes, rCNV dosage type (deletions or duplications), and locus features.
Design, Setting, and Participants
This genetic association study is an analysis of data from the Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) case-cohort sample of individuals born in Denmark in 1981-2008 and followed up until 2015, including (1) all individuals (n = 92 531) with a hospital discharge diagnosis of attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder, major depressive disorder (MDD), or schizophrenia spectrum disorder (SSD) and (2) a subcohort (n = 50 625) randomly drawn from the source population. Data were analyzed from January 2021 to August 2023.
Exposures
Carrier status of deletions and duplications at 27 autosomal rCNV loci was determined from neonatal blood samples genotyped on single-nucleotide variant microarrays.
Main Outcomes and Measures
Population-based rCNV prevalence was estimated with a survey model using finite population correction to account for oversampling of cases. Hazard ratio (HR) estimates and 95% CIs for psychiatric disorders were derived using weighted Cox proportional hazard models. Risks were compared across outcomes, dosage type, and locus features using generalized estimating equation models.
Results
A total of 3547 rCNVs were identified in 64 735 individuals assigned male at birth (53.8%) and 55 512 individuals assigned female at birth (46.2%) whose age at the end of follow-up ranged from 7.0 to 34.7 years (mean, 21.8 years). Most observed increases in rCNV-associated risk for ADHD, ASD, or SSD were moderate, and risk estimates were highly correlated across these disorders. Notable exceptions included high ASD-associated risk observed for Prader-Willi/Angelman syndrome duplications (HR, 20.8; 95% CI, 7.9-55). No rCNV was associated with increased MDD risk. Also, rCNV-associated risk was positively correlated with locus size and gene constraint but not with dosage type. Comparison with published case-control and community-based studies revealed a higher prevalence of deletions and lower associated increase in risk for several rCNVs in iPSYCH2015.
Conclusions and Relevance
This study found that several rCNVs were more prevalent and conferred less risk of psychiatric disorders than estimated previously. Most case-control studies overestimate rCNV-associated risk of psychiatric disorders, likely because of selection bias. In an era where genetics is increasingly being clinically applied, these results highlight the importance of population-based risk estimates for genetics-based predictions.
Recurrent copy number variants (rCNVs) arise through nonallelic homologous recombination mediated by low-copy repeat sequences and are maintained at stable population frequencies through high de novo mutation rates, despite often being associated with substantial pathogenicity and reduced reproductive fitness.1-3 Around 40 loci in the human genome harbor rCNVs conferring risk of intellectual disability, developmental delay, and congenital malformations.4 Many of these have also been associated with increased risk of neuropsychiatric disorders; autism spectrum disorder (ASD),1 attention-deficit/hyperactivity disorder (ADHD),5 and schizophrenia spectrum disorder (SSD),6-8 while evidence for bipolar disorder (BPD)9 and major depressive disorder (MDD)10 is limited.
Most rCNV studies have estimated disease-associated risk from case-control samples, often based on highly selected cases (eg, with long-term illness) and controls (eg, screened for lack of any family history of mental illness). While this design increases power for detecting associations between genetic exposures and disease outcomes, it yields exaggerated estimates of the prevalence of the exposure and the associated risk of the outcome. Studies involving community-based samples (eg, the UK Biobank [UKB]) benefit from larger sample sizes and a multitude of available study outcomes.10,11 However, the UKB suffers a “healthy volunteer” bias, which affects association estimates in genetic studies.12 Two health care system–based biobank studies from the US estimated population-based rCNV risk for a range of conditions but could not estimate risk for specific psychiatric disorders such as ADHD, ASD, and SSD because of small case samples.13,14 This lack of population-based risk estimates for rCNVs compromises implementation of genetically informed clinical care and genetic counselling.
The Lundbeck Foundation Initiative for Integrative Psychiatric Research 2012 (iPSYCH2012) case-cohort sample leverages Danish registers and biobanks with full population coverage to maximize study power for psychiatric disorders within a population-based setting.15 Our previous studies of rCNVs at 7 loci in iPSYCH2012 found that in many instances, the population-based prevalence was higher, and associated risk of psychiatric disorders lower, than typically reported in case-control studies.16,17
Here, we expand our population-based analyses by (1) systematically considering deletions and duplications across all rCNV loci, including many that have been associated with genomic disorders and clinical phenotypes11; (2) comparing the effects on risk across psychiatric outcomes and rCNVs; and (3) doing so with a substantially increased sample size (iPSYCH2015),18 affording unprecedented precision and power to characterize population-unbiased rCNV risk profiles with respect to psychiatric outcomes. The study follows the Strengthening the Reporting of Observational Studies in Epidemiology () reporting guideline on case-cohort studies.19,20
The study is based on the iPSYCH2015 case-cohort sample,18 which is an update and expansion of iPSYCH201215 and includes 140 116 individuals from 1 657 449 singletons born between May 1, 1981, and December 31, 2008, in Denmark, who were residents in Denmark at 1 year of age and have a mother registered in the Danish Civil Registration System.21 The case-cohort includes 2 components: first, the cases were all individuals with a clinical diagnosis of ADHD, ASD, BPD, MDD, or SSD until December 31, 2015, obtained from the Psychiatric Central Research Register.22 Second, the cohort was 50 615 individuals randomly drawn from the same source population, corresponding to 3% (see the eMethods in Supplement 1 for details).
In addition to the psychiatric diagnosis groups specifically targeted in the iPSYCH2015 case-cohort design, the same information (date of diagnosis) was obtained through Psychiatric Central Research Register22 and the Danish National Patient Registry23 for intellectual disability and epilepsy (eMethods in Supplement 1). For this study, the Danish Scientific Ethics Committee has, in accordance with Danish legislation, waived the need for informed consent in biomedical research based on existing biobanks.
DNA was extracted from neonatal blood spots retrieved from the Danish Neonatal Screening Biobank (DNSB).24 Genotyping of samples is described elsewhere.15,18 Extraction of B-allele frequency and probe intensities (through log-R-ratio; LRR) was done with GenomeStudio software (Illumina). Samples with a genotyping call rate less than 95% were excluded.18 Access to the data and its use for research purposes was granted by the Danish Scientific Ethics Committee, Danish Health Data Authority, Danish data protection agency, and DNSB Steering Committee.
CNV Detection and Quality Control
We selected rCNV loci based on CNVs studied in the UKB11 that are also listed as recurrent CNVs in ClinGen25 with some indication of either haploinsufficiency or triplosensitivity26 (eMethods and eTable 1 in Supplement 1). PennCNV27 calls from these loci were validated using software developed in-house.28 We excluded samples with an LRR-SD of 0.35 or greater , B-allele frequency drift of 0.005 or greater, and/or GC-wave factor of 0.02 or greater (eFigure 1 in Supplement 1). We assessed the quality of visual inspection of calls through intrarater and interrater reliability (eTable 2 in Supplement 1) and tested whether intrarater and interrater reliability and the fraction of calls with undetermined carrier status correlated with LRR-SD and/or locus, array, or dosage type (eTable 3 in Supplement 1). We then performed a series of analyses to verify that sensitivity in rCNV detection (the ability of PennCNV to detect true CNVs) and specificity of visual validation (our ability to correctly classify all PennCNV-detected calls as true or false) was not affected by the level of LRR-SD or by array type (eMethods, eFigure 2, and eTables 4-6 in Supplement 1).
We calculated population-based rCNV prevalence with 95% CIs using the R survey package, with finite population correction to account for oversampling of cases, and compared deletion-vs-duplication prevalence in iPSYCH2015, as well as iPSYCH2015-vs-UKB prevalence, with a Welch test assuming unequal variance and false discovery rate (FDR) to account for multiple testing (eMethods in Supplement 1). We used sex-stratified, weighted Cox proportional hazard (CPH) models from the R survival package to estimate the rCNV-associated risk of the psychiatric disorders targeted by the case-cohort design, with age at first hospital diagnosis as outcome and rCNV carrier status as exposure (eMethods in Supplement 1). Data were analyzed from January 2021 to August 2023.
To test for differences in rCNV-associated risk between pairs of diagnoses (ASD, ADHD, and SSD) we applied generalized estimating equation (GEE) models using the R glmtoolbox package, accounting for the possibility of being diagnosed with both disorders. Details, including of analyses of rCNV dosage type (deletion [del] vs duplication [dup]) and locus features, also involving use of GEE models, are provided in eMethods in Supplement 1.
As a sensitivity analysis, we performed a logistic regression analysis corresponding to each CPH analysis. Only individuals with the corresponding outcome or from the subcohort were included in each analysis, with the diagnosis as a binary outcome and sex, age at end of follow-up, and genotyping array as covariates (eMethods in Supplement 1).
We processed PennCNV27 calls through the QCtreeCNV pipeline28 and further quality control (eMethods in Supplement 1), yielding 3547 verified rCNV calls in 120 247 samples: 79 535 cases and 43 311 cohort samples, including an overlap (cases-in-cohort) of 2599 samples. Sample genotype success rate was 85.8% and did not differ across case, cohort, and cases-in-cohort groups (Fisher exact test, simulated P = .12). The sample included 64 735 male (53.8%) and 55 512 female (46.2%) individuals (recorded at birth). Age at end of follow-up ranged from 7.0 to 34.7 years (mean [SD], 21.8 [7.0] years). All participants were born in Denmark, but information on self-declared ethnicity was not available. Nine loci were excluded because of low rCNV prevalence (eTable 1 in Supplement 1).
Consistent with previous studies,4,11 the overall prevalence of duplications was higher than that of deletions (1.29% vs 0.99%, respectively; P = 2.2 ×ĉ10−5). However, the prevalence disparity differed substantially across loci (Figure 1A), with 6 and 2 loci showing significantly higher prevalence of duplications and deletions, respectively (eTable 7 in Supplement 1). When comparing rCNV prevalence with the UKB11 for the 17 loci with estimates in both studies, the overall population-based prevalence of deletions was higher in iPSYCH2015 (0.98% vs 0.70%, respectively; P = 7.7 ×ĉ10−12) whereas overall duplication prevalence did not differ significantly (1.27% vs 1.23%, respectively; P = .42). In the per-locus comparison, 9 deletions and 1 duplication had higher prevalence (FDR P < .05) in iPSYCH2015, 1 duplication was more prevalent in UKB, and 23 rCNVs showed no significant differences (Figure 1B and eTable 7 in Supplement 1).
Per-diagnosis results for the main diagnosis groups (ASD, ADHD, MDD, and SSD) are shown in Figure 2, while a per-locus version of all estimates, including for schizophrenia (International Statistical Classification of Diseases and Related Health Problems, Tenth Revision, code F20), BPD, and the broader diagnosis groups (any affective disorder and any iPSYCH2015 disorder) are provided in eFigure 3 in Supplement 1. Risk estimates exceeded 1.5 for a majority of rCNVs tested for ASD (18/34), ADHD (21/32), and SSD (19/31) but only for 6 of 32 tested for MDD (Figure 2). In fact, we observed no significantly increased HRs for MDD (Figure 2) or the broader category of any affective disorder (eTable 8 in Supplement 1). The greatest risk increase was observed for Prader-Willi/Angelman Syndrome (PWAS)-dup with ASD (HR, 20.7; 95% CI, 7.8-54.9) and SSD (HR, 9.1; 95% CI, 1.7-49.2) and for 17q12-del with ADHD (HR, 4.2; 95% CI, 1.2-13.9) (Figure 2).
We compared risk estimates for the main diagnosis groups with the largest available published case-control studies. The number of rCNVs with available estimates in those studies ranged from 11 for ASD1 to 13 for ADHD,5 20 for MDD,10 and 22 for SSD.6-8 Risk estimates were significantly lower in our study for 1q21.1-del, 15q11.2-del, 16p11.2-dup, 17p12-del, and 22q11.2-del for SSD; 16p11.2-del and 16p11.2-dup for ASD; and 22q11.2-del for ADHD (eFigure 4 and eTable 9 in Supplement 1). The only significantly higher risk estimate in iPSYCH2015 was for 22q11.2-dup in SSD, for which we find no indication of the protective effect reported in the comparison study.6 Our sensitivity analysis found very high correlation across all diagnosis groups between the CPH-derived HRs and corresponding odds ratios (ORs) derived from logistic regression models of iPSYCH2015 data (Pearson correlation test: r = 0.94; P &; 2.2 ×ĉ10−16) (eTable 10 in Supplement 2).
As no individual rCNV was significantly associated with increased risk of MDD, we assessed the omnibus effect of rCNV carrier status on MDD diagnosis and found no effect of rCNV carrier status on MDD prediction (χ2 = 40.32; P = .28); a stark contrast to corresponding results for ADHD, ASD, and SSD (χ2 > 104.68 and P &; 1.27 ×ĉ10−8 for each). Therefore, we excluded MDD from the following comparative analyses of rCNV-associated risk across diagnosis groups, between deletions and duplications, and across loci.
We tested if rCNV-associated risk varied across outcomes through complementary approaches made possible by having information on multiple outcomes in the same population-based sample. First, we plotted CNV-associated HR estimates for each pair of the remaining diagnosis groups (ADHD, ASD, and SSD) and observed high pairwise correlations (rڴ-ٱٱ = 0.68, P = 2 ×ĉ10−5; rڴ-ٱ = 0.71, P = 7 ×ĉ10−6; rڴٱ-ٱ = 0.75, P = 2 ×ĉ10−6) (Figure 3). Then, we fitted nested GEE models with and without an interaction term between the diagnosis group and overall rCNV carrier status. This revealed a significant overall difference in rCNV-associated effects between ASD and SSD (χ2 = 7.47, P = .006), and ADHD and SSD (χ2 = 10.08, P = .002), but not between ASD and ADHD (χ2 = 0.21, P = .64). Finally, we used the same approach but replaced the overall rCNV carrier variable with a categorical rCNV variable. This analysis showed a significant effect of adding the rCNV-diagnosis interaction term in all 3 comparisons, suggesting that for each pair of disorders at least some rCNVs differ in effect size (χ2ڴ-ٱٱ = 94.6, P = 4.2 ×ĉ10−8; χ2ڴ-ٱ = 72.9, P = 3.1 ×ĉ10−5; χ2ڴٱ-ٱ = 50.1, P = .01). Post hoc analyses of the interaction coefficients from the GEE analyses revealed no consistent trend for rCNVs conferring higher or lower risk of any of the 3 diagnoses compared with either of the other 2 (eTable 11 and eFigure 5 in Supplement 1), but among rCNVs with significant interaction coefficients (FDR P < .05) HRs for both PWAS-dup and 16p11.2-del were higher for ASD than for ADHD or SSD, and HRs for 15q11.2-del and 22q11.2-dup were higher for ASD and ADHD than SSD (Figure 3).
The distribution of HR estimates for deletions vs duplications across loci for the 3 main diagnosis groups revealed no indication of a consistent dosage-dependent effect (Figure 4), and we found no evidence of an omnibus effect of dosage type on rCNV-associated risk in a formal test across all 3 diagnosis groups (χ2 = 0.11, P = .73). Also, a comparison of the GEE models with and without an interaction term between rCNV dosage type and locus found no indication of significant locus-specific rCNV dosage type effects (χ2 = 20.66, P = .14).
Hypothesizing that rCNV-associated risk correlated with gene content or locus size, we plotted estimated HRs for the 3 disorders against locus size and locus-wide loss-of-function observed/expected upper-bound fraction (LOEUF) score29 across loci (Figure 5 in Supplement 1 and eTable 12 in Supplement 2). We also hypothesized that rCNV-associated risk correlated with their population-based prevalence or the difference in prevalence between iPSYCH2015 and UKB11 and plotted rCNV prevalence and the iPSYCH2015/UKB prevalence ratio (eFigure 6 in Supplement 1) by fitting log-linear trend lines for deletion- and duplication-associated estimates separately (excluding PWAS because of outlying locus size and HRs). While differing somewhat across diagnoses, the plots indicate that HRs tend to increase with locus size (Figure 5A), LOEUF score (Figure 5B), and iPSYCH/UKB rCNV prevalence ratio (eFigure 6B in Supplement 1), while decreasing with rCNV prevalence itself (eFigure 6A in Supplement 1). A formal test of these trends with nested GEE models confirmed that both locus size and LOEUF score were significantly associated with overall case prediction (size β= = 0.09, P = .01; LOEUF β = 0.10, P = .002) as was the iPSYCH/UKB rCNV prevalence ratio (β = 0.12, P = .009), whereas iPSYCH2015 prevalence was not (β = −0.02, P = .31).
We applied weighted CPH models to estimate the risk of epilepsy and intellectual disability (ID) associated with rCNVs in our dataset (eMethods in Supplement 1). These disorders were not targeted specifically in the iPSYCH2015 case ascertainment, and the risk estimates are therefore not representative of the source population. Nine rCNVs were significantly associated with increased risk of both epilepsy and ID (out of 28 and 33 tested, respectively), and a further 10 rCNVs were associated with increased risk of one or the other disorder (eFigure 7 and eTable 13 in Supplement 1).
Our assessment of rCNV-associated risk across common psychiatric disorders extends beyond current knowledge because of 2 features of the study design. First, the systematic assessment of deletions and duplications across most known rCNV loci enabled us to fill a gap in the existing literature, especially for ADHD and ASD where our samples are larger and involve more rCNVs than earlier studies.1,5 Second, the unbiased population-based study design18 allowed critical revising of estimates of rCNV prevalence and associated risk of psychiatric disorders.
iPSYCH2015 includes 5 times as many ASD cases and 3 times as many ADHD cases as the largest published rCNV case-control studies,1,5 and risk estimates for a majority of the examined rCNV are novel for both disorders. We observed 6 and 7 rCNVs associated with increased risk of ASD and ADHD, respectively, that had either not been included or not found significant in previous studies. For ASD, we confirmed high risk associated with PWAS-dup and 17q12-del, while 16p11.2 (del and dup) were associated with less increased risk than earlier reported.1 For ADHD-associated risk, increases were modest, with HRs ranging from 1.4 at 15q11.2 to 4.2 at 17q12.
While our SSD analysis includes fewer cases than published case-control meta-analyses,6-8 we provide risk estimates for 9 rCNVs not included in those studies. We confirm high risk associated with PWAS-dup and moderate risk (HRs, 2.0-4.8) associated with 7 other rCNVs. For 6 rCNVs, our estimates differed significantly from previous studies,6-8 most drastically so for 22q11.2-del (previously reported with OR >506 but with an associated HR of only 3.3 in iPSYCH2015) and 22q11.2-dup (previously reported with strong protective effect against schizophrenia6 but not significantly affecting SSD risk in iPSYCH2015).
The comparison of rCNV risk profiles across ASD, ADHD, and SSD suggests that PWAS-dup and 16p11.2-del are associated with a greater increase in risk of ASD than of ADHD or SSD and that 15q11.2-del and 22q11.2-dup are associated with less (if at all any) increase in risk of SSD than of ADHD or ASD. Notwithstanding these interesting signs of variability in risk across diagnosis groups, we observed high correlation of risk estimates for rCNVs across these disorders.
The main exception from this pleiotropic trend was MDD, for which we found no evidence of an overall rCNV-associated increase in risk. This contrasts with the results of 2 US health care system population studies, where overall rCNV carriage associated with a significant increase (OR = 1.5) of depression and/or anxiety risk13,14 A limitation to such comparison of overall rCNV effect on risk is that we were unable to include some of the rarest (and potentially most pathogenic) rCNVs that we originally set out to test. A study of the UKB also reported modest increases in MDD risk for a few rCNVs.10 The BPD case sample in our study is small and does not add meaningfully to the limited evidence for rCNV involvement in BPD.9
Although meiotic nonallelic homologous recombination produces de novo rCNV deletions more often than duplications,30 rCNV studies consistently report higher population prevalence of duplications.4,11 The general view is that this prevalence disparity reflects a greater effect on viability and reproductive fitness associated with deletions than duplications, a part of which has been attributed to psychiatric disorders. However, we do not observe greater risk associated with deletions than duplications for ASD, ADHD, or SSD. Across the 17 rCNV loci included in our analysis, both locus size and locus-wide LOEUF scores were independently associated with increased disease risk while rCNV dosage type was not. This is consistent with a recent study that found the summed constraint of overlapped genes to best explain CNV-associated risk with ASD.31
While the case-cohort design of iPSYCH2015 is optimal to assess rCNV-associated risk of its targeted psychiatric disorders, the relatively young age of participants limited the study power for late-onset psychiatric illness, such as BPD, and the quality of the probe intensity data used to infer CNV carrier status was not as high in the dried blood spot samples used in our study as when using whole blood samples. Also, because of low prevalence, we could not derive estimates for rCNVs at 9 of 27 loci, including some with high reported risk of psychiatric disorders in previous studies (such as 3q29-del and WBS-dup). Third, while having nationwide coverage of hospital-based diagnoses, iPSYCH2015 does not include individuals diagnosed and treated solely outside the hospital regimen, and we additionally note that the diagnostic criteria for specific disorders in some instances differ between iPSYCH2015 and the respective previously published case-control study that we draw comparison to. Finally, iPSYCH2015 was designed to address population-valid risk associated with its targeted common psychiatric disorders, and our observations of less increased risk of these disorders associated with, eg, 22q11.2-del, do not mean that the same is necessarily true of other disorders, such as ID. In fact, from hospital records on ID and epilepsy in our study, we find evidence of high risk of both disorders associated with 22q11.2-del.
This study found that several rCNVs were more prevalent and conferred less risk of psychiatric disorders than estimated previously. To summarize, we propose that (1) overall, rCNVs are associated with increased risk of psychiatric disorders to a degree that varies according to locus size and gene content but not according to gene dosage; (2) the risk increase is overall similar across ADHD, ASD, and SSD but much less (if at all) pronounced for MDD; (3) some rCNVs, such as PWAS-dup, may not adhere to these general statements; and (4) case-control studies have in some instances overestimated rCNV-associated risk of psychiatric disorders due to selection bias and screening of controls. In an era where genetics is increasingly being used in clinical practice, our results highlight the importance of accurate population-based risk estimates for genetic exposures such as rCNVs.
Accepted for Publication: April 17, 2024.
Published Online: June 26, 2024. doi:10.1001/jamapsychiatry.2024.1453
Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2024 Vaez M et al. JAMA Psychiatry.
Corresponding Author: Andrés Ingason, PhD, Forskningsinstitut bygn 134, Boserupvej 2, DK-4000 Roskilde, Denmark (andres.ingason@regionh.dk).
Author Contributions: Dr Ingason had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: Vaez, Montalbano, Mortensen, Geschwind, Schork, Helenius, Raznahan, Thompson, Werge, Ingason.
Acquisition, analysis, or interpretation of data: Vaez, Montalbano, Calle Sánchez, Georgii Hellberg, Dehkordi, Krebs, Meijsen, Shorter, Bybjerg-Grauholm, Børglum, Hougaard, Nordentoft, Buil, Helenius, Raznahan, Werge, Ingason.
Drafting of the manuscript: Vaez, Montalbano, Georgii Hellberg, Dehkordi, Raznahan, Thompson, Ingason.
Critical review of the manuscript for important intellectual content: Vaez, Montalbano, Calle Sánchez, Georgii Hellberg, Krebs, Meijsen, Shorter, Bybjerg-Grauholm, Mortensen, Børglum, Hougaard, Nordentoft, Geschwind, Buil, Schork, Helenius, Raznahan, Werge, Ingason.
Statistical analysis: Vaez, Montalbano, Calle Sánchez, Dehkordi, Krebs, Meijsen, Shorter, Helenius, Raznahan, Thompson, Werge, Ingason.
Obtained funding: Vaez, Calle Sánchez, Mortensen, Børglum, Nordentoft, Werge, Ingason.
Administrative, technical, or material support: Vaez, Meijsen, Bybjerg-Grauholm, Mortensen, Børglum, Hougaard, Werge, Ingason.
Supervision: Geschwind, Buil, Schork, Helenius, Werge, Ingason.
Conflict of Interest Disclosures: Dr Børglum reported speaking fees from Lundbeck Foundation outside the submitted work. Dr Werge reported being a member of the board of directors of the Lundbeck Foundation (considered independent). No other disclosures were reported.
Funding/Support: This study was supported by grants from the Lundbeck Foundation (R165-2013-15320, R155-2014-1724, and R248-2017-2003), National Institutes of Health (R01-MH124789-01 and 1R01-MH124851-01), the UK Research and Innovation Medical Research Council (MR/W014416/1), and EU Horizon Europe (grant agreement 101057385 [R2D2-MH]).
Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Group Information: The iPSYCH investigators are listed in Supplement 3.
Data Sharing Statement: See Supplement 4.
1.Malhotra
D, Sebat
J. CNVs: harbingers of a rare variant revolution in psychiatric genetics. . 2012;148(6):1223-1241. doi:
2.Rees
E, Kirov
G. Copy number variation and neuropsychiatric illness. Curr Opin Genet Dev. 2021;68:57-63. doi:
3.Stefansson
H, Meyer-Lindenberg
A, Steinberg
S,
et al. CNVs conferring risk of autism or schizophrenia affect cognition in controls. ٳܰ. 2014;505(7483):361-366. doi:
4.Coe
BP, Witherspoon
K, Rosenfeld
JA,
et al. Refining analyses of copy number variation identifies specific genes associated with developmental delay. Nat Genet. 2014;46(10):1063-1071. doi:
5.Gudmundsson
OO, Walters
GB, Ingason
A,
et al. Attention-deficit hyperactivity disorder shares copy number variant risk with schizophrenia and autism spectrum disorder. Transl Psychiatry. 2019;9(1):258. doi:
6.Marshall
CR, Howrigan
DP, Merico
D,
et al; Psychosis Endophenotypes International Consortium; CNV and Schizophrenia Working Groups of the Psychiatric Genomics Consortium. Contribution of copy number variants to schizophrenia from a genome-wide study of 41,321 subjects. Nat Genet. 2017;49(1):27-35. doi:
7.Rees
E, Kendall
K, Pardiñas
AF,
et al. Analysis of intellectual disability copy number variants for association with schizophrenia. JAMA Psychiatry. 2016;73(9):963-969. doi:
8.Rees
E, Walters
JT, Georgieva
L,
et al. Analysis of copy number variations at 15 schizophrenia-associated loci. Br J Psychiatry. 2014;204(2):108-114. doi:
9.Green
EK, Rees
E, Walters
JTR,
et al. Copy number variation in bipolar disorder. Mol Psychiatry. 2016;21(1):89-93. doi:
10.Kendall
KM, Rees
E, Bracher-Smith
M,
et al. Association of rare copy number variants with risk of depression. JAMA Psychiatry. 2019;76(8):818-825. doi:
11.Crawford
K, Bracher-Smith
M, Owen
D,
et al. Medical consequences of pathogenic CNVs in adults: analysis of the UK Biobank. J Med Genet. 2019;56(3):131-138. doi:
12.Schoeler
T, Speed
D, Porcu
E, Pirastu
N, Pingault
JB, Kutalik
Z. Participation bias in the UK Biobank distorts genetic associations and downstream analyses. Nat Hum Behav. 2023;7(7):1216-1227. doi:
13.Martin
CL, Wain
KE, Oetjens
MT,
et al. Identification of neuropsychiatric copy number variants in a health care system population. JAMA Psychiatry. 2020;77(12):1276-1285. doi:
14.Birnbaum
R, Mahjani
B, Loos
RJF, Sharp
AJ. Clinical characterization of copy number variants associated with neurodevelopmental disorders in a large-scale multiancestry biobank. JAMA Psychiatry. 2022;79(3):250-259. doi:
15.Pedersen
CB, Bybjerg-Grauholm
J, Pedersen
MG,
et al. The iPSYCH2012 case-cohort sample: new directions for unravelling genetic and environmental architectures of severe mental disorders. Mol Psychiatry. 2018;23(1):6-14. doi:
16.Olsen
L, Sparsø
T, Weinsheimer
SM,
et al. Prevalence of rearrangements in the 22q11.2 region and population-based risk of neuropsychiatric and developmental disorders in a Danish population: a case-cohort study. Lancet Psychiatry. 2018;5(7):573-580. doi:
17.Calle Sánchez
X, Helenius
D, Bybjerg-Grauholm
J,
et al. Comparing copy number variations in a Danish case cohort of individuals with psychiatric disorders. JAMA Psychiatry. 2022;79(1):59-69. doi:
18.Bybjerg-Grauholm
J, Pedersen
CB, Bækvad-Hansen
M,
et al. The iPSYCH2015 case-cohort sample: updated directions for unravelling genetic and environmental architectures of severe mental disorders. 澱. 2020;11(30):20237768. doi:
19.Vandenbroucke
JP, von Elm
E, Altman
DG,
et al; STROBE Initiative. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE): explanation and elaboration. PLoS Med. 2007;4(10):e297. doi:
20.Sharp
SJ, Poulaliou
M, Thompson
SG, White
IR, Wood
AM. A review of published analyses of case-cohort studies and recommendations for future reporting. PLoS One. 2014;9(6):e101176. doi:
21.Pedersen
CB. The Danish Civil Registration System. Scand J Public Health. 2011;39(7)(suppl):22-25. doi:
22.Mors
O, Perto
GP, Mortensen
PB. The Danish Psychiatric Central Research Register. Scand J Public Health. 2011;39(7)(suppl):54-57. doi:
23.Schmidt
M, Schmidt
SAJ, Sandegaard
JL, Ehrenstein
V, Pedersen
L, Sørensen
HT. The Danish National Patient Registry: a review of content, data quality, and research potential. Clin Epidemiol. 2015;7:449-490. doi:
24.Hollegaard
MV, Grove
J, Grauholm
J,
et al. Robustness of genome-wide scanning using archived dried blood spot samples as a DNA source. BMC Genet. 2011;12:58. doi:
25.Rehm
HL, Berg
JS, Brooks
LD,
et al; ClinGen. ClinGen–the Clinical Genome Resource. N Engl J Med. 2015;372(23):2235-2242. doi:
26.Riggs
ER, Nelson
T, Merz
A,
et al. Copy number variant discrepancy resolution using the ClinGen dosage sensitivity map results in updated clinical interpretations in ClinVar. Hum Mutat. 2018;39(11):1650-1659. doi:
27.Wang
K, Li
M, Hadley
D,
et al. PennCNV: an integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Res. 2007;17(11):1665-1674. doi:
28.Montalbano
S, Sánchez
XC, Vaez
M, Helenius
D, Werge
T, Ingason
A. Accurate and effective detection of recurrent copy number variants in large SNP genotype datasets. Curr Protoc. 2022;2(12):e621. doi:
29.Karczewski
KJ, Francioli
LC, Tiao
G,
et al; Genome Aggregation Database Consortium. The mutational constraint spectrum quantified from variation in 141,456 humans. ٳܰ. 2020;581(7809):434-443. doi:
30.Turner
DJ, Miretti
M, Rajan
D,
et al. Germline rates of de novo meiotic deletions and duplications causing several genomic disorders. Nat Genet. 2008;40(1):90-95. doi:
31.Douard
E, Zeribi
A, Schramm
C,
et al. Effect sizes of deletions and duplications on autism risk across the genome. Am J Psychiatry. 2021;178(1):87-98. doi: