Skip to main content

Mitochondrial DNA sequence divergence and diversity of Glossina fuscipes fuscipes in the Lake Victoria basin of Uganda: implications for control



Glossina fuscipes fuscipes is the main vector of African Trypanosomiasis affecting both humans and livestock in Uganda. The human disease (sleeping sickness) manifests itself in two forms: acute and chronic. The Lake Victoria basin in Uganda has the acute form and a history of tsetse re-emergence despite concerted efforts to control tsetse. The government of Uganda has targeted the basin for tsetse eradication. To provide empirical data for this initiative, we screened tsetse flies from the basin for genetic variation at the mitochondrial DNA cytochrome oxidase II (mtDNA COII) gene with the goal of investigating genetic diversity and gene flow among tsetse, tsetse demographic history; and compare these results with results from a previous study based on microsatellite loci data in the same area.


We collected 429 Gff tsetse fly samples from 14 localities in the entire Ugandan portion of the Lake Victoria coast, covering 40,000 km2. We performed genetic analyses on them and added data collected for 56 Gff individuals from 4 additional sampling sites in the basin. The 529pb partial mitochondrial DNA cytochrome oxidase II (mtDNA COII) sequences totaling 485 were analysed for genetic differentiation, structuring and demographic history. The results were compared with findings from a previous study based on microsatellite loci data from the basin.


The differences within sampling sites explained a significant proportion of the genetic variation. We found three very closely related mtDNA population clusters, which co-occurred in multiple sites. Although Φ ST (0 – 0.592; P < 0.05) and Bayesian analyses suggest some level of weak genetic differentiation, there is no correlation between genetic divergence and geographic distance (r = 0.109, P = 0.185), and demographic tests provide evidence of locality-based demographic history.


The mtDNA data analysed here complement inferences made in a previous study based on microsatellite data. Given the differences in mutation rates, mtDNA afforded a look further back in time than microsatellites and revealed that Gff populations were more connected in the past. Microsatellite data revealed more genetic structuring than mtDNA. The differences in connectedness and structuring over time could be related to vector control efforts. Tsetse re-emergence after control interventions may be due to re-invasions from outside the treated areas, which emphasizes the need for an integrated area-wide tsetse eradication strategy for sustainable removal of the tsetse and trypanosomiasis problem from this area.


Tsetse flies (Diptera: Glossinidae) are the major vectors of Human African Trypanosomiasis (HAT) and Animal African Trypanosomoses (AAT) in sub-Saharan Africa [1, 2]. Approximately 70 million people in 1.55 million km2 are estimated to be at risk of HAT caused by two species of trypanosomes [3]: Trypanosoma brucei gambiense (Tbg), responsible for the chronic form of the disease, and Trypanosoma brucei rhodesiense (Tbr), which causes the acute form [4, 5]. There is evidence that tsetse have influenced food production, urbanization, and institutional development dating back to historical Africa [6]. AAT is a major obstacle to the development of more efficient and sustainable livestock production systems, and thus one of the most important causes of hunger and poverty [7, 8]. There are currently no vaccines for the above diseases, and the available drugs are expensive, toxic, and logistically difficult to administer.

Since reducing host/vector contact can rapidly slow human trypanosomiasis transmission [9], controlling the tsetse fly remains the most efficient and sustainable way of managing African trypanosomiasis. Available environmentally-friendly tsetse control techniques include the sequential aerosol technique (SAT), which is an aerial application of ultra-low-volume non-residual insecticides [10], the use of insecticide-impregnated targets and traps that can be odour-baited [11], the application of residual insecticides on livestock, referred to as the live bait technique [12], and the sterile insect technique (SIT) [13].

In 2001, the African Union established the Pan African Tsetse and Trypanosomiasis Eradication Campaign (PATTEC) with a view of using an integrated area-wide approach to control HAT and AAT with the available methods. A prerequisite to any vector control campaign aiming at eradication is to identify and target isolated populations to minimize the risk of reinvasion. If not already isolated, populations could be isolated by creating physical obstacles, such as the insecticide-impregnated biconical trap barriers. Such a method has been used to effectively control Glossina palpalis gambiensis and Glossina tachinoides in a 3000 km2 area in an agro-pastoral zone of Sideradougou, in the Guinea savannah in Cameroon [14].

Population genetic techniques can help understand and quantify gene flow between populations, which can be used as a proxy for dispersal [15]. Dispersal rates for Glossina fuscipes fuscipes (Gff) based on mark-release-recapture (MRR) studies are about 14.2 km per generation, given a movement estimate of 338 m/day [9].

Fine-scale genetic analysis based on microsatellites confirmed that Gff disperse up to 14 km per generation [16], Gff appear to be genetically homogeneous over 1–5 km2.

Information about dispersal derived using population genetic techniques can be used to support vector control decision-making [17, 18] at various spatial levels and ecological settings. For example, regional studies such as the one on riverine Glossina palpalis palpalis in west and central Africa [19] have provided information that is useful for control of riverine palpalis tsetse group in cross-boundary projects. Studies of tsetse in Burkina Faso, Guinea and Senegal have identified riverine tsetse populations that are sufficiently isolated to warrant attempts at complete eradication [20, 21]. In the morsitans or savannah tsetse group, population genetic studies have indicated high gene flow among Glossina morsitans morsitans populations separated at geographic scales of 12–917 km in East and Southern Africa [22, 23].

In Uganda, Gff, a riverine subspecies in the palpalis group, is the major vector of HAT. The acute form of HAT (T. b. rhodesiense) previously had its historical focus along the shores of Lake Victoria, but has recently extended its range northwards into central Uganda [4, 24]. If this distribution continues extending, the range might overlap with that of the chronic form of HAT (T. b. gambiense) found in northwestern Uganda, thereby complicating diagnosis, treatment, and providing new challenges, as recombination between the two trypanosome forms can occur and could lead to unforeseen pathologies [25, 26].

In an effort to eliminate the acute form of the disease and to prevent potential challenges associated with overlap of the two forms of HAT in Uganda, in 2008 Pan African Tsetse and Trypanosomiasis Eradication Campaign (PATTEC) activities were initiated against Gff in the Lake Victoria basin; an area with a history of tsetse re-emergence despite concerted tsetse control efforts [27]. Tsetse re-emergence is a major obstacle to elimination of the tsetse fly vector in Africa [28]. Understanding the population genetics of Gff in the Lake Victoria basin may elucidate the factors influencing re-emergence. Indeed, genetic tools have revealed genetic structuring among localities north, south and west of Lake Kyoga in Uganda, occurrence of gene flow among genetic clusters [29], and temporal stability of these genetic patterns [30]. We previously screened for genetic variation at 15 microsatellite loci using tsetse flies from 14 sampling sites from continental and island locations along Lake Victoria in Uganda [16]. That study identified four genetically distinct clusters and showed that gene flow occurred at varying levels between these clusters.

In this study, we followed up on the work of our group [16] by screening 485 tsetse flies from 18 sampling sites (Fig. 1) for genetic variation in a fragment of the mtDNA COII gene (526 bp). In contrast to the bi-parentally inherited microsatellites, mtDNA is maternally inherited and lacks recombination [31, 32]. Given these differences, as well as the slower mutation rate in mtDNA than microsatellites [33, 34], we can compare differences in genetic variation among different timescales. The insight about temporal dynamics that the comparison of mtDNA and microsatellite data affords, could further inform the ongoing PATTEC control and monitoring efforts in the area and possibly beyond.

Fig. 1

Map showing the location of the 18 sampling sites and the distribution of the 23 COII mtDNA haplotypes of G. f. fuscipes recovered from the analysis of 485 individuals in the Lake Victoria Basin, Uganda. Blue dots represent sampled localities, pie charts indicate frequencies of the haplotypes in the sampled localities and each colour in a pie chart represents a haplotype. The inset in the upper right corner shows the location of sampling sites with reference to the whole of Uganda and neighboring countries



We obtained 429 Gff tsetse fly samples from 14 localities in the entire Ugandan portion of the Lake Victoria coast, covering 40,000 km2 (Fig. 1). The samples were collected from continental (Masaka, MA; Entebbe, EB; Budondo, BD; Okame, OK and Busime, BU) and offshore islands on Lake Victoria (Buvuma Islands: Buvuma, BV; Bugaya, BY; Buziri, BZ and Lingira, LI; Koome islands: Damba, DB; Nsazi, NS and Koome, KO; Ssese islands: Kalangala, KG and Ssese, SS) as shown in Table 1a. A maximum of 30 individuals per location were collected between October 2009 and March 2011, using biconical traps [35], following standard protocols. Whole tsetse samples were stored individually in 90 % ethanol and kept at 4 °C.

Table 1 Sampling localities and genetic diversity statistics for the mitochondrial COII sequences from 18 Gff localities in the lake Victoria basin in Uganda. N = Number of individuals analyzed, Nh = Number of haplotypes, Hd = haplotype diversity and π = nucleotide diversity. (1a) New sampling localities for this study. (1b) Sampling localities added from previous studies (Echodu et al. 2013)

DNA extraction, Amplification and Sequencing

Total genomic DNA was extracted from legs of individual tsetse flies using the PrepGEM™ Insect kit (ZYGEM Corp. Ltd) as per the manufacturer’s protocol. A 570 bp fragment of mtDNA COII gene was PCR-amplified using the primers COIF1 (5’ – CCT CAA CAC TTT TTA GGT TTA G – 3’) and COIIR1 (5’ – GGT TCT CTA ATT TCA TCA AGT A – 3’), as described by [29]. Reactions contained 1–10 ng of template DNA, 2.6 μl (5X) buffer (GoTaq colorless, Promega), 1.1 μl (10 mM) dNTPs, 0.5 μl (10 mM) primers, 1.1 μl (25 mM) MgCl2, and 0.1 μl (U/μL) GoTaq polymerase, and 6.9 μl of water for a total volume of 13 μl. Amplification involved an initial denaturation step at 95 °C for 5 min, followed by 95 °C for 30 s of denaturation, 40 cycles each for 30 s at 50 °C for annealing, 45 s at 72 °C for extension and a final extension step at 72 °C for 20 min. The PCR products were purified using ExoSAP-IT (Affymetrix, Inc.) as per the manufacturer’s protocol. Sequencing was carried out for both forward and reverse strands at the DNA Analysis Facility on Science Hill at Yale University (

Chromatograms were visually inspected and sequences trimmed to remove poor quality data using the CLC Workbench (CLC Bio Denmark). The forward and reverse strands were used to create a consensus sequence for each sample. In addition to the newly sequenced 429 samples, mtDNA COII gene sequences for 56 Gff individuals from 4 additional sampling sites in the basin [36] (Table 1b) were added to the dataset. Thus, making the final number of analysed sequences from the same sampling sites where previous microsatellite data were collected [16] 485. The total length of these sequences, prior to analysis, was 570 bp. This fragment was trimmed to a 526 bp long fragment common to all the samples.

Genetic diversity, network and population structure analysis

We analyzed the data for haplotype diversity (Hd) and nucleotide diversity (π) using DnaSP version 5.10 [37]. Significance was assessed with 1000 permutations. The partitioning of the genetic diversity within and among sampling sites was evaluated using the analysis of molecular variance (AMOVA) as implemented in Arlequin 3.5 [38]. We performed a nested analysis of variance (AMOVA) framework to partition the total amount of genetic differentiation between hierarchical levels of population subdivision [39], and produced Φ-statistics that measure the similarity of pairs of haplotypes in each hierarchical level of the analysis, relative to pairs drawn from the pool of sequences in the higher hierarchical level. Significance of the Φ-statistics was tested by permuting haplotypes among the corresponding hierarchical levels, and recalculating the statistics to obtain their null distributions [40].

To understand the evolutionary relationship of the mtDNA haplotypes, we constructed a median-joining haplotype network [41], where individual sequences were collapsed into haplotypes using the default settings in the NETWORK 4.6.1 software ( This program implements the median-joining method in the absence of recombination. The method, which provides an estimation of the haplotype genealogical relationships, is a more powerful method than bifurcating trees, when studying phylogenetic relationships at the intraspecific level, because it allows for the inclusion of multi-furcations and reticulations [42]. The program GenGIS [43] was used to visualize haplotype diversity and its relationship between geographical localities.

Genetic differentiation among the 18 sampling sites was evaluated with and without spatial information as a priori [44], using the Bayesian approach implemented in BAPS 6 [45]. We employed the spatial model option in BAPS, using local populations inhabiting discrete habitat patches (localities) with known geographical coordinates as the population units to be clustered. All molecular data collected from a particular local population were used to obtain the posterior distribution of haplotype frequencies for that population. Under the spatial model, the genetic structure is calculated assuming a priori that the structure within a particular area depends on the neighbouring areas. This program uses a statistical genetic model that treats nucleotide frequencies and K (the number of genetically diverged groups in a population) as random variables. The best K was determined using posterior probabilities. The best partition was visualized using a Voronoi tessellation as implemented in BAPS.

To obtain pairwise estimates of genetic differentiation we computed Φ ST values among sampling sites using Arlequin 3.5 with 1000 random permutations. We used Φ ST because it also accounts for the evolutionary relatedness of the mtDNA haplotypes. To test the correlation between these pairwise genetic distances and pairwise geographic distances, we used Mantel’s test [46] with 9,999 permutations, as implemented in GenAlEx 6.5. Pairwise geographic (Euclidean) distances were generated using the coordinates of the sampling localities in GenAlEx 6.5 [47].

Demographic history

We used mismatch distributions (number of pairwise mutational differences) [48] to determine if the mtDNA data showed signatures of population expansion and calculated the raggedness statistic to analyse the goodness of fit of the population expansion model to evaluate the extent to which the distribution followed the smooth unimodal curve, which one would expect under a population growth scenario. However, as this approach does not use all the information in the sequence data, we also used Tajima’s D [49] and Fu’s F S [50] statistics to test for deviations from neutral expectations. Positive values indicate an excess of intermediate-frequency haplotypes, which might result from balancing selection or bottlenecks, while negative values reflect an excess of rare polymorphisms, which might result from population growth but also genetic hitchhiking, selective sweeps, or background selection. For all these tests we used DnaSP version 5.10 [37] and significance was evaluated by comparing observed and expected statistics to a distribution of values generated with 5000 coalescent simulations.


Genetic diversity

A 526 bp fragment of the mtDNA cytochrome oxidase II (COII) was analysed from 485 individuals from 18 localities around the Lake Victoria basin (Table 1, Fig. 1). The collection of sequences was comprised of 23 haplotypes and 29 polymorphic sites (Table 2). The number of haplotypes within each sampling site varied considerably (from 1 to 10 haplotypes per sampling location) despite equal sample sizes. Similarly, both haplotype diversities ranged widely from 0 in BU to 0.774 in BY (Table 1). On the contrary, nucleotide diversity was very low ranging from 0 in BU a coastal site and SS an island site, to 0.008 in KG in Ssese islands. These low levels of nucleotide diversity may be due to relatively recent reduction in population size or recent colonization events, as sampling effort was the same for every site (Table 1). However, the fact that for some sites we recovered high haplotypic diversity suggests differences in demographic dynamics among sites.

Table 2 Haplotype distributions among the 18 G. f. fuscipes studied, based on mitochondrial CO II sequence data: 1st column: Haplotype code name (Hap1-Hap23); 2nd column: segregating sites in each haplotype, numbers on top of 2nd column are the variable sites in the reference sequence JFJR01006635.1, dots represent identical nucleotides to the ones for Hap1. The location code names (column 3 to 20) are those shown in Table 1. The last column shows the frequency of each haplotype in the whole mitochondrial CO II sequence data

Table 3 shows the results for the AMOVA analysis on the 18 sampling sites; overall genetic variation within sampling sites was much larger (85.21 %) than the variation among sampling sites (14.79 %), which is indicative of shallow levels of genetic divergence among sampling sites. This is further supported by the distribution of haplotypes among the sampled localities (Fig. 1, Tables 1 and 2) and their evolutionary relationships (Fig. 2).

Table 3 Results of AMOVA (Excoffier et al. 1992) on 485 mitochondrial COII sequences from 18 localities in the Lake Victoria Basin, Uganda, computed using the Arlequin program (Excoffier et al. 2009). Significance was tested using 1000 random permutations
Fig. 2

Median-Joining network [41] for 23 COII mtDNA haplotypes of G. f. fuscipes from 485 individuals in the Lake Victoria Basin, Uganda. Each colour represents a haplotype and the size of the circle is proportional to the number of individuals with that haplotype. Each line represents one mutational step, colour coding is the same as that in Fig. 1 and a white circle represents an inferred missing haplotype

Figure 1 shows that Haplotype 1 (HAP1), the most common haplotype (72.4 %; Table 2), is ubiquitous. The second most common haplotype, HAP2 (Table 2) was by far much less frequent (3.9 %) than HAP1 and was found in only 8 localities. Six other haplotypes occurred in two or more localities. These eight haplotypes represented 95.7 % of the sample. The other fifteen haplotypes (4.3 % of the sample) were unique to specific localities. The high percentage of shared haplotypes with the most common haplotype found at all sampling sites suggests high connectivity of Gff in the past. However, some haplotypes were retrieved from only geographically proximate areas, suggesting the occurrence of some genetic structuring. For example, HAP4 (Table 2; Light-blue in Fig. 1) was retrieved from BD, BV, LI, BZ, MG and BY, all geographically proximate localities; HAP7 (Green in Fig. 1) appears only in the extreme west of the basin (SS and MA), and HAP8 was retrieved exclusively from OK, a sampling site at the eastern edge of the Gff belt (Fig. 1; Table 2). Interestingly, HAP5 occurred exclusively on islands, particularly sites KG, NS, KO, DB and BY, some of which are located more than 100 km apart.

Figure 2 shows the evolutionary relationships among the 23 haplotypes. The network shows two haplogroups separated by five mutational steps. The most common haplotype (HAP1) is located internally in the larger haplogroup, with the other haplotypes arising from it, suggesting that HAP1 is the ancestral haplotype of this haplogroup. In addition, a star-like polytomy separated from HAP1 by two mutation steps was found in this haplogroup. The second haplogroup has only two haplotypes, HAP5 and HAP18, each separated by one mutation step from an unknown haplotype. Overall the network shows very low levels of sequence divergence among haplotypes and a high frequency of singletons (i.e., haplotypes seen only once in a group of samples), a pattern suggesting recent divergence and possibly population expansion.

Demographic history

To investigate demographic history and explore evidence of recent population expansions or reductions, we carried out mismatch distribution analyses by combining all the sampling sites (Additional file 1: Figure S1). Harpending’s Raggedness index rejected the null hypothesis of exponential growth (r > 0.05, P > 1.000). The observed distributions suggest a unimodal pattern, indicating a signal of past population expansion. Tajima’s D and Fu’s F S (Table 4) were both negative and significant for the study area (D = −1.661; P = 0.014; F S  = −10.787, P = 0.009), confirming population expansion of the Gff population in this part of the basin. At locality level, however, F S and D statistics (Table 4) confirmed demographic dynamics being different among localities as the values were negative for some sites and positive for the others.

Table 4 Neutrality and Demographic parameters: Tajima’s D, Fu’s Fs, Harpending’s raggedness index (r) based on mitochondrial COII sequence data of 18 localities of G. f. fuscipes belonging to the Lake Victoria Basin as implemented in the program DnaSP (Librado and Rozas 2009) for population size changes. In bold are statistically significant values at 0.05 Significance level

Population differentiation patterns

Figure 3 shows results of the BAPS analyses. The analysis, which incorporates spatial information of sampling sites as prior information inferred existence of three (K = 3) genetic clusters. In agreement with the shallow genetic divergence and haplotypic distribution shown above, these clusters do not group entirely according to geographical location of tsetse samples. For instance, cluster 1 (red in Fig. 3), the cluster that groups the majority of individuals (67.8 %) includes tsetse flies from all sampling sites regardless of their geographic proximity. On the other hand there is some evidence of genetic structuring, because cluster 2 (blue in Fig. 3) includes only individuals from each of the Buvuma archipelago sites (LI, BY, BV, BZ) as well as samples from BD, a mainland site about 50 km away from the Buvuma islands. However, cluster 3 (green in Fig. 3) includes individuals from two different island groups (KG and KO) on the west side of the study area and OK, located at the opposite end of the Gff distribution in the Lake Victoria basin.

Fig. 3

Genetic clustering of local populations in the Lake Victoria basin inferred with the program BAPS [44] using mtDNA COII marker. Locality codes are those described in Table 1 (a) Mixture clustering graphical output for K = 3, where K is the optimal number of clusters identified. Each vertical block is a sampling site, colour indicates membership of its individuals to population clusters (red - cluster 1, blue – cluster 2, green – cluster 3). Localities are ordered geographically from west to east across the basin. b Spatial clustering model for K = 3, each bordered cell represents a sampling site and colour indicates membership of its individuals to the same three population clusters as in A. X and Y-axes are spatial coordinates of the localities

Similar conclusions in terms of overall levels of genetic divergence can be inferred from the pairwise Φ ST values (Table 5). Among localities these values ranged from zero between populations from Damba Island (DB) and Nsazi island (NS), located about 4 km apart in the Koome archipelago, to relatively high and statistically significant values between Budondo (BD) and Busime (BU; Φ ST  = 0.592, P ≤ 0.05), which are continental sites about 100 km apart. Samples from Lingira (LI), an island site in Buvuma islands and Nkumba (NA), a continental site about 40 km away, were not genetically distinct (Φ ST  = 0.05, P < 0.05), suggesting that there has been gene flow between islands and continental sampling sites. Surprisingly, samples from KO, a site only 5 km from DB and NS in the Koome archipelago, were genetically distinct from all samples including those from DB and NS which are only 5 km within the same archipelago, but similar to samples from KG, an island site more than 100 km away in Kalangala islands, suggesting possible long-range dispersal among the islands’ Gff. Additional file 2: Figure S2 shows the results of the Mantel test, which suggests no correlation between genetic and geographic distances among localities (R =0.109, P = 0.185), confirming the findings from the Φ ST and the BAPS analyses.

Table 5 Pairwise differentiation estimates of mtDNA Φst between the 18 localities arranged from West to East across the basin: Computed in Arlequin 3.5 (Excoffier et al. 2009), bold numbers show statistically significant comparisons at 0.05 Significance level


Lack of mtDNA structure in Lake Victoria basin Gff

Sequence analysis of the COII mitochondrial DNA fragment from Gff populations across the Lake Victoria basin revealed very little genetic structuring. Most of the genetic variation at this locus was found within rather than between sampling sites (Table 3). Bayesian clustering inferred three spatially overlapping clusters, which do not group according to geographical origin of the samples. The overlapping spatial clustering could be a result of stochasticity in the process of lineage sorting of haplotypes followed by introgression due to gene flow from continental sites not included in this study, resulting in spatial mixing of the haplotype groups. A previous study indeed showed high levels of gene flow among different continental sampling sites separated by hundreds of kilometers in both Southern and Northern Uganda [29], reinforcing this hypothesis. Given the data at hand, it is not possible to distinguish between ancestral polymorphisms or recent introgression, as both could produce the observed patterns [51]. On the other hand the influence of reproductively inherited symbionts such as Wolbachia [52] could be investigated.

Regardless of the very little genetic structuring that we detected among sampling sites, we found relatively high levels of genetic diversity, as 14 of the 23 haplotypes recovered in this study are singletons (Table 2). Although this could reflect technical artifacts rather than the actual diversity of this mtDNA fragment, we feel that this is unlikely for a variety of reasons. The observed mtDNA sequence diversity is unlikely to be due to the presence of transcriptionally inactive mtDNA fragments inserted in the nuclear genome, numts [53]. We did not find evidence of mixed templates when sequencing the PCR products, or stop codons when the DNA sequences were translated into amino acids. Moreover, numts were never observed in any of previous studies of Gff mtDNA polymorphism, which included samples from a larger spatial scale than the current study [29, 30, 36, 54]. It is therefore unlikely that the patterns observed in the mtDNA data could be attributed to accidental cross-contamination or sample mixing, given that we checked for cross-contamination at each step, including negative controls. Indeed data were collected for both markers at the same time from the same DNA extractions and the microsatellite markers did not show any evidence of cross-contamination [16]. Additionally, several samples were genotyped and sequenced in duplicate and yielded identical results.

Genetic drift and gene flow equilibrium in Lake Victoria basin Gff

The Mantel test (Additional file 2: Figure S2) detected no significant correlation between geographic and genetic distance, and pairwise Φ ST comparisons showed higher differentiation between geographically close localities than distant localities, which suggests the existence of a complex and locality-dependent population. This could be facilitated by local environmental conditions, which would allow both genetic drift and gene flow to occur concurrently. Gff are found in highly fragmented habitats where genetic drift could be the predominant force. However, Gff also occur in contiguous riverine habitats along Lake Victoria and the Nile River, which can facilitate gene flow by acting as a corridor for individual dispersal among localities across the basin. The role of contiguous riverine habitat in facilitating long-range dispersal in tsetse has been previously discussed for the same species in Uganda but at a larger geographic scale [36, 55] and also for another riverine tsetse species G. tachinoides in Burkina Faso [56].

The haplotype network depicts a frequent haplotype (HAP1), with the majority of the haplotypes (91.3 %) in the network originating from it. This haplotype has a range-wide distribution across the basin, and more than 95 % of all the haplotypes are shared among the localities, suggesting long-range gene flow across the basin. Despite the long-range gene flow, some haplogroups were retrieved from only geographically proximate localities, which, coupled with the presence of private haplotypes at some localities further supports the importance of both gene flow and genetic drift in shaping the observed genetic patterns. Localities around the source of the Nile, such as BD, BV, LI, BY and MG (Fig. 1) had the highest haplotype diversities, confirming the role of contiguous riverine habitats in facilitating gene flow and the importance of the river Nile in facilitating gene flow between the lake Victoria basin and the northern Gff lineage, as previously suggested [29]. Interestingly, one haplogroup was exclusively retrieved from islands, some of which are more than 100 km apart. This haplogroup is five mutation steps from the dominant haplogroup, indicating higher connectivity among the Gff populations across islands than to the coastal area, suggesting that the islands could have been connected in the past.

Localized demographic dynamics

Both mismatch distributions and neutrality tests indicated demographic expansion for the study area, but the difference in demographic dynamics exhibited by the neutrality tests at locality level is further evidence for population sub-division rather than panmixia of Gff in the basin. The positive Fs and D values indicate that Gff experienced localized population reductions or re-colonization events at some localities as opposed to the expansions at the other localities showing negative values. This could be a result of unsustained small-scale tsetse control projects that register temporally successes at those localities, but are followed by re-infestations from adjacent un-treated areas when the projects end.

Comparison between mitochondrial and nuclear DNA markers

In another study of Gff [16] in the Lake Victoria basin, frequency-based analysis of microsatellite genotypic data revealed a complex genetic structure with four distinct meta-populations, which, although genetically distinct and spatially separated, also showed considerable amounts of gene flow. The microsatellite analyses also revealed existence of isolation by distance (IBD) within and between the distinct genetic clusters. Genetically derived dispersal distances varied between clusters ranging from about 2.5 to 14 km and matched reasonably well with dispersal rates predicted from mark–release–recapture (MRR) data for Gff and other riverine species [9]. Hierarchical F ST and individual assignment tests indicated that there were four genetic clusters, and that flies in clusters 3 and 4 shared many migrants, while clusters 1 and 2 were more isolated. The difference in gene flow among these clusters was attributed to heterogeneity in human influence. Clustering of Gff from island sites with Gff from mainland sites led to a conclusion that the Lake Victoria does not act as a barrier to fly movement and gene flow, possibly due to passive dispersal mediated by boat traffic.

The results presented in this study show both agreements and disagreements with previous results [16]. Both studies recorded high gene flow between islands and adjacent mainland sites; however, they differed in the level of genetic structuring that was identified. Unlike mitochondrial DNA, microsatellite data indicated the presence of four distinct genetic clusters in a small area, with different degrees of isolation from the rest. Additionally, in contrast to mitochondrial DNA, which indicated population expansion throughout Gff demographic history, microsatellites pointed to population stability over several generations in the Lake Victoria region [16, 30] as well as other areas in Uganda [36]. Since mitochondrial DNA has lower mutation rates than microsatellites [33, 34], does not recombine [57], and has a smaller effective population size because of its maternal inheritance, it provides insights on older evolutionary events than microsatellite data [58]. So, by revealing patterns further back in the demographic history of Gff in the Lake Victoria region, the mtDNA results in this study complement inferences based on microsatellites [16]. One mtDNA haplotype was present in all sampling sites, suggesting a higher degree of connectivity between these sites in the past. It is possible that due to human activity, especially vector control efforts and human development, Gff populations have become more and more fragmented, which is why the microsatellites reveal more genetic structuring. The Gff structuring revealed by microsatellites in this study is also in line with recent work modeling predictions of Gff distributions in southern Uganda [59].


Results of gene structuring and connectivity based on partial mtDNA sequences alone may underestimate current levels of genetic differentiation. As revealed by microsatellite data, lack of significant partitioning among groups or populations based on mtDNA data may not necessarily be indicative of current panmixia, but instead reflects historical events. This study has revealed the demographic history of Gff in the Lake Victoria basin, enabling us to better understand the factors behind the observed tsetse re-emergences after successful control interventions in the basin.

In terms of tsetse and trypanosomiasis control, interventions implemented at local scales are unlikely to produce long-lasting results due to re-invasion(s) from adjacent areas and/or residual tsetse pockets. As such, the high levels of genetic mixing between Gff in the island and mainland sites suggests that island and the mainland populations should be handled at the same time when implementing interventions. These findings support the need for an integrated area-wide elimination strategy for tsetse and trypanosomiasis from Uganda.


  1. 1.

    Leak SGA. Tsetse biology and ecology: their role in the epidemiology and control of trypanosomosis. Wallingford: CABI; 1998.

  2. 2.

    Van den Bossche P, de La Rocque S, Hendrickx G, Bouyer J. A changing environment and the epidemiology of tsetse-transmitted livestock trypanosomiasis. Trends Parasitol. 2010;26:236–43.

    PubMed  Article  Google Scholar 

  3. 3.

    Simarro PP, Cecchi G, Franco JR, Paone M, Diarra A, Ruiz-Postigo JA, et al. Estimating and Mapping the Population at Risk of Sleeping Sickness. PLoS Negl Trop Dis. 2012. doi:10.1371/journal.pntd.0001859

  4. 4.

    Matovu E, Stewart ML, Geiser F, Brun R, Mäser P, Wallace LJ, Burchmore RJ, Enyaru JC, Barrett MP, Kaminsky R, Seebeck T, de Koning HP. Mechanisms of arsenical and diamidine uptake and resistance in Trypanosoma brucei. Eukaryot Cell. 2003;2:1003–8.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  5. 5.

    Cecchi G, Paone M, Franco JR, Fèvre EM, Diarra A, Ruiz JA, et al. Towards the Atlas of human African trypanosomiasis. Int J Health Geogr. 2009;8:15.

    PubMed Central  PubMed  Article  Google Scholar 

  6. 6.

    Alsan BM, Alesina A, Bates R, et al. The Effect of the TseTse Fly on African Development. American Economic Review. 2015;105(1):382–410.

  7. 7.

    Cecchi G, Mattioli RC, Slingenbergh J, La Rocque SD, Feldmann U. Standardizing land cover mapping for tsetse and trypanosomiasis decision making. PAAT Tech Sci Ser. 2008;8:1–97.

    Google Scholar 

  8. 8.

    Molyneux D, Hallaj Z, Keusch GT, McManus DP, Ngowi H, Cleaveland S, Ramos-Jimenez P, Gotuzzo E, Kar K, Sanchez A, Garba A, Carabin H, Bassili A, Chaignat CL, Meslin FX, Abushama HM, Willingham AL, Kioy D. Zoonoses and marginalised infectious diseases of poverty: where do we stand? Parasit Vectors. 2011;4:106.

    PubMed Central  PubMed  Article  Google Scholar 

  9. 9.

    Rogers D. Study of a Natural Population of Glossina fuscipes fuscipes Newstead and a Model of Fly Movement. J Anim Ecol. 1977;46:309.

    Article  Google Scholar 

  10. 10.

    Jordan AM, Curtis CF. Productivity of Glossina morsitans morsitans Westwood maintained in the laboratory, with particular reference to the sterile-insect release method. Bull World Health Organ. 1972;46:33–8.

    CAS  PubMed Central  PubMed  Google Scholar 

  11. 11.

    Green CH. Advances in Parasitology Volume 34. Adv Parasitol. 1994;34:229–91.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Thomson PC, Marlow NJ, Rose K, Kok NE. The effectiveness of a large-scale baiting campaign and an evaluation of a buffer zone strategy for fox control. Wildl Res. 2000;27:465.

    Article  Google Scholar 

  13. 13.

    Vreysen MJB, Saleh KM, Ali MY, Abdulla AM, Zhu ZR, Juma KG, Dyck VA, Msangi AR, Mkonyi PA, Feldmann HU. Glossina austeni (Diptera: Glossinidae) Eradicated on the Island of Unguja, Zanzibar, Using the Sterile Insect Technique. J Econ Entomol. 2000;93:123–35.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Mamoudou A, Zoli A, Delespaux V, Cuisance D, Geerts S, van den Bossche P. Half a century of tsetse and animal trypanosomosis control on the Adamawa plateau in Cameroon. Rev Elev Med Vet Pays Trop. 2009;62:33–8.

    Google Scholar 

  15. 15.

    Gooding RH, Krafsur ES. Tsetse genetics: contributions to biology, systematics, and control of tsetse flies. Annu Rev Entomol. 2005;50:101–23.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  16. 16.

    Hyseni C, Kato AB, Okedi LM, Masembe C, Ouma JO, Aksoy S, et al. The population structure of Glossina fuscipes fuscipes in the Lake Victoria basin in Uganda: implications for vector control. Parasit Vectors. 2012;5:222.

    PubMed Central  PubMed  Article  Google Scholar 

  17. 17.

    Solano P, Kaba D, Ravel S, Dyer NA, Sall B, Vreysen MJ, Seck MT, Darbyshir H, Gardes L, Donnelly MJ, De Meeûs T, Bouyer J. Population genetics as a tool to select tsetse control strategies: Suppression or eradication of Glossina palpalis gambiensis in the niayes of senegal. PLoS Negl Trop Dis. 2010;4:1–11.

    Article  Google Scholar 

  18. 18.

    Aksoy S, Caccone A, Galvani AP, Okedi LM. Glossina fuscipes populations provide insights for human African trypanosomiasis transmission in Uganda. Trends Parasitol. 2013;29:394–406.

    PubMed Central  PubMed  Article  Google Scholar 

  19. 19.

    Melachio TTT, Simo G, Ravel S, De Meeûs T, Causse S, Solano P, et al. Population genetics of Glossina palpalis palpalis from central African sleeping sickness foci. Parasit Vectors. 2011;4:140.

    PubMed Central  PubMed  Article  Google Scholar 

  20. 20.

    Kagbadouno MS, Camara M, Bouyer J, Courtin F, Onikoyamou MF, Schofield CJ, et al. Progress towards the eradication of Tsetse from the Loos islands, Guinea. Parasit Vectors. 2011;4:18.

    PubMed Central  PubMed  Article  Google Scholar 

  21. 21.

    Koné N, Bouyer J, Ravel S, Vreysen MJB, Domagni KT, Causse S, et al. Contrasting population structures of two vectors of African Trypanosomoses in Burkina Faso: Consequences for control. PLoS Negl Trop Dis. 2011;5:1–10.

    Article  Google Scholar 

  22. 22.

    Ouma JO, Marquez JG, Krafsur ES. Patterns of genetic diversity and differentiation in the tsetse fly Glossina morsitans morsitans Westwood populations in East and southern Africa. Genetica. 2007;130:139–51.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  23. 23.

    Krafsur ES, Marquez JG, Ouma JO. Structure of some East African Glossina fuscipes fuscipes populations. Med Vet Entomol. 2008;22:222–7.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  24. 24.

    Waiswa C, Picozzi K, Katunguka-Rwakishaya E, Olaho-Mukani W, Musoke RA, Welburn SC. Glossina fuscipes fuscipes in the trypanosomiasis endemic areas of south eastern Uganda: apparent density, trypanosome infection rates and host feeding preferences. Acta Trop. 2006;99:23–9.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Hao Z, Kasumba I, Lehane MJ, Gibson WC, Kwon J, Aksoy S. Tsetse immune responses and trypanosome transmission: implications for the development of tsetse-based strategies to reduce trypanosomiasis. Proc Natl Acad Sci U S A. 2001;98:12648–53.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  26. 26.

    Hamilton PB, Gibson WC, Stevens JR. Patterns of co-evolution between trypanosomes and their hosts deduced from ribosomal RNA and protein-coding gene phylogenies. Mol Phylogenet Evol. 2007;44:15–25.

    CAS  PubMed  Article  Google Scholar 

  27. 27.

    Luyimbazi F: Detailed work plan/action plan for the collection of entomological baseline data. Integrated area-wide program for the creation of sustainable tsetse and trypanosomiasis free areas in the Lake Victoria basin. Entebbe: Ministry of Agriculture, Animal Industry and Fisheries; 2006.

  28. 28.

    De La Rocque S, Augusseau X, Guillobez S, Michel V, De Wispelaere G, Bauer B, et al. The changing distribution of two riverine tsetse flies over 15 years in an increasingly cultivated area of Burkina Faso. Bull Entomol Res. 2001;91:157–66.

    Google Scholar 

  29. 29.

    Beadell JS, Hyseni C, Abila PP, Azabo R, Enyaru JCK, Ouma JO, Mohammed YO, Okedi LM, Aksoy S, Caccone A (2010) Phylogeography and population structure of Glossina fuscipes fuscipes in Uganda: Implications for control of tsetse. PLoS Negl Trop Dis. doi:10.1371/journal.pntd.0000636

  30. 30.

    Echodu R, Beadell JS, Okedi LM, Hyseni C, Aksoy S, Caccone A. Temporal stability of Glossina fuscipes fuscipes populations in Uganda. Parasit Vectors. 2011;4:19.

    PubMed Central  PubMed  Article  Google Scholar 

  31. 31.

    Gillham NW. Organelle genes and genomes. USA: Oxford University Press; 1994.

    Google Scholar 

  32. 32.

    Rokas A, Williams BL, King N, Carroll SB. Genome-scale approaches to resolving incongruence in molecular phylogenies. Nature. 2003;425:798–804.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Whittaker JC, Harbord RM, Boxall N, Mackay I, Dawson G, Sibly RM. Likelihood-based estimation of microsatellite mutation rates. Genetics. 2003;164:781–7.

    PubMed Central  PubMed  Google Scholar 

  34. 34.

    Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, Brandon M, Easley K, Chen E, Brown MD, Sukernik RI, Olckers A, Wallace DC. Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci U S A. 2003;100:171–6.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  35. 35.

    Challier A, Laveissière C. Un nouveau piège pour la capture des glossines (Glossina: Diptera, Muscidae): description et essais sur le terrain. Cah ORSTOMSérie Entomol Médicale Parasitol. 1973;11:251–62.

    Google Scholar 

  36. 36.

    Echodu R, Sistrom M, Hyseni C, Enyaru J, Okedi L, Aksoy S, Caccone A. Genetically distinct Glossina fuscipes fuscipes populations in the lake Kyoga region of Uganda and its relevance for human African trypanosomiasis. Biomed Res Int. 2013. doi:10.1155/2013/614721

  37. 37.

    Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19:2496–7.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Excoffier L, Lischer H. Arlequin 3.5: An Integrated Software Package for Population Genetics Data Analysis. 2011. doi:10.1111/j.1755-0998.2010.02847.x

  39. 39.

    Weiss KM. Genetic data analysis: Method for discrete population genetic data. By B. S. Weir. xii + 337 pp. Sunderland, MA: Sinauer associates, 1990, $27.00 (paper), $48.00 (cloth). Am J Hum Biol. 1991;3:212–3.

    Article  Google Scholar 

  40. 40.

    Excoffier L, Smouse P, Quattro J. Analysis of molecular variance infered from metric distances among DNA haplotypes: application to human mitochondrial DNA restricyion data. Genetics. 1992;131:479–91.

    CAS  PubMed Central  PubMed  Google Scholar 

  41. 41.

    Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    David Posadaand KeithA. Crandall (2001) Intraspecific gene genealogies: trees grafting into networks. In: TRENDS Ecol. Evol. Vol.16 No.1 January 2001.

  43. 43.

    Parks DH, Porter M, Churcher S, Wang S, Blouin C, Whalley J, et al. GenGIS: A geospatial information system for genomic data. Genome Res. 2009;19:1896–904.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  44. 44.

    Corander J, Sirén J, Arjas E. Bayesian spatial modeling of genetic population structure. Comput Stat. 2008;23:111–29.

    Article  Google Scholar 

  45. 45.

    Corander J, Corander J, Marttinen P, Marttinen P, Tang J, Tang J. BAPS: Bayesian Analysis of Population Structure. Analysis. 2007;1–28.

  46. 46.

    Smouse PE, Long JC, Sokal RR. Multiple Regression and Correlation Extensions of the Mantel Test of Matrix Correspondence. Syst Zool. 1986;35:627.

    Article  Google Scholar 

  47. 47.

    Peakall R, Smouse PE. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.

    Article  Google Scholar 

  48. 48.

    Harpending HC. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Hum. Biol. an Int. Rec. Res. Hum Biol. 1994;66-591-600.

  49. 49.

    Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

    CAS  PubMed Central  PubMed  Google Scholar 

  50. 50.

    Fu YX. Statistical Tests of Neutrality of Mutations Against Population Growth, Hitchhiking and Background Selection. Genetics. 1997;147:915–25.

    CAS  PubMed Central  PubMed  Google Scholar 

  51. 51.

    Kvie KS, Hogner S, Aarvik L, Lifjeld JT, Johnsen A. Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata). Ecol Evol. 2013;3:126–44.

    PubMed Central  Article  Google Scholar 

  52. 52.

    Hurst GDD, Jiggins FM. Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts. Proc Biol Sci. 2005;272:1525–34.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  53. 53.

    Richly E, Leister D. NUMTs in sequenced eukaryotic genomes. Mol Biol Evol. 2004;21:1081–4.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Abila PP, Slotman MA, Parmakelis A, Dion KB, Robinson AS, Muwanika VB, Enyaru JC, Okedi LM, Aksoy S, Caccone A. High levels of genetic differentiation between Ugandan Glossina fuscipes fuscipes populations separated by Lake Kyoga. PLoS Negl Trop Dis. 2008;2:e242.

    PubMed Central  PubMed  Article  Google Scholar 

  55. 55.

    Beadell JS, Hyseni C, Abila PP, Azabo R, Enyaru JCK, Ouma JO, et al. Phylogeography and population structure of Glossina fuscipes fuscipes in Uganda: implications for control of tsetse. PLoS Negl Trop Dis. 2010;4:e636.

    PubMed Central  PubMed  Article  Google Scholar 

  56. 56.

    Bouyer J, Balenghien T, Ravel S, Vial L, SidibÉ I, ThÉvenon S, Solano P, De MeeÛs T. Population sizes and dispersal pattern of tsetse flies: Rolling on the river? Mol Ecol. 2009. doi:10.1111/j.1365-294X.2009.04233.x

  57. 57.

    Buburuzan L, Gorgan L, Bara I. Types of Dna Used in Speciation and Phylogeny Studies. Analele Ştiinţifice ale Universităţii "Alexandru Ioan Cuza”, Secţiunea Genetică şi Biologie Moleculară, TOM VIII. 2007.

  58. 58.

    Dyer RJ, Nason JD, Garrick RC. Landscape modelling of gene flow: improved power using conditional genetic distance derived from the topology of population networks. Mol Ecol. 2010;19(17):3746–59.

    PubMed  Article  Google Scholar 

  59. 59.

    Albert M, Wardrop NA, Atkinson PM, Torr SJ, Welburn SC. Tsetse Fly (G.f. fuscipes) Distribution in the Lake Victoria Basin of Uganda. (2015). PLoS Negl Trop Dis. 2015; 9(4):e0003705. doi:10.1371/journal.pntd.0003705.

Download references


This study was supported by grants from NIH (R01 AI068932 and D43 TW007391) to SA, AC and LMO, and WHO-TDR (A80132) to JOO, AC and LMO. The research was accomplished while ABK was a Fogarty Research Fellow at Yale University. We are thankful for the support of Drs. Vincent Muwanika and Anne Akol (Makerere University, Uganda), PATTEC Uganda’s STATFA Project and Mukono District Administration. We are grateful to the technical staff of NaLIRRI for excellent assistance with field sampling.

Author information



Corresponding author

Correspondence to Charles Masembe.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AC, SA, LMO, CM and JOO designed the study. ABK and LMO collected samples in Uganda. CH and ABK performed the lab work. ABK provided background and ecological information, carried out the statistical analyses and wrote the initial draft of the manuscript. AC, CH, SA, and CM revised the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1: Figure S1.

Mismatch distributions plot [48] obtained using pairwise differences in mitochondrial COII sequence nucleotides for Glossina f. fuscipes in the lake Victoria Basin, Uganda. On the X-axis are the pairwise nucleotide differences, Y-axis are the number of pairs (Frequency). The solid grey lines show observed frequency distribution while the dotted black lines show the distribution expected under constant growth. The data were obtained using DNASP version 5.10 [37].

Additional file 2: Figure S2.

Mantel Test plot of genetic distance (Φst /(1- Φst)) versus geographic distance for pairwise comparisons among 18 localities of G. f. fuscipes in the lake Victoria Basin, Uganda. Blue dots represent pairwise comparisons of localities and the black line is the linear correlation of genetic and geographic distances across the basin. There is no isolation by distance (R = 0.109, P value = 0.185).

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kato, A.B., Hyseni, C., Okedi, L.M. et al. Mitochondrial DNA sequence divergence and diversity of Glossina fuscipes fuscipes in the Lake Victoria basin of Uganda: implications for control. Parasites Vectors 8, 385 (2015).

Download citation


  • Glossina fuscipes fuscipes
  • Tsetse fly
  • Trypanosomiasis
  • Vector control
  • Gene flow
  • Demographic dynamics