Deletion and duplication of 16p11.2 are associated with opposing effects on visual evoked potential amplitude
© The Author(s). 2016
Received: 11 February 2016
Accepted: 21 June 2016
Published: 27 June 2016
Duplication and deletion of the chromosomal region 16p11.2 cause a broad range of impairments, including intellectual disability, language disorders, and sensory symptoms. However, it is unclear how changes in 16p11.2 dosage affect cortical circuitry during development. The aim of this study was to investigate whether the visual evoked potential (VEP) could be used as a noninvasive quantitative measure of cortical processing in children with 16p11.2 copy number variation.
Pattern-reversal VEPs were successfully recorded in 19 deletion carriers, 9 duplication carriers, and 13 typically developing children between the ages of 3 and 14 years. The stimulus was a black and white checkerboard (60’) that reversed contrast at 2 Hz. VEP responses were extracted from continuous EEG recorded using a high-density elasticized electrode net.
Quantitative analysis of the VEP waveform revealed that, relative to controls, deletion carriers displayed increased amplitude and duplication carriers displayed diminished amplitude. Latencies of the VEP waveform components were unaffected by 16p11.2 status. P1 amplitude did not correlate with age, IQ, or head circumference.
The results of this study suggest that recording VEP is a useful method to assay cortical processing in children with 16p11.2 copy number variation. There is a gene dosage-dependent effect on P1 amplitude that merits further investigation. The VEP is directly translatable to animal models, offering a promising way to probe the neurobiological mechanisms underlying cortical dysfunction in this developmental disorder.
Keywords16p11.2 copy number variation Visual evoked potential Visual cortex
Copy number variation of the 16p11.2 chromosomal region (16p CNV) is associated with developmental delay and a wide range of physical abnormalities, behavioral problems, and intellectual deficits . Some features are similarly manifested in deletion and duplication carriers, such as seizures and cognitive impairment , and others are affected in opposite ways, including head circumference and body weight [2, 3]. This begs the question of whether brain function is affected in a gene dosage-dependent manner. Although progress is being made on characterizing the disorder and parsing the neurobiological mechanisms, there remains a gap in our understanding of how these genetic changes ultimately produce the varied cognitive and behavioral outcomes in this population.
Abnormal processing of sensory information is emerging as a common problem in neurodevelopmental disorders [4, 5], and testing sensory modalities like vision offers an easily accessible window into cortical processing. Abnormalities of basic sensory processing may contribute to and be indicative of global cortical deficits that would be expressed as intellectual and behavioral deficits. A recent study reported delayed auditory evoked responses using magnetoencephalography in deletion but not duplication carriers . We were interested in whether this deficit would also pertain to the processing of visual information. Vision is an ideal system because it is evolutionarily conserved between rodents and humans, well characterized, easily tested in a standardized manner, and relevant to socialization, communication, and cognition.
Recording pattern-reversal visual evoked potentials (VEPs) is a powerful method to noninvasively assess the integrity of cortical processing in babies and children, even those with profound disabilities [4, 6]. Eliciting VEP involves a quick and passive procedure that only requires fixation on a monitor, and the evoked response is stable and robust enough to produce a recognizable waveform after only a few trials. The waveform shape elicited by this particular large check stimulus matures within the first 2 years of life [7–11], and group differences in latency and amplitude can be quantified [4, 12].
Here, we recorded pattern-reversal VEP in children carrying 16p deletions and duplications. Our goals were to assess if the VEP could detect differences in cortical processing between typically developing children and those with 16p CNV and between those with deletions and duplications. We predicted that the latency and/or amplitude of visual responses would either be similarly affected by 16p deletion and duplication due to a common effect of disruption of this chromosomal region on synapses and circuitry, or would be differentially affected by 16p copy number due to a gene dosage effect on the wiring of the visual cortex [5, 13, 14].
All procedures were reviewed and approved by the Boston Children’s Hospital (BCH) Office of Clinical Investigation and the BCH Institutional Review Board, and written informed consent was obtained from each participant or guardian prior to testing. Individuals carrying 16p11.2 deletions and duplications were recruited to participate in the study through the Simons Variation in Individuals Project (SVIP). Eligibility and exclusionary criteria are described previously . Typically developing children were recruited through the BCH participant registry and did not have any neurological or developmental disorders, seizures, and ophthalmological problems. Data were collected onsite at BCH and at SVIP family meetings in Orlando, Florida, and Chicago, Illinois.
Description of study participants
Age (months) median, range
Phenotypic information was obtained from the Simons Foundation Autism Research Initiative (SFARI) Base (www.sfari.org/resources/sfari-base). Full-scale intelligence quotient (IQ) scores were derived without adjusting for prematurity and were only used if taken within 6 months of the VEP session date. IQ was not assessed for typically developing participants because they were recruited independently through BCH and did not undergo the SFARI phenotypic battery. However, because the typically developing sample was drawn from a database of nearly 20,000 families and we have previously assessed IQ in a subset of these children as part of another project, we expect the IQ of our current group is average to above average. Autism spectrum disorder status was determined based on the diagnosis given by experienced, licensed clinicians. Seizure history was available for 11/19 deletion carriers and 6/9 duplication carriers through the SFARI neurologic record review form. Of these individuals, 4 deletion carriers and no duplication carriers had a history of seizures. Current medication information was available for 7/19 deletion carriers and 4/9 duplication carriers through the SFARI medication questionnaire. Of these individuals, 2 deletion carriers and no duplication carriers were currently taking antiepileptic drugs at the time of the VEP session.
Participants were not given an ophthalmological exam at the time of the VEP session due to budget and time constraints, but information related to eye surgeries, cataract, coloboma, glaucoma, and corrected vision was available on the SFARI developmental and medical history form for 17/19 deletion carriers and 8/9 duplication carriers. A total of 7 deletion carriers and 3 duplication carriers were reported to have corrected vision but no other uncorrected visual problems.
Visual evoked potential recordings
Pattern VEPs were recorded as described previously . Participants sat 60 cm from a Tobii Eye Tracker monitor (Tobii Technology, Sweden). Phase reversing 99 % contrast black and white checkerboard patterns were presented on the monitor using ePrime software (Psychology Software Tools Inc., Pittsburgh, PA, USA). The diagonal check size was 0.5 cycles/degree (60’), and the reversal rate was once every 500 ms. Approximately 50–100 trials were presented to each participant, and VEP sessions lasted no longer than several minutes. Stimulus presentation was experimenter-driven based on verification of binocular fixation on the monitor.
Continuous electroencephalogram (EEG) was recorded from a 128-channel HydroCel Geodesic Sensor Net (Electrical Geodesics Inc., Eugene, OR, USA), referenced to Cz, amplified with a NetAmps 300 amplifier, and digitized at 500 Hz. Data were analyzed offline using NetStation 4.5.4 software. The signal was filtered (0.3–30 Hz), segmented into 300 ms post-stimulus recording periods, and baseline-corrected against the mean voltage during the 100 ms pre-stimulus period. A 34-ms offset that was consistent for all participants was accounted for at this point in the analysis (16 ms between the stimulus trigger and appearance on the monitor and 18 ms due to anti-aliasing filters within the Electrical Geodesics Inc. NetAmps 300 amplifier).
Quantification of VEP components and comparison between groups
H(2): 0.22, p = 0.894
Latency measurements (ms)
H(2): 6.85, p = 0.033
Ctr/Del: p = 0.053
Ctr/Dup: p = 0.099
Del/Dup: p = 1.000
H(2): 1.72, p = 0.422
H(2): 0.13, p = 0.937
H(2): 0.37, p = 0.858
Amplitude measurements (uV)
H(2): 3.97, p = 0.137
H(2): 11.92, p = 0.003
Ctr/Del: p = 0.088
Ctr/Dup: p = 0.628
Del/Dup: p = 0.003
H(2): 1.55, p = 0.460
H(2): 6.81, p = 0.033
Ctr/Del: p = 1.000
Ctr/Dup: p = 0.226
Del/Dup: p = 0.028
H(2): 10.48, p = 0.005
Ctr/Del: p = 0.207
Ctr/Dup: p = 0.455
Del/Dup: p = 0.005
Visual evoked potential analysis
We analyzed activity over the midline occipital electrode Oz. The N1, P1, and N2 components were identified for each participant’s average waveform, which are equivalent to the N75, P100, and N135  or the C1, P1, and N1 , respectively. The P1 was identified visually as the first prominent positive inflection point closest to 100 ms and occurred within 56–132 ms post-stimulus onset for all participants. The N1 was the negative inflection point immediately preceding the P1 and occurred within 0–70 ms post-stimulus onset for all participants. The N2 was the first negative inflection point following the P1 (or the last peak in a series of multiple peaks, as described previously ) and occurred within 108–266 ms post-stimulus onset in all participants. The absolute amplitude and latency of each component, as well as the relative amplitudes and times between components, were quantified for each individual’s average waveform.
Statistical analysis was performed with IBM SPSS Statistics v.21 software. All data were analyzed blind to group and based on biological replicates (each data point is from one participant). Nonparametric tests were used due to small sample sizes and non-normal distribution of data. Results were considered significant if p < 0.05. The three participant groups (control, deletion, and duplication) were compared using a Kruskal-Wallis test, reported as H(df) where H is the test statistic and df is the degrees of freedom. If significant, the Kruskal-Wallis test was followed by paired comparisons with adjusted p values. Spearman correlations were used to compare the relationships between continuous variables.
When comparing the three groups, no differences in latencies of any component were apparent, with the exception of a difference in N1 latency that did not remain significant for any paired comparisons (Table 2). Instead, there was a consistent effect on the amplitude of the principal positive component (P1), including absolute P1 amplitude and relative N1-P1 and P1-N2 amplitudes. This effect was due to a tendency towards large amplitudes in deletion carriers and small amplitudes in duplication carriers (Table 2, Fig. 1).
Correlations with age and IQ
Correlation with age
Correlation with IQ
Head size is known to be differentially affected by 16p11.2 copy number variation  and this could impact VEP latency and/or amplitude [17, 18]. Head circumference was measured at the beginning of the session in order to select an appropriate size net for each participant. While there was no statistically significant difference in head size between groups (H(2):4.06, p = 0.131), children carrying 16p duplications did tend to have smaller heads (TD: median = 52.0, mean = 52.4, range = 48.5–55.5; deletion: median = 51.5, mean = 52.4, range = 48.0–58.0; duplication: median = 49.8, mean = 50.1, range = 46.0-53.0; unit = cm). There was also no correlation within any group between head circumference and any VEP measurement, including N1-P1 amplitude (TD: r = 0.47, p = 0.127; deletion: r = 0.17, p = 0.494; duplication: r = 0.26; p = 0.108). Taken together, differences in head size are unlikely to explain the group differences in N1-P1 amplitude.
16p11.2 CNV syndrome increases the risk for autism spectrum disorder (ASD), particularly for those carrying deletions [19, 20]. Only a small number of individuals in our dataset had a confirmed diagnosis of ASD (n = 5 deletion carriers, n = 1 duplication carrier), and we designated these individuals as open circles in the scatter plots in Fig. 1c, d in order to visualize any trends. All individuals with ASD had P1 latencies close to the median value for their respective group (Fig. 1c). Individuals with ASD tended to have P1 amplitudes above the median value for their respective group. This result indicates that individuals carrying 16p deletions and a diagnosis of ASD may have larger VEP amplitudes, but future analysis of larger datasets is necessary.
The differing effect of 16p11.2 deletion and duplication on VEP P1 amplitude is intriguing and may offer insight into the neural differences between the two genetic conditions. It is known that deletion and duplication have opposite effects on head circumference (macrocephaly and microcephaly, respectively) , brain volume , and white matter microstructure . Little is understood about how brain size and volume may affect VEP amplitude measured from the scalp, but this could be a contributing factor. It must also be considered that reduced cortical thickness has been reported in both deletion and duplication carriers , a finding that may be particularly relevant for a cortical recording method like VEP.
A reduction in P1 amplitude of the VEP has been found in psychoses [23–25], and 16p11.2 CNVs are associated with schizophrenia and bipolar disorder [20, 26]. This potential association could be investigated further in an older cohort of individuals with 16p11.2 CNVs.
A similar increase in VEP P1 amplitude has been found in mice carrying deletion of 16p11.2 (S. Cooke and M. Bear, personal communication). If this holds true, this parallel result could validate the translatable relevance of the mouse model to the human condition and will ultimately help us to better understand the synaptic and circuit mechanisms underlying cortical dysfunction in 16p CNV syndrome. Some recent in-depth analysis of cortical organization has been done in mice carrying 16p deletions. Specifically, these mice display fewer upper layer pyramidal cells and increased layer 6 corticothalamic projection neurons due to ERK dysregulation . Such changes in cortical microcircuitry could contribute to the changes in VEP amplitude observed in both mice and humans.
Interestingly, auditory responses measured using magnoencephalography were reported to be delayed in children with 16p deletions but not duplications, and there was no contribution of age or IQ to this result . The representative traces displayed in the first figure of that paper appear to show increased signal amplitude in deletions and decreased amplitude in duplications; however, quantification of amplitude was not reported so this may not be a significant finding. Here, we used EEG to assess visual processing in our study and found differences in amplitude and not latency. It is unclear if our findings differ due to the sensory system that was being assessed or the method used.
Limitations of our study
One potential limitation of our study is the broad age range of our sample that spans childhood. However, the VEP is mature by 2 years of age [7–11], which is well below the age of our youngest participant, and we saw no effect of age on any VEP component within any group. Another limitation is the small sample size and variability of phenotype and comorbid conditions, including seizure history, medication use, and ASD diagnosis. This is a challenge inherent in studying this population and means that any differences that emerge from group comparison are then potentially specific to the genetic condition. Cognitive ability is one phenotypic variable that we were not able to adequately address in this study because IQ was not assessed in the typically developing group due to the nature of the study design. IQ did not differ between the 16p deletion and duplication groups, however, and because there was no correlation between IQ and VEP amplitude in the deletion group, IQ is not likely to explain the differences in P1 amplitude in our dataset.
Finally, because visual information must first travel through the retina and thalamus before reaching the cortex where the VEP is recorded, we cannot exclude that ophthalmological abnormalities could contribute to any group differences we see. However, no participant was reported to have had eye surgery, cataract, coloboma, glaucoma, or any other uncorrected visual problems. In addition, the latency of the components was unaffected, suggesting that the arrival of visual information from the periphery to the cortex was intact.
Our study has shown that it is possible to successfully record VEP in children with 16p11.2 CNV disorder who present with a range of intellectual, behavioral, and physical disabilities. We have shown that VEP responses are robust, easily quantifiable, and unaffected by age, head circumference, or IQ. Furthermore, the amplitude of the principal P1 component is enhanced in deletion carriers and diminished in duplication carriers, providing a biomarker to further probe cortical processing in the two conditions.
ASD, autism spectrum disorder; BCH, Boston Children’s Hospital; CNV, copy number variation; EEG, electroencephalogram; IQ, intelligence quotient; SFARI, Simons Foundation Autism Research Initiative; SVIP, Simons Variation in Individuals Project; VEP, visual evoked potential
This work was supported by a grant from the Simons Foundation (SFARI #206787, C.A.N). We would like to thank the research staff members who coordinated the study and collected the data, including V. Vogel-Farley, M. Crossman, T. Clarkson, G. Holt, and M-K. Driscoll. We would also like to thank the participating children and families who generously contributed their time.
JJL performed the data analysis, interpreted the results, and wrote the manuscript. CAN designed the study, interpreted the results, and edited the manuscript. All authors read and approved the manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
All procedures were reviewed and approved by the Boston Children’s Hospital (BCH) Office of Clinical Investigation and the BCH Institutional Review Board, and written informed consent was obtained from each participant or guardian prior to testing.
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.
- Shinawi M, Liu P, Kang SL, Shen J, Belmont JW, Scott DA, Probst FJ, Craigen WJ, Graham BH, Clark G, Lee J, Proud M, Stocco A, Rodriguez DL, Kozel BA, Sparagana S, Roeder ER, Mcgrew SG, Thaddeus W. Recurrent reciprocal 16p11.2 rearrangements associated with global developmental delay, behavioural problems, dysmorphism, epilepsy, and abnormal head size. J Med Genet. 2010;47:332–41.View ArticlePubMedGoogle Scholar
- Jacquemont S, Reymond A, Zufferey F, Harewood L, Walters RG, Kutalik Z, Martinet D, Shen Y, Valsesia A, Beckmann ND, Thorleifsson G, Belfiore M, Bouquillon S, Campion D, de Leeuw N, de Vries BBA, Esko T, Fernandez BA, Fernández-Aranda F, Fernández-Real JM, Gratacòs M, Guilmatre A, Hoyer J, Jarvelin M-R, Frank Kooy R, Kurg A, Le Caignec C, Männik K, Platt OS, Sanlaville D, et al. Mirror extreme BMI phenotypes associated with gene dosage at the chromosome 16p11.2 locus. Nature. 2011;478:97–102.View ArticlePubMedPubMed CentralGoogle Scholar
- Zufferey F, Sherr EH, Beckmann ND, Hanson E, Maillard AM, Hippolyte L, Mace A, Ferrari C, Kutalik Z, Andrieux J, Aylward E, Barker M, Bernier R, Bouquillon S, Conus P, Delobel B, Faucett WA, Goin-Kochel RP, Grant E, Harewood L, Hunter JV, Lebon S, Ledbetter DH, Martin CL, Mannik K, Martinet D, Mukherjee P, Ramocki MB, Spence SJ, Steinman KJ, et al. A 600 kb deletion syndrome at 16p11.2 leads to energy imbalance and neuropsychiatric disorders. J Med Genet. 2012;49:660–8.View ArticlePubMedPubMed CentralGoogle Scholar
- LeBlanc JJ, DeGregorio G, Centofante E, Vogel-Farley VK, Barnes K, Kaufmann WE, Fagiolini M, Nelson CA. Visual evoked potentials detect cortical processing deficits in Rett syndrome. Ann Neurol. 2015;78:775–86.View ArticlePubMedGoogle Scholar
- Jenkins J, Chow V, Blaskey L, Kuschner E, Qasmieh S, Gaetz L, Edgar JC, Mukherjee P, Buckner R, Nagarajan SS, Chung WK, Spiro JE, Sherr EH, Berman JI, Roberts TPL. Auditory evoked M100 response latency is delayed in children with 16p11.2 deletion but not 16p11.2 duplication. Cereb Cortex. 2015;26:1–8.Google Scholar
- Varcin KJ, Nelson CA, Ko J, Sahin M, Wu JY, Jeste SS. Visual evoked potentials as a readout of cortical function in infants with tuberous sclerosis complex. J Child Neurol. 2015;31(2):195–202.View ArticlePubMedGoogle Scholar
- Sokol S. Measurement of infant visual acuity from pattern reversal evoked potentials. Vision Res. 1978;18:33–9.View ArticlePubMedGoogle Scholar
- Sokol S, Jones K. Implicit time of pattern evoked potentials in infants: an index of maturation of spatial vision. Vision Res. 1979;19:747–55.View ArticlePubMedGoogle Scholar
- McCulloch DL, Skarf B. Development of the human visual system: monocular and binocular pattern VEP latency. Investig Ophthalmol Vis Sci. 1991;32:2372–81.Google Scholar
- Lippe S, Roy M-S, Perchet C, Lassonde M. Electrophysiological markers of visuocortical development. Cereb Cortex. 2007;17(January):100–7.PubMedGoogle Scholar
- Lenassi E, Likar K, Stirn-Kranjc B, Brecelj J. VEP maturation and visual acuity in infants and preschool children. Doc Ophthalmol. 2008;117:111–20.View ArticlePubMedGoogle Scholar
- Odom JV, Bach M, Brigell M, Holder GE, McCulloch DL, Tormene AP, Vaegan. SCEV standard for clinical visual evoked potentials (2009 update). Doc Ophthalmol. 2009;2010(120):111–9.Google Scholar
- Qureshi AY, Mueller S, Snyder AZ, Mukherjee P, Berman JI, Roberts TPL, Nagarajan SS, Spiro JE, Chung WK, Sherr EH, Buckner RL. Opposing brain differences in 16p11.2 deletion and duplication carriers. J Neurosci. 2014;34:11199–211.View ArticlePubMedPubMed CentralGoogle Scholar
- Pucilowska J, Vithayathil J, Tavares EJ, Kelly C, Karlo JC, Landreth GE. The 16p11.2 deletion mouse model of autism exhibits altered cortical progenitor proliferation and brain cytoarchitecture linked to the ERK MAPK pathway. J Neurosci. 2015;35:3190–200.View ArticlePubMedGoogle Scholar
- Di Russo F, Martínez A, Sereno MI, Pitzalis S, Hillyard SA. Cortical sources of the early components of the visual evoked potential. Hum Brain Mapp. 2002;15:95–111.View ArticlePubMedGoogle Scholar
- Jones SJ, Turano G, Kriss A, Shawkat F, Kendall B, Thompson AJ. Visual evoked potentials in phenylketonuria: association with brain MRI, dietary state, and IQ. J Neurol Neurosurg Psychiatry. 1995;59:260–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Erwin C, Erwin A, Hartwell J, Wilson W. P100 latency as a function of head size. Am J EEG Technol. 1991;31:279–88.Google Scholar
- Guthkelch AN, Bursick D, Sclabassi RJ. The relationship of the latency of the visual P100 wave to gender and head size. Electroencephalogr Clin Neurophysiol Evoked Potentials. 1987;68:219–22.View ArticlePubMedGoogle Scholar
- Kumar RA, Karamohamed S, Sudi J, Conrad DF, Brune C, Badner JA, Gilliam TC, Nowak NJ, Cook EH, Dobyns WB, Christian SL. Recurrent 16p11.2 microdeletions in autism. Hum Mol Genet. 2008;17:628–38.View ArticlePubMedGoogle Scholar
- McCarthy SE, Makarov V, Kirov G, Addington AM, McClellan J, Yoon S, Perkins DO, Dickel DE, Kusenda M, Krastoshevsky O, Krause V, Kumar RA, Grozeva D, Malhotra D, Walsh T, Zackai EH, Kaplan P, Ganesh J, Krantz ID, Spinner NB, Roccanova P, Bhandari A, Pavon K, Lakshmi B, Leotta A, Kendall J, Lee Y, Vacic V, Gary S, Iakoucheva LM, et al. Microduplications of 16p11.2 are associated with schizophrenia. Nat Genet. 2009;41:1223–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Owen JP, Chang YS, Pojman NJ, Bukshpun P, Wakahiro MLJ, Marco EJ, Berman JI, Spiro JE, Chung WK, Buckner RL, Roberts TPL, Nagarajan SS, Sherr EH, Mukherjee P. Aberrant white matter microstructure in children with 16p11.2 deletions. J Neurosci. 2014;34:6214–23.View ArticlePubMedGoogle Scholar
- Maillard AM, Ruef A, Pizzagalli F, Migliavacca E, Hippolyte L, Adaszewski S, Dukart J, Ferrari C, Conus P, Männik K, Zazhytska M, Siffredi V, Maeder P, Kutalik Z, Kherif F, Hadjikhani N, Beckmann JS, Reymond A, Draganski B, Jacquemont S. The 16p11.2 locus modulates brain structures common to autism, schizophrenia and obesity. Mol Psychiatry. 2015;20:140–7.View ArticlePubMedGoogle Scholar
- Yeap S, MRCPsych, Kelly S, Sehatpour P, Magno E, Javitt DC, Garavan H, Thakore JH, Foxe JJ. Early visual sensory deficits as endophenotypes for schizophrenia. Arch Gen Psychiatry. 2006;63(Nov):1180–8.View ArticlePubMedGoogle Scholar
- Yeap S, Kelly SP, Sehatpour P, Magno E, Garavan H, Thakore JH, Foxe JJ. Visual sensory processing deficits in Schizophrenia and their relationship to disease state. Eur Arch Psychiatry Clin Neurosci. 2008;258:305–16.View ArticlePubMedGoogle Scholar
- Lalor EC, Yeap S, Reilly RB, Pearlmutter BA, Foxe JJ. Dissecting the cellular contributions to early visual sensory processing deficits in schizophrenia using the VESPA evoked response. Schizophr Res. 2008;98:256–64.View ArticlePubMedGoogle Scholar
- Green EK, Rees E, Walters JTR, Smith K, Forty L, Grozeva D, Moran JL, Sklar P, Ripke S, Chambert KD, Genovese G. Copy number variation in bipolar disorder. Mol Psychiatry. 2016;21:89–93.View ArticlePubMedGoogle Scholar