Pioneer study of population genetics of Rhodnius ecuadoriensis (Hemiptera: Reduviidae) from the central coastand southern Andean regions of Ecuador

Bacterial microbiota from the gut of Rhodnius ecuadoriensis, a vector of Chagas disease in Ecuador's Central Coast and Southern Andes

Introduction

Chagas disease is a neglected tropical disease caused by the parasite Trypanosoma cruzi that is transmitted mainly by the feces of infected Triatomines. In Ecuador the main vector is Rhodnius ecuadoriensis which is distributed in several provinces of the country. More than 40% of these insects in the wild have T. cruzi as part of their intestinal microbiota. For this reason, the objective of this research was to characterize the intestinal bacterial microbiota of R. ecuadoriensis.

Methods

The methodology used was based on the DNA extraction of the intestinal contents from the wild collected insects (adults and nymphs V), as well as the insects maintained at the insectary of the CISeAL. Finally, the samples were analyzed by metagenomics extensions based on the different selected criteria.

Results

The intestinal microbiota of R. ecuadoriensis presented a marked divergence between laboratory-raised and wild collected insects. This difference was observed in all stages and was similar between insects from Loja and Manabí. A large loss of microbial symbionts was observed in laboratory-raised insects.

Discussion

This study is a crucial first step in investigating microbiota interactions and advancing new methodologies.

Evaluation of the Effectiveness of Chemical Control for Chagas Disease Vectors in Loja Province, Ecuador

The objective of this study was to evaluate the effectiveness of selective and community-wide house insecticide spraying in controlling triatomines in the subtropical areas of Loja Province, Ecuador. We designed a quasi-experimental pre-post-test without a control group to compare entomological levels before and after spraying. The baseline study was conducted in 2008. Second, third, and fourth visits were conducted in 2010, 2011, and 2012 in three rural communities. Out of the 130 domestic units (DU) visited, 41 domestic units were examined in each of the four visits. Selective and community-wide insecticide interventions included spraying with 5% deltamethrin at 25 mg/m2 active ingredient. At each visit, a questionnaire was administered to identify the characteristics of households, and DUs were searched for triatomine bugs. In addition, parasitological analysis was carried out in life triatomines. One and two rounds of selective insecticide spraying decreased the probability of infestation by 62% (pairwise odds ratios [POR] 0.38, 95% confidence interval [CI] 0.17-0.89, p = 0.024) and 51% (POR 0.49, 95% CI 0.23-1.01, p = 0.054), respectively. A similar effect was observed after one round of community-wide insecticide application in Chaquizhca and Guara (POR 0.55, CI 0.24-1.25, p = 0.155) and Bellamaria (POR 0.62, CI 0.22-1.79, p = 0.379); however, it was not statistically significant. Trypanosoma cruzi infection in triatomines (n = 483) increased overtime, from 2008 (42.9% and 8.5% for Rhodnius ecuadoriensis and Panstrongylus chinai, respectively) to 2012 (79.5% and 100%). Neither of the two spraying methodologies was effective for triatomine control in this area and our results point to a high likelihood of reinfestation after insecticide application. This underscores the importance of the implementation of physical barriers that prevent invasion and colonization of triatomines in households, such as home improvement initiatives, accompanied by a concerted effort to address the underlying socioeconomic issues that keep this population at risk of developing Chagas disease.

Population genomics and geographic dispersal in Chagas disease vectors: Landscape drivers and evidence of possible adaptation to the domestic setting

Accurate prediction of vectors dispersal, as well as identification of adaptations that allow blood-feeding vectors to thrive in built environments, are a basis for effective disease control. Here we adopted a landscape genomics approach to assay gene flow, possible local adaptation, and drivers of population structure in Rhodnius ecuadoriensis, an important vector of Chagas disease. We used a reduced-representation sequencing technique (2b-RADseq) to obtain 2,552 SNP markers across 272 R. ecuadoriensis samples from 25 collection sites in southern Ecuador. Evidence of high and directional gene flow between seven wild and domestic population pairs across our study site indicates insecticide-based control will be hindered by repeated re-infestation of houses from the forest. Preliminary genome scans across multiple population pairs revealed shared outlier loci potentially consistent with local adaptation to the domestic setting, which we mapped to genes involved with embryogenesis and saliva production. Landscape genomic models showed elevation is a key barrier to R. ecuadoriensis dispersal. Together our results shed early light on the genomic adaptation in triatomine vectors and facilitate vector control by predicting that spatially-targeted, proactive interventions would be more efficacious than current, reactive approaches.

Re-infestation of recently insecticide-treated houses by wild/secondary triatomine, their potential adaptation to this new environment and capabilities to geographically disperse across multiple human communities jeopardise sustainable Chagas disease control. This is the first study in Chagas disease vectors that identifies genomic regions possibly linked to adaptations to the built environment and describes landscape drivers for accurate prediction of geographic dispersal. We sampled multiple domestic and wild Rhodnius ecuadoriensis population pairs across a mountainous terrain in southern Ecuador. We evidenced that triatomine movement from forest to built enviroments does occur at a high rate. In these highly connected population pairs we detected loci possibly linked to local adaptation among the genomic makers we evaluated and in doing so we pave the way for future triatomine genomic research. We highlighted that current haphazardous vector control in the zone will be hindered by reinfestation of triatomines from the forest. Instead, we recommend frequent and spatially-targeted vector control and provided a landacape genomic model that identifies highly connected and isolated triatomine populations to facilitate efficient vector control.

Do the new triatomine species pose new challenges or strategies for monitoring Chagas disease? An overview from 1979-2021

Chagas disease persists as one of the most important, and yet most neglected, diseases in the world, and several changes in its epidemiological aspects have been recorded since its discovery. Currently, some of the most relevant changes are related to: (i) the reduction in the incidence of the endemic due to the control of the most important vectors, Triatoma infestans and Rhodnius prolixus, in many countries; (ii) the migration of human populations spreading cases of the disease throughout the world, from endemic to non-endemic areas, transforming Chagas disease into a global threat; and (iii) new acute cases and deaths caused by oral transmission, especially in the north of Brazil. Despite the reduction in the number of cases, new challenges need to be responded to, including monitoring and control activities aiming to prevent house infestation by the secondary vectors from occurring. In 1979, Lent & Wygodzinsky(1) published the most complete review of the subfamily Triatominae, encompassing 111 recognised species in the taxon. Forty-two years later, 46 new species and one subspecies have been described or revalidated. Here we summarise the new species and contextualise them regarding their ecology, epidemiologic importance, and the obstacles they pose to the control of Chagas disease around the world.

Under pressure: phenotypic divergence and convergence associated with microhabitat adaptations in Triatominae

Background

Triatomine bugs, the vectors of Chagas disease, associate with vertebrate hosts in highly diverse ecotopes. It has been proposed that occupation of new microhabitats may trigger selection for distinct phenotypic variants in these blood-sucking bugs. Although understanding phenotypic variation is key to the study of adaptive evolution and central to phenotype-based taxonomy, the drivers of phenotypic change and diversity in triatomines remain poorly understood.

Methods/results

We combined a detailed phenotypic appraisal (including morphology and morphometrics) with mitochondrial cytb and nuclear ITS2 DNA sequence analyses to study Rhodnius ecuadoriensis populations from across the species’ range. We found three major, naked-eye phenotypic variants. Southern-Andean bugs primarily from vertebrate-nest microhabitats (Ecuador/Peru) are typical, light-colored, small bugs with short heads/wings. Northern-Andean bugs from wet-forest palms (Ecuador) are dark, large bugs with long heads/wings. Finally, northern-lowland bugs primarily from dry-forest palms (Ecuador) are light-colored and medium-sized. Wing and (size-free) head shapes are similar across Ecuadorian populations, regardless of habitat or phenotype, but distinct in Peruvian bugs. Bayesian phylogenetic and multispecies-coalescent DNA sequence analyses strongly suggest that Ecuadorian and Peruvian populations are two independently evolving lineages, with little within-lineage phylogeographic structuring or differentiation.

Conclusions

We report sharp naked-eye phenotypic divergence of genetically similar Ecuadorian R. ecuadoriensis (nest-dwelling southern-Andean vs palm-dwelling northern bugs; and palm-dwelling Andean vs lowland), and sharp naked-eye phenotypic similarity of typical, yet genetically distinct, southern-Andean bugs primarily from vertebrate-nest (but not palm) microhabitats. This remarkable phenotypic diversity within a single nominal species likely stems from microhabitat adaptations possibly involving predator-driven selection (yielding substrate-matching camouflage coloration) and a shift from palm-crown to vertebrate-nest microhabitats (yielding smaller bodies and shorter and stouter heads). These findings shed new light on the origins of phenotypic diversity in triatomines, warn against excess reliance on phenotype-based triatomine-bug taxonomy, and confirm the Triatominae as an informative model system for the study of phenotypic change under ecological pressure.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04647-z.

Would tropical climatic variations impact the genetic variability of triatomines: Rhodnius ecuadoriensis, principal vector of Chagas disease in Ecuador?

Rhodnius ecuadoriensis is one of the most important vector species of Chagas disease in Ecuador. This species is distributed in the Central coast region and in the south Andean region, and an incipient speciation process between these geographical populations was previously proposed. The current population genetics study only focused on the Central coast region and analyzed 96 sylvatic specimens of R. ecuadoriensis associated with Phytelephas aequatorialis palm trees. We used Cytb and 16S-rRNA sequences and a Cytb-16S-rRNA concatenated set to explore (i) the genetic variability, spatial structuring, and demographic history of R. ecuadoriensis, and to determine (ii) the relationship between the genetic and climatic variabilities. A particularly high genetic variability was observed without detectable general genetic structure; only some terminal genetic clusters were observed. We did not observe isolation by geographical distance (IBD), and it is likely that ancient expansion occurred, according to Fs index and mismatch distribution for Cytb-16S-rRNA concatenated sequences. Hierarchical clustering showed that the current locality origins of the bugs were grouped into four bioclimatic clusters. Genetic and bioclimatic distances were not correlated, but some genetic clusters were associated with bioclimatic ones. The results showed an ancient evolution of the species in the region with a possible old expansion. The absence of spatial genetic structure could be due to climatic conditions (possible selection of singular genotypes) and to passive transportation of palms tree materials where R. ecuadoriensis are living.

Population genetics of two chromatic morphs of the Chagas disease vector Rhodnius pallescens Barber, 1932 in Panamá

Rhodnius pallescens is the principal vector of Chagas disease in Panama. Recently a dark chromatic morph has been discovered in the highlands of Veraguas Province. Limited genetic studies have been conducted with regards to the population structure and dispersal potential of Triatominae vectors, particularly in R. pallescens. Next generation sequencing methods such as RADseq and complete mitochondrial DNA (mtDNA) genome sequencing have great potential for examining vector biology across space and time. Here we utilize a RADseq method (3RAD), along with complete mtDNA sequencing, to examine the population structure of the two chromatic morpho types of R. pallescens in Panama. We sequenced 105 R. pallescens samples from five localities in Panama. We generated a 2216 SNP dataset and 6 complete mtDNA genomes. RADseq showed significant differentiation among the five localities (F_CT = 0.695; P = .004), but most of this was between localities with the dark vs. light chromatic morphs (Veraguas vs. Panama Oeste). The mtDNA genomes showed a 97-98% similarity between dark and light chromatic morphs across all genes and a 502 bp insert in light morphs. Thus, both the RADseq and mtDNA data showed highly differentiated clades with essentially no gene flow between the dark and light chromatic morphs from Veraguas and central Panama respectively. We discuss the growing evidence showing clear distinctions between these two morpho types with the possibility that these are separate species, an area of research that requires further investigation. Finally, we discuss the cost-effectiveness of 3RAD which is a third of the cost compared to other RADseq methods used recently in Chagas disease vector research.

Remarkable genetic diversity of Trypanosoma cruzi and Trypanosoma rangeli in two localities of southern Ecuador identified via deep sequencing of mini-exon gene amplicons

BACKGROUND

Trypanosoma cruzi, the causative agent of Chagas disease, and T. rangeli are kinetoplastid parasites endemic to Latin America. Although closely related to T. cruzi and capable of infecting humans, T. rangeli is non-pathogenic. Both parasite species are transmitted by triatomine bugs, and the presence of T. rangeli constitutes a confounding factor in the study of Chagas disease prevalence and transmission dynamics. Trypanosoma cruzi possesses high molecular heterogeneity: seven discrete typing units (DTUs) are currently recognized. In Ecuador, T. cruzi TcI and T. rangeli KP1(-) predominate, while other genetic lineages are seldom reported.

METHODS

Infection by T. cruzi and/or T. rangeli in different developmental stages of triatomine bugs from two communities of southern Ecuador was evaluated via polymerase chain reaction product size polymorphism of kinetoplast minicircle sequences and the non-transcribed spacer region of the mini-exon gene (n = 48). Forty-three mini-exon amplicons were also deep sequenced to analyze single-nucleotide polymorphisms within single and mixed infections. Mini-exon products from ten monoclonal reference strains were included as controls.

RESULTS

Trypanosoma cruzi genetic richness and diversity was not significantly greater in adult vectors than in nymphal stages III and V. In contrast, instar V individuals showed significantly higher T. rangeli richness when compared with other developmental stages. Among infected triatomines, deep sequencing revealed one T. rangeli infection (3%), 8 T. cruzi infections (23.5%) and 25 T. cruzi + T. rangeli co-infections (73.5%), suggesting that T. rangeli prevalence has been largely underestimated in the region. Furthermore, deep sequencing detected TcIV sequences in nine samples; this DTU had not previously been reported in Loja Province.

CONCLUSIONS

Our data indicate that deep sequencing allows for better parasite identification/typing than amplicon size analysis alone for mixed infections containing both T. cruzi and T. rangeli, or when multiple T. cruzi DTUs are present. Additionally, our analysis showed extensive overlap among the parasite populations present in the two studied localities (c.28 km apart), suggesting active parasite dispersal over the study area. Our results highlight the value of amplicon sequencing methodologies to clarify the population dynamics of kinetoplastid parasites in endemic regions and inform control campaigns in southern Ecuador.

Chagas vectors Panstrongylus chinai (Del Ponte, 1929) and Panstrongylus howardi (Neiva, 1911): chromatic forms or true species?

Background

Chagas disease is a parasitic infection transmitted by “kissing bugs” (Hemiptera: Reduviidae: Triatominae) that has a huge economic impact in Latin American countries. The vector species with the upmost epidemiological importance in Ecuador are Rhodnius ecuadoriensis (Lent & Leon, 1958) and Triatoma dimidiata (Latreille, 1811). However, other species such as Panstrongylus howardi (Neiva, 1911) and Panstrongylus chinai (Del Ponte, 1929) act as secondary vectors due to their growing adaptation to domestic structures and their ability to transmit the parasite to humans. The latter two taxa are distributed in two different regions, they are allopatric and differ mainly by their general color. Their relative morphological similarity led some authors to suspect that P. chinai is a melanic form of P. howardi.

Methods

The present study explored this question using different approaches: antennal phenotype; geometric morphometrics of heads, wings and eggs; cytogenetics; molecular genetics; experimental crosses; and ecological niche modeling.

Results

The antennal morphology, geometric morphometrics of head and wing shape and cytogenetic analysis were unable to show distinct differences between the two taxa. However, geometric morphometrics of the eggs, molecular genetics, ecological niche modeling and experimental crosses including chromosomal analyses of the F1 hybrids, in addition to their coloration and current distribution support the hypothesis that P. chinai and P. howardi are separate species.

Conclusions

Based on the evidence provided here, P. howardi and P. chinai should not be synonymized. They represent two valid, closely related species.

2b-RAD genotyping for population genomic studies of Chagas disease vectors: Rhodnius ecuadoriensis in Ecuador

BACKGROUND

Rhodnius ecuadoriensis is the main triatomine vector of Chagas disease, American trypanosomiasis, in Southern Ecuador and Northern Peru. Genomic approaches and next generation sequencing technologies have become powerful tools for investigating population diversity and structure which is a key consideration for vector control. Here we assess the effectiveness of three different 2b restriction site-associated DNA (2b-RAD) genotyping strategies in R. ecuadoriensis to provide sufficient genomic resolution to tease apart microevolutionary processes and undertake some pilot population genomic analyses.

METHODOLOGY/PRINCIPAL FINDINGS

The 2b-RAD protocol was carried out in-house at a non-specialized laboratory using 20 R. ecuadoriensis adults collected from the central coast and southern Andean region of Ecuador, from June 2006 to July 2013. 2b-RAD sequencing data was performed on an Illumina MiSeq instrument and analyzed with the STACKS de novo pipeline for loci assembly and Single Nucleotide Polymorphism (SNP) discovery. Preliminary population genomic analyses (global AMOVA and Bayesian clustering) were implemented. Our results showed that the 2b-RAD genotyping protocol is effective for R. ecuadoriensis and likely for other triatomine species. However, only BcgI and CspCI restriction enzymes provided a number of markers suitable for population genomic analysis at the read depth we generated. Our preliminary genomic analyses detected a signal of genetic structuring across the study area.

CONCLUSIONS/SIGNIFICANCE

Our findings suggest that 2b-RAD genotyping is both a cost effective and methodologically simple approach for generating high resolution genomic data for Chagas disease vectors with the power to distinguish between different vector populations at epidemiologically relevant scales. As such, 2b-RAD represents a powerful tool in the hands of medical entomologists with limited access to specialized molecular biological equipment.