Rapid and Inexpensive Screening of Genomic Copy Number Variations Using a Novel Quantitative Fluorescent PCR Method
Detection of human microdeletion and microduplication syndromes poses significant burden on public healthcare systems in developing countries. With genome-wide diagnostic assays frequently inaccessible, targeted low-cost PCR-based approaches are preferred. However, their reproducibility depends on equally efficient amplification using a number of target and control primers. To address this, the recently described technique called Microdeletion/Microduplication Quantitative Fluorescent PCR (MQF-PCR) was shown to reliably detect four human syndromes by quantifying DNA amplification in an internally controlled PCR reaction. Here, we confirm its utility in the detection of eight human microdeletion syndromes, including the more common WAGR, Smith-Magenis, and Potocki-Lupski syndromes with 100% sensitivity and 100% specificity. We present selection, design, and performance evaluation of detection primers using variety of approaches. We conclude that MQF-PCR is an easily adaptable method for detection of human pathological chromosomal aberrations.
Repetitive DNA accounts for more than 50% of human genome  and is one of the primary driving forces of evolution of all eukaryotes. While mostly harmless, stretches of highly similar sequences flanking genes with essential functions create hotspots prone to undergo unequal recombination during meiosis. The resulting pathological copy number variations (CNVs)  include microscopic and submicroscopic deletions and duplications (microdeletions and microduplications) that have been implicated as etiologic factors in many human clinical syndromes. Molecular diagnosis of these aberrations has significantly evolved in resolution from simple metaphase karyotyping to genome-wide molecular screens [3, 4]. However, there are serious limitations to the widespread use of complex diagnostic assays, especially the high cost and rapid obsolescence of the necessary equipment, expensive supplies, and the need for high technical expertise. As a result, pangenomic diagnostic approaches remain scarce in developing countries, where often microdeletion and microduplication syndromes remain undiagnosed. We have described recently the development of a simple, rapid, and inexpensive targeted PCR-based CNV detection method called Microdeletion/Microduplication Quantitative Fluorescent PCR (MQF-PCR) . In our original report, we developed assays for two common microdeletion syndromes, the Williams-Beuren syndrome (7q11.23 microdeletion) and the 22q11.2 microdeletion syndromes and discovered that MQF-PCR could detect both with 100% sensitivity and 100% specificity. Additionally, we demonstrated that the same principle could be reliably used for detection of microduplication syndromes, by using samples from patients with the Lubs (MECP2 duplication) syndrome and the 17q11.2 microduplication involving the NF1 gene .
Now, we present the development of MQF-PCR primer pairs to detect several other microdeletion syndromes: Wilms tumour, aniridia, genitourinary anomalies, and mental retardation syndrome (WAGR; OMIM: 194072), Smith-Magenis (OMIM: 182290), Potocki-Lupski (OMIM: 610883), microdeletions at Xp22.3 (OMIM: 308100, 308700), Langer-Giedion (OMIM: 150230), Greig cephalopolysyndactyly (OMIM: 175700), Miller-Dieker lissencephaly (OMIM: 247200), and Saethre-Chotzen (OMIM: 101400) syndromes.
2. Materials and Methods
2.1. Ethics Statement
DNA samples previously diagnosed with known microdeletion syndromes at GENE-Núcleo de Genética Médica, Belo Horizonte, MG, were anonymized and used in testing the diagnostic procedure. The participants did not provide written or verbal consent and no other formal documented measures were taken. This procedure was approved by the Ethics Committee of Santa Casa de Misericórdia in Belo Horizonte with the registration number CEP 009/2008. WAGR patients were recruited through the International WAGR Syndrome Association. The study was approved by the institutional review board of the National Institute of Child Health and Human Development. Written informed consent was obtained from adults who were competent to provide consent and from the parents or legal guardians of children and adults with cognitive impairment. Assent was obtained from patients who were unable to give consent. Patient samples with Smith-Magenis and Potocki-Lupski syndromes were collected after informed consent was obtained and in accordance with local institutional review board approved protocols from Virginia Commonwealth University and/or Michigan State University.
Genomic DNA from normal male and female individuals and from individuals previously diagnosed with microdeletion syndromes diagnosed at GENE Núcleo de Genética Médica, Belo Horizonte, MG, Brazil, was isolated as described . WAGR sample GM06803 was obtained from Coriell Cell Repositories (Camden, New Jersey, USA). The remaining WAGR samples diagnosed by genome-wide aCGH were described previously . DNA samples of Smith-Magenis and Potocki-Lupski patients diagnosed by FISH or metaphase karyotyping were published previously . In order to simulate Xp22.3 microdeletions, we used normal female and male samples as control individuals and affected samples, respectively. The Saethre-Chotzen syndrome patient 10051 had microdeletion at 7p21.1 diagnosed by aCGH. The Langer-Giedion syndrome patient 12550 showed cytogenetically visible partial deletion of the long arm of chromosome 8 and was diagnosed at our clinic. The remaining samples were obtained from Coriell Cell Repositories: GM09888 (Langer-Giedion), GM10925 (Greig cephalopolysyndactyly), and GM09208 (Miller-Dieker lissencephaly).
2.3. MQF-PCR Primer Selection, Assay Design, and Analysis
MQF-PCR takes advantage of the presence of highly similar sequences in distinct parts of human genome, which differ in length by a few base pairs (Figure 1(a)). A detailed literature search of previously published molecular genetic analyses was performed for each of the syndromes in order to determine minimal critical regions suitable for MQF-PCR primer design (Table 1, GRCh37 human genome assembly) [6–12, 14]. These regions were used as queries in Basic Local Alignment Search Tool (http://blast.ncbi.nlm.nih.gov/). A single PCR primer pair (arrows in Figure 1(a)) was then designed, capable of amplification of two fragments, one that locates to the critical genomic region of disease and another, shorter or longer, that serves as CNV-free control fragment (Figure 1(a); Table 1). In order to reliably quantify the differently sized amplicons, we employed a cost-efficient method of fluorescent labeling of PCR products with the common M13-40-NED primer (Figure 1(a)) . PCR amplification, amplicon separation, and data collection and analysis were performed as described in Stofanko et al. . We used Receiver Operating Characteristic (ROC) analysis of the Z or UPZ score statistics [15, 16] using MedCalc v12 (MedCalc Software, Mariakerke, Belgium). Box-and-Whisker plots were also constructed using MedCalc . 95% confidence intervals of true positives, false positives, true negatives and false negatives were calculated using the Fisher’s exact test  in GraphPad Prism 6 (GraphPad software, CA, USA).
3. Results and Discussion
In affected individuals, because of the deletions or duplications, the amount of template DNA corresponding to the syndrome’s critical region is inherently halved or doubled, respectively, in comparison with the control region when equal amounts of DNA from the affected and control individuals are dispensed in the assay system. Following PCR, this translates into an approximately 50% reduction or increase in the area under the electrophoretic peak of the amplicon from the critical microdeletion region (Figures 1(b)–1(i)). For example, in the detection of the WAGR syndrome (Figure 1(b)), we observed in normal controls a ratio close to 1 : 1 of the peak that corresponds to the WAGR critical region at 11p13 (black) and the control peak at 17q22 (white). In a hemizygous patient, the peak area ratio changed to approximately 1 : 2 (black/white), which reflects the template DNA dosage change. Similar reduction in peak area ratios was noted for Smith-Magenis, Langer-Giedion, Greig cephalopolysyndactyly, Miller-Dieker lissencephaly, and Saethre-Chotzen syndromes even if the peak area ratio in controls was not exactly the expected 1 : 1 ratio, presumably because amplicon-specific factors influence the efficiency of PCR (Figures 1(c) and 1(e)–1(i)). Figure 1(d) shows detection of autosomal duplication of 17p11.2 that causes Potocki-Lupski syndrome. In this patient, the peak area ratio changed from about 1 : 1 to about 3 : 2 (black/white), which reflects the template DNA dosage change. Since the number of patients available for testing is limited in rare syndromes, in the case of the Xp22.3 deletions (Figure 1(e)) we used female samples (XX) as controls and male DNA (XY) as a simulated model for affected patients. This way we could obtain a better picture of the performance of our method on a larger scale and extend its use to include diagnosis of deletions in the Xp22.3 region in male patients for whom only the control fragment would be detected. Finally, in order to correct for the observed uneven amplification efficiency, we normalized our data to a reference sample of mixed DNA from 100 individuals that was processed at the same time . Standard deviation Z-scores or pseudo-Z scores  were then calculated for each sample (Table 2), and the performance of each assay was evaluated by ROC analysis  (Figures 2(a)–2(d)). Instances with fewer than five positive patients were assessed using Box-and-Whisker plot  (Figures 2(e)–2(h)). We achieved unbiased and reproducible PCR amplification in WAGR and Smith-Magenis syndromes and the simulated Xp22.3 deletions with 100% detection sensitivity and 100% exclusion specificity. Similarly, the results for patients with rarer syndromes appeared as outliers in corresponding Box-and-Whisker plots. The 95% confidence intervals for the sensitivity and specificity of MQF-PCR tests for WAGR, Smith-Magenis, Potocki-Lupski, and the simulated Xp22.3 syndromes are shown in Table 3.
It should be noted that a limitation of this protocol for diagnosing patients with 11p13 region deletions is that it is 100% sensitive only in patients with aniridia because the targeted critical region lies between PAX6 (the gene responsible for aniridia) and WT1 (the gene responsible for Wilms tumor and genitourinary anomalies); this protocol would fail to detect 11p13 deletions that involve WT1, but spare PAX6. Such patients would not be considered to have classic WAGR syndrome however, and this MQF-PCR could be easily modified to target other critical regions to diagnose proximal 11p deletion and/or Potocki-Shaffer syndromes.
Here we report development of detection primers for 8 human microdeletion syndromes using the recently described MQF-PCR. We showed that the method is easily adaptable to detect any pathological CNV with high accuracy. With the expansion of the capability to diagnose even more disease-related CNVs, MQF-PCR appears to be a reliable, low cost, and rapid detection method, especially for use in developing countries.
The Microdeletion/Microduplication Quantitative Fluorescent PCR (MQF-PCR) test, being described herein, is the subject of a patent application (No. BR102012021281-1) by GENE-Nucleo de Genetica Médica in the Instituto Nacional de Propriedade Intelectual (INPI) of Brazil, entitled “Método e kit para identificação de microdeleções cromossômicas humanas pela quantificação fluorescente de PCR (MQF-PCR)”.
Conflict of Interests
Joan C. Han, Sarah H. Elsea, Heloísa B. Pena and Higgor Gonçalves-Dornelas declare that they have no conflict of interests.
The authors thank Conselho Nacional de Desenvolvimento Científico e Tecnológico Grant no. 401983/2010-2 for providing financial funding for this work. The research was also supported in part by the Intramural Research Program of the Eunice Kennedy Shriver National Institute of Child Health and Human Development (ZIAHD008898). We thank PRISMS (SMS support group) and the patients and their families for their participation and support of research.
M. Stofanko, H. Goncalves-Dornelas, P. S. Cunha, H. B. Pena, A. M. Vianna-Morgante, and S. D. Pena, “Simple, rapid and inexpensive quantitative fluorescent PCR method for detection of microdeletion and microduplication syndromes,” PloS ONE, vol. 8, Article ID e61328, 2013.View at: Google Scholar
L. M. B. Oliveira, S. B. Seminara, M. Beranova et al., “The importance of autosomal genes in Kallmann syndrome: genotype-phenotype correlations and neuroendocrine characteristics,” The Journal of Clinical Endocrinology and Metabolism, vol. 86, no. 4, pp. 1532–1538, 2001.View at: Publisher Site | Google Scholar
J. Hou, J. Parrish, H.-J. Ludecke et al., “A 4-megabase YAC contig that spans the Langer-Giedion syndrome region on human chromosome 8q24.1: use in refining the location of the trichorhinophalangeal syndrome and multiple exostoses genes (TRPS1 and EXT1),” Genomics, vol. 29, no. 1, pp. 87–97, 1995.View at: Publisher Site | Google Scholar
J. J. Johnston, I. Olivos-Glander, J. Turner et al., “Clinical and molecular delineation of the greig cephalopolysyndactyly contiguous gene deletion syndrome and its distinction from acrocallosal syndrome,” The American Journal of Medical Genetics A, vol. 123, no. 3, pp. 236–242, 2003.View at: Google Scholar
D. Johnson, S. W. Horsley, D. M. Moloney et al., “A comprehensive screen for TWIST mutations in patients with craniosynostosis identifies a new microdeletion syndrome of chromosome band 7p21.1,” The American Journal of Human Genetics, vol. 63, no. 5, pp. 1282–1293, 1998.View at: Publisher Site | Google Scholar
J. R. Lanzante, “Resistant, robust and non-parametric techniques for the analysis of climate data: theory and examples, including applications to historical radiosonde station data,” International Journal of Climatology, vol. 16, no. 11, pp. 1197–1226, 1996.View at: Google Scholar
M. H. Zweig and G. Campbell, “Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine,” Clinical Chemistry, vol. 39, no. 4, pp. 561–577, 1993.View at: Google Scholar
J. V. Tukey, Exploratory Data Analysis, Addison-Wesley, Reading, Mass, USA, 1977.
R. A. Fisher, Statistical Methods for Research Workers, Oliver and Boyd, Edinburgh, UK, 1925.