Establishing the sources of reinfestation after residual insecticide spraying is crucial for vector elimination programs. Triatoma infestans, traditionally considered to be limited to domestic or peridomestic (abbreviated as D/PD) habitats throughout most of its range, is the target of an elimination program that has achieved limited success in the Gran Chaco region in South America.
During a two-year period we conducted semi-annual searches for triatomine bugs in every D/PD site and surrounding sylvatic habitats after full-coverage spraying of pyrethroid insecticides of all houses in a well-defined rural area in northwestern Argentina. We found six low-density sylvatic foci with 24 T. infestans in fallen or standing trees located 110–2,300 m from the nearest house or infested D/PD site detected after insecticide spraying, when house infestations were rare. Analysis of two mitochondrial gene fragments of 20 sylvatic specimens confirmed their species identity as T. infestans and showed that their composite haplotypes were the same as or closely related to D/PD haplotypes. Population studies with 10 polymorphic microsatellite loci and wing geometric morphometry consistently indicated the occurrence of unrestricted gene flow between local D/PD and sylvatic populations. Mitochondrial DNA and microsatellite sibship analyses in the most abundant sylvatic colony revealed descendents from five different females. Spatial analysis showed a significant association between two sylvatic foci and the nearest D/PD bug population found before insecticide spraying.
Our study shows that, despite of its high degree of domesticity, T. infestans has sylvatic colonies with normal chromatic characters (not melanic morphs) highly connected to D/PD conspecifics in the Argentinean Chaco. Sylvatic habitats may provide a transient or permanent refuge after control interventions, and function as sources for D/PD reinfestation. The occurrence of sylvatic foci of T. infestans in the Gran Chaco may pose additional threats to ongoing vector elimination efforts.
Triatoma infestans, a highly domesticated species and historically the main vector of Trypanosoma cruzi, is the target of an insecticide-based elimination program in the southern cone countries of South America since 1991. Only limited success has been achieved in the Gran Chaco region due to repeated reinfestations. We conducted full-coverage spraying of pyrethroid insecticides of all houses in a well-defined rural area in northwestern Argentina, followed by intense monitoring of house reinfestation and searches for triatomine bugs in sylvatic habitats during the next two years, to establish the putative sources of new bug colonies. We found low-density sylvatic foci of T. infestans in trees located within the species' flight range from the nearest infested house detected before control interventions. Using multiple methods (fine-resolution satellite imagery, geographic information systems, spatial statistics, genetic markers and wing geometric morphometry), we corroborated the species identity of the sylvatic bugs as T. infestans and found they were indistinguishable from or closely related to local domestic or peridomestic bug populations. Two sylvatic foci were spatially associated to the nearest peridomestic bug populations found before interventions. Sylvatic habitats harbor hidden foci of T. infestans that may represent a threat to vector suppression attempts.
Citation: Ceballos LA, Piccinali RV, Marcet PL, Vazquez-Prokopec GM, Cardinal MV, Schachter-Broide J, et al. (2011) Hidden Sylvatic Foci of the Main Vector of Chagas Disease Triatoma infestans: Threats to the Vector Elimination Campaign? PLoS Negl Trop Dis 5(10): e1365. doi:10.1371/journal.pntd.0001365
Editor: Jorge A. Huete-Pérez, Universidad Centroamericana, Nicaragua
Received: May 20, 2011; Accepted: September 5, 2011; Published: October 25, 2011
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Funding: This work was supported by the National Institutes of Health/National Science Foundation Ecology of Disease program through NIH Research Grant R01 TW05836 funded by the Fogarty International Center and the National Institute of Environmental Health Sciences (to U.K., R.E.G. and Joel E. Cohen), Agencia Nacional de Promoción Científica y Técnica (Argentina) and by the University of Buenos Aires (R.E.G.). The participation of R.E.G. was also supported by UNDP/World Bank/WHO/TDR (Grant No. A70596) and International Development Research Centre (Grant No. 103696-009). 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.
Disease eradication or elimination programs depend on time-limited intensive campaigns and are likely to fail if resistance to insecticides or drugs (i.e., malaria) or sylvatic transmission cycles (i.e., yellow fever) occur. Chagas disease is the most important vector-borne disease in Latin America in terms of disability-adjusted lost years, with an estimated 10–18 million people infected with Trypanosoma cruzi . Elimination of domestic or peridomestic (hereafter abbreviated D/PD) populations of the insect vectors of T. cruzi through residual spraying with insecticides has shown varying degrees of success depending on the species and the occurrence of sylvatic foci. Several vector species occupy sylvatic habitats and show different degrees of domestication, such as T. dimidiata in Central America, Panstrongylus megistus, T. brasiliensis and T. pseudomaculata in Brazil, Rhodnius ecuadoriensis in northern Peru and Ecuador, and T. pallidipennis and related species in Mexico –. Species of sylvatic or peridomestic triatomines that were not recognized as control targets have emerged as primary vectors of T. cruzi in geographically defined areas over the last two decades [e.g., 5]. For species such as R. prolixus, house reinfestations may also be driven by invasion from peridomestic or sylvatic foci .
Triatoma infestans historically is the main vector of human T. cruzi infection. In 1991, this species was the target of a regional elimination program (the Southern Cone Initiative) that interrupted vector- and blood-borne transmission to humans in Chile, Uruguay, Brazil, eastern Paraguay and parts of Argentina . However, only limited success in the elimination of T. infestans and interruption of vector-borne transmission has been achieved in the Gran Chaco region due to repeated reinfestations even in areas under intensive professional vector control . The Gran Chaco, an ecoregion of 1.3 million km2 mainly spanning northern Argentina, Bolivia and Paraguay, has high levels of poverty and is hyperendemic for Chagas disease . Recurrent reinfestation after residual spraying with insecticides and lack of a sustainable vector surveillance program result in renewed parasite transmission 3–5 years after community-wide vector control campaigns –. The obstacles to the elimination of T. infestans in the Gran Chaco may stem from different processes yet to be identified conclusively.
The Southern Cone Initiative for the elimination of T. infestans was based on two major assumptions with wide consensus and limited supporting evidence , : (i) the species was restricted to D/PD habitats –, with true sylvatic foci only occurring in rock piles associated with wild guinea pigs in the Cochabamba and Sucre Andean valleys in Bolivia –, and (ii) T. infestans had low genetic variability and therefore was very unlikely to develop resistance to modern pyrethroid insecticides. Rare findings of T. infestans in sylvatic habitats up to the early 1980 s were judged to be of little relevance by several investigators (reviewed in , ). The surprising finding of melanic forms (“dark morphs”) in isolated dry forests in the Bolivian ,  and Argentine Chaco , and more recently in the Paraguayan Chaco , combined with the discovery of sylvatic foci with normal phenotypes in Chile and Bolivia – challenged the highly domesticated status of T. infestans. In addition, recent evidence showed T. infestans had richer genetic variability than previously assumed –, with strong chromosomal and DNA content differences between T. infestans from different sources , whereas pyrethroid resistance emerged in northwestern Argentina and throughout Bolivia since the late 1990 s , . Understanding the ecological dynamics of reinfestation in insecticide-treated villages and untangling the mechanisms underlying the observed patterns is crucial for devising improved vector control tactics and the eventual elimination of T. infestans and other major triatomine vectors , . Genetic  and phenetic , markers combined with carefully georeferenced bug samples collected before and after control interventions, a geographic information system (GIS) and spatial statistics  provide the means to better understand reinfestation dynamics. Here we first integrate the use of all these tools to investigate house reinfestation dynamics in the context of control interventions.
As part of a longitudinal project on the eco-epidemiology and control of Chagas disease in a well-defined rural area in the dry Argentine Chaco , we detected isolated findings of adult T. infestans and recently established, very low-density D/PD colonies during two years after a community-wide residual spraying of pyrethroid insecticides of all houses. To identify the putative sources for such occurrences and the sylvatic vectors of T. cruzi , we conducted intensive surveys for triatomine bugs in diverse sylvatic habitats after interventions and surprisingly found various sylvatic foci of T. infestans. Using fine-resolution satellite imagery, GIS, spatial statistics, genetic markers and wing geometric morphometry, we investigated the relatedness between sylvatic and D/PD populations of T. infestans and the threat that they may represent to vector control and elimination attempts in the Argentinean Chaco. Based on previous findings of sylvatic T. infestans in the Bolivian Chaco  and of an isolated adult specimen of T. infestans infected with T. cruzi in semi-sylvatic habitats of our study area in the mid-1980 s , we speculated that similar foci might exist in the Argentinean Chaco and that Triatoma guasayana was a likely candidate sylvatic vector of T. cruzi given its high abundance, widespread occurrence and occasional infection with the parasite , , .
Materials and Methods
Field studies were carried out in Amamá (27° 12′ 30″S, 63° 02′ 30″W) and neighboring rural villages in a 650 km2 area situated in the Moreno Department, Province of Santiago del Estero, Argentina (Figure 1). This area is located in the dry Chaco ecoregion  and its history of infestation since the mid-1980 s has been described elsewhere . Based on the history of control interventions, the study area was subdivided into core (5 villages, 143 domiciles and 790 peridomestic sites) and peripheral (7 villages, 132 houses and 709 peridomestic sites) areas with all sites georeferenced. In April 2004, community-wide residual spraying with 2.5% deltamethrin (K-Othrin, Bayer) of nearly all houses was conducted by professional vector-control personnel using a standard insecticide dose in domiciles (25 mg/m2) and standard or double dose in peridomestic sites for enhanced impact. Here we only report results from the core area (villages of Amamá, Trinidad, Mercedes, Villa Matilde and Pampa Pozo; Figure 1) because no systematic searches for bugs were performed in sylvatic habitats around the peripheral communities.
Red triangles indicate the position of T. infestans-positive mouse-baited traps. Inset shows the location of the study area (black square) within the Gran Chaco region.
Timed manual searches for triatomine bugs with a dislodging spray (0.2% tetramethrin, Espacial) were conducted in all domestic (0.5 person-hour) and peridomestic sites (one person-hour per house compound) from all study villages in October 2004, April and December 2005, and November 2006 as described before . In the core area, 143 domiciles and 764 peridomestic sites were inspected for triatomine bugs at least once between 2004 and 2006. All detected foci were immediately sprayed with deltamethrin using the same procedures. As part of an ongoing monitoring program, discriminant dose assays demonstrated that no pyrethroid resistance occurred in local populations of T. infestans (María Inés Picollo, unpublished results).
We conducted four intensive surveys of triatomine bugs in sylvatic habitats using mouse-baited (Noireau) traps fitted with adhesive tape (Plasto®, Brazil)  in October and November 2005, April and November-December 2006 as described before . Mean temperatures varied between 23°C and 26°C in October-December (spring) surveys, and were below 20°C in April (fall). Searches for sylvatic triatomine foci were conducted in 15 sampling areas that included representative forest sections with different degrees of disturbance (i.e., degraded forest under logging operations, cleared sections, ecotones, and implanted grasslands preceded by selective deforestation) and in all sorts of refuges potentially suitable for triatomine bugs. The total capture effort was 598 trap-nights (range per survey, 129 to 169). Traps were usually placed far from houses in holes of fallen or standing trees (live or dead), trunks or tree stumps and in between terrestrial bromeliads (Bromelia serra and Bromelia hieronymi), cacti (Opuntia quimilo and Opuntia ficus-indica) or piles of shrubs (Figure S1). Traps were deployed when the weather was warm and not rainy approximately between 17.00–18.00 hs and retrieved before 10.00 hs to protect mice from exposure to extreme temperatures. All trap locations were georeferenced using a GPS (Garmin, Etrex Legend C). All sylvatic sites surveyed in October and November 2005 were different except one, and 98% of them were re-inspected with mouse-baited traps on April 2006 to assess bug occurrence, persistence and invasion. The survey conducted in November-December 2006 only included sites that had not been surveyed previously.
Flight-dispersing triatomine bugs were collected using black-light traps  placed in 36 georeferenced sylvatic sites where concurrent searches with mouse-baited traps were made (i.e., in the same areas). Light traps were deployed away from houses in habitats where there was a wide opening in the forest that allowed at least a 100 m visibility. Light traps were operated from approximately 19:45 (i.e., 15 min before sunset) to 22:00–23:00 hs because the flight activity of T. infestans peaks during the first hour after sunset, and is more likely to occur when air temperature exceeds 20°C and wind speed is <5 km/h –. Suitable conditions for flight initiation of T. infestans occurred during the surveys conducted in October-November 2005 but not in April 2006.
All collected bugs were kept alive in plastic vials with folded filter paper, identified to species following Lent and Wygodzinsky  and counted. Species identification of very small first- or second-instar nymphs sometimes was considered tentative depending on the integrity of the material; no such doubts remained for third-instars or later stages. Fourth- or fifth-instar nymphs and adult bugs collected in 2005 were individually weighed on an electronic balance (OHAUS, precision, 0.1 mg) and total body length (L) measured from the end of the clipeus to the end of the abdomen with a vernier caliber (precision, 0.02 mm) to estimate a weight-to-length ratio (W∶L) –a quantitative index of nutritional status. The qualitative nutritional status of nymphs was determined by a cross-sectional view of the abdomen and cuticle distension and classified into four categories that ranged from unfed to large blood contents . Feces from live third-instars and larger stages were examined microscopically for T. cruzi infection at 400× magnification.
DNA from bugs assigned to T. infestans (based on morphological evidence) was obtained, PCR-amplified, and sequenced for a 661 bp fragment of the mitochondrial genes cytochrome oxidase I (mtCOI)  and a 572 bp fragment of the cytochrome B (mtcytB) gene . Sequences from sylvatic T. infestans were compared with Triatoma spp sequences available at Genbank and from previous surveys on the instraspecific variability of T. infestans , –.
Sylvatic T. infestans mtCOI plus mtcytB composite haplotypes were compared with previously recorded haplotypes of D/PD T. infestans from the study villages (collected in 2001–2002), from other more distant (40 km) localities within Santiago del Estero Province (Quilumpa, Km 40, La Loma and Invernada Norte, collected in 2003–2004), and from other Argentinean Provinces more than 300 km apart (Salta, La Rioja, Tucumán and Formosa, collected in 2000–2005). A detailed description of the source localities was published elsewhere . Genetic variability was estimated as the mean number of pairwise differences per site (π), Watterson's estimator (θW) and the haplotype diversity (Hd) with DnaSP 5.0  and a statistical parsimony haplotype network was built with TCS 1.21 .
For higher resolution of the relationships between sylvatic and D/PD populations of T. infestans, the multilocus (ML) genotype for 10 microsatellite loci was obtained for sylvatic T. infestans using primers and PCR conditions previously described . ML genotypes were compared with those from T. infestans captured in D/PD sites from Amamá and neighboring villages in October 2002 and April 2004 before full-coverage insecticide spraying . Inter-individual genetic distance based on the complement of the proportion of shared alleles  was estimated with MICROSAT 1.5d (http://hpgl.stanford.edu/projects/microsat/), and a neighbor-joining (NJ) tree was built with MEGA 3.0 .
Using the genotypes of local D/PD T. infestans as reference populations, we applied the Bayesian based assignment-exclusion test implemented in GENECLASS 2  to individually assign sylvatic individuals to the local pre-spraying D/PD populations (defined as the total gene pool at a given community in each capture date). No post-spraying reference groups could be formed because after community-wide insecticide spraying (2004–2006) most bug collections contained one or a few insects per site that were sparsely distributed throughout the communities (i.e., no established populations of T. infestans were detected). Reference populations were not excluded as the putative origin of the sylvatic insects when the marginal probability exceeded 0.05. We used 100,000 replications and a simulation algorithm .
Sibship of T. infestans bugs collected in traps with more than one individual (TN-92 and TN-139) was inferred with the maximum likelihood approach implemented in COLONY 2.0  performing two independent runs and assuming a probability of null alleles of 0.05 in loci ms42, ms64 and ms65 due to departures from Hardy-Weinberg expectations.
The wing geometric morphometry of the only sylvatic T. infestans male collected was compared with T. infestans males captured in D/PD sites from Amamá and neighboring study villages in October 2002 (n = 87) and April 2004 (n = 74) as described elsewhere . The geometric coordinates of 11 type-I landmarks (venation intersections) from all right wings were digitized by the same user (JSB). After performing the generalized Procrustes superposition (GPA, ), the residual coordinates of the total sample (including the sylvatic specimen) were transformed into partial warps (PW). These shape variables allow standard statistical analyses such as principal component (PCA) or discriminant analyses (DA). To cope with small sample sizes in some villages, the first nine principal components of the PW were used as input for a DA performed on the village samples (excluding the sylvatic specimen). These principal components are also called relative warps (RW). The sylvatic specimen was then used as supplementary data and its position in the morphospace examined in terms of Mahalanobis distances. Digitization, GPA, PCA and DA were performed using the corresponding modules of the CLIC package .
Global positioning system readings from all sampling sites (with mouse-baited and light traps) were integrated into a Geographic Information System (ArcGIS 9.1, ESRI, Redlands, CA, U.S.A.) of the study communities containing a georeferenced satellite image (Ikonos2, Space Imaging Inc., Atlanta, GA, U.S.A.) and the position of all houses and peridomestic sites sprayed with insecticides in 2004. Cartesian coordinates (Universal Transverse Mercator, UTM, Zone 20S) were calculated for each D/PD site and trapping location in order to perform spatial analysis. A focal spatial statistic (Gi(d))  was used to determine the presence and extent of spatial clustering of T. infestans D/PD abundance (average of timed manual catches of bugs per site in 2002 and 2004] around each T. infestans-positive sylvatic focus (point i). This local statistic is additive in the sense that it focuses on the sum of the j values in the vicinity of point i. Hence, we took each T. infestans-positive sylvatic focus, one at a time, and searched the nearby area for occurrences of more or fewer D/PD T. infestans bugs collected before full-coverage insecticide spraying than expected by random. This procedure identified specific trap locations as members or non-members of infestation clusters. We used a binary weight wij based on a distance threshold (d) scheme. Clustering of D/PD T. infestans abundance around a positive sylvatic site occurred when the observed Gi was higher than 2.32 (the expected value at P<0.01). We evaluated the value of Gi up to 3 km from each sylvatic site with T. infestans –a tentative upper bound of the flight range of T. infestans. Analyses were performed using the software Point Pattern Analysis .
Humane care and use of laboratory animals were performed according to Institutional Animal Care and Use Committee (IACUC, CICUAL in Spanish) guidelines at UBA's Faculty of Exact and Natural Sciences. Animal care and use is guided by the International Guiding Principles for Biomedical Research Involving Animals developed by the Council for International Organizations of Medical Sciences.
Collection and nutritional state of triatomine bugs
A total of 13 (9.1% of 143 domiciles) domestic foci of T. infestans with 23 bugs and 38 (5.0% of 764 sites) peridomestic foci with 223 bugs were detected between 2004 and 2006 after full-coverage spraying with deltamethrin. Nearly 25% of all collected T. infestans were adult bugs.
Only 30 (5%) of 598 mouse-baited traps deployed overnight in sylvatic habitats were positive for triatomine bugs (Table 1). Six sylvatic foci of T. infestans with normal chromatic characters (totaling 23 nymphs and 1 male; range per site, 1–17) were found in tree holes or trunks (Figures S1 and S2). Another probable sylvatic foci of this species with two first- or second-instar nymphs was conservatively excluded because the morphological identification of these stages was uncertain and mtDNA markers did not amplify; this probable focus occurred in the vicinity of the largest sylvatic colony of T. infestans (trap TN-139, Figure S2). The apparent density of sylvatic T. infestans was 4 per 100 trap-nights (24 bugs in 598 trap-nights; mean ± SD, 3.8±6.4 bugs per site). One sylvatic focus located west of Amamá (trap TN-139) was infested both in October (1 male) and November 2005 (14 first- or second-instar nymphs and 2 fourth-instars) and was taken as one colony. No T. infestans bugs were collected with mouse-baited traps in April or November 2006.
T. guasayana occurred more frequently (3.0% of mouse-baited traps in all surveyed habitats) than T. infestans (1.2%, Table 1). Feces and hairs of Didelphis opossums were found in one T. guasayana focus. All first- or second-instars of Triatoma sp. not identified to species level most likely were T. guasayana based on morphology, size and type of habitat. Light-trap collections yielded 110 adult T. guasayana, one specimen of T. garciabesi (female) and one of T. platensis (male), and no T. infestans in 41 light-trap-nights (Table 1). Of the 41 light-trap nights, 28 (68.3%) were positive for triatomine bugs. The adult sex ratio in T. guasayana was 1∶2.2 (male to female).
Sylvatic foci of T. infestans occurred at 5 sampling areas located 2.0–11.5 km apart (Figure 2). Most triatomines (17 or 70.8% of 24 T. infestans and 18 or 64.3% of 28 T. guasayana) caught with mouse-baited traps occurred in areas that had been deforested selectively (totalling 40 bugs at 11 sites); the other seven T. infestans were caught in secondary forest with medium-sized or a few large-sized trees. The only three T. garciabesi found were caught in mature forest under active deforestation. The remaining triatomine bugs were caught in secondary forest with medium- or large-sized trees. The main identified micro-habitats of T. infestans were in holes of fallen trees and decaying tree trunks lying on the ground (21 or 87% of 24 bugs collected), a tree stump and a live standing tree. These ecotopes included 4 ‘quebracho colorado’ (Schinopsis lorentzii) and 2 ‘mistol’ (Zizyphus mistol) trees (Figure S1).
Plot of (Gi(d)) values estimated for the peridomestic abundance of T. infestans (2002–2004) as a function of distance to each sylvatic focus within 3 km of Amamá (A) and Mercedes (B). Dotted lines represent 99% confidence intervals.
Nearly all triatomine bugs caught with mouse-baited traps and examined for qualitative nutritional status (n = 36) were unfed (61.1%) or had very little remnants of a blood meal (33.3%) and very low W/L ratios (Table S1). Of 140 sylvatic triatomine bugs examined microscopically (10 T. infestans, 21 T. guasayana and 3 T. garciabesi caught with mouse-baited traps and 106 T. guasayana collected with light traps) none was found microscope-positive for T. cruzi.
mtDNA analyses of T. infestans
The morphological identification of 20 sylvatic bugs as T. infestans was confirmed by DNA sequencing of mtCOI and/or mtcytB fragments; DNA from six other bugs (all first- to third-instars identified as T. infestans based on morphological characters) could not be amplified. The two third-instar nymphs not amplified were taken as T. infestans because a morphological misidentification (relative to the locally known species) was considered very unlikely. None of the sylvatic T. infestans bugs carried the T_C change at position 556, which is characteristic of T. platensis and is absent in a large sample of T. infestans from Argentina, Bolivia, Peru, and Uruguay .
Sylvatic T. infestans with mtCOI and mtcytB composite haplotypes (n = 16, Table S2) exhibited high nucleotide variability (θW = 0.006, π = 0.007) and haplotype diversity (Hd = 0.901). No shared haplotypes were found among bugs from different traps, whereas traps with more than one bug had one (TN-92, n = 3) and five (TN-139, n = 11) different haplotypes (Table S2). Of eight sylvatic haplotypes identified, six were exclusive of sylvatic bugs whereas two haplotypes were recorded in local peridomestic populations of T. infestans and elsewhere in Argentina (Figure 3). Sylvatic haplotypes were spread along the entire statistical parsimony network; they did not form a unique cluster separated from the rest and were more closely related to D/PD than to other sylvatic haplotypes (Figure 3). One sylvatic haplotype was highly divergent (am-XIV) but also was closely connected to an Amamá peridomestic haplotype (haplotype b-XIV).
Each line represents a mutational step and the small empty circles are unobserved haplotypes.
Microsatellite and wing morphometry analyses
The multilocus (ML) genotype for 10 microsatellite loci was obtained for 21 sylvatic T. infestans. We identified a total of 86 different alleles for the 10 loci, of which only 15 (17.5%) and 17 (19.8%) were private alleles not detected in the local D/PD populations in 2002 and 2004, respectively. Sylvatic T. infestans clustered among D/PD bugs with no sharp discontinuity (Figure 4). T. infestans bugs captured concurrently at trap TN-139 clustered together whereas bugs collected there at different times were more closely related to different clusters of Amamá peridomestic bugs (i.e., the closest village). In addition, insects from trap TN-139 had five different mtCOI-mtcytB haplotypes (Table S2). Sibship microsatellite analyses showed that bugs that shared a mitochondrial haplotype (or that had consistent haplotypes because of missing data for mtCOI or mtcytB) were most likely full- or half-sibs whereas bugs with different haplotypes were not (Tables S3 and S4). Bugs from trap TN-92 clustered together and closely to bugs from Mercedes village (where the trap was located) and from another village at ∼5 km (Pampa Pozo). These three bugs were full- or half-sibs and shared the same mitochondrial haplotype (Tables S3 and S4). The bug from site trap TN-182 was grouped with bugs from the nearest village (Mercedes) located at ∼8 km. The bug collected at trap TN-101 (close to Villa Matilde, Fig. 1) clustered with bugs from Amamá and Pampa Pozo.
Comparison of sylvatic with domestic or peridomestic T. infestans populations captured before full-coverage insecticide spraying in 2002.
The Bayesian-based assignment-exclusion test indicated that 18 of 21 sylvatic ML genotypes were not excluded from one or more of the D/PD reference populations (Table 2). D/PD populations were excluded as putative sources for three sylvatic insects captured in two different sites (traps TN-182 and -139). The mtCOI-mtcytB haplotype from the bug in trap TN-182 (al VII, Figure 4) was also genetically distant from the local D/PD populations and was closely related to D/PD populations from La Rioja, more than 400 km far from the study area (Figure 4).
Wing geometric morphometry was used to compare the only sylvatic T. infestans male collected (trap TN-139) with T. infestans males captured in local PD sites in 2002 and 2004. The factorial map showed that the sylvatic bug clearly overlapped with 2002 PD bugs from Amamá –the closest village to its capture site (Figure 5) and it was also assigned to 2004 PD bugs from Amamá (not shown).
All sylvatic foci of T. infestans were located 110–2,300 m from the nearest D/PD sites ever found to be infested by this species after full-coverage insecticide spraying (i.e., detected during the preceding 18 months) (Figure 2). Trap location TN-182 included two T. infestans-positive sites (TN-180 and TN-182) that were analyzed together because their separation (13 m) was smaller than the distance resolution of the Gi(d) test (50 m). The distance between traps positive for T. infestans to the nearest house varied from 125 to 1,900 m. Spatial analysis showed a statistically significant association (Gi(d)>2.32, P<0.01) between two sylvatic foci of T. infestans found within three km of a D/PD site and the average timed-manual catch of bugs before insecticide spraying (Figure 2). Significant clustering occurred up to 1.2 km in Amamá (trap TN-139, with 17 insects) and up to 150 m in Mercedes (trap TN-101, with one third-instar nymph) (Figure 2). The remaining three sylvatic foci of T. infestans (TN-182, TN-180 and TN-92) were located at 430–1,846 m from the nearest infested house, but did not appear to be significantly associated with any of them (Gi(d)>1.96; P>0.05).
We report here the first finding of multiple sylvatic foci of T. infestans: i) with normal chromatic characters (not “dark morphs”) in the Gran Chaco region outside Bolivia; ii) with morphological identification confirmed by DNA sequence information –ruling out taxonomic misdiagnosis of nymphs, and iii) with a genetic makeup indistinguishable from their local D/PD conspecifics in nearly all cases. The discovery of sylvatic foci of T. infestans was made possible by the extensive deployment of mouse-baited sticky traps in a wide diversity of habitats potentially suitable for the species as suggested by surveys in the Bolivian Chaco . Although mouse-baited sticky traps may not achieve perfect detection of sylvatic foci , , the alternatives of using timed manual collections with a dislodging spray or habitat destruction are even less satisfactory or feasible . Easy-to-use, more sensitive sampling methods for triatomine bugs in sylvatic habitats are crucially needed. Therefore, the actual prevalence of sylvatic foci of T. infestans as determined with mouse-baited traps was most likely underestimated.
Mitochondrial and microsatellite DNA markers coupled with wing geometric morphometry consistently indicated the occurrence of unrestricted gene flow between local D/PD and sylvatic T. infestans populations. In spite of the occurrence of private mitochondrial haplotypes and microsatellite alleles in sylvatic bugs, analyses suggest a strong genetic relationship with D/PD bugs. In the phylogenetic network, mtDNA sylvatic haplotypes were more frequently connected to peridomestic haplotypes rather than to other sylvatic variants ―which indicates that they did not form a population that evolved in isolation for a long period of time. A limitation here is that mtDNA only allows the estimation of historical female-based gene flow. However, microsatellite-based analyses ―a more suitable tool for detecting current gene flow― failed to reject that neighboring villages were the putative sources of sylvatic bugs in most cases. In addition, the small level of differentiation between sylvatic and D/PD specimens fell within the observed levels of within-population diversity , . Therefore, there is no sufficient evidence to support restriction of gene flow between sylvatic and D/PD populations of T. infestans from the surrounding villages except in one case (trap TN-182).
Sibship analyses coupled with mitochondrial haplotype information at trap TN-139 over two trapping sessions separated by one month suggest that descendents from five different females were found at this rather remote site. Microsatellite data corroborated the heterogeneous genetic composition of TN-139 bugs as the insects were assigned to different reference populations. In the context of rare, light D/PD infestations after the insecticide spraying campaign, the finding of multiple haplotypes at a defined site was surprising. This sylvatic colony of T. infestans (the largest) was located 1.1 km away from the nearest infested house in an isolated habitat with no signs of current or past human use over the previous two decades (Figure S2). Moreover, another probable sylvatic foci of T. infestans with early-stage nymphs was detected in the vicinity of the largest sylvatic colony. Passive transport of T. infestans in the belongings of rural workers at a transitory camp may have provided an additional means of disseminating bugs within and between communities or the surrounding landscape. This alternative is worth considering because long-distance passive bug transport beyond its distribution range is well known and still occurs . Thus, genetic and morphological evidence combined with the past history of denser D/PD infestations  suggests that the sylvatic specimens of T. infestans may have been feral derivatives (“spill-over”) of D/PD populations. Lack of sampling in sylvatic habitats before full-coverage insecticide spraying unfortunately does not allow establishing whether the sylvatic foci of T. infestans existed before or were formed as a consequence of flight dispersal of D/PD adult bugs or human-assisted passive transport of bugs.
These findings question the widely held notion of an unlikely continuous exchange of T. infestans bugs between wild and domestic habitats in the Chaco. Earlier studies using allozymes or morphometrics ,  and mitochondrial DNA ,  did not detect differences between sylvatic and domestic populations of T. infestans in the Andean Bolivian valleys, neither could mitochondrial markers in Chile . In the allozyme-based study, the findings were interpreted as suggesting that sylvatic foci could be recent derivatives from nearby D/PD bug populations or vice versa –a pattern that was also consistent with unrestricted gene flow between domestic and sylvatic T. cruzi . Intense gene flow between both types of bug populations (abundant at that time) could have generated the same patterns. Using head morphometry in the same study area in the Andean Bolivian valleys, reinfestant specimens of T. infestans found six months after house spraying with pyrethroids were considered survivors of the original domestic bug population unrelated to local sylvatic specimens . Microsatellite data comparing the genetic makeup of sylvatic and D/PD populations of T. infestans showed restricted gene flow between sylvatic and peridomestic populations separated by only 300–650 m at 2,700 m altitude in the Andean Bolivian valleys , whereas they were highly structured and with evidence of low, asymmetric gene flow in a remote, well-preserved dry forest in the Argentinean Chaco . Our data collected in highly-disturbed dry forest with more scattered houses show a different pattern and suggest that the occurrence of sylvatic foci of T. infestans may explain at least some of the new D/PD foci detected after full-coverage residual spraying of insecticides.
All sylvatic foci of T. infestans were 110–2,300 m from the nearest house or infested D/PD site detected after full-coverage insecticide spraying. These distances are within the estimated flight range of this species (1.5 km) derived from direct and indirect observations , , , , . Because T. infestans may sustain tethered flights for >20 min at speeds of 2 m/s , the upper bound of its flight range may reach 3 km and remains uncertain. Therefore, the significant spatial associations detected combined with the range of distances between sylvatic and D/PD foci of T. infestans suggest that these habitats were probably connected through flight dispersal of adult bugs.
The identified habitats of sylvatic T. infestans in our study area were nearly all associated with trees at ground level, in fallen trees or tree stumps. No rocky outcrops were available. Potential bug refuges at higher altitude in the canopy ―difficult to spot and sample― were much less represented in our surveys. Compared with other sylvatic foci investigated with mouse-baited traps, the local apparent density of T. infestans (4 bugs per 100 trap-nights) was slightly higher than that recorded in remote dry forest in the Argentine Chaco (1.2 bugs per 100 trap-nights) , and substantially lower than in the Bolivian Chaco (17 bugs per 100 trap-nights)  or the Andean valleys (8–123 bugs per 100 trap-nights) . The finding of small, malnourished sylvatic colonies with immature stages of T. infestans indicates that despite extensive deforestation and land-use change, the degraded forest still maintained suitable conditions and resources for bug development but at reduced levels: the apparent abundance and availability of blood-meal sources (not identified yet) in local sylvatic habitats was poorer and more unstable than in D/PD habitats. Some of the sylvatic bug foci in our study could be considered “semi-sylvatic”, in the sense that these habitats were intermediate between peridomestic ecotopes (such as pig or goat corrals made with piled thorny shrubs) and sylvatic habitats in terms of resident host species and abundance . Semi-sylvatic habitats also tend to be less used and modified by regular human activities than peridomestic ecotopes. These findings suggest the possibility of sylvatic foci of T. infestans in almost any rural area within its geographic range. The domestication process T. infestans underwent in the past does not prevent the species from surviving at low density in a wide diversity of sylvatic or semi-sylvatic habitats, even after community-wide insecticide spraying.
Unlike previous reports in Argentina where the presence of sylvatic T. infestans may be a result of spill over from heavy D/PD infestations [reviewed in 24], the sylvatic T. infestans of this study occurred in sampling areas around villages under vector surveillance and selective control activities that only allowed the establishment of very few low density D/DP foci for limited time periods between surveys. A relevant question is whether the small-sized sylvatic bug populations we found are viable in the absence of immigration from D/PD sources (i.e., ‘rescue effects’) or they simply are temporary sinks. Our second-year follow-up data raise doubts about their viability over a longer time horizon in the absence of immigration, although removal of bugs may have contributed to apparent local extinctions. However, it is noteworthy that “dark morph” populations of T. infestans in the Bolivian and Argentine Chaco were viable despite having very low density and remote locations excluding them from D/PD ‘rescue effects’ , , .
T. guasayana was far more abundant than T. infestans in sylvatic habitats, and light-trap collections demonstrated the large number of flight-dispersing adult T. guasayana, as was found in the Bolivian and Paraguayan Chaco , . Previous studies showed that T. guasayana also colonized peridomestic structures and semi-sylvatic ecotopes where it was associated positively with the local abundance of goats and the density of cacti and bromeliads . Householders frequently collected adult bugs of this species when invading human habitations at sunset but this species was not able to colonize domestic premises before or after apparent suppression of T. infestans , , . In the present study, the concurrent finding of T. guasayana in a fallen tree with fresh signs of Didelphis opossums suggests a close association with the main local sylvatic reservoir of T. cruzi typically infected with discrete typing unit I , . The widespread occurrence and large abundance of T. guasayana combined with its ocassional infection, opportunistic blood-feeding behavior and dispersal ability implicate it as a secondary vector of T. cruzi in the peridomestic environment  and sylvatic habitats in the Argentine Chaco.
Implications for vector control and elimination
A long-standing, key scientific question with vast implications for vector control is what is the source of the triatomine bugs appearing after community-wide insecticide spraying , , . Are they (i) survivors or the offspring of previously existing bugs; (ii) immigrants from untreated D/PD or sylvatic foci; or (iii) migrants brought by passive transport from other villages or elsewhere? This issue is applicable to all major triatomine vector control programs throughout Latin America and the responses may differ between settings and even within the same species, as with T. dimidiata in Central America and Mexico or T. brasiliensis and P. megistus in Brazil –all of which display substantial within-species differences in habitat distribution, invasive capacity and other relevant traits. As with other species of triatomine bugs, T. infestans adults and nymphs are attracted to lights , . Sylvatic populations of T. infestans are much more widespread than assumed in the past – and have recently been discovered in the Paraguayan Chaco . Because sylvatic habitats are not targeted for vector control operations, they may provide hidden refuges for T. infestans from which they may reinvade houses in search of more suitable conditions and resources. Our results suggest that in areas with recurrent reinfestation, vector control programs should consider the potential occurrence of external sources (semi-sylvatic or sylvatic) around the target community. The role that sylvatic populations of T. infestans (either with melanic or normal phenotype) play in the process of recolonization of insecticide-treated villages and their invasive capacity needs to be more widely investigated to evaluate the risk they pose to effective vector control and eventual elimination in the Gran Chaco and elsewhere.
Ecotopes where sylvatic foci of T. infestans were searched for and eventually detected. A) holes of standing trees, B) dry cacti (Opuntia quimilo and Opuntia ficus-indica), C) terrestrial bromeliads (Bromelia serra and Bromelia hieronymi), D) piles of shrubs, E) tree trunks or stumps, F) holes of fallen trees.
Two of the sylvatic foci where T. infestans was detected. A) TN-139. B) TN-182.
Weight-to-length ratios of T. infestans and T. guasayana. Mean for adult bugs, medians for nymphs, minimum and maximum values are reported according to collection site. Data for peridomestic bugs collected in Amamá, Trinidad and Mercedes (October 2000–August 2001) were taken from Ceballos et al. 2005 and L. A. Ceballos, unpublished data).
Mitochondrial haplotypes and microsatellite genotypes of sylvatic T. infestans. NA: no PCR amplification. mtCOI Genbank accession numbers: EF451012-4, FJ811845, GQ478993, GQ478995, GQ478993. mtcytB Genbank accession numbers: AY062165, JN006793-9.
Sibship maximum likelihood analyses in traps TN-92 and TN-139. Only values with probabilities greater than 0.6 are shown.
Reconstructed full- and half-sib families in TN-92 and TN-139. A question mark means an unknown haplotype.
Translation of the abstract into Spanish by author Ricardo E. Gürtler.
In memory of François Noireau. For support and advice we thank Gustavo Azzimonti, M. Carla Cecere, Juan M. Gurevitz, Jonathan B. Miranda, Walter R. Escalada, “Uchi” Escalada, and Laura Tomassone. R. E. G. and J. P. D. thank the ECLAT network for helpful comments. R. E. G., R. V. P. and M. V. C. are members of CONICET Researcher's Career, Argentina. CDC core facilities provided the microsatellite and mtcytB oligonucleotides. The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the views of the Centers for Disease Control and Prevention.
Conceived and designed the experiments: REG UK JPD EMD RVP LAC. Performed the experiments: LAC GMVP PLM MVC RVP JSB. Analyzed the data: LAC RVP PLM GMVP MVC JSB. Contributed reagents/materials/analysis tools: LAC GMVP PLM RVP JSB EMD JPD UK REG. Wrote the paper: REG LAC GMVP PLM RVP JPD EMD UK.
- 1. World Health Organization (2004) World Health Report 2004: Changing history (World Health Organization, Geneva). Available: www.who.int/whr/2004/en/report04_en.pdf. Accessed 2011 April 11.
- 2. Diotaiuti L, Faria Filho OF, Carneiro FCF, Dias JCP, Pires HHR, et al. (2000) Aspectos operacionais do controle do Triatoma brasiliensis. Cad Saúde Públ 16: 7–14.
- 3. Ramsey JM, Cruz-Celis A, Salgado L, Espinosa L, Ordonez R, et al. (2003) Efficacy of pyrethroid insecticides against domestic and peridomestic populations of Triatoma pallidipennis and Triatoma barberi (Reduviidae: Triatominae) vectors of Chagas' disease in Mexico. J Med Entomol 40: 912–920.
- 4. Dumonteil E, Ruiz-Pina H, Rodriguez-Felix E, Barrera-Perez M, Ramirez-Sierra MJ, et al. (2004) Re-infestation of houses by Triatoma dimidiata after intradomicile insecticide application in the Yucatan Peninsula, Mexico. Mem Inst Oswaldo Cruz 99: 253–266.
- 5. Matias A, de la Riva J, Martínez E, Torrez M, Dujardin JP (2003) Domiciliation process of Rhodnius ståli (Hemiptera: Reduviidae) in Alto Beni, La Paz, Bolivia. Trop Med Int Health 8: 264–268.
- 6. Sánchez-Martin MJ, Feliciangeli MD, Campbell-Lendrum D, Davies CR (2006) Could the Chagas disease elimination programme in Venezuela be compromised by reinvasion of houses by sylvatic Rhodnius prolixus bug populations? Trop Med Int Health 11: 1585–1593.
- 7. Silveira AC (2002) O controle da doença de Chagas nos países do Cone Sul da America. História de uma iniciativa internacional. 1991–2001. In: Silveira AC, Rojas de Arias A, Guillén G, Russomando G, Schenone H, et al., editors. El control de la enfermedad de Chagas en los países del Cono Sur de América. pp. 15–43. Historia de una iniciativa internacional. 1991–2001. Facultad de Medicina do Triangulo Mineiro, Uberaba, Brazil.
- 8. Gürtler RE, Kitron U, Cecere MC, Segura EL, Cohen JE (2007) Sustainable vector control and management of Chagas disease in the Gran Chaco, Argentina. Proc Natl Acad Sci USA 104: 16194–16199.
- 9. Gürtler RE (2009) Sustainability of vector control strategies in the Gran Chaco region: current challenges and possible approaches. Mem Inst Oswaldo Cruz 104: Suppl 152–59.
- 10. Gürtler RE, Cecere MC, Canale D, Castañera MB, Chuit R, et al. (1999) Monitoring house reinfestation by vectors of Chagas disease: a comparative trial of detection methods during a four-year follow-up. Acta Trop 72: 213–234.
- 11. Gürtler RE, Cecere MC, Lauricella MA, Petersen RM, Canale D, et al. (2005) Incidence of Trypanosoma cruzi infection among children following domestic reinfestation after insecticide spraying in rural northwestern Argentina. Am J Trop Med Hyg 73: 95–103.
- 12. Cecere MC, Vazquez-Prokopec GM, Gürtler RE, Kitron U (2006) Reinfestation sources for Chagas disease vector, Triatoma infestans, Argentina. Emerg Infect Dis 12: 1096–1102.
- 13. Gorla DE, Porcasi X, Hrellac H, Catalá SS (2009) Spatial stratification of house infestation by Triatoma infestans in La Rioja, Argentina. Am J Trop Med Hyg 80: 405–409.
- 14. Schmunis GA, Zicker F, Moncayo A (1996) Interruption of Chagas' disease transmission through vector elimination. Lancet 348: 1171.
- 15. Schofield CJ, Dias JC (1999) The Southern Cone Initiative against Chagas disease. Adv Parasit 42: 1–27.
- 16. Abalos JW, Wygodzinsky P (1951) Las Triatominae Argentinas (Reduviidae, Hemiptera). Tucumán: Universidad Nacional de Tucumán. Publicación 601. Monografía N° 2: 179.
- 17. Usinger RL, Wygodzinsky P, Ryckman RE (1966) The biosystematics of Triatominae. Annu Rev Entomol 11: 309–330.
- 18. Schofield CJ (1988) Biosystematics of the Triatominae. In: Service MW, editor. Biosystematics of Haematophagous Insects. Oxford: Clarendon Press. pp. 284–312.
- 19. Lent H, Wygodzinsky P (1979) Revision of the Triatominae (Hemiptera: Reduviidae) and their significance as vectors of Chagas' disease. Bull Am Mus Nat Hist 163: 185–199.
- 20. Dujardin JP, Tibayrenc M, Venegas E, Maldonado L, Desjeux P, et al. (1987) Isozyme evidence of lack of speciation between wild and domestic Triatoma infestans (Heteroptera: Reduviidae) in Bolivia. J Med Entomol 24: 40–45.
- 21. Dujardin JP, Bermudez H, Schofield CJ (1997) The use of morphometrics in entomological surveillance of silvatic foci of Triatoma infestans in Bolivia. Acta Trop 66: 145–153.
- 22. Dujardin JP, Schofield CJ, Tibayrenc M (1998) Population structure of Andean Triatoma infestans: an isozymic study and its epidemiological relevance. Med Vet Entomol 12: 20–29.
- 23. Noireau F, Cortez MR, Monteiro FA, Jansen AM, Torrico F (2005) Can wild Triatoma infestans foci in Bolivia jeopardize Chagas disease control efforts? Trends Parasitol 21: 7–10.
- 24. Ceballos LA, Piccinali RV, Berkunsky I, Kitron U, Gürtler RE (2009) First finding of melanic sylvatic Triatoma infestans (Hemiptera: Reduviidae) in the Argentine Chaco. J Med Entomol 46: 1195–1202.
- 25. Noireau F, Flores R, Gutierrez T, Dujardin JP (1997) Detection of sylvatic dark morphs of Triatoma infestans in the Bolivian Chaco. Mem Inst Oswaldo Cruz 92: 583–584.
- 26. Rolón M, Vega MC, Román F, Gómez A, Rojas de Arias A (2011) First report of colonies of sylvatic Triatoma infestans (Hemiptera: Reduviidae) in the Paraguayan Chaco, using a trained dog. PLoS Negl Trop Dis 5: e1026.
- 27. Bacigalupo BA, Segura MJA, Garcia CA, Hidalgo CJ, Galuppo GS, et al. (2006) Primer hallazgo de vectores de la enfermedad de Chagas asociados a matorrales silvestres en la Región Metropolitana, Chile. Rev Méd Chile 134: 1230–1236.
- 28. Bacigalupo A, Torres-Pérez F, Segovia V, García A, Correa JP, et al. (2010) Sylvatic foci of the Chagas disease vector Triatoma infestans in Chile: description of a new focus and challenges for control programs. Mem Inst Oswaldo Cruz 105: 633–641.
- 29. Buitrago R, Waleckx E, Bosseno MF, Zoveda F, Vidaurre P, et al. (2010) First report of widespread wild populations of Triatoma infestans (Reduviidae, Triatominae) in the valleys of La Paz, Bolivia. Am J Trop Med Hyg 82: 574–579.
- 30. Bargues MD, Klisiowicz DR, Panzera F, Noireau F, Marcilla A, et al. (2006) Origin and phylogeography of the Chagas disease main vector Triatoma infestans based on nuclear rDNA sequences and genome size. Infect Genet Evol 6: 46–62.
- 31. Pérez de Rosas AR, Segura EL, García BA (2007) Microsatellite analysis of genetic structure in natural Triatoma infestans (Hemiptera: Reduviidae) populations from Argentina: its implication in assessing the effectiveness of Chagas' disease vector control programmes. Mol Ecol 16: 1401–1412.
- 32. Piccinali RV, Marcet PL, Noireau F, Kitron U, Gürtler RE, et al. (2009) Molecular population genetics and phylogeography of the Chagas disease vector Triatoma infestans in South America. J Med Entomol 46: 796–809.
- 33. Piccinali RV, Marcet PL, Ceballos LA, Kitron U, Gürtler RE, et al. (2011) Genetic variability, phylogenetic relationships and gene flow in Triatoma infestans dark morphs from the Argentinean Chaco. Infect Genet Evol 11: 895–903.
- 34. Panzera F, Dujardin JP, Nicolini P, Caraccio MN, Rose V, et al. (2004) Genomic changes of Chagas disease vector, South America. Emerg Infect Dis 10: 438–446.
- 35. Picollo MI, Vassena C, Santo Orihuela P, Barrios S, Zaidenberg M, et al. (2005) High resistance to pyrethroid insecticides associated with ineffective field treatments in Triatoma infestans (Hemiptera: Reduviidae) from Northern Argentina. J Med Entomol 42: 637–642.
- 36. Lardeux F, Depickère S, Duchon S, Chavez T (2010) Insecticide resistance of Triatoma infestans (Hemiptera, Reduviidae) vector of Chagas disease in Bolivia. Trop Med Int Health 15: 1037–1048.
- 37. Miles MA, Feliciangeli MD, Rojas de Arias A (2003) Science, medicine, and the future: American trypanosomiasis (Chagas' disease) and the role of molecular epidemiology in guiding control strategies. BMJ 326: 1444–1448.
- 38. Fitzpatrick S, Feliciangeli MD, Sanchez-Martin MJ, Monteiro FA, Miles MA (2008) Molecular genetics reveal that silvatic Rhodnius prolixus do colonise rural houses. PLoS Negl Trop Dis 2: e210.
- 39. Dujardin JP, Beard CB, Ryckman R (2007) The relevance of wing geometry in entomological surveillance of Triatominae, vectors of Chagas disease. Infect Genet Evol 7: 161–167.
- 40. Feliciangeli MD, Sanchez-Martin M, Marrero R, Davies C, Dujardin JP (2007) Morphometric evidence for a possible role of Rhodnius prolixus from palm trees in house re-infestation in the State of Barinas (Venezuela). Acta Trop 101: 169–177.
- 41. Kitron U, Clennon JA, Cecere MC, Gürtler RE, King CH, et al. (2006) Upscale or downscale: applications of fine scale remotely sensed data to Chagas disease in Argentina and schistosomiasis in Kenya. Geospatial Health Applic 1: 49–58.
- 42. Ceballos LA, Cardinal MV, Vazquez-Prokopec GM, Lauricella MA, Orozco MM, et al. (2006) Long-term reduction of Trypanosoma cruzi infection in sylvatic mammals following deforestation and sustained vector surveillance in northwestern Argentina. Acta Trop 98: 286–296.
- 43. Noireau F, Flores R, Gutierrez T, Abad-Franch F, Flores E, et al. (2000) Natural ecotopes of Triatoma infestans dark morph and other wild triatomines in the Bolivian Chaco. Trans R Soc Trop Med Hyg 94: 23–27.
- 44. Wisnivesky-Colli C, Schweigmann NJ, Alberti A, Pietrokovsky SM, Conti O, et al. (1992) Sylvatic American trypanosomiasis in Argentina. Trypanosoma cruzi infection in mammals from the Chaco forest in Santiago del Estero. Trans R Soc Trop Med Hyg 86: 38–41.
- 45. Cecere MC, Castañera MB, Canale DM, Chuit R, Gürtler RE (1999) Trypanosoma cruzi infection in Triatoma infestans and other triatomines: long-term effects of a control program in a rural area of northwestern Argentina. Rev Panam Salud Publ 5: 392–399.
- 46. Vazquez-Prokopec GM, Cecere MC, Kitron U, Gürtler RE (2008) Environmental and demographic factors determining the spatial distribution of Triatoma guasayana in peridomestic and semi-sylvatic habitats of rural northwestern Argentina. Med Vet Entomol 22: 273–282.
- 47. Noireau F, Flores R, Vargas F (1999) Trapping sylvatic Triatominae (Reduviidae) in hollow trees. Trans R Soc Trop Med Hyg 93: 13–14.
- 48. Vazquez-Prokopec GM, Ceballos LA, Marcet PL, Cecere MC, Cardinal MV, et al. (2006) Seasonal variations in active dispersal of natural populations of Triatoma infestans in rural north-western Argentina. Med Vet Entomol 20: 1–6.
- 49. Schofield CJ, Lehane MJ, McEwen P, Catala SS, Gorla DE (1992) Dispersive flight by Triatoma infestans under natural climatic conditions in Argentina. Med Vet Entomol 6: 51–56.
- 50. Gurevitz JM, Kitron U, Gürtler RE (2007) Flight muscle dimorphism and heterogeneity in flight initiation of field-collected Triatoma infestans (Hemiptera: Reduviidae). J Med Entomol 44: 186–191.
- 51. Ceballos LA, Vazquez-Prokopec GM, Cecere MC, Marcet PL, Gürtler RE (2005) Feeding rates, nutritional status and flight dispersal potential of peridomestic populations of Triatoma infestans in rural northwestern Argentina. Acta Trop 95: 149–159.
- 52. Monteiro FA, Barrett TV, Fitzpatrick S, Cordon-Rosales C, Feliciangeli D, et al. (2003) Molecular phylogeography of the Amazonian Chagas disease vectors Rhodnius prolixus and R. robustus. Mol Ecol 12: 997–1006.
- 53. García BA, Powell JR (1998) Phylogeny of species of Triatoma (Hemiptera: Reduviidae) based on mitochondrial DNA sequences. J Med Entomol 35: 232–238.
- 54. Monteiro FA, Perez R, Panzera F, Dujardin JP, Galvão C, et al. (1999) Mitochondrial DNA variation of Triatoma infestans populations and its implication on the specific status of T. melanosoma. Mem Inst Oswaldo Cruz 94: 229–238.
- 55. Giordano R, Pizarro Cortez JC, Paulk S, Stevens L (2005) Genetic diversity of Triatoma infestans (Hemiptera: Reduviidae) in Chuquisaca, Bolivia based on the mitochondrial cytochrome b gene. Mem Inst Oswaldo Cruz 100: 753–760.
- 56. Cortez MR, Emperaire L, Piccinali RV, Gürtler RE, Torrico F, et al. (2007) Sylvatic Triatoma infestans (Reduviidae, Triatominae) in the Andean valleys of Bolivia. Acta Trop 102: 47–54.
- 57. Librado P, Rozas J (2009) DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452.
- 58. Clement M, Posada D, Crandall KA (2000) TCS: a computer program to estimate gene genealogies. Mol Ecol 9: 1657–1660.
- 59. Marcet PL, Lehmann T, Groner G, Gürtler RE, Kitron U, et al. (2006) Identification and characterization of microsatellite markers in the Chagas disease vector Triatoma infestans (Heteroptera: Reduviidae). Infect Genet Evol 6: 32–37.
- 60. Marcet PL, Mora MS, Cutrera AP, Jones LJ, Gürtler RE, et al. (2008) Genetic structure of Triatoma infestans populations in rural villages of Santiago del Estero, northern Argentina. Infect Genet Evol 8: 835–846.
- 61. Bowcock AM, Ruiz-Linares A, Tomfohrde J, Minch E, Kidd JR, et al. (1994) High resolution of human evolutionary trees with polymorphic microsatellites. Nature 368: 455–457.
- 62. Kumar S, Tamura K, Nei M (2004) MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinformatics 5: 150–163.
- 63. Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, et al. (2004) GENECLASS2: a software for genetic assignment and first-generation migrant detection. J Hered 95: 536–539.
- 64. Paetkau D, Slade R, Burden M, Estoup A (2004) Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Mol Ecol 13: 55–65.
- 65. Jones O, Wang J (2009) COLONY: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour 10: 551–555.
- 66. Schachter-Broide J, Dujardin JP, Kitron U, Gürtler RE (2004) Spatial structuring of Triatoma infestans (Hemiptera, Reduviidae) populations from northwestern Argentina using wing geometric morphometry. J Med Entomol 41: 643–649.
- 67. Rohlf FJ (1996) Morphometric spaces, shape components and the effects of linear transformations. In Advances in Morphometrics. Proceedings of the 1993 NATO-ASI on Morphometrics. In: Marcus LF, Corti M, Loy A, Naylor G, Slice D, editors. pp. 117–129. New York: Plenum Publ. NATO ASI, ser. A, Life Sciences.
- 68. Dujardin JP, Kaba D, Henry AB (2010) The exchangeability of shape. BMC Research Notes 3: 266.
- 69. Getis A, Ord JK (1996) Local spatial statistics: an overview. In: Longley P, Batty M, editors. Spatial Analysis: Modelling in a GIS Environment. Cambridge: Geoinformation International. pp. 261–277.
- 70. Chen D, Getis A (1998) Point Pattern Analysis. San Diego State University, San Diego, CA. Available: http://www.nku.edu/~longa/cgi-bin/cgi-tcl-examples/generic/ppa/ppa.cgi. Accessed: 2011, April 11.
- 71. Abad-Franch F, Ferraz G, Campos C, Palomeque FS, Grijalva MJ, et al. (2010) Modeling disease vector occurrence when detection is imperfect: Infestation of Amazonian palm trees by triatomine bugs at three spatial scales. PLoS Negl Trop Dis 4: e620.
- 72. Piccinali RV, Canale DM, Sandoval AE, Cardinal MV, Jensen O, et al. (2010) Triatoma infestans bugs in southern Patagonia. Emerg Inf Dis 16: 887–889.
- 73. Torres-Pérez F, Acuna-Retamar M, Cook JA, Bacigalupo A, García A, et al. (2010) Statistical phylogeography of Chagas disease vector Triatoma infestans: Testing biogeographic hypotheses of dispersal. Infect Genet Evol 11: 167–174.
- 74. Richer W, Kengne P, Cortez MR, Perrineau MM, Cohuet A, et al. (2007) Active dispersal by wild Triatoma infestans in the Bolivian Andes. Trop Med Int Health 12: 759–764.
- 75. Lehane MJ, Schofield CJ (1981) Field experiments of dispersive flight by Triatoma infestans. Trans R Soc Trop Med Hyg 75: 399–400.
- 76. Schweigmann N, Vallvé S, Muscio O, Guillini M, Alberti A, et al. (1988) Dispersal flight by Triatoma infestans in an arid area of Argentina. Med Vet Entomol 2: 401–404.
- 77. Ward JP, Baker PS (1982) The tethered flight performance of a laboratory population of Triatoma infestans Klug (Hemiptera: Reduviidae). Bull Entomol Res 72: 17–28.
- 78. Noireau F, Gutierrez T, Flores R, Breniere F, Bosseno MF, et al. (1999b) Ecogenetics of Triatoma sordida and Triatoma guasayana (Hemiptera: reduviidae) in the Bolivian chaco. Mem Inst Oswaldo Cruz 94: 451–457.
- 79. Yeo M, Acosta N, Llewellyn M, Sanchez H, Adamson S, et al. (2005) Origins of Chagas disease: Didelphis species are natural hosts of Trypanosoma cruzi I and armadillos hosts of Trypanosoma cruzi II, including hybrids. Int J Parasitol 35: 225–233.
- 80. 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 3: 1533–1543.
- 81. Romaña CA (1963) Enfermedad de Chagas. López Libreros Editores, Buenos Aires.
- 82. Minoli SA, Lazzari CR (2006) Take-off activity and orientation of triatomines (Heteroptera: Reduviidae) in relation to the presence of artificial lights. Acta Trop 97: 324–330.