- Research article
- Open Access
A method for accurate detection of genomic microdeletions using real-time quantitative PCR
BMC Genomics volume 6, Article number: 180 (2005)
Quantitative Polymerase Chain Reaction (qPCR) is a well-established method for quantifying levels of gene expression, but has not been routinely applied to the detection of constitutional copy number alterations of human genomic DNA. Microdeletions or microduplications of the human genome are associated with a variety of genetic disorders. Although, clinical laboratories routinely use fluorescence in situ hybridization (FISH) to identify such cryptic genomic alterations, there remains a significant number of individuals in which constitutional genomic imbalance is suspected, based on clinical parameters, but cannot be readily detected using current cytogenetic techniques.
In this study, a novel application for real-time qPCR is presented that can be used to reproducibly detect chromosomal microdeletions and microduplications. This approach was applied to DNA from a series of patient samples and controls to validate genomic copy number alteration at cytoband 22q11. The study group comprised 12 patients with clinical symptoms of chromosome 22q11 deletion syndrome (22q11DS), 1 patient trisomic for 22q11 and 4 normal controls. 6 of the patients (group 1) had known hemizygous deletions, as detected by standard diagnostic FISH, whilst the remaining 6 patients (group 2) were classified as 22q11DS negative using the clinical FISH assay. Screening of the patients and controls with a set of 10 real time qPCR primers, spanning the 22q11.2-deleted region and flanking sequence, confirmed the FISH assay results for all patients with 100% concordance. Moreover, this qPCR enabled a refinement of the region of deletion at 22q11. Analysis of DNA from chromosome 22 trisomic sample demonstrated genomic duplication within 22q11.
In this paper we present a qPCR approach for the detection of chromosomal microdeletions and microduplications. The strategic use of in silico modelling for qPCR primer design to avoid regions of repetitive DNA, whilst providing a level of genomic resolution greater than standard cytogenetic assays. The implementation of qPCR detection in clinical laboratories will address the need to replace complex, expensive and time consuming FISH screening to detect genomic microdeletions or duplications of clinical importance.
Array comparative genomic [1–3] and FISH-based methods  have been widely used for the detection of DNA copy number changes. However, the resolution of commercially available DNA arrays can be too low to detect microdeletions or microduplications [5, 6], whilst FISH is generally only useful when the regions of interest have been previously defined. Currently, DNA arrays providing full coverage of the human genome are not widely available and too expensive to diagnostically screen large numbers of patients. Moreover, the findings that are emerging from recent array comparative genomic hybridization studies indicate that significant validations of both controls and patient populations will be required to make phenotype-genotype interpretations in a clinical context. Similarly, standard FISH methods are time consuming, costly and suffer the significant limitation that some patients with uniquely localized microdeletions or duplications may yield normal clinical FISH findings because the probe set used does not map precisely to the entire region of deletions/duplication.
In this study we have chosen to validate the use of qPCR technology for detection of microdeletions or microduplications using Velocardiofacial Syndrome (VCFS) or chromosome 22q11 deletion syndrome (22q11DS) as a test model. The frequency of the causative deletion for 22q11DS in the general population is 1 in 3000 live births, making it one of the most common constitutional genomic alterations found in humans . 22q11DS is suspected in individuals with characteristic clinical findings and is confirmed in most cases by detection of a sub-microscopic deletion using FISH.
The currently accepted clinical laboratory assay for 22q11DS uses the TUPLE1 FISH probe, which is located within a typically deleted region of approximately 3 Mb. Although this assay can detect the majority of affected patients (85–90%), many patients with phenotypic features of 22q11DS have no deletion detectable by FISH testing. As a consequence these patients will go undiagnosed due to the presence of atypical deletions that map outside the area covered by the TUPLE1 probe [8–12]. In addition, there have been reports of individuals with some features of 22q11DS with microduplications of 22q11.2 . Unfortunately, clinical FISH assays are not usually capable of detecting such duplications, so alternative methods, such as FISH analysis of interphase nuclei, are required . However, such techniques require advanced optical instrumentation, presently only used in specialized research laboratories. Alternative molecular technologies that could potentially be used to screen deletions in the diagnose of 22q11DS is multiplex ligation-dependent probe amplification (MLPA)  or microsatellites marker analysis .
The novel approach validated in this study utilized qPCR rather than FISH to detect copy number alterations (microdeletions and microduplications) in patient DNA. This approach had several advantages. Primers were selected within regions of unique sequence utilizing publicly available sequence databases [16–24]. This method can allow for the production of a high-resolution map of any region of interest; the 22q11DS region is used in this example. The qPCR technique provides a quantitative measurement of DNA copy number and accurately characterizing chromosomal breakpoints. This method will therefore permit the identification of individuals who would otherwise go undetected by the currently available clinical FISH methods. In addition, qPCR provides greater flexibility and adaptability, whilst being less technically challenging than FISH, thus making it more appropriate for use in a large number of laboratories. Furthermore, the qPCR technique only takes a fraction of the time usually required for FISH, which allows for multiple samples and multiple primer sets to be studied in parallel, using convenient and cost effective high throughput analysis. The method described in this paper has been evaluated using patient and controls samples with known copy number changes on 22q11. The approach can be readily adapted for molecular diagnostics of any region of the genome suffering recurrent constitutional genomic deletion or duplication.
The twelve patients utilized in this study have previously been analysed in a clinical laboratory by FISH using the TUPLE 1 probe. Using this test, six of the patients (group 1) were identified as 22q11 deletion positive whilst the other six (group 2) and the four controls showed normal cytogenetic results. The DNA sample known to have three copies of 22q11 (trisomy for chromosome 22) had been previously analysed by FISH .
The most critical factor for successful detection of micro-alterations using qPCR was primer design. To guarantee optimal primer design a high-resolution 22q11.2 physical map was constructed using information available from published reports [8–12, 26–34] and online databases and repositories [16–24]. This allowed for the identification of unique sequences within the 22q11DS affected region whilst also avoiding the complex repetitive regions. Figure 1 shows a schematic representation of the 22q11.2 region studied in this work, the location of previously reported deletions and deletions identified from our study are shown as are the location of primers, repeat sequences, known genes and pseudogenes.
Sequences were selected to lie, mostly, within exonic and/or intervening regions from known or putative genes: UniSTS marker D22S181 (UniSTS is a comprehensive database of sequence tagged sites) , Proline Dehydrogenase/Proline Oxidase (PRODH), TUP-like enhancer of SPLIT 1 (TUPLE1), Catechol-O-Methyltransferase (COMT), Zinc Finger Protein 74 (ZNF74), Phosphatidylinositol 4-kinase (PIK4CA), Leucine-zipper-like Transcriptional Regulator 1 (LZTR1), Cationic Amino acid Transporter-4 (CAT4), D22S936 and Similar to Pre-B lymphocyte gene 1 (VPREB1). The selected sequences were aligned against the human genome using the BLAT program  to ensure that only contiguous sequences with 100% homology to one unique location were selected. A single primer set was then designed from each of these unique sequences using the Primer Express v 2.0 program (Applied Biosystems). We designed reference primers for each of the "housekeeping genes" Glucose-6-Phosphate Dehydrogenase (G6PDH) and Hydroxymethylbilane Synthase (HEM3) which were selected using published guidelines ; moreover we insured that they were single copy genomic sequences according to the BLAT alignment , in addition agarose gel visualization confirmed a single band of the expected size. All primer sequences are shown in Table 1. The use of reference primers was to control for varying input amounts of DNA from each separate patient. Thus any differences in the qPCR values obtained for test primers/markers would correspond to differences in the amount of the primers target sequence.
Primer concentrations were optimised over a range of final concentrations, 100 nM to 900 nM at 100 nM intervals. The optimal concentration was that which obtained the lowest threshold cycle (Ct) and maximum ΔRn while minimizing non-specific amplification. The results indicated that for all of the 22q11.2 primer sets the optimal final primer concentration was 800 nM, whilst for the G6PDH and HEM3 primer sets this value was 400 nM. The specificity of amplification for each qPCR product was confirmed by determining that the melting curve, showed a single dissociation peak corresponding to the melting temperature of the analysed amplicons (See Additional File 3, Figure 2 b. Dissociation Curve for PRODH -14 DNA samples).
To allow for comparisons between each primer sets all had to amplify with comparable efficiency. This was assessed by analysis of the slopes from the standard curves, which were generated using a log10 dilution series of input genomic DNA (range 102 nM to 10-2 nM). If all conditions are optimal and reactions are 100% efficient, it will take approximately 3.32 cycles for ten fold amplification (log2 of 10 = 3.321928) of product, a value that is equal to the slope of the standard curve. This translates to 1 cycle to copy 1 molecule into 2; a second cycle to copy 2 molecules into 4; a third cycle to copy 4 into 8 and 0.32 cycles to copy 8 into 10. There was a linear relationship between the amount of input DNA and the threshold cycle (Ct) values for the various reactions. Regression analyses of the Ct values generated by the log10 dilution series gave R2 values for all reactions in excess of 0.99 (Table 1).
At the optimal primer concentrations all of the primer sets gave slope values of 3.32 ± 0.25, indicating that the reactions were occurring with similar efficiencies. Following primer optimisation a baseline Ct was identified for each primer set, which was used when analysing subsequent data (R2, slope and threshold cycle values are shown in Table 1).
To control for differences in sampling, DNA preparation, reaction efficiency (the varying PCR efficiencies between patient samples and the calibrator)  and other variables such as the Ct values for each primer pair  from 22q11.2 were corrected/normalized using the Ct value of the G6PDH and HEM3 products for the same sample. Although the input of template was standardised at 10 ng of DNA, the Ct values for the housekeeping genes differed slightly from patient to patient and from group to group (Table 3), thus demonstrating the need for correction of the raw data prior to further analysis.
Correction was performed using a method described by Moody et al. (2000) . Once the corrected Ct values (KCt) for each of the test markers had been determined it was then possible to identify fold copy number change (ΔKCt) for each of the markers from the 22q11.2 region, using a formula described by Sijben et al. (2003) . In the context of 22q11DS this approach calculates a ratio; by comparing the Ct value obtained for each primer pair between the normal (control) DNA samples and the patient (affected) DNA samples. This value is then translated into fold changes (copy number gain or loss) per sample. For all calculations please refer to the Methods section.
A summary of the qPCR results is presented in Table 2; light grey shading denotes loss, whilst dark grey shading denotes gain. We obtained ΔKCt values of either 0 ± 0.35 indicating an equal ratio of the target and reference, which corresponds to no loss and therefore no genetic abnormality, or -1 ± 0.35, indicating loss of one copy (microdeletion), for the affected samples. For the trisomy 22 patient we obtained a ratio of 1 ± 0.35, indicating gain of one copy. Each experiment was performed in triplicate, with replicates performed on different days. The inter-assay (same assay repeated on different days) Ct variation and the intra-assay Ct variation (the triplicates) was less than or equal to ± 0.35 cycles.
In our analysis we have used a series of 12 qPCR primer sets to analyze twelve patients with clinical symptoms of 22q11DS and four controls. Ten of the primer sets (test primers) amplify markers localized within and around the chromosome 22q11.2-deleted region (D22S181, PRODH, TUPLE1, COMT, ZNF74, PIK4CA, LZTR1, CAT4, D22S936 and VPREB1) and two amplify "housekeeping" genes/markers (reference primers) (G6PDH and HEM3). The 22q11DS locus contains approximately 50 genes or pseudogenes and is characterized by an unusual genomic architecture comprising a large polymorphic chromosome-specific low copy repeats (Figure 1). The repetitive nature of this region of the genome is thought to increase the frequency of deletions and duplication events associated with clinical disease. The repetitive nature of the region of 22q11.2 under study required bioinformatics to identify unique regions for primer design.
In the context of 22q11DS we were trying to discriminate between 2 copies of a product (normal) versus 1 copy of the same product (deleted). In order to do this we had to implement data correction  to control for variations such as the input DNA concentration or reaction efficiencies . To perform the correction we used the reference genes. The copy number of the reference genes will be the same in all of the samples under investigation. Any variation in copy number for the reference between samples will be the result of differences in initial template concentration, as long as the same DNA is sampled for both the control and target primers. Once the Ct values for the reference were determined for each sample, it was then possible to use these to correct the values of the test markers for variation in initial template concentration. This correction permitted determination of copy number differences between the samples under investigation.
Individuals who had previously shown a deletion of TUPLE1 by FISH also showed deletion by qPCR (group 1). Our data demonstrates that the ΔKCt values, -1 ± 0.35, for the primers PRODH, TUPLE1, COMT, ZNF74, PIK4CA, LZTR1, CAT4 and D22S936 are indicative of deletion (Table 2). A finding that is in 100% concordance with the FISH results for the TUPLE1 probe. The region of deletion spanning PRODH, to D22S936, represents an interval of 2,502,410 base pairs (bp). Furthermore, the implementation of qPCR has allowed for the identification of breakpoint within the 22q11.2 region. The markers D22S181 and VPREB1 show values indicative of no deletion (0 ± 0.35), suggesting that the proximal deletion breakpoint occurs between the markers D22S181 and PRODH and the distal breakpoint between D22S936 and VPREB1 (see Table 2 and Figure 1).
The individuals that did not demonstrate loss using the TUPLE1 FISH probe were also deletion negative by qPCR, showing ΔKCt values similar to the normal controls (0 ± 0.35) thus indicative of no loss. For the 22q11 trisomy sample the qPCR results showed ΔKCt values indicative of duplication (1 ± 0.35) for PRODH, TUPLE1, COMT, ZNF74, PIK4CA, LZTR1, CAT4 and D22S936. The markers D22S181 and VPREB1 again showed values indicative of no copy number change (0 ± 0.35).
To accurately calculate fold changes, to the high precision required for genomic DNA quantitation, the real-time data analysis used here makes the most of the existing methods. The slope values determined from the standard curve allows for correction for variations in the primers efficiency and reaction kinetics, whilst the relative abundance ratio, calculated after the samples are normalized using the reference genes, allows for determination of gene/marker copy number. This rationale reduces, as much as possible, the unavoidable approximations introduced with any method of data interpretation.
The use of qPCR to detect and refine copy number differences in patients suffering from 22q11DS provides a further novel application for qPCR methodology. When used in a research setting, this type of analysis has proven to be very useful when comparing levels of a transcript or genomic markers [43, 44] between different groups, and has been widely used in the study of human malignancies. However, such an approach has rarely been applied in a clinical setting to study constitutional deletions. The advantage of this technique over standard FISH assays is that qPCR provides a quantitative measurement of DNA copy number.
Here we demonstrate the application of a robust, fast and accurate real-time quantitative PCR based assay using SYBR ® Green I dye, that is capable of screening for copy-number alterations in genomic DNA.
Although qPCR detection methods have previously been used in 22q11.2 deletion analysis [45–47], these reports have only used a small number of primers/markers and have not been able to refine the region of deletion as has been done here. The utility of the approach outlined in this paper is the ease with which one can increase resolution by increasing the number of primers in the 22q11 deleted region thus facilitating accurate mapping of deletion breakpoints. The fine structure of qPCR mapping of deletions will reveal important clues into the mechanism by which the deletion occurs and thus will offer insights into the "at risk" factors predictive of deletions or other rearrangements.
The implementation of qPCR for genomic copy number profiling will provide a valuable tool for detection of atypical microdeletions and/or microduplications in individuals who go undiagnosed by the current available FISH methods. Such data will be useful for phenotype correlation studies. In addition, this methodology has the advantage of providing greater flexibility and adaptability than the currently available cytogenetic methods and will be beneficial in molecular classification and diagnosis.
Primers (Table 1) were designed using Primer Express v2.0 (Applied Biosystems). The parameters for primer design were as follows; amplicons of 100–250 bp with a penalty score no higher than 10–12; primer melting temperature (Tm) range 58°C to 60°C; primer length range 18 to 25 bp, optimal 20 bp; primer G/C content range 20% to 60%; amplicon maximum Tm of 85°C. Runs of more than three identical nucleotides were avoided as polyG or polyC stretches can promote non-specific annealing, whilst runs of polyA and polyT can potentially open up stretches of the primer-template complex. Primers were required to be free of self-complementary sequence in order to avoid hairpin loop formation. In addition, no more than two G and/or C bases were permitted in the last five nucleotides at the 3' end to avoid GC clamp formation. Minimal deviation from the parameters was allowed only when there were no other options due to the complex nature of the 22q11.2 sequence. Primer and amplicon sequences were compared to the human genome using the BLAT program to guarantee that they showed 100% homology to only the sequence from which they were designed and also to guarantee that the forward and reverse primers were free of single nucleotide polymorphisms.
Isolation of DNA
The DNA used in this study was obtained from the lymphocytes of four normal controls and 12 patients displaying clinical features of 22q11DS using phenol-chloroform extraction . DNA was quality tested. Suitable sample for the qPCR reactions were of high molecular weight (un-sheared band of undigested DNA visible on a 0.5% agarose gel) and as clean as possible (an OD 260/280 ranging from 1.8 to 2.0).
Reactions were performed using SYBR Green I PCR Master Mix (Applied Biosystems), which includes the internal reference (ROX). Each qPCR reaction comprised 12.5 μl 2× SYBR Green PCR Master Mix, forward and reverse primer at optimized concentrations of 800 nM (final concentration) for the 22q11.2 test primers and 400 nM (final concentration) for the reference primers, 10 ng/μl genomic DNA template and sterile water up to a final volume of 25 μl. The qPCR reactions were performed using the ABI Prism 7900 high-throughput sequence detection system. The reaction profile was: initial step, 50°C for 2 min, denaturation, 95°C for 10 min, then 40 cycles of denaturing at 95°C for 15 sec and combined annealing and extension at 60°C for 60 sec.
Generating the standard curve
To generate standard curves for the selected primers and the reference primers a log10 dilution series of genomic DNA was prepared at concentrations ranging from 102 nM to 10-2 nM. Each dilution was tested in triplicate. When analyzed by qPCR, the dilution series produced a set of standard curves, which were used to calculate the slope value with the aid of the SDS software version 2.1 Applied Biosystems (values are shown in Table 1). See Additional File 3, Figure 2 e. and f. for an example of SDS output report of standard and amplification plots.
DNA quantification data analysis
Each qPCR experiment contained triplicates of the no-template-controls and patient samples for all of the primers tested. On the same reaction plate all DNA samples were tested with the test and reference primers. When any particular sample was being tested, the qPCR using reference primers and that sample were always included on the reaction plate. Each experiment was performed in triplicate, with replicates being performed on different days. Quantification was based on the increased fluorescence, which was measured and recorded using the ABI Prism 7900 sequence detection system and associated SDS software version 2.1 (Applied Biosystems). Results were expressed in terms of the threshold cycle value (Ct; the cycle at which the change in fluorescence for the SYBR dye passes a significance threshold). The threshold values are shown in Table 1. The output of the results was exported in tab-delimited text file format. Further calculations were performed using Microsoft Excel. PCR products were resolved by agarose gel electrophoresis to confirm the presence of a single band of the expected size (See Additional File 3, Figure 2 a. showing SDS output of amplification plot for 14 DNA samples for the PRODH and G6PDH primer sets along with images of corresponding gel bands c. and d.).
The qPCR data was normalized adapting a method devised by Moody et al. (2000)  and also described by Sijben et al. (2003) . (See Additional File 1: Derivation of the formula)
KCti = 'Corrected C t ' (KC t ) of the test primer (T) against the reference (R)
ACtR = the Average Ct value for Reference primer set for all the samples included in one qPCR run (control and patient).
CtRi = Ct value for Reference primer set for the sample to be corrected.
SR = slope value (from the standard curve) for reference primer set.
ST = slope value (from the standard curve) for test primer set.
CtTi = Ct value for test primer set.
Copy number calculation
Fold copy number (ΔKCt) change for each of the markers from the 22q11.2 region, was obtained using the formula:
ΔKC t = KCt/control- KCt/affected
ΔKCt = fold change (copy number gain or loss)
KCt/control = "Corrected Ct" of the test primer for the control samples.
KCt/affected = "Corrected Ct" of the test primer for the affected sample.
(Of note when multiple controls were used in the same reaction run the KCt/control was obtained by averaging all the controls'
"Corrected Ct" – KCt/control's). (See Additional File 2: Working example)
- qPCR –:
Real time Quantitative Polymerase Chain Reaction
- G6PDH –:
- Reference (R) –:
refers to the sequence we use to standardize against – the sequence which we expect to remain unchanged in all the DNA samples (patient and control), and is the equivalent of the housekeeping gene use in expression studies. In our experiments we use G6PDH and HEM3.
- Target (T) –:
refers to the sequence originating from the region of interest which we test and which we expect to find modifications (deletion or duplication).
Ishkanian AS, Malloff CA, Watson SK, DeLeeuw RJ, Chi B, Coe BP, Snijders A, Albertson DG, Pinkel D, Marra MA, Ling V, MacAulay C, Lam WL: A tiling resolution DNA microarray with complete coverage of the human genome. Nat Genet. 2004, 36: 299-303. 10.1038/ng1307.
Koolen DA, Vissers LE, Nillesen W, Smeets D, van Ravenswaaij CM, Sistermans EA, Veltman JA, de Vries BB: A novel microdeletion, del(2)(q22.3q23.3) in a mentally retarded patient, detected by array-based comparative genomic hybridization. Clin Genet. 2004, 65: 429-432. 10.1111/j.0009-9163.2004.00245.x.
Shaw-Smith C, Redon R, Rickman L, Rio M, Willatt L, Fiegler H, Firth H, Sanlaville D, Winter R, Colleaux L, Bobrow M, Carter NP: Microarray based comparative genomic hybridisation (array-CGH) detects submicroscopic chromosomal deletions and duplications in patients with learning disability/mental retardation and dysmorphic features. J Med Genet. 2004, 41: 241-248. 10.1136/jmg.2003.017731.
Babovic-Vuksanovic D, Jenkins SC, Ensenauer R, Newman DC, Jalal SM: Subtelomeric deletion of 18p in an adult with paranoid schizophrenia and mental retardation. Am J Med Genet A. 2004, 124: 318-322. 10.1002/ajmg.a.20391.
Mantripragada KK, Buckley PG, de Stahl TD, Dumanski JP: Genomic microarrays in the spotlight. Trends Genet. 2004, 20: 87-94. 10.1016/j.tig.2003.12.008.
Mantripragada KK, Tapia-Paez I, Blennow E, Nilsson P, Wedell A, Dumanski JP: DNA copy-number analysis of the 22q11 deletion-syndrome region using array-CGH with genomic and PCR-based targets. Int J Mol Med. 2004, 13: 273-279.
Saitta SC, Harris SE, Gaeth AP, Driscoll DA, McDonald-McGinn DM, Maisenbacher MK, Yersak JM, Chakraborty PK, Hacker AM, Zackai EH, Ashley T, Emanuel BS: Aberrant interchromosomal exchanges are the predominant cause of the 22q11.2 deletion. Hum Mol Genet. 2004, 13: 417-428. 10.1093/hmg/ddh041.
Kurahashi H, Tsuda E, Kohama R, Nakayama T, Masuno M, Imaizumi K, Kamiya T, Sano T, Okada S, Nishisho I: Another critical region for deletion of 22q11: a study of 100 patients. Am J Med Genet. 1997, 72: 180-185. 10.1002/(SICI)1096-8628(19971017)72:2<180::AID-AJMG10>3.0.CO;2-J.
Amati F, Conti E, Novelli A, Bengala M, Diglio MC, Marino B, Giannotti A, Gabrielli O, Novelli G, Dallapiccola B: Atypical deletions suggest five 22q11.2 critical regions related to the DiGeorge/velo-cardio-facial syndrome. Eur J Hum Genet. 1999, 7: 903-909.
Lu JH, Chung MY, Betau H, Chien HP, Lu JK: Molecular characterization of tetralogy of fallot within Digeorge critical region of the chromosome 22. Pediatr Cardiol. 2001, 22: 279-284. 10.1007/s002460010219.
Rauch A, Pfeiffer RA, Leipold G, Singer H, Tigges M, Hofbeck M: A novel 22q11.2 microdeletion in DiGeorge syndrome. Am J Hum Genet. 1999, 64: 659-666. 10.1086/302235.
Saitta SC, McGrath JM, Mensch H, Shaikh TH, Zackai EH, Emanuel BS: A 22q11.2 deletion that excludes UFD1L and CDC45L in a patient with conotruncal and craniofacial defects. Am J Hum Genet. 1999, 65: 562-566. 10.1086/302514.
Ensenauer RE, Adeyinka A, Flynn HC, Michels VV, Lindor NM, Dawson DB, Thorland EC, Lorentz CP, Goldstein JL, McDonald MT, Smith WE, Simon-Fayard E, Alexander AA, Kulharya AS, Ketterling RP, Clark RD, Jalal SM: Microduplication 22q11.2, an emerging syndrome: clinical, cytogenetic, and molecular analysis of thirteen patients. Am J Hum Genet. 2003, 73: 1027-1040. 10.1086/378818.
Ikeuchi T: Extending the capabilities of human chromosome analysis: from high-resolution banding to chromatin fiber-FISH. Hum Cell. 1997, 10: 121-134.
Sellner LN, Taylor GR: MLPA and MAPH: new techniques for detection of gene deletions. Hum Mutat. 2004, 23: 413-419. 10.1002/humu.20035.
UCSC Genome Bioinformatics. [http://genome.ucsc.edu/]
Roe Lab Human Genomic DNA Sequencing Progress. [http://www.genome.ou.edu/hum_totals.html]
Human BAC Resource. [http://www.ncbi.nlm.nih.gov/genome/cyto/hbrc.shtml]
Human Organized Whole Genome Database. [http://www-alis.tokyo.jst.go.jp/HOWDY/]
YACs on Chromosome 22. [http://www.biologia.uniba.it/rmc/2-YAC-BAC/YAC/YAC-22.html]
Human Genome Segmental Duplication Database. [http://projects.tcag.ca/humandup/]
Yobb TM, Somerville MJ, Willatt L, Firth HV, Harrison K, Mackenzie J, Gallo N, Morrow BE, Shaffer LG, Babcock M, Chernos J, Bernier F, Sprysak K, Christiansen J, Haase S, Elyas B, Lilley M, Bamforth S, McDermid HE: Microduplication and triplication of 22q11.2: a highly variable syndrome. Am J Hum Genet. 2005, 76: 865-876. 10.1086/429841.
Carlson C, Sirotkin H, Pandita R, Goldberg R, McKie J, Wadey R, Patanjali SR, Weissman SM, Anyane-Yeboa K, Warburton D, Scambler P, Shprintzen R, Kucherlapati R, Morrow BE: Molecular definition of 22q11 deletions in 151 velo-cardio-facial syndrome patients. Am J Hum Genet. 1997, 61: 620-629.
Demczuk S, Aledo R, Zucman J, Delattre O, Desmaze C, Dauphinot L, Jalbert P, Rouleau GA, Thomas G, Aurias A: Cloning of a balanced translocation breakpoint in the DiGeorge syndrome critical region and isolation of a novel potential adhesion receptor gene in its vicinity. Hum Mol Genet. 1995, 4: 551-558.
Funke B, Edelmann L, McCain N, Pandita RK, Ferreira J, Merscher S, Zohouri M, Cannizzaro L, Shanske A, Morrow BE: Der(22) syndrome and velo-cardio-facial syndrome/DiGeorge syndrome share a 1.5-Mb region of overlap on chromosome 22q11. Am J Hum Genet. 1999, 64: 747-758. 10.1086/302284.
Kurahashi H, Akagi K, Inazawa J, Ohta T, Niikawa N, Kayatani F, Sano T, Okada S, Nishisho I: Isolation and characterization of a novel gene deleted in DiGeorge syndrome. Hum Mol Genet. 1995, 4: 541-549.
Lindsay EA: Chromosomal microdeletions: dissecting del22q11 syndrome. Nat Rev Genet. 2001, 2: 858-868. 10.1038/35098574.
Lindsay EA, Vitelli F, Su H, Morishima M, Huynh T, Pramparo T, Jurecic V, Ogunrinu G, Sutherland HF, Scambler PJ, Bradley A, Baldini A: Tbx1 haploinsufficieny in the DiGeorge syndrome region causes aortic arch defects in mice. Nature. 2001, 410: 97-101. 10.1038/35065105.
McQuade L, Christodoulou J, Budarf M, Sachdev R, Wilson M, Emanuel B, Colley A: Patient with a 22q11.2 deletion with no overlap of the minimal DiGeorge syndrome critical region (MDGCR). Am J Med Genet. 1999, 86: 27-33. 10.1002/(SICI)1096-8628(19990903)86:1<27::AID-AJMG6>3.0.CO;2-7.
Perez E, Sullivan KE: Chromosome 22q11.2 deletion syndrome (DiGeorge and velocardiofacial syndromes). Curr Opin Pediatr. 2002, 14: 678-683. 10.1097/00008480-200212000-00005.
Shaikh TH, Kurahashi H, Saitta SC, O'Hare AM, Hu P, Roe BA, Driscoll DA, McDonald-McGinn DM, Zackai EH, Budarf ML, Emanuel BS: Chromosome 22-specific low copy repeats and the 22q11.2 deletion syndrome: genomic organization and deletion endpoint analysis. Hum Mol Genet. 2000, 9: 489-501. 10.1093/hmg/9.4.489.
Johnson A, Minoshima S, Asakawa S, Shimizu N, Shizuya H, Roe BA, McDermid HE: A 1.5-Mb contig within the cat eye syndrome critical region at human chromosome 22q11.2. Genomics. 1999, 57: 306-309. 10.1006/geno.1999.5757.
Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664. 10.1101/gr.229202. Article published online before March 2002.
Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun. 2004, 313: 856-862. 10.1016/j.bbrc.2003.11.177.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.
Meijerink J, Mandigers C, van de Locht L, Tonnissen E, Goodsaid F, Raemaekers J: A novel method to compensate for different amplification efficiencies between patient DNA samples in quantitative real-time PCR. J Mol Diagn. 2001, 3: 55-61.
Tichopad A, Dilger M, Schwarz G, Pfaffl MW: Standardized determination of real-time PCR efficiency from a single reaction set-up. Nucleic Acids Res. 2003, 31: e122-10.1093/nar/gng122.
Moody A, Sellers S, Bumstead N: Measuring infectious bursal disease virus RNA in blood by multiplex real-time quantitative RT-PCR. J Virol Methods. 2000, 85: 55-64. 10.1016/S0166-0934(99)00156-1.
Sijben JW, Klasing KC, Schrama JW, Parmentier HK, van der Poel JJ, Savelkoul HF, Kaiser P: Early in vivo cytokine genes expression in chickens after challenge with Salmonella typhimurium lipopolysaccharide and modulation by dietary n--3 polyunsaturated fatty acids. Dev Comp Immunol. 2003, 27: 611-619. 10.1016/S0145-305X(03)00031-4.
Parr C, Jiang WG: Quantitative analysis of lymphangiogenic markers in human colorectal cancer. Int J Oncol. 2003, 23: 533-539.
Cairns JP, Chiang PW, Ramamoorthy S, Kurnit DM, Sidransky D: A comparison between microsatellite and quantitative PCR analyses to detect frequent p16 copy number changes in primary bladder tumors. Clin Cancer Res. 1998, 4: 441-444.
Ivanov D, Kirov G, Norton N, Williams HJ, Williams NM, Nikolov I, Tzwetkova R, Stambolova SM, Murphy KC, Toncheva D, Thapar A, O'Donovan MC, Owen MJ: Chromosome 22q11 deletions, velo-cardio-facial syndrome and early-onset psychosis. Molecular genetic study. Br J Psychiatry. 2003, 183: 409-413. 10.1192/bjp.183.5.409.
Kariyazono H, Ohno T, Ihara K, Igarashi H, Joh-o K, Ishikawa S, Hara T: Rapid detection of the 22q11.2 deletion with quantitative real-time PCR. Mol Cell Probes. 2001, 15: 71-73. 10.1006/mcpr.2000.0340.
Shi YR, Hsieh KS, Wu JY, Lee CC, Tsai CH, Yu MT, Chang JS, Tsai FJ: Genetic analysis of chromosome 22q11.2 markers in congenital heart disease. J Clin Lab Anal. 2003, 17: 28-35. 10.1002/jcla.10062.
Nicholas C. Dracopoli, Jonathan L. Haines, Bruce R. Korf, Cynthia C. Morton, Christine E. Seidman, J.G. Seidman, Smith DR: Short Protocols in Human Genetics. Edited by: Leonard SA. 1994, , John Wiley & Sons, Inc.
We would like to thank Dr. H. McDermid from University of Alberta, Canada for providing a cell line carrying a duplication of the 22q11.2 region for this study . We would like to recognize the invaluable input from Dr. I. Teshima and for their technical assistance the researchers at The Hospital for Sick Children, Toronto, Canada and also Dr. Sandy Der's lab at The University of Toronto. Also we would like to thank Dr. G. Underwood (Bioproperties, Australia) Dr. C. Jones (CAB International, UK) and Dr. A. Moody (Astra Zeneca, UK) for their critical review of this manuscript. This research was supported by the National Cancer Institute of Canada with funds from the Canadian Cancer Society.
JS and RW conceived the study, and participated in its design, coordination and final writing; SH resourced the mathematical model and the methodology, revised the data and the manuscript for publication; LM developed the required bioinformatics, carried out the molecular genetics study, performed the data interpretation adapting the mathematical model, additional materials and, together with SH drafted the manuscript; ASB and EWCC contributed microdeletion cell lines, pertinent clinical information and contributed to overall study design; all authors revised and approved the final manuscript.
Rosanna Weksberg, Simon Hughes contributed equally to this work.
Electronic supplementary material
Additional File 1: MS Word 2000 document demonstrating the mathematical derivation of the formula for copy number transition (fold change). (DOC 42 KB)
Additional File 2: MS Word 2000 document describing the mathematical formula applied to a practical example (working example) where fold change is calculated from the Ct values of SDS output in excel. (DOC 96 KB)
Additional File 3: Figure 2 (Fig. 2) Example of SDS output report and agarose gel for samples run on real-time qPCR. a. Amplification plot for PRODH (14 DNA samples) and G6PDH (14 DNA samples). b. Dissociation Curve PRODH -14 DNA samples. c. and d. Agarose gel electrophoresis images showing unique bands for 14 DNA samples qPCR products corresponding to 101 bp fragment (PRODH primer set) and respectively 202 bp – G6PDH. e. COMT Standard Curve Plot and f. Amplification Plot for standard dilutions. (PDF 724 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Weksberg, R., Hughes, S., Moldovan, L. et al. A method for accurate detection of genomic microdeletions using real-time quantitative PCR. BMC Genomics 6, 180 (2005). https://doi.org/10.1186/1471-2164-6-180
- Copy Number Change
- Copy Number Gain
- Chromosomal Microdeletions
- Copy Number Difference
- Log10 Dilution Series
View archived comments (1)