Skip to main content

Evidence of temporal stability in allelic and mitochondrial haplotype diversity in populations of Glossina fuscipes fuscipes (Diptera: Glossinidae) in northern Uganda



Glossina fuscipes fuscipes is a tsetse species of high economic importance in Uganda where it is responsible for transmitting animal African trypanosomiasis (AAT) and both the chronic and acute forms of human African trypanosomiasis (HAT). We used genotype data from 17 microsatellites and a mitochondrial DNA marker to assess temporal changes in gene frequency for samples collected between the periods ranging from 2008 to 2014 in nine localities spanning regions known to harbor the two forms of HAT in northern Uganda.


Our findings suggest that the majority of the studied populations in both HAT foci are genetically stable across the time span sampled. Pairwise estimates of differentiation using standardized FST and Jost’s DEST between time points sampled for each site were generally low and ranged between 0.0019 and 0.1312 for both sets of indices. We observed the highest values of FST and DEST between time points sampled from Kitgum (KT), Karuma (KR), Moyo (MY) and Pader (PD), and the possible reasons for this are discussed. Effective population size (Ne) estimates using Waple’s temporal method ranged from 103 (95 % CI: 73–138) in Kitgum to 962 (95 % CI: 669–1309) in Oculoi (OC). Additionally, evidence of a bottleneck event was detected in only one population at one time point sampled; Aminakwach (AM-27) from December 2014 (P < 0.03889).


Findings suggest general temporal stability of tsetse vectors in foci of both forms of HAT in northern Uganda. Genetic stability and the moderate effective population sizes imply that a re-emergence of vectors from local residual populations missed by control efforts is an important risk. This underscores the need for more sensitive sampling and monitoring to detect residual populations following control activities.


Tsetse flies (Diptera: Glossinidae) are vectors of human African trypanosomiasis (HAT) and animal African trypanosomiasis (AAT), two diseases that exert a significant constraint on human health, animal production and agricultural livelihood in vast areas of rural sub-Saharan Africa. The occurrence of the two diseases follows the restricted distribution of the tsetse flies across 37 countries, covering more than nine million km2 between 14° North and 20° South [1, 2]. In 2000, the African Union recognized trypanosomiasis as “one of Africa’s greatest constraints to socioeconomic development” [3].

HAT is among the most debilitating diseases on the continent, with an estimated 70 million people at risk [4]. The HAT disease presents in two different forms; the Rhodesian form that is restricted to eastern and southern sub-Saharan Africa and is caused by Trypanosoma brucei rhodesiense, and the Gambian form that is restricted to central and western Africa and is caused by T. b. gambiense. The two trypanosome subspecies are identical morphologically but display stark differences in epidemiological features, and require different diagnosis and treatment methods [5, 6]. AAT, on the other hand, stands as a major obstacle to the development of more efficient and sustainable livestock production systems in tsetse-infested areas [7]. The Programme on African Animal Trypanosomiasis (PAAT) estimates that AAT causes approximately 3 million cattle deaths per year with farmers using approximately 35 million doses of costly trypanocidal drugs [8]. Overall, annual losses in agriculture alone have been estimated at over US $5 billion [9].

In Uganda, about two thirds of the total land area is infested with tsetse flies [10], and this is the only country known to sustain active transmission of both forms of human pathogenic trypanosomes; T. b. gambiense in the northwest and T. b. rhodesiense in the southeast [11]. There is a growing body of evidence that the disease ranges are expanding [12], and thereby narrowing a disease-free belt of less than 150 km just north of Lake Kyoga [13, 14]. Overlap of the two diseases in north-central districts of Uganda [15] will complicate diagnosis and treatment activities, and provide new challenges, as recombination between the two trypanosome forms can occur in tsetse flies and could lead to novel pathologies [16, 17]. This highlights the need for increased understanding and implementation of vector control and disease prevention measures in this region.

The primary vector of both human and animal forms of trypanosomiasis in Uganda is the tsetse species Glossina fuscipes fuscipes, where it is responsible for an estimated 90 % of all disease cases. Glossina f. fuscipes belongs to the Palpalis group of tsetse, which inhabits low bushes or forests at the edges of riverine habitats, and like other tsetse species, has a unique reproductive strategy where larvae develop in utero. Populations appear to respond to seasonal weather patterns, often disappearing during the bi-annual dry season from sites where they were previously abundant [18], which can cause seasonal bottlenecks and thereby result in large temporal changes in gene frequencies due to genetic drift. However, it has also been hypothesized that larval development in utero can stabilize population sizes during dry periods [19], which would result in stable gene frequencies due to low rates of genetic drift. Thus, our understanding of population fluctuations remains incomplete. Likewise, knowledge of tsetse dispersal is also inadequate. Ecological studies using capture-release-recapture data show that tsetse species have a high capacity for dispersal, but genetic data indicate surprisingly high differentiation among populations [20]. Studying temporal patterns of genetic variation can provide insight into some of these knowledge gaps and the apparent contradictions between ecological and genetic data [21]. Temporal population genetics data can also provide estimates of effective population sizes [22], and probability of population bottlenecks, and thereby aid in the choice of spatial and temporal parameters in vector control programmes [23].

Previous work on G. f. fuscipes temporal population dynamics has focused on southern, central and southeastern Uganda and has demonstrated stability, with little evidence of seasonal variation in population size [24, 25]. However, regions north of Lake Kyoga have not yet been investigated, and remain a high priority because tsetse from this region are genetically distinct [2628], and it is where the two forms of HAT will likely merge in the near future [13, 14]. This region also differs significantly in climate, especially in annual precipitation [29]. Southern Uganda has a somewhat cooler climate and is less humid, with mean annual rainfall near Lake Victoria often exceeding 2100–3000 mm; the high temperature varies by 2–3 °C over the year, with a mean daily high being around 26 °C. In the north, the rainfall is between 1000 and 2000 mm, and temperature varies by 5 °C over the year, with the mean daily high being 29 °C [30]. The unique genetic and climatic background of northern Uganda may therefore create different tsetse population dynamics, and thus change the best strategy of vector control for the region.

In this study, we investigate temporal changes in nuclear allele and mitochondrial haplotype frequencies of G. f. fuscipes populations in areas north of Lake Kyoga by examining temporal samples collected across multiple locations in both the T. b. gambiense and T. b. rhodesiense disease belts. We do this by evaluating genetic stability at one mitochondrial marker and 17 microsatellite loci so as to contribute reliable scientific data around which sustainable vector control programs can be planned and implemented.


Tsetse sampling and DNA extraction

Tsetse flies were sampled using biconical traps [31] in nine localities spanning two disease foci in northern Uganda during March-December 2014. Populations sampled included Arua (AR), Moyo (MY), Kitgum (KT), Pader (PD), Kole (KO), Apac (AP), Aminakwach (AM), Oculoi (OC) and Karuma (KR) (Fig. 1). We chose sites with genotype data from previous collections done between 2008 and 2011. All captured tsetse flies were preserved individually in cryotubes containing 95 % ethanol. Assuming that G. f. fuscipes goes through approximately 8 generations per year [32, 33], our temporal separation ranged from 22 to 52 generations. DNA was extracted from tsetse legs using PrepGEM Insect DNA Extraction kit (ZyGEM New Zealand, 2013), following the manufacturer’s protocols.

Fig. 1
figure 1

Map of Uganda showing sampling sites. The temporal sites are shown with population codes, as well as the two HAT disease belts in Uganda (T. b. gambiense and Tbr is T. b. rhodesiense)

Nuclear microsatellite amplification and genotyping

We evaluated patterns of nuclear DNA (nDNA) genetic diversity using 17 microsatellite loci extensively evaluated for minimal allele dropouts and null alleles and optimized for use in G. f. fuscipes in previous studies [24, 25, 27]. Amplifications were performed with fluorescently labeled forward primers (6-FAM, HEX and NED) using a touchdown PCR (10 cycles of annealing at progressively lower temperatures from 60 °C to 51 °C, followed by 35 cycles at 50 °C) in 13.0 μl reaction volumes containing 2.6 μl 5× PCR buffer, 1.1 μl 10 mM dNTPs, 1.1 μl 25 mM MgCl2 and 0.1 μl 5U/μl GoTaq (Promega, USA), 0.1 μl 100X BSA (New England Biolabs, USA), 0.5 μl 10 μM fluorescently-labeled M13 primer, 0.5 μl 10 μM reverse primer, 0.3 μl 2 μM M13-tailed forward primer. For loci C7b and GmL11, 0.5 units of Taq Gold polymerase (Life Technologies, USA) were used instead of Promega GoTaq. PCR products were multiplexed in groups of two or three and genotyped on an ABI 3730xL automated sequencer (Life technologies, USA) at the DNA Analysis Facility on Science Hill at Yale University ( Alleles were scored using the program Genemarker v 2.4.0 (Soft Genetics, State College, PA, USA) with manual editing of the automatically scored peaks.

Microsatellite marker validation and summary statistics

We evaluated microsatellite genotypes with Genepop v 4.2 [34] to test for deviation from Hardy Weinberg equilibrium (HWE), linkage disequilibrium (LD), and to estimate inbreeding coefficients (FIS) using the Markov chain method [35] with 100,000 dememorizations, 1,000 batches and 10,000 iterations per batch. Arlequin v 3.5 [36] was used to calculate observed (HO) and expected (HE) heterozygosity for each population. Significance values for all multiple tests and comparisons were adjusted using the Benjamini-Hochberg method in preference to the Bonferroni method because of lower incidence of false negatives [37]. We used the program FSTAT v 3.1 [38] to calculate allelic richness and to provide locus and population-specific estimates of microsatellite allele frequencies.

Mitochondrial DNA amplification, sequencing, and summary statistics

A 570 bp fragment of the region spanning across the mitochondrial DNA (mtDNA) COI and COII gene was PCR-amplified using the primers COII-F1 (5'- CCT CAA CAC TTT TTA GGT TTA G -3') and COI-R1 (5'- GGT TCT CTA ATT TCA TCA AGT A -3') as described by [27] with an initial denaturation step at 95 °C for 5 min, followed by 40 cycles of annealing at 50 °C, and a final extension step at 72 °C for 20 min. We used a reaction volume of 13.0 μl containing 1 μl of template genomic DNA, 2.6 μl 5× PCR buffer, 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 (Promega, USA). The PCR products were purified using ExoSAP-IT (Affymetrix Inc., USA) as per the manufacturer’s protocol. Sequencing was carried out for both forward and reverse strands on the ABI 3730xL automated sequencer at the DNA Analysis Facility on Science Hill at Yale University ( Sequence chromatograms were inspected by eye and sequences trimmed to remove poor quality data using Geneious v 6.1.8 (Biomatters, New Zealand). The forward and reverse strands were used to create a consensus sequence for each sample, and the sequences trimmed to a length of 404 bp and aligned in Geneious v 6.1.8. DnaSP v 5.10.01 [39] was used to calculate mtDNA haplotype diversity (Hd) and nucleotide diversity (π).

Temporal genetic analyses


We determined the overall proportion of the variance attributable to differences in sampling dates for all 18 sampled points using the analysis of molecular variance (AMOVA) implemented in Arlequin v 3.5; first, we ran AMOVA with two groups containing all sites at generation zero pooled together to form one group and sites with multiple generations forming another. Then, we performed AMOVA with multiple groups by having each temporal sample (generation 0, 22, 27, 29, 35, 38, 48 and 52) forming a separate group.

Genetic differentiation between temporal samples from the same population was quantified by computing pairwise FST values in Arlequin v 3.5, and Jost’s DEST [40] using DEMEtics [41] in R [42]. It is well documented that high-allelic diversity markers like microsatellites can lead to underestimates of FST [40, 43, 44]. To account for this potential bias, we standardized FST values using the formulae developed by Hedricks [43], and used Jost's DEST. P-values and confidence intervals for Jost’s DEST were also calculated based on 1000 bootstrap resamplings. To test for correlation of differentiation indices and time since first sampling, we ran a linear regression of standardized FST and DEST against number of generations separating time points sampled in JMP® v11.0 (SAS Institute Inc., Cary, NC, USA, 1989–2007).


We employed the same AMOVA strategy outlined above to determine the overall proportion of the variance attributable to differences in sampling dates for the mitochondrial data. We also quantified genetic differentiation between temporal samples from the same population by computing pairwise Φ ST values between temporal samples in Arlequin v3.5, and used the same program to perform Fisher’s exact test of sample differentiation based on haplotype frequencies with 10,000 iterations and 1000 dememorization steps.

Effective population size (Ne) estimates and tests for bottlenecks

To provide an estimate of effective populations (Ne) of each site, we employed two methods using the nuclear microsatellite markers; the modified temporal method [45] based on [46], and the linkage disequilibrium (LD) method, as implemented in NeEstimator v 2.01 [47]. The temporal method relates the standardized variance of allele frequencies across generations to the effective population size [48, 49] while the LD method uses the level of non-random associations among alleles at different loci [50, 51] to estimate the action of genetic drift, and thus Ne.

Evidence of population bottleneck events was tested using two methods, both of which are implemented in the program BOTTLENECK [52]. The first method tested for an excess of heterozygosity relative to observed allelic diversity [53], and was performed separately for each sample and microsatellite locus. Simulation of heterozygosity at mutation-drift equilibrium distributions assumed the two-phase mutation model (TPM) with 70 % single-step mutations and 30 % of multiple-step mutation, as recommended for microsatellite loci [54], and significance was assessed using Wilcoxon’s signed rank test, as recommended for fewer than 20 loci [52]. The second method tested for a bottleneck-induced mode shift in allele frequency distributions that were based on equal increments of 0.1 [55].


Summary statistics


A total of 404 G. f. fuscipes tsetse flies were genotyped using 17 microsatellite loci. Allelic richness ranged from 3.45 in population AM to 4.50 in KR and was generally similar for samples of the same locality analyzed from different time points, except for KR which differed slightly between generation 0 (4.50) and generation 35 (4.08) (Table 1). None of the pairs of loci analyzed exhibited significant LD after correction for multiple testing, which confirmed previous work that these loci are unlinked. However, locus Pg17 exhibited significant departures from Hardy-Weinberg Equilibrium after correction for multiple testing in most of the populations (results not shown), and was therefore dropped from all further analyses. HE and HO values remained almost constant between temporal samples, and showed no evidence of violating Hardy-Weinberg expectations. Site AM (for both generations 0 and 27) had the lowest observed (HO) and expected heterozygosity (HE) values, while site AR (generations 0 and 52) had the highest HO and HE values (Table 1). None of the FIS values was significantly greater than zero after correction for multiple tests; values were all positive except for one population (OC at generation 0).

Table 1 Sample sizes and genetic diversity statistics for 16 nDNA microsatellite loci and an mtDNA marker


For mitochondrial DNA sequences, we recovered a total of 26 haplotypes from the 18 sampled points. Haplotype diversity (Hd) was moderately high across samples; the number of haplotypes (h) ranged from 2 to 6. Nucleotide diversity (π) was highest in samples from KR (at generations 0 and 35), and lowest in samples from AM (at generation 27) and KT (at generation 22). However, both Hd and π remained similar for temporal samples except for KT and PD (Table 1).

Temporal genetic variation and population stability

Results of AMOVA using both microsatellite and mitochondrial DNA frequencies suggested that differences between samples from different time points do not explain a significant amount of the overall genetic variation. On the other hand, differences among sites contributed significantly to the overall variation (Table 2). STRUCTURE assignment shows that genetic structure remains homogenous from one time point to the next (Fig. 2), and that genetic distance is correlated with spatial distance rather than temporal distance.

Table 2 Results of AMOVA testing for temporal genetic structure
Fig. 2
figure 2

Structure plots of temporal samples. Plots show clustering of samples, organized north to south, with the first time-point sampled followed by the second time-point sampled


Overall, the allele frequency changes were minimal for most of the sites and appear homogenous among temporal samples from the same locations (Fig. 3). Standardized pairwise FST (F’ST) for microsatellites between time points from the same site were generally low (Table 3), with the highest values found between time points at KT (0.0354; P < 0.000); MY (0.0692; P = 0.001), PD (0.0614; P = 0.015) and KR (0.0713; P < 0.0001). The moderately high values of F’ST estimates were in concordance with DEST estimates. Number of generations separating sampled time points was not clearly correlated with differentiation indices (R 2 = 0.1416, P = 0.318; Fig. 4) suggesting that time between our samplings was sufficient to make accurate estimates. As an example, the two temporal populations from AR (AR-0 and AR-52), which were separated by 52 generations, had F’ST and DEST of 0.0385 (P = 0.016) and 0.0236 (P = 0.003) respectively, whereas other populations separated by fewer generations had similar or even higher values of differentiation indices to AR (Table 3).

Fig. 3
figure 3

Microsatellite allele frequencies from two time points in nine populations of G. f. fuscipes in northern Uganda. Numbers after location code indicate the time interval (in generations) since the first sampling. Vertical bars represent allele frequencies, and each color represents a single allele. The heading of each chart lists the microsatellite locus (A03b, B05, C7b, D101, D05, GmmL03, GmmD15, Gmm8, GpB20b, GmmA06, GmmB20, GmmL11, GpC10, CAG29, Pg28, CAG133)

Table 3 Pairwise estimates of genetic differentiation between temporal samples
Fig. 4
figure 4

Correlation of genetic differentiation with time between sampling efforts. Linear regression of standardized pairwise FST (F’ST) and DEST against time interval (generations) since the first sampling


The same trend was apparent for the mitochondrial DNA haplotype frequencies; some sampled localities gained and others lost haplotypes across time points (Fig. 5). More critically, though, the frequencies of the most common haplotypes remained fairly homogenous for pairs of samples from the same site. The Φ ST estimates between time points from the same site were generally low, with the highest values found between time points at MY (0.04715; P < 0.0001) and KR (0.01494; P < 0.0001).

Fig. 5
figure 5

Mitochondrial haplotype frequencies from two time points in nine populations of G. f. fuscipes in northern Uganda. Numbers after location codes indicate the time interval (in generations) since the first sampling. Vertical bars represent haplotype frequencies, and each color represents a haplotype. Generally, frequencies of the most common haplotypes remained similar across time points

Effective population size (Ne) estimates and tests for bottlenecks

Estimates of effective population size (Ne) using the temporal method were all bound by a 95 % confidence interval and ranged from 103 (95 % CI: 73–138) in KT to 962 (95 % CI: 669–1309) in OC. Based on the linkage disequilibrium method, only two populations (AM-27 and KR-35) exhibited 95 % confidence intervals that excluded infinity. Ne estimates using the two methods generally agreed for most populations except for KR, OC and AP. Ne estimates for both KR and OC using the LD method were much lower than from the temporal method (five and eight times respectively). In contrast, Ne estimate for AP using the LD method was higher than the temporal method, and indistinguishable from infinity. For both methods, Ne estimates were not proportional to the number of generations separating temporal samples; for example, OC had the highest Ne estimate using the temporal method at 962 and was based on samples separated by 22 generations, while MY had a Ne estimate close to the overall mean at 313 (95 % CI: 227–412) and was based on samples separated by 38 generations (Table 4).

Table 4 Effective population size estimates and tests for population bottleneck events

Following the Benjamini-Hochberg FDR correction, only one site (AM at generation 27) tested positive for bottleneck events at P < 0.0389 under the TPM model (Table 4); however, OC at generation 0 and KO at generation 29 were significant at P < 0.05 (Table 4). Additionally, tests for recent bottlenecks using the mode-shift indicator approach indicated that allele frequency distributions in all populations except AM-27 (Table 4) approximated the expected normal L-shape, thus showing no loss of rare alleles via drift as might be expected in a population which has undergone a recent bottleneck event.


We carried out an evaluation of changes in gene frequencies at microsatellite loci and at an mtDNA marker to provide insight into the temporal stability of G. f. fuscipes populations that harbor the two forms of HAT in northern Uganda. Generally, temporal population genetics studies are valuable in evaluating population stability and persistence, especially under changes in environment or following a perturbation [56]. Tsetse flies are an interesting study system because they exist in highly structured meta-populations that are sensitive to external environmental change, and because they act as vectors of dangerous human and animal diseases. The presence of the two forms of HAT in Uganda north of Lake Kyoga and the narrowing gap between the two disease belts provides added impetus for monitoring population dynamics of tsetse flies in the region. Although previous work on G.f. fuscipes from southern and central Uganda indicated population stability [24, 25], there has not been any work in the genetically unique tsetse populations found north of Lake Kyoga. Temporal stability was also demonstrated for other closely related tsetse fly species like G. pallidipes in Kenya [21]. One limitation in previous studies has been short time intervals between samples of only 1–2 years. In this study, we sampled at comparatively longer time intervals, and have demonstrated similar stability of G. f. fuscipes in northern Uganda, the region where the two forms of HAT will likely merge in the near future.

Estimates of genetic diversity in both nDNA and mtDNA indicate stable populations that are currently under migration-drift equilibrium. Microsatellite results indicate relatively stable molecular diversity indices (AR, HO and HE), while MtDNA sequences show stable h, Hd, and π except for some populations such as PD and KT (Table 2). These findings indicate stable populations that currently might be under migration-drift equilibrium. Generally low FIS estimates indicate mostly random mating within these populations. The pattern of reduced FIS values in the second sample at all sites except MY suggests a general increase in outbreeding in recent years, potentially driven by migration.

Pairwise estimates of genetic divergence and AMOVA results in both the nDNA and the mtDNA analyses indicated temporal stability of sampled populations, with low estimates of genetic differentiation between temporal samples from the same locality. Exceptions included sites KT, MY and PD, which showed moderate differentiation using both standardized FST and Jost’s DEST. The moderate differentiation observed in the temporal sites KT and PD might be partially attributable to the geographic distance between trapping locations. Trapping of the generation 22 samples from KT was done 36.7 km from the original sample. Similarly, trapping of the generation 52 samples from PD was done 51.2 km from the original sample. Tsetse populations can be genetically distinct even at geographic scales of only 4–5 km, for example, as has been observed in G. pallidipes [57], G.palpalis gambiensis [58], and even in G. f. fuscipes in Uganda [27]. Thus the differentiation we observed in KT and PD may indicate micro-geographic structure between sampling sites rather than temporal instability. On the other hand, two sites (MY and KR) represent potential regions of instability. Both F’ST and Φ ST values for these sites were highly significant (See Table 3). Temporal sampling efforts for these two sites were not geographically highly separated (5.86 km for MY and 1.22 km for KR), and therefore, these may represent the only areas of instability observed in our study.

Microsatellite-based population size estimates were variable and generally associated with large confidence intervals. Our Ne estimates were lower than most estimates from Uganda obtained by [24] but similar to some localities most proximal to our study sites, such as Masindi (MS). However, our estimates were slightly higher than estimates from southern Uganda around the Lake Victoria basin obtained by [25], and similar to estimates obtained by [59] in other riverine species like G. palpalis palpalis in West Africa. Relatively large Ne estimates coupled with evidence of low temporal differentiation indicates genetic drift does not strongly impact these populations [60]. Although Ne is generally difficult to estimate and is affected by many possible biases, the temporal method we have used here is suitable because high allelic richness (Table 1) should outweigh the upward bias in temporal method estimates of Ne [61]. Furthermore, the large sampling intervals used decreases the bias caused by overlapping generations and age structure [62]. The apparently low Ne estimated at some sites, such as KT, may be due to micro-geographic population structure rather than to a true shift in allele frequencies by genetic drift through time, and therefore may represent a false signal of low Ne.

Despite the tendency of G. f. fuscipes to exist in discrete patches and the potential for seasonal fluctuations in population size presented by the climate of northern Uganda, there is little evidence of genetic bottlenecks in the populations we studied; only population AM at generation 27, sampled in December 2014, showed evidence of a bottleneck. Our finding of limited bottleneck events is congruent with previous studies in southern Uganda [24, 25, 33]. The apparent lack of extensive seasonal variation in abundance supports the hypothesis that larval development in utero may help to relieve tsetse from harsh environmental conditions during reproduction, and hence stabilizes populations [19]. Furthermore, the networks of rivers, streams and other semi-permanent water bodies common in this region may drive overall stability. We hypothesize that waterways may be facilitating connectivity and ‘rescue effects’ between populations. Glossina f. fuscipes are known to generally disperse along waterways, following riverbeds or the edges of gallery forests, where they are able to survive low humidity conditions during dry seasons. Additionally, tsetse flies in northern Uganda have not been subjected to extensive control measures. A prolonged period of civil unrest in the region led to a disruption of control efforts and a breakdown in social infrastructure [63], and later control activities were decentralized and are currently managed by underfunded district health and veterinary authorities [64]. At the height of the insurgencies, the populace was displaced into Internally Displaced Persons’ (IDPs) camps leaving tsetse populations in natural conditions for approximately two decades. This time frame was also free of major drought [65], which potentially stabilized population sizes and distribution in the region.

In summary, we find temporal stability in gene diversities of most tsetse populations sampled, which lends further credence to the hypothesis that larval development in utero helps to stabilize populations during dry periods, and as a result, decreases seasonal variation in tsetse number. Additionally, large populations of pupa, which develop in the ground over a period of weeks, may also help to ensure the continuity of tsetse populations and would contribute to reducing the variance in genetic changes over time [24]. For northern Uganda, population stability is also probably aided by the presence of many year-round rivers and streams, and historical happenstance, which left tsetse populations in natural conditions free from control efforts for the last several decades.


The findings point to general temporal stability of tsetse vectors in both foci of the two forms of HAT in northern Uganda. Additionally, migration and gene flow is probably important in shaping the observed genetic stability, which suggests that area-wide control strategies are more desirable than localized efforts. Furthermore, genetic stability and moderate effective population sizes suggest there is a significant risk of re-emergence from local residual populations missed by control efforts. This underscores the need for more sensitive sampling techniques to detect residual populations when tsetse densities decrease due to vector control to avoid a population rebound effect when control and monitoring activities are relaxed. Our results, especially the Ne estimates, may also be useful in evaluating the effectiveness of future control efforts, for instance, by estimating reduction in Ne resulting from control.


  1. Jordan A. Control of tsetse flies (Diptera: Glossinidae) with the aid of attractants. J Am Mosq Control Assoc. 1995;11(2 Pt 2):249–55.

    CAS  PubMed  Google Scholar 

  2. FAO. Training Manual for Tsetse Control Personnel, Volume 4: Use of Attractive Devices for Tsetse Survey and Control. Information Resources and Training Manuals, Agriculture Department, Animal Production and Health Division. FAO, Rome. 1992:196. 

  3. Kabayo JP. Aiming to eliminate tsetse from Africa. Trends Parasitol. 2002;18:473–5.

    Article  PubMed  Google Scholar 

  4. Simarro PP, Cecchi G, Franco JR, Paone M, Diarra A, Ruiz-Postigo JA, Fèvre EM, Mattioli RC, Jannin JG. Estimating and mapping the population at risk of sleeping sickness. PLoS Negl Trop Dis. 2012;6:e1859.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Gibson WC, de C Marshall TF, Godfrey DG. Numerical analysis of enzyme polymorphism: new approach to the epidemiology and taxonomy of trypanomsomes of the subgenus Typanozoon. Adv Parasitol. 1980;18:175–246.

  6. Africa E. Pan African Tsetse and Trypanosomosis Eradication Campaign (PATTEC) Enhancing Africa’ S Prosperity Plan of Action Organization of African Unity. 2001.

    Google Scholar 

  7. Itard J, Cuisance D, Tacher G. Trypanosomoses: Historique Repartition Géographique. Principales Maladies Infectieuses et Parasitaires Du Bétail,” in Europe et Régions Chaudes. Volume 2. Paris: Lavoisier; 2003.

  8. PAAT. Programme Against African Trypanosomiasis, Newsletter No. 8. 2000.

    Google Scholar 

  9. Budd L. Economic analysis. In: DFID-Funded Tsetse and Trypanosomiasis Research and Development since 1980. Chatham: Dep Int Dev Livest Prod Program Anim Heal Program Resour Syst Program; 1999.

    Google Scholar 

  10. Chizyuka G. FAO Liaison Officers Summary Report. Harare: Food and Agricultural Organisation; 1998.

    Google Scholar 

  11. Hutchinson O, Fèvre E, Carrington M, Welburn SC. Lessons learned from the emergence of a new Trypanosoma brucei rhodesiense sleeping sickness ofcus in Uganda. Lancet. 2003;3:42–5.

  12. Fèvre EM, Coleman PG, Odiit M, Magona JW, Welburn SC, Woolhouse MEJ. The origins of a new Trypanosoma brucei rhodesiense sleeping sickness outbreak in eastern Uganda. Lancet. 2001;358:625–8.

  13. Simarro PP, Diarra A, Ruiz Postigo JA, Franco JR, Jannin JG. The human African trypanosomiasis control and surveillance programme of the World Health Organization 2000–2009: the Way forward. PLoS Negl Trop Dis. 2011;5:e1007.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Picozzi K, Fèvre EM, Odiit M, Carrington M, Eisler MC, Maudlin I, Welburn SC. Sleeping sickness in Uganda: a thin line between two fatal diseases. BMJ. 2005;331:1238–41.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Welburn SCOM. Recent developments in human African trypanosomiasis. Curr Opin Infect Dis. 2002;15:477–84.

    Article  CAS  PubMed  Google Scholar 

  16. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Hamilton PB, Adams ER, Malele II, Gibson WC. A novel, high-throughput technique for species identification reveals a new species of tsetse-transmitted trypanosome related to the Trypanosoma brucei subgenus, Trypanozoon. Infect Genet Evol. 2008;8:26–33.

  18. Rayaisse JB, Tirados I, Kaba D, Dewhirst SY, Logan JG, Diarrassouba A, Salou E, Omolo MO, Solano P, Lehane MJ, Pickett JA, Vale GA, Torr SJ, Esterhuizen J. Prospects for the development of odour baits to control the tsetse flies Glossina tachinoides and G. palpalis s.l. PLoS Negl Trop Dis. 2010;4:e632.

  19. Hargrove JW. Tsetse population dynamics. In: Maudlin I, Holmes PH MM, editors. The Trypanosomiases. Wallingford: CABI Publishing; 2004.

    Google Scholar 

  20. Krafsur ES, Griffiths N, Brockhouse CLBJ. Breeding structure of Glossina pallidipes (Diptera: Glossinidae) populations in East and southern Africa. Bull Entomol Res. 1997;87:67–73.

  21. Ouma JO. Population genetic studies of the tsetse fly Glossina pallidipes (Diptera: Glossinidae). PhD thesis, Iowa State University, USA. 2004.

  22. Richards CLP. Temporal changes in allele frequencies and population’s history of severe bottlenecks. Conserv Biol. 1996;10:832–9.

    Article  Google Scholar 

  23. Solano P, Kaba D, Ravel S, Dyer NA, Sall B, Vreysen MJB, 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 

  24. 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.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  26. Abila PP, Slotman MA, Parmakelis A, Dion KB, Robinson AS, Muwanika VB, Enyaru JCK, Lokedi 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.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Beadell JS, Hyseni C, Abila PP, Azabo R, Enyaru JCK, Ouma JO, Mohammed YO, Okedi LM, Aksoy S, Caccone A. Phylogeography and population structure of Glossina fuscipes fuscipes in Uganda: Implications for control of tsetse. PLoS Negl Trop Dis. 2010;4, e636.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 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;2013:614721.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Basalirwa CP. Delineation of Uganda into climatological rainfall zones using the method of principal component analysis. Int J Climatol. 1995;15:1161–77.

  31. Challier A, Laveissiere C. A new trap for capturing Glossina flies (Diptera: Muscidae), description and field trials. Cah ORSTOM Entomol Med Parasitol. 1973;11:251–62.

    Google Scholar 

  32. Hargrove JW. Extinction probabilities and times to extinction for populations of tsetse flies Glossina spp. (Diptera: Glossinidae) subjected to various control measures. Bull Entomol Res. 2005;95:13–21.

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.

    Article  PubMed  Google Scholar 

  35. TE Guo SW. Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992;48:361–72.

    Article  CAS  PubMed  Google Scholar 

  36. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.

    Article  PubMed  Google Scholar 

  37. Benjamini Y, Hochberg Y. Benjamini and Y FDR.pdf. J R Stat Soc Ser B Methodol. 1995;57(No. 1):289–300.

    Google Scholar 

  38. Goudet J, Raymond M, de Meeüs TRF. Testing differentiation in diploid populations. Genetics. 1996;144:1933–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Rojas J. DNA sequence polymorphism analysis using DnaSP. In: Bioinformatics for DNA Sequence Analysis; Methods in Molecular Biology Series. Hum Press NJ. 2009;537:337–50.

    Google Scholar 

  40. JOST L. G ST and its relatives do not measure differentiation. Mol Ecol. 2008;17:4015–26.

    Article  PubMed  Google Scholar 

  41. Gerlach G, Jueterbock A, Kraemer P, Deppermann J, Harmand P. Calculations of population differentiation based on GST and D: Forget GST but not all of statistics. Mol Ecol. 2010;19:3845–52.

    Article  PubMed  Google Scholar 

  42. R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2011.

    Google Scholar 

  43. Hendricks P. A standardized genetic differentiation measure. Evolution (N Y). 2005;59:1633–8.

    Google Scholar 

  44. Meirmans PG. The trouble with isolation by distance. Mol Ecol. 2012;21:2839–46.

    Article  PubMed  Google Scholar 

  45. Waples R. A generalized approach for estimating effective population size. Genetics. 1989;379–391.

  46. Jorde PERN. Unbiased estimator for genetic drift and effective population size. Genetics. 2007;177(2):927–35.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Do C, Waples RS, Peel D, Macbeth GM, Tillet BJOJ. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol Ecol Resour. 2013;14:209–14.

    Article  PubMed  Google Scholar 

  48. Nei M, Tajima F. Genetic drift and estimation of effective population size. Genetics. 1981;98:625–40.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. Pollak E. Population size from allele frequency changes. 1983:531–548.

  50. Hill WG. Estimation of effective population size from data on linkage disequilibrium. Genet Res (Camb). 1981;38:209–16.

    Article  Google Scholar 

  51. Waples RS. Genetic methods for estimating effective population size of Cetacean populations. In: Hoelzel AR, editor. Genetic ecology of whales and dolphins, vol. Special Issue 13. Cambridge: International Whaling Commission; 1991. p. 279–300.

  52. Piry S, Luikart G, Cornuet JM. BOTTLENECK: A computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999;90:502–3.

    Article  Google Scholar 

  53. Marie J, Luikartt G, Cornuet JM, Luikartt G. Power analysis of two tests for detecting recent population bottlenecks from allele frequency data. 1996.

  54. Di Rienzo A, Peterson CA, Garza JC, Valdes AM, Slatkin M, Freimer NB. Mutational processes of simple-sequence repeat loci in human populations. Proc Natl Acad Sci U S A. 1994;91:3166–70.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Luikart G, Allendorf FW, Sherwin WB. Distortion of allele frequency distributions bottlenecks. 1998:238–47.

  56. Ozerov MY, Veselov AE, Lumme J, Primmer CR. Temporal variation of genetic composition in Atlantic salmon populations from the Western White Sea Basin: influence of anthropogenic factors? BMC Genet. 2013;14:88.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Odulaja A, Baumgärtner J, Mihok SAI. Spatial and temporal distribution of tsetse fly trap catches at Nguruman, southwest Kenya. Bull Entomol Res. 2001;91:213–20.

    CAS  PubMed  Google Scholar 

  58. Camara M, Caro-riaño H, Ravel S, Dujardin J-p, Hervouet J-p, De MeEüs T, Kagbadouno MS, Bouyer JPS. No Title. J Med Entomol. 2006;43:853–60.

    Article  PubMed  Google Scholar 

  59. Dyer NA, Furtado A, Cano J, Ferreira F, Odete Afonso M, Ndong-Mabale N, Ndong-Asumu P, Centeno-Lima S, Benito A, Weetman D, Donnelly MJ, Pinto J. Evidence for a discrete evolutionary lineage within Equatorial Guinea suggests that the tsetse fly Glossina palpalis palpalis exists as a species complex. Mol Ecol. 2009;18:3268–82.

    Article  CAS  PubMed  Google Scholar 

  60. Palstra FP, Ruzzante DE. Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence? Mol Ecol. 2008;17:3428–47.

    Article  PubMed  Google Scholar 

  61. Turner TF, Salter LA, Gold JR. Temporal-method estimates of N e from highly polymorphic loci. Conserv Genet. 2001;2:297–308.

    Article  Google Scholar 

  62. Luikart G, Ryman N, Tallmon DA, Schwartz MK, Allendorf FW. Estimation of census and effective population sizes : the increasing usefulness of DNA-based approaches. 2010:355–373.

  63. PRDP. Peace, Recovery and Development Plan for Northern Uganda (PDRP)-2007-2010. 2007.

    Google Scholar 

  64. Menon S, Rossi R, Nshimyumukiza L, Zinszer. Revisiting zoonotic human African trypanosomiasis control in Uganda. Journal of Public Health Policy. 2016;37(1):51.

  65. Trends PR. A Climate Trend Analysis of Uganda. 2010.

    Google Scholar 

Download references


This work received financial support from Fogarty Global Infectious Diseases Training Grant D43TW007391, NIH R01 award AI068932 and 5T32AI007404-24. We are grateful to Carol Mariani and Augustine Dunn of the Yale Caccone laboratory for help with sample processing and laboratory analysis. Special thanks also goes to Mary Burak for assistance with ArcGIS and for producing the map in Fig. 1. The research was accomplished while RO was a Fogarty Research Fellow at Yale University.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Robert Opiro.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Collected field samples: RO, EAO and RE. Performed the lab work: RO, KD and AH. Analyzed the data: RO and NS. Wrote the paper: RO, NS and AC. Critically revised manuscript: SA, EAO and RE. Contributed reagents/materials/analysis tools: AC, SA. Conceived and designed study: AC, SA, RE, RO and EAO. All authors read and approved the final version of the manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Opiro, R., Saarman, N.P., Echodu, R. et al. Evidence of temporal stability in allelic and mitochondrial haplotype diversity in populations of Glossina fuscipes fuscipes (Diptera: Glossinidae) in northern Uganda. Parasites Vectors 9, 258 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: