- Short Report
- Open Access
Genetic screening of Fabry patients with EcoTILLING and HRM technology
BMC Research Notesvolume 4, Article number: 323 (2011)
Anderson-Fabry disease (FD) is caused by a deficit of the α-galactosidase A enzyme which leads to the accumulation of complex sphingolipids, especially globotriaosylceramide (Gb3), in all the cells of the body, causing the onset of a multi-systemic disease with poor prognosis in adulthood. In this article, we describe two alternative methods for screening the GLA gene which codes for the α-galactosidase A enzyme in subjects with probable FD in order to test analysis strategies which include or rely on initial pre-screening.
We analyzed 740 samples using EcoTILLING, comparing two mismatch-specific endonucleases, CEL I and ENDO-1, while conducting a parallel screening of the same samples using HRM (High Resolution Melting). Afterwards, all samples were subjected to direct sequencing. Overall, we identified 12 different genetic variations: -10C>T, -12G>A, -30G>A, IVS2-76_80del5, D165H, C172Y, IVS4+16A>G, IVS4 +68 A>G, c.718_719delAA, D313Y, IVS6-22C>T, G395A. This was consistent with the high genetic heterogeneity found in FD patients and carriers. All of the mutations were detected by HRM, whereas 17% of the mutations were not found by EcoTILLING. The results obtained by EcoTILLING comparing the CEL I and ENDO-1 endonucleases were perfectly overlapping.
On the basis of its simplicity, flexibility, repeatability, and sensitivity, we believe that HRM analysis of the GLA gene is a reliable presequencing screening tool. This method can be applied to any genomic feature to identify known and unknown genetic alterations, and it is ideal for conducting screening and population studies.
Fabry disease (FD) is a lysosomal storage disease caused by a congenital error in glycosphingolipid metabolism resulting from the deficient or absent activity of the lysosomal enzyme α-galactosidase A . This enzyme is essential for the sphingolipid recycling process that occurs within cells. It is able to remove the galactose bound in the alpha position in complex sphingolipids, allowing their metabolism. The enzyme deficiency causes the interruption of the process of waste product demolition, leading to the intracellular accumulation of complex sphingolipids, especially globotriaosylceramide (Gb3) . The accumulation of these products is the basis of FD symptomology that includes acroparesthesias, angiokeratoma, corneal and lenticular opacities, gastrointestinal problems and diseases of the kidney, heart and central nervous system [3–5]. The GLA gene (12 kb) contains seven exons (92 to 291 bp), is located on the X chromosome, and encodes a polypeptide of 429 amino acids .
In the literature, more than 500 mutations in the GLA gene have been associated with FD . Of these, only 5% are localized in exons and are able to cause the formation of a truncated protein. Instead, for the most part, they involve missense mutations generating null alleles or changes in specific domains of the protein; others, such as those that occur in splicing sites, can lead to aberrant splicing and can be responsible for altered forms of the enzyme [8, 9]. Polymorphisms, small deletions, and insertions in the GLA gene have been detected by direct sequencing , through single-strand conformation polymorphism analysis , and by denaturing high performance liquid chromatography (DHPLC) . However, these investigative techniques may be limited in the analysis of large numbers of samples, so it is therefore necessary to move towards new analytical approaches which include or rely on pre-screening.
In our study, we performed genetic analyses of the GLA gene in 740 subjects with symptoms related to FD using two different approaches to genetic screening. The first was a method based on the enzymatic cleavage of mismatches, known by the acronym EcoTILLING (Targeting Induced Local Lesions in Genomes), originally used for studying the plant genome and adapted to researching human genetic variations . We conducted these analyses using a screening strategy based on pools of 8 samples using two mismatch-specific endonucleases, CEL I (SURVEYOR Mutation Detection Kit-Transgenomic) and ENDO-1 (Serialgenetic). The second method was based on HRM which is able to detect genetic mutations based on the dissociation temperature of the two strands (Tm) and the analysis of the melting curves (LightCycler480, Roche). We applied these techniques to analyze the coding regions and intron junction regions of the GLA gene in order to investigate mutations in the coding regions and splicing regulatory sites. Finally, to verify the results obtained, each sample was analyzed by direct sequencing.
Materials and methods
DNA extraction and quantification
DNA was extracted by column extraction (GenElute Blood Genomic DNA Kit, Miniprep - Sigma Aldrich) from 740 peripheral blood samples of subjects with clinical manifestations related to FD symptomatology. The study was approved by University Hospital Ethics Committee and written informed consent was obtained from all participants.
The concentration of the DNA was determined by spectrophotometer (Victor X3 - PerkinElmer), and the concentration was adjusted to 100 ng/microliter in TE (10 mM Tris-HCl and 1 mM EDTA, pH 7.4) or water.
Seven pairs of primers were designed for the analyses of the seven target regions containing the seven exons of the GLA gene and the splicing regulatory regions flanking them, and they were obtained using the GeneTools program. The gene sequence was obtained from the site http://www.ensembl.org (ENSG00000102393). The primers used for the EcoTILLING and HRM analyses were identical but for marking with IRDye700 in forward and IRDye800 in reverse (Table 1).
For the simultaneous analysis of a large number of samples, we chose a strategy of pooling, as described by Till et al. . For a high-throughput screening, the 740 DNA samples were divided into groups of 64 samples each, arranged in 8 × 8 grids and then assembled in pools of 8 samples (Figure 1). Each pool was obtained by adding all the samples present in a row or column of the grid, repeating the procedure for the 8 rows and 8 columns of each grid. The 16 pools derived from each grid were arranged in a second plate and subjected to the EcoTILLING method. The presence of the fragments in two pools containing the same sample allowed us to trace the exact location of the sample in the grid and then to identify the subject.
Each PCR was performed in a final volume of 50 μl (Hot Start Master Mix - GE Healthcare) with 1 μM of primers, 50 ng of DNA. The amplification was performed using the following parameters: 95°C for 2 minutes, 40 cycles with 95°C for 30 minutes, 58°C for 30 minutes, 72°C for 1 minute, and a final extension at 72°C for 10 minutes.
Mixing 8 different samples ensures the formation of a heteroduplex in the presence of a mutation. The samples thus obtained were denatured and renaturated according to the following parameters: 95°C for 2 minutes, followed by a decrease in temperature of 0.5°C for 15 minutes, until 45°C. The heteroduplexes were subjected to enzymatic digestion by the endonucleases CEL I (SURVEYOR Mutation Detection Kit - Transgenomic) and, in parallel, ENDO-1 (Serialgenetics). The reaction was incubated at 42°C for 20 minutes and stopped by adding 1.5 μl stop buffer, or at 37°C for 40 minutes and stopped with 5 μl stop buffer, respectively. All reagents were supplied by the Surveyor kit (SURVEYOR Mutation Detection Kit - Transgenomic) and by the ENDO-1 Endonuclease Kit (Serialgenetics).
The fragments were separated by electrophoretic migration (DNAanalyzer 4300 - Licor) on polyacrylamide gel (KB-plus 6.5% gel-matrix - Licor). The scanning of gels subjected to electrophoresis permits the detection of the presence of mutations and their location in the analyzed fragment, deduced on the basis of the size of the fragments obtained. The gels were analyzed using the GelBuddy program in manual mode.
PCR-High Resolution Melting (HRM)
Genetic screening for the presence of mutations in the GLA gene was conducted by Real-Time PCR and subsequent HRM (Figure 2). The PCR was carried out on a final volume of 20 μl, according to the following amplification program: 95°C for 10' minutes, 50 cycles at 95°C for 10 minutes, 62°C for 15 minutes, 72°C for 18 minutes. The amplification in Real-Time PCR and the dissociation of the two DNA strands were performed in the presence of a fluorescent marker (Light Cycler® 480 High Resolution Melting Dye). In order to ensure standardization of the method, male subjects were analyzed by adding a wild type DNA counterpart, in a 1:1 volume ratio.
To test the reliability of the results obtained, all samples were subjected to direct sequencing according to the following program: 94°C for 2 minutes, 30 cycles at 96°C 30 minutes, 58°C for 30 minutes, 72°C for 1 minute, holding at 4°C.
The need to perform simultaneous analyses of many samples led to the adoption of a high-throughput screening strategy based on pools of 8 samples, as reported by Till et al. . In this study, we have carried out genetic screening of the seven chosen target regions for the analysis of the GLA gene on 740 samples (Figure 3). The analysis was done by comparing two enzymes with endonuclease activity (ENDO-1 and CEL I), and their recognition efficiency and heteroduplex cut fragments were identical. The mutations and polymorphisms identified by EcoTILLING are shown in Table 2.
High Resolution Melting
In order to evaluate the results obtained by EcoTILLING, we conducted another screening with High Resolution Melting. The ten mutations/polymorphisms identified by EcoTILLING were confirmed, and two further gene variations (G395A; IVS6-22c>t), undetected by enzymatic digestion by the CEL I and ENDO-1 endonucleases, were identified (Table 3). The mutations which were undetected by EcoTILLING were analyzed separately, applying variations of the protocols: increases in the enzymatic digestion times or changes in analysis of the relationship between mutated and wild type samples in the pool (Figure 4).
After direct sequencing of the 740 samples, no discrepancies were found in the HRM results (Table 4). According to the data, this technique appears to be the optimal method for pre-sequencing screening of the GLA gene. Indeed, the robustness of this method ensures the identification of all possible genomic variants.
The use of a screening system permits the identification of samples to be subjected to sequencing, the performance of the simultaneous analysis of numerous genetic traits, and the undertaking of population studies.
EcoTILLING technology is recognized by the scientific community as a method for the identification and characterization of single nucleotide polymorphisms (SNPs) and small and large deletions (indels) [16–18]. In this study we analyzed pools of eight samples to reduce screening time and/or the cost associated with the analysis of 740 samples, identifying ten genetic variants across the seven exons and their flanking regions of the GLA gene. After HRM analysis, we discovered two other previously undetected mutations.
The detection of mutations with EcoTILLING shows a number of disadvantages that limit its use in analyses of human genetic diseases and in clinical medicine investigations. In particular, mismatch recognition sensitivity is reduced when mutations are present at the ends of the target sequence [15, 19–23] and as the number of polymorphisms identified for each fragment increases.
The acquisition of melting profiles in HRM is not dependent on the location and number of mutations present in the analyzed fragment, but only on changes in the dissociation kinetics of the amplified fragments with respect to the wild type [24–27]. Because of its high reliability and flexibility, it is currently employed as a method of genotyping in clinical and diagnostic settings [28, 29]. In genetics, this method has been applied to BRCA1/2 hereditary breast cancer , cystic fibrosis , haemophilia , the instability of microsatellites , the rearrangement of heavy chain genes in lymphomas , the prenatal diagnosis of B-thalassaemia , analysis of DNA methylation , the identification of bacterial species [37, 38], and other applications. Considering the polymorphic nature of the GLA gene, we sustain that HRM analysis appears to be an ideal investigative tool which permits the search for all possible gene variations with a high level of reliability.
Finally, in the cases reported, some mutations seemed to occur in haplotype patterns, as studies of genealogical trees of the proband families confirmed. Indeed, 12% of the subjects analyzed showed polymorphisms in the promoter region of the GLA gene and, of these, 99% exhibited several simultaneous polymorphisms spread throughout the gene. In particular, the -10 c>t, IVS2-76_80 del5, IVS4-16 A>G, IVS6-22C>T and -12 g>a, IVS4+68 A>G, IVS6-22C>T polymorphisms occured simultaneously in 8.9% and 3.7% of the subjects, and the significance of this haplotype in FD pathology remains unknown. In a pilot study of patients with small fiber neuropathy of unknown etiology, published in August 2010, Tanislav et al. reported the presence of a complex intronic haplotype within the GLA gene (-10 c>t, IVS2-76_80 del5, IVS4-16 A>G, IVS6-22C>T) in four patients which is identical to the first of the haplotypes mentioned above for our subjects . The -12 g>a, IVS4+68 A>G, IVS6-22C>T haplotype identified in approximately 4% of our samples has not been reported in the literature.
Hoffman Bjorn: Fabry disease: recent avances in pathology, diagnosis, treatment and monitoring. Orphanet J Rare Dis. 2009, 4: 21-10.1186/1750-1172-4-21.
Aerts Johannes: Elevated globotriaosylsphingosine is a hallmark of Fabry disease. Proc Natl Acad Sci USA. 2008, 105 (8): 2812-2817. 10.1073/pnas.0712309105.
Shen Jin-Song: Globotriaosylceramide induces oxidative stress and up-regulates cell adhesion molecule expression in Fabry disease endothelial cells. Mol Genet Metab. 2008, 95 (3): 163-168. 10.1016/j.ymgme.2008.06.016.
Linhart Aleš, Elliott Perry: The heart in AndersonıFabry disease and other lysosomal storage disorders. Heart. 2007, 93 (4): 528-535. 10.1136/hrt.2005.063818.
Schiffmann Raphael: Fabry disease: progression of nephropathy, and prevalence of cardiac and cerebrovascular events before enzyme replacement therapy. Nephrol Dial Transplant. 2009, 24 (7): 2102-2111. 10.1093/ndt/gfp031.
Shows TB: Assignment of a-Galactosidase (aGal) to the q22-qter region of the X chromosome in man. Cytogenet Cell Genet. 1978, 22: 541-4. 10.1159/000131020.
Eng CM: Nature and frequency of mutations in the alpha-galactosidase A gene that cause Fabry disease. Am J Hum Genet. 1993, 53 (6): 1186-1197.
Filoni C: Functional studies of new GLA gene mutations leading to conformational fabry disease. Biochimica et Biophysica Acta. 2010, 1802: 247-252.
Lin HY: Novel human pathological mutations. Gene symbol: GLA. Disease: Fabry disease. Hum Genet. 2010, 127 (1): 124.-
Eng CM: Nature and frequency of mutations in the α-galactosidase A gene that cause Fabry disease. Am J Hum Genet. 1993, 53: 1186-97.
Davies JP: Sequence variations in the first exon of α-galactosidase A. J Med Genet. 1993, 30: 658-63. 10.1136/jmg.30.8.658.
Shabbeer J: Detection of α-galactosidase a mutations causing Fabry disease by denaturing high performance liquid chromatography. Hum Mutat. 2005, 25: 299-305. 10.1002/humu.20144.
Barkley NA, Wang M: Application of TILLING and EcoTILLING as Reverse Genetic Approaches to Elucidate the Function of Genes in Plants and Animals. Curr Genomics. 2008, 9 (4): 212-26. 10.2174/138920208784533656.
Montgomery JL: High-resolution DNA melting analysis in clinical research and diagnostics. Expert Rev Mol Diagn. 2010, 10 (2): 219-40. 10.1586/erm.09.84.
Till BJ: High-throughput discovery of rare human nucleotide polymorphisms by Ecotilling. Nucleic Acids Res. 2006, 34 (13): e99.-
Coassin S: An optimized procedure for the design and evaluation of Ecotilling assays. BMC Genomics. 2008, 9: 510-10.1186/1471-2164-9-510.
Till BJ: High-throughput TILLING for functional genomics. Methods Mol Biol. 2003, 236: 205-20.
McCallum CM: Targeted screening for induced mutations. Nat Biotechnol. 2000, 18 (4): 455-7. 10.1038/74542.
Greene EA: Spectrum of chemically induced mutations from a large-scale reversegenetic screen in Arabidopsis. Genetics. 2003, 164: 731-740.
Till BJ: Discovery of chemically induced mutations in rice by TILLING. BMC Plant Biol. 2007, 7: 19-10.1186/1471-2229-7-19.
Raghavan C: Rapid method for detecting SNPs on agarose gels and its application in candidate gene mapping. Mol Breed. 2007, 19: 87-101. 10.1007/s11032-006-9046-x.
Comai L: Efficient discovery of DNA polymorphisms in natural populations by Ecotilling. Plant J. 2004, 37: 778-786. 10.1111/j.0960-7412.2003.01999.x.
Gilchrist EJ: Use of Ecotilling as an efficient SNP discovery tool to survey genetic variation in wild populations of Populus trichocarpa. Mol Ecol. 2006, 15: 1367-1378. 10.1111/j.1365-294X.2006.02885.x.
Wittwer CT: High-resolution genotyping by amplicon melting analysis using LCGreen. Clin Chem. 2003, 49: 853-860. 10.1373/49.6.853.
Erali M: SNP genotyping by unlabeled probe melting analysis. Methods Mol Biol. 2008, 429: 199-206. 10.1007/978-1-60327-040-3_14.
Reed GH: Sensitivity and specificity of single-nucleotide polymorphism scanning by high-resolution melting analysis. Clin Chem. 2004, 50: 1748-1754. 10.1373/clinchem.2003.029751.
Liew M: Genotyping of single-nucleotide polymorphisms by high-resolution melting of mall amplicons. Clin Chem. 2004, 50: 1156-1164. 10.1373/clinchem.2004.032136.
Erali Maria: High Resolution Melting Applications for Clinical Laboratory Medicine. Exp Mol Pathol. 2008, 85 (1): 50-58. 10.1016/j.yexmp.2008.03.012.
Krypuy M: High resolution melting analysis for the rapid and sensitive detection of mutations in clinical samples: KRAS codon 12 and 13 mutations in non-small cell lung cancer. BMC Cancer. 2006, 6: 295-10.1186/1471-2407-6-295.
Takano EA: Rapid detection of carriers with BRCA1 and BRCA2 mutations using high resolution melting analysis. BMC Cancer. 2008, 8: 59-10.1186/1471-2407-8-59.
Krenková P: Evaluation of high-resolution melting (HRM) for mutation scanning of selected exons of the CFTR gene. Folia Biol (Praha). 2009, 55 (6): 238-42.
Laurie AD: Detection of factor VIII gene mutations by high-resolution melting analysis. Clin Chem. 2007, 53 (12): 2211-4. 10.1373/clinchem.2007.093781.
Janavicius R: Microsatellite instability detection by high-resolution melting analysis. Clin Chem. 2010, 56 (11): 1750-7. 10.1373/clinchem.2010.150680.
Kummalue T: Detection of monoclonal immunoglobulin heavy chain gene rearrangement (FR3) in Thai malignant lymphoma by High Resolution Melting curve analysis. Diagn Pathol. 2010, 5: 31-10.1186/1746-1596-5-31.
Pornprasert S, Sukunthamala K: SYTO9 and SYBR GREEN1 with a high-resolution melting analysis for prenatal diagnosis of β-thalassemia/hemoglobin-E. Eur J Haematol. 2010, 85 (5): 424-9. 10.1111/j.1600-0609.2010.01512.x.
Wong EM, Dobrovic A: Assessing gene-specific methylation using HRM-based analysis. Methods Mol Biol. 2011, 687: 207-17. 10.1007/978-1-60761-944-4_14.
Castellanos E: Rapid identification and differentiation of Mycobacterium avium subspecies paratuberculosis types by use of real-time PCR and high-resolution melt analysis of the MAP1506 locus. J Clin Microbiol. 2010, 48 (4): 1474-7. 10.1128/JCM.02484-09.
Ong DC: Rapid detection of rifampicin- and isoniazid-resistant Mycobacterium tuberculosis by high-resolution melting analysis. J Clin Microbiol. 2010, 48 (4): 1047-54. 10.1128/JCM.02036-09.
Tanislav C: Frequency of Fabry disease in patients with small-fibre neuropathy of unknown aetiology: a pilot study. Eur J Neurol. 2010
This work was made possible by a contribution from Shire Italia.
The study was approved by University Hospital Ethics Committee and written informed consent was obtained from all participants.
The authors declare that they have no competing interests.
CB drafted the manuscript and conducted EcoTILLING experiments, DN design the study and conducted HRM experiments, GA carried out the molecular genetic study, CZ performed the statistical analysis and participate in the sequence alignment, DF participate in the sequence alignment, FI participated in acquisitions of a data and interpretations of data, ES participated in acquisitions of a data and interpretations of data, GD participated in the coordination of the study obtaining the informed consent and ethics approval to conduct the study approval. All authors read and approved the final manuscript.
Caterina Bono, Domenico Nuzzo contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.