Evidence of temporal stability in allelic and mitochondrial haplotype diversity in populations of Glossina fuscipes fuscipes (Diptera: Glossinidae) in northern Uganda
© Opiro et al. 2016
Received: 6 January 2016
Accepted: 20 April 2016
Published: 3 May 2016
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” .
HAT is among the most debilitating diseases on the continent, with an estimated 70 million people at risk . 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 . 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 . Overall, annual losses in agriculture alone have been estimated at over US $5 billion .
In Uganda, about two thirds of the total land area is infested with tsetse flies , 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 . There is a growing body of evidence that the disease ranges are expanding , 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  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 , 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 , 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 . Studying temporal patterns of genetic variation can provide insight into some of these knowledge gaps and the apparent contradictions between ecological and genetic data . Temporal population genetics data can also provide estimates of effective population sizes , and probability of population bottlenecks, and thereby aid in the choice of spatial and temporal parameters in vector control programmes .
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 [26–28], 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 . 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 . 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
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 (http://dna-analysis.yale.edu/). 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  to test for deviation from Hardy Weinberg equilibrium (HWE), linkage disequilibrium (LD), and to estimate inbreeding coefficients (FIS) using the Markov chain method  with 100,000 dememorizations, 1,000 batches and 10,000 iterations per batch. Arlequin v 3.5  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 . We used the program FSTAT v 3.1  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  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 (http://dna-analysis.research.yale.edu/). 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  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  using DEMEtics  in R . 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 , 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  based on , and the linkage disequilibrium (LD) method, as implemented in NeEstimator v 2.01 . 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 . The first method tested for an excess of heterozygosity relative to observed allelic diversity , 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 , and significance was assessed using Wilcoxon’s signed rank test, as recommended for fewer than 20 loci . The second method tested for a bottleneck-induced mode shift in allele frequency distributions that were based on equal increments of 0.1 .
Sample sizes and genetic diversity statistics for 16 nDNA microsatellite loci and an mtDNA marker
Date of sampling
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 testing for temporal genetic structure
Sum of squares
2 temporal groups (generation 0 & X)
Among temporal groups
Among sites within groups
8 temporal groups (generation 0, 22, 27, 29, 35, 38, 48, 52)
Among sites within groups
2 temporal groups (generation 0 & X)
Among temporal groups
Among sites within groups
8 temporal groups (generation 0, 22, 27, 29, 35, 38, 48, 52)
Among temporal groups
Among sites within groups
Pairwise estimates of genetic differentiation between temporal samples
KT-0 vs KT-22
OC-0 vs OC-27
AM-0 vs AM-27
KR-0 vs KR-35
AP-0 vs AP-48
KO-0 vs KO-29
AR-0 vs AR-52
PD-0 vs PD-52
MY-0 vs MY-38
Effective population size (Ne) estimates and tests for bottlenecks
Effective population size estimates and tests for population bottleneck events
Effective population size
95 % CI
95 % CI
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 . 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 . 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 , G.palpalis gambiensis , and even in G. f. fuscipes in Uganda . 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  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 , and similar to estimates obtained by  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 . 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 . Furthermore, the large sampling intervals used decreases the bias caused by overlapping generations and age structure . 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 . 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 , and later control activities were decentralized and are currently managed by underfunded district health and veterinary authorities . 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 , 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 . 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.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- 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.PubMedGoogle Scholar
- 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. Google Scholar
- Kabayo JP. Aiming to eliminate tsetse from Africa. Trends Parasitol. 2002;18:473–5.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- Africa E. Pan African Tsetse and Trypanosomosis Eradication Campaign (PATTEC) Enhancing Africa’ S Prosperity Plan of Action Organization of African Unity. 2001.Google Scholar
- 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.Google Scholar
- PAAT. Programme Against African Trypanosomiasis, Newsletter No. 8. 2000.Google Scholar
- 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
- Chizyuka G. FAO Liaison Officers Summary Report. Harare: Food and Agricultural Organisation; 1998.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Welburn SCOM. Recent developments in human African trypanosomiasis. Curr Opin Infect Dis. 2002;15:477–84.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- Hargrove JW. Tsetse population dynamics. In: Maudlin I, Holmes PH MM, editors. The Trypanosomiases. Wallingford: CABI Publishing; 2004.Google Scholar
- 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.Google Scholar
- Ouma JO. Population genetic studies of the tsetse fly Glossina pallidipes (Diptera: Glossinidae). PhD thesis, Iowa State University, USA. 2004.Google Scholar
- Richards CLP. Temporal changes in allele frequencies and population’s history of severe bottlenecks. Conserv Biol. 1996;10:832–9.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Basalirwa CP. Delineation of Uganda into climatological rainfall zones using the method of principal component analysis. Int J Climatol. 1995;15:1161–77.Google Scholar
- 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
- 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.Google Scholar
- Krafsur ES, Marquez JGOJ. Structure of some East African Glossina fuscipes fuscipes populations. Med Vet Entomol. 2008;22:222–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.View ArticlePubMedGoogle Scholar
- TE Guo SW. Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992;48:361–72.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y. Benjamini and Y FDR.pdf. J R Stat Soc Ser B Methodol. 1995;57(No. 1):289–300.Google Scholar
- Goudet J, Raymond M, de Meeüs TRF. Testing differentiation in diploid populations. Genetics. 1996;144:1933–40.PubMedPubMed CentralGoogle Scholar
- 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
- JOST L. G ST and its relatives do not measure differentiation. Mol Ecol. 2008;17:4015–26.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2011.Google Scholar
- Hendricks P. A standardized genetic differentiation measure. Evolution (N Y). 2005;59:1633–8.Google Scholar
- Meirmans PG. The trouble with isolation by distance. Mol Ecol. 2012;21:2839–46.View ArticlePubMedGoogle Scholar
- Waples R. A generalized approach for estimating effective population size. Genetics. 1989;379–391.Google Scholar
- Jorde PERN. Unbiased estimator for genetic drift and effective population size. Genetics. 2007;177(2):927–35.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Nei M, Tajima F. Genetic drift and estimation of effective population size. Genetics. 1981;98:625–40.PubMedPubMed CentralGoogle Scholar
- Pollak E. Population size from allele frequency changes. 1983:531–548.Google Scholar
- Hill WG. Estimation of effective population size from data on linkage disequilibrium. Genet Res (Camb). 1981;38:209–16.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- Marie J, Luikartt G, Cornuet JM, Luikartt G. Power analysis of two tests for detecting recent population bottlenecks from allele frequency data. 1996.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Luikart G, Allendorf FW, Sherwin WB. Distortion of allele frequency distributions bottlenecks. 1998:238–47.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Turner TF, Salter LA, Gold JR. Temporal-method estimates of N e from highly polymorphic loci. Conserv Genet. 2001;2:297–308.View ArticleGoogle Scholar
- 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.Google Scholar
- PRDP. Peace, Recovery and Development Plan for Northern Uganda (PDRP)-2007-2010. 2007.Google Scholar
- Menon S, Rossi R, Nshimyumukiza L, Zinszer. Revisiting zoonotic human African trypanosomiasis control in Uganda. Journal of Public Health Policy. 2016;37(1):51.Google Scholar
- Trends PR. A Climate Trend Analysis of Uganda. 2010.Google Scholar