Molecular Vision 2012; 18:528-536 <http://www.molvis.org/molvis/v18/a58>
Received 23 November 2011 | Accepted 27 February 2012 | Published 1 March 2012

A Chinese family with Oguchi’s disease due to compound heterozygosity including a novel deletion in the arrestin gene

Lingli Huang,1 Wen Li,1,2 Weilin Tang,2 Xiaohua Zhu,3 Pingbo Ou-yang,3 Guangxiu Lu1,2

1Institute of Reproductive and Stem Cell Engineering, Central South University, Changsha, P.R. China; 2Reproductive and Genetic Hospital of Citic-Xiangya, Changsha, P.R. China; 3Department of Ophthalmology, the Second Affiliated Xiangya Hospital, Central South University, Changsha, P.R. China

Correspondence to: Guangxiu Lu, Institute of Reproductive and Stem Cell Engineering, Central South University, Changsha 410078, P.R. China; Phone: +86-731-84805319; FAX: +86-731-84497661; email: lugxdirector@yahoo.com.cn.

Abstract

Purpose: Oguchi’s disease is a rare autosomal recessive disease and known to be caused by mutations in the rhodopsin kinase (GRK1) gene or the arrestin (SAG) gene. SAG contains 16 exons and encodes a protein with 405 amino acids. This study was to identify the underlying genetic defects in a non-consanguineous Chinese family with Oguchi’s disease.

Methods: Ophthalmologic examinations including fundus photography and electroretinography (ERG) were performed on all family members. All exons of the GRK1 gene and the SAG gene were amplified with PCR and directly sequenced. Quantitative real-time PCR (qPCR) was performed to screen heterozygous deletions/duplications in the SAG gene. Long-range PCR and direct sequencing were further performed to define the breakpoints.

Results: The patient had characteristic clinical features of Oguchi’s disease, including night blindness, normal vision fields, typical fundus appearance with the Mizuo-Nakamura phenomenon, nearly undetectable rod b waves in the scotopic 0.01 ERGs, and nearly “negative” scotopic 3.0 ERGs. No mutations were found in the GRK1 gene. A heterozygous nonsense Arg193stop (R193X) mutation was found in the SAG gene in the patient and the unaffected mother. No pathogenic SAG mutations were found in the unaffected father. qPCRs showed a heterozygous deletion encompassing exon 2 of the SAG gene in the patient and the unaffected father. Long-range PCR and direct sequencing verified the deletion and revealed the breakpoints of the deletion, skipping a 3,224-bp fragment of the SAG gene. The deletion was not detected in 96 unrelated healthy controls. This deletion was predicted to eliminate the exon 2 and the AUG initiate codon in the mature SAG mRNA and cause no production of the SAG protein or low-level production of a non-functional truncated protein lacking 134 amino acids in the NH2 terminus.

Conclusions: Compound heterozygosity of a nonsense R193X mutation and a heterozygous deletion of 3,224 bp encompassing exon 2 in the SAG gene is the cause of Oguchi’s disease in this Chinese family. qPCR analysis should be performed if there is a negative result of the mutation screening of the SAG gene in patients with Oguchi’s disease.

Introduction

Oguchi’s disease (OMIM 258100) is a rare form of stationary night blindness with autosomal recessive inheritance, characterized by a typical clinical feature called the Mizuo-Nakamura phenomenon in which the golden-yellow discoloration of the fundus disappears in the dark-adapted condition and reappears shortly after exposure to light [1]. In addition, patients with Oguchi’s disease usually have night blindness but normal color vision, visual acuity, and visual field [2]. Electroretinographic examinations showed reduced or no rod functions with normal cone functions [2,3].

Two causative genes have been reported for Oguchi’s disease: the arrestin (SAG; OMIM 181031) gene [1] and the rhodopsin kinase (GRK1; OMIM 180381) gene [4]. Beginning in 1995 with the first report on a homozygous SAG frameshift 926delA (formerly referred to as 1147delA) [5-7] mutation in five of six unrelated Japanese patients [1], more SAG mutations have been found in Oguchi’s disease, including an additional homozygous frameshift 926delA mutation in seven other Japanese families [8-12], a homozygous nonsense Arg193stop (R193X) mutation in an Indian family [13], a compound heterozygous mutation of a nonsense R175X mutation plus a frameshift 926delA mutation in a Japanese family [14], and a homozygous nonsense R292X mutation in a Japanese family [14] and in a South Asian family [15]. Recently, an additional heterozygous frameshift 926delA mutation without identification of the other mutation in another allele during mutation screening of the SAG gene was reported in a Japanese patient [16], suggesting that other types of mutations, undetectable with mutation screening of the SAG gene, may exist.

In this study, we report a heterozygous nonsense R193X mutation and a novel heterozygous deletion of 3,224 bp encompassing exon 2 in the SAG gene in a Chinese family. The deletion was identified using quantitative real-time PCR (qPCR). To our knowledge, this is the first case found in a Chinese family and is the first report of a novel heterozygous deletion in the SAG gene.

Methods

A Chinese family with a 13-year-old female patient and her unaffected parents participated in this study. This study conformed to the tenets of the Declaration of Helsinki, and the research protocol was approved by the Ethics Committee of the Reproductive and Genetic Hospital of Citic-Xiangya (Changsha, P.R. China). Informed consent was obtained from all family members after an explanation of the purpose of this study was provided.

Clinical examinations

Family members were clinically examined at the Department of Ophthalmology, the Second Affiliated Xiangya Hospital of Central South University. Ophthalmologic examinations included best-corrected visual acuity (BCVA), slit lamp biomicroscopy, fundus photography, kinetic perimetry, and electroretinography. Kinetic perimetry was performed on a Twinfield perimeter (Oculus Inc., Wetzlar, Germany) using defined stimuli according to the Goldmann standard. Standard full-field electroretinograms (ERGs) were elicited with Ganzfeld stimuli after 30 min of dark adaptation using the commercial ERG system (RetiPort32; Roland Consult Systems, Bradenburg, Germany) according to the guidelines of the International Society for Clinical Electrophysiology of Vision (ISCEV)(standard flash, 3 cd·s·m−2) [17]. The stimuli were 0.0095 cd·s·m−2 for rod stimulation (Scotopic 0.01 ERG) and 3 cd·s·m−2 for other standard responses (including scotopic 3.0 ERG, photopic 3.0 ERG, and photopic 3.0 flicker). For photopic ERG, the background luminance was set at 34 cd·m−2.

Mutation screening

All exons of the GRK1 gene, including exon/intron boundaries, were amplified with PCR and directly sequenced as previously described [18]. All exons of the SAG gene with exon/intron boundaries were amplified with PCR using primers (Table 1) designed by Primer 3 and checked by NCBI BLAST for specificity followed by direct sequencing for mutation screening.

Quantitative real-time polymerase chain reaction (qPCR)

qPCR was performed to detect exonic deletions/duplications and fine mapping of the breakpoints. Primers were designed with an amplicon size of less than 300 bp (Table 2). For fine mapping of the heterozygous deletion encompassing exon 2 identified in qPCR, additional primer pairs were designed and named In1 (located in intron 1), In2–1, In2–2, In2–3, In2–4, and In2–5 (located in intron 2; Table 2). The albumin (ALB) gene was used as the reference gene [19]. The qPCR reactions were performed in 10 µl volume containing 1×iQTM SYBR® Green Supermix (Bio-Rad Laboratories, Inc., Hercules, CA), primers at a final concentration of 0.25 μM each and genomic DNA at a concentration of 5 ng/μl. The qPCR reactions were run in the 7500 Real-time PCR system (Applied Biosystems, Foster City, CA). Cycling conditions were 95 °C for 10 min, and 40 cycles of 95 °C for 15 s and 60 °C for 1 min. A melting curve analysis of products was performed routinely following the amplification to test the specificity of the PCR products. Genomic DNA from a normal individual was used as a deletion-negative control (the calibrator sample) to normalization. Two independent experiments were performed. Quantitative PCR data were analyzed with the comparative threshold cycle method [20,21]. A relative quantity of about 1 in a normal sample, about 0.5 in a sample with a heterozygous SAG exon deletion, and about 1.5 in a sample with a heterozygous SAG exon duplication was expected.

Long-range polymerase chain reaction

To validate the heterozygous deletion detected with qPCR analysis, a long-range PCR approach was performed using the primers In1-F and In2–4-R (Table 2). The PCR conditions were 95 °C for 5 min, and 35 cycles of 94 °C for 30 s, 60 °C for 30 s, and 68 °C for 1 min, followed by the final extension at 68 °C for 5 min. The purified PCR products were directly sequenced.

Results

Clinical findings

The proband was a 13-year-old girl who had night blindness since early childhood. Her BCVA was 1.5 in each eye. The refractive errors were −3.50 −0.50×180° in the right eye and −3.75 −0.50×180° in the left eye. Fundus examinations showed the characteristic golden-yellow discoloration in both eyes and the Mizuo-Nakamura phenomenon, in which the golden-yellow discoloration of the fundus (Figure 1A) disappeared after prolonged dark adaptation (Figure 1B). Neither vascular attenuation nor retinal degeneration was seen throughout the retina, and no maculopathy was observed. The visual fields of both eyes were within normal limits. The full-field scotopic 0.01 ERGs showed undetectable rod b waves in the patient. The scotopic 3.0 ERGs showed a “negative” configuration with a significantly reduced a wave and a nearly absent b wave in both eyes. The photopic 3.0 ERGs and photopic 3.0 flickers seemed almost normal (Figure 2). The girl was the first offspring of an unaffected non-consanguineous couple, and no family history was found in the family pedigree (Figure 3A).

Molecular genetic findings

No mutations in the GRK1 gene were found in the patient. Mutation screening of the SAG gene revealed a heterozygous nonsense mutation R193X (c.577C>T; Figure 3B,C), and three heterozygous single nucleotide polymorphisms (SNPs, recorded in the SNP database [build 130]) including IVS6–18G>C, IVS9+31T>G, and c.1207G>A (V430I; Figure 3B) in the patient. Among four nucleotide changes, the heterozygous nucleotide change IVS9+31T>G was seen in both parents while the other three changes were found only in the mother, suggesting those changes were inherited from the mother (Figure 3B). No putative pathogenic mutation was identified in the father.

The existence of four heterozygous nucleotide changes in the patient suggested that exon 7, exon 8, exon 9, and exon 16 were deletion-negative. qPCR analysis of the other 12 exons of the SAG gene was performed and revealed a heterozygous deletion of the SAG exon 2 in the patient and her unaffected father (Figure 4A). Fine mapping of the exon 2 deletion with qPCR analysis using the primers including In1, In2–1, In2–2, In2–3, In2–4, and In2–5 revealed the heterozygous deletion involving In2–1, In2–2, and In2–3 in the patient and her father (Figure 4B, Figure 5A).

Long-range PCR using the primers In1-F and In2–4-R amplified a shortened fragment of about 1,500 bp in the patient and her father (Figure 5B). Direct sequencing of the long-range PCR products revealed the breakpoints of the deletion, skipping a 3,224-bp fragment of the SAG gene. The 5′-end of the breakpoint was located in the intron 1/exon 2 boundary, and the 3′-end of the breakpoint was located in the intron 2, 3,122 bp downstream of exon 2 (IVS2+3122; Figure 5C). No shortened fragment was seen in 96 unrelated healthy controls.

Discussion

Oguchi’s disease is a rare autosomal recessive inherited disease. Patients usually have typical clinical features including night blindness, typical fundus changes called the Mizuo-Nakamura phenomenon, and characteristic ERGs revealing reduced or no rod function with normal cone function [1-3]. The patient described here could be clinically diagnosed with Oguchi’s disease according to the typical ophthalmologic features revealed by clinical examinations.

Two causative genes, SAG [1] and GRK1 [4], have been reported for Oguchi’s disease. To date, only four SAG mutations including 926delA, R193X, R175X, and R292X have been reported in 16 Japanese families [1,8-12,14,16], one Indian family [13], and one South Asian family [15] with Oguchi’s disease. The most recent report of Oguchi’s disease was identified in a Japanese patient with only one allele of the frameshift 926delA mutation [16], suggesting that other types of mutation, undetectable by mutation screening of the coding regions of the SAG gene, may exist in the other allele of the SAG gene in the patient.

In our study, only one heterozygous nonsense R193X mutation was identified in the patient and her mother by mutation screening of the coding regions of the SAG gene. The homozygous nonsense R193X mutation has been reported in only one Indian family with Oguchi’s disease. qPCR analysis followed by long-range PCR revealed a novel heterozygous deletion of 3,224 bp, encompassing exon 2 and partial intron 2, in the patient and her father. To our knowledge, this is the first report of a novel heterozygous deletion in the SAG gene identified using qPCR and the first case of Oguchi’s disease found in a Chinese family.

This deletion was predicted to eliminate the entire exon 2 of 103 nucleotides (nt) with the AUG initiation site in the mature mRNA of SAG and cause the use of downstream alternative AUG codons (Figure 6). The eukaryotic ribosomes are loaded on the 5′-cap of the mutant mRNA, scan for the translation initiation signal (TIS), and initiate the first AUG codon they encounter [22]. The first downstream AUG codon out of frame located in the exon 3 is recognized by the ribosome, initiating the translation of its open reading frame (ORF). A short peptide of 19 amino acids is produced, which is expected to stall ribosomes to its stop codon [23,24]. In addition, the premature termination codon leads to the rapid mRNA degradation by the pathway of the nonsense-mediated mRNA decay (NMD) [25]. If possible, the mRNAs if not all degraded may remain connected to the un-disassociated ribosomes stalled at the stop codon of the ORF, and thus, these ribosomes may resume scanning and reinitiate the next downstream TIS [26,27]. In this case, there would be three additional downstream out-of-frame AUG codons to be scanned by the ribosomes in the same way. And if the fifth in-frame AUG codon has a chance to be initiated, a predicted protein lacking 134 amino acids in the NH2 terminus would be expressed in a dramatically reduced level. Even if the expression level is high enough, the predicted truncated protein could be non-functional with impaired conformation of the NH2-terminal half of the SAG protein and disrupted interaction between the N domain and the C domain [28-31].

Interestingly, no missense mutations have been reported so far. All the previously reported mutations, nonsense or frameshift, were predicted to produce truncated proteins lacking the COOH terminus if the transcripts of the gene were not all degraded by NMD [1,8-16]. Even if the expression level is high enough, the predicted truncated proteins could be non-functional [1,8-16]. Taken together with the novel deletion reported in this study, these results led us to suggest that the predicted aberrant SAG gene products with a loss of either the COOH-terminus or the NH2-terminus caused by nonsense mutations, frameshift mutations, or large fragment deletions eliminating natural AUG codon, having great potential to impair or abolish the function of the SAG protein and thus cause its transcripts to be exposed to the NMD [25], will cause Oguchi’s disease. Missense mutations that are not the target of NMD may remain some basic function of the SAG protein and therefore cause mild Oguchi’s disease without significant night blindness. However, more reports of Oguchi’s disease may help to elucidate the underlying mechanism.

Recently, two recorded copy number variants in the Database of Genomic Variants (DGV), Gain_79017 [32] and Loss_53177 [33], have been identified in normal individuals with incidence of 2 in 90 (>2.2%) and 2 in 1,854 (1,064 plus 790; approximately 1‰), respectively. Those two variants could be potential hot-spot rearrangement sites for Oguchi’s disease. The duplication of the exon 3 (Gain_79017) is predicted to be pathogenic by inserting 61 nucleotides in the coding region of the mRNA transcripts, causing a frameshift effect and triggering the NMD. The loss of the exon 11 to 16 (Loss_53177) is also predicted to be pathogenic by lacking the function of the C-terminus. Therefore, the incidence of Oguchi’s disease caused by the SAG gene mutations could be higher than expected. The relative low number of reports on SAG gene mutations could partly be explained by the putative pathogenic SAG deletions/duplications that escape the most common mutation screening methods for detecting point or minute mutations. The identification of the first deletion in the SAG gene described in this study will suggest the idea that the search for possible rearrangements in the SAG gene should be routinely performed in all patients with Oguchi’s disease if there is a negative result of the mutation screening of the SAG gene in one or two alleles.

In summary, the compound heterozygosity for a heterozygous nonsense mutation and a heterozygous deletion of 3,224 bp encompassing exon 2 in the SAG gene would putatively explain Oguchi’s disease in a Chinese family. Here, we emphasize that qPCR analysis should be performed if the mutation screening of the SAG gene fails to detect the aberrant variants in patients with Oguchi’s disease. The novel heterozygous deletion in the SAG gene identified with the qPCR method in this study will expand the spectrum of SAG mutations associated with Oguchi’s disease, and this will help to elucidate further the role of this gene in the etiology of Oguchi’s disease.

References

  1. Fuchs S, Nakazawa M, Maw M, Tamai M, Oguchi Y, Gal A. A homozygous 1-base pair deletion in the arrestin gene is a frequent cause of Oguchi disease in Japanese. Nat Genet. 1995; 10:360-2. [PMID: 7670478]
  2. Carr RE, Gouras P. Oguchi's Disease. Arch Ophthalmol. 1965; 73:646-56. [PMID: 14281981]
  3. Miyake Y, Horiguchi M, Suzuki S, Kondo M, Tanikawa A. Electrophysiological findings in patients with Oguchi's disease. Jpn J Ophthalmol. 1996; 40:511-9. [PMID: 9130055]
  4. Yamamoto S, Sippel KC, Berson EL, Dryja TP. Defects in the rhodopsin kinase gene in the Oguchi form of stationary night blindness. Nat Genet. 1997; 15:175-8. [PMID: 9020843]
  5. den Dunnen JT, Antonarakis SE. Mutation nomenclature extensions and suggestions to describe complex mutations: a discussion. Hum Mutat. 2000; 15:7-12. [PMID: 10612815]
  6. Saga M, Mashima Y, Kudoh J, Oguchi Y, Shimizu N. Gene analysis and evaluation of the single founder effect in Japanese patients with Oguchi disease. Jpn J Ophthalmol. 2004; 48:350-2. [PMID: 15295660]
  7. Yoshida S, Yamaji Y, Yoshida A, Ikeda Y, Yamamoto K, Ishibashi T. Rapid detection of SAG 926delA mutation using real-time polymerase chain reaction. Mol Vis. 2006; 12:1552-7. [PMID: 17200654]
  8. Nakazawa M, Wada Y, Fuchs S, Gal A, Tamai M. Oguchi disease: phenotypic characteristics of patients with the frequent 1147delA mutation in the arrestin gene. Retina. 1997; 17:17-22. [PMID: 9051837]
  9. Nakamachi Y, Nakamura M, Fujii S, Yamamoto M, Okubo K. Oguchi disease with sectoral retinitis pigmentosa harboring adenine deletion at position 1147 in the arrestin gene. Am J Ophthalmol. 1998; 125:249-51. [PMID: 9467455]
  10. Yoshii M, Murakami A, Akeo K, Nakamura A, Shimoyama M, Ikeda Y, Kikuchi Y, Okisaka S, Yanashima K, Oguchi Y. Visual function and gene analysis in a family with Oguchi's disease. Ophthalmic Res. 1998; 30:394-401. [PMID: 9731122]
  11. Yamada T, Matsumoto M, Kadoi C, Nagaki Y, Hayasaka Y, Hayasaka S. 1147 del A mutation in the arrestin gene in Japanese patients with Oguchi disease. Ophthalmic Genet. 1999; 20:117-20. [PMID: 10420197]
  12. Hayashi T, Tsuzuranuki S, Kozaki K, Urashima M, Tsuneoka H. Macular Dysfunction in Oguchi Disease with the Frequent Mutation 1147delA in the SAG Gene. Ophthalmic Res. 2011; 46:175-80. [PMID: 21447990]
  13. Maw M, Kumaramanickavel G, Kar B, John S, Bridges R, Denton M. Two Indian siblings with Oguchi disease are homozygous for an arrestin mutation encoding premature termination. Hum Mutat. 1998; Suppl 1:S317-9. [PMID: 9452120]
  14. Nakamura M, Yamamoto S, Okada M, Ito S, Tano Y, Miyake Y. Novel mutations in the arrestin gene and associated clinical features in Japanese patients with Oguchi's disease. Ophthalmology. 2004; 111:1410-4. [PMID: 15234147]
  15. Sergouniotis PI, Davidson AE, Sehmi K, Webster AR, Robson AG, Moore AT. Mizuo-Nakamura phenomenon in Oguchi disease due to a homozygous nonsense mutation in the SAG gene. Eye (Lond). 2011; 25:1098-101. [PMID: 21494281]
  16. Fujinami K, Tsunoda K, Nakamura M, Oguchi Y, Miyake Y. Oguchi Disease With Unusual Findings Associated With a Heterozygous Mutation in the SAG Gene. Arch Ophthalmol. 2011; 129:1375-6. [PMID: 21987685]
  17. Marmor MF, Fulton AB, Holder GE, Miyake Y, Brigell M, Bach M. ISCEV Standard for full-field clinical electroretinography (2008 update). Doc Ophthalmol. 2009; 118:69-77. [PMID: 19030905]
  18. Zhang Q, Zulfiqar F, Riazuddin SA, Xiao X, Yasmeen A, Rogan PK, Caruso R, Sieving PA, Riazuddin S, Hejtmancik JF. A variant form of Oguchi disease mapped to 13q34 associated with partial deletion of GRK1 gene. Mol Vis. 2005; 11:977-85. [PMID: 16319817]
  19. Laurendeau I, Bahuau M, Vodovar N, Larramendy C, Olivi M, Bieche I, Vidaud M, Vidaud D. TaqMan PCR-based gene dosage assay for predictive testing in individuals from a cancer family with INK4 locus haploinsufficiency. Clin Chem. 1999; 45:982-6. [PMID: 10388473]
  20. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Methods. 2001; 25:402-8. [PMID: 11846609]
  21. Hoebeeck J, Speleman F, Vandesompele J. Real-time quantitative PCR as an alternative to Southern blot or fluorescence in situ hybridization for detection of gene copy number changes. Methods Mol Biol. 2007; 353:205-26. [PMID: 17332643]
  22. Kozak M. Structural features in eukaryotic mRNAs that modulate the initiation of translation. J Biol Chem. 1991; 266:19867-70. [PMID: 1939050]
  23. Crowe ML, Wang XQ, Rothnagel JA. Evidence for conservation and selection of upstream open reading frames suggests probable encoding of bioactive peptides. BMC Genomics. 2006; 7:16 [PMID: 16438715]
  24. Oyama M, Itagaki C, Hata H, Suzuki Y, Izumi T, Natsume T, Isobe T, Sugano S. Analysis of small human proteins reveals the translation of upstream open reading frames of mRNAs. Genome Res. 2004; 14:2048-52. [PMID: 15489325]
  25. Mendell JT, Sharifi NA, Meyers JL, Martinez-Murillo F, Dietz HC. Nonsense surveillance regulates expression of diverse classes of mammalian transcripts and mutes genomic noise. Nat Genet. 2004; 36:1073-8. [PMID: 15448691]
  26. Kochetov AV. Alternative translation start sites and hidden coding potential of eukaryotic mRNAs. Bioessays. 2008; 30:683-91. [PMID: 18536038]
  27. Wang XQ, Rothnagel JA. 5′-untranslated regions with multiple upstream AUG codons can support low-level translation via leaky scanning and reinitiation. Nucleic Acids Res. 2004; 32:1382-91. [PMID: 14990743]
  28. Gurevich VV, Benovic JL. Visual arrestin interaction with rhodopsin. Sequential multisite binding ensures strict selectivity toward light-activated phosphorylated rhodopsin. J Biol Chem. 1993; 268:11628-38. [PMID: 8505295]
  29. Gurevich VV, Chen CY, Kim CM, Benovic JL. Visual arrestin binding to rhodopsin. Intramolecular interaction between the basic N-terminus and acidic C-terminus of arrestin may regulate binding selectivity. J Biol Chem. 1994; 269:8721-7. [PMID: 8132602]
  30. Vishnivetskiy SA, Paz CL, Schubert C, Hirsch JA, Sigler PB, Gurevich VV. How does arrestin respond to the phosphorylated state of rhodopsin? J Biol Chem. 1999; 274:11451-4. [PMID: 10206946]
  31. Hirsch JA, Schubert C, Gurevich VV, Sigler PB. The 2.8 A crystal structure of visual arrestin: a model for arrestin's regulation. Cell. 1999; 97:257-69. [PMID: 10219246]
  32. Matsuzaki H, Wang PH, Hu J, Rava R, Fu GK. High resolution discovery and confirmation of copy number variants in 90 Yoruba Nigerians. Genome Biol. 2009; 10:R125 [PMID: 19900272]
  33. Itsara A, Cooper GM, Baker C, Girirajan S, Li J, Absher D, Krauss RM, Myers RM, Ridker PM, Chasman DI, Mefford H, Ying P, Nickerson DA, Eichler EE. Population analysis of large copy number variants and hotspots of human genetic disease. Am J Hum Genet. 2009; 84:148-61. [PMID: 19166990]