Polymorphic microsatellite markers for the tsetse fly Glossina fuscipes fuscipes (Diptera: Glossinidae), a vector of human African trypanosomiasis

Development and characterization of microsatellite markers for the tsetse species Glossina brevipalpis and preliminary population genetics analyses

Tsetse flies, the vectors of African trypanosomes are of key medical and economic importance and one of the constraints for the development of Africa. Tsetse fly control is one of the most effective and sustainable strategies used for controlling the disease. Knowledge about population structure and level of gene flow between neighbouring populations of the target vector is of high importance to develop appropriate strategies for implementing effective management programmes. Microsatellites are commonly used to identify population structure and assess dispersal of the target populations and have been developed for several tsetse species but were lacking for Glossina brevipalpis. In this study, we screened the genome of G. brevipalpis to search for suitable microsatellite markers and nine were found to be efficient enough to distinguish between different tsetse populations. The availability of these novel microsatellite loci will help to better understand the population biology of G. brevipalpis and to assess the level of gene flow between different populations. Such information will help with the development of appropriate strategies to implement the sterile insect technique (SIT) in the framework of an area-wide integrated pest management (AW-IPM) approach to manage tsetse populations and ultimately address the trypanosomoses problem in these targeted areas.

Population Structure and Migration Patterns of the Tsetse Fly Glossina fuscipes in Congo-Brazzaville

Tsetse flies of the palpalis group, particularly Glossina fuscipes, are the main vectors of human African trypanosomiasis or sleeping sickness in Congo-Brazzaville. They transmit the deadly human parasite, Trypanosoma brucei gambiense and other trypanosomes that cause animal trypanosomiasis. Knowledge on diversity, population structure, population size, and gene flow is a prerequisite for designing effective tsetse control strategies. There is limited published information on these parameters including migration patterns of G. fuscipes in Congo-Brazzaville. We genotyped 288 samples of G. fuscipes from Bomassa (BMSA), Bouemba (BEMB), and Talangai (TLG) locations at 10 microsatellite loci and determined levels of genetic diversity, differentiation, structuring, and gene flow among populations. We observed high genetic diversity in all three localities. Mean expected heterozygosity was 0.77 ± 0.04, and mean allelic richness was 11.2 ± 1.35. Deficiency of heterozygosity was observed in all populations with positive and significant F IS values (0.077-0.149). Structure analysis revealed three clusters with genetic admixtures, evidence of closely related but potentially different taxa within G. fuscipes. Genetic differentiation indices were low but significant (F ST = 0.049, P < 0.05), indicating ongoing gene flow countered with a stronger force of drift. We recorded significant migration from all the three populations, suggesting exchange of genetic information between and among locations. Ne estimates revealed high and infinite population sizes in BEMB and TLG. These critical factors should be considered when planning area-wide tsetse control interventions in the country to prevent resurgence of tsetse from relict populations and/or reinvasion of cleared habitats.

Developing and quality testing of microsatellite loci for four species of Glossina

Microsatellite loci still represent valuable resources for the study of the population biology of non-model organisms. Discovering or adapting new suitable microsatellite markers in species of interest still represents a useful task, especially so for non-model organisms as tsetse flies (genus Glossina), which remain a serious threat to the health of humans and animals in sub-Saharan Africa. In this paper, we present the development of new microsatellite loci for four species of Glossina: two from the Morsitans group, G. morsitans morsitans (Gmm) from Zimbabwe, G. pallidipes (Gpalli) from Tanzania; and the other two from the Palpalis group, G. fuscipes fuscipes (Gff) from Chad, and G. palpalis gambiensis (Gpg) from Guinea. We found frequent short allele dominance and null alleles. Stuttering could also be found and amended when possible. Cryptic species seemed to occur frequently in all taxa but Gff. This explains why it may be difficult finding ecumenical primers, which thus need adaptation according to each taxonomic and geographic context. Amplification problems occurred more often in published old markers, and Gmm and Gpg were the most affected (stronger heterozygote deficits). Trinucleotide markers displayed selection signature in some instances (Gmm). Combining old and new loci, for Gmm, eight loci can be safely used (with correction for null alleles); and five seem particularly promising; for Gpalli, only five to three loci worked well, depending on the clade, which means that the use of loci from other species (four morsitans loci seemed to work well), or other new primers will need to be used; for Gff, 14 loci behaved well, but with null alleles, seven of which worked very well; and for G. palpalis sl, only four loci, needing null allele and stuttering corrections seem to work well, and other loci from the literature are thus needed, including X-linked markers, five of which seem to work rather well (in females only), but new markers will probably be needed. Finally, the high proportion of X-linked markers (around 30%) was explained by the non-Y DNA quantity and chromosome structure of tsetse flies studied so far.

The population genomics of multiple tsetse fly (Glossina fuscipes fuscipes) admixture zones in Uganda

Understanding the mechanisms that enforce, maintain or reverse the process of speciation is an important challenge in evolutionary biology. This study investigates the patterns of divergence and discusses the processes that form and maintain divergent lineages of the tsetse fly Glossina fuscipes fuscipes in Uganda. We sampled 251 flies from 18 sites spanning known genetic lineages and the four admixture zones between them. We apply population genomics, hybrid zone and approximate Bayesian computation to the analysis of three types of genetic markers: 55,267 double-digest restriction site-associated DNA (ddRAD) SNPs to assess genome-wide admixture, 16 microsatellites to provide continuity with published data and accurate biogeographic modelling, and a 491-bp fragment of mitochondrial cytochrome oxidase I and II to infer maternal inheritance patterns. Admixture zones correspond with regions impacted by the reorganization of Uganda's river networks that occurred during the formation of the West African Rift system over the last several hundred thousand years. Because tsetse fly population distributions are defined by rivers, admixture zones likely represent both old and new regions of secondary contact. Our results indicate that older hybrid zones contain mostly parental types, while younger zones contain variable hybrid types resulting from multiple generations of interbreeding. These findings suggest that reproductive barriers are nearly complete in the older admixture zones, while nearly absent in the younger admixture zones. Findings are consistent with predictions of hybrid zone theory: Populations in zones of secondary contact transition rapidly from early to late stages of speciation or collapse all together.

Nuclear and Wolbachia-based multimarker approach for the rapid and accurate identification of tsetse species

Background

Tsetse flies (Diptera: Glossinidae) are solely responsible for the transmission of African trypanosomes, causative agents of sleeping sickness in humans and nagana in livestock. Due to the lack of efficient vaccines and the emergence of drug resistance, vector control approaches such as the sterile insect technique (SIT), remain the most effective way to control disease. SIT is a species-specific approach and therefore requires accurate identification of natural pest populations at the species level. However, the presence of morphologically similar species (species complexes and sub-species) in tsetse flies challenges the successful implementation of SIT-based population control.

Results

In this study, we evaluate different molecular tools that can be applied for the delimitation of different Glossina species using tsetse samples derived from laboratory colonies, natural populations and museum specimens. The use of mitochondrial markers, nuclear markers (including internal transcribed spacer 1 (ITS1) and different microsatellites), and bacterial symbiotic markers (Wolbachia infection status) in combination with relatively inexpensive techniques such as PCR, agarose gel electrophoresis, and to some extent sequencing provided a rapid, cost effective, and accurate identification of several tsetse species.

Conclusions

The effectiveness of SIT benefits from the fine resolution of species limits in nature. The present study supports the quick identification of large samples using simple and cost effective universalized protocols, which can be easily applied by countries/laboratories with limited resources and expertise.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1295-4) contains supplementary material, which is available to authorized users.

Cryptic diversity within the major trypanosomiasis vector Glossina fuscipes revealed by molecular markers

BACKGROUND

The tsetse fly Glossina fuscipes s.l. is responsible for the transmission of approximately 90% of cases of human African trypanosomiasis (HAT) or sleeping sickness. Three G. fuscipes subspecies have been described, primarily based upon subtle differences in the morphology of their genitalia. Here we describe a study conducted across the range of this important vector to determine whether molecular evidence generated from nuclear DNA (microsatellites and gene sequence information), mitochondrial DNA and symbiont DNA support the existence of these taxa as discrete taxonomic units.

PRINCIPAL FINDINGS

The nuclear ribosomal Internal transcribed spacer 1 (ITS1) provided support for the three subspecies. However nuclear and mitochondrial sequence data did not support the monophyly of the morphological subspecies G. f. fuscipes or G. f. quanzensis. Instead, the most strongly supported monophyletic group was comprised of flies sampled from Ethiopia. Maternally inherited loci (mtDNA and symbiont) also suggested monophyly of a group from Lake Victoria basin and Tanzania, but this group was not supported by nuclear loci, suggesting different histories of these markers. Microsatellite data confirmed strong structuring across the range of G. fuscipes s.l., and was useful for deriving the interrelationship of closely related populations.

CONCLUSION/SIGNIFICANCE

We propose that the morphological classification alone is not used to classify populations of G. fuscipes for control purposes. The Ethiopian population, which is scheduled to be the target of a sterile insect release (SIT) programme, was notably discrete. From a programmatic perspective this may be both positive, given that it may reflect limited migration into the area or negative if the high levels of differentiation are also reflected in reproductive isolation between this population and the flies to be used in the release programme.

Polyandry is a common event in wild populations of the Tsetse fly Glossina fuscipes fuscipes and may impact population reduction measures

BACKGROUND

Glossina fuscipes fuscipes is the main vector of human and animal trypanosomiasis in Africa, particularly in Uganda. Attempts to control/eradicate this species using biological methods require knowledge of its reproductive biology. An important aspect is the number of times a female mates in the wild as this influences the effective population size and may constitute a critical factor in determining the success of control methods. To date, polyandry in G.f. fuscipes has not been investigated in the laboratory or in the wild. Interest in assessing the presence of remating in Ugandan populations is driven by the fact that eradication of this species is at the planning stage in this country.

METHODOLOGY/PRINCIPAL FINDINGS

Two well established populations, Kabukanga in the West and Buvuma Island in Lake Victoria, were sampled to assess the presence and frequency of female remating. Six informative microsatellite loci were used to estimate the number of matings per female by genotyping sperm preserved in the female spermathecae. The direct count of the minimum number of males that transferred sperm to the spermathecae was compared to Maximum Likelihood and Bayesian probability estimates. The three estimates provided evidence that remating is common in the populations but the frequency is substantially different: 57% in Kabukanga and 33% in Buvuma.

CONCLUSIONS/SIGNIFICANCE

The presence of remating, with females maintaining sperm from different mates, may constitute a critical factor in cases of re-infestation of cleared areas and/or of residual populations. Remating may enhance the reproductive potential of re-invading propagules in terms of their effective population size. We suggest that population age structure may influence remating frequency. Considering the seasonal demographic changes that this fly undergoes during the dry and wet seasons, control programmes based on SIT should release large numbers of sterile males, even in residual surviving target populations, in the dry season.

Phylogeography and population structure of Glossina fuscipes fuscipes in Uganda: implications for control of tsetse

Background

Glossina fuscipes fuscipes, a riverine species of tsetse, is the main vector of both human and animal trypanosomiasis in Uganda. Successful implementation of vector control will require establishing an appropriate geographical scale for these activities. Population genetics can help to resolve this issue by characterizing the extent of linkage among apparently isolated groups of tsetse.

Methodology/Principal Findings

We conducted genetic analyses on mitochondrial and microsatellite data accumulated from approximately 1000 individual tsetse captured in Uganda and neighboring regions of Kenya and Sudan. Phylogeographic analyses suggested that the largest scale genetic structure in G. f. fuscipes arose from an historical event that divided two divergent mitochondrial lineages. These lineages are currently partitioned to northern and southern Uganda and co-occur only in a narrow zone of contact extending across central Uganda. Bayesian assignment tests, which provided evidence for admixture between northern and southern flies at the zone of contact and evidence for northerly gene flow across the zone of contact, indicated that this structure may be impermanent. On the other hand, microsatellite structure within the southern lineage indicated that gene flow is currently limited between populations in western and southeastern Uganda. Within regions, the average F_ST between populations separated by less than 100 km was less than ∼0.1. Significant tests of isolation by distance suggested that gene flow is ongoing between neighboring populations and that island populations are not uniformly more isolated than mainland populations.

Conclusions/Significance

Despite the presence of population structure arising from historical colonization events, our results have revealed strong signals of current gene flow within regions that should be accounted for when planning tsetse control in Uganda. Populations in southeastern Uganda appeared to receive little gene flow from populations in western or northern Uganda, supporting the feasibility of area wide control in the Lake Victoria region by the Pan African Tsetse and Trypanosomiasis Eradication Campaign.

Glossina fuscipes fuscipes is the most common species of tsetse in Uganda, where it transmits human sleeping sickness and nagana, a related disease of cattle. A consortium of African countries dedicated to controlling these diseases is poised to begin area wide control of tsetse, but a critical question remains: What is the most appropriate geographical scale for these activities? To address this question, we used population genetics to determine the extent of linkage between populations of tsetse confined to discrete patches of riverine habitat. Our results suggest that Uganda was colonized by two distinct lineages of G. f. fuscipes, which now co-occur only in a narrow band across central Uganda. Evidence for interbreeding at the zone of contact and movement of genes from the south to the north suggest that this historical genetic structure may dissolve in the future. At smaller scales, we have demonstrated that exchange of genes among neighboring populations via dispersal is at equilibrium with the differentiating force of genetic drift. Our results highlight the need for investment in vector control programs that account for the linkage observed among tsetse populations. Given its genetic isolation and its location at the far edge of G. fuscipes' range, the Lake Victoria region appears to be an appropriate target for area wide control.