Microsatellite variation revealed panmictic pattern for Triatoma brasiliensis (Triatominae: Reduviidae) in rural northeastern Brazil: the control measures implications

Population genetics and genomics of Triatoma brasiliensis (Hemiptera, Reduviidae) in an area of high pressure of domiciliary infestation in Northeastern Brazil

Understanding the population dynamics of vectors is crucial for effective control of vector-borne diseases. In the Northeastern Brazilian semi-arid region, Triatoma brasiliensis persists as the most significant Chagas disease vector, frequently displaying recurrent domiciliary infestations. This situation raises relevant public health concerns in the municipality of Currais Novos in the state of Rio Grande do Norte. This area has experienced a high prevalence of peridomiciliary re-infestations by T. brasiliensis, coupled with elevated rates of Trypanosoma cruzi infection. Therefore, we assessed the distribution of genetic variation via mitochondrial Cytochrome b gene (MT-CYB) sequencing (n = 109) and single nucleotide polymorphisms (SNPs, n = 86) to assess the gene flow among distinct populations distributed in varied geographic spots and environments, mainly sylvatic and peridomiciliary. Insects were collected from rural communities at Currais Novos, enclosed within a 16 km radius. Sampling included 13 populations: one intradomiciliary, eight peridomiciliary, and four sylvatic. Furthermore, an external population located 220 km from Currais Novos was also included in the study. The method employed to obtain SNP information relied on ddRAD-seq genotyping-by-sequencing (GBS), enabling a genome-wide analysis to infer genetic variation. Through AMOVA analysis of MT-CYB gene variation, we identified four distinct population groups with statistical significance (FCT= 0.42; p<0.05). We identified a total of 3,013 SNPs through GBS, with 11 loci showing putative signs of being under selection. The variation based on 3,002 neutral loci evidenced low genetic structuration based on low FST values (p>0.05), indicating local panmixia. However, resampling algorithms pointed out that three samples from the external population were assigned (>98 %) in a cluster contrasting from the ones putatively under local panmixia - validating the newly applied genome-wide marker for studies on the population genetics at finer-scale resolution for T. brasiliensis. The presence of population structuring in some of the sampled points, as suggested by the mitochondrial marker, leads us to assume that infestations were probably initiated by small populations of females - demographic event poses a risk for rapid re-infestations. The local panmictic pattern revealed by the GBS marker poses a challenge for vector control measures, as re-infestation foci may be distributed over a wide geographical and ecological range. In such instances, vectors exhibit reduced susceptibility to conventional insecticide spraying operations since sylvatic populations are beyond the reach of these interventions. The pattern of infestation exhibited by T. brasiliensis necessitates integrating innovative strategies into the existing control framework, holding the potential to create a more resilient and adaptive vector control program. In our dataset, the results demonstrated that the genetic signals from both markers were complementary. Therefore, it is essential to consider the nature and inheritance pattern of each marker when inferring the pattern of re-infestations.

Experimental Hybrids of the Triatoma brasiliensis Species Complex Show Higher Susceptibility to the Trypanosoma cruzi Infection Than Their Parentals

The Triatoma brasiliensis species complex is a monophyletic group encompassing two subspecies and six species. Recently, a hybrid zone of members of this complex was recorded in the state of Pernambuco. Questions concerning the capability of the hybrids to become infected with Trypanosoma cruzi have been raised. This study aimed to compare the susceptibility of Triatoma b. brasiliensis, Triatoma juazeirensis, and their experimental hybrids to infection with T. cruzi. We infected the parentals and their experimental hybrids (obtained through reciprocal crosses) through artificial feeding with citrated rabbit blood, to which the TcI 0354 strain of T. cruzi had been added. The insects were weighed before and after feeding on the rabbit blood, and then they were dissected on the 10th, 20th, and 30th day after infection. Both the hybrids and the parentals remained infected throughout the experiment. The parasite was mostly found in the epimastigote form. The number of epimastigotes was significantly lower in the stomach and small intestine of T. juazeirensis than in the hybrids or in T. b. brasiliensis. A significantly higher percentage of metacyclic trypomastigotes was detected in the small intestine and rectum of the hybrids. Hybrids demonstrated higher susceptibility to the TcI 0354 strain than their parentals, opening up new avenues to be investigated.

Dynamics of Panstrongylus megistus infestation,the primary vector of Trypanosoma cruzi in Minas Gerais,Brazil

Panstrongylus megistus is considered one of the primary species of epidemiological importance for the transmission of Chagas disease in Brazil due to its wide geographical distribution throughout the national territory, good ability to invade and colonize houses, and high rates of natural infection by the Trypanosoma cruzi. The importance of this species in Minas Gerais has been recognized since the 80s. It is responsible for the high prevalence rates of Chagas disease in the west of the state. Studies conducted in the municipality of Jaboticatubas show that P. megistus is still the most captured vector in the region, even after 40 years of uninterrupted actions of the Chagas Disease Control Program in the municipality. Despite the importance of the species, its population dynamics is little studied. Consequently, crucial genetic information such as genetic diversity and gene flow among environments have not been well characterized yet. In this context, this work presents a population genetic analysis at the microgeographic level using microsatellite markers in samples of P. megistus obtained from Jaboticatubas to better understand the infestation dynamics of the primary vector species of T. cruzi in the region. The observed and expected heterozygosity ranged from 0.26 to 0.47 and 0.47 to 0.65, respectively. Most loci presented Hardy-Weinberg disequilibrium due to the excess of homozygotes. The pairwise Fst ranged from 0.05 to 0.35, with the p-value significant for all comparisons, indicating the absence of gene flow between them. The values of Fis found ranged from 0.25 to 0.52, all values of p ≤ 0.05, probably due to structured populations, inbreeding, or null alleles. The results suggest an extended stay of this species in the colonized environments, with rare dispersal to other locations. Such results differ from that observed for Triatoma brasiliensis and Triatoma dimidiata, species that present constant movement, even after colonization of the artificial environment, and are similar to the behavior of Triatoma infestans. Thus, the study reinforces the importance of continuous entomological surveillance in the areas of occurrence of P. megistus to avoid the installation and formation of new foci of household infestation in the region.

Morphological Stasis in Time? A Triatoma brasiliensis brasiliensis Study Using Geometric Morphometrics in the Long Run

Simple Summary

Triatomines are vector insects capable of transmitting the protozoan that causes Chagas disease, thus representing a health risk in several countries, especially in Central and South America. Triatoma brasiliensis brasiliensis, the main triatomine vector in northeastern Brazil, needs frequent monitoring as it is able to colonize various natural and artificial ecotopes as well as to infest domiciles. This research uses geometric morphometrics as a tool to evaluate changes in the morphology and analyze a large temporal dataset of 102 years of collections of T. b. brasiliensis. Finding novelty results which indicate a morphological stasis in the morphology.

Abstract

Triatoma brasiliensis brasiliensis Neiva, 1911 is one of the most important vectors of Chagas disease in the Brazilian semiarid regions in the north-east. The risk imposed by T. b. brasiliensis to the human populations, due to frequent invasions and/or colonization of the domiciles, demands constant monitoring and control actions as well as an understanding of its evolutionary process. In this context, the following research studies the pattern of shape adaptation over time using a large dataset from 102 years of specimen collections in order to identify the morphological plasticity of this vector in Brazil. This dataset was analyzed using geometric morphometrics tools and the timescale was divided into eight different groups, containing specimens from 1912 to 2014. Geometric morphometrics analysis showed an interesting morphological stasis in the wing shape of T. b. brasiliensis, which allowed us to understand the high capacity of adaptation to changes in climate condition through time, and the invasive status which Triatoma species have around the world. Moreover, these results showed novel findings as an interesting phenotypic pattern, with no modifications in more than 100 years, leading us to understand the shape evolution in Triatominae as a vector species of diseases.

Identification and characterization of microsatellite markers for population genetic studies of Panstrongylus megistus (Burmeister, 1835) (Triatominae: Reduviidae)

Background

Panstrongylus megistus is the most important vector of Chagas disease in Brazil. Studies show that the principal factor hindering the control of triatomines is reinfestation of houses previously treated with insecticides. Studies at the microgeographic level are therefore necessary to better understand these events. However, an efficient molecular marker is not yet available for carrying out such analyses in this species. The aim of the present study was to identify and characterize microsatellite loci for future population genetic studies of P. megistus.

Methods

This study work consisted of five stages: (i) sequencing of genomic DNA; (ii) assembly and selection of contigs containing microsatellites; (iii) validation of amplification and evaluation of polymorphic loci; (iv) standardization of the polymorphic loci; and (v) verification of cross-amplification with other triatomine species.

Results

Sequencing of males and females generated 7,908,463 contigs with a total length of 2,043,422,613 bp. A total of 2,043,690 regions with microsatellites in 1,441,091 contigs were obtained, with mononucleotide repeats being the most abundant class. From a panel of 96 loci it was possible to visualize polymorphisms in 64.55% of the loci. Of the 20 loci genotyped, the number of alleles varied from two to nine with an average of 4.9. Cross-amplification with other species of triatomines was observed in 13 of the loci.

Conclusions

Due to the high number of alleles encountered, polymorphism and the capacity to amplify from geographically distant populations, the microsatellites described here show promise for utilization in population genetic studies of P. megistus.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04771-w.