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. Author summary: Understanding Chagas disease vector (triatomine) population dispersal is key for the design of control measures tailored for the epidemiological situation of a particular region. In Ecuador, Rhodnius ecuadoriensis is a cause of concern for Chagas disease transmission, since it is widely distributed from the central coast to southern Ecuador. Here, a genome-wide sequencing (2b-RAD) approach was performed in 20 specimens from four communities from Manabí (central coast) and Loja (southern) provinces of Ecuador, and the effectiveness of three type IIB restriction enzymes was assessed. The findings of this study show that this genotyping methodology is cost effective in R. ecuadoriensis and likely in other triatomine species. In addition, preliminary population genomic analysis results detected a signal of population structure among geographically distinct communities and genetic variability within communities. As such, 2b-RAD shows significant promise as a relatively low-tech solution for determination of vector population genomics, dynamics, and spread

Limitations of selective deltamethrin application for triatomine control in central coastal Ecuador

<p>Abstract</p> <p>Background</p> <p>This year-long study evaluated the effectiveness of a strategy involving selective deltamethrin spraying and community education for control of Chagas disease vectors in domestic units located in rural communities of coastal Ecuador.</p> <p>Results</p> <p>Surveys for triatomines revealed peridomestic infestation with <it>Rhodnius ecuadoriensis </it>and <it>Panstrongylus howardi</it>, with infestation indices remaining high during the study (13%, 17%, and 10%, at initial, 6-month, and 12-month visits, respectively), which indicates a limitation of this strategy for triatomine population control. Infestation was found 6 and 12 months after spraying with deltamethrin. In addition, a large number of previously vector-free domestic units also were found infested at the 6- and 12-month surveys, which indicates new infestations by sylvatic triatomines. The predominance of young nymphs and adults suggests new infestation events, likely from sylvatic foci. In addition, infection with <it>Trypanosoma cruzi </it>was found in 65%, 21% and 29% at initial, 6-month and 12-month visits, respectively. All parasites isolated (n = 20) were identified as TcI.</p> <p>Conclusion</p> <p>New vector control strategies need to be devised and evaluated for reduction of <it>T. cruzi </it>transmission in this region.</p

Culture-free genome-wide locus sequence typing (GLST) provides new perspectives on Trypanosoma cruzi dispersal and infection complexity.

Analysis of genetic polymorphism is a powerful tool for epidemiological surveillance and research. Powerful inference from pathogen genetic variation, however, is often restrained by limited access to representative target DNA, especially in the study of obligate parasitic species for which ex vivo culture is resource-intensive or bias-prone. Modern sequence capture methods enable pathogen genetic variation to be analyzed directly from host/vector material but are often too complex and expensive for resource-poor settings where infectious diseases prevail. This study proposes a simple, cost-effective 'genome-wide locus sequence typing' (GLST) tool based on massive parallel amplification of information hotspots throughout the target pathogen genome. The multiplexed polymerase chain reaction amplifies hundreds of different, user-defined genetic targets in a single reaction tube, and subsequent agarose gel-based clean-up and barcoding completes library preparation at under 4 USD per sample. Our study generates a flexible GLST primer panel design workflow for Trypanosoma cruzi, the parasitic agent of Chagas disease. We successfully apply our 203-target GLST panel to direct, culture-free metagenomic extracts from triatomine vectors containing a minimum of 3.69 pg/μl T. cruzi DNA and further elaborate on method performance by sequencing GLST libraries from T. cruzi reference clones representing discrete typing units (DTUs) TcI, TcIII, TcIV, TcV and TcVI. The 780 SNP sites we identify in the sample set repeatably distinguish parasites infecting sympatric vectors and detect correlations between genetic and geographic distances at regional (< 150 km) as well as continental scales. The markers also clearly separate TcI, TcIII, TcIV and TcV + TcVI and appear to distinguish multiclonal infections within TcI. We discuss the advantages, limitations and prospects of our method across a spectrum of epidemiological research

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

Ecological factors related to the widespread distribution of sylvatic Rhodnius ecuadoriensis populations in southern Ecuador

<p>Abstract</p> <p>Background</p> <p>Chagas disease transmission risk is a function of the presence of triatomines in domestic habitats. <it>Rhodnius ecuadoriensis </it>is one of the main vectors implicated in transmission of <it>Trypanosoma cruzi </it>in Ecuador. This triatomine species is present in domestic, peridomestic and sylvatic habitats in the country. To determine the distribution of sylvatic populations of <it>R. ecuadoriensis </it>and the factors related to this distribution, triatomine searches were conducted between 2005 and 2009 in southern Ecuador.</p> <p>Methods</p> <p>Manual triatomine searches were conducted by skilled bug collectors in 23 communities. Sylvatic searched sites were selected by a) directed sampling, where microhabitats were selected by the searchers and b) random sampling, where sampling points where randomly generated. Domiciliary triatomine searches were conducted using the one man-hour method. Natural trypanosome infection was determined by microscopic examination and PCR. Generalized linear models were used to test the effect of environmental factors on the presence of sylvatic triatomines.</p> <p>Results</p> <p>In total, 1,923 sylvatic individuals were collected representing a sampling effort of 751 man-hours. Collected sylvatic triatomines were associated with mammal and bird nests. The 1,219 sampled nests presented an infestation index of 11.9%, a crowding of 13 bugs per infested nest, and a colonization of 80% of the nests. Triatomine abundance was significantly higher in squirrel (<it>Sciurus stramineus</it>) nests located above five meters from ground level and close to the houses. In addition, 8.5% of the 820 examined houses in the same localities were infested with triatomines. There was a significant correlation between <it>R. ecuadoriensis </it>infestation rates found in sylvatic and synanthropic environments within communities (<it>p </it>= 0.012). Parasitological analysis revealed that 64.7% and 15.7% of the sylvatic bugs examined (n = 300) were infected with <it>Trypanosoma cruzi </it>and <it>T. rangeli </it>respectively, and 8% of the bugs presented mixed infections.</p> <p>Conclusions</p> <p>The wide distribution of sylvatic <it>R. ecuadoriensis </it>populations may jeopardize the effectiveness of control campaigns conducted to eliminate domestic populations of this species. Also, the high <it>T. cruzi </it>infection rates found in sylvatic <it>R. ecuadoriensis </it>populations in southern Ecuador could constitute a risk for house re-infestation and persistent long-term Chagas disease transmission in the region.</p

Sex, Subdivision, and Domestic Dispersal of Trypanosoma cruzi Lineage I in Southern Ecuador

Trypanosoma cruzi is transmitted by blood sucking insects known as triatomines. This protozoan parasite commonly infects wild and domestic mammals in South and Central America. However, triatomines also transmit the parasite to people, and human infection with T. cruzi is known as Chagas disease, a major public health concern in Latin America. Understanding the complex dynamics of parasite spread between wild and domestic environments is essential to design effective control measures to prevent the spread of Chagas disease. Here we describe T. cruzi genetic diversity and population dynamics in southern Ecuador. Our findings indicate that the parasite circulates in two largely independent cycles: one corresponding to the sylvatic environment and one related to the domestic/peridomestic environment. Furthermore, our data indicate that human activity might promote parasite dispersal among communties. This information is the key for the design of control programmes in Southern Ecuador. Finally, we have encountered evidence of a sexual reproductive mode in the domestic T. cruzi population, which constitutes a new and intriguing finding with regards to the biology of this parasite

Repeat-Driven Generation of Antigenic Diversity in a Major Human Pathogen, Trypanosoma cruzi.

Trypanosoma cruzi, a zoonotic kinetoplastid protozoan parasite, is the causative agent of American trypanosomiasis (Chagas disease). Having a very plastic, repetitive and complex genome, the parasite displays a highly diverse repertoire of surface molecules, with pivotal roles in cell invasion, immune evasion and pathogenesis. Before 2016, the complexity of the genomic regions containing these genes impaired the assembly of a genome at chromosomal level, making it impossible to study the structure and function of the several thousand repetitive genes encoding the surface molecules of the parasite. We here describe the genome assembly of the Sylvio X10/1 genome sequence, which since 2016 has been used as a reference genome sequence for T. cruzi clade I (TcI), produced using high coverage PacBio single-molecule sequencing. It was used to analyze deep Illumina sequence data from 34 T. cruzi TcI isolates and clones from different geographic locations, sample sources and clinical outcomes. Resolution of the surface molecule gene distribution showed the unusual duality in the organization of the parasite genome, a synteny of the core genomic region with related protozoa flanked by unique and highly plastic multigene family clusters encoding surface antigens. The presence of abundant interspersed retrotransposons in these multigene family clusters suggests that these elements are involved in a recombination mechanism for the generation of antigenic variation and evasion of the host immune response on these TcI strains. The comparative genomic analysis of the cohort of TcI strains revealed multiple cases of such recombination events involving surface molecule genes and has provided new insights into T. cruzi population structure

Trypanosoma cruzi population dynamics in the Central Ecuadorian Coast.

Chagas disease is the most important parasitic disease in Latin America. The causative agent, Trypanosoma cruzi, displays high genetic diversity and circulates in complex transmission cycles among domestic, peridomestic and sylvatic environments. In Ecuador, Rhodnius ecuadoriensis is known to be the major vector species implicated in T. cruzi transmission. However, across vast areas of Ecuador, little is known about T. cruzi genetic diversity in relation to different parasite transmission scenarios. Fifty-eight T. cruzi stocks from the central Ecuadorian coast, most of them derived from R. ecuadoriensis, were included in the study. All of them were genotyped as T. cruzi discrete typing unit I (DTU TcI). Analysis of 23 polymorphic microsatellite loci through neighbor joining and discriminant analysis of principal components yielded broadly congruent results and indicate genetic subdivision between sylvatic and peridomestic transmission cycles. However, both analyses also suggest that any barriers are imperfect and significant gene flow between parasite subpopulations in different habitats exists. Also consistent with moderate partition and residual gene flow between subpopulations, the fixation index (FST) was significant, but of low magnitude. Finally, the lack of private alleles in the domestic/peridomestic transmission cycle suggests the sylvatic strains constitute the ancestral population. The T. cruzi population in the central Ecuadorian coast shows moderate tendency to subdivision according to transmission cycle. However, connectivity between cycles exists and the sylvatic T. cruzi population harbored by R. ecuadoriensis vectors appears to constitute a source from which the parasite invades human domiciles and their surroundings in this region. We discuss the implications these findings have for the planning, implementation and evaluation of local Chagas disease control interventions

From roads to biobanks: roadkill animals as a valuable source of genetic data

To protect biodiversity we must understand its structure and composition including the bacteria and microparasites associated with wildlife, which may pose risks to human health. However, acquiring this knowledge often presents challenges, particularly in areas of high biodiversity where there are many undescribed and poorly studied species and funding resources can be limited. A solution to fill this knowledge gap is sampling roadkill (animals that die on roads as a result of collisions with circulating vehicles). These specimens can help characterize local wildlife and their associated parasites with fewer ethical and logistical challenges compared to traditional specimen collection. Here we test this approach by analyzing 817 tissue samples obtained from 590 roadkill vertebrate specimens (Amphibia, Reptilia, Aves and Mammalia) collected in roads within the Tropical Andes of Ecuador. First, we tested if the quantity and quality of recovered DNA varied across roadkill specimens collected at different times since death, exploring if decomposition affected the potential to identify vertebrate species and associated microorganisms. Second, we compared DNA stability across taxa and tissues to identify potential limitations and offer recommendations for future work. Finally, we illustrate how these samples can aid in taxonomic identification and parasite detection. Our study shows that sampling roadkill can help study biodiversity. DNA was recovered and amplified (allowing species identification and parasite detection) from roadkill even 120 hours after death, although risk of degradation increased overtime. DNA was extracted from all vertebrate classes but in smaller quantities and with lower quality from amphibians. We recommend sampling liver if possible as it produced the highest amounts of DNA (muscle produced the lowest). Additional testing of this approach in areas with different environmental and traffic conditions is needed, but our results show that sampling roadkill specimens can help detect and potentially monitor biodiversity and could be a valuable approach to create biobanks and preserve genetic data

2b-RAD raw data of Rhodnius ecuadoriensis

Sequencing A) Clustering was done by ‘onboard clustering’ and samples were sequenced on MiSeq (MSC 2.5.0.5/RTA 1.18.54) with a 1x51 setup using ‘Version2’ chemistry. The Bcl to FastQ conversion was performed using bcl2fastq-1.8.4 from the CASAVA software suite. The quality scale used is Sanger / phred33 / Illumina 1.8+