| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |



From the Department of Pathology,
*
Brigham and Womens Hospital, Boston; the Department of Medical Oncology,
Dana-Farber Cancer Institute, Boston; Harvard Medical School,
Boston, Massachusetts; and the Department of Pathology and Laboratory Medicine,
University of Pennsylvania Medical Center, Philadelphia, Pennsylvania
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Approximately 94% of clinically typical SMA patients lack both copies of SMN1 exon 7.8
SMN gene dosage analysis, a method to determine the copy number of SMN1, can be used to identify SMA carriers. Exon 7 of SMN1 and SMN2 are co-amplified with genomic and internal standards. The PCR products are then digested with DraI, which cuts only SMN2 exon 7 PCR products, followed by quantification of the PCR products.9, 10, 11
Other methods for SMN gene dosage analysis have been described.12, 13, 14
A single copy of SMN1 by gene dosage analysis confirms carrier status; this analysis is therefore of clinical importance. A single-copy result also supports the diagnosis of SMN1-related SMA in an affected individual, who may have one deleted allele and one allele with a small intragenic mutation. However, the final diagnosis depends largely on the index of clinical suspicion.15
This is because the frequency of single-copy carriers in the general population (
2%) approaches the frequency of individuals affected with SMN1-related SMA who have a single-copy test result (
3.6%12
).16
The copy number of SMN2 correlates inversely with disease severity.9, 10, 12, 13, 14 Feldkötter et al14 found that SMN2 copy number also correlates directly with length of survival. Potential therapies for SMA include approaches to increase the expression of full-length transcripts from SMN2. Full-length SMN2 transcripts are increased in vitro and in vivo by sodium butyrate,17 and in vitro by aclarubicin.18 Hence, in the future, accurate determination of SMN2 copy number may have both prognostic and therapeutic significance.
We described previously the effect of SMN1/SMN2 heteroduplex formation on SMN gene dosage analysis.11 We calculated that unless SMN2 is absent, apparent SMN1 peaks contain between approximately 2.9% and approximately 14% SMN2 PCR products (depending on the genotype) due to DraI-undigestable SMN1/SMN2 heteroduplexes.11 However, in our experience, there seemed to be less SMN2 PCR product than SMN1 PCR product in almost all samples, even after correcting for heteroduplex formation. The hypothesis that incomplete DraI digestion falsely increased the SMN1 signal and decreased the SMN2 signal seemed unlikely because of assay controls lacking SMN1, in which undigested SMN2 signal has never been detected.11 We hypothesized that there might be a considerable difference in PCR efficiency (PCR bias) between SMN1 and SMN2. Using a large number of samples in our SMN gene dosage analysis, a robust quantitative PCR assay, we quantify herein consistent PCR bias caused by the single nucleotide polymorphism (SNP) between SMN1 exon 7 and SMN2 exon 7. We also validate methods to determine SMN2 copy number. The precision and accuracy of our SMN gene dosage analysis are high because our assay design and controls take advantage of the consistency of the PCR bias.
PCR bias caused by an SNP, not in primer binding sites, can significantly affect the accuracy and precision of quantitative PCR assays. To assure high precision and accuracy of quantitative PCR assays, standards and controls must be chosen judiciously, and the signal intensities of the PCR products must be calculated and normalized appropriately. In addition, close monitoring of results on clinical samples and controls/standards is essential for quality assurance in any clinical molecular diagnostic laboratory that performs quantitative PCR. As additional clinically significant SNPs are discovered, assessment of PCR bias will be increasingly important.
| Materials and Methods |
|---|
|
|
|---|
SMN1 and SMN2 Copy Number Assay (SMN Gene Dosage Analysis)
SMN1 gene dosage analysis was originally developed by McAndrew et al9
and modified as a non-radioisotopic assay as described previously.10
The assay has since been modified further, using 23 cycles of PCR and the ABI Prism 310 Genetic Analyzer (Applied Biosystems, Foster City, CA).11
All test samples were analyzed in duplicate using the same PCR master mix. The assay takes advantage of the SNP in exon 7 to distinguish SMN1 from SMN2 after DraI digestion. The copy number of SMN1 per cell (or, more precisely, per diploid genome) was determined as described previously.11
Briefly, we first normalize the SMN1 signal of each sample, using both a genomic standard (CFTR exon 4), as well as internal standards for SMN1 and CFTR that are added to the PCR reaction. We then normalize the result to the mean of five control samples, each with two copies of SMN1, to obtain the SMN1 value designated as "C(SMN1)," (Table 1)
which stands for "calculated SMN1,"11
as described in the Appendix. Theoretically, C(SMN1) should be close to an integer number, and should indicate the copy number of SMN1. The use of control samples with two copies of SMN1 were validated as described.10, 11
|
The coefficient of variation (CV) between the two C(SMN1) values for each sample in the duplicate testing of our 272 samples ranged from 0% to 12%, with a mean of 2.5% and a median of 1.9%.11 The CV between the two N(SMN2) values for each sample in the duplicate testing of our 259 samples (excluding the 2:0 genotype) ranged from 0% to 20%, with a mean of 3.4% and a median of 2.6%. The average values of C(SMN1), C(SMN2) (see below) and N(SMN2) from the two runs for each sample were used for further analyses.
Quantification of PCR Bias between SMN1 and SMN2
To quantify PCR bias between SMN1 and SMN2, we calculated the SMN2 signal relative to that of SMN1 ["C(SMN2)"], using the same set of external quantification standards as C(SMN1) (see Appendix for details). To calculate PCR bias between SMN1 and SMN2 accurately, we normalized the SMN1 and SMN2 signals to the same standards. We refer to the means of C(SMN1), C(SMN2), and N(SMN2) for a given SMN1:SMN2 genotype X:Y as MC(SMN1)X:Y, MC(SMN2)X:Y, and MN(SMN2)X:Y, respectively. Because SMN1/SMN2 heteroduplexes cannot be digested by DraI, they falsely increase the SMN1 signal and falsely decrease the SMN2 signal.11
Thus, we first needed to correct our values of MC(SMN1)X:Y and MC(SMN2)X:Y for the false increase of the SMN1 signal, and the false decrease of the SMN2 signal, caused by SMN1/SMN2 heteroduplex formation. Because SMN1/SMN2 heteroduplexes cannot form in samples of the 2:0 genotype, we quantified SMN1/SMN2 heteroduplex formation by comparing the MC(SMN1)X:Y to the MC(SMN1)2:0.11
The amount of SMN1/SMN2 heteroduplexes is expressed as the percentage of the SMN1 signal MC(SMN1)X:Y that is contributed by SMN2 (which we refer to as "%SMN2X:Y " for a genotype X:Y). The extent of heteroduplex formation depends on the ratio of SMN1 copies to SMN2 copies, and therefore must be quantified separately for each genotype.11
We re-analyzed the data of Ogino et al11
after changing the genotype assignments from 1:3 to 1:4 for two samples (see below), resulting in minor changes in the values for %SMN2X:Y in these genotypes (Table 2)
. After taking into account heteroduplex formation, we quantified PCR bias between SMN1 and SMN2 (referred to as "PBX:Y " for a given genotype X:Y) for each genotype, by dividing the corrected value for SMN2 per SMN2 copy by the corrected value for SMN1 per SMN1 copy. Our methods for calculating the PCR bias (PBX:Y) are described in detail in the Appendix.
|
| Results |
|---|
|
|
|---|
|
To determine the consistency of SMN1/SMN2 heteroduplex formation and PCR bias between SMN1 and SMN2, we calculated the uncorrected ratio of C(SMN2)/Y to C(SMN1)/X for each sample in a given genotype X:Y. The mean ratios for each genotype were (mean ± SD): 0.651 ± 0.083 for 2:3, 0.729 ± 0.016 for 3:2, 0.723 ± 0.061 for 1:3, 0.694 ± 0.065 for 2:2, 0.681 ± 0.067 for 3:1, 0.698 ± 0.056 for 1:2, 0.687 ± 0.075 for 2:1, and 0.712 ± 0.050 for 1:1. Thus, heteroduplex formation and PCR bias were reproducible between samples.
| Discussion |
|---|
|
|
|---|
Quantification of PCR bias due to an SNP is complicated by heteroduplex formation between the two sequences, and by the need to quantify nearly identical PCR products independently. We present herein methods to overcome these difficulties. We describe the quantification of PCR bias due to an SNP using a large number of samples in our robust SMN gene dosage analysis. The consistency of our PCR bias measurements between samples and between genotypes reinforces the validity of our results. We used samples with various copy numbers of SMN1 and of SMN2, which were present in precise integer ratios. In addition, the simultaneous amplification of a genomic reference sequence (CFTR exon 4) allowed us to normalize initial template input.
We also evaluated the effect of PCR bias on our SMN gene dosage analysis,10, 11 which is based on the method of McAndrew et al9 A priori, we assumed that the amplification efficiencies for SMN1 and SMN2 would be nearly identical since the single nucleotide difference in the segment amplified is not in the primer binding sites. The data presented herein demonstrate that this assumption was incorrect. The PCR bias between SMN1 and SMN2, and any other PCR bias that may occur between the various genomic and internal-standard sequences in our assay, was consistent between samples and between runs, which allows us to maintain assay precision. Because SMN2 amplifies approximately 20% less efficiently than SMN1 in our assay, we normalize SMN2 signals using SMN2 signals from controls of known SMN2 copy number. For SMN gene dosage analysis, the laboratory should verify that the ratio of apparent SMN2 signal to apparent SMN1 signal is consistent in each SMN genotype. In our assay, this ratio is approximately 70% because both PCR bias and heteroduplex formation cause an apparent increase in the SMN1 signal and an apparent decrease in the SMN2 signal.
The cause of PCR bias due to an SNP or a point mutation not in primer binding sites is poorly understood. Bernard et al19 hypothesized that the sequence CXGG might cause PCR bias. The segment of SMN1 and SMN2 amplified in our assay contains the sequence ... CAGGGTTT(C or T)A(G to A)ACAA...(where the reverse-primer binding site is italicized, "C or T " refers to the SMN exon 7 polymorphism, and "G to A " refers to the nucleotide change generated by primer mismatch to create a DraI site in SMN2 PCR product). The CXGG sequence is present, though not at the site of the polymorphism. One may hypothesize that a difference in the exact locations of SMN1 and SMN2 on 5q13 might cause PCR bias. However, if that were true, converted telomeric SMN2 (or centromeric SMN1, if it exists) would amplify with an efficiency similar to that of native telomeric SMN1 (or centromeric SMN2, respectively). Samples in our study with converted telomeric SMN2 would have had less or no PCR bias, since at least one SMN2 copy would have been amplified with a similar efficiency to that of the native SMN1. The consistency and reproducibility of our PCR bias data do not support this hypothesis. Moreover, it seems unlikely that a particular chromosomal structure would be present in our purified DNA samples. Alternatively, because the polymorphism lies only one nucleotide from the 3' end of the reverse-primer binding site, it might affect the initial interaction of the DNA polymerase with the template and dNTP. The SNP might also affect initial primer binding.
Heteroduplex formation should diminish as PCR cycle number decreases because it depends on the amount of PCR products formed.11 In contrast, PCR bias may be present even after the initial cycles of amplification, which could affect quantification in real-time PCR assays. Recently, Feldkötter et al14 described a real-time PCR assay for the quantification of SMN1 and SMN2 copy numbers. They used allele-specific PCR, with slightly different primer pairs for SMN1 and SMN2. In their assay, SMN2 amplified somewhat better than SMN1. The precision and accuracy of the method of Feldkötter et al14 are similar to, or slightly lower than, those of our method (10, 11 and data herein), which is a modification of the original method of McAndrew et al9 All of these methods,9, 10, 11, 14 in addition to that of Gérard et al,13 determine SMN1 and SMN2 copy numbers reliably. Gérard et al13 also demonstrated an efficiency bias in their primer-extension assay, slightly in favor of SMN2, even though they generated larger products from SMN2 (27 bp) than those from SMN1 (23 bp). In contrast, they found less or no PCR bias between the SMN sequence and its 3 bp-smaller SMN internal standard, and between their genomic reference (PBGD) and its 5 bp-larger internal standard.13 Thus PCR bias appears both primer- and template-specific.
In conclusion, even a single nucleotide difference, not in primer binding sites, can cause reproducible PCR bias. The precision and accuracy of our SMN gene dosage analysis are high because our assay design and controls take advantage of the consistency of the PCR bias. As additional clinically significant SNPs in the human genome are discovered, assessment of PCR bias, and judicious selection of standards and controls, will be increasingly important for quantitative PCR assays.
| Appendix 1 |
|---|
|
|
|---|
A1, SMN internal standard PCR product;
A2, SMN2 PCR product;
A3, SMN1 PCR product;
A4, CFTR internal standard PCR product;
A5, CFTR PCR product.
As described previously by Ogino et al,11
![]() | (1) |
![]() | (2) |
From Equation (1)
![]() | (3) |
![]() | (4) |
is normalized as follows:
![]() | (5) |
We obtained a weighted average of k4, designated as k4*, from Equation 5
in seven control samples comprising three with the 2:2 genotype, two with the 2:1 genotype, one with the 1:2 genotype, and one with the 1:1 genotype: ie, the sum of
was divided by 11, the total number of copies of SMN2 in the seven samples.
From Equation 4
![]() | (6) |
Equation 6 is designated herein as the "normalized SMN2 signal" or N(SMN2).
2. Measurement of PCR Bias between SMN1 and SMN2 (PBX:Y)
We defined the "calculated SMN2 signal relative to SMN1" or "C(SMN2)" as follows {using k3*, which is defined above for the calculation of C(SMN1)}:
![]() | (7) |
We define the mean C(SMN1) values and the mean C(SMN2) values for each genotype X:Y as MC(SMN1)X:Y and MC(SMN2)X:Y, respectively. We quantified the fraction of SMN2 products derived from SMN1/SMN2 heteroduplexes in an SMN1 signal in the genotype X:Y as described previously,11
and designate this fraction as %SMN2X:Y. The corrected (for heteroduplex formation) signal intensity for SMN1, per copy of SMN1, is:
![]() | (8) |
Likewise, the corrected (for heteroduplex formation) signal intensity for SMN2, per copy of SMN2, is:
![]() | (9) |
Thus, the difference in amplification efficiency (PCR bias) between SMN1 and SMN2 in a genotype X:Y, which we designate as "PBX:Y", is defined as the corrected (for heteroduplex formation) SMN2 signal intensity per copy of SMN2, divided by the corrected (for heteroduplex formation) SMN1 signal intensity per copy of SMN1, or (from Equations 8
and 9
):
![]() | (10) |
| Acknowledgments |
|---|
| Footnotes |
|---|
Accepted for publication August 29, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Ogino, T. Kawasaki, M. Brahmandam, M. Cantor, G. J. Kirkner, D. Spiegelman, G. M. Makrigiorgos, D. J. Weisenberger, P. W. Laird, M. Loda, et al. Precision and Performance Characteristics of Bisulfite Conversion and Real-Time PCR (MethyLight) for Quantitative DNA Methylation Analysis J. Mol. Diagn., May 1, 2006; 8(2): 209 - 217. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ogino, T. Kawasaki, M. Brahmandam, L. Yan, M. Cantor, C. Namgyal, M. Mino-Kenudson, G. Y. Lauwers, M. Loda, and C. S. Fuchs Sensitive Sequencing Method for KRAS Mutation Detection by Pyrosequencing J. Mol. Diagn., August 1, 2005; 7(3): 413 - 421. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H. Watterson, S. Raha, C. C. Kotoris, C. C. Wust, F. Gharabaghi, S. C. Jantzi, N. K. Haynes, N. H. Gendron, U. J. Krull, A. E. Mackenzie, et al. Rapid detection of single nucleotide polymorphisms associated with spinal muscular atrophy by use of a reusable fibre-optic biosensor Nucleic Acids Res., January 23, 2004; 32(2): e18 - e18. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |