Open Access

A replication study of schizophrenia-related rare copy number variations in a Han Southern Chinese population

  • Jianmin Yuan1,
  • Jianlin Hu2,
  • Zhiqiang Li3,
  • Fuquan Zhang1,
  • Dexiang Zhou1 and
  • Chunhui Jin1Email author
Contributed equally
Hereditas2017154:2

https://doi.org/10.1186/s41065-016-0025-x

Received: 5 July 2016

Accepted: 26 December 2016

Published: 14 January 2017

Abstract

Background

Schizophrenia (SCZ) is a common, complex and severe psychiatric disorder associated with many different genetic and environmental risk factors. Evidence from genetic studies has revealed the role of genome structural variations, specifically copy number variants (CNVs), in the etiology of SCZ. Nevertheless, the occurrence of CNVs and their relation to SCZ has remained relatively unstudied in the diverse Han Chinese population.

Results

We used a case/control paradigm, including 476 cases and 1023 controls. All samples were genotyped using the Axiom® Exome Genotyping Arrays. Four CNVs, including two deletions and two duplications, were detected in this study. Notably, the 16p11.2 duplication from 29.3 Mb to 29.6 Mb was detected in four cases (0.84%) and one control (0.098%) (p = 0.0377).

Conclusions

The results highlight the potential role of these deletions and duplications in the development of SCZ. Clearly, larger sample sized studies are needed for a careful localization of these CNVs and to possibly detect more deletions and/or duplications, associated with the development of SCZ in the Han Chinese population.

Keywords

SchizophreniaCopy number variantsHan Chinese

Background

The role of genome structural variations, specifically copy number variants (CNVs), in the development of mental disorders is a growing area of study. Due to the potential polygenic effects and the capacity to alter gene dosage, CNVs may have dramatic consequences on expression levels and therefore phenotype. Szatkiewicz et al hypothesize that CNV deletions are more severe compared to CNV duplications, due to the loss of function effect [1]. Furthermore, large CNVs are more severe than small CNVs, which affect fewer genes and regulatory regions. Because of their direct impact on gene modification, CNVs play a major role in evolution and have the potential to disturb some genomically unstable areas. For instance, the 1q21.1 region is associated with neurological development and thus could have impacted brain evolution. Therefore, deletion of 1q21.1 could explain the neurological and psychiatric conditions observed in patients carrying this CNV. The larger expression of the genome by brain tissue compared to other tissues explains the rationale for selecting CNVs as likely candidates for modulating psychiatric illnesses. As a result of the cost associated with whole genome sequencing, exome arrays have emerged as a viable option for the detection of both single nucleotide polymorphisms (SNPs) and CNVs. Only protein-coding DNA will be examined, thereby providing an approximate portrayal of the genomic expression. Exome arrays have both confirmed SNPs previously established by genome wide association studies (GWAS), as well as identified novel SNPs associated with a wide-range of diseases. Exome sequencing has been utilized for CNV analysis in the genome analysis of Swedish patients with schizophrenia (SCZ) and has confirmed some of the previously determined CNVs from GWAS [1]. Nevertheless, the diverse Han Chinese population represents a major group that has remained relatively unstudied with regards to the occurrence of CNVs and their relation to mental health disorders, such as SCZ. Several previous studies implicated numerous rare CNVs in SCZ, including the 1q21.1 deletion [2], 15q11.2 deletion [3], Williams-Beuren syndrome (WBS) region duplication [4] and 16p11.2 duplication [1]. These new findings raise fundamental clinical and scientific questions concerning classification of major neuropsychiatric disorders, modes of inheritance, diagnostics, and genetic counseling. None of these CNVs have been investigated in Chinese Han population. In this paper we analyze the exome SNP array data in 476 SCZ cases and 1023 healthy controls taken from southern China for the presence of CNVs, particularly focusing on the four CNVs listed above.

Results and Discussion

The 16p11.2 duplication from 29.3 Mb to 29.6 Mb was detected in four cases (0.84%) and one control (0.098%) (p = 0.0377, Fisher’s exact test) (Table 1). A previous investigation by McCarthy et al [5] reported an association of an adjacent duplication from 29.6 Mb to 30.2 Mb within 16p11.2 with SCZ (p = 1.2 × 10−5 for the initial cohort, p = 0.022 in the replication cohort). A subsequent meta-analysis from a follow-up study [5] involving 8,590 individuals with schizophrenia, 2,172 with developmental delay or autism spectrum disorders (ASD), 4,822 with bipolar disorder, and 30,492 controls reported a p value of 4.8 × 10−7 for association with SCZ. The authors also reported an association between the 16p11.2 duplication and both bipolar disorder and autism. Alternatively, the 16p11.2 deletion has been associated with autism and developmental disorders, suggesting an overlapping etiology between SCZ, ASD and other developmental disorders. The 16p11.2 duplication is also associated with a decrease in brain size compared to 16p11.2 deletion, and it is worth noting the possible connection between autism and increased brain volumes, not only in the content of McCarthy et al, but also in other studies.
Table 1

The significant level of the identified rare CNVs

Locus

Position in Mb (hg19)

Cases (N = 476)

Control (N = 1023)

OR

Fisher's P

Cases from previous studiesa (%)

Controls from previous studiesa (%)

Frequency in Chinese populationb (%)

1q21.1 del

chr1:146.5–147.8

1

0

NA

0.3175

0.170

0.016

0.092

15q11.2 del

chr15:22.7–23.1

1

1

2.15

0.5344

0.640

0.410

0.254

WBS dup

chr7:72.7–74.3

2

0

NA

0.1007

0.044

0.016

0.016

16p11.2 dup

chr16:29.3–29.6

4

1

8.66

0.0377

0.059

0.000

0.087

aResults from previous reported European data sets (Br J Psychiatry. 2014;204(2):108–114). bFrequencies of the CNVs in Chinese population, which were established from the Chinese study (Biol Psychiatry. 2016;80(4):331–337)

One out of 476 cases (0.21%) carried the 1q21.1 deletion, which was not expressed in any of the controls (p = 0.3175, Fisher’s exact test) (Table 1). The 1q21.1 deletion spans 1.3 Mb (chromosome 1: 146.5 Mb- 147.8 Mb) in the distal region of 1q21.1, not containing the TAR region. This deletion is also well known to be associated with SCZ, although the deletion region in our results indicated a much smaller area of deletion than Stephansson et al. described [68] The associated region is gene rich, and contains the gene BCL9, and when mutated, BCL9 has been shown by Li et al. to be associated with schizophrenia, as well as bipolar and major depression in the Chinese Han population [2]. However, the gene function for BCL9 is largely unknown. Besides schizophrenia, the highly unstable 1q21.1 region has also been implicated in a host of other human disorders, among which are idiopathic mental retardation [9, 10], autism/ASD [9], congenital heart disease [11], microcephaly/macrocephaly [9, 12] and neuroblastoma [13]. Dumas and Sikela indicated that the 1q21.1 CNV region contains the vast majority of DUF1220 protein domains, which portray a striking human lineage with a specific increase in copy number relative to the other great apes [14]. The increase in DUF1220 roughly correlates to increased brain size, and thus it is relevant to point out that deletions in 1q21.1 may also be associated with microcephaly, possibly due to a decrease in DUF1220. This trait is a common phenotype for schizophrenia as opposed to autism [15].

The 15q11.2 deletion from 22.7 Mb to 23.1 Mb occurred in one case (0.21%) and one control (0.098%) (p = 0.5344, Fisher’s exact test) (Table 1). The 15q11.2 deletion is one of the regions deleted in the Prader-Willi syndrome, whose symptoms include obsessive compulsive behavior, adaptive skills difficulty, and low intellectual ability. Of note, the Prader-Willi deletion in the Prader-Willi syndrome is much larger than that we identified. In addition, the 15q11.2 deletion has been reported in three different studies to be associated with SCZ [4]. The most relevant gene affected in the 15q11.2 region, CYFIP1 interacts with the fragile X mental retardation protein (FMRP) as well as the Rho GTPase Rac1, which is involved in the maintenance of neuronal structures [16].

We also observed that the Williams-Beuren syndrome (WBS) region from 22.7 Mb to 23.1 Mb on chromosome 7 was duplicated in two cases (0.42%), but not in controls (Table 1), for a p-value of 0.1007. The association between the 7q11.23 duplication and schizophrenia was previously reported by Mulle et al., [17] as well as Kirov et al [4]. Despite the fact that the specific biological pathways are unknown, the SCZ patients with WBS duplication also all had social anxiety and language delay prior to the onset of SCZ.

Take all the four highlighted CNVs as a whole, 1.68% of our cases have a pathogenic CNV, it is about eight times higher than that in the control group (0.20%). The detailed locations are given in Fig. 1 and the Additional file 1: Table S2.
Fig. 1

Schematic illustration of the chromosomal locations of the rare CNVs using UCSC Genome Browser. a 1q21.1 del; b 15q11.2 del; c 16p11.2 dup; d WBS dup. Note: In custom tracks, red color indicates deletion in case, dark red indicates deletion in control. For duplication bright blue color indicates duplication in case, while dark blue indicates duplication in control

Conclusions

To summarize, this study reported four CNV findings, including two deletions and two duplications. The results highlight the potential role of these deletions and duplications in the development of SCZ. Clearly, larger sample sized studies are needed for a careful localization of these CNVs and to possibly detect more deletions and/or duplications, associated with the development of SCZ in the Han Chinese population.

Methods

We used a case/control paradigm, including 476 cases and 1023 controls (Additional file 1: Table S1). Before study enrollment, all SCZ patients and healthy subjects were required to sign a consent form. Diagnosis of SCZ was confirmed by interviews with at least two experienced psychiatrists in accordance with criteria in the Diagnostic and Statistical Manual of Mental Disorders, Fourth edition (DSM-IV). Exclusion criteria included the presence of other mood or neurodevelopmental disorders, epilepsy, or mental retardation. Control subjects were recruited by advertisement and selected based on the results from a questionnaire, concerning history of obstetric complications, substantial head injury, seizures, neurological or psychiatric disease, hypertension, diabetes, and substance abuse, as well as family medical histories. Subjects with any personal or family history of psychiatric illness among their first-degree relatives were excluded. The patients were recruited from several hospitals in Jiangsu province. All patients and healthy controls were age-matched and recruited from the cities of Wuxi and Nanjing in the Jiangsu province. Individuals not born in Jiangsu or whose family members were not born in Jiangsu were excluded from the study based on self-reported data. In addition, there was no blood relation between patients. All participants are Han Chinese. All samples were genotyped using the Axiom® Exome Genotyping Arrays. We first calculated the quality control measures (Dish-QC, DQC), and arrays with DQC < 0.82 (n = 26) signified poor quality, and could result in false-positive CNVs, so they were excluded from the study. We then called CNV for the left samples using PennCNV [18]. The genotyping calls and signal intensity data were generated using Affymetrix Power Tools, and then converts into Log R Ratio (LRR) and B Allele Frequency (BAF) values. To minimize the noise of genome-wide intensity signals, only samples within the s.d. of the normalized intensity (LRR < 0.35) were included. Wave artifacts roughly correlating with GC content resulting from hybridization bias of low full-length DNA quantity are known to interfere with the accurate inference of CNVs. Only samples where |GC base pair wave factor (GCWF)| < 0.05 were accepted. PLINK's identity-by-descent (IBD) analysis was used to detect cryptic relatedness based on the common variants (minor allele frequency > 1%) [19]. When a pair of individuals exhibited a PI_HAT > 0.2, the member of the pair with the lower call rate was excluded from the analysis (n = 3). The clean_cnv.pl program was adopted to merge adjacent calls of the same individual using the fraction threshold (0.5, calculated as base pair length). CNVs were then excluded if they were covered by < 5 probes or had a probe density > 1 probe/40 kb. If the count of CNV calls made by PennCNV exceeds 70, the DNA quality was assumed to be poor. Therefore, only samples with CNV call count <70 were included. After quality control, 476 cases and 1023 controls were kept for further analysis.

More than ten CNVs have been convincingly associated with SCZ [20, 21], and most of the previous studies are mainly carried out in Caucasian populations. We try to replicate the reported CNVs identified in Caucasian populations using a Han Chinese sample.

Four CNVs were found to be overrepresented in our cases and further evaluated as a result of our study, the 1q21.1 deletion, the 15q11.2 deletion, the WBS duplication and the 16p11.2 duplication. We compared the frequencies in the cases and controls groups for each particular region, and a Fisher's exact test was used to test the significance of the difference. All of these have been previously found to be associated with schizophrenia by multiple studies.

Abbreviations

ASD: 

Autism spectrum disorder

CNV: 

Copy number variant

GWAS: 

Genome wide association study

SCZ: 

Schizophrenia

SNP: 

Single nucleotide polymorphism

WBS: 

Williams-Beuren syndrome

Declarations

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No.81301147), the Medical Research Pharmaceutical Program of Wuxi Hospital Management Center (YGZXY1306), and the Science & Technology Development Foundation of Wuxi Municipal Science & Technology Bureau (CSZON1316).

Availability of data and materials

All data generated or analyzed during this study are included in this published article [and its Additional file 1].

Authors’ contributions

JY and CJ conceived and designed this study. JY, JH, ZL, FZ and DZ performed the experiments. ZL and FZ contributed reagents/materials/analysis tools. JY and CJ wrote the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The clinical study was in compliance with the World Medical Association (WMA) Declaration of Helsinki - Ethical Principles for Medical Research Involving Human Subjects. The study was approved by the Ethics Committees of the Wuxi Health Mental Center.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Wuxi Mental Health Center, Nanjing Medial University
(2)
Wuxi Second Traditional Chinese Medicine Hospital
(3)
Bio-X Institutes, Key Laboratory for the Genetics of Developmental and Neuropsychiatric Disorders (Ministry of Education), Shanghai Jiao Tong University

References

  1. Szatkiewicz JP, Neale BM, O'Dushlaine C, Fromer M, Goldstein JI, Moran JL, et al. Detecting large copy number variants using exome genotyping arrays in a large Swedish schizophrenia sample. Mol Psychiatry. 2013;18(11):1178–84.View ArticlePubMedPubMed CentralGoogle Scholar
  2. Li J, Zhou G, Ji W, Feng G, Zhao Q, Liu J, et al. Common variants in the BCL9 gene conferring risk of schizophrenia. Arch Gen Psychiatry. 2011;68(3):232–40.View ArticlePubMedGoogle Scholar
  3. Zhao Q, Li T, Zhao X, Huang K, Wang T, Li Z, et al. Rare CNVs and tag SNPs at 15q11.2 are associated with schizophrenia in the Han Chinese population. Schizophr Bull. 2013;39(3):712–9.View ArticlePubMedGoogle Scholar
  4. Kirov G, Pocklington AJ, Holmans P, Ivanov D, Ikeda M, Ruderfer D, et al. De novo CNV analysis implicates specific abnormalities of postsynaptic signalling complexes in the pathogenesis of schizophrenia. Mol Psychiatry. 2012;17(2):142–53.View ArticlePubMedGoogle Scholar
  5. McCarthy SE, Makarov V, Kirov G, Addington AM, McClellan J, Yoon S, et al. Microduplications of 16p11.2 are associated with schizophrenia. Nat Genet. 2009;41(11):1223–7.View ArticlePubMedPubMed CentralGoogle Scholar
  6. Stone JL, O'Donovan MC, Gurling H, Kirov GK, Blackwood DH, Corvin A, et al. Rare chromosomal deletions and duplications increase risk of schizophrenia. Nature. 2008;455(7210):237–41.View ArticleGoogle Scholar
  7. Levinson DF, Shi J, Wang K, Oh S, Riley B, Pulver AE, et al. Genome-wide association study of multiplex schizophrenia pedigrees. Am J Psychiatry. 2012;169(9):963–73.View ArticlePubMedGoogle Scholar
  8. Stefansson H, Rujescu D, Cichon S, Pietilainen OP, Ingason A, Steinberg S, et al. Large recurrent microdeletions associated with schizophrenia. Nature. 2008;455(7210):232–6.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Mefford HC, Sharp AJ, Baker C, Itsara A, Jiang Z, Buysse K, et al. Recurrent rearrangements of chromosome 1q21.1 and variable pediatric phenotypes. N Engl J Med. 2008;359(16):1685–99.View ArticlePubMedPubMed CentralGoogle Scholar
  10. Sharp AJ, Hansen S, Selzer RR, Cheng Z, Regan R, Hurst JA, et al. Discovery of previously unidentified genomic disorders from the duplication architecture of the human genome. Nat Genet. 2006;38(9):1038–42.View ArticlePubMedGoogle Scholar
  11. Christiansen J, Dyck JD, Elyas BG, Lilley M, Bamforth JS, Hicks M, et al. Chromosome 1q21.1 contiguous gene deletion is associated with congenital heart disease. Circ Res. 2004;94(11):1429–35.View ArticlePubMedGoogle Scholar
  12. Brunetti-Pierri N, Berg JS, Scaglia F, Belmont J, Bacino CA, Sahoo T, et al. Recurrent reciprocal 1q21.1 deletions and duplications associated with microcephaly or macrocephaly and developmental and behavioral abnormalities. Nat Genet. 2008;40(12):1466–71.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Vandepoele K, Van Roy N, Staes K, Speleman F, van Roy F. A novel gene family NBPF: intricate structure generated by gene duplications during primate evolution. Mol Biol Evol. 2005;22(11):2265–74.View ArticlePubMedGoogle Scholar
  14. Dumas L, Sikela JM. DUF1220 domains, cognitive disease, and human brain evolution. Cold Spring Harb Symp Quant Biol. 2009;74:375–82.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Crespi B, Badcock C. Psychosis and autism as diametrical disorders of the social brain. Behav Brain Sci. 2008;31(3):241–61. discussion 61-320.PubMedGoogle Scholar
  16. Schenck A, Bardoni B, Moro A, Bagni C, Mandel JL. A highly conserved protein family interacting with the fragile X mental retardation protein (FMRP) and displaying selective interactions with FMRP-related proteins FXR1P and FXR2P. Proc Natl Acad Sci U S A. 2001;98(15):8844–9.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Mulle JG, Pulver AE, McGrath JA, Wolyniec PS, Dodd AF, Cutler DJ, et al. Reciprocal duplication of the Williams-Beuren syndrome deletion on chromosome 7q11.23 is associated with schizophrenia. Biol Psychiatry. 2014;75(5):371–7.View ArticlePubMedGoogle Scholar
  18. Wang K, Li M, Hadley D, Liu R, Glessner J, Grant SFA, et al. PennCNV: An integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Resh. 2007;17(11):1665–74.View ArticleGoogle Scholar
  19. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Rees E, Walters JTR, Georgieva L, Isles AR, Chambert KD, Richards AL, et al. Analysis of copy number variations at 15 schizophrenia-associated loci. Bri J Psychiatry. 2014;204(2):108–14.View ArticleGoogle Scholar
  21. Li Z, Chen J, Xu Y, Yi Q, Ji W, Wang P, et al. Genome-wide Analysis of the Role of Copy Number Variation in Schizophrenia Risk in Chinese. Biol Psychiatry. 2016;80(4):331–7.View ArticlePubMedGoogle Scholar

Copyright

© The Author(s) 2017

Advertisement