- Short report
- Open Access
Identification of Trypanosoma cruzi Discrete Typing Units (DTUs) through the implementation of a High-Resolution Melting (HRM) genotyping assay
Parasites & Vectorsvolume 6, Article number: 112 (2013)
Chagas disease, caused by Trypanosoma cruzi, is a geographically widespread anthropozoonosis that is considered a major public health problem in Latin America. Because this parasite presents high genetic variability, a nomenclature has been adopted to classify the parasite into six discrete typing units (DTUs): TcI, TcII, TcIII, TcIV, TcV, and TcVI, which present different eco-epidemiological, clinical, and geographic associations. Currently, the available genotyping methods present a series of drawbacks that implies the need for developing new methods for characterizing T. cruzi DTU’s. The aim of this work was to genotype reference populations from T. cruzi by means of a High-Resolution Melting (HRM) genotyping assay.
To genotype the DTUs of 38 strains and 14 reference clones of T. cruzi from diverse sources, real-time PCR (qPCR) was coupled to high-resolution melting (HRM) based on the amplification of two molecular markers—the divergent domain of the 24 sα rRNA gene and the intergenic region of the mini-exon gene.
Amplification of the mini-exon gene allowed the genotyping of three distinct groups: TcI, TcII- TcIV- TcV, and TcIII-TcVI, while amplification of the 24sα gene generated non-overlapping melting temperature ranges for each DTU that were used to identify the groups in the six existing DTUs of Trypanosoma cruzi.
The proposed genotyping assay allowed discrimination of the six genetic groups by obtaining specific melting curves for each DTU. The application of this technique is proposed because of its specificity, sensitivity, high performance, and low cost compared with other previously described characterization methods.
Chagas disease is an anthropozoonosis caused by the parasite Trypanosoma cruzi and affects an estimated 10 million people, posing a serious public health problem . Trypanosoma cruzi has a high genetic variability with evidence of six Discrete Typing Units (DTU’s): T. cruzi I (TcI), T. cruzi II (TcII), T. cruzi III (TcIII), T. cruzi IV (TcIV), T. cruzi V (TcV), and T. cruzi VI (TcVI), each of which presents characteristics based on the geographic distribution, eco-epidemiological associations, and clinical manifestations of the disease . Currently, several methods have been reported for characterizing and genotyping Trypanosoma cruzi populations, including the PCR amplification of the SL-IR (spliced-leader intergenic region of the mini-exon gene), 24sα, and 18s rDNA gene regions. The comparison criterion for this technique is based on the size differences in the resulting amplicons [3, 4]. However, the characterization obtained using this technique is limited by the absence of amplification products or the need for multiple amplification procedures [2, 5]. Another proposed methodology is based on genetic polymorphism analysis using multilocus PCR-RFLP; however, this technique requires the development of combinations and multiple assays, increasing its complexity and limiting its use [6, 7]. In 2009, Llewellyn and colleagues implemented an approach based on the Multilocus Microsatellite Typing (MLMT) of a panel of 48 polymorphic microsatellite markers. This technique is a quick and feasible alternative but requires the inclusion of sequencing procedures that lead to increased costs . A similar initiative was proposed to develop a Multilocus Sequence Typing (MLST) scheme based on the analysis of a panel of four molecular targets: Met-III, RB19, TcGPXII, and DHFR-TS; however, the development of this process requires both bioinformatic analysis and sequencing steps . Recently, the analysis of the polymorphisms of the 24sα rDNA gene, the subunit II of the cytochrome oxidase, and four microsatellite loci has been shown to be an alternative method . This technique is highly efficient but requires multiple amplification steps that can lead to an increased risk of contamination.
Thus, a low-cost genotyping technique with reliable and robust results for the identification of the six genetic groups in a single assay is needed to facilitate the characterization of the circulating strains in a specific geographic area. A technique such as High-Resolution Melting (HRM) offers a viable alternative for genotyping T. cruzi given its robustness, speed, cost, and the elimination of sequencing or hybridization procedures . This molecular analysis has recently been used to genotype different populations of parasitic microorganisms [12–15]. The aim of this study was to genotype 36 reference strains and 18 reference clones of Trypanosoma cruzi previously genotyped by different molecular markers from different biological and geographical origins via HRM.
We employed 36 strains and 18 reference clones of Trypanosoma cruzi previously characterized via MLST, MLEE, MLMT, and RFLP-PCR (Table 1). Additionally, in order to test specificity we included DNA from reference strains of T. rangeli and Leishmania panamensis. A real-time PCR assay coupled with HRM analysis was developed to target SL-IR and the 24sα rRNA gene using a 7500 Fast Real-Time PCR System (Applied Biosystems, Inc., CA, USA). These genes were selected because they had previously been used as molecular markers for the PCR genotyping of T. cruzi[4, 16]. The amplification reactions were carried out in a total volume of 21 μL. The reaction mixture consisted of 1× MeltDoctor™ HRM Master Mix (Applied Biosystems, Inc., CA, USA); a 5 μM solution of each of the primers TCC (5′CCC CCC TCC CAG GCC ACA CTG 3) and TC2 (5′CCT GCA GGC ACA cGT GTG TGT G 3) for the amplification of SL-IR and 24sα, D71 (5′AAG GTG CGT CGA CAG TGT GG 3) and D72 (5′TTT TCA GAA TGG CCG AAC AGT 3); 6.6 μL of water; and 20 ng/μL DNA template. The cycle conditions for the qPCR were as follows: 95°C for 10 minutes (1 cycle), followed by amplification for 40 cycles at 95°C for 15 seconds (denaturation) and 60°C for 1 minute. The amplicon dissociation was immediately initiated by a melting step. The thermal profile consisted of denaturation at 95°C for 10 seconds, 60°C for 1 minute (annealing), 95°C for 30 seconds (high-resolution melting), and a final annealing stage at 60°C for 15 seconds. The melting curve analysis was carried out by supplying raw fluorescence (F) vs. temperature (T) values or already transformed -dF/dT values. The raw value,differentiation was performed and the first derivative curve is displayed with the melting peaks. Peaks and Tm values were automatically identified, peak areas calculated and peaks flagged as ‘negative’ if the peak area is below the user-defined threshold value. During this process, the PCR amplicons were denatured prior to developing the melting curves, and the changes in fluorescence with respect to the changes in temperature (dF/dT) were recorded by the instrument with a ramp of 0.5°C/seg. Each DNA sample was analyzed in duplicate. Finally, the normalized and derivative melting curve profiles were obtained. The analysis was performed using a blind methodology to reduce the possibility of bias. To evaluate the reproducibility of the assay, the genotyping technique was performed on five independent days under the same conditions. The specificity was assessed by evaluating the concordance between the pre-established genotypes of the DNA samples used and the genotypes characterized by HRM. The consistency of the method was determined by correlating the previously quantified reproducibility and specificity . Finally, a randomization and subsequent non-parametric Kruskal Wallis test (p < 0.05) were conducted to detect any significant differences between the melting temperature data for each DTU. The analyses were performed using the R platform and GraphPad software package.
A specific melting temperature was obtained for each of the reference strains and clones included in this study, both for amplifying SL-IR and for the 24sα. The plots were created for the melting temperature distribution of each DTU (Figure 1). The 24sα gene exhibited higher discriminatory power than the SL-IR gene, producing non-overlapping ranges of the melting temperature for each DTU. The normalized fluorescence curve and the derivative melting curve yielded a different and easily distinguishable profile for each DTU. The representative profiles were collected for the derivative, the normalized melting curves for the SL-IR gene and the 24sα amplicons (Figure 1). Finally, an algorithm was constructed (Figure 2) which allowed the discrimination of the DTU’s after the amplification of the SL-IR gene and the subsequent amplification of the 24sα gene. To validate the method, three parameters were determined based on the methodology of Moncayo and Luquetti (1990) ; the reproducibility was determined for both the SL-IR gene amplification and the 24sα rRNA and yielded a value of 100% for both tests. The specificity and consistency values were obtained only for the amplification of the 24sα rRNA, which was the only marker capable of generating non-overlapping melting temperature ranges that was specific for the six DTU’s and the DNAs from T. rangeli and L. panamensis.
The use of the HRM methodology allowed the adequate and congruent discrimination of the T. cruzi genotypes for the strains and reference clones used. Many problems, such as the need for sequencing, the degree of complexity, the design of multiple pairs of primers, and genotyping by the absence of a band, among others, have been reported with the existing genotyping methods. Thus, new alternatives are needed to simplify the process and minimize the difficulties. For example, methods based on the absence of an amplification product are considered problematic because no actual evidence is generated that may relate to the presence of a specific DTU . To obtain a correct identification, other genotyping methods for which high efficiency, sensitivity, and robustness have been reported require multiple amplifications of a single sample using different markers or coupled with RAPD’s and RFLP analysis, which can significantly increase the economic cost and lead to sample contamination . The HRM technique does not present any of the aforementioned difficulties and, given the as-obtained results, could be considered a potential alternative strategy for genotyping T. cruzi.
The markers were selected based on two criteria: the fact that both the SL-IR gene and the 24sα rRNA have previously been used to characterize T. cruzi and because they exhibit a high number of copies within the parasite genome, which increases the probability of detection [6, 19, 20]. The results of this study indicated that the SL-IR gene did not allow the absolute discrimination of the six DTUs in a single assay, a fact that may be due to the recently confirmed variation in the copy number of the gene and the chromosomal rearrangement at the level of the different strains . To examine the reproducibility of the method, the assays were conducted in duplicate, and the tests were repeated over several days, maintaining the same test conditions, reagent quantities, and DNA concentrations. The technique generated a reproducibility value of 100% for both molecular markers. This result is of great significance because genotyping techniques must be reproducible and robust to avoid producing conflicting results and to accurately establish the genotype of each sample analyzed .
The use of a real-time PCR assay coupled with HRM analysis is proposed for genotyping the DTU’s of the T. cruzi population for the development of research projects that require genotyping of the circulating strains because it provides an easy, quick, and low (1.5 USD per assay compared to 2.5 USD for genotyping) alternative compared to the existing methods. Among the advantages over other currently used methods are its speed and simplicity (given that this method does not require gel electrophoresis), the reduction of the processing time, and the specificity and sensitivity of the method . Therefore, HRM coupled to qPCR provides a cost-effective alternative for the future genotyping and identification of the six DTU’s of T. cruzi which fulfills one of the main research targets in Chagas disease.
Discrete Typing Unit
Trypanosoma cruzi I
Trypanosoma cruzi II
Trypanosoma cruzi III
Trypanosoma cruzi IV
Trypanosoma cruzi V
Trypanosoma cruzi VI
Spliced Leader Intergenic Region.
Rassi A, Rassi A, de Rezende M: American Trypanosomiasis Chagas Disease. Infect Dis Clin North Am. 2012, 262: 275-291.
Zingales B, Miles MA, Campbell DA, Tibayrenc M, Macedo AM, Teixeira MM, Schijman AG, Llewellyn MS, Lages-Silva E, Machado CR, Andrade SG, Sturm NR: The revised Trypanosoma cruzi subspecific nomenclature: Rationale, epidemiological relevance and research applications. Infect Genet Evol. 2012, 122: 240-253.
Brisse S, Verhoef J, Tibayrenc M: Characterization of large and small subunit rRNA and mini-exon genes further supports the distinction of six Trypanosoma cruzi lineages. Int J Parasitol. 2001, 311: 1218-1226.
Souto RP, Fernandes O, Macedo AM, Campbell DA, Zingales B: DNA markers define two major phylogenetic lineages of Trypanosoma cruzi. Mol Biochem Parasitol. 1996, 83: 141-152. 10.1016/S0166-6851(96)02755-7.
Burgos JM, Diez M, Vigliano C, Bisio M, Risso M, Duffy T, Cura C, Brusses B, Favaloro L, Leguizamon MS, Lucero RH, Laguens R, Levin MJ, Favaloro R, Schijman AG: Molecular Identification of Trypanosoma cruzi Discrete Typing Units in End-Stage Chronic Chagas Heart Disease and Reactivation after Heart Transplantation. Clin Inf Dis. 2010, 515: 485-495.
Rozas M, Doncker S, Adaui V, Coronado X, Barnabé C, Tibayrenc M, Solari A, Dujardin JC: Multilocus polymerase chain reaction restriction fragment-length polymorphism genotyping of Trypanosoma cruzi (Chagas disease): taxonomic and clinical applications. J Infect Dis. 2007, 195: 1381-8. 10.1086/513440.
Lewis MD, Ma J, Yeo M, Carrasco HJ, Llewellyn MS, Miles MA: Genotyping of Trypanosoma cruzi: systematic selection of assays allowing rapid and accurate discrimination of all known lineages. AmJTrop Med Hyg. 2009, 81: 1041-1049. 10.4269/ajtmh.2009.09-0305.
Llewellyn MS, Miles MA, Carrasco HJ, Lewis MD, Yeo M, Vargas J, Torrico F, Diosque P, Valente SA, Gaunt MW: Genome-scale multilocus microsatellite typing of Trypanosoma cruzi discrete typing unit I reveals phylogeographic structure and specific genotypes linked to human infection. PLoS Pathog. 2009, 5: e1000410-10.1371/journal.ppat.1000410.
Yeo M, Mauricio IL, Messenger LA, Lewis MD, Llewellyn MS, Acosta N, Bhattachaeyya T, Diosque P, Carrasco HJ, Miles MA: Multilocus Sequence Typing MLST for Lineage Assignment and High Resolution Diversity Studies in Trypanosoma cruzi. PLoS Negl Trop Dis. 2011, 56: e1049-
Valadares HMS, Pimenta JR, Segatto M, Veloso VM, Gomes ML, Chiari E, Gollob KJ, Bahia M, de Lana M, Franco GR, Machado CR, Pena SD, Macedo AM: Unequivocal Identification of Subpopulations in Putative Multiclonal Trypanosoma cruzi Strains by FACs Single Cell Sorting and Genotyping. PLoS Negl Trop Dis. 2012, 67: e1722-10.1371/journal.pntd.0001722.
Schütz E, von Ahsen N: Influencing factors of dsDNA dye high.resolution melting curves and improved genotype call based on thermodynamic considerations. Anal Biochem. 2009, 3851: 143-152.
Zhang P, Liu Y, Alsarakibi M, Li J, Liu T, Li Y, Li G: Application of HRM assays with EvaGreen dye for genotyping Giardia duodenalis zoonotic assemblages. Parasitol Res. 2012, 111: 2157-2163. 10.1007/s00436-012-3064-x.
Al-Mohammed H: Genotypes of Giardia intestinalis clinical isolates of gastrointestinal symptomatic and asymptomatic Saudi children. Parasitol Res. 2011, 108: 1375-1381. 10.1007/s00436-010-2033-5.
Ngui R, Lim YA, Chua KH: Rapid Detection and Identification of Human Hookworm Infections through High Resolution Melting HRM Analysis. PLoS One. 2012, 77: e41996-
Costa JM, Cabaret O, Moukoury S, Bretagne S: Genotyping of the protozoan pathogen Toxoplasma gondii using high-resolution melting analysis of the repeated B1 gene. J Microbiol Met. 2011, 863: 357-363.
Ramírez JD, Guhl F, Rendón LM, Rosas F, Marin-Neto JA, Morillo CA: Chagas Cardiomyopathy Manifestations and Trypanosoma cruzi Genotypes Circulating in Chronic Chagasic Patients. PLoS Neg Trop Dis. 2010, 41: e899-
Moncayo A, Luquetti AO: Multicentre double blind study for evaluation of Trypanosoma cruzi defined antigens as diagnostic reagents. Mem Inst Oswaldo Cruz. 1990, 85: 489-495. 10.1590/S0074-02761990000400020.
Peña VH, Fernández GJ, Gómez-Palacio AM, Mejía-Jaramillo AM, Cantillo O, Triana-Chávez O: High-Resolution Melting HRM of the Cytochrome B Gene: A Powerful Approach to Identify Blood-Meal Sources in Chagas Disease Vectors. PLoS Neg Trop Dis. 2012, 62: e1530-
Cardinal MV, Lauricella MA, Ceballos LA, Lanati L, Marcet PL, Levin MJ, Kitron U, Gürtler RE, Schijman AG: Molecular epidemiology of domestic and sylvatic Trypanosoma cruzi infection in rural northwestern Argentina. Int J Parasitol. 2008, 381: 1533-1543.
Messenger LA, Llewellyn MS, Bhattacharyya T, Franzén O, Lewis MD, Ramírez JD, Carrasco HJ, Andersson B, Miles MA: Multiple Mitochondrial Introgression Events and Heteroplasmy in Trypanosoma cruzi Revealed by Maxicircle MLST and Next Generation Sequencing. PLoS Neg Trop Dis. 2012, 64: e1584-
Minning TA, Weatherly DB, Flibotte S, Tarleton RL: Widespread, focal copy number variations CNV and whole chromosome aneuploidies in Trypanosoma cruzi strains revealed by array comparative genomic hybridization. BMC Genom. 2011, 121: 139-146.
Hong H, Xu L, Liu J, Jones WD, Su Z, Ning B, Perkins R, Ge W, Miclaus K, Zhang L, Park K, Green B, Han T, Fang H, Lambert CG, Vega SC, Lin SM, Jafari N, Czika W, Wolfinger RD, Goodsaid F, Tong W, Shi L: Technical Reproducibility of Genotyping SNP Arrays Used in Genome-Wide Association Studies. PLoS One. 2012, 79: e44483-
We thank IRD and LSHTM for providing the DNA of the reference strains and clones in the framework of the Chagas EpiNet project. Financial support was provided by the Project Chagas EpiNet from The European Union Seventh Framework Programme, contract no.223034.
The authors declare that they have no competing interest.
JDR designed the experiments. SLH and JDR standardized the HRM assay and developed the experiments to calculate the parameters of reproducibility and specificity. SLH, JDR and FG wrote the manuscript. All authors read and approved the final version of the manuscript.
Sonia L Higuera and Juan David Ramírez contributed equally to this work.