Individual Copy Number Variants Research Articles

Increased copy-number variant load of associated risk genes in sporadic cases of amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is an age-related neurodegenerative disease characterized by selective loss of motor neurons in the brainstem and spinal cord. Several genetic factors have been associated to ALS, ranging from causal genes and potential risk factors to disease modifiers. The search for pathogenic variants in these genes has mostly focused on single nucleotide variants (SNVs) while relatively understudied and not fully elucidated is the contribution of structural variants, such as copy number variations (CNVs). Here, we applied an exon-centric aCGH method to investigate, in sporadic ALS patients, the load of CNVs in 131 genes previously associated to ALS. Our approach revealed that CNV load, defined as the total number of CNVs or their size, was significantly higher in ALS cases than controls. About 87% of patients harbored multiple CNVs in ALS-related genes, and 75% structural variants compromised genes directly implicated in ALS pathogenesis (C9orf72, CHCHD10, EPHA4, FUS, HNRNPA1, KIF5A, NEK1, OPTN, PFN1, SOD1, TARDBP, TBK1, UBQLN2, UNC13A, VAPB, VCP). CNV load was also associated to higher onset age and disease progression rate. Although the contribution of individual CNVs in ALS is still unknown, their extensive load in disease-related genes may have relevant implications for the diagnostic, prognostic and therapeutical management of this devastating disorder.

Genome-wide identification and phenotypic characterization of seizure-associated copy number variations in 741,075 individuals

Copy number variants (CNV) are established risk factors for neurodevelopmental disorders with seizures or epilepsy. With the hypothesis that seizure disorders share genetic risk factors, we pooled CNV data from 10,590 individuals with seizure disorders, 16,109 individuals with clinically validated epilepsy, and 492,324 population controls and identified 25 genome-wide significant loci, 22 of which are novel for seizure disorders, such as deletions at 1p36.33, 1q44, 2p21-p16.3, 3q29, 8p23.3-p23.2, 9p24.3, 10q26.3, 15q11.2, 15q12-q13.1, 16p12.2, 17q21.31, duplications at 2q13, 9q34.3, 16p13.3, 17q12, 19p13.3, 20q13.33, and reciprocal CNVs at 16p11.2, and 22q11.21. Using genetic data from additional 248,751 individuals with 23 neuropsychiatric phenotypes, we explored the pleiotropy of these 25 loci. Finally, in a subset of individuals with epilepsy and detailed clinical data available, we performed phenome-wide association analyses between individual CNVs and clinical annotations categorized through the Human Phenotype Ontology (HPO). For six CNVs, we identified 19 significant associations with specific HPO terms and generated, for all CNVs, phenotype signatures across 17 clinical categories relevant for epileptologists. This is the most comprehensive investigation of CNVs in epilepsy and related seizure disorders, with potential implications for clinical practice.

Impact of Copy Number Variants and Polygenic Risk Scores on Psychopathology in the UK Biobank

BackgroundOur understanding of the impact of copy number variants (CNVs) on psychopathology and their joint influence with polygenic risk scores (PRSs) remains limited. MethodsThe UK Biobank recruited 502,534 individuals ages 37 to 73 years living in the United Kingdom between 2006 and 2010. After quality control, genotype data from 459,855 individuals were available for CNV calling. A total of 61 commonly studied recurrent neuropsychiatric CNVs were selected for analyses and examined individually and in aggregate (any CNV, deletion, or duplication). CNV risk scores were used to quantify intolerance of CNVs to haploinsufficiency. Major depressive disorder and generalized anxiety disorder PRSs were generated for White British individuals (N = 408,870). Mood/anxiety factor scores were generated using item-level questionnaire data (N = 501,289). ResultsCNV carriers showed higher mood/anxiety scores than noncarriers, with the largest effects seen for intolerant deletions. A total of 11 individual deletions and 8 duplications were associated with higher mood/anxiety. Carriers of the 9p24.3 (DMRT1) duplication showed lower mood/anxiety. Associations remained significant for most CNVs when excluding individuals with psychiatric diagnoses. Nominally significant CNV × PRS interactions provided preliminary evidence that associations between select individual CNVs, but not CNVs in aggregate, and mood/anxiety may be modulated by PRSs. ConclusionsCNVs associated with risk for psychiatric disorders showed small to large effects on dimensional mood/anxiety scores in a general population cohort, even when excluding individuals with psychiatric diagnoses. CNV × PRS interactions showed that associations between select CNVs and mood/anxiety may be modulated by PRSs.

The contribution of copy number variants to psychiatric symptoms and cognitive ability.

Copy number variants (CNVs) are deletions and duplications of DNA sequence. The most frequently studied CNVs, which are described in this review, are recurrent CNVs that occur in the same locations on the genome. These CNVs have been strongly implicated in neurodevelopmental disorders, namely autism spectrum disorder (ASD), intellectual disability (ID), and developmental delay (DD), but also in schizophrenia. More recent work has also shown that CNVs increase risk for other psychiatric disorders, namely, depression, bipolar disorder, and post-traumatic stress disorder. Many of the same CNVs are implicated across all of these disorders, and these neuropsychiatric CNVs are also associated with cognitive ability in the general population, as well as with structural and functional brain alterations. Neuropsychiatric CNVs also show incomplete penetrance, such that carriers do not always develop any psychiatric disorder, and may show only mild symptoms, if any. Variable expressivity, whereby the same CNVs are associated with many different phenotypes of varied severity, also points to highly complex mechanisms underlying disease risk in CNV carriers. Comprehensive and longitudinal phenotyping studies of individual CNVs have provided initial insights into these mechanisms. However, more work is needed to estimate and predict the effect of non-recurrent, ultra-rare CNVs, which also contribute to psychiatric and cognitive outcomes. Moreover, delineating the broader phenotypic landscape of neuropsychiatric CNVs in both clinical and general population cohorts may also offer important mechanistic insights.

Widespread occurrence of copy number variants and fixation of pyrethroid target site resistance in Anopheles gambiae (s.l.) from southern Côte d’Ivoire

Resistance to pyrethroid and organophosphate insecticides in the malaria vector Anopheles gambiae (s.l.) is conferred by a variety of genetic mutations, including single nucleotide polymorphisms (SNPs) and copy number variants (CNVs). Knowledge of the distribution of these mutations in mosquito populations is a prerequisite for establishing better strategies for their management. In this study, a total of 755 Anopheles gambiae (s.l.) from southern Côte d’Ivoire were exposed to deltamethrin or pirimiphos-methyl insecticides and were screened to assess the distribution of SNPs and CNVs known or believed to confer resistance to one or other of the insecticide classes. Most individuals from the An. gambiae (s.l.) complex were identified by molecular tests as Anopheles coluzzii. Survival to deltamethrin (from 94% to 97%) was higher than to pirimiphos-methyl (from 10% to 49%). In An. gambiae (s.s.), the SNP in the Voltage Gated Sodium Channel (Vgsc) at the 995F locus (Vgsc-995F) was fixed, while other target site mutations were rare or absent (Vgsc-402L: 0%; Vgsc-1570Y: 0%, Acetylcholinesterase Acel-280S: 14%). In An. coluzzii, Vgsc-995F was the target site SNP found at highest frequency (65%) followed by other target site mutations (Vgsc-402L: 36%; Vgsc-1570Y: 0.33%; Acel-280S: 45%). The Vgsc-995S SNP was not present. The presence of the Ace1-280S SNP was found to be significantly linked to the presence of the Ace1-CNV, Ace1_AgDup. Significant association was found between the presence of the Ace1_AgDup and pirimiphos-methyl resistance in An. gambiae (s.s.) but not in An. coluzzii. The deletion Ace1_Del97 was found in one specimen of An. gambiae (s.s.). Four CNVs in the Cyp6aa/Cyp6p gene cluster, which contains genes of known importance for resistance, were detected in An. coluzzii, the most frequent being Dup 7 (42%) and Dup 14 (26%). While none of these individual CNV alleles were significantly associated with resistance, copy number in the Cyp6aa gene region in general was associated with increased resistance to deltamethrin. Elevated expression of Cyp6p3 was nearly associated with deltamethrin resistance, although there was no association of resistance with copy number. Use of alternative insecticides and control methods to arrest resistance spread in An. coluzzii populations is merited.

Whole exome sequencing reveals copy number variants in individuals with disorders of sex development

Whole exome sequencing reveals copy number variants in individuals with disorders of sex development

Whole exome sequencing reveals copy number variants in individuals with disorders of sex development

Complete androgen insensitivity syndrome (CAIS), where 46,XY individuals present as female, is caused by variants in the androgen receptor gene (AR). We analyzed the DNA of a patient with suspected CAIS using a targeted gene sequencing panel and whole exome sequencing (WES) but did not detect any small nucleotide variants in AR. Analysis of WES data using our bioinformatics pipeline designed to detect copy number variations (CNV) uncovered a rare duplication of exon 2 of AR. Using array comparative genomic hybridization, the duplication was found to span 43.6 kb and is predicted to cause a frameshift and loss of AR protein. We confirmed the power of our WES-CNV detection protocol by identifying pathogenic CNVs in FSHR and NR5A1 in previously undiagnosed patients with disorders of sex development. Our findings illustrate the usefulness of CNV analysis in WES data to detect pathogenic genomic changes that may go undetected using only standard analysis protocols.

16. Copy number variants at Chr16p11.2 and relevant clinical presentations

Recurrent microdeletions and microduplications on chromosome 16p11.2 (MIM: 611913 and 614671) confer susceptibility to autism spectrum disorder (ASD) in ∼1% of individuals with ASD, in addition to other clinical features. To further analyze copy number variants (CNVs) in this region, we studied 41 individuals with 16p11.2 duplications and 53 individuals with 16p11.2 deletions at Greenwood Genetic Center. They share common clinical features including: developmental delay, intellectual disability, autism spectrum disorder, seizures, and some features that are unique to either deletions or duplications. We summarize unique CNVs in individuals with either deletion (n=41) or duplication (n=35) of these regions. By overlapping the del/dup regions, we identify six subtypes of CNVs and classify them by potential relevant breakpoints based on size and location. Types 1, 2 and 3 are more common with both deletions and duplications found in these regions, and they overlap with BP1-BP3, BP1-BP2 and BP4-BP5 regions, respectively, in the literature. Types 4, 5 and 6 are less frequent, locating at 31.5-32.5 Mb, 29.6-47.5 Mb and 28.48-28.49 Mb regions on chromosome 16. Type 5 is found only in individuals with duplications and type 6 only in individuals with deletions. Type 4 variants are more frequently observed in individuals of African American decent, yet no clear association could be determined due to small sample size in this study. Genotype-phenotype correlation analysis was also performed to narrow down the potentially causative genes for specific clinical features. It is also worth noting that de novo CNVs are observed with much higher frequency in individuals with 16p11.2 deletions (68.2%, n=22) compared to those with duplications (28.6%, n=14). In conclusion, our analysis explores 16p11.2 CNVs in more detail for rarely observed regions, and could broaden the understanding of the clinical significance of genes and non-coding regions in this chromosomal region.

ClassifyCNV: a tool for clinical annotation of copy-number variants

Copy-number variants (CNVs) are an important part of human genetic variation. They can be benign or can play a role in human disease by creating dosage imbalances and disrupting genes and regulatory elements. Accurate identification and clinical annotation of CNVs is essential, however, manual evaluation of individual CNVs by clinicians is challenging on a large scale. Here, we present ClassifyCNV, an easy-to-use tool that implements the 2019 ACMG classification guidelines to assess CNV pathogenicity. ClassifyCNV uses genomic coordinates and CNV type as input and reports a clinical classification for each variant, a classification score breakdown, and a list of genes of potential importance for variant interpretation. We validate ClassifyCNV’s performance using a set of known clinical CNVs and a set of manually evaluated variants. ClassifyCNV matches the pathogenicity category for 81% of manually evaluated variants with the significance of the remaining pathogenic and benign variants automatically determined as uncertain, requiring a further evaluation by a clinician. ClassifyCNV facilitates the implementation of the latest ACMG guidelines in high-throughput CNV analysis, is suitable for integration into NGS analysis pipelines, and can decrease time to diagnosis. The tool is available at https://github.com/Genotek/ClassifyCNV.

Association of CNVs with methylation variation

Germline copy number variants (CNVs) and single-nucleotide polymorphisms (SNPs) form the basis of inter-individual genetic variation. Although the phenotypic effects of SNPs have been extensively investigated, the effects of CNVs is relatively less understood. To better characterize mechanisms by which CNVs affect cellular phenotype, we tested their association with variable CpG methylation in a genome-wide manner. Using paired CNV and methylation data from the 1000 genomes and HapMap projects, we identified genome-wide associations by methylation quantitative trait locus (mQTL) analysis. We found individual CNVs being associated with methylation of multiple CpGs and vice versa. CNV-associated methylation changes were correlated with gene expression. CNV-mQTLs were enriched for regulatory regions, transcription factor-binding sites (TFBSs), and were involved in long-range physical interactions with associated CpGs. Some CNV-mQTLs were associated with methylation of imprinted genes. Several CNV-mQTLs and/or associated genes were among those previously reported by genome-wide association studies (GWASs). We demonstrate that germline CNVs in the genome are associated with CpG methylation. Our findings suggest that structural variation together with methylation may affect cellular phenotype.

Estimating the effects of copy-number variants on intelligence using hierarchical Bayesian models.

It is challenging to estimate the phenotypic impact of the structural genome changes known as copy-number variations (CNVs), since there are many unique CNVs which are nonrecurrent, and most are too rare to be studied individually. In recent work, we found that CNV-aggregated genomic annotations, that is, specifically the intolerance to mutation as measured by the pLI score (probability of being loss-of-function intolerant), can be strong predictors of intellectual quotient (IQ) loss. However, this aggregation method only estimates the individual CNV effects indirectly. Here, we propose the use of hierarchical Bayesian models to directly estimate individual effects of rare CNVs on measures of intelligence. Annotation information on the impact of major mutations in genomic regions is extracted from genomic databases and used to define prior information for the approach we call HBIQ. We applied HBIQ to the analysis of CNV deletions and duplications from three datasets and identified several genomic regions containing CNVs demonstrating significant deleterious effects on IQ, some of which validate previously known associations. We also show that several CNVs were identified as deleterious by HBIQ even if they have a zero pLI score, and the converse is also true. Furthermore, we show that our new model yields higher out-of-sample concordance (78%) for predicting the consequences of carrying known recurrent CNVs compared with our previous approach.

Increasing the Clinical Psychiatric Knowledge Base About Pathogenic Copy Number Variation.

Specific copy number variants (CNVs) have been robustly associated with intellectual disability, autism, and schizophrenia. Most of the literature focus has been on documenting the existence of these phenomena. There are few data to guide therapeutic choices for these “orphan” diseases. We call for systematic and longitudinal case reports which, if carefully conducted, may provide crucial initial knowledge to guide therapeutics. We provide a step-by-step overview, a tailored set of consensus criteria for high-quality case reports, and a specific set of learning resources.

Critical exon indexing improves clinical interpretation of copy number variants in neurodevelopmental disorders.

ObjectiveTo evaluate a new tool to aid interpretation of copy number variants (CNVs) in individuals with neurodevelopmental disabilities.MethodsCritical exon indexing (CEI) was used to identify genes with critical exons (CEGs) from clinically reported CNVs, which may contribute to neurodevelopmental disorders (NDDs). The 742 pathogenic CNVs and 1,363 variants of unknown significance (VUS) identified by chromosomal microarray analysis in 5,487 individuals with NDDs were subjected to CEI to identify CEGs. CEGs identified in a subsequent random series of VUS were evaluated for relevance to CNV interpretation.ResultsCEI identified a total of 2,492 unique CEGs in pathogenic CNVs and 953 in VUS compared with 259 CEGs in 6,965 CNVs from 873 controls. These differences are highly significant (p < 0.00001) whether compared as frequency, average, or normalized by CNV size. Twenty-one percent of VUS CEGs were not represented in Online Mendelian Inheritance in Man, highlighting limitations of existing resources for identifying potentially impactful genes within CNVs. CEGs were highly correlated with other indices and known pathways of relevance. Separately, 136 random VUS reports were reevaluated, and 76% of CEGs had not been commented on. In multiple cases, further investigation yielded additional relevant literature aiding interpretation. As one specific example, we discuss GTF2I as a CEG, which likely alters interpretation of several reported duplication VUS in the Williams-Beuren region.ConclusionsApplication of CEI to CNVs in individuals with NDDs can identify genes of potential clinical relevance, aid laboratories in effectively searching the clinical literature, and support the clinical reporting of poorly annotated VUS.

Rare copy number variants in individuals at clinical high risk for psychosis: Enrichment of synaptic/brain-related functional pathways.

Schizophrenia is a complex and chronic neuropsychiatric disorder, with a heritability of around 60-80%. Large (>100 kb) rare (<1%) copy number variants (CNVs) occur more frequently in schizophrenia patients compared to controls. Currently, there are no studies reporting genome-wide CNVs in clinical high risk for psychosis (CHR-P) individuals. The aim of this study was to investigate the role of rare genome-wide CNVs in 84 CHR-P individuals and 124 presumably healthy controls. There were no significant differences in all rare CNV frequencies and sizes between CHR-P individuals and controls. However, brain-related CNVs and brain-related deletions were significantly more frequent in CHR-P individuals than controls. In CHR-P individuals, significant associations were found between brain-related CNV carriers and attenuated positive symptoms syndrome or cognitive disturbances (OR = 3.07, p = .0286). Brain-related CNV carriers experienced significantly higher negative symptoms (p = .0047), higher depressive symptoms (p = .0175), and higher disturbances of self and surroundings (p = .0029) than noncarriers. Furthermore, enrichment analysis of genes was performed in the regions of rare CNVs using three independent methods, which confirmed significant clustering of predefined genes involved in synaptic/brain-related functional pathways in CHR-P individuals. These results suggest that rare CNVs might affect synaptic/brain-related functional pathways in CHR-P individuals.

CNVRanger: association analysis of CNVs with gene expression and quantitative phenotypes.

Copy number variation (CNV) is a major type of structural genomic variation that is increasingly studied across different species for association with diseases and production traits. Established protocols for experimental detection and computational inference of CNVs from SNP array and next-generation sequencing data are available. We present the CNVRanger R/Bioconductor package which implements a comprehensive toolbox for structured downstream analysis of CNVs. This includes functionality for summarizing individual CNV calls across a population, assessing overlap with functional genomic regions, and genome-wide association analysis with gene expression and quantitative phenotypes. http://bioconductor.org/packages/CNVRanger.

Copy number variations in individuals with conotruncal heart defects reveal some shared developmental pathways irrespective of 22q11.2 deletion status.

Over 50% of patients with 22q11.2 deletion syndrome (DS) have a conotruncal or related cardiac defect (CTRD). We hypothesized that similar genetic variants, developmental pathways and biological functions, contribute to disease risk for CTRD in patients without a 22q11.2 deletion (ND-CTRD) and with a 22q11.2 deletion (DS-CTRD). To test this hypothesis, we performed rare CNV (rCNV)-based analyses on 630 ND-CTRD cases and 602 DS-CTRD cases with comparable cardiac lesions separately and jointly. First, we detected a collection of heart development related pathways from Gene Ontology and Mammalian Phenotype Ontology analysis. We then constructed gene regulation networks using unique genes collected from the rCNVs found in the ND-CTRD and DS-CTRD cohorts. These gene networks were clustered and their predicted functions were examined. We further investigated expression patterns of those unique genes using publicly available mouse embryo microarray expression data from single-cell embryos to fully developed hearts. By these bioinformatics approaches, we identified a commonly shared gene expression pattern in both the ND-CTRD and DS-CTRD cohorts. Computational analysis of gene functions characterized with this expression pattern revealed a collection of significantly enriched terms related to cardiovascular development. By our combined analysis of rCNVs in the ND-CTRD and DS-CTRD cohorts, a group of statistically significant shared pathways, biological functions, and gene expression patterns were identified that can be tested in future studies for their biological relevance.

Possibility of early diagnosis in a fetus affected by Prader‑Willi syndrome with maternal hetero‑UPD15: A lesson to be learned.

Prader-Willi syndrome (PWS), a complicated neurodevelopmental disorder arising from errors in genomic imprinting, is characterized by evident hypotonia along with feeding difficulties and the absence of crying in early infancy. Hyperphagia and obesity are not uncommon in patients with PWS, usually accompanied by intellectual disability, cognitive impairment, short stature, small hands and feet, as well as hypogonadism and typical facial features. Due to the severe complications associated with PWS, a thorough understanding of its features and an early diagnosis, preferably in the fetal period, are important for clinical management. According to previous studies, prenatal diagnosis has been confirmed in only a few cases of PWS, using ultrasound, or as an accidental finding by cytogenetic molecular techniques, as no precise fetal phenotype has been defined. In this present study, an infant with PWS arising from maternal heterodisomy of chromosome 15 is described. This is a typical case of missed diagnosis by fetal ultrasound examination, chromosome karyotype analysis and chromosome microarray (CMA) conducted during the pregnancy. To delineate the complex prenatal characteristics of a fetus with PWS, prenatally-diagnosed cases of PWS described in the literature were reviewed. This present study indicated that although prenatal signs are not sufficient for a diagnosis to be confirmed, a comprehensive consideration of these signs is important in leading to a diagnosis of suspected PWS, and thus prompts further prenatal investigations using molecular genetic tools. Furthermore, this present study also suggested that CMA can lead to a missed diagnosis of PWS/Angelman syndrome and other imprinting disorders despite its high value in the detection of copy-number variants in individuals with developmental delay. If clinical signs strongly suggest PWS, other prenatal molecular genetic investigations, including methylation tests and short tandem repeat-based linkage analysis for uniparental disomy, are recommended as an additional tool to aid diagnosis.

Copy number variants detection by microarray and multiplex ligation-dependent probe amplification in congenital heart diseases

Congenital heart diseases (CHDs) are the most common birth defects among life births, which could be presented as isolated or syndromic with other congenital malformations. The etiology of CHD largely unknown, genetic and environmental factors contribute to the disease. Recurrent copy number variants (CNVs) have been reported in the pathogenesis of CHD. The aim of this study was to evaluate the clinical utility of multiplex ligation-dependent probe amplification (MLPA) and microarray analyses on isolated and syndromic CHD cases and to explore the relationship between identified CNVs and CHD. Eighteen prenatal samples, 16 isolated and 33 syndromic patients with mild to severe CHD phenotype were tested. Prenatal and isolated CHD cases did not show pathogenic CNVs. Clinically significant CNVs were detected in 7/33 (21%) syndromic CHD patients: del 22q11.2 (n = 2), 8p23.1 duplication (n = 2), deletion 5p (n = 1), deletion 6q21q22 (n = 1), unbalanced translocation causing partial deletion of 4q34.3 and duplication of 6q25.1 (n = 1). These genomic imbalances contain genes that has been associated with human CHD before. The present study demonstrates that using microarray and MLPA analysis increase the detection rate of causal CNVs in individuals with syndromic CHD.

Association of Rare Copy Number Variants With Risk of Depression

The role of large, rare copy number variants (CNVs) in neuropsychiatric disorders is well established, but their association with common psychiatric disorders, such as depression, remains unclear. To examine the association of a group of 53 CNVs associated with neurodevelopmental disorders and burden of rare CNVs with risk of depression. This case-control study used data from the UK Biobank study sample, which comprised 502 534 individuals living in the United Kingdom. Individuals with autism spectrum disorder, intellectual disability, attention-deficit/hyperactivity disorder, schizophrenia, or bipolar affective disorder diagnoses were excluded. Analyses were further restricted to individuals of European genetic ancestry (n = 407 074). The study was conducted from January 2017 to September 2018. CNV carrier status. For the primary outcome, individuals who reported that a physician had told them they had a depression diagnosis were defined as cases. Analyses were repeated using 2 alternative depression definitions: self-reported lifetime depression with current antidepressant prescription at the time of visit 1, and hospital discharge diagnosis of depression. Copy number variants were identified in 488 366 individuals aged 37 to 73 years. In total, 407 074 individuals with European genetic ancestry (220 201 female [54.1%]; mean [SD] age of 56.9 [8.0] years) were included in the study. Of these individuals, 23 979 (5.9%) had self-reported lifetime depression and 383 095 (94.1%) reported no lifetime depression. The group of 53 neurodevelopmental CNVs was associated with self-reported depression (odds ratio [OR], 1.34; 95% CI, 1.19-1.49, uncorrected P = 1.38 × 10-7), and these results were consistent when using 2 alternative definitions of depression. This association was partially explained by physical health, educational attainment, social deprivation, smoking status, and alcohol consumption. A strong independent association remained between the neurodevelopmental CNVs and depression in analyses that incorporated these other measures (OR, 1.26; 95% CI, 1.11-1.43; P = 2.87 × 10-4). Eight individual CNVs were nominally associated with risk of depression, and 3 of these 8 CNVs (1q21.1 duplication, Prader-Willi syndrome duplication, and 16p11.2 duplication) survived Bonferroni correction for the 53 CNVs tested. After the exclusion of carriers of neurodevelopmental CNVs, no association was found between measures of CNV burden and depression. Neurodevelopmental CNVs appear to be associated with depression, extending the spectrum of clinical phenotypes that are associated with CNV carrier status.

Genome-wide Burden of Rare Short Deletions Is Enriched in Major Depressive Disorder in Four Cohorts

BackgroundMajor depressive disorder (MDD) is moderately heritable, with a high prevalence and a presumed high heterogeneity. Copy number variants (CNVs) could contribute to the heritable component of risk, but the two previous genome-wide association studies of rare CNVs did not report significant findings. MethodsIn this meta-analysis of four cohorts (5780 patients and 6626 control subjects), we analyzed the association of MDD to 1) genome-wide burden of rare deletions and duplications, partitioned by length (<100 kb or >100 kb) and other characteristics, and 2) individual rare exonic CNVs and CNV regions. ResultsPatients with MDD carried significantly more short deletions than control subjects (p = .0059) but not long deletions or short or long duplications. The confidence interval for long deletions overlapped with that for short deletions, but long deletions were 70% less frequent genome-wide, reducing the power to detect increased burden. The increased burden of short deletions was primarily in intergenic regions. Short deletions in cases were also modestly enriched for high-confidence enhancer regions. No individual CNV achieved thresholds for suggestive or significant association after genome-wide correction. p values < .01 were observed for 15q11.2 duplications (TUBGCP5, CYFIP1, NIPA1, and NIPA2), deletions in or near PRKN or MSR1, and exonic duplications of ATG5. ConclusionsThe increased burden of short deletions in patients with MDD suggests that rare CNVs increase the risk of MDD by disrupting regulatory regions. Results for longer deletions were less clear, but no large effects were observed for long multigenic CNVs (as seen in schizophrenia and autism). Further studies with larger sample sizes are warranted.