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

Abstract

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.