|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular Diagnostics and Genetics |
1 Dipartimento di Patologia Generale, Seconda Università degli Studi di Napoli, Naples; 2 Telethon Institute of Genetics and Medicine (TIGEM), Naples; 3 Dipartimento di Medicina Sperimentale, Servizio di Cardiomiologia e Genetica Medica, Seconda Università degli Studi di Napoli, Naples; 4 Centro di Eccellenza per le Malattie Cardiovascolari, Seconda Università degli Studi di Napoli, Naples, Italy.
aAddress correspondence to this author at: Dipartimento di Patologia Generale, Seconda Università degli Studi di Napoli, S. Andrea delle Dame, via L. De Crecchio 7-80138 Napoli, Italy. Fax 39 0815665704; e-mail vincenzo.nigro{at}unina2.it, nigro{at}tigem.it.
| Abstract |
|---|
|
|
|---|
Methods: We approached the challenge of multiplexing by modifying the PCR chemistry. We set up a rapid protocol that analyzes all dystrophin exons and flanking introns (57.5 kb). We grouped exons according to their effect on the reading frame and ran 2 PCR reactions for DMD mutations and 2 reactions for BMD mutations under the same conditions. The PCR products are evenly spaced logarithmically on the gel (Log-PCR) in an order that reproduces their chromosomal locations. This strategy enables both simultaneous mapping of all the mutation borders and distinguishing between DMD and BMD. As a proof of principle, we reexamined samples from 506 patients who had received a DMD or BMD diagnosis.
Results: We observed gross rearrangements in 428 of the patients (84.6%; 74.5% deletions and 10.1% duplications). We also recognized a much broader spectrum of mutations and identified 14.6% additional cases.
Conclusions: This study is the first exhaustive investigation of this subject and has made possible the development of a cost-effective test for diagnosing a larger proportion of cases. The benefit of this approach may allow more focused efforts for discovering small or deep-intronic mutations among the few remaining undiagnosed cases. The same protocol can be extended to set up Log-PCRs for other high-throughput applications.
| Introduction |
|---|
|
|
|---|
Because DMD is a lethal X-linked disorder, one third of all mutated alleles are removed each generation because they are carried by affected males, who rarely have children. According to Haldanes rule(4), the rate of the appearance of mutated alleles must equal the rate of their removal(5). Thus, the frequency of new mutations has to be very high (up to 10–4), and the mutation spectrum is extremely heterogeneous. Both DMD and BMD are caused by mutations in the dystrophin gene locus [DMD,2 dystrophin (muscular dystrophy, Duchenne and Becker types)](6), which consists of 79 exons and 8 tissue-specific promoters. Dystrophin is not detectable in DMD patients, whereas in BMD patients some dystrophin—albeit abnormal—is visible on a western blot. The dystrophin gene spans 2.22 Mb, encoding the longest known primary transcript of the human genome. The most common mutations are large intragenic deletions or duplications that encompass one or more exons(7). The effects on the transcript depend less on the extent of a deletion (or duplication) than on whether it disrupts the translational reading frame. In general, mutations that maintain the reading frame are associated with BMD, whereas deletions and duplications tend to disrupt the reading frame in DMD(8)(9). A survey of 4700 mutations in the Leiden database (http://www.dmd.nl) indicates that the reading-frame rule holds for 91% of cases(10). Deletion/duplication breakpoints may occur anywhere in the gene. There are, however, large differences in intron size, from 54 bp (intron 14) to 248 342 bp (intron 44). In addition, 20 consecutive exons (exons 23–42) are symmetric (0,0), and each one can be deleted without causing a frameshift. This feature accounts for the locations of 2 major disease-linked hot spots, one within exons 44–52 (710 kb) and the other within exons 2–19 (530 kb). The analysis of a limited number of exons has been thought to detect most of the deletions found in patients; thus, genetic testing for deletions has relied worldwide on different methods, which were originally based on the multiplex PCR technique. This technique initially was set up to analyze 18 DMD exons(11)(12). Other PCR protocols were subsequently developed to cover some untested exons and to define deletion borders(13)(14), and the descriptions of several additional multiplex primer sets have been published over the years. Possibly because of the different sensitivities and efficacies of all the multiplex PCR methods and/or the methods used for clinical diagnosis, deletions have been detected in different countries [Germany(15), Greece(16), Mexico(17), Egypt(18), Morocco(19), Saudi Arabia(20), India(21), and China(22)] at frequencies that range from 50% to 65% of DMD cases. In the majority of laboratories, molecular diagnosis is performed via multiplex PCR alone. Other mutation types, such as atypical deletions, duplications, small mutations, deep-intronic deletions, and the insertion of repetitive sequences, remain undetected.
Additional diagnostic approaches, such as quantitative PCR(23)(24), single-strand conformation polymorphism analysis(25), denaturing HPLC(26)(27), SCAIP ("single condition amplification/internal primer") sequencing(28), multiplex amplifiable probe hybridization (MAPH)(29), or multiplex ligation-dependent probe amplification (MLPA)(30)(31) are required to resolve these cases; however, the latest, more sophisticated techniques require specialized equipment and expertise that are very often unsuitable for a common diagnostic laboratory, especially in developing countries.
We have developed a novel tool, a multiplex PCR assay in which the PCR products are evenly spaced logarithmically on a gel (Log-PCR), that can detect deletions and duplications via analysis of all DMD exons in 4 multiplex PCRs that run under the same conditions.
We validated the method performing the pivotal analysis on a group of 506 DMD/BMD patients.
We found that this rapid, simple and inexpensive tool enabled a definitive molecular diagnosis in 85% of DMD/BMD patients.
| Materials and Methods |
|---|
|
|
|---|
DNA from peripheral blood leukocytes was used in accordance with the standard operating procedures adopted by the EuroBioBank network and was stored at the Naples Human Mutation Gene Bank (Cardiomyology and Medical Genetics) or at the Telethon Institute of Genetics and Medicine.
primer design
Dystrophin amplimers were divided into 2 groups. The first group consists of 40 amplicons of phase 1 or phase 2 asymmetric exons [(0,1), (0,2), (1,0), (1,2), (2,0), (2,1)] that produce a frameshift when absent or duplicated. This amplicon group is composed of 2 PCR sets, A and B. Set A includes exons 1, 3, 5, 7, 11, 17, 19, 21, 43, 45, 50, 52, 54, 56, 58, 61, 63, 66, 68, and 75; set B includes exons 2, 4, 6, 8, 12, 18, 20, 22, 44, 46, 51, 53, 55, 57, 59, 62, 65, 67, 70, and 76.
The second dystrophin amplimer group consists of 38 amplicons of symmetric exons [(0,0), (1,1), (2,2)]. Deletion or duplications involving these exons do not produce a frameshift. This group is also composed of 2 sets, C and D. Set C includes exons 9, 13, 16, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 48, 60, 69, 72, and 74, plus an amplicon for the AMELX gene [amelogenin (amelogenesis imperfecta 1, X-linked)] when required. Set D includes exons 10, 14–15, 23, 25, 27, 29, 31, 33, 37, 39, 41, 47, 49, 64, 71, 73, and 77, plus an amplicon for the AMELY gene (amelogenin, Y-linked) when required. The presence or absence of an amelogenin amplicon helps to distinguish female and male DNA and to check for substantial maternal contamination in samples of DNA extracted from chorionic villi. Consecutive exons are never included in the same set. We excluded exons 78 and 79, which have never been found to be deleted in DMD and BMD(32). For each exon of the DMD gene (NM_004006.1), we based our design of primer pairs so that we (a) produced a predetermined spacing of the PCR products on agarose gels after electrophoresis and (b) retained the chromosomal order of the exons. To establish the length of each fragment, we created constant spacing between electrophoretic bands by using the equation: log Xn = log X(n – 1) + 0.047, where X is amplicon length in bp. We empirically determined that this increment produced the optimal distance with conventional agarose gels. Fragment lengths are alternated between sets to retain the option of running sets A and B together and sets C and D together with higher-resolution electrophoresis methods. Moreover, we designed each primer to fit the following requirements: (a) 28–32 bp in length, (b) at least 13 C or G nucleotides, and (c) at least 3 C or G nucleotides at the 3' end. We used BLASTn (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi) to check primer sequences to avoid matching with repeated human sequences and single-nucleotide polymorphisms. AMELX and AMELY sequences were obtained from GenBank (NM_182680 and NM_001143, respectively). We aligned and compared sequences with ClustalW software to select 2 pairs of primers, the first specific for AMELX alone and the second specific for AMELY. To avoid PCR artifacts due to polymorphic variation, we designed the primer sets to have 4 internal amplification controls: Exon 31 is present both in set D and within the amplimer of exon 32 in set C, exon 11 is present both in set A and within the amplimer of exon 10 in set D, exon 25 is present both in set D and within the amplimer of exon 24 in set C, and exon 34 is contained both in set C and in the amplimer of exon 35 in set D. These redundancies are very useful. For instance, we have found a polymorphic 84-bp deletion at the 5' end of intron 30 (frequency, 0.015) that causes a failure to amplify exon 31, but the internal control helps us to distinguish between this polymorphism and true deletion of exon 31.
multiplex pcr chemistry
Each multiplex PCR was performed in a final volume of 15–20 µL containing 60–100 ng of DNA template, 10x buffer (final concentrations: 15 mmol/L Tris base, 20 mmol/L HEPES free acid, 25 mmol/L KCl, 10 mmol/L MgCl2, pH 8.20), 2x polyol solution (final concentrations: 150 mmol/L maltitol, 350 mmol/L sorbitol), a primer mixture (total final concentration of primers, 4 µmol/L; each primer, 0.1 µmol/L), deoxynucleoside triphosphates (final concentration of each, 0.25 mmol/L), and 1.25 U of AmpliTaq DNA polymerase (Applied Biosystems).
The cycling conditions of all 4 multiplex PCRs were identical, so we were able to perform all 4 reactions simultaneously with the same thermal cycler. The multiplex PCRs were carried out with a hot-start method on an MWG AG Biotech thermocycler. The cycling conditions consisted of a first step at 99 °C for 30 s, a pause at 85 °C (for Taq addition), and 21 cycles of 20 s at 97 °C and 7 min plus 20 s/cycle at 66 °C (total time, 4 h 30 min). Two microliters of each sample were run on a DNase- and RNase-free agarose gel (16 g/L SeaKem LE; Lonza), containing GelStar stain (10 000x; Cambrex) and 0.5x TTE buffer (final concentrations: 45 mmol/L Tris, 15 mmol/L taurine, and 0.3 mmol/L EDTA, pH 8.0). By using GelStar as the DNA-staining dye (10-fold more sensitive than ethidium bromide), we were able to reduce the number of PCR cycles to 21, thereby maintaining the PCR reaction within the logarithmic phase. Gel electrophoresis was performed for 50 min at 160 V (Fig. 1
).
|
| Results |
|---|
|
|
|---|
|
|
We first developed reaction and cycling conditions for set A by testing various multiplex PCR protocols with different deoxynucleoside triphosphate, primer, and Taq concentrations and at various annealing/extension temperatures and times. We used 0.1 µmol/L of each primer (4 µmol/L total) and 60 ng of DNA template. Because we were initially unable to amplify all of the fragments homogeneously with one protocol, we investigated various modifications to the PCR chemistry. We found that all products were amplifiable with a buffer consisting of 15 mmol/L Tris base, 20 mmol/L HEPES free acid, and 25 mmol/L KCl. One or more products were not amplified with different salts (e.g., ammonium sulfate). We also investigated the addition of polyols or sugar alcohols to improve band homogeneity and consistency by testing different concentrations of 5 compounds, either alone or in pairs: D-sorbitol, D-mannitol, xylitol, maltitol, and meso-erythritol. A mixture of maltitol and sorbitol was the most effective combination. We found the optimal MgCl2 concentration to be 10 mmol/L, a concentration that is at least 6-fold higher than typically used in PCR reactions. We then successfully applied these conditions to the other 3 multiplex PCRs containing the B, C, and D primer sets.
We validated the Log-PCR method by retesting samples from a large cohort of DMD and BMD patients (n = 506). The patients were divided into 2 groups The first group consisted of 152 patients who had no apparent deletions according to the tests of Chamberlain and Gibbs(11) and Beggs et al.(12). Our Log-PCR method revealed that 74 (48.7%) of the DMD and BMD patients in this group had mutations (23 deletions and 51 duplications) that were previously undetected. The second group consisted of 354 patients who had documented deletions. We verified the deletion endpoints in all of these patients and found larger deletions in 170 (48%) of the patients in this group. In brief, Log-PCR identified deletions in 74.5% of cases and duplications in 10.1% (see the figures in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol54/issue6). Our approach improved the detection rate of deletions and duplications by 14.6%. We identified 377 deletions with 127 different borders and 51 duplications with 43 different borders. Thus, each mutation is unique in almost all cases. We tested all deletions and duplications with the "reading frame checker"(32) and identified 21 exceptions to the reading-frame rule (see supplementary data in the online Data Supplement).
deletions
We detected 66 deletions of single exons, 9 of which are undetected with standard tests. All single-exon deletions were confirmed via a single PCR with different primer pairs. Multiple deletions spanned 2 to 43 exons. We documented 41 deletions (11%) in the 5' hot spot (exons 2–19), and 242 deletions (64%) in the 3' hot spot (exons 43–52).
We refined the locations of the deletion endpoints in 44% of the cases with deletions in the 5' hot spot and in 34% of the deletions in the 3' hot spot. The largest deletions spanned exons 17–59 (1100 kb), exons 19–52 (820 kb), exons 10–44 (670 kb), exons 20–47 (615 kb), exons 20–45 (550 kb), and exons 8–42 (400 kb). We found 22 deletions outside the hot spots (see Fig. 4 in the online Data Supplement).
duplications
The Log-PCR also detects duplications when a small number of cycles is used. In the patient group with no previously documented deletions (n = 152), we detected 51 duplications, all of which were confirmed by real-time PCR analysis. Nine of these duplications involved a single exon. In the 5' hot spot, the frequencies of deletions and duplications were 11% and 43%, respectively. We found the reverse situation in the 3' hot spot, where the frequencies of deletions and duplications were 66% and <20%, respectively, suggesting different mutation susceptibilities in the 2 regions. The maximum duplication size was 1 Mb, spanning exons 33–60 (see Fig. 5 in the online Data Supplement).
negative results
We found large deletions and duplications to be absent in 78 patients (15.4%). Such a finding is usually explained by a failure to detect small or deep-intronic mutations. We have analyzed all of the exons and flanking intronic regions in these patients via high-throughput denaturing HPLC followed by direct sequencing, and we have identified a causative point mutation in 56 of these patients (unpublished data). This result improved the rate of detection of causative mutations to 95.7% (see Fig. 6 in the online Data Supplement), which is very close to the rate obtained by Zeng et al.(35). In 4.3% of the cases, no mutations in the dystrophin gene could be identified after DNA analysis. We cannot exclude the presence of deep-intronic or atypical mutations. Alternatively, milder phenotypes may have been misdiagnosed as BMD, a situation that might occur especially when muscular dystrophy is diagnosed at a young age.
| Discussion |
|---|
|
|
|---|
An appreciable increase in the detection rate of deletions appears possible only if PCR tests are coupled with independent methods of analysis(33)(34)(35).
MAPH(29) and MLPA(30)(31) are 2 new techniques that have been developed to detect deletions and are able to screen many target sequences simultaneously. MLPA is also commercially available (MRC-Holland) as 8 separate reactions for detection on agarose gels. This method relies on hybridization of a sequence-specific probe to genomic DNA, subsequent amplification of the hybridized probes with a common primer set, and analysis of the resulting PCR products. As with conventional PCR, polymorphisms or single-base mutations in the probe-binding regions may affect results. The short length of the region identified by the specific probe (21 nucleotides) implies that mismatches at the binding site may prevent probe hybridization and hence prevent ligation and, ultimately, detection. Consequently, single-base changes may appear as exon deletions. In cases of exon rearrangements, such as translocations or the insertion of long interspersed elements, MLPA can miss mutations. Moreover, some true deletions randomly remain undetected with MLPA(33). Thus, MLPA results should be confirmed with an independent method.
We propose Log-PCR as a new tool for complete screening of dystrophin exons and for sex testing. This method uses only 4 quantitative multiplex PCRs, which are run under the same reaction and cycling conditions. Protocols for quantitative PCR(13) and simultaneously amplifying all of the dystrophin exons(28) have already been described. The novelty of the Log-PCR approach is the achievement of both results in a single-step assay. It is applicable to genomic DNA extracted from blood or chorionic villi.
Multiplex PCRs are generally characterized by low product yield, product dropout, and nonspecific amplification. We carried out a long series of experiments to develop a robust quantitative multiplex PCR method that produces yields that are comparable for all amplicons and proportional to the initial amount of DNA template.
The idea of testing polyols and cosolvents originated from our review of early studies that found that the presence of sugars and polyols increased the thermal stability of proteins(36). A later study demonstrated the use of trehalose for the thermostabilization of reverse transcriptase(37). Including such compounds is crucial for the simultaneous and quantitative amplification of 57.5 kb of DNA fragments. The Log-PCR procedure substantially shortens reaction setup times and reduces reagent consumption, and these features support its use as a simple but universal test. We designed primers of 28–32 nucleotides in length with GC-rich 3' ends that are never affected by latent polymorphic nucleotides. In addition, we used denaturing HPLC with 600 control samples to scan for rare variants in the annealing region. Thus, the Log-PCR was able to minimize the occurrence of false single-exon deletions, and we confirmed 66 of 66 single-exon deletions. Moreover, because most of the primers anneal to sequences far from exonic sequences, flanking intronic deletions can also be detected. For example, the Log-PCR method discovered a case of exon 22 deletion, which was not detectable with exonic primers because the missing nucleotides were in the flanking intron sequences and produced exon skipping.
To verify the reliability of the technique, we used the Log-PCR method with samples from a group of 506 patients with DMD or BMD that had been diagnosed at a single center according to homogeneous clinical criteria. We excluded samples from other sources to avoid introducing any bias. To our knowledge, this study is the most extensive single-center study of this kind ever done. Thus, the statistics for the dystrophin mutations we have described (506 cases) can be considered a very good approximation to the real frequency distribution of mutations. In fact, these data are not biased by the detection method because we used a homogeneous procedure to study all of the exons of every patient. In contrast, current databases include data produced with a heterogeneous group of methods and incomplete analyses. Such data sets are likely to overrepresent easy-to-discover mutations.
We hypothesize that geographic and ethnic differences are not important in X-linked lethal disorders, because all of the mutations have a short life in the pedigrees; however, because our data are derived from a Caucasian population, this distribution should be checked in other populations.
clinical use of the log-pcr
Complete visualization of all dystrophin exons allows mutation boundaries to be precisely defined. Such a comprehensive approach may be important when diagnostic studies are performed at a young age, when distinguishing between DMD and BMD is difficult. The identification of mutation endpoints may be important for future patient therapy via antisense oligonucleotide exon skipping(38).
Log-PCR is a noninvasive, sensitive, and specific laboratory test for the diagnosis of DMD and BMD. The assay requires approximately 6 h and is cost-effective. It produces a report that is direct and easy to interpret, because all of the bands are ordered and evenly spaced. The assay requires the development of a reagent set that can be used in routine diagnostic laboratories. Compared with the MLPA reagents from MRC-Holland, the reagents should be 5 times less expensive, but labor requirements also should be addressed. The Log-PCR approach can also be used on microfluidics-based platforms (such as the Agilent 2100 Bioanalyzer) that reproducibly and completely quantify each fragment.
Log-PCR can also be useful for assessing carrier status (data not shown) when the mutation is known. We detected deletions in all of the carriers we tested. We believe that the development of the Log-PCR method also provides a proof of principle for higher-throughput multiplex PCR methods. In the present study, we quantitatively amplified 57 kb of DNA sequence, and we expect this approach to be of use in resequencing large genomic regions (enrichment step). The method is therefore directly applicable to such resequencing without any further development.
Results obtained for a preliminary test series of samples sequenced on an ABI 3130XL DNA Sequencer (Applied Biosystems) have shown that direct sequencing of Log-PCR products is possible with common purification procedures and specific primers for sequencing.
| Acknowledgments |
|---|
Financial Disclosures: None declared.
Acknowledgments: We thank Alessandra Ferlini for helpful suggestions and Maria Grazia Esposito at Cardiomyology and Medical Genetics for DNA processing. We acknowledge the SUN-Naples Human Mutation Gene Bank (Cardiomyology and Medical Genetics), which is a partner of the EuroBioBank network.
| Footnotes |
|---|
2 Human genes: DMD, dystrophin (muscular dystrophy, Duchenne and Becker types); AMELX, amelogenin (amelogenesis imperfecta 1, X-linked); AMELY, amelogenin, Y-linked. ![]()
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
A. Torella, A. Trimarco, F. Del Vecchio Blanco, A. Cuomo, S. Aurino, G. Piluso, C. Minetti, L. Politano, and V. Nigro One Hundred Twenty-One Dystrophin Point Mutations Detected from Stored DNA Samples by Combinatorial Denaturing High-Performance Liquid Chromatography J. Mol. Diagn., January 1, 2010; 12(1): 65 - 73. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |