# Wild populations of Triatoma infestans are highly connected to intra-peridomestic conspecific populations in the Bolivian Andes

Simone Frédérique Brenière, Renata Salas, Rocio Buitrago, Philippe Brémond, V. Sosa, Marie-France Bosseno, Etienne Waleckx, Stéphanie Depickere, Christian Barnabé

<figure><img src="image_1.png" /></figure>

## To cite this version:

Simone Frédérique Brenière, Renata Salas, Rocio Buitrago, Philippe Brémond, V. Sosa, et al. Wild populations of Triatoma infestans are highly connected to intra-peridomestic conspecific populations in the Bolivian Andes. PLoS ONE, 2013, 8 (11), pp.e80786. 10.1371/journal.pone.0080786. hal-01256052

HAL Id: hal-01256052

https://hal.science/hal-01256052v1

Submitted on 21 Sep 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

# Wild Populations of *Triatoma infestans* Are Highly Connected to Intra-Peridomestic Conspecific Populations in the Bolivian Andes

Simone Frédérique Brenière¹,²*, Renata Salas¹,², Rosio Buitrago¹,², Philippe Brémond¹, Victor Sosa¹,³, Marie-France Bosseno¹,², Etienne Waleckx¹,⁴, Stéphanie Depickère¹,², Christian Barnabé¹,²

1 Health Department, MIVEGEC (Université de Montpellier 1 et 2, CNRS 5290, IRD 224), Maladies Infectieuses et Vecteurs : Ecologie, Génétique, Evolution et Contrôle, Institut de recherche pour le développement (IRD), Montpellier, France, 2 Instituto Nacional de Laboratorios de Salud (INLASA), Laboratorio de Entomología Médica, La Paz, Bolivia, 3 Instituto de Investigaciones Biomédicas y de Interacción Social (IIBISMED), Facultad de Medicina, Universidad Mayor de San Simón, Cochabamba, Bolivia, 4 Laboratorio de Parasitología, Centro de Investigaciones Regionales Dr. Hideyo Noguchi, Mérida, Yucatán, México

## Abstract

*Triatoma infestans*, the major vector of Chagas disease south of the Amazon in South America, has a large distribution of wild populations, contrary to what has previously been stated. These populations have been suspected of being the source of reinfestation of human habitats and could impede the full success of vector control campaigns. This study examined gene flow between intra-peridomestic populations and wild populations collected in the surround areas in three Andean localities in Bolivia. The populations were defined according to temporal, ecological, and spatial criteria. After DNA extraction from the legs of each insect, the samples were analyzed using seven microsatellite markers. First, the analysis of molecular variance (AMOVA) detected an absence of differentiation between wild and intra-peridomestic populations, although strong structuring was observed between the populations within each environment. Then for some populations, the Bayesian method of assignment to inferred populations showed very similar assignment patterns of the members of wild or intra-peridomestic populations in each locality. Finally, the detection of the first-generation migrants within the different populations provided evidence of insect displacement from the wild to the intra-peridomestic environment. This result indicates that, after control campaigns in the Andes, controlling this new paradigm of vector transmission risk stemming from the invasion of human habitats by wild populations of *T. infestans* requires long-term maintenance of public monitoring to keep the risk at a minimal level. Since wild populations of *T. infestans* have also been detected elsewhere in Argentina, Paraguay, and Chile, there is an urgent need to take these populations into account in future monitoring of Chagas disease transmission.

**Citation:** Brenière SF, Salas R, Buitrago R, Brémond P, Sosa V, et al. (2013) Wild Populations of *Triatoma infestans* Are Highly Connected to Intra-Peridomestic Conspecific Populations in the Bolivian Andes. PLoS ONE 8(11): e80786. doi:10.1371/journal.pone.0080786

**Editor:** Pedro Lagerblad Oliveira, Universidade Federal do Rio de Janeiro, Brazil

**Received** July 17, 2013; **Accepted** October 7, 2013; **Published** November 20, 2013

**Copyright:** © 2013 Brenière 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:** This work was funded by the National Agency for Research (ANR No 3624, France), the IRD Institute, and the European Union Seventh Program Grant 223034 (ChagasEpiNet). 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.

* E-mail: Frederique.Brenière@ird.fr

## Introduction

In the Southern Cone countries of South America and in Bolivia, *Triatoma infestans* (Reduviidae, Triatominae) remains the main and most widespread vector of *Trypanosoma cruzi*, the causative agent of Chagas disease. This species has long been considered almost exclusively domestic, but the increase in wild population discoveries in the Andes, in the Gran Chaco lowland region, and in Chile shows that wild *T. infestans* populations are much more widespread than previously thought [1-7]; the potential map of the distribution of wild *T. infestans*, based on an earlier sample of Bolivian and Argentinean wild foci, supports a very wide distribution [8].

In 1991, the regional program (INCOSUR-Chagas) of the Southern Cone countries was initiated, aiming to control Chagas disease and eliminate *T. infestans* based on house spraying with residual insecticides. This program has led to the interruption of vector-borne transmissions to humans in Chile, Uruguay, Brazil, eastern Paraguay, and parts of Argentina [9,10]. However, *T. infestans* persists primarily in the Bolivian Andean valleys and in different locations of the Gran Chaco ecoregion. Interestingly, these occurrences match the distribution where wild populations have been found. Determining whether the success in eliminating vectors could be disrupted by repeated reinfestations of dwellings and peridomestic areas by the nearby wild populations remains an
open question. In this case, the current vigilance system, based on selective insecticide spraying of dwellings where vectors are resurgent, remains of limited efficacy and new control strategies must be developed [4,11-13].

Studies in the Gran Chaco region have highlighted that the main source of reinfestation was from habitats within a village rather than from nearby villages [14-17], emphasizing the role of reemerging peridomestic populations as the primarily source of dispersion [18]. Nevertheless, studies of populations at capture sites within villages, using different genetic markers (allozymes, microsatellites, etc.), depicted a high degree of microgeographical genetic structure among populations in the Andean as well as Gran Chaco regions [16,19,20]; a substantial differentiation between populations, even separated by short distances (around 10 m), has been observed. In conclusion, the genetic structure of *T. infestans* in human habitats may be based on infrequent events of bug displacements between ecotopes, giving rise to reduced genetic flow, population isolation, and genetic drift at the local level (colony). Since the assumption of movement of *T. infestans* populations between the wilderness environment and human habitats is recent, few studies are available. Given that the reinfestation phenomenon has been recurring in the Argentinean Chaco for several years, researchers have strengthened the search for wild populations and discovered their existence. In order to assess the origin of reinfesting bugs, they have used several methods, including the analysis of genetic markers (microsatellites and mitochondrial genes), and for the first time assessed reinvasion of domestic and peridomestic structures by nearby wild *T. infestans* colonies [5].

As part of the "TiBo" project founded by the ANR (French Research National Agency), aimed to depict the epidemiological role of the wild *T. infestans* population in Bolivia, unexpectedly, several new sites of wild populations were discovered in the Inter-Andean Dry Forest and Prepuna Andean ecoregions [21]. These bugs exhibited very high infection rates by *T. cruzi*, up to 85.7% in adults captured in the Andean valleys of La Paz Department [4].

Insofar as (i) *T. infestans* is one of the triatomine species that is best adapted to human habitats, (ii) it is an efficient vector, (iii) there is no genetic argument that has questioned the specific or subspecific status of the wild populations, and (iv) their infection rate can be extremely high, the invasion of human habitats by wild Bolivian populations is an extreme danger, potentially leading to resurgence of transmission in areas where domestic populations of triatomines have been controlled.

To establish the population source of reinfesting bugs (residual specimens after the chemical spray of the houses or migrants of domestic or wild origin) observed in three Bolivian Andean areas (two rural and one urban), populations were defined at a small geographical scale (capture site level) according to their habitat (indoors, peridomestic structures, and natural structures found in agropastoral or wild areas), and genetic variation was examined at seven microsatellite DNA loci. Briefly, the results support the displacement of bugs from wild to intra-peridomestic ecotopes as well as displacements between sites with the same type of ecotope.

## Materials and Methods

The field work was realized in collaboration and coordination with the SEDES (Health Departmental Service) of the Bolivian La Paz and Cochabamba Departments, and with the agreement of local authorities in each location (see material and methods). Search for triatomines in private spaces was done with the consent of every inhabitant who accompanied the active search for triatomines by the team members. The field studies did not involve endangered or protected species.

## Geographical origin of *T. infestans*

*T. infestans* were collected in three areas (localities) located in the Eastern Andean Cordillera and the individuals were distributed into wild, peridomestic, and intra-peridomestic subsamples defined according to space and ecotopes (Figure 1). In wild environments, the insects were captured with mice bait traps [22], and in intra-peridomestic sites using a manual active search during the day by the research team. The three areas belong to high valleys of the Inter-Andean Dry Forest ecoregion characterized by mild weather, wet summers, and cool dry winters; most rainfall occurs during summer from December through March. The dry season is substantially colder and lasts from April through August. Depending on the altitude, maximal and minimal temperatures vary slightly.

The first locality is Sapini (16°48'47"S, 67°42'10"W, 1880 m), a small village along the Sapahaqui River in La Paz Department, composed of 30 dwellings where approximately one-third are inhabited or abandoned. At the beginning of the study the team took contact to the authority of the village which has permit to search for triatomines around and in the village. The village is surrounded by high mountains and the lands near the river are used for growing vegetables using an irrigation system. Since 2003 the village has been under entomological vigilance by La Paz Department's Departmental Health Service, visiting the village every 6 months. They found a persistence of dwelling infestation in peridomestic and/or intradomestic sites, ranging from 3.9% to 25% (2003-2010), and was 9.1% at the beginning 2010. The current bugs, collected during October 2010, were grouped into six subsamples, two from intra-peridomestic sites (piles of stones, indoors) and four from wilderness sites: two of them were captured in the cliffs near the village (around 50-100 m away) and two others farther away (300-500 m) in a field of stones and thorns and a pile of stones.

The second area, Quillacollo (17°25'22.03"S 66°17'30.38"W, 2600 m), is located in Cochabamba department at the urban periphery of Quillacollo city, 10 km from Cochabamba city, in the same valley. This suburb was under entomological vigilance by the local municipality authorities, but insecticide spraying of infested/reinfested dwellings was not applied during the study period due to lack of resources from the municipality. The work was developed in collaboration with the authorities of Quillacollo municipality that has provided staff to guide the research team to collect triatomines in the inhabited part. Four populations were collected in peridomestic structures (piles of river stones and piles of bricks where hens often sleep) belonging to different dwellings located at the periphery of the
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Geographical localization of wild and intra-peridomestic populations of *T. infestans*. a – General localization of the three study areas in Bolivia, blue arrows and names. b – Rural locality of Sapini in La Paz Department, inhabited area delimited by blue lines. c – Periurban area of Quillacollo city in Cochabamba Department, blue line delimiting the urban area from the sylvatic area. d – Rural locality of Thago Thago in Potosi Department, blue lines delimiting the grouped dwellings. In a, b, and d, green dots indicate specific sites where the wild *T. infestans* were collected; red dots represent specific sites where intra-peridomestic *T. infestans* were collected; the names of each population in black. See also Table 2 for the ecotope origin.</figcaption></figure>

doi: 10.1371/journal.pone.0080786.g001

urbanized area; four others were in the wilderness in rock outcrops. The periphery of the urban area is surrounded by hills where there are rock outcrops always infested with *T. infestans*. All these populations were collected from October 2009 to June 2010.

The third area is a very small village, Thago Thago (18° 0'41.53"S 65°48'32.68"W, 2000 m), about 80 km southeast of Quillacollo in Potosi Department. As in the village of Sapini, the study was done in accordance with local authorities and also the community participation. The village has ten dwellings along a small road bordered by two hills. The dwellings were grouped (no more than 200 m apart), except for two that were set apart from the others, about 300 and 760 m away. The

dwellings were treated at least once with insecticide before the study, but the village was not monitored by the health authorities during 2010. Seven populations were collected in peridomestic and/or indoor structures belonging to nine different dwellings, and four in the wilderness in rocky areas in the two hills on either side of the group of dwellings.

**Microsatellite genotyping**

DNA was extracted from two bug legs using the slightly modified CTAB-chloroform method [23]. Briefly, legs were crushed in 200 µl of 2% CTAB (cetyl trimethyl ammonium bromide) and incubated overnight at 37°C after the addition of
30 µl of proteinase K (20 mg/ml). The extraction was then performed as described previously [23]. The final DNA pellet was suspended in 20 µl of distilled water. Seven sets of primers identified for *T. infestans* microsatellite loci (TiA02, TiE02, TiC09, TiF11, TiF03, TiD09, TiE12 [24] were used in the polymerase chain reaction (PCR). The forward primer of each set was the 5' end labeled with one of the three fluorescent dyes, FAM, NED, or VIC (Applied Biosystems, Villebon sur Yvette, France). PCR amplifications were performed in a 25-µl solution containing the following products at the concentrations of 1.5 mM for MgCl2, 200 µM for each dNTP, 0.4 µM for each primer (forward-labeled), 1 U of Taq polymerase (Quiagen, Courteboeuf, France), 1X Taq Polymerase buffer, and 20-25 ng template DNA. Amplifications were carried out in a Mastercycler (Eppendorf, Hamburg, Germany) under the following conditions: a starting step of 94°C for 5 min, 35 cycles at 94°C for 30 s, 50°C for 30 s, and 72°C for 30 s, and a final step at 72°C for 15 min. Amplified products were visualized after electrophoresis on a 2.5% agarose gel. Electrophoreses of fluorescent labeled PCR products, diluted and denatured in 20 µl of HiDi formamide, were carried out on an ABI PRISM310 Genetic Analyzer (Applied Biosystems, Carlsbad, CA, USA), with Genescan 500 LIZ as the internal size standard. GeneMapper® software (Applied Biosystems, Carlsbad, CA, USA) was used to characterize alleles.

## Data analysis

The significance of genotype association between pairs of loci (linkage disequilibrium) was tested using Fstat 2.9.3.2 software [25]. The *P*-values of the test for all pairs of loci was obtained as follows: genotypes at two loci were associated at random 2100 times, and the log-likelihood G-test statistic was recalculated on the randomized data. The *P*-value was estimated as the proportion of statistics from randomized data sets that are larger or equal to the proportion observed.

Analysis of genetic variability per loci, per sample, and overall consisted in general descriptive indices estimated with Fstat [26]. Allelic richness (aRich) measured the number of alleles independently of sample size, gene diversity per locus, and the sample used as an unbiased estimator [27], and $F_{is}$ index [28] measured the heterozygote deficit or excess. The $F_{is}$ index was calculated for each sample and the Hardy Weinberg equilibrium was tested within the samples by allele permutations among individuals (1000 times) using Fstat. The 95% confidence interval of the $F_{is}$ values per sample was estimated by bootstrapping the alleles within each sample using the GENETIX v 4.05.2 package [29]. Tests for differences among groups of populations for allelic richness, gene diversity, and $F_{is}$ between localities were carried out using Fstat.

The level of genetic structuring was estimated through different hierarchical analyses of molecular variance (Analysis of Molecular Variance, AMOVA), with Arlequin software v. 3.5.1.2 [30]. The pairwise $F_{ST}$ indices can be used as short-term genetic distances between samples and were computed with Arlequin. Differentiation between samples was tested by different types of permutation for each covariance component. The genetic relationships between samples were inferred with

the neighbor-joining method (NJ) [31] using Cavalli Sforza's genetic distances between pairwise samples. This analysis was conducted in Population v. 1.2.30. software (Olivier Langella, CNRS UPR9034; http://bioinformatics.org/~tryphon/ populations/).

A Bayesian approach was implemented to determine the number of genetic clusters within the entire data set using STRUCTURE v. 2.2 [32,33]. Runs with a 50,000-iteration burn-in period and 50,000 iterations in length, up to 100,000 for Thago Thago populations, with variable genetic clusters ($k = 2–10$) were carried out to identify the best assignment. This procedure places individuals in $k$ clusters chosen in advance and membership coefficients are calculated for each individual in inferred clusters. The “admixture” model for the ancestry of individuals (individuals may have ancestry from various populations) was applied to the data sets with the “Allele Frequencies Correlated” option (meaning that allele frequencies in the different populations are most probably similar), with the other parameters left at their default values. For these runs, the LOCPRIOR option was used, which allows clustering by using sampling locations of individuals as prior information. This option is more appropriate for data sets with few markers and small sample sizes. Five independent runs for each $k$ were produced to assess the consistency of the results across runs, one of the criteria to estimate the best $k$ [32]. Further analysis with STRUCTURE HARVESTED website program (http://taylor0.biology.ucla.edu/structureHarvester/) was applied to visualize likelihood values across the multiple runs performed for the overall data set and the locality data sets and to help determine the best number of clusters corresponding to the data set [34].

GENECLASS2 software [35] was used to detect the putative first-generation migrants within the data set in each locality. To test the probability of being a first-generation migrant, the probability for each bug of being a resident (not a migrant) was expressed using Bayesian assignment criteria [36] and the Monte-Carlo resampling method [37]. Also, the probability that the putative population source of the first-generation migrant was among the subsamples within each locality was tested using the assignment-exclusion test implemented in GENECLASS2 with the same model as for the detection of migrants.

## Results

### Genetic variability overall and within each locality

Data from seven previously described microsatellite loci were analyzed for a total of 277 *T. infestans* specimens collected in three localities in wild and domestic environments (Sapini, Quillacollo and Thago Thago). The overall amplification success of the loci was 97.7%, ranging from 92.8% (locus TiD09) to 99.3% (loci TiE02, TiC09, TiE12). All loci were polymorphic. No significant association between the pairs of locus genotypes within the overall sample and within each sampled locality was found at the nominal 1% level, a result supporting the statistical independence of loci. The number of alleles per locus in the entire sample ranged from three (TiF11) to 19 (TiF03), and the mean allelic richness was 8.0 ± 3.9.
<table><caption>**Table 1.** Genetic indices evaluated among seven microsatellite loci per locality, Sapini, Quillacollo and Thago Thago.</caption><thead><tr><th rowspan="2">Genetic indices</th><th colspan="3">Locality</th></tr><tr><th>Sapini</th><th>Quillacollo</th><th>Thago Thago</th></tr></thead><tbody><tr><td>Av. N ± SD</td><td>37.1 ± 1.6</td><td>115.6 ± 1.7</td><td>116.7 ± 5.0</td></tr><tr><td>Av. aRich ± SD</td><td>2.13 ± 1.33</td><td>5.97 ± 2.58</td><td>7.13 ± 3.3</td></tr><tr><td>Av. GD ± SD</td><td>0.28 ± 0.25</td><td>0.69 ± 0.18</td><td>0.72 ± 0.14</td></tr><tr><td>F<sub>IS</sub></td><td>0.042</td><td>0.064</td><td>0.181</td></tr><tr><td>P-value for F<sub>IS</sub></td><td>0.27</td><td>0.0010*</td><td>0.0005*</td></tr></tbody></table>

Av., average among loci; SD, standard deviation; N, number of individuals; aRich, allelic richness; GD, genetic diversity; P-value for F<sub>IS</sub> was based on 2100 randomisations the nominal level for 5% was 0.00238; * significant value.
doi: 10.1371/journal.pone.0080786.t001

Genetic variability indices in the three samples composed of all analyzed individuals in each locality are presented in Table 1. Significant deviations from HW expectations due to heterozygote deficit (positive F<sub>IS</sub> values significantly different from random predictions) were observed in Quillacollo and Thago Thago but not in Sapini (Table 1). The locality pairwise comparison showed that Sapini differed for all indices (aRich, observed proportion of heterozygotes, gene diversity, and F<sub>IS</sub>) from Quillacollo and Thago Thago (P < 0.02), while only the F<sub>IS</sub> value was significantly different between Quillacollo and Thago Thago (P < 0.01). Samples from each locality were divided into subsamples according to space and ecotopes (intra-peridomestic or wild) and the summary of genetic variability, per locus and per subsample, is presented in Table 2 for Sapini, Quillacollo, and Thago Thago. Remarkably, in Sapini two loci were monomorphic among all the subsamples (TiE02 and TiF11) and up to six loci in the Dom2 and IP11 subsamples, while in Quillacollo and Thago Thago, out of 19 subsamples, only one was monomorphic at one locus, TiF11 (CV119, Quillacollo). Within the 25 subsamples, no significant deviations from HW expectations were observed except for one subsample in Thago Thago (CV78) composed of bugs collected in a single dwelling (significant positive F<sub>IS</sub>).

## Genetic structure and differentiation between intra-peridomestic and wild populations

The genetic structure of the subsamples clustered according to different spatial levels was explored with the analysis of molecular variance (AMOVA) using all loci for distance computation and all data, with 0.12% threshold for missing data (Table 3). Among the 25 subsamples, most of the genetic variation was assigned within individual (78.7% of the total variation, P < 10<sup>-4</sup>, see Table 3), but a significant 19.5% variation was attributed to differences between subsamples, showing some genetic differentiation among them (P < 10<sup>-4</sup>). The first hierarchical AMOVA, grouping subsamples according to their geographic origin (locality), showed a higher variance component among localities (18%) than among subsamples within localities (6.7%), and in both cases the F<sub>CT</sub> (among localities) and F<sub>SC</sub> (among subsamples within localities) indices

were significant (P < 10<sup>-4</sup>), showing geographical structuring between localities and also some possible structuring within localities. This geographical structure was fully confirmed by the topology of the neighbor-joining tree built from Cavalli Sforza's genetic distances between wild and intra-peridomestic populations from the three localities, where the bootstrap values at the nodes grouping each locality were highly significant (Figure 2). To further explore a possible structure between ecotopes (wild/intra-peridomestic), additional hierarchical AMOVA were performed within each locality. In the three localities, the variability between the two kinds of ecotope accounted for a very low total variability rate, 0.30% in Quillacollo, 1.48% in Sapini, and 1.66% in Thago Thago, and the differentiation between the two ecotopes was not significant for Sapini and Quillacollo and borderline for Thago Thago (Table 3). However, in each locality the variance component attributed to differences between subsamples (hence among wild or domestic groups of subsamples) was significant, showing some genetic differentiation between subsamples within ecotopes. Moreover, pairwise comparisons between subsamples using the F<sub>ST</sub> estimator showed significant differentiation between most of the subsamples (94% of the 300 pairwise comparisons, data not shown).

The Bayesian method of the assignment implemented in STRUCTURE detected two clusters in the overall sample (Figure 3a) whether the “admixture” or the “no admixture” model was applied. One genetic cluster was composed of all the Sapini subsamples and the second one all the Quillacollo and Thago Thago samples, indicating a clear geographical differentiation. To explore the population structure in each locality in greater detail, the corresponding subsamples were analyzed separately. In Sapini and Quillacollo, two genetic clusters were detected (Figure 3b and c), and no correlation was observed with the bugs’ capture area, intra-peridomestic vs. wild. The individual patterns of assignment to the two genetic clusters were generally homogeneous within subsamples that were defined based on spatial criteria (grouping collected bugs in the same place or nearby). This pattern may depict the existence of discrete populations, in domestic as well as in wild ecotopes. Similar individual genetic assignments in wild and intra-peridomestic subsamples were also observed in Sapini and Quillacollo: some subsamples were composed of members assigned to a single cluster, while others were composed of individuals that presented an admixed origin. In Thago Thago, there were five significant clusters (Figure 3d, Table 4). Only 50% of the individuals were assigned to one of the five inferred clusters, with a probability greater than 70%. The other individuals came from two or three clusters. Only two populations contained the majority of its individuals assigned to a single cluster (C12 wild and CV78 intra-peridomestic).

## Detection of putative first-generation migrants within localities

In each locality, the putative first-generation migrants were detected on the basis of the highest likelihood ratio calculated between the subsample from which the individual was sampled and among all other subsamples in the locality including its
<table><caption>**Table 2.** Summary of the genetic variability for seven microsatellite loci in intra-peridomestic and wild subsamples of *T. infestans* in Sapini, Quillacollo and Thago localities.</caption><thead><tr><th rowspan="2">Subsample</th><th rowspan="2">Ecotope</th><th colspan="7">Locus</th><th rowspan="2">Overall loci<br>Av. ± SD</th></tr><tr><th>TiA02</th><th>TiE02</th><th>TiC09</th><th>TiF11</th><th>TiF03</th><th>TiD09</th><th>TiE12</th></tr></thead><tbody><tr><td colspan="10"><strong>Sapini</strong></td></tr><tr><td>Dom1</td><td>Intra-peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>6</td><td>6</td><td>6</td><td>4</td><td>6</td><td>6</td><td>6</td><td>5.7 ± 0.7</td></tr><tr><td>aRich</td><td></td><td>1.98</td><td>1</td><td>1.98</td><td>1</td><td>2.91</td><td>2.00</td><td>2.00</td><td>1.84 ± 0.66</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.250</td><td>NA</td><td>0.615</td><td>NA</td><td>0.524</td><td>-0.667</td><td>-0.667</td><td>-0.053</td></tr><tr><td>Dom2</td><td>Peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>6</td><td>6</td><td>6</td><td>6</td><td>5</td><td>6</td><td>6</td><td>5.8 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>2.96</td><td>1</td><td>1</td><td>1.28 ± 0.7</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>-0.333</td><td>NA</td><td>NA</td><td>-0.333</td></tr><tr><td>Syl1</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>7</td><td>7</td><td>7</td><td>6</td><td>7</td><td>7</td><td>7</td><td>6.9 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>1.99</td><td>1</td><td>1</td><td>1</td><td>3.39</td><td>2.00</td><td>2.00</td><td>1.77 ± 0.87</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.500</td><td>NA</td><td>NA</td><td>NA</td><td>-0.263</td><td>-0.714</td><td>-0.714</td><td>-0.525</td></tr><tr><td>IP11</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>7</td><td>7</td><td>7</td><td>7</td><td>7</td><td>7</td><td>7</td><td>7.0 ± 0.0</td></tr><tr><td>aRich</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>2.890</td><td>1</td><td>1</td><td>1.27 ± 0.71</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>-0.412</td><td>NA</td><td>NA</td><td>-0.412</td></tr><tr><td>Z11</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>7</td><td>7</td><td>7</td><td>6</td><td>7</td><td>7</td><td>7</td><td>6.8 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>1.98</td><td>1</td><td>1</td><td>1</td><td>3.24</td><td>1.99</td><td>1.99</td><td>1.74 ± 0.82</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.333</td><td>NA</td><td>NA</td><td>NA</td><td>0.020</td><td>0.143</td><td>0.143</td><td>0.006</td></tr><tr><td>Z10</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>5</td><td>5</td><td>5</td><td>5</td><td>4</td><td>5</td><td>5</td><td>4.8 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>3.00</td><td>1.98</td><td>1.98</td><td>1.56 ± 0.78</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>0.294</td><td>-0.143</td><td>-0.143</td><td>0.077</td></tr><tr><td colspan="10"><strong>Quillacollo</strong></td></tr><tr><td>CV119</td><td>Peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>15</td><td>16</td><td>16</td><td>16</td><td>16</td><td>16</td><td>16</td><td>15.8 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>2.975</td><td>2.659</td><td>2.920</td><td>1.000</td><td>3.980</td><td>3.751</td><td>2.779</td><td>2.86 ± 0.96</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.125</td><td>-0.098</td><td>0.125</td><td>NA</td><td>-0.254</td><td>0.174</td><td>-0.099</td><td>0.002</td></tr><tr><td>CV125</td><td>Peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>9</td><td>9</td><td>9</td><td>8</td><td>8</td><td>5</td><td>7</td><td>7.8 ± 1.5</td></tr><tr><td>aRich</td><td></td><td>2.554</td><td>2.817</td><td>2.817</td><td>1.992</td><td>2.000</td><td>3.000</td><td>1.989</td><td>2.45 ± 0.45</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.474</td><td>-0.655</td><td>-0.067</td><td>-0.273</td><td>0.000</td><td>-0.333</td><td>0.625</td><td>-0.208</td></tr><tr><td>CV134</td><td>Peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>11</td><td>11</td><td>11</td><td>11</td><td>11</td><td>10</td><td>11</td><td>10.8 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>3.688</td><td>3.385</td><td>3.026</td><td>2.387</td><td>5.118</td><td>2.747</td><td>2.883</td><td>3.32 ± 0.90</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.104</td><td>-0.184</td><td>-0.188</td><td>0.059</td><td>0.011</td><td>0.063</td><td>0.101</td><td>-0.046</td></tr><tr><td>CV59</td><td>Peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>18</td><td>18</td><td>18</td><td>18</td><td>18</td><td>18</td><td>18</td><td>18.0 ± 0.0</td></tr><tr><td>aRich</td><td></td><td>4.089</td><td>2.957</td><td>3.287</td><td>1.987</td><td>5.385</td><td>5.100</td><td>3.619</td><td>3.62 ± 1.19</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.215</td><td>-0.172</td><td>0.182</td><td>0.370</td><td>0.056</td><td>0.004</td><td>-0.245</td><td>0.039</td></tr><tr><td>P1</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>15</td><td>15</td><td>15</td><td>15</td><td>15</td><td>15</td><td>15</td><td>15.0 ± 0.0</td></tr><tr><td>aRich</td><td></td><td>3.264</td><td>2.718</td><td>3.766</td><td>1.988</td><td>5.176</td><td>4.261</td><td>3.530</td><td>3.53 ± 1.03</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.131</td><td>0.234</td><td>0.149</td><td>-0.077</td><td>-0.130</td><td>-0.277</td><td>-0.171</td><td>-0.026</td></tr><tr><td>P2</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>14</td><td>14</td><td>14</td><td>14</td><td>14</td><td>14</td><td>14</td><td>14.0 ± 0.0</td></tr><tr><td>aRich</td><td></td><td>3.650</td><td>2.592</td><td>2.899</td><td>1.997</td><td>5.137</td><td>4.181</td><td>4.999</td><td>3.64 ± 1.21</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.004</td><td>-0.020</td><td>-0.156</td><td>0.409</td><td>-0.095</td><td>-0.143</td><td>0.010</td><td>-0.018</td></tr><tr><td>P3</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table 2 (continued).</caption><thead><tr><th rowspan="2">Subsample</th><th rowspan="2">Ecotope</th><th colspan="7">Locus</th><th rowspan="2">Overall loci<br>Av. ± SD</th></tr><tr><th>TiA02</th><th>TiE02</th><th>TiC09</th><th>TiF11</th><th>TiF03</th><th>TiD09</th><th>TiE12</th></tr></thead><tbody><tr><td>N</td><td></td><td>22</td><td>22</td><td>22</td><td>22</td><td>22</td><td>22</td><td>22</td><td>22.0 ± 0.0</td></tr><tr><td>aRich</td><td></td><td>3.577</td><td>3.437</td><td>2.969</td><td>1.744</td><td>3.943</td><td>4.370</td><td>4.958</td><td>3.57 ± 1.03</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.016</td><td>-0.112</td><td>-0.141</td><td>0.784</td><td>0.088</td><td>-0.050</td><td>0.003</td><td>0.007</td></tr><tr><td>P4</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>12</td><td>12</td><td>12</td><td>12</td><td>12</td><td>12</td><td>12</td><td>12.0 ± 0.0</td></tr><tr><td>aRich</td><td></td><td>4.214</td><td>3.190</td><td>2.416</td><td>1.417</td><td>5.706</td><td>5.853</td><td>4.028</td><td>3.83 ± 1.63</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.025</td><td>-0.329</td><td>-0.197</td><td>0.000</td><td>0.143</td><td>-0.052</td><td>-0.111</td><td>-0.063</td></tr><tr><td>Thago Thago</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CV7687</td><td>Intra-peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>13</td><td>12</td><td>12</td><td>13</td><td>13</td><td>11</td><td>13</td><td>12.4 ± 0.8</td></tr><tr><td>aRich</td><td></td><td>3.66</td><td>3.04</td><td>2.41</td><td>2.15</td><td>3.91</td><td>3.73</td><td>3.35</td><td>3.18 ± 0.68</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.032</td><td>0.017</td><td>0.492</td><td>0.231</td><td>0.253</td><td>-0.029</td><td>0.065</td><td>0.130</td></tr><tr><td>CV85</td><td>Peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>4</td><td>4</td><td>4</td><td>4</td><td>3</td><td>4</td><td>4</td><td>3.9 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>2.50</td><td>2.71</td><td>1.75</td><td>2.00</td><td>2.00</td><td>1.96</td><td>4.21</td><td>2.45 ± 0.85</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.500</td><td>-0.286</td><td>0.000</td><td>0.571</td><td>-0.333</td><td>-0.200</td><td>0.455</td><td>0.156</td></tr><tr><td>SN0102</td><td>Intra-peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>17</td><td>19</td><td>18</td><td>18</td><td>18</td><td>16</td><td>19</td><td>17.9 ± 1.1</td></tr><tr><td>aRich</td><td></td><td>3.21</td><td>2.54</td><td>2.36</td><td>1.97</td><td>3.71</td><td>4.02</td><td>3.84</td><td>3.09 ± 0.81</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.074</td><td>0.179</td><td>0.526</td><td>0.452</td><td>0.012</td><td>-0.047</td><td>-0.038</td><td>0.121</td></tr><tr><td>CV82</td><td>Intra-peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>8</td><td>8</td><td>8</td><td>8</td><td>8</td><td>4</td><td>8</td><td>7.4 ± 1.5</td></tr><tr><td>aRich</td><td></td><td>2.94</td><td>2.95</td><td>1.62</td><td>1.97</td><td>3.99</td><td>4.39</td><td>3.00</td><td>2.98 ± 0.99</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.646</td><td>0.041</td><td>-0.077</td><td>0.517</td><td>0.263</td><td>0.182</td><td>-0.029</td><td>0.247</td></tr><tr><td>CV74</td><td>Intra-peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>10</td><td>10</td><td>10</td><td>10</td><td>9</td><td>8</td><td>10</td><td>9.6 ± 0.8</td></tr><tr><td>aRich</td><td></td><td>3.35</td><td>3.02</td><td>3.02</td><td>1.96</td><td>3.90</td><td>4.30</td><td>3.43</td><td>3.28 ± 0.75</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.087</td><td>0.211</td><td>0.211</td><td>0.386</td><td>-0.094</td><td>-0.167</td><td>0.031</td><td>0.066</td></tr><tr><td>CV78</td><td>Intra-peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>13</td><td>13</td><td>13</td><td>12</td><td>13</td><td>13</td><td>13</td><td>12.9 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>3.52</td><td>3.14</td><td>3.47</td><td>1.71</td><td>3.85</td><td>4.21</td><td>3.15</td><td>3.29 ± 0.79</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.114</td><td>0.439</td><td>0.306</td><td>0.436</td><td>0.446</td><td>0.284</td><td>0.208</td><td>0.274*</td></tr><tr><td>CV68</td><td>Peridomestic</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>6</td><td>6</td><td>6</td><td>6</td><td>5</td><td>6</td><td>6</td><td>5.8 ± 0.4</td></tr><tr><td>aRich</td><td></td><td>3.04</td><td>2.82</td><td>1.50</td><td>1.99</td><td>4.50</td><td>3.81</td><td>3.23</td><td>2.98 ± 1.02</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.026</td><td>-0.250</td><td>0.000</td><td>0.706</td><td>0.111</td><td>0.184</td><td>-0.136</td><td>0.081</td></tr><tr><td>C11</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>12</td><td>15</td><td>15</td><td>14</td><td>14</td><td>13</td><td>15</td><td>14.0 ± 1.1</td></tr><tr><td>aRich</td><td></td><td>3.48</td><td>2.18</td><td>2.69</td><td>1.95</td><td>4.23</td><td>3.68</td><td>3.86</td><td>3.15 ± 0.88</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.107</td><td>0.279</td><td>0.590</td><td>0.103</td><td>0.259</td><td>-0.082</td><td>0.099</td><td>0.186</td></tr><tr><td>C12</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>15</td><td>14</td><td>15</td><td>14</td><td>14</td><td>14</td><td>15</td><td>14.4 ± 0.5</td></tr><tr><td>aRich</td><td></td><td>2.77</td><td>2.48</td><td>3.14</td><td>2.11</td><td>3.69</td><td>3.76</td><td>3.78</td><td>3.10 ± 0.67</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.116</td><td>0.027</td><td>0.357</td><td>0.406</td><td>-0.273</td><td>-0.087</td><td>-0.087</td><td>0.010</td></tr><tr><td>C13</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>12</td><td>14</td><td>14</td><td>14</td><td>13</td><td>13</td><td>14</td><td>13.4 ± 0.8</td></tr><tr><td>aRich</td><td></td><td>3.61</td><td>1.64</td><td>2.55</td><td>2.76</td><td>3.47</td><td>4.47</td><td>3.80</td><td>3.18 ± 0.94</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>-0.063</td><td>-0.040</td><td>0.090</td><td>0.266</td><td>-0.109</td><td>0.040</td><td>-0.068</td><td>0.015</td></tr><tr><td>C2</td><td>Wild</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>N</td><td></td><td>5</td><td>5</td><td>5</td><td>5</td><td>5</td><td>5</td><td>5</td><td>5.0 ± 0.0</td></tr><tr><td>aRich</td><td></td><td>3.07</td><td>3.93</td><td>2.00</td><td>1.97</td><td>3.43</td><td>4.17</td><td>4.50</td><td>3.29 ± 1.01</td></tr><tr><td>F<sub>IS</sub></td><td></td><td>0.077</td><td>0.294</td><td>0.667</td><td>0.600</td><td>0.250</td><td>-0.176</td><td>-0.143</td><td>0.180</td></tr></tbody></table>
<table><caption>**Table 2 (continued).**</caption><thead><tr><th rowspan="2">Subsample</th><th rowspan="2">Ecotope</th><th colspan="8">Locus</th></tr><tr><th>TiA02</th><th>TiE02</th><th>TiC09</th><th>TiF11</th><th>TiF03</th><th>TiD09</th><th>TiE12</th><th>Overall loci</th></tr></thead><tbody><tr><td colspan="10">Av. ± SD</td></tr></tbody></table>

N, number of individuals; aRich, allelic richness; Av., average among loci; SD, standard deviation; *P*-value for *F*<sub>IS</sub> was based on 4200 randomizations, the nominal level for 5% was 0.00119, there was no significant value; NA, not available.

doi: 10.1371/journal.pone.0080786.t002

<table><caption>**Table 3.** Analysis of Molecular Variance (AMOVA) for *T. infestans* samples according to different hierarchical levels.</caption><thead><tr><th>Hierarchical level</th><th>Source of variation</th><th>d. f.</th><th>% of variation</th><th>F-index</th><th>P-value*</th></tr></thead><tbody><tr><td rowspan="3">Subsamples (25)</td><td>Among subsamples</td><td>24</td><td>19.50</td><td>F<sub>ST</sub> = 0.195</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td>Among individuals within subsamples</td><td>252</td><td>1.82</td><td>F<sub>IS</sub> = 0.022</td><td>0.15</td></tr><tr><td>Within individuals</td><td>277</td><td>78.70</td><td>F<sub>IT</sub> = 0.213</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td rowspan="4">Localities (Sapini, Quillacollo, Thago Thago)</td><td>Among localities</td><td>2</td><td>17.98</td><td>F<sub>CT</sub> = 0.180</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td>Among subsamples within localities</td><td>22</td><td>6.74</td><td>F<sub>SC</sub> = 0.082</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td>Among individuals within subsamples</td><td>252</td><td>1.70</td><td>F<sub>IS</sub> = 0.023</td><td>0.15</td></tr><tr><td>Within individuals</td><td>277</td><td>73.59</td><td>F<sub>IT</sub> = 0.264</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td rowspan="4">Ecotopes in Sapini (wild / peri-intra domestic)</td><td>Among ecotopes</td><td>1</td><td>1.48</td><td>F<sub>CT</sub> = 0.015</td><td>0.40</td></tr><tr><td>Among subsamples within ecotopes</td><td>4</td><td>21.46</td><td>F<sub>SC</sub> = 0.218</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td>Among individuals within subsamples</td><td>32</td><td>-16.81</td><td>F<sub>IS</sub> = -0.218</td><td>0.98</td></tr><tr><td>Within individuals</td><td>38</td><td>93.86</td><td>F<sub>IT</sub> = 0.061</td><td>0.42</td></tr><tr><td rowspan="4">Ecotopes in Quillacollo (wild / peridomestic)</td><td>Among ecotopes</td><td>1</td><td>0.30</td><td>F<sub>CT</sub> = 0.003</td><td>0.31</td></tr><tr><td>Among subsamples within ecotopes</td><td>6</td><td>8.59</td><td>F<sub>SC</sub> = 0.086</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td>Among individuals within subsamples</td><td>109</td><td>-2.57</td><td>F<sub>IS</sub> = -0.028</td><td>0.87</td></tr><tr><td>Within individuals</td><td>117</td><td>93.69</td><td>F<sub>IT</sub> = 0.063</td><td>0.01</td></tr><tr><td rowspan="4">Ecotopes in Thago Thago (wild / peri-intra domestic)</td><td>Among ecotopes</td><td>1</td><td>1.66</td><td>F<sub>CT</sub> = 0.017</td><td>0.048</td></tr><tr><td>Among subsamples within ecotopes</td><td>9</td><td>5.61</td><td>F<sub>SC</sub> = 0.057</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td>Among individuals within subsamples</td><td>111</td><td>11.82</td><td>F<sub>IS</sub> = 0.127</td><td>&lt; 10<sup>-4</sup></td></tr><tr><td>Within individuals</td><td>122</td><td>80.91</td><td>F<sub>IT</sub> = 0.190</td><td>&lt; 10<sup>-4</sup></td></tr></tbody></table>

d.f., degree of freedom; * *P*-value for *F*-index based on 10,000 permutations.

doi: 10.1371/journal.pone.0080786.t003

own subsample (Table 5). Putative first-generation migrants were detected in 12 of 25 samples, the migrant rates were 5.3% (two migrants) in Sapini, 2.6% (three migrants) in Quillacollo, and 7.4% (nine migrants) in Thago Thago. Of the 14 putative migrants, four were males, four females, and six third-, fourth-, and fifth-stage nymphs. Eight were captured in intra-peridomicile and seven in wild ecotopes, six could have originated from one or several other subsamples of the corresponding locality (the assignment-exclusion test in GENECLASS software provided a significant probability P<sub>2</sub> > 0.70, Table 5). The results supported the displacement of bugs from wild to intra-peridomestic ecotopes in Sapini and Thago Thago, and also displacements between intra-peridomestic ecotopes in Thago Thago. In Quillacollo, the origin of the migrant remained unknown.

## Discussion

### Genetic diversity in *T. infestans*

The overall genetic diversity of *T. infestans* observed in Quillacollo and Thago Thago is high and similar to that previously observed in the same area analyzing sylvatic populations with the same set of microsatellite loci [38]. Strong genetic variability has also been found in several Argentinean populations, each composed of insects from several houses in the same village [39] and even at the micro-geographical scale among house compound populations [16]. Therefore, Sapini seems to be an exception with significantly much lower genetic diversity than in Quillacollo and Thago Thago (0.28 vs. around 0.70). In Sapini two loci were monomorphic, TiF11 and TiE02, analyzed for 34 and 38 insects, respectively. Low genetic variability was also previously detected in Sapini upon analysis of the polymorphism of cytochrome b sequences; indeed, a single haplotype belonging to a “special group of haplotypes” composed of very divergent haplotypes was identified [40]. This group of divergent haplotypes was sampled only in
<figure><img src="image_5.png" /><figcaption>**Figure 2. Neighbor-joining between wild and intra-peridomestic populations.** The tree was built from Cavalli Sforza's genetic distances between the six populations that grouped, in each locality (bugs collected in the wild and the intra-peridomestic environments). The bootstrap values are expressed at the nodes of the branches.</figcaption></figure>

doi: 10.1371/journal.pone.0080786.g002

northern Andean valleys of La Paz Department, but not in the Cochabamba region. The authors suggested the intervention of isolation and bottleneck phenomena due to the particular feature of these valleys, highly fragmented, with extreme weather conditions and where many endemic species have evolved [21]. The low genetic diversity currently observed for microsatellite polymorphisms in these populations is in agreement with this hypothesis.

## Panmictic unit

The search for the size of the natural panmictic unit contributes to determining the geographical scale where random mating among members of a species is potentially possible. The putative panmictic unit is the most suitable to define populations for population genetics analysis to derive the relationships that exist between the members of these populations at different geographical scales, which is needed for control strategies. The allozyme polymorphism analyses first showed that the *T. infestans* members within a village evolve in agreement with the assumption of panmixia [41]; then it was proposed that the panmictic unit can be smaller, such as populations limited to individual houses or peridomestic structures in a village [20]. Other studies with microsatellite markers generally agree with this last assumption for *T*.

*infestans* [16,39,41] as well as for other species of triatomine such as *Triatoma longipennis* for which the panmictic unit corresponds to a collection of insects captured in a single peridomestic structure [42]. However, for *Triatoma dimidiata* from Guatemalan populations, the results suggest that the size of the panmictic unit is not the same in different regions, probably due to variable migration rates between populations [43]. Here, Hardy-Weinberg equilibrium was observed within all populations studied except one. The above-mentioned overall results show that, to develop genetic studies of triatomine populations, the best way is to consider a group of insects collected in a single capture site or sites that are extremely close (a few meters) in the same period, as a population. This approach better assesses the genetic relationships between populations in each environment and can infer the dispersal ability of insects in a given area, a frequently raised question that is directly related to the epidemiology of the disease.

## Genetic structure in *T. infestans*

In the present case, a strong geographic structure was observed between the three localities (the $F_{ST}$ value is 0.179, $P < 10^{-5}$), up to 242 km away (Sapini to Thago Thago) as observed between populations collected in different localities in Argentina over long distances [44]. However, high genetic
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Plots of ancestry estimates of *T. infestans* individuals using STRUCTURE. software. a - Plot of the individuals belonging to the 25 study populations, one to four intra-peridomestic populations from Quillacollo, five to eight wild populations from Quillacollo, nine and ten intra-peridomestic populations from Sapini, 11 to 14 wild populations from Sapini, 15 to 22 intra-peridomestic populations from Thago Thago, 23 to 25 wild populations from Thago Thago. b - Plot of the individuals belonging to the six populations in Sapini (one and two intra-peridomestic, three to six wild). c - Plot of the individuals belonging to the eight populations in Quillacollo (one to four intra-peridomestic, five to eight wild). d - Plot of the individuals belonging to the 11 populations in Thago Thago (one to eight intra-peridomestic, nine to 11 wild). Each individual in the data set is represented in the plots by a single vertical line, which is partitioned into different colored segments representing the individual membership estimates in each of the inferred clusters (k).</figcaption></figure>

doi: 10.1371/journal.pone.0080786.g003
<table><caption>**Table 4.** Assignement of individual genotypes of *T. infestans* belonging to 11 sub samples of Thago Thago to five inferred clusters.</caption><thead><tr><th rowspan="3">Subsample</th><th rowspan="3">No. of T. infestans</th><th colspan="5">Number of individuals assigned to each inferred cluster with a probability > 0.70</th><th colspan="2">No assigned</th></tr><tr><th colspan="5">Five inferred clusters (k = 5)</th><th rowspan="2">4</th><th rowspan="2">5</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr></thead><tbody><tr><td>CV7687</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>5</td><td>8</td></tr><tr><td>CV85</td><td>4</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td></tr><tr><td>SN0102</td><td>19</td><td>0</td><td>0</td><td>12</td><td>0</td><td>0</td><td>0</td><td>7</td></tr><tr><td>CV82</td><td>8</td><td>0</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>5</td></tr><tr><td>CV74</td><td>10</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>10</td></tr><tr><td>CV78</td><td>13</td><td>0</td><td>9</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td>CV68</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>3</td></tr><tr><td>C11</td><td>15</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>4</td><td>11</td></tr><tr><td>C12</td><td>15</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td></tr><tr><td>C13</td><td>14</td><td>0</td><td>0</td><td>2</td><td>0</td><td>7</td><td>0</td><td>5</td></tr><tr><td>C2</td><td>5</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3</td></tr></tbody></table>

doi: 10.1371/journal.pone.0080786.t004

<table><caption>**Table 5.** Putative first generation migrants and possible origin.</caption><thead><tr><th rowspan="2">T. infestans code</th><th rowspan="2">Locality</th><th rowspan="2">Stage</th><th colspan="3">Population where the individual was sampled</th><th colspan="3">Putative original population</th><th rowspan="2">Distance between populations (m)</th><th rowspan="2">Putative ecotope displacement**</th></tr><tr><th>Ecotope</th><th>Name</th><th>P<sub>1</sub></th><th>Ecotope</th><th>Name</th><th>P<sub>2</sub></th></tr></thead><tbody><tr><td>Sap 398</td><td>Sapini</td><td>Female</td><td>Intra-peri</td><td>Dom1</td><td>&lt; 0.05</td><td>Wild</td><td>Syl1</td><td>0.95</td><td>65</td><td>Wild to domestic</td></tr><tr><td>Sap 471</td><td>Sapini</td><td>5<sup>th</sup> nymph</td><td>Wild</td><td>Z10</td><td>&lt; 0.05</td><td>Wild</td><td>Z11</td><td>0.76</td><td>320</td><td>Wild to wild</td></tr><tr><td>Qui 253</td><td>Quillacollo</td><td>4<sup>th</sup> nymph</td><td>Peri</td><td>CV125</td><td>&lt; 0.01</td><td>Wild</td><td>P1</td><td>0.21</td><td>250</td><td>uk</td></tr><tr><td>Qui 379</td><td>Quillacollo</td><td>4<sup>th</sup> nymph</td><td>Peri</td><td>CV59</td><td>&lt; 0.01</td><td>Wild</td><td>P1</td><td>0.27</td><td>175</td><td>uk</td></tr><tr><td>Qui 659</td><td>Quillacollo</td><td>Male</td><td>Wild</td><td>P2</td><td>&lt; 0.001</td><td>Wild</td><td>P3</td><td>0.45</td><td>233</td><td>uk</td></tr><tr><td>Tor 622</td><td>Thago Thago</td><td>5<sup>th</sup> nymph</td><td>Intra-peri</td><td>CV7687</td><td>&lt; 0.01</td><td>Wild</td><td>C13</td><td>&lt; 0.01</td><td>150</td><td>uk</td></tr><tr><td>Tor 369</td><td>Thago Thago</td><td>Female</td><td>Intra-peri</td><td>SN0102</td><td>&lt; 0.001</td><td>Intra-peri</td><td>CV78</td><td>&lt; 0.006</td><td>120</td><td>uk</td></tr><tr><td>Tor 404</td><td>Thago Thago</td><td>Female</td><td>Intra-peri</td><td>CV82</td><td>&lt; 0.01</td><td>Wild</td><td>C12</td><td>0.72</td><td>125</td><td>Wild to domestic</td></tr><tr><td>Tor 149</td><td>Thago Thago</td><td>Male</td><td>Intra-peri</td><td>CV78</td><td>&lt; 0.001</td><td>Intra-peri</td><td>CV82 or SN0102*</td><td>0.86</td><td>130</td><td>Domestic to domestic</td></tr><tr><td>Tor 246</td><td>Thago Thago</td><td>Male</td><td>Peri</td><td>CV68</td><td>&lt; 0.001</td><td>Wild</td><td>C11</td><td>0.42</td><td>1000</td><td>uk</td></tr><tr><td>Tor 494</td><td>Thago Thago</td><td>Male</td><td>wild</td><td>C11</td><td>&lt; 0.01</td><td>Intra-peri</td><td>CV7687</td><td>0.39</td><td>40</td><td>uk</td></tr><tr><td>Tor 435</td><td>Thago Thago</td><td>1<sup>st</sup> nymph</td><td>Intra-peri</td><td>CV7687</td><td>&lt; 0.01</td><td>Wild</td><td>C13</td><td>0.49</td><td>150</td><td>uk</td></tr><tr><td>Tor 442</td><td>Thago Thago</td><td>Female</td><td>Wild</td><td>C13</td><td>&lt; 0.001</td><td>Wild</td><td>C12</td><td>0.84</td><td>140</td><td>Wild to wild</td></tr><tr><td>Tor 502</td><td>Thago Thago</td><td>4<sup>th</sup> nymph</td><td>Wild</td><td>C13</td><td>&lt; 0.01</td><td>Wild</td><td>C2</td><td>0.80</td><td>200</td><td>Wild to wild</td></tr></tbody></table>

P<sub>1</sub> is the probability to be a resident that is to say to be no a first generation migrant; P<sub>2</sub> is the probability to originate from the putative population detected by the test; * both subsamples had the same probability; ** domestic includes peri or intra domestic origins, putative origin was considered when P<sub>2</sub> was > 0.70; uk, unknown; GENECLASS software was applied for these analysis.

doi: 10.1371/journal.pone.0080786.t005

structure was also depicted at the locality scale among wild and intra-peridomestic populations. It is increasingly common and expected to see strong structuring among populations from peridomestic collection sites because the bugs may be less inclined to disperse in this medium than in sylvatic areas, since food host mammals are locally more stable. Few data exist concerning wild sites. The previous data regarding wild populations collected near Quillacollo, whose maximal distance is about 1 km, showed differentiation between only some of

them [38], while in the current analyses of the four wild populations of Quillacollo, most of the pairwise F<sub>ST</sub> value comparisons showed a significant differentiation (of 18 comparisons, 14 cases were significant). These results also support restricted gene flow among wild populations at a small geographic scale that can be explained by the conjunction of two phenomena: an intrinsically low ability of *T. infestans* to move actively by walking or flying in the case of adults and the
foundation of colonies mostly created by few founder
specimens (founder effect).

Even in this context, the dispersion remains an important
parameter that leads to major epidemiological consequences,
such as the colonization of new ecotopes. The factors
influencing dispersion are still poorly known. In addition to the
physical parameters, the search for blood meals is mandatory
for triatomine survival. This feeding habit links the insects to
their guests, whose dynamics, whatever the medium, is
strongly dependent on the level of anthropization. The
environmental conditions and their changes may influence the
dispersal patterns, which in fact could be very different
depending on the environmental pressures.

Genetic flow of *T. infestans* between wilderness and human habitat

The analysis of genetic flow between *T. infestans* collected in
wild or domestic areas is a very new issue in understanding
reinfestation events in areas treated by insecticide spraying.
Indeed, *T. infestans* has long been considered exclusively
domestic except in the valleys of Cochabamba, Bolivia;
therefore, the problem was geographically limited and not fully
treated. This is not the case of the vector species, which plays
an epidemiological role in other endemic areas and whose wild
populations are recognized as threatening because of their
dispersion to inhabited areas [45-47].

In the case of *T. infestans*, the first work of morphometry applied to domestic and wild populations from the Cochabamba valley highlighted differences that have been interpreted as incipient separation [48] because they had not been observed by enzymatic analysis before [49]. Moreover, in this same area, the comparison of reinfesting populations with domestic populations before treatment and the surrounding wild populations, using the morphometric marker, have led to the conclusion that reinfesting specimens came from residual populations [50]. In contrast, other analyses, based on mitochondrial sequence comparison, showed shared haplotypes between wild and domestic bugs, indicating a common origin [5,17,40,51]. In addition, major variability was observed in the wild environment, one argument in favor of the origin of dispersion from the wild to human habitats [40]. At a small geographic scale, microsatellite markers appeared to be more informative of bug movement between populations of different ecotopes, but few data are currently available for *T. infestans*. In Argentina, a microsatellite polymorphism study argues for genetic differentiation between a wild and a peridomestic population, 7 km apart [51], while a significant link between two sylvatic foci and the nearest domestic/peridomestic bug populations (100 m and 1.2 km away) was found in an Argentinean area with recurrent reinfestation, raising the hypothesis of sylvatic contributors to reinfestation [5].

In the current study, the analytical effort to estimate gene
flow between different ecotopes is based on (1) the sampling
strategy that grouped the bugs in potentially panmictic
populations using the spatial criterion and (2) the use of
microsatellites, which are highly polymorphic markers.

The absence of significant differentiation between the insects
collected in intra-peridomestic and nearby wild environments
favors the hypothesis of bug dispersion between ecotopes and
the putative sylvatic origin of the reinfesting bugs. The analyses
conducted with the STRUCTURE program confirmed this
result, showing that close intra-peridomestic and wild
populations may be composed of members that are genetically
very similar.

Epidemiological implications of dispersing from
wilderness to human habitat

The first study by [5] and the present one strongly suggest
that wild populations of *T. infestans* can reach intra-
peridomestic areas. For *T. infestans*, it seems that this flow of
insects between the two environments has previously existed
in different areas because a higher level of genetic relatedness
exists between wild and intra-peridomestic in each
geographical area than between areas [40]. This process could
originate from multiple geographical domestication events,
more likely than the occurrence of a single event, in the
Bolivian Andes, followed by a recent expansion using passive
transport associated with human migrations [13,52,53]. The
consequences are the recurrent arrival of new bugs after
villages have been sprayed with insecticide and recolonization
primarily of peridomestic areas [54]. To date, there is no data
showing specific properties of wild populations that can be
related to a lower adaptive capacity to the intra-peridomestic
environment compared to the domestic environment.
Remarkably, we have shown that these wild populations,
discovered in abundance in the Andean valleys of Bolivia, also
feed on humans when the latter conduct outdoor activities [55].
More generally, the feeding behavior of numerous triatomine
species is eclectic, based on a wide range of mammal species
blood meal sources. Interestingly, they can also feed on birds
and even reptiles. Several studies have shown that triatomine
species sporadically found indoors (incursive species) easily
feed on humans when they enter dwellings [56,57]. It now
seems clear that the wild populations of *T. infestans* represent
a health danger, but few data are available on the kinetics of
the recolonization process of the human habitat by these
populations.

The Southern Cone Initiative to Control/Eliminate Chagas
(INCOSUR, http://www1.paho.org/english/hcp/hct/dch/
incosur.htm. Accessed 2013 Oct 21), launched in 1991 to
control *T. infestans*, has met extensive success in Brazil, Chile,
Paraguay, and Uruguay, but control failures have been
declared in the Gran Chaco of Argentina, Bolivia, and
Paraguay [17,58,59]. In accordance with the objectives of the
program, transmission interruption has been officially declared,
region by region, wherever mean dwelling infestation rates
below 3% are reported. In 2011, transmission interruption was
achieved in La Paz Department. This certification was based
on the indisputable effectiveness of the work carried out for
several years in this department where the rate of dwelling
infestation in the endemic municipalities concerned has
decreased from 40% to less than 1% (http://www.paho.org/
Southern Cone Initiative to Control/Eliminate Chagas Disease
(INCOSUR), Annual meeting 2011). The certification process is
a tool that validates the efforts made and the effectiveness of the eradication actions, but monitoring by health authorities often breaks down thereafter, and these attitudes should be precisely examined to prevent transmission risk from persisting especially in areas (Bolivian Andes) where population flows between wild and domestic environments are strongly suspected. This issue was recently fully addressed by emphasizing the role of native vector populations, which occur in natural ecotopes and reinvade/reinfest human habitats [60]. Native populations can be composed of members of the same species as those combatted in the intra-peridomestic environment or by other species. In fact, endemic regions where triatomines do not exist in the wilderness are probably uncommon. Also, public monitoring of Chagas disease should be maintained. This inventory, most particularly the discovery of a far wider distribution of wild *T. infestans* populations than initially described, requires further research to understand the dynamics of wild population movements and their ability to adapt to new ecotopes. In monitoring, emphasis on inhabitant-based participation rather than the systematic search for triatomines by professional teams, which is very costly, could be considered by the health authorities, as well as prioritizing awareness and education campaigns for the human populations at risk [57,61]. Sustained action over time is therefore the major challenge to fight against this new paradigm of Chagas disease.

It would be incomplete not to mention that the importance of the persistence of dwelling infestation by *T. infestans* in the

Gran Chaco, is more likely related to other phenomena than reinfestation by wild populations; one of the involved mechanisms is the development of a pyrethroid resistance [62-64]; however, other factors such as the reduced effectiveness of pyrethroids in peridomestic environment (lack of insecticide stability caused by climatic conditions, intricate peridomestic structures) or the incomplete covering of the village leaving uncontrolled/untreated dwellings, may play a role [11]. In fact, these other phenomena also need to be explored in any area where the reinfestation persists.

## Acknowledgements

We are particularly grateful to the leaders of the Inlasa (Instituto de Laboratorios de Salud, La 390 Paz Bolivia), Dr. Walter Agreda, and Dr. Tamara Chavez, director of the Department of Entomology, for having hosted this work. We dedicate this work to the memory of our colleague François Noireau who worked on the epidemiological impact of wild *T. infestans* populations.

## Author Contributions

Conceived and designed the experiments: SFB MFB EW CB. Performed the experiments: SFB RS RB PB VS MFB EW SD CB. Analyzed the data: SFB CB. Contributed reagents/ materials/analysis tools: RS PB VS MFB. Wrote the manuscript: SFB CB.

## References

1. Noireau F, Bastrenta B, Catalá S, Dujardin JP, Panzera F et al. (2000) Sylvatic population of *Triatoma infestans* from the Bolivian Chaco: from field collection to characterization. Mem Inst Oswaldo Cruz 95 Suppl 1: 119-122. doi:10.1590/S0074-02762000000700020. PubMed: 11142700.

2. 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(5): 633-641. doi:10.1590/S0074-02762010000500006. PubMed: 20835609.

3. Ceballos LA, Piccinali RV, Berkunsky I, Kitron U, Gürtler RE (2009) First finding of melanic sylvatic *Triatoma infestans* (Hemiptera: Reduviidae) colonies in the Argentine Chaco. J Med Entomol 46(5): 1195-1202. doi:10.1603/033.046.0530. PubMed: 19769054.

4. 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(4): 574-579. doi:10.4269/ajtmh.2010.09-0325. PubMed: 20348501.

5. Ceballos LA, Piccinali RV, Marcet PL, Vazquez-Prokopec GM, Cardinal MV et al. (2011) Hidden sylvatic foci of the main vector of Chagas disease *Triatoma infestans*: threats to the vector elimination campaign? PLoS Negl Trop. Drosophila Inf Service 5(10): e1365.

6. Rolón M, Vega MC, Román F, Gómez A, Rojas de Arias A et al. (2011) First Report of Colonies of Sylvatic *Triatoma infestans* (Hemiptera: Reduviidae) in the Paraguayan Chaco, Using a Trained Dog. PLoS Negl Trop. Drosophila Inf Service 5(5): e1026.

7. Waleckx E, Depickère S, Salas R, Aliaga C, Monje M et al. (2012) New discoveries of sylvatic *Triatoma infestans* (Hemiptera: Reduviidae) throughout the Bolivian Chaco. Am J Trop Med Hyg 86(3): 455-458. doi:10.4269/ajtmh.2012.11-0205. PubMed: 22403316.

8. Gorla D, Noireau F (2010) Geographic distribution of Triatominae vectors in America. In: J TelleriaM Tibayrenc. American Trypanosomiasis Chagas disease. Burlington, USA: Elsevier. pp. 209-231.

9. Schofield CJ, Jannin J, Salvatella R (2006) The future of Chagas disease control. Trends Parasitol 22(12): 583-588. doi:10.1016/j.pt. 2006.09.011. PubMed: 17049308.

10. Dias JC (2007) Southern Cone Initiative for the elimination of domestic populations of *Triatoma infestans* and the interruption of transfusional Chagas disease. Historical aspects, present situation, and perspectives. Mem Inst Oswaldo Cruz 102 Suppl 1: 11-18. doi:10.1590/ S0074-02762007005000092. PubMed: 17891281.

11. Gürtler RE (1999) Control campaigns against *Triatoma infestans* in a rural community of northwestern Argentina. Medicina 59 (Suppl 2): 47-54.

12. Noireau F, Cortez MG, Monteiro FA, Jansen AM, Torrico F (2005) Can wild *Triatoma infestans* foci in Bolivia jeopardize Chagas disease control efforts? Trends Parasitol 21(1): 7-10. doi:10.1016/j.pt. 2004.10.007. PubMed: 15639733.

13. Cortez MR, Monteiro FA, Noireau F (2010) New insights on the spread of *Triatoma infestans* from Bolivia-implications for Chagas disease emergence in the southern cone. Infect Genet Evol 10(2): 350-353. doi: 10.1016/j.meegid.2009.12.006. PubMed: 20060504.

14. Cecere MC, Gürtler RE, Canale D, Chuit R, Cohen JE (1997) The role of the peridomestic area in the elimination of *Triatoma infestans* from rural Argentine communities. Rev Panam Salud Publica 1(4): 273-279. doi:10.1590/S1020-49891997000400003. PubMed: 9149523.

15. Cecere MC, Vázquez-Prokopec GM, Gürtler RE, Kitron U (2006) Reinfestation sources for Chagas disease vector, *Triatoma infestans*, Argentina. Emerg Infect Dis 12(7): 1096-1102. doi:10.3201/ eid1207.051445. PubMed: 16836826.

16. Marcet PL, Mora MS, Cutrera AP, Jones L, Gürtler RE et al. (2008) Genetic structure of *Triatoma infestans* populations in rural communities of Santiago del Estero, northern Argentina. Infect Genet Evol 8(6): 835-846. doi:10.1016/j.meegid.2008.08.002. PubMed: 18773972.

17. Quisberth S, Waleckx E, Monje M, Chang B, Noireau F et al. (2011) "Andean" and "non-Andean" ITS-2 and mtCytB haplotypes of *Triatoma infestans* are observed in the Gran Chaco (Bolivia): population genetics and the origin of reinfestation. Infect Genet Evol 11(5): 1006-1014. doi: 10.1016/j.meegid.2011.03.014. PubMed: 21457795.
18. Cecere MC, Gürtler RE, Canale DM, Chuit R, Cohen JE (2002) Effects of partial housing improvement and insecticide spraying on the reinfestation dynamics of *Triatoma infestans* in rural northwestern Argentina. *Acta Trop* 84(2): 101-116. doi:10.1016/S0001-706X(02)00183-3. PubMed: 12429427.

19. Breniere SF, Bosseno MF, Vargas F, Yaksic N, Noireau F et al. (1998) Smallness of the panmictic unit of *Triatoma infestans* (Hemiptera: Reduviidae). J Med Entomol 35(6): 911-917. PubMed: 9835679.

20. Pizarro JC, Gilligan LM, Stevens L (2008) Microsatellites reveal a high population structure in *Triatoma infestans* from Chuquisaca, Bolivia. PLoS Neglec Trop. Drosophila Inf Service 2(3): e202.

21. Ibisch PL, Beck SG, Gerckmann B, Carretero A (2008) La Diversidad biológica. In: PL IbischG Mérida. Biodiversidad: la riqueza de Bolivia Estado de conocimiento y conservación. Santa Cruz de la Sierra: FAN Bolivia

22. Noireau F, Abad-Franch F, Valente SA, Dias-Lima A, Lopes CM et al. (2002) Trapping Triatominae in silvatic habitats. Mem Inst Oswaldo Cruz 97(1): 61-63. doi:10.1590/S0074-02762002000100009. PubMed: 11992149.

23. Edwards JK (1998) Miniprep procedures for the isolation of plant DNA. London: Chapman & Hall.

24. García BA, Zheng L, Pérez de Rosas AR, Segura EL (2004) Isolation and characterization of polymorphic microsatellite loci in the Chagas'disease vector *Triatoma infestans* (Hemiptera: Reduviidae). Mol Ecol Notes 4: 568-571. doi:10.1111/j.1471-8286.2004.00735.x.

25. Goudet J (2001) FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Available: http://www2.unil.ch/popgen/softwares/fstat.htm; Updated from Goudet (1995). Accessed 2013 Oct 21

26. Goudet J (1995) FSTAT (Version 1.2): A Computer Program to Calculate F-Statistics. J Hered 86: 485-486.

27. Nei M (1987) Molecular Evolutionary Genetics. Columbia University Press, New York.

28. Weir BS, Cockerham CC (1984) Estimating F-Statistics for the Analysis of Population-Structure. Evolution 38(6): 1358-1370. doi: 10.2307/2408641.

29. Belkhir K, Borsa P, Goudet J, Chikhi L, Bonhomme F et al. (2001) GENETIX 4.02, logiciel sous WindowsTM pour la génétique des populations. Montpellier Laboratoire Génome, Populations, Interactions CNRS UMR 5000, Université de Montpellier II (France).

30. Excoffier LG, Laval G, Schneider S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics 1: 47-50.

31. Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4(4): 406-425. PubMed: 3447015.

32. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155(2): 945-959. PubMed: 10835412.

33. Falush D, Stephens M, Pritchard JK (2007) Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol Ecol Notes 7(4): 574-578. doi:10.1111/j.1471-8286.2007.01758.x. PubMed: 18784791.

34. Earl DA, vonHoldt BM (2012) Structure Havester: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour 4(2): 359-361. doi:10.1007/s12686-011-9548-7.

35. 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(6): 536-539. doi:10.1093/jhered/esh074. PubMed: 15475402.

36. Rannala B, Mountain JL (1997) Detecting immigration by using multilocus genotypes. Proc Natl Acad Sci U S A 94(17): 9197-9201. doi:10.1073/pnas.94.17.9197. PubMed: 9256459.

37. 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(1): 55-65. doi:10.1046/j.1365-294X.2004.02008.x. PubMed: 14653788.

38. 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(6): 759-764. doi:10.1111/j.1365-3156.2007.01846.x. PubMed: 17550473.

39. 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(7): 1401-1412. doi:10.1111/j.1365-294X.2007.03251.x. PubMed: 17391265.

40. Waleckx E, Salas R, Huamán N, Buitrago R, Bosseno MF et al. (2011) New insights on the Chagas disease main vector *Triatoma infestans* (Reduviidae, Triatominae) brought by the genetic analysis of Bolivian sylvatic populations. Infect Genet Evol 11(5): 1045-1057. doi:10.1016/j.meegid.2011.03.020. PubMed: 21463708.

41. Dujardin JP, Schofield CJ, Tibayrenc M (1998) Population structure of Andean *Triatoma infestans*: allozyme frequencies and their epidemiological relevance. Med Vet Entomol 12(1): 20-29. doi: 10.1046/j.1365-2915.1998.00076.x. PubMed: 9513935.

42. Brenière SF, Waleckx E, Magallón-Gastélum E, Bosseno MF, Hardy X et al. (2012) Population genetic structure of *Meccus longipennis* (Hemiptera, Reduviidae, Triatominae), vector of Chagas disease in West Mexico. Infect Genet Evol 12(2): 254-262. doi:10.1016/j.meegid.2011.11.003. PubMed: 22142488.

43. Dorn PL, Monroy C, Curtis A (2007) *Triatoma dimidiata* (Latreille, 1811): a review of its diversity across its geographic range and the relationship among populations. Infect Genet Evol 7(2): 343-352. doi: 10.1016/j.meegid.2006.10.001. PubMed: 17097928.

44. Pérez de Rosas AR, Segura EL, Fichera L, García BA (2008) Macrogeographic and microgeographic genetic structure of the Chagas' disease vector *Triatoma infestans* (Hemiptera: Reduviidae) from Catamarca, Argentina. Genetica 133(3): 247-260. doi:10.1007/s10709-007-9208-8. PubMed: 17885811.

45. Guhl F, Pinto N, Aguilera G (2009) Sylvatic triatominae: a new challenge in vector control transmission. Mem Inst Oswaldo Cruz 104(4): 71-75. PubMed: 19753461.

46. Barbu C, Dumonteil E, Gourbière S (2011) Evaluation of spatially targeted strategies to control non-domiciliated *Triatoma dimidiata* vector of Chagas disease. PLoS Negl Trop. Drosophila Inf Service 5(5): e1045.

47. Gourbière S, Dorn P, Tripet F, Dumonteil E (2012) Genetics and evolution of triatomines: from phylogeny to vector control. Heredity (Edinb) 108(3): 190-202. doi:10.1038/hdy.2011.71. PubMed: 21897436.

48. Dujardin JP, Bermudez H, Schofield CJ (1997a) The use of morphometrics in entomological surveillance of sylvatic foci of *Triatoma infestans* in Bolivia. Acta Trop 66(3): 145-153. doi:10.1016/S0001-706X(97)00038-7. PubMed: 9210965.

49. 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(1): 40-45. PubMed: 3546693.

50. Dujardin JP, Bermudez H, Casini C, Schofield CJ, Tibayrenc M (1997b) Metric differences between silvatic and domestic *Triatoma infestans* (Heteroptera: Reduviidae) in Bolivia. J Med Entomol 34(5): 544-551. PubMed: 9379460.

51. Piccinini 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(5): 895-903. doi:10.1016/j.meegid.2011.02.013. PubMed: 21352954.

52. Usinger RL, Wygodzinsky P, Ryckman RE (1966) The biosystematics of Triatominae. Annu Rev Entomol 11: 309-330. doi:10.1146/annurev.en.11.010166.001521. PubMed: 5321579.

53. Schofield CJ (1988) Biosystematics of the Triatominae. Oxford: In: M.W. Service, Biosystematics of Haematophagous Insects, Clarendon Press. 55: 284-312.

54. Gurevitz JM, Ceballos LA, Gaspe MS, Alvarado-Otegui JA, Enriquez GF et al. (2011) Factors affecting infestation by *Triatoma infestans* in a rural area of the humid Chaco in Argentina: a multi-model inference approach. PLoS Negl Trop. Drosophila Inf Service 5(10): e1349.

55. Buitrago R, Bosseno MF, Waleckx E, Brémond P, Vidaurre P et al. (2013) Risk of transmission of *Trypanosoma cruzi* by wild *Triatoma infestans* (Hemiptera: Reduviidae) in Bolivia supported by the detection of human blood meals. Infect Genet Evol 19: 141-144.

56. Canals M, Cruzat L, Molina MC, Ferreira A, Cattan PE (2001) Blood host sources of *Mepraia spinolai* (Heteroptera: Reduviidae), wild vector of chagas disease in Chile. J Med Entomol 38(2): 303-307. doi: 10.1603/0022-2585-38.2.303. PubMed: 11296839.

57. Brenière SF, Bosseno MF, Gastélum EM, Gutierrez Soto, de Jesus; Kasten Monges M et al. (2010) Community participation and domiciliary occurrence of infected *Meccus longipennis* in two Mexican villages in Jalisco state. Am J Trop Med Hyg 83(2): 382-387. doi:10.4269/ajtmh.2010.10-0080. PubMed: 20682887. Available online at: doi:10.4269/ajtmh.2010.10-0080 Available online at: PubMed: 20682887

58. Gürtler RE (2007) Combining residual insecticide spraying campaigns with targeted detection and specific chemotherapy for *Trypanosoma cruzi* infection in children. PLoS Negl Trop. Drosophila Inf Service1(3): e168.
59. Gaspe MS, Gurevitz JM, Gürtler RE, Dujardin JP (2013) Origins of house reinfestation with *Triatoma infestans* after insecticide spraying in the Argentine Chaco using wing geometric morphometry. *Infect Genet Evol* 17: 93-100. doi:10.1016/j.meegid.2013.03.044. PubMed: 23557838.

60. Abad-Franch F, Diotaiuti L, Gurgel-Gonçalves R, Gürtler RE (2013) Certifying the interruption of Chagas disease transmission by native vectors: cui bono? *Mem Inst Oswaldo Cruz* 108(2): 251-254. PubMed: 23579810.

61. Abad-Franch F, Vega MC, Rolón MS, Santos WS, Rojas de Arias A (2011) Community participation in Chagas disease vector surveillance: systematic review. *PLoS Negl Trop. Drosophila Inf Service* 5(6): e1207.

62. Picollo MI, Vassena CV, 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. Available online at: doi: 10.1603/0022-2585(2005)042[0637:HRTPIA]2.0.CO;2. PubMed: 16119553

63. 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. PubMed: 20545921.

64. Depickère S, Buitrago R, Siñani E, Baune M, Monje M et al. (2012) Susceptibility and resistance to deltamethrin of wild and domestic populations of *Triatoma infestans* (Reduviidae: Triatominae) in Bolivia: new discoveries. *Mem Inst Oswaldo Cruz* 107(8): 1042-1047. doi: 10.1590/S0074-02762012000800013. PubMed: 23295756.