|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular Diagnostics and Genetics |
1 Translational Research Laboratory, IDIBELL-Institut Català dOncologia, Barcelona, Spain; 2 Bioinformatics and Biostatistics Unit, Department of Epidemiology, IDIBELL-Institut Català dOncologia, Barcelona, Spain; 3 Human Genetics Group, Spanish National Cancer Centre (CNIO), and Centre for Biomedical Research on Rare Diseases, Instituto de Salud Carlos III, CIBERER; 4 Genetic Counselling Unit, IDIBELL-Institut Català dOncologia, Barcelona, Spain; 5 Department of Genetics, INSERM U614, Faculty of Medicine, Rouen, France.
aAddress correspondence to this author at: Laboratori de Recerca Translacional (LRT1), Institut Català dOncologia, Hospital Duran i Reynals, Gran via s/n Km 2.7, 08907 LHospitalet de Llobregat, Spain. Fax +34 93-260-74-66; e-mail gcapella{at}iconcologia.net.
| Abstract |
|---|
|
|
|---|
Methods: We set up and validated APC-gene QMPSF using 23 negative and 1 positive control and examined 45 (13 FAP and 32 AFAP) unrelated members of APC/MUTYH mutation-negative families for copy number alterations. We confirmed the results using multiplex ligation-dependent probe amplification (MLPA). We used different approaches such as sequencing, quantitative real time-PCR (QRT-PCR), and fluorescence in situ hybridization (FISH) to further characterize the identified deletions.
Results: APC QMPSF was capable of detecting deletions with an acceptable variability, as shown by mean values (SD) of allele dosage for the deleted control obtained from intra- and interexperimental replicates [0.52 (0.05) and 0.45 (0.10)]. We detected 3 gross deletions in 13 (23%) of the classic FAP cases analyzed (1 complete gene deletion and 2 partial deletions encompassing exons 9 and 10 and exons 11–15, respectively). No rearrangements were detected in the 32 AFAP cases.
Conclusions: QMPSF is able to detect rearrangements of the APC gene. Our findings highlight the importance of using a copy number alteration methodology as a first step in the routine genetic testing of FAP families in the clinical setting.
| Introduction |
|---|
|
|
|---|
Cytogenetically detectable interstitial deletions at 5q22 have been reported in FAP patients exhibiting some degree of mental retardation and dysmorphism (6). Submicroscopic deletions have also been detected in 6% to 17% of FAP families (6) by Southern blot analysis, apparent nonparental segregation of intragenic or flanking polymorphic marker alleles, fluorescence in situ hybridization (FISH) analysis, quantitative real-time PCR (QRT-PCR), comparative genomic hybridization, RNA analysis, and multiplex ligation-dependent probe amplification (MLPA) (6)(7)(8)(9)(10)(11).
A simple, semiquantitative procedure based on quantitative multiplex PCR of short fluorescent fragments (QMPSF) has been shown to be useful for detecting genomic deletions or duplications of mismatch repair and BRCA genes (12)(13)(14). We adapted this technique for the study of copy number alterations in the APC gene in a series of unrelated FAP patients with no previously detected point mutation of the APC or MUTYH genes. We compared the results to those obtained with other techniques for determining the relative copy number of DNA sequences such as MLPA (15) and characterized the molecular nature of the alterations.
| Materials and Methods |
|---|
|
|
|---|
apc qmpsf
Between 8 and 10 exons can be amplified together in quantitative multiplex PCR assays (14)(16). The APC gene (
100 kb) is composed of 14 small exons and exon 15, which covers more than three quarters of the coding sequence. Thus, 18 amplicons of between 130 and 330 bp were designed to cover the entire 8.5-kb coding region and distributed across 3 multiplex reactions (Table 1
). As a negative control, we included in each reaction a fragment of a different gene in which deletion was not expected (HMBS, hydroxymethylbilane synthase at 11q23.3). Primers, designed using Primer Premier 5 (Premier Biosoft International), shared common melting temperatures of between 58 and 60 °C. At least 1 of the primers was located in the exonic region, and no primer covered polymorphic regions. All primers were chimeric and carried 10-nt extensions on their 5' ends to homogenize annealing temperatures and relative primer concentrations (available on request). One primer from each pair was 5'-labeled with 6-FAM fluorochrome (Table 1
). We set the distribution of the primers across the 3 multiplex reactions using WinMultiPLX v1.00 (BioData Ltd.). We amplified 100 ng genomic DNA in a final volume of 25 µL including 0.4–1 µmol/L primers, 200 µmol/L dNTPs, 1.5 mmol/L MgCl2, 1 U Taq polymerase (Thermoprime Plus DNA Polymerase; ABgene), 1 µL DMSO, and 1.25 µL triethylammonium acetate. After an initial denaturing step, samples underwent 23 cycles (10 s at 94 °C, 15 s at 50 °C, 20 s at 72 °C), which ensured that the amplification ended during the exponential phase. After the multiplex reactions, the DNA fragments were separated on an ABI Prism 3130 DNA sequencer and analyzed using GeneScan Analysis Software v3.1.2 (Applied Biosystems).
|
Three different methods were used to calculate allele dosage in this study: visual sample-to-control comparison, numerical sample-to-control comparison, and numerical cumulative comparison. For the visual sample-to-control comparison, we estimated allele dosage by superimposing the electropherogram of the tested sample onto the corresponding image for a control DNA sample after adjusting the vertical scale of the different amplicons. Allelic losses of one or more amplicons are represented by a 2-fold reduction in the intensity (peak heights) of an amplicon of the sample analyzed. In the numerical sample-to-control comparison, we calculated numerical normalized ratios (R) using the formula: R = (peak intensityamplicon x sample/peak intensityamplicon x control)/(peak intensityamplicon ref sample/peak intensityamplicon ref control). An R value close to 0.5 (0.4–0.6) represented a 2-fold reduction. For numerical cumulative comparison, we adapted the analytical method used in MLPA (www.mrc-holland.com) to QMPSF. After normalization of all peak intensities to the intensity of the reference (HMBS) peak, we obtained the numerical normalized ratio (R') of a given exon and sample using the formula R' = normalized peak intensity of analyzed sample/mean of normalized peak intensities of all samples. As above, R' values close to 1 and 0.5 indicate whole and half allele dosage, respectively.
mlpa analysis
MLPA analysis, carried out using an APC-specific MLPA test (Salsa MLPA Kit P043 APC, MRC-Holland) according to the manufacturers instructions, used 300 ng genomic DNA. After hybridization, ligation, and amplification, we calculated the PCR product ratios by dividing each measured peak area by the sum of all peak areas for the sample. The ratio of each peaks individual relative probe area was normalized to the mean obtained with 2 control samples. We analyzed the products using the ABI 3700 DNA sequencer and Peak Scanner Software v1.0 (Applied Biosystems).
molecular characterization of deletions
We performed PCR to characterize the deletion of exons 9 and 10 (forward flanking primer 5'-ACCTATAGTCTAAATTATACCATC-3'; reverse 5'-CTGAGCTATCTTAAGAAATACATG-3') in a final volume of 50 µL containing 0.2 µmol/L each primer, 1 U Taq polymerase (Thermoprime Plus DNA Polymerase; ABgene), 200 µmol/L dNTPs, and 1.5 mmol/L MgCl2. After denaturation at 94 °C for 10 min, 30 cycles of PCR were performed (94 °C for 30 s; 55 °C for 30 s; 72 °C for 5 min). After using the PCR Product Purification Spin Kit/250 (Genomed), 1 µL of the amplified product was sequenced using the BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) according to the manufacturers instructions. The primers used for sequencing the breakpoint were 5'-TTTCCGGTTTCTAATCTGCTTCT-3' as the forward primer and 5'-CCTTTTTATGCCAGTTTCTCCAC-3' as the reverse primer.
We performed QRT-PCR with the Roche LightCycler® 2.0 System according to the SYBR Green Detection methodology. We used 8 sets of primer pairs, including a pair that amplified a fragment used as a reference for the presence of the 2 alleles (APC_e5) and a pair that amplified a deleted fragment (APC_e15) (see Supplemental Table 1 in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol54/issue7 ). Reactions were performed in a volume of 10 µL containing 1 µL LightCycler Fast Start DNA Master SYBR Green I (Roche), 1.5 mmol/L MgCl2, 0.5 µmol/L of each primer, and 50 ng genomic DNA for 40 cycles (95 °C for 6 s; melting temperature for 16 s; 72 °C for 20 s). The samples were then subjected to melting curve analysis from 60 to 95 °C. Crossing points (the cycle number when the amplification becomes exponential) were determined using LightCycler software. We calculated the gene copy number using the relative standard curve method (17).
We performed FISH on chromosome spreads obtained from peripheral blood of the relevant patients. Cultures were set up a few hours after blood collection in heparinized tubes. We added 500 µL whole blood to a flask containing 9 mL RPMI 1640 (Gibco) with 10% FCS, 5 mg/L phytohemagglutinin (Gibco), 100 000 U/L penicillin, and 100 mg/L streptomycin. Cultures were set up in duplicate and incubated at 37 °C for 72 h. We harvested cells after colcemid treatment (0.4 mg/L) for 20 min and then performed FISH. Briefly, we used 7 different probes from chromosome 5q22 to cytogenetically clone the deletion breakpoints RP11-809J10, RP11-619D6, RP11-616C24, RP11-389E12, and RP11-145D17 (ordered from the centromere to the telomere and covering approximately 1 Mb). All probes were obtained from the 32K BAC Rearray of the Childrens Hospital Oakland Research Institute (www.chori.org/bccpcc) and kindly provided by Luis Pérez-Jurado. Probes (1 µg of each) were directly labeled by a standard nick translation reaction, hybridized in pairs, and detected according to standard protocols (18). The FISH analysis was performed under an Olympus BX60 fluorescence microscope equipped with the appropriate filter set, and images were captured using the Cytovision software (Applied Imaging).
| Results |
|---|
|
|
|---|
|
We assayed 10 replicates of both a negative control and the deleted control in the same experiment to assess the intraexperimental variability (Supplemental Data Fig. 2A). Using the cumulative comparison, mean R' values of 1 (0.09; 0.73–1.25) and 0.52 (0.05; 0.38–0.71) were obtained for the normal and the deleted controls, respectively. As in the previous assay, the variability was concentrated in a few amplicons, and more extreme values appeared in only 1 or 2 amplicons of a replicate, thus reducing the likelihood of interpretation errors. Finally, interexperimental variability was studied in 13 independent QMPSF assays using the deleted control. Numerical cumulative comparison with the 23 HNPCC patients yielded R' values of 0.45 (0.10; 0.14–0.69) (see Supplemental Fig. 2B in the online Data Supplement). There was no overlap between any of the wt and deleted control ratios.
Thus, QMPSF was able to detect the loss of 1 copy of the gene for all fragments analyzed, with an acceptable amount of intra- and interexperimental variability.
deletion screening in patients
We used the above data to establish the following experimental design incorporating 2 duplicates of 2 wt controls in each run and a single sample of the deleted control. The run was disregarded when the ratio between replicates of wt controls was >1.2 or <0.8. Once a copy number variation was observed, the experiments were independently replicated at least twice using independent DNA aliquots. Deletions were detected in 3 of the 13 (23%) selected APC/MUTYH mutation-negative classic FAP families (see Supplemental Table 2 and Supplemental Fig. 3 in the online Data Supplement). Two of these deletions were partial, and the other encompassed the entire gene. The partial deletions encompassed exons 9 and 10 and from exon 11 to all amplicons of exon 15, respectively. The MLPA analyses performed for the 13 classic FAP families confirmed the QMPSF results. No gross rearrangements were detected in the 32 AFAP mutation-negative families.
Regarding genotype-phenotype correlation in classic FAP families, carriers of gross deletions (n = 4, including the cytogenetically deleted control) were diagnosed at an earlier age [23.25 (7.09) years] than APC/MUTYH mutation-negative index cases [36.9 (7.84) years; n = 10; P = 0.028, Wilcoxon rank-sum test] (see Supplemental Table 2 and Supplemental Fig. 3 in the online Data Supplement). Differences were detected with respect to other APC gene mutation carriers [29.28 (13.39) years; n = 57], although the differences were not significant, possibly owing to the low number of deletion carriers. In contrast, no differences were observed in the number of colorectal polyps, the presence of extracolonic manifestations, or family history.
molecular characterization of the partial apc deletions
We confirmed all 3 of the observed deletions by replicating QMPSF assays and using MLPA. Further characterization was carried out in partial deletion carriers, which provided insight into the molecular nature of these aberrations.
In an index case, QMPSF analysis showed a deletion removing exons 9 and 10 (Fig. 2
, A and B). Genomic DNA PCR encompassing exons 8–11 revealed the amplification of an abnormal fragment of approximately 4 kb in the affected patient, corresponding to the junction fragment comprising the deleted region (Fig. 2
, A and C). Sequence analysis revealed that the first 3007 bp of intron 8 were fused to the last 496 bp of intron 10, deleting a 8011-bp fragment of the gene (g.80715_88725del8011; NC_000005.8; GI:51511721) (Fig. 2D
). Alignment of intronic sequences using the RepeatMasker program revealed that the recombination did not involve any flanking repeat elements, suggesting a nonhomologous recombination process (Fig. 2D
).
|
Fig. 3
shows a schematic representation of molecular studies done in an index case carrying a deletion of exon 11 to all fragments of exon 15. A QMPSF pattern showing 2-fold reduction of intensities of these exons was detected in the index case (Fig. 4A
) and in another affected relative (data not shown). To further characterize the 3' extension of the deletion, genomic DNA QRT-PCR was performed on eight fragments up to 400 kb downstream of the deleted area using exon 5 of the APC gene as a reference for normalization (Fig. 3
, Fig. 4A
; see Supplemental Table 1 in the online Data Supplement). Samples from an affected member and a healthy relative were compared with control DNA. The 3' breakpoint was located within a region of approximately 83 kb encompassed by fragments REEP5_e3 and DCP2_e2 in the affected subject (Fig. 3
, Fig. 4B
). FISH analysis using 5 probes covering 835 kb also confirmed the approximate location of the breakpoint (Fig. 3
, Fig. 4C
). Probes corresponding to RP11-809J10, RP11-619D06, RP11-389E12, and RP11-145D17 clones yielded a double signal concordant with the QMPSF and QRT-PCR results. Although approximately half of the area covered by RP11-619D06 was deleted, the nondeleted area was sufficient to hybridize and identify the corresponding fluorescent signal. Only the RP11-616C24 probe yielded a single signal. No attempt was made to sequence the 5' breakpoint.
|
|
| Discussion |
|---|
|
|
|---|
Copy number variation is known to be an important factor in determining genetic predisposition to disease, in particular cancer. Germline deletions in cancer genes have been associated with a number of hereditary cancer syndromes including neurofibromatosis, HNPCC, hereditary breast cancer, and FAP. Duplications have been detected, albeit at a lower frequency (19). Various techniques have simplified this type of analysis: MLPA is one widely used method (15), real-time PCR methods have also been successfully developed (20), and QMPSF is a simple semiquantitative method that had already been used in the analysis of germline and somatic copy number variations (12)(13)(14)(21). Whereas primer design is a common critical step in all of these methods, the quality of DNA may also influence the final outcome. To minimize the putative impact of this source of variability, the entire procedure for sample collection and processing must be highly standardized.
In this study, APC QMPSF consistently detected the loss of 1 copy of all gene fragments. In this approach, it is not necessary to use the large number of control probes included in each MLPA reaction. On the other hand, while MLPA is performed in a single PCR reaction, QMPSF requires 3 reactions. The even distribution of exons across the 3 multiplex reactions diminishes the possibility of artifacts. In contrast with the ready-to-use MLPA kit, our in-house QMPSF approach may be more prone to variability when used in other centers. Finally, the flexibility of QMPSF allows the exploration of any genomic region of interest without the need for a specific commercial kit.
The method used to analyze allele dosage is also an important factor. We prefer the numerical cumulative approach, as it provides more accurate calculations when multiple negative samples are considered. This is often the case, since most of the tested samples will be negative. A strategy for maintaining accuracy when analyzing a smaller number of samples is to include negative controls. Alternatively, the numerical sample-to-control comparison can be used. In the case of MLPA, the latter approach is not routinely used, so it is necessary to analyze a relatively large number of samples in each assay. Although visual inspection can be more intuitive, it should always be used in combination with the sample-to-control numerical calculation.
Robustness is vital if the results of the assay will be used in the routine clinical setting. All deletions detected in the present study were confirmed by performing independent replicates of the QMPSF analysis and by MLPA, which validated the reliability of our results. Although molecular characterization of deletions is strictly necessary only when only 1 or a few exons are altered (22)(23), its importance should not be underestimated in the case of aberrations involving larger numbers of exons. Therefore, we attempted to further characterize the 2 partial deletions found. We used QRT-PCR using a relative quantification approach (relative standard curve method) and FISH analyses to obtain a partial characterization of the deletion encompassing exon 11 to the 3' end of the gene, in which both boundaries were completely unknown. Finally, the exon 9 and 10 deletion was detected by QMPSF and carefully characterized using PCR and breakpoint sequencing. We believe that the diagnosis of gross rearrangements should not be made unless it has been confirmed by more than one analytical method.
Gross deletions were detected in 23% (3 of 13) of the APC/MUTYH mutation-negative classic FAP families analyzed. When all classic FAP families were considered, gross deletions were found in approximately 5% (4 of 75), including the cytogenetically detected deletion used as a control. This prevalence is consistent with previous reports (6)(8)(9) and corroborates that comprehensive mutation scanning of APC must include gene dosage analysis in mutation-negative classic FAP families. QMPSF is a straightforward, rapid, and low-cost screening technique—similar to MLPA—that should be used before mutational analysis of coding sequences, which is consistent with the conclusions of previous reports (24). Although other more advanced methods such as array-based comparative genomic hybridization may eventually provide better resolution, their current cost makes them unaffordable as a diagnostic screening tool. As in most previous research, no gross rearrangements were detected in the 30 AFAP families. Nonetheless, deletions have occasionally been associated with an attenuated phenotype (25)(26)(27), making it unwise to rule out the use of deletion screening in this setting. In the limited set of cases analyzed in this study, deletion carriers had an earlier age at diagnosis, although this association was not found in a recent review (25). As reported, extracolonic features commonly occur in these cases.
In summary, QMPSF is a robust method for APC copy number alteration analysis either as the principal screening method or as a complementary tool for confirming MLPA results. Our results also confirm the relevance of gross germline deletions in APC/MUTYH mutation-negative classic FAP families and support the use of QMPSF in the diagnostic algorithm for FAP patients. Importantly, this technique is not restricted to germline FAP/AFAP cases. It has been shown recently that the use of QMPSF to assess somatic copy number variations may be of prognostic value in colorectal cancer (21).
| Acknowledgments |
|---|
Financial Disclosures: None declared.
Acknowledgments: We thank Jesús del Valle, Conxi Lázaro, Marta Pineda, and Mireia Menéndez for critical assessment and helpful comments on the manuscript, and Olaya Villa for providing the BAC clones.
| Footnotes |
|---|
2 Human genes: APC, adenomatosis polyposis coli; MUTYH (or MYH), mutY homolog (E. coli); BRCA, breast cancer; MSH2, mutS homolog 2, colon cancer, nonpolyposis type 1 (E. coli); MLH1, mutL homolog 1, colon cancer, nonpolyposis type 2 (E. coli); HMBS, hydroxymethylbilane synthase (formerly PBGD); REEP5, receptor accessory protein 5 (formerly C5orf18); DCP2, DCP2 decapping enzyme homolog (S. cerevisiae). ![]()
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |