Nested clade and phylogeographic analyses of the Chagas disease vector Triatoma brasiliensis in Northeast Brazil

Triatoma brasiliensis (Hemiptera: Reduviidae: Triatominae) is the most important Chagas disease vector in the semiarid areas of Northeast Brazil. We analyzed mitochondrial cytochrome b sequence variation among 136 individuals representing 16 populations from across the species' distribution. Neighbor-joining and parsimony tree-building methods were used in conjunction with nested clade analysis to describe the systematics and phylogeography of this species. Our results indicate that T. brasiliensis is composed of four genetically distinct chromatic forms (referred to as brasiliensis, macromelasoma, juazeiro, and melanica) that present inter-population divergence values (0.027-0.119, corrected K2-p) and a pattern of haplotype geographic distribution compatible with the existence of a species complex. As a consequence, such forms can be treated as isolated targets in vector control programs. We were unable to infer what is shaping the population structure of the brasiliensis form as we obtained mutually exclusive causes of structure, namely a barrier to gene flow caused by past population fragmentation, and isolation by distance between populations (which would permit gene flow). We found indication of mitochondrial DNA introgression occurring among forms in putative hybrid zones.