Genome-wide Copy Number Variation Study Research Articles

Exonic deletions in IMMP2L in schizophrenia with enhanced glycation stress subtype.

We previously identified a subtype of schizophrenia (SCZ) characterized by increased plasma pentosidine, a marker of glycation and oxidative stress (PEN-SCZ). However, the genetic factors associated with PEN-SCZ have not been fully clarified. We performed a genome-wide copy number variation (CNV) analysis to identify CNVs associated with PEN-SCZ to provide an insight into the novel therapeutic targets for PEN-SCZ. Plasma pentosidine was measured by high-performance liquid chromatography in 185 patients with SCZ harboring rare CNVs detected by array comparative genomic hybridization. In three patients with PEN-SCZ showing additional autistic features, we detected a novel deletion at 7q31.1 within exons 2 and 3 of IMMP2L, which encodes the inner mitochondrial membrane peptidase subunit 2. The deletion was neither observed in non-PEN-SCZ nor in public database of control subjects. IMMP2L is one of the SCZ risk loci genes identified in a previous SCZ genome-wide association study, and its trans-populational association was recently described. Interestingly, deletions in IMMP2L have been previously linked with autism spectrum disorder. Disrupted IMMP2L function has been shown to cause glycation/oxidative stress in neuronal cells in an age-dependent manner. To our knowledge, this is the first genome-wide CNV study to suggest the involvement of IMMP2L exons 2 and 3 in the etiology of PEN-SCZ. The combination of genomic information with plasma pentosidine levels may contribute to the classification of biological SCZ subtypes that show additional autistic features. Modifying IMMP2L functions may be useful for treating PEN-SCZ if the underlying biological mechanism can be clarified in further studies.

P.739 High-throughput proteomics for suicide biomarker discovery in patients with bipolar disorder

The genetic basis of suicide attempts (SA) remains unclear. Especially the role of copy number variations (CNVs) remains to be elucidated. The present study aimed to identify susceptibility variants associated with SA among Chinese with major depressive disorder (MDD), covering both CNVs and single-nucleotide polymorphisms (SNPs).We conducted a genome-wide association study (GWAS) on MDD patients with and without SA and top results were tested in a replication study. A genome-wide CNV study was also performed. Subsequently, a validation assay using qRT-PCR technology was performed to confirm any associated CNVs and then applied to the entire cohort to examine the association.Although GWAS did not identify any SNPs reaching genome-wide significance, we identified TPH2 as the top susceptibility gene (p-value = 2.75e−05) in gene-based analysis, which is a strong biological candidate for its role in the serotonergic system. As for CNV analysis, we found that the global rate of CNV was higher in SA than that in non-SA subjects (p-value = 0.023). Genome-wide CNV study revealed an SA-associated CNV region that achieved genome-wide significance (corrected p-value = 0.014). The associated CNV was successfully validated with a more rigorous qRT-PCR assay and identified to be a common variant in this cohort. Its deletion rate was higher in SA subjects [OR = 2.05 (1.02–4.12), adjusted p-value = 0.045]. Based on the GTEx database, genetic variants that probed this CNV were significantly associated with the expression level of ZNF33B in two brain regions (p-value < 4.2e−05). In stratified analysis, the CNV showed a significant effect [OR = 2.58 (1.06–6.27), p-value = 0.039] in those with high neuroticism but not in those with average or low neuroticism.We identified a new common CNV likely involved in the etiology of SA. This finding sheds light on an important role of common CNVs in the pathophysiology of SA, suggesting a new promising avenue for investigating its genetic architecture.

Genome-wide copy number variation-, validation- and screening study implicates a new copy number polymorphism associated with suicide attempts in major depressive disorder.

The genetic basis of suicide attempts (SA) remains unclear. Especially the role of copy number variations (CNVs) remains to be elucidated. The present study aimed to identify susceptibility variants associated with SA among Chinese with major depressive disorder (MDD), covering both CNVs and single-nucleotide polymorphisms (SNPs). We conducted a genome-wide association study (GWAS) on MDD patients with and without SA and top results were tested in a replication study. A genome-wide CNV study was also performed. Subsequently, a validation assay using qRT-PCR technology was performed to confirm any associated CNVs and then applied to the entire cohort to examine the association. Although GWAS did not identify any SNPs reaching genome-wide significance, we identified TPH2 as the top susceptibility gene (p-value=2.75e-05) in gene-based analysis, which is a strong biological candidate for its role in the serotonergic system. As for CNV analysis, we found that the global rate of CNV was higher in SA than that in non-SA subjects (p-value=0.023). Genome-wide CNV study revealed an SA-associated CNV region that achieved genome-wide significance (corrected p-value=0.014). The associated CNV was successfully validated with a more rigorous qRT-PCR assay and identified to be a common variant in this cohort. Its deletion rate was higher in SA subjects [OR=2.05 (1.02-4.12), adjusted p-value=0.045]. Based on the GTEx database, genetic variants that probed this CNV were significantly associated with the expression level of ZNF33B in two brain regions (p-value<4.2e-05). In stratified analysis, the CNV showed a significant effect [OR=2.58 (1.06-6.27), p-value=0.039] in those with high neuroticism but not in those with average or low neuroticism. We identified a new common CNV likely involved in the etiology of SA. This finding sheds light on an important role of common CNVs in the pathophysiology of SA, suggesting a new promising avenue for investigating its genetic architecture.

Copy number variation in the susceptibility to systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component and etiology characterized by chronic inflammation and autoantibody production. The purpose of this study was to ascertain copy number variation (CNV) in SLE using a case-control design in an admixed Brazilian population. The whole-genome detection of CNV was performed using Cytoscan HD array in SLE patients and healthy controls. The best CNV candidates were then evaluated by quantitative real-time PCR in a larger cohort or validated using droplet digital PCR. Logistic regression models adjusted for sex and ancestry covariates was applied to evaluate the association between CNV with SLE susceptibility. The data showed a synergistic effect between the FCGR3B and ADAM3A loci with the presence of deletions in both loci significantly increasing the risk to SLE (5.9-fold) compared to the deletion in the single FCGR3B locus (3.6-fold). In addition, duplications in these genes were indeed more frequent in healthy subjects, suggesting that high FCGR3B/ADAM3A gene copy numbers are protective factors against to disease development. Overall, 21 rare CNVs were identified in SLE patients using a four-step pipeline created for identification of rare variants. Furthermore, heterozygous deletions overlapping the CFHR4, CFHR5 and HLA-DPB2 genes were described for the first time in SLE patients. Here we present the first genome-wide CNV study of SLE patients in a tri-hybrid population. The results show that novel susceptibility loci to SLE can be found once the distribution of structural variants is analyzed throughout the whole genome.

Association between genome-wide copy number variation and arsenic-induced skin lesions: a prospective study

BackgroundExposure to arsenic in drinking water is a global health problem and arsenic-induced skin lesions are hallmark of chronic arsenic toxicity. We and others have reported germline genetic variations as risk factors for such skin lesions. The role of copy number variation (CNV) in the germline DNA in this regard is unknown.MethodsFrom a large prospectively followed-up cohort, exposed to arsenic, we randomly selected 2171 subjects without arsenic-induced skin lesions at enrollment and genotyped their whole blood DNA samples on Illumina Cyto12v2.1 SNP chips to generate DNA copy number. Participants were followed up every 2 years for a total of 8 years, especially for the development of skin lesions. In Cox regression models, each CNV segment was used as a predictor, accounting for other potential covariates, for incidence of skin lesions.ResultThe presence of genomic deletion(s) in a number of genes (OR5J2, GOLGA6L7P, APBA2, GALNTL5, VN1R31P, PHKG1P2, SGCZ, ZNF658) and lincRNA genes (RP11-76I14.1, CTC-535 M15.2, RP11-73B2.2) were associated with higher risk [HR between 1.67 (CI 1.3-2.1) and 2.15 (CI 1.5-2.9) for different CNVs] for development of skin lesions independent of gender, age, and arsenic exposure. Some deletions had stronger effect in a specific gender (ZNF658 in males, SGCZ in females) and some had stronger effect in higher arsenic exposure (lincRNA CTD-3179P9.1) suggesting a possible gene-environment interaction.ConclusionThis first genome-wide CNV study in a prospectively followed-up large cohort, exposed to arsenic, suggests that DNA deletion in several genes and lincRNA genes may predispose an individual to a higher risk of development of arsenic-induced skin lesions.

Identification of a Copy Number Variation on Chromosome 20q13.12 Associated with Osteoporotic Fractures in the Korean Population.

Osteoporotic fractures (OFs) are critical hard outcomes of osteoporosis and are characterized by decreased bone strength induced by low bone density and microarchitectural deterioration in bone tissue. Most OFs cause acute pain, hospitalization, immobilization, and slow recovery in patients and are associated with increased mortality. A variety of genetic studies have suggested associations of genetic variants with the risk of OF. Genome-wide association studies have reported various single-nucleotide polymorphisms and copy number variations (CNVs) in European and Asian populations. To identify CNV regions associated with OF risk, we conducted a genome-wide CNV study in a Korean population. We performed logistic regression analyses in 1,537 Korean subjects (299 OF cases and 1,238 healthy controls) and identified a total of 8 CNV regions significantly associated with OF (p < 0.05). Then, one CNV region located on chromosome 20q13.12 was selected for experimental validation. The selected CNV region was experimentally validated by quantitative polymerase chain reaction. The CNV region of chromosome 20q13.12 is positioned upstream of a family of long non-coding RNAs, LINC01260. Our findings could provide new information on the genetic factors associated with the risk of OF.

Genome-wide copy number variation study reveals KCNIP1 as a modulator of insulin secretion

Copy number variations (CNVs) have emerged as another important genetic marker in addition to SNP for understanding etiology of complex diseases. In light of this, we performed a genome-wide CNV study to identify type 2 diabetes (T2D)-associated CNV using an array comparative genomic hybridization from 3180 subjects for T2D cases (n=863) and controls (n=2,317). Thus, five CNV regions having a p-value threshold ≤0.05 were identified and evaluated by validation with quantitative PCR and comparison with previously reported CNV regions in the Database of Genomic Variants. Furthermore, we performed a functional experiment to assess the biological significance of a gene encompassing a CNV region. The inhibition of KCNIP1 led to increased insulin secretion in a glucose-dependent manner, but had no effect on insulin gene transcription as well as cell apoptosis. Taken together, these data indicate that KCNIP1 from CNV study might function as a T2D-susceptibility gene whose dysregulation alters insulin production.

Adult neuropsychiatric expression and familial segregation of 2q13 duplications.

New genomic disorders associated with large, rare, recurrent copy number variations (CNVs) are being discovered at a rapid pace. Detailed phenotyping and family studies are rare, however, as are data on adult phenotypic expression. Duplications at 2q13 were recently identified as risk factors for developmental delay/autism and reported in the prenatal setting, yet few individuals (all children) have been extensively phenotyped. During a genome-wide CNV study of schizophrenia, we identified two unrelated probands with 2q13 duplications. In this study, detailed phenotyping and genotyping using high-resolution microarrays was performed for 12 individuals across their two families. 2q13 duplications were present in six adults, and co-segregated with clinically significant later-onset neuropsychiatric disorders. Convergent lines of evidence implicated GABAminergic dysfunction. Analysis of the genic content revealed promising candidates for neuropsychiatric disease, including BCL2L11, ANAPC1, and MERTK. Intrafamilial genetic heterogeneity and "second hits" in one family may have been the consequence of assortative mating. Clinical genetic testing for the 2q13 duplication and the associated genetic counseling was well received. In summary, large rare 2q13 duplications appear to be associated with variable adult neuropsychiatric and other expression. The findings represent progress toward clinical translation of research results in schizophrenia. There are implications for other emerging genomic disorders where there is interest in lifelong expression.

Genome-wide copy number variation study and gene expression analysis identify ABI3BP as a susceptibility gene for Kashin–Beck disease

Kashin-Beck disease (KBD) is a chronic osteochondropathy. In this study, we conducted the first genome-wide copy number variation study (GCNVS) of KBD totally involving 2,743 Chinese Han adults. GCNVS was first performed using Affymetrix Human SNP6.0 Arrays. The identified copy number variations (CNVs) were then replicated in an independent Chinese Han sample containing 1,026 subjects. SNP genotyping, CNV identification and quality control were implemented by Birdsuite. STRUCTURE and EIGENSTRAT were applied for controlling potential population stratification in the GCNVS. Association analysis was conducted using PLINK. Microarray and qRT-PCR were also conducted to compare the expression levels of the genes overlapping with identified CNVs between KBD patients and healthy controls. GCNVS found that CNV452 (P value = 7.78 × 10(-5)) overlapping with ABI3BP gene was significantly associated with KBD. Replication association study observed that rs9850273 (P value = 0.008) and rs7613610 (P value = 0.021) in ABI3BP gene were significantly associated with KBD. Gene expression analysis also found that ABI3BP was up-regulated in KBD patients compared to healthy controls. Our results suggest that ABI3BP was a novel susceptibility gene for KBD.

Large scale variation in DNA copy number in chicken breeds

BackgroundDetecting genetic variation is a critical step in elucidating the molecular mechanisms underlying phenotypic diversity. Until recently, such detection has mostly focused on single nucleotide polymorphisms (SNPs) because of the ease in screening complete genomes. Another type of variant, copy number variation (CNV), is emerging as a significant contributor to phenotypic variation in many species. Here we describe a genome-wide CNV study using array comparative genomic hybridization (aCGH) in a wide variety of chicken breeds.ResultsWe identified 3,154 CNVs, grouped into 1,556 CNV regions (CNVRs). Thirty percent of the CNVs were detected in at least 2 individuals. The average size of the CNVs detected was 46.3 kb with the largest CNV, located on GGAZ, being 4.3 Mb. Approximately 75% of the CNVs are copy number losses relatively to the Red Jungle Fowl reference genome. The genome coverage of CNVRs in this study is 60 Mb, which represents almost 5.4% of the chicken genome. In particular large gene families such as the keratin gene family and the MHC show extensive CNV.ConclusionsA relative large group of the CNVs are line-specific, several of which were previously shown to be related to the causative mutation for a number of phenotypic variants. The chance that inter-specific CNVs fall into CNVRs detected in chicken is related to the evolutionary distance between the species. Our results provide a valuable resource for the study of genetic and phenotypic variation in this phenotypically diverse species.

Genome-wide copy number variation study in anorectal malformations

Anorectal malformations (ARMs, congenital obstruction of the anal opening) are among the most common birth defects requiring surgical treatment (2-5/10 000 live-births) and carry significant chronic morbidity. ARMs present either as isolated or as part of the phenotypic spectrum of some chromosomal abnormalities or monogenic syndromes. The etiology is unknown. To assess the genetic contribution to ARMs, we investigated single-nucleotide polymorphisms and copy number variations (CNVs) at genome-wide scale. A total of 363 Han Chinese sporadic ARM patients and 4006 Han Chinese controls were included. Overall, we detected a 1.3-fold significant excess of rare CNVs in patients. Stratification of patients by presence/absence of other congenital anomalies showed that while syndromic ARM patients carried significantly longer rare duplications than controls (P = 0.049), non-syndromic patients were enriched with both rare deletions and duplications when compared with controls (P = 0.00031). Twelve chromosomal aberrations and 114 rare CNVs were observed in patients but not in 868 controls nor 11 943 healthy individuals from the Database of Genomic Variants. Importantly, these aberrations were observed in isolated ARM patients. Gene-based analysis revealed 79 genes interfered by CNVs in patients only. In particular, we identified a de novo DKK4 duplication. DKK4 is a member of the WNT signaling pathway which is involved in the development of the anorectal region. In mice, Wnt disruption results in ARMs. Our data suggest a role for rare CNVs not only in syndromic but also in isolated ARM patients and provide a list of plausible candidate genes for the disorder.

Genome-Wide Copy Number Variation Association Analyses for Age at Menarche

Menarche is a significant physiological event for women. Age at menarche (AAM) is a heritable trait associated with many common female diseases. The genetic basis and the mechanism for AAM are largely unknown. Copy number variation (CNV) is a common type of genetic variation underlying human complex traits. The importance of CNV to AAM variation is unclear. The objective of the study was to identify CNV important to AAM variation. We performed the first genome-wide CNV study of AAM in 1654 Caucasian females using Affymetrix human single-nucleotide polymorphism 6.0 array. We also replicated our findings in another Chinese cohort containing 752 women. We identified a CNV, variation_38399, in the 2q14.2 region, for association with AAM (P = 1.03 × 10(-3)). The CNV has two variants (one copy and two copy), with a mean AAM of 14.00 yr and 12.90 yr, respectively. Interestingly, in a Chinese sample containing 752 women, this CNV has been replicated both with a marginally significant P = 0.090 and with a same direction of effect (a lower copy number for a later AAM). The CNV is located approximately 75 kb upstream of the diazepam binding inhibitor (DBI), a gene known to regulate estrogen levels, a key factor for menarche. Our findings for the first time identified a novel CNV and suggested the DBI-mediated endocrinological pathway as a potential mechanism for AAM regulation.

Genome-wide copy number variation study associates metabotropic glutamate receptor gene networks with attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a common, heritable neuropsychiatric disorder of unknown etiology. We performed a whole-genome copy number variation (CNV) study on 1,013 cases with ADHD and 4,105 healthy children of European ancestry using 550,000 SNPs. We evaluated statistically significant findings in multiple independent cohorts, with a total of 2,493 cases with ADHD and 9,222 controls of European ancestry, using matched platforms. CNVs affecting metabotropic glutamate receptor genes were enriched across all cohorts (P = 2.1 × 10(-9)). We saw GRM5 (encoding glutamate receptor, metabotropic 5) deletions in ten cases and one control (P = 1.36 × 10(-6)). We saw GRM7 deletions in six cases, and we saw GRM8 deletions in eight cases and no controls. GRM1 was duplicated in eight cases. We experimentally validated the observed variants using quantitative RT-PCR. A gene network analysis showed that genes interacting with the genes in the GRM family are enriched for CNVs in ∼10% of the cases (P = 4.38 × 10(-10)) after correction for occurrence in the controls. We identified rare recurrent CNVs affecting glutamatergic neurotransmission genes that were overrepresented in multiple ADHD cohorts.