High-Resolution Melting (HRM) of the Cytochrome B Gene: A Powerful Approach to Identify Blood-Meal Sources in Chagas Disease Vectors

Abstract

Methods to determine blood-meal sources of hematophagous Triatominae bugs (Chagas disease vectors) are serological or based on PCR employing species-specific primers or heteroduplex analysis, but these are expensive, inaccurate, or problematic when the insect has fed on more than one species. To solve those problems, we developed a technique based on HRM analysis of the mitochondrial gene cytochrome B (Cyt b). This technique recognized 14 species involved in several ecoepidemiological cycles of the transmission of Trypanosoma cruzi and it was suitable with DNA extracted from intestinal content and feces 30 days after feeding, revealing a resolution power that can display mixed feedings. Field samples were analyzed showing blood meal sources corresponding to domestic, peridomiciliary and sylvatic cycles. The technique only requires a single pair of primers that amplify the Cyt b gene in vertebrates and no other standardization, making it quick, easy, relatively inexpensive, and highly accurate.

Abstract

Methods to determine blood-meal sources of hematophagous Triatominae bugs (Chagas disease vectors) are serological or based on PCR employing species-specific primers or heteroduplex analysis, but these are expensive, inaccurate, or problematic when the insect has fed on more than one species. To solve those problems, we developed a technique based on HRM analysis of the mitochondrial gene cytochrome B (Cyt b). This technique recognized 14 species involved in several ecoepidemiological cycles of the transmission of Trypanosoma cruzi and it was suitable with DNA extracted from intestinal content and feces 30 days after feeding, revealing a resolution power that can display mixed feedings. Field samples were analyzed showing blood meal sources corresponding to domestic, peridomiciliary and sylvatic cycles. The technique only requires a single pair of primers that amplify the Cyt b gene in vertebrates and no other standardization, making it quick, easy, relatively inexpensive, and highly accurate.

Author Summary

Chagas disease is one of the most important tropical diseases in America. This disease is caused by the parasite Trypanosoma cruzi and transmitted through the feces of blood-sucking insects known as triatomines. Different species of insects have different habits and food sources that confer variable degrees of epidemiological importance. In this paper, we propose the use of High Resolution Melting (HRM) analysis of cytochrome b (cyt b) gene PCR products to identify blood-food sources in triatomines. This tool can effectively differentiate blood-meal sources of insects collected from the field. Such data allows for targeted investigations of insect species that are likely to be involved in the transmission of the parasite to humans in different regions. This knowledge is very important because it allows establishing and designing vector control and surveillance strategies according to each geographical area and to stop the transmission of the parasite to human populations by insects.

Citation: Peña VH, Fernández GJ, Gómez-Palacio AM, Mejía-Jaramillo AM, Cantillo O, et al. (2012) High-Resolution Melting (HRM) of the Cytochrome B Gene: A Powerful Approach to Identify Blood-Meal Sources in Chagas Disease Vectors. PLoS Negl Trop Dis 6(2): e1530. doi:10.1371/journal.pntd.0001530

Editor: Daniel K. Masiga, International Centre of Insect Physiology and Ecology, Kenya

Received: October 20, 2011; Accepted: January 2, 2012; Published: February 28, 2012

Copyright: © 2012 Peña et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: Financial support was provided by Colciencias, project number 111549326149. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Trypanosoma cruzi, the etiological agent of Chagas disease, is transmitted by hematophagous Reduviidae insects belonging to Triatominae subfamily, which are mostly distributed on the American continent and comprises 141 species in 15 genera [1], [2]. T. cruzi can infect more than 180 species of terrestrial and arboreal mammals belonging to nine orders and 25 subfamilies, which are natural blood-meal sources of triatomine bugs [3].

Depending on the habitats of triatomines, at least two transmission cycles of T. cruzi have been identified: domestic and sylvatic. The first one includes humans, domestic animals, and domiciliated vectors. The second one involves sylvatic insects and wild animals [4]. A third cycle associated with anthropic environments can occur in several places where domestic animals are kept in sheds near dwellings, called by some authors a peridomestic cycle [5]. However, the epidemiological scenario of Chagas disease is more complex because it can include overlapping cycles [6]-[10], intrusions of wild vectors into human dwellings [11], [12], colonization and re-infestation of new sylvatic vector species after spraying houses [4], [13]-[15], the regular intrusion of wild mammals inside houses [16], house infestation by two vector species [17], etc.

The identification of T. cruzi's transmission cycle has been considered an essential issue to establish and design vector control and surveillance strategies [18]. In this sense, the determination of blood-meal sources in triatomine insects is essential to identify infection foci, risk factors, and new vector and mammal species involved in the transmission cycles [19]. In triatomines, so far, the methods used for this purpose include serologic and PCR-based molecular techniques. ELISA-based serological techniques are the most frequently used in blood-meal determination of Chagas disease vectors [20]-[24]. However, this method lacks specificity because it cross-reacts with species within the same family and is time-consuming and very expensive because it requires the preparation of specific antisera for each vertebrate host involved [25].

Unlike serological tools, PCR-based molecular methods are more specific and usually easier to perform, although the design of specific primers for each species involved in transmission scenarios makes it difficult to identify the vector and mammals from the wild cycle where a wide variety of species exist. Additionally, this technique is wasteful because each sample has to be amplified by PCR many times with different molecular markers [26]-[28]. A more powerful approach [29]-[31] is the heteroduplex assay of the mitochondrial cytochrome B (Cyt b) gene (reviewed in [19]). This technique can detect differences in the DNA sequences of this gene among species by differences on electrophoretic profiles [29]. Although it is increasingly used in triatomines, standardization is often difficult, and when insects have fed on more than one species the application and the analysis could be complex.

Recently, high-resolution melting (HRM) has become a sensitive genotyping method, based on the characteristics of thermal denaturation of the amplicons. This method has a much higher performance information never before achieved by classical DNA melting curve analysis [32]. HRM is performed using a fluorescent double-stranded DNA dye that can be used in fully saturating conditions [33], [34]. The amplicon is analyzed by gradual denaturation through increasing temperature and decreased fluorescence caused by the release of intercalating dye from DNA. The melting temperature (Tm) and specific shape of the melting curve result from the DNA sequence, GC content, and amplicon length [35], [36], so it is a useful technique to obtain species-specific genotypes. Due to the increased demand for rapid, economic, easy, and high-throughput genotyping analyses, there has been a considerable focus on HRM, which can detect sequence variants without sequencing or hybridization procedures [37]-[40]. Here we present a fast, inexpensive, accurate, and sensitive technique to identify blood-meal sources of triatomines based on Cyt b gene HRM analysis of host-specific genotypes and its application in field studies.

Materials and Methods

Ethics statement

All Animals were handled in strict accordance with good animal practice as defined by the Colombian code of practice for the care and use of animals for scientific purposes. Ethical approval (Act N° 53, 30/06/2009) for analyzing animal specimens was obtained from the animal ethics Committee of the University of Antioquia, Medellin, Colombia.

Identification of HRM-Cyt b profiles in standard samples

To differentiate among species-specific genotypes of natural blood-meal sources of triatomine bugs, a 383-bp fragment from the Cyt b gene was amplified from DNA samples extracted from tissue or blood of 36 individuals corresponding to 14 species. Species were chosen according to their epidemiological relevance as reservoirs of T. cruzi or by their proximity to anthropic environments. The species, number of individuals per species, and type of sample used to obtain DNA are shown in Table 1.

Table 1. Source of the sample of the standard species used in the HRM analysis.

Validation of HRM-Cyt b profiles in samples from triatomines

To validate HRM-Cyt b profiles, six fifth-stage nymphs per species of Rhodnius prolixus, R. colombiensis, and Triatoma maculata maintained under controlled laboratory conditions of temperature (28°C) and relative humidity (70%) were fed with chicken blood until satiated. Feces from three individuals of each species and intestinal content from three other individuals were collected 5 days post-feeding for posterior DNA extraction and amplification of the Cyt b gene. Additionally, the sensitivity of HRM Cyt b over time was analyzed in two groups of five R. prolixus adults fed with chicken or human blood. The intestinal content and feces were collected 1, 5, 15, and 30 days post-feeding and processed as below.

To evaluate the capability of this technique to detect HRM-Cyt b profiles in mixed feeds, five fifth-instar R. prolixus individuals were fed with both mice and chicken blood and the intestinal content and feces were collected at 5 days post-feeding for subsequent HRM analysis.

Application of the HRM Cyt b technique to identify blood-meal sources from field triatomines

To apply the HRM technique to triatomines from the field, a total of 20 insects were collected in four localities of Colombia's Caribbean region, which display different eco-epidemiological transmission cycles of T. cruzi, as follows: (1) six domiciliary R. prolixus and three peridomiciliary T. dimidiata were collected by timed manual capture in an indigenous community with high epidemiological risk from the Sierra Nevada de Santa Marta area [41]; (2) three peridomiciliary T. maculata and one sylvatic Eratyrus cuspidatus (light-trap captured) were collected in the town of Talaigua Nuevo, located in the Bolivar department where a moderate epidemiological risk exists [42]; (3) three sylvatic T. dimidiata were collected manually in a hole in a fallen tree in Turbo, near the border between Colombia and Panama in the Antioquia department; this region is considered a nonendemic zone and a few reports on triatomine visiting houses have been released [43]; (4) finally, four sylvatic R. pallescens were collected in a palm tree in Aguachica in the Cesar department, where nondomiciliary triatomines have been reported [43]; they were caught using live chicken-bait traps (Table S1).

All insects were placed in plastic bottles, marked, transported to the laboratory, and identified according to the classification proposed by Lent and Wygodzinsky [44]. Feces samples from insects were collected and analyzed by HRM of the Cyt b gene to identify the blood-meal sources of each insect.

DNA extraction

DNA was extracted using the phenol-chloroform method [45], with some modifications. Briefly, 200 µL of serum, 25 mg of tissue, 25 µL of intestinal content, or 30 µL of feces were incubated at 37°C for 4 h with 200 µg of proteinase K, 1% sodium dodecyl sulfate, 2.5 mM disodium EDTA, and 25 mM sodium acetate. DNA was extracted with phenol and chloroform, precipitated with 0.3 M sodium acetate and absolute ethanol, washed with 70% ethanol, vacuum dried, and then dissolved in 50 µL of distilled water. Quality and DNA concentration were measured using a Nanodrop 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA).

Amplification of the Cyt b gene using real-time PCR

A 383-bp fragment from the Cyt b gene was real time PCR-amplified from genomic DNA samples using the 5'-CCCCTCAGAATGATATTTGTCCTCA-3' and 5'- CCATCCAACATCTCAGCATGATGAAA-3' primers [29]. Real time PCR was performed at a final volume of 10 µL, containing 0.5 µM of each primer, 0.25 mM of each dNTP, 3.5 mM MgCl2, 0.5 U TrueStart Hot Start polymerase, 1.5× SYBR green dye (Invitrogen), and 10 ng of genomic DNA. Each sample was amplified in triplicate. The thermal cycling conditions were as follows: pre-heating at 95°C for 5 min, 35 cycles at 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s in a thermal cycler Rotor Gene Q system (Qiagen). Following the real time PCR, melting-curve analysis of amplicons was conducted in the thermocycler (Rotor-Gene Q) by increasing the temperature from 72°C to 92°C at ramping increments from 0.1°C/s, recording changes in fluorescence with changes in temperature (dF/dT), and plotting against changes in temperature. DNA extracted from triatomine legs and T. cruzi were used as negative control.

High-resolution melt curve analyses

HRM analysis was carried out using the Rotor Gene Q software v2.2 with normalization regions between 76.15-78.65°C and 89.50-91.00°C. Genotypes were defined by selecting a sample from each standard species as a reference control to identify unknown samples. The software then auto-called the genotype and melting temperatures of each amplicon and provided a confidence percentage based on the square root of the correlation coefficient between samples and the reference genotypes. Any specimen generating melt curves but having a dF/dT less than 1.0 was considered not subjected to genotyping. Averages of melting temperatures and the confidence percentage of the specimen replicates were assigned to a representative genotype, and the standard deviation, confidence intervals, and variation coefficient were calculated using the Prism GraphPad v4.0 software (GraphPad software, Inc).

Sequencing

To confirm the sequence identity of those samples tested by HRM analysis, a 358-bp fragment of Cyt b was sequenced for the standard samples and ten test samples chosen randomly (2, 3, 4, 5, 10, 12, 15, 16, 19, and 20 (Table S1)). For each specimen, both forward and reverse sequences were used to generate a consensus sequence using Bioedit v. 7.0.5 [46], and then positional nucleotide sites were compared after multiple alignment done using the ClustalW algorithm [47] implemented in Bioedit v. 7.0.5 [48]. Additionally, a nucleotide Basic Local Alignment Search (Blastn) was performed to estimate the matched hit scores, identity percentages and e-values of the Cyt b tested. Finally, the distance tree based on net genetic distances (p-distances) was performed using the neighbor-joining algorithm with 1000 bootstrap replicates using MEGA 5.05 [49].

Results

Genotyping of standard samples using HRM analysis of the Cyt b gene

All samples analyzed amplified a 383-bp product, as previously reported to the Cyt b gene [29] (data not shown) and no product was obtained when the insect DNA or T. cruzi DNA were used. The amplicon melt curve and Tm showed low intraspecies variability (Table 2, Figure 1). A specific profile for each species analyzed was observed (Figures 2 and S1), with confidence percentages (%C) ranging from 78.81% and 98.69% (Table 2), indicating that samples were correctly identified in each of the experiments conducted in triplicate.

[Figure omitted, see PDF]

Figure 1. Tm distribution of 14 genotypes included in the study.

Tm: melting temperature, Sp: species.

[Figure omitted, see PDF]

Figure 2. Dissociation curve of 13 of the 14 genotypes included in the study.

Table 2. Fourteen standard species included in the study and its respective Tm and HRM analysis.

Although some species exhibited the same Tm value, HRM profiles in these species were clearly discriminated and recognized as different genotypes, such as dog and sheep or pig and horse (Table 2, Figure S1). On the other hand, although the cow melt curve showed two peaks, the genotype was recognized with a high confidence percentage (Table 2). These two peaks were found in different types of samples (muscle and blood).

Validation of the HRM technique to identify blood-meal sources in triatomines

Analysis of intestinal contents and feces from different species of insects.

The HRM profile of the Cyt b gene obtained from feces and intestinal content from R. colombiensis, R. prolixus, and T. maculata fed with chicken blood were analyzed and classified based on the genotypes defined with the standard samples. All samples were correctly classified as blood-fed chicken (C) and just one of the samples obtained from feces (R. prolixus) showed low confidence percentage. Table S2 summarizes the results of the amplifications obtained for each of these species.

Analysis of intestinal content and feces at different times post-feeding and after mixed feeding.

All samples were successfully amplified by qPCR 1, 5, 15, and 30 days post-feeding (Table 3). Intestinal content and feces samples showed successful amplification and detection of the blood-meal source after 30 days post-feeding, indicating the high sensitivity of Cyt b-HRM at this time. Interestingly, the samples from insects fed on humans showed one peak classified as a chicken genotype at 1 and 5 days post-feeding. This result is not surprising considering that insects were maintained in laboratory conditions with chicken blood meal before feeding with human blood. However, after 15 days post-feeding, a second peak appeared which corresponded to the human genotype (Table 3).

Table 3. Samples from intestinal content and feces at different times after feeding.

Therefore, we evaluated the possibility of detecting at least two different blood meal sources. In this case, the melt curve showed two peaks belonging to mouse and chicken amplicons, with their respective Tm. In addition, the dissociation curve had a shape with two falls (Figure S2). These types of curves were analyzed by modifying the normalization area to every fall of the curve to recognize genotypes separately with their respective standard genotype.

Identification of blood sources in triatomines from the field

The HRM analysis of Cyt b in R. prolixus collected inside dwellings indicated that five of the six individuals evaluated had at least a human blood-meal source, and the remaining one a dog blood-meal source (Table 4). Two peridomiciliary T. dimidiata showed mixed blood-meal sources (human and dog; human and a nonidentified source) (Table 4). Sylvatic T. dimidiata showed a mixed blood-meal source of opossum, human, and a nonidentified source. Of two T. maculata caught in a house, one was found to have fed on a human blood-meal source, while the other had a mixed human and nonidentified source. Finally, the sylvatic E. cuspidatus and R. pallescens were identified as being blood-fed from human, chicken, and nonidentified sources (Table 4, Figure 3).

[Figure omitted, see PDF]

Figure 3. Field sample showing mixed feeding including opossum and human with their respective controls.

Table 4. Amplification results from field samples (Table S1) showing one or two peaks in the melting curve.

As mentioned above, four mixed samples, including one intradomiciliary insect, showed a peak with a Tm = 84.83±0.05, which was inconsistent with any of our standards (Table 4), indicating another additional blood-meal source that was not considered in this study. Five samples did not amplify.

Sequence analysis of blood sources in triatomines confirmed the HRM result

Blastn results showed significant identity values greater than 95% for all samples except for sample number 20 (Table 5). Multiple alignment of the test sample and standards showed great nucleotide similarity and coherence with HRM results previously obtained. The percentage of identity among the field-collected samples and their respective standards ranged between 95% and 99% in samples from human and opossum, and 57% in the sample from chicken (Table 5). Finally, the distance tree showed well-supported groups according to HRM results (Figure S3). It is worth mentioning that the field-collected samples sequenced included a sample (number 12) that was positive for mixed feeding for HRM analysis. In this case, the sequence analysis also identified the two species (human (12_1) and opossum (12_2)) (Table 5).

Table 5. Blastn results of ten insects collected from the field.

Discussion

The epidemiological scenario of Chagas disease has become increasingly complex over the years. The natural habitats of some human populations within the forest and deforestation caused by humans are but two of the reasons that may complicate this scenario. Thus, the classical separation of transmission cycles defined for this disease could be different for many places, making it difficult to determine the epidemiological characteristics of particular regions. With this unclear scenario, the determination of blood-meal sources in hematophagous vectors has become essential to surveillance and prevention of potential infection foci.

In this study, HRM analysis of the Cyt b gene made it possible to identify 14 species successfully even when some of them had the same Tm values. Each species was well recognized under a variety of species reaching high confidence percentage values, thus showing the power of this classification using gene amplification and HRM analysis. It is worth noting that this recognition is possible with a single PCR, making it a quick and inexpensive technique.

It is important to highlight that no studies with this number of species standards have been conducted with Chagas disease vectors. Pizarro and Stevens [28] designed primers against different targets from 11 species and Mota et al. [26] also identified 11 species in their study. ELISA-based techniques could identify the same number of species, but it is very expensive and, in most cases, anti-serum for many species is not commercially available. On the other hand, studies based on PCR methods may also increase the number of species, but it is necessary to standardize PCR conditions for each species. Additionally, every DNA sample must be submitted to a number of PCRs to identify the blood sources, making the technique very expensive. Instead, the technique used herein only needs a single pair of primers, making it easy to include other species according to every epidemiological scenario with no other standardizations of PCR conditions.

A similar power is reached when DNA extracted from intestinal content and feces is analyzed with this technique, demonstrating that both can be used for the identification of the blood source from insects. Although no misidentification was recorded with DNA from insect feces, low confidence percentages were reached and other samples did not amplify, showing that the quality of the sample is a limitation of the technique, as reported for real-time PCR [50]. Nonetheless, work with feces has substantial advantages because it is sometimes hard to extract intestinal content from insects, especially when they have been starved for a long time. Another benefit of feces is the possibility of making a diagnosis of T. cruzi infection and the determination of the blood-meal source; both tests that can be done with a single sample of DNA without sacrificing the insect.

The HRM technique proved to work in both types of samples (intestinal and fecal) at least until 30 days after insects were fed. This showed the applicability of the technique even when the blood-meal sources were limited to the triatomines. Pizarro and Stevens [28] reported on the feasibility of their technique, based on the primer design for each species, for 2 months. It is probable that our technique has identical results, but this was not verified. However, many factors have an influence on the presence of blood in the intestinal tract over time, as described by Lee et al. [30]. Chicken blood has nucleated erythrocytes, so a large amount of DNA should remain in the intestine for a long time. This could explain the presence of unspecific peaks in the melt curve when insects were fed with human blood because the insects are usually fed with chicken blood in laboratory conditions. Thus, the present results suggest that chicken blood remains in the intestinal tract longer than mouse and human blood (data not shown).

On the other hand, the cow profile showed two melting peaks, suggesting the amplification of nonspecific fragments. This could stem from nuclear copies of mtDNA (Numts). It is well known that Numts can be amplified in genetic studies based on mtDNA [51], [52]. Some of these Numts have been described in cattle [53]. However, this phenomenon does not produce mistakes in classification because the shape remains the same across different specimens included.

Unspecific peaks led to testing mixed feedings, showing that the technique detects both Tm sources. This allows identification of many species from which the insects have been fed, even when there are no standards for those species. This is impossible with methods such as ELISA or techniques based on PCR-specific primers. Heteroduplex analysis identifies mixed feeding as well, but the electrophoretic profile becomes very complex, making it difficult to analyze, and sometimes it does not recognize even a single species within the mixture [30]. Sequencing of DNA fragments can resolve mixed feedings well [54], but is quite difficult to analyze and time-consuming. However, we suggest using sequencing of amplicons only when unidentified samples by HRM analysis are present.

We processed 20 samples from the field but only 15 amplified. The other five samples were probably of insufficient quality, which is not surprising because it was feces DNA.

In general, the expected genotypes were recognized by the analysis. Those insects captured inside houses showed human feeding and one dog feeding, suggesting a domestic transmission cycle involving humans and to a lesser extent dogs, which has been reported to have epidemiological relevance [55]-[58]. Samples from the peridomestic or sylvatic cycle also showed human feeding, in agreement with other studies showing that R. prolixus and T. dimidiata fed mostly on human blood [23]. Sometimes, mixed feeding with human blood showed a low peak in the melt curve (Figure 3). It is possible that the other techniques were not sensitive enough to detect human blood because of the small amount of DNA in the sample [31].

We included R. pallescens caught with live-bait trap, which showed chicken-blood feeding as expected by the type of bait used. This control from the field confirms the capacity and accuracy of the technique even when field samples are analyzed, which is the aim of technique.

E. cuspidatus showed human feeding, which is surprising because is a strictly sylvatic species. However, Dib et al. [59] incriminated E. cuspidatus in transmission of parasite when comparing RAPD profiles from T. cruzi strain isolates from human and vector. These two results showed that although sylvatic this species probably visits dwellings looking for feeding sources. Additionally, Cortés and Suárez [42] reported E. cuspidatus inside houses in the study area, suggesting that the insect occasionally colonized intra- and peridomestic environments, which is confirmed by the present result. More studies must be conducted to determine the role of this species in transmission of T. cruzi in this area.

There was only one genotype that could not be identified with a Tm = 84.83°C±0.05, indicating that the number of species standards included in the study was adequate. However, this could be improved knowing the fauna and diversity of species of a particular study area. When the sample is not recognized, the next step is sequencing, as reported by other authors [29]-[31]. Based on these results, we recommend the use of Tm values as a first approach to identify a sample from insects collected in natural populations. Then the species selected based on the Tm could be used as a standard sample for HRM analysis, and finally the species not identified could be sequenced.

In conclusion, we believe that epidemiological studies involving vectorial incrimination and transmission dynamics must identify the blood-meal source to cover the entire panorama of transmission. HRM analysis of the Cyt b gene is the most powerful technique in this type of study because it can accurately identify the species even when the vector has mixed feeding, it has a high resolution power, and it is fast, easy, and inexpensive. However, it is important to obtain high-quality DNA and be mindful of the fauna of the study area to have an adequate number of standard species.

Supporting Information

Figure S1.

Discrimination of genotypes with similar Tm. A: sheep and dog genotypes with Tm of 82.29 for both. B: horse and pig genotypes with Tm of 82.56 and 82.58, respectively.

(TIF)

Figure S2.

Dissociation curve of DNA sample from insect with mixed feeding with mouse (Tm = 81.63) and chicken (Tm = 86.27). Arrows indicate falls of the curve indicating each genotype.

(TIF)

Figure S3.

Distance tree based on net genetic distances (p-distances) with neighbor-joining algorithm with 1000 bootstrap replicates. Gavialis gangeticus was used as outgroup.

(PDF)

Table S1.

Field samples included in this study and geographic distribution. Location corresponds to locality followed by the department. Ecotope corresponds to the place where the bugs were collected. SNSM: Sierra Nevada de Santa Marta.

(DOCX)

Table S2.

Tm and HRM analysis from DNA samples from feces and intestinal content of three triatominae species. IC: Intestinal content, F: Feces; Tm: Melting Temperature, %C: confidence percentage, %C SD: Standard deviation of confidence percentage, C: Chicken.

(DOCX)

Author Contributions

Conceived and designed the experiments: AMGP AMMJ GJF OTC. Performed the experiments: VHP GJF. Analyzed the data: VHP GJF AMGP AMMJ OC OTC. Contributed reagents/materials/analysis tools: AMMJ OTC. Wrote the paper: VHP GJF AMGP AMMJ OC OTC.

References

Schofield CJ, Galvão C (2009) Classification, evolution, and species groups within the Triatominae. Acta Trop 110: 88-100. doi: 10.1016/j.actatropica.2009.01.010. Find this article online

Jurberg J, Rocha DdS, Galvão C (2009) Rhodnius zeledoni sp. nov. similar to Rhodnius paraensis Sherlock, Guitton & Miles, 1977 (Hemiptera, Reduviidae, Triatominae). Biota Neotrop 9: 123-128. Find this article online

Herrera L (2010) Una revisión sobre reservorios de Trypanosoma (Schizotrypanum) cruzi (Chagas, 1909), agente etiológico de la Enfermedad de Chagas. Bol Mal Salud Amb 50: 3-15. Find this article online

Miles MA, Feliciangeli MD, de Arias AR (2003) American trypanosomiasis (Chagas' disease) and the role of molecular epidemiology in guiding control strategies. BMJ 326: 1444-1448. doi: 10.1136/bmj.326.7404.1444. Find this article online

Zingales B, Souto RP, Mangia RH, Lisboa CV, Campbell DA, et al. (1998) Molecular epidemiology of American trypanosomiasis in Brazil based on dimorphisms of rRNA and mini-exon gene sequences. Int J Parasitol 28: 105-112. doi: 10.1016/S0020-7519(97)00178-1. Find this article online

Fernandes AJ, Diotaiuti L, Dias JC, Romanha AJ, Chiari E (1994) [Relationships between Trypanosoma cruzi transmission cycles in the county of Bambuí, Minas Gerais, Brazil]. Cad Saude Publica 10: 473-480. Find this article online

Diotaiuti L, Pereira AS, Loiola CF, Fernandes AJ, Schofield JC, et al. (1995) Inter-relation of sylvatic and domestic transmission of Trypanosoma cruzi in areas with and without domestic vectorial transmission in Minas Gerais, Brazil. Mem Inst Oswaldo Cruz 90: 443-448. doi: 10.1590/S0074-02761995000400002. Find this article online

Téllez-Meneses J, Mejía-Jaramillo AM, Triana-Chávez O (2008) Biological characterization of Trypanosoma cruzi stocks from domestic and sylvatic vectors in Sierra Nevada of Santa Marta, Colombia. Acta Trop 108: 26-34. doi: 10.1016/j.actatropica.2008.08.006. Find this article online

Cardinal MV, Lauricella MA, Ceballos LA, Lanati L, Marcet PL, et al. (2008) Molecular epidemiology of domestic and sylvatic Trypanosoma cruzi infection in rural northwestern Argentina. Int J Parasitol 38: 1533-1543. doi: 10.1016/j.ijpara.2008.04.010. Find this article online

Rodríguez IB, Botero A, Mejía-Jaramillo AM, Marquez EJ, Ortiz S, et al. (2009) Transmission dynamics of Trypanosoma cruzi determined by low-stringency single primer polymerase chain reaction and southern blot analyses in four indigenous communities of the Sierra Nevada de Santa Marta, Colombia. Am J Trop Med Hyg 81: 396-403. Find this article online

Costa J (1999) The synanthropic process of Chagas disease vectors in Brazil, with special attention to Triatoma brasiliensis Neiva, 1911 (Hemiptera, reduviidae, triatominae) population, genetical, ecological, and epidemiological aspects. Mem Inst Oswaldo Cruz 94: Suppl 1239-241. doi: 10.1590/S0074-02761999000700038. Find this article online

Dias-Lima AG, Sherlock IA (2000) Sylvatic vectors invading houses and the risk of emergence of cases of Chagas disease in Salvador, State of Bahia, Northeast Brazil. Mem Inst Oswaldo Cruz 95: 611-613. doi: 10.1590/S0074-02762000000500004. Find this article online

Noireau F, Dujardin JP (2001) Flight and nutritional status of sylvatic Triatoma sordida and Triatoma guasayana. Mem Inst Oswaldo Cruz 96: 385-389. doi: 10.1590/S0074-02762001000300018. Find this article online

Costa J, Almeida CE, Dotson EM, Lins A, Vinhaes M, et al. (2003) The epidemiologic importance of Triatoma brasiliensis as a chagas disease vector in Brazil: a revision of domiciliary captures during 1993-1999. Mem Inst Oswaldo Cruz 98: 443-449. doi: 10.1590/S0074-02762003000400002. Find this article online

Guhl F, Pinto N, Aguilera G (2009) Sylvatic triatominae: a new challenge in vector control transmission. Mem Inst Oswaldo Cruz 104: Suppl 171-75. doi: 10.1590/S0074-02762009000900012. Find this article online

Roque AL, Xavier SC, da Rocha MG, Duarte AC, D'Andrea PS, et al. (2008) Trypanosoma cruzi transmission cycle among wild and domestic mammals in three areas of orally transmitted Chagas disease outbreaks. Am J Trop Med Hyg 79: 742-749. Find this article online

Feliciangeli MD, Carrasco H, Patterson JS, Suarez B, Martínez C, et al. (2004) Mixed domestic infestation by Rhodnius prolixus Stal, 1859 and Panstrongylus geniculatus Latreille, 1811, vector incrimination, and seroprevalence for Trypanosoma cruzi among inhabitants in El Guamito, Lara State, Venezuela. Am J Trop Med Hyg 71: 501-505. Find this article online

Silveira A, Vinhaes M (1999) Elimination of vector-borne transmission of Chagas disease. Mem Inst Oswaldo Cruz 94: Suppl 1405-411. doi: 10.1590/S0074-02761999000700080. Find this article online

Kent RJ (2009) Molecular methods for arthropod bloodmeal identification and applications to ecological and vector-borne disease studies. Mol Ecol Resour 9: 4-18. doi: 10.1111/j.1755-0998.2008.02469.x. Find this article online

Zeledón R, Solano G, Zúñiga A, Swartzwelder JC (1973) Biology and ethology of Triatoma dimidiata (Latreille, 1811). 3. Habitat and blood sources. J Med Entomol 10: 363-370. Find this article online

Jirón LF, Zeledón R (1982) [Feeding preferences in 3 species of Triatominae (Hemiptera: Reduviidae) in experimental conditions]. Rev Biol Trop 30: 151-159. Find this article online

Gurtler RE, Cohen JE, Cecere MC, Chuit R (1997) Shifting Host Choices of the Vector of Chagas Disease, Triatoma Infestans, in Relation to the Availability of Host in Houses in North-West Argentina. Journal of Applied Ecology 34: 699-715. doi: 10.2307/2404917. Find this article online

Sasaki H, Rosales R, Tabaru Y (2003) Host feeding profiles of Rhodnius prolixus and Triatoma dimidiata in Guatemala (Hemiptera: Reduviidae: Triatominae). Medical Entomology and Zoology 54: 283-289. Find this article online

Pineda V, Montalvo E, Alvarez D, Santamaría AM, Calzada JE, et al. (2008) Feeding sources and trypanosome infection index of Rhodnius pallescens in a Chagas disease endemic area of Amador County, Panama. Rev Inst Med Trop Sao Paulo 50: 113-116. doi: 10.1590/S0036-46652008000200009. Find this article online

Chow E, Wirtz RA, Scott TW (1993) Identification of blood meals in Aedes aegypti by antibody sandwich enzyme-linked immunosorbent assay. J Am Mosq Control Assoc 9: 196-205. Find this article online

Mota J, Chacon JC, Gutiérrez-Cabrera AE, Sánchez-Cordero V, Wirtz RA, et al. (2007) Identification of blood meal source and infection with Trypanosoma cruzi of Chagas disease vectors using a multiplex cytochrome b polymerase chain reaction assay. Vector Borne Zoonotic Dis 7: 617-627. Find this article online

Pizarro JC, Lucero D, Stevens L (2007) A method for the identification of guinea pig blood meal in the Chagas disease vector, Triatoma infestans. Kinetoplastid Biol Dis 6: 1. Find this article online

Pizarro JC, Stevens L (2008) A new method for forensic DNA analysis of the blood meal in chagas disease vectors demonstrated using Triatoma infestans from Chuquisaca, Bolivia. PLoS One 3: e3585. doi: 10.1371/journal.pone.0003585. Find this article online

Boakye DA, Tang J, Truc P, Merriweather A, Unnasch TR (1999) Identification of bloodmeals in haematophagous Diptera by cytochrome B heteroduplex analysis. Med Vet Entomol 13: 282-287. doi: 10.1046/j.1365-2915.1999.00193.x. Find this article online

Lee JH, Hassan H, Hill G, Cupp EW, Higazi TB, et al. (2002) Identification of mosquito avian-derived blood meals by polymerase chain reaction-heteroduplex analysis. Am J Trop Med Hyg 66: 599-604. Find this article online

Bosseno MF, García LS, Baunaure F, Gastelúm EM, Gutierrez MS, et al. (2006) Identification in triatomine vectors of feeding sources and Trypanosoma cruzi variants by heteroduplex assay and a multiplex miniexon polymerase chain reaction. Am J Trop Med Hyg 74: 303-305. Find this article online

Wilhelm J, Pingoud A, Hahn M (2003) Validation of an algorithm for automatic quantification of nucleic acid copy numbers by real-time polymerase chain reaction. Anal Biochem 317: 218-225. doi: 10.1016/S0003-2697(03)00167-2. Find this article online

Chateigner-Boutin AL, Small I (2007) A rapid high-throughput method for the detection and quantification of RNA editing based on high-resolution melting of amplicons. Nucleic Acids Res 35: e114. doi: 10.1093/nar/gkm640. Find this article online

Wittwer CT, Reed GH, Gundry CN, Vandersteen JG, Pryor RJ (2003) High-resolution genotyping by amplicon melting analysis using LCGreen. Clin Chem 49: 853-860. doi: 10.1373/49.6.853. Find this article online

Wittwer CT (2009) High-resolution DNA melting analysis: advancements and limitations. Hum Mutat 30: 857-859. doi: 10.1002/humu.20951. Find this article online

Vossen RH, Aten E, Roos A, den Dunnen JT (2009) High-resolution melting analysis (HRMA): more than just sequence variant screening. Hum Mutat 30: 860-866. doi: 10.1002/humu.21019. Find this article online

Reed GH, Wittwer CT (2004) Sensitivity and specificity of single-nucleotide polymorphism scanning by high-resolution melting analysis. Clin Chem 50: 1748-1754. doi: 10.1373/clinchem.2003.029751. Find this article online

Mackay JF, Wright CD, Bonfiglioli RG (2008) A new approach to varietal identification in plants by microsatellite high resolution melting analysis: application to the verification of grapevine and olive cultivars. Plant Methods 4: 8. doi: 10.1186/1746-4811-4-8. Find this article online

Wu SB, Wirthensohn MG, Hunt P, Gibson JP, Sedgley M (2008) High resolution melting analysis of almond SNPs derived from ESTs. Theor Appl Genet 118: 1-14. doi: 10.1007/s00122-008-0870-8. Find this article online

Tindall EA, Petersen DC, Woodbridge P, Schipany K, Hayes VM (2009) Assessing high-resolution melt curve analysis for accurate detection of gene variants in complex DNA fragments. Hum Mutat 30: 876-883. doi: 10.1002/humu.20919. Find this article online

Montilla M, Soto H, Parra E, Torres M, Carrillo P, et al. (2011) Infestation by triatomine bugs in indigenous communities of Valledupar, Colombia. Rev Saude Publica 45: 773-780. Find this article online

Cortés LA, Suárez HA (2005) [Triatomines (Reduviidae: Triatominae) in a Chagas disease focus in Talaigua Nuevo (Bolívar, Colombia)]. Biomedica 25: 568-574. Find this article online

Guhl F, Aguilera G, Pinto N, Vergara D (2007) [Updated geographical distribution and ecoepidemiology of the triatomine fauna (Reduviidae: Triatominae) in Colombia]. Biomedica 27: Suppl 1143-162. Find this article online

Lent H, Wigodzinski P (1979) Revision of the Triatominae (Hemiptera, Reduviidae) and their significance as vector of Chagas diseae. Bulletin of the American Museum of Natural History 163: 123-520. Find this article online

Sambrook J, Fritsch EF, Maniatis T (1989) Molecular cloning: A laboratory manual. Second edition. New York: Cold Spring Harbor Laboratory Press, US

Hall T (2004) BioEdit, Biological sequence alignment editor for Win95/98/NT/2K/XP. Department of microbiology, North Carolina State University. Available: http://www.mbio.ncsu.edu/BioEdit/bioedit.html. Accesed 2012 Jan 30

Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673-4680. doi: 10.1093/nar/22.22.4673. Find this article online

Hall T (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp ser 41: 95-98. Find this article online

Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol 28: 2731-2739. doi: 10.1093/molbev/msr121. Find this article online

Aarskog NK, Vedeler CA (2000) Real-time quantitative polymerase chain reaction. A new method that detects both the peripheral myelin protein 22 duplication in Charcot-Marie-Tooth type 1A disease and the peripheral myelin protein 22 deletion in hereditary neuropathy with liability to pressure palsies. Hum Genet 107: 494-498. Find this article online

Zhang DX, Hewitt GM (1996) Nuclear integrations: challenges for mitochondrial DNA markers. Trends Ecol Evol 11: 247-251. doi: 10.1016/0169-5347(96)10031-8. Find this article online

Bensasson D, Zhang D, Hartl DL, Hewitt GM (2001) Mitochondrial pseudogenes: evolution's misplaced witnesses. Trends Ecol Evol 16: 314-321. doi: 10.1016/S0169-5347(01)02151-6. Find this article online

Liu Y, Zhao X (2007) Distribution of nuclear mitochondrial DNA in cattle nuclear genome. J Anim Breed Genet 124: 264-268. doi: 10.1111/j.1439-0388.2007.00674.x. Find this article online

Alcaide M, Rico C, Ruiz S, Soriguer R, Muñoz J, et al. (2009) Disentangling vector-borne transmission networks: a universal DNA barcoding method to identify vertebrate hosts from arthropod bloodmeals. PLoS One 4: e7092. Find this article online

Pineda V, Saldaña A, Monfante I, Santamaría A, Gottdenker NL, et al. (2011) Prevalence of trypanosome infections in dogs from Chagas disease endemic regions in Panama, Central America. Vet Parasitol 178: 360-363. doi: 10.1016/j.vetpar.2010.12.043. Find this article online

Jiménez-Coello M, Guzmán-Marín E, Ortega-Pacheco A, Acosta-Viana KY (2010) Serological survey of American trypanosomiasis in dogs and their owners from an urban area of Mérida Yucatàn, México. Transbound Emerg Dis 57: 33-36. doi: 10.1111/j.1865-1682.2010.01130.x. Find this article online

Jimenez-Coello M, Ortega-Pacheco A, Guzman-Marin E, Guiris-Andrade DM, Martinez-Figueroa L, et al. (2010) Stray dogs as reservoirs of the zoonotic agents Leptospira interrogans, Trypanosoma cruzi, and Aspergillus spp. in an urban area of Chiapas in southern Mexico. Vector Borne Zoonotic Dis 10: 135-141. Find this article online

Rowland ME, Maloney J, Cohen S, Yabsley MJ, Huang J, et al. (2010) Factors associated with Trypanosoma cruzi exposure among domestic canines in Tennessee. J Parasitol 96: 547-551. doi: 10.1645/GE-2299.1. Find this article online

Dib J, Barnabe C, Tibayrenc M, Triana O (2009) Incrimination of Eratyrus cuspidatus (Stal) in the transmission of Chagas' disease by molecular epidemiology analysis of Trypanosoma cruzi isolates from a geographically restricted area in the north of Colombia. Acta Trop 111: 237-242. doi: 10.1016/j.actatropica.2009.05.004. Find this article online