Search tips
Search criteria 


Logo of plosntdPLoS Neglected Tropical DiseasesView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS Negl Trop Dis. 2010 May; 4(5): e692.
Published online 2010 May 25. doi:  10.1371/journal.pntd.0000692
PMCID: PMC2876113

Population Genetics as a Tool to Select Tsetse Control Strategies: Suppression or Eradication of Glossina palpalis gambiensis in the Niayes of Senegal

Joseph Mathu Ndung'u, Editor



The Government of Senegal has initiated the “Projet de lutte contre les glossines dans les Niayes” to remove the trypanosomosis problem from this area in a sustainable way. Due to past failures to sustainably eradicate Glossina palpalis gambiensis from the Niayes area, controversies remain as to the best strategy implement, i.e. “eradication” versus “suppression.” To inform this debate, we used population genetics to measure genetic differentiation between G. palpalis gambiensis from the Niayes and those from the southern tsetse belt (Missira).

Methodology/Principal Findings

Three different markers (microsatellite DNA, mitochondrial CO1 DNA, and geometric morphometrics of the wings) were used on 153 individuals and revealed that the G. p. gambiensis populations of the Niayes were genetically isolated from the nearest proximate known population of Missira. The genetic differentiation measured between these two areas (θ = 0.12 using microsatellites) was equivalent to a between-taxa differentiation. We also demonstrated that within the Niayes, the population from Dakar – Hann was isolated from the others and had probably experienced a bottleneck.


The information presented in this paper leads to the recommendation that an eradication strategy for the Niayes populations is advisable. This kind of study may be repeated in other habitats and for other tsetse species to (i) help decision on appropriate tsetse control strategies and (ii) find other possible discontinuities in tsetse distribution.

Author Summary

Tsetse flies transmit trypanosomes to humans (sleeping sickness) and animals (nagana). Controlling these vectors is a very efficient way to contain these diseases. There are several strategies and methods that can be used for control, each being more or less efficient depending on several factors. The Government of Senegal wants to sustainably eliminate trypanosomosis from the Niayes region by controlling the tsetse vector, Glossina palpalis gambiensis. To reach this objective, two different strategies may be used: suppression (decrease in tsetse densities) or eradication (remove all the tsetse in the region until last one). For eradication, the approach has to be area-wide, i.e. the control effort targets an entire pest population within a circumscribed area, to avoid any possible reinvasion. Three different tools (microsatellite DNA, mitochondrial DNA and morphometrics) were used, and all showed an absence of gene flow between G. p. gambiensis from the Niayes and from the nearest known population in the south east of the country (Missira). This genetic isolation of the target population leads to the recommendation that an eradication strategy for the Niayes populations is advisable. This kind of study may be extended to other areas on other tsetse species.


The Niayes of Senegal harbours the most northern and western population of Glossina palpalis gambiensis Vanderplank, which is a major vector of the debilitating diseases Human African Trypanosomosis (HAT) or sleeping sickness, and African Animal Trypanosomosis (AAT) or nagana (reviewed in [1]). Particular meteorological and ecological characteristics of this area provide great potential for agricultural development in general and animal production (cattle, donkeys, horses, small ruminants, pigs and poultry) in particular. Most of these animals are however susceptible to AAT which is seriously limiting the development of efficient and productive, sustainable livestock systems. The socio-economic impact of the disease is therefore dramatic and very often underestimated [2]. In the 1970s and 1980s, it was attempted to eliminate the G. p. gambiensis population from the Niayes mainly using ground spraying of residual insecticides [3]. The tsetse and trypanosomosis problem seemed to have disappeared until flies were detected again in 1998 (unpublished report of the Direction de l'Elevage - DIREL). In 2005, the DIREL initiated a control campaign called «Projet de lutte contre les glossines dans les Niayes» with the objective of developing a sustainable solution to the tsetse and trypanosomosis problem in the Niayes. The programme is funded by the Government of Senegal and technically and financially supported by the Food and Agriculture Organization of the United Nations (FAO) and the International Atomic Energy Agency (IAEA). The project is implemented in the context of the African Union - Pan African Tsetse and Trypanosomiasis Eradication Campaign (PATTEC), a political initiative of the African heads of state that calls for increased efforts to manage the tsetse and trypanosomosis problem.

Tsetse populations may be reduced using a variety of techniques, including insecticide impregnated traps and targets, live-baits, sequential aerial spraying, and sterile male releases [4][9]. In the past, most control efforts were not implemented according to area-wide principles [10], [11], and as a consequence, when the control effort was reduced or stopped, the tsetse populations tended to recover – due to either flies surviving the initial interventions, or migrant flies coming from untreated regions, or both [12]. This has fuelled a debate as to whether in some instances “eradication”, defined by FAO [13], [14] as the creation of a tsetse free zone, may be more cost effective than “suppression” where tsetse densities are reduced to a level minimizing the risk of disease transmission. A sound decision whether to select an eradication or suppression strategy will be facilitated when the population structure within the target region, in particular the degree of genetic isolation of the target population from its adjacent populations is clearly understood. For isolated populations, eradication may be the most cost-effective strategy, as reported for Glossina austeni Newstead in Unguja Island, Zanzibar [5]. But for most mainland populations of tsetse, the geographical limits of target tsetse populations are less easily defined. Application of population genetics techniques can quantify rates of gene flow between sub-populations [15][19], and guide decisions on the choice of control strategies [20].

The level of isolation of the targeted tsetse populations will be an important parameter to guide the Government of Senegal to select the most optimal control strategy. Here we report population genetic analyses of microsatellite and mtDNA markers combined with morphometrics of G. p. gambiensis populations sampled from the Niayes area and from the nearest population in the south-eastern part of the country (fig. 1) to assess their degree of isolation by measuring gene flow among the different populations. The genetic differentiation of the various G. p. gambiensis populations within the Niayes was also assessed to determine if the different populations of the Niayes can be targeted at the same time (if it is a single panmictic unit), or if a sequential control strategy can be contemplated (if substantial genetic differentiation between populations is found)., which will also depend on their respective history, including effective population sizes and possible bottlenecks.

Figure 1
Geographic location of the study sites in Senegal.


Study sites and sampling

In the Niayes region, four tsetse populations were sampled using Vavoua traps: Dakar Hann which is a swamp forest harbouring an animal park within the city of Dakar, Diacsaw Peul, an area of riparian thicket where tsetse and cattle are in intense contact, Sebikotan and Pout, which are mango and citrus-tree plantations where tsetse and people are in close contact. The tsetse flies collected in the Niaye area and analyzed in the present study are the sole property of the Senegalese authorities. They were collected by the national veterinary services through official mission orders, in one wildlife park (Parc de Hann) and three private sites, with the oral consent of the owners. No written consent is mandatory for tsetse fly collection in Senegal. In the south-eastern part of the country, the area of Missira was sampled: it is the nearest known infested area from the Niayes, according to a detailed tsetse survey implemented as part of the baseline data collection of the Niayes tsetse Control project (see Fig. 1). Areas between the Niayes and Missira are not favourable for tsetse, which was confirmed by zero tsetse catches (JB, BS unpublished data) despite intensive trapping efforts.

In total 153 tsetse individuals were analysed originating from Diacsaw Peul (22 females (F), 8 males (M)), Dakar Hann (23F, 6M), Sebikotan (21F, 11M), Missira (23F, 12M), and Pout (13F, 14M).

Microsatellite DNA markers

Markers used and PCR conditions

A total of 10 microsatellite markers were used (preceded by “X” for X-linked loci): X55.3 [21], XPgp11, Pgp1, XPgp13, Pgp24 [22], A10, XB104, XB110, C102 (A. Robinson, FAO/IAEA, pers. com.) and GPCAG [23]. The samples were then processed for Polymerase Chain Reaction (PCR) and genotyping on a 4300 DNA Analysis System from LI-COR (Lincoln, NE) exactly as described in [19].

Population structure analyses

Wright's F-statistics [24], the parameters most widely used to describe population structure [25], were initially defined for a three-levels hierarchical population structure (individuals, sub-populations and total). In such a structure, three fixation indices or F-statistics can be defined. F IS is a measure of the inbreeding of individuals (hence I) resulting from non random union of gametes within each sub-population (hence S). F ST is a measure of the relatedness between individuals resulting from non-random distribution of individuals among sub-populations, relative to the total population; F ST quantifies the differentiation between sub-populations in the total population (hence S and T). F IT is a measure of the inbreeding of individuals resulting both from non-random union of gametes within sub-populations and from population structure (deviation from panmixia of all individuals of the total population, hence I and T). These F-statistics are classically estimated by Weir and Cockerham 's unbiased estimators f (for F IS), θ (for F ST) and F (for F IT) [26]. When appropriate, these statistics were estimated with Fstat (updated from [27]). However, more than two levels (i.e., individuals, sub-populations and total) may exist. This is the case for the different subsamples belonging to the Niayes region when to be compared to tsetse from Missira. Hierfstat version 0.03–2 [28] is a package for the statistical software R [29]. This package computes hierarchical F-statistics from any number of hierarchical levels [28]. A user-friendly description of the software is presented elsewhere [30].

The significance of the F-statistics was tested by randomization (10000 permutations in each case). The significance of F IS was tested by randomizing alleles between individuals within sub-samples. The significance of F ST was tested by randomizing individuals among sub-samples. These tests were performed with Fstat. The significance of F SR, i.e., the homozygosity due to subdivision into subsamples within the Niayes region, was tested by randomising individuals among subsamples of the Niayes. The significance of F RT, which measures the relative homozygosity due to the geographical separation between the Niayes region and Missira was tested by randomizing subsamples (with all individuals contained) between the two regions. These two last tests were performed with Hierfstat. For F IS the statistic used was directly the f (unbiased estimator of F IS). For other differentiation measures (F ST, F SR, F RT), the statistic used was the maximum likelihood ratio G [28], [31].

Non random association of alleles within individuals (F IS>0) may be due to the presence of null alleles. We used Micro-Checker 2.2.3 [32] to detect null alleles and estimate their frequency pn at each locus according to Brookfield's second method [33]. For X-linked loci, males were not included in such analyses and null frequency was also directly estimated by the proportion of missing genotypes. The global expected frequency of blanks (null homozygotes) was then compared to the observed one with an exact unilateral (H1: there are less blanks than expected) binomial test with R.

Linkage disequilibrium between pairs of loci was tested using Fstat by randomising loci combinations (free recombination) across loci with a G-based statistics permitting a global test across sub-samples for each pair of loci. Because there are as many tests as locus pairs tested, we used the exact binomial test to check if the proportion of tests found significant at the 5% level was significantly above 0.05 with the “Test alpha' and look for k'” option of Multitest V.1.2 [34]. To optimise power we disregarded loci displaying one allele at a frequency above 0.9, which was the case of locus XB110.

Some loci are X-linked. For F IS based analyses males were coded as missing data at these loci, otherwise these loci were coded as homozygous (differentiation and LD based analyses).

In order to visualize the partition in genetic differentiation among all samples, a NJTree dendrogram was constructed with the software MEGA 3.1 [35] from a chord distance matrix (noted DC-S&E) [36], as recommended by [37], and computed using Genetix V4 [38].

Effective population sizes

Effective population sizes were computed with three methods. Estim 1.2 [39] uses the connection between migration and effective population size with heterozygosity and linkage disequilibrium between loci. Linkage disequilibrium is indeed connected to effective population size, a property exploited by the method from Bartley et al. [40], which is in fact coming from Hill [41] and modified by Waples [42]. From equation 2 of [40], modified as in NeEstimator Help file [43], allelic correlations (r) can be estimated as An external file that holds a picture, illustration, etc.
Object name is pntd.0000692.e001.jpg, where p and q are frequencies of allele A at locus 1 and allele B at locus 2, respectively and D is the so called (e.g. [40], [44]) Burrow's composite measure of disequilibrium [45]. The r 2-values across all pairs of alleles are averaged to yield a single r 2 for each pair of loci. Finally, an arithmetic mean of the r 2 values for all pairs of loci is used to obtain a single correlation coefficient and to obtain an Ne estimate using the equation An external file that holds a picture, illustration, etc.
Object name is pntd.0000692.e002.jpg where An external file that holds a picture, illustration, etc.
Object name is pntd.0000692.e003.jpg is the harmonic mean of the sample sizes of each pairwise comparison between loci [40], [46]. Males were coded homozygous for X-linked loci, which should not be a problem for a composite based linkage disequilibrium measure. The method was implemented with NeEstimator [43]. Heterozygote excess method from Pudovkin et al. [47] (see also [48]) corrected by Balloux [49], uses the fact that, in dioecious (or self incompatible) populations, alleles from females can only combine with alleles contained in males and a heterozygote excess is expected as compared to Hardy-Weinberg expectations, and this excess is proportional to the effective population size. This method was implemented using Weir and Cockerham estimator of F IS in the equation Ne = 1/(−2F IS)−F IS/(1+F IS) [49] and was only applicable in subsamples with heterozygote excess.

For intra locus based methods (Estim and Balloux's methods), loci showing evidence of stuttering and null alleles were removed.

Bottleneck detection

Signatures of bottleneck events were investigated by comparing the expected heterozygosity for a sample (HE) with the heterozygosity that would be expected for a sample taken in a population at mutation/drift equilibrium with the same size and allele number (HEQ). As allele number decreases faster than heterozygosity, a bottleneck is signed by HE > HEQ in subsequent generations [50]. This analysis was performed with Bottleneck v1.2 software [51] assuming that mutations of microsatellite loci followed either an IAM (infinite allele model), a SMM (stepwise mutation model), or a TPM (two phase model), in the last case of which we assumed that 70% of mutations consist of one step and 30% consist of multistep change with a variance of 30 (default values). Tests were performed using unilateral Wilcoxon tests as recommended [50].

The method of bottleneck detection described by Cornuet and Luikart [50] allows a rough estimate of the effective population size right after the bottleneck event (Ne -pb). From the Figure 3A of [50], with 10 loci and mean sample sizes of 10–30 individuals, a bottleneck can be detected if it occurred between τ 12Ne -pb and τ22Ne -pb generations before sampling with τ 1 = 0.1 and τ 2 = 2.5.

Figure 3
Unrooted Cavalli-Sforza and Edward's chord distance based Neighbour-Joining-Tree.

The difference in genetic diversity, as measured by Hs (Nei's unbiased estimator [52]), between the Niayes and Missira was evaluated with a bilateral Wilcoxon signed rank test for paired data. The different microsatellite loci were used as the pairing factor. This test was undertaken with R.

The Bonferronni procedure [53] was used each time multiple testing was done and individual tests significance required.

User-friendly descriptions of most of the tests and procedures used in the genetic data analyses can be found in a recent review [54].

Mitochondrial DNA markers

A portion of the 5′ end of the mitochondrial gene COI was amplified, purified and sequenced using the primers CI-J-2195 TTGATTTTTTGGTCATCCAGAAGT [55] and CULR TGAAGCTTAAATTCATTGCACTAATC using the same conditions reported by [56]. The following statistics were calculated using DNAsp version 4.50.3 [57]: Hd Haplotype diversity, was calculated using equations 8.4 and 8.12 in [58], Pi, the nucleotide diversity, which is the average number of nucleotide differences per site between two sequences, and its sampling variance was calculated using equations 10.5–10.7 in [58], K, the average number of nucleotide differences and the total variance of K is (sampling plus stochastic), assuming no recombination were calculated using equations from [59]. FST was calculated according to equation 3 in [60], HST was calculated according to equation 2–4, and K*ST according to equations 7–11 in [60]. HST is a haplotype frequency based genetic differentiation statistic that does not take into account the number of differences separating different haplotypes. HST = 1−(HS/HT), where HS is the weighted average of subpopulation genetic diversities and H T is the estimated haplotype diversity of the total population. K*ST = 1−(KS/KT), where K*S is a weighted average of the log corrected average number of sequence difference in the populations being compared, and KT is the average number of difference between sequences. A permutation test, in which haplotypes or sequences were randomly assigned to the different localities 10000 times, was used to test the significance of HST and K*ST) [60]. The average number of nucleotide difference (equation A3, [61]) and its variance (sampling plus stochastic) were also calculated using DNAsp version 4.50.3. Haplotype trees for 738 nucleotides of COI from G. p. gambiensis (data from a total of 148 individuals in this study) were generated using the algorithm of [62]. The TCS 1.21 programme was used to estimate the haplotype tree, with the connection limit (probability of parsimony) at 95% [63]. The maximum number of connection steps at 95% was 11.

Geometric morphometrics

In total, the number of analysed wings was 20, 21, 18, 34 and 18 for Diacsaw Peul, Dakar Hann, Sebikotan, Missira and Pout respectively. These analysed individuals were all also analysed by the molecular markers.

Wings were dry-mounted between two microscope slides and scanned with a scanner. From this picture, 10 landmarks defined by vein intersections were recorded as previously described [64]. Each landmark has X and Y coordinates, and the 10 LM defined per wing represent a polygon. After scaling, translating and rotating all these polygons so that they can be compared, data were subjected to generalized Procrustes analysis (GPA) [65], [66] allowing to implement shape variables, here represented by 16 “partial warps” (PW) (including uniform component of shape). These PW were used to conduct a discriminant analysis to allow for individual reclassification based on Mahalanobis distances (noted DM), which were calculated between populations. The statistical significance of Mahalanobis distances was estimated by 1,000-runs permutation tests [67]. A re-classification score was computed where individuals are assigned to each group and the percentage of good classification was then calculated.

Evaluating correlation between genetic distances and morphometric distance matrices

To evaluate the correlation between distance matrices we undertook three Mantel tests [68] with the “Mantelize it” option of Fstat 2.9.4 [27] between DC-S&E and Kst*, between DC-S&E and DM, and between Kst* and DM.


Microsatellite loci

Within subsamples analyses

The per-locus number of alleles was 3.4, 3.7, 3.9 and 4 for Diacsaw Peul, Dakar Hann, Sebikotan and Pout respectively (all coming from the Niayes), but it was 8.7 for Missira. The Wilcoxon test confirmed that mean genetic diversity in the Niayes (Hs = 0.5) was lower than the one found in Missira (Hs = 0.7) (P-value = 0.0059).

There was a global and significant homozygote excess over all loci (F IS = 0.17, P-value = 0.0001). Figure 2 shows that this is mainly due to certain loci (i.e. pGp1, pGp24, A10 and XB110). This strongly suggests locus specific technical problems. Microchecker analysis confirmed this view as only A10 and XB110 did not display enough blanks (not enough male blanks for XB110) as compared to expected ones if F IS were only due to null alleles. However, both loci seemed to suffer from stuttering. Without these four loci, the six remaining loci did not result in rejection of the hypothesis of panmixia (F IS = 0.035, P-value = 0.0983).

Figure 2
F IS across loci and overall loci.

Among the 36 possible tests, four pairs of loci displayed a significant linkage at the 5% level of significance. This proportion (4/36) is not different from the proportion expected under the null hypothesis (exact binomial test, P-value = 0.1037).

Hierarchical population structure

The overall mean F ST value was estimated at θ = 0.11 (CI95: 0.088<θ<0.143). The mean θ value between Missira and all the Niayes populations was high and significant, θ = 0.12 (CI95: 0.08<θ<0.14, p<0.001). Within the Niayes, the mean θ value between Dakar Hann and the three others was also high and significant θ = 0.12 (CI95: 0.05<θ<0.20, p<0.001), which is of the same order of magnitude as the one between Missira and the Niayes. This is illustrated in Fig. 3, and this explains why HierFstat analysis reveals a substantial and highly significant contribution of sites (F SR = 0.069, P-value = 0.001) and an absence of significant contribution of region (P-value = 0.209). This latter point is due to the fact that Dakar-Hann is as different from the other Niayes subsamples, than Niayes are from the remote Missira subsample.

Effective population sizes

A constant seemed to be that Diacsaw Peul, Dakar-Hann and Sebikotan displayed very small effective population sizes (3 to 28), and Pout and Missira appeared much bigger (51, 83 respectively). However, some of the computations led to infinite population sizes and sometimes the different methods did not converge to a reliable estimator (see Table S1), which suggests bigger effective population sizes than those found.


Only two subsamples provided some evidence for a bottleneck: Dakar-Hann with P-values of 0.001, 0.001 and 0.116 for IAM, TPM and SMM mutation models respectively; and Pout with P-values 0.012, 0.348 and 0.935 for IAM, TPM and SMM (see Table S2). These values were still significant after Bonferronni procedure for multiple testing in Dakar Hann only: P = 0.015 and 0.022 for IAM and TPM respectively. It is thus probable that a bottleneck affected the population from Dakar-Hann, between 0.2Ne -pb and 5Ne -pb generations ago. Using 16 as a proxy for the effective population size in this site (obtained with Bartley's method) ranges the bottleneck between 3.2 and 80 generations ago. With the 95% confidence interval upper limits (i.e. 23) the ranges increases to [5–115]. This bottleneck provides a good explanation for the genetic isolation noticed for this site during the population differentiation analysis.

Mitochondrial DNA

The statistics of genetic differentiation between the Niayes population samples from Sebikotan, Pout, Dakar Hann and Diasca Peuhl and the south-eastern Senegal site Missira based on 738 nucleotides of mitochondrial cytochrome oxidase 1 (COI) are shown in table 1. Statistics of genetic differentiation based on nucleotide sequence and haplotype frequency both indicate a very high level of genetic differentiation between Niayes populations and Missira, since there were no common haplotypes between the two regions (see Fig 4). Genetic differentiation between Diacsaw Peul and Pout was also significant (0.05>P>0.01). Otherwise, there was no significant evidence for differentiation between populations within the Niayes region. The haplotype and nucleotide diversity of the Niayes populations was very low in all cases (see table 2).

Figure 4
Haplotype trees for 738 nucleotides of COI from G. p. gambiensis.
Table 1
Population differentiation statistics for mitochondrial cytochrome oxidase 1 (COI) sequences (738 sites in alignment, no gaps).
Table 2
Diversity statistics for the studied populations based on COI sequences.


An analysis of the shape of the wings was conducted by using Procustes analysis to compare wing shape among the five different populations. All Mahalanobis distances between populations were significant except between Pout and Sebikotan. This is illustrated by reclassification scores which were over 70% indicating differentiation between populations (see Table 3). Missira sample showed the best one (85%) confirming this sample was easily distinguished from the others.

Table 3
Geometric morphometrics: reclassification scores of individuals in their population of origin.

Correlation between distance matrices

The correlation between the genetic distances and morphometrric distances were all strong and significant (see Figure 5). The strong correlation between microsatellites and COI distance matrices was mainly due to the fact that there are two sets of points: those from the Niayes sites and those differentiating the Niayes from Missira (very strong distances for both markers). For the correlation between genetic and morphometric data, Figure 5 shows that around 50% of the variance in Mahalanobis distance can be explained by genetics, meaning that the other 50% are probably explained by environmental parameters.

Figure 5
Mantel correlation tests (P-values) between the different markers used in this study.


The results of the genetic and morphometric analyses indicated limited gene flow between the G. p. gambiensis population of the Niayes and those of the main tsetse belt in the south-eastern part of Senegal. Using three different kinds of markers, i.e. microsatellite DNA, mitochondrial DNA and geometric morphometrics, the data led to the conclusion that the G. p. gambiensis population from the Niayes can be considered isolated with very little risk of re-invasion should the population eventually be eradicated. These data corroborate the results of the entomological baseline data collection and observations in the field indicating the absence of G. p. gambiensis in the 120 km long area between Missira and the Niayes (JB, BS unpublished results.). On this basis, eradication of G. p. gambiensis from the Niayes can therefore be recommended as an appropriate control strategy. The data from this study, which was part of a comprehensive baseline data collection effort, confirmed that eradication can be recommended as an appropriate control strategy, and as such the study greatly assisted in the decision making on which strategy to select. Although the isolation of the target population is not an absolute prerequisite for AW-IPM, tackling a continuous pest populations is actually more complex, requiring more resources and a long-term, regional commitment (for details, see [10], [11], [14].

Looking at F ST values reported in previous studies on G. p. gambiensis, the one observed in the present work between Missira and the Niayes populations was ten times higher than those observed along 260 km on the Mouhoun river in Burkina Faso [17], and two times higher than the values observed in Guinea between Loos islands and the continent [64]. These values are of the same order of magnitude as those observed between the two different taxa, G. p. gambiensis and G. palpalis palpalis (DK, PS, unpublished data). It is also noteworthy that some of the microsatellite loci used on the individuals in this study amplified poorly, a behaviour not recorded in earlier studies of G. p. gambiensis [17], [19], [64]. This may be an additional argument for genetic divergence, since it is known that in different taxa polymorphisms in sequences flanking the microsatellite may occur [69] leading then to mismatches in the primer binding sites. If gene flow between the sites occurs then shared haplotypes of the mitochondrial gene COI would be expected. In fact, we observed no haplotypes shared between Missira and the Niayes region. Future work should include more detailed examination of potential cryptic species within the distribution area of G. p. gambiensis since such cryptic species have been suspected and shown to exist in closely related taxa [70], [71].

In case the Government of Senegal decides to integrate the sterile insect technique (SIT) as part of the AW-IPM (area-wide integrated pest management) approach, it would be advisable that prior to the release of sterile males, experiments be conducted to assess mating compatibility between the G. p. gambiensis from the Niayes and those which should be released if they do not originate from the same area (e.g. from Burkina Faso). Although previous studies with G. p. gambiensis originating from Mali and Burkina Faso revealed no mating barriers between these populations (G. Mutika, personal communication), nor did another study between G. p. gambiensis and G. p. palpalis [72], these tests would provide evidence that genetic differentiation of the Niayes population has not been accompanied by pre-mating barriers, that would threaten the success of the programme. Similarly, it may also be recommended to release sterile flies coming from Burkina Faso in the Niayes to assess their behaviour and performance (dispersal, dispersion, mobility, lifespan, mating frequency etc.) in the natural habitat. Indeed, the particular eco-system of the Niayes, the fragmented nature of the preferred vegetation types of G. p. gambiensis (mango-tree plantations, Euphorbia spp. hedges) and an annual precipitation below 400–500 mm might have induced the development of some level of xerotolerance of these relatively small populations [73], [74].

The second question that has relevance for a future AW-IPM programme was to know whether within the Niayes region, the four populations showed any genetic differentiation, or if they constituted a single, panmictic (i.e random mating) population. Low haplotype diversity of the populations of this area is in agreement with a previous Single Strand Conformation Polymorphism-based study on mitochondrial DNA haplotypes in the Niayes region of Senegal [75]. The Niayes population is probably a remnant population, possibly of small size and therefore likely to lose rare haplotypes by genetic drift more rapidly than larger populations. Many ancestral G. p. gambiensis haplotypes have probably been lost from the Niayes region, leaving just two haplotypes remaining today. The low diversity was also observed at microsatellite loci with lower genetic diversity than in Missira. The three different markers used in this study generally showed good agreement in the differentiations observed, as can be seen by the high and significant correlations coefficients computed between them.

The only apparent discrepancy between microsatellites on one hand and COI based and morphometrics based distances on the other hand probably comes from the variance expected for paired F ST estimates and the small sample sizes, the effect of which is probably increased by historical events such as the probable bottleneck that affected the Dakar-Hann population. This may explain why only microsatellite data suggest that the population from Dakar Hann is so much differentiated from all the others, including the geographically closest ones. From an operational point of view however, the population of Dakar Hann is probably also isolated from the three others, since the bottleneck signature was still visible, suggesting very limited exchanges (if any) with the other populations. Field data also corroborate this since the population of Dakar Hann is located within an animal park of this huge city which can be seen as an isolated refuge for this tsetse population. This further suggests that an elimination operation may probably be implemented on this population.

The results also suggest that the samples from Diacsaw Peul, Sebikotane and Pout, although showing some genetic differentiation, are not completely isolated from each other. This is also consistent with ecological data since it appears that the maximum distance between the forest patches of this area is less than 2km (data not shown.).

To conclude, the use of genetic and morphometric markers has been instrumental in the decision-making process of selecting and developing of an appropriate intervention strategy to create a sustainably zone free of G. p. gambiensis and Trypanosomosis in the Niayes region of Senegal. In the near future, it should be encouraged to carry out such studies prior to the selection of target areas or the choice of control strategies, and these should be part of the overall collection of baseline data (see ref [76] for recent review). In addition, the results obtained here suggest that efforts should be made to look for other genetic discontinuities in G. palpalis s.l. distribution that may be indicative of the presence of cryptic species.

Supporting Information

Alternative Language Abstract S1

Translation of the abstract into French by PS.

(0.03 MB DOC)

Table S1

Effective population size estimates according to the different methods described.

(0.02 MB XLS)

Table S2

Probabilities of detecting bottleneck in each population according to microsatellite mutation model.

(0.02 MB XLS)


The authors acknowledge the Direction de Services Vétérinaires for their logistic support, and the technicians of the DSV who assisted with the sampling of the flies, especially Abdoulaye Mbaye, Alphonse Manga, Aliou Babel and Babacar Camara. TDM thank Vince Obre for his strong support during data analyses.


The authors have declared that no competing interests exist.

The authors acknowledge the support of the EU-funded INCO-DEV TFCASS project, IAEA Technical Cooperation project SEN5029 for financing the mission of JB to guide the sampling effort, the FAO/IAEA Joint Programme of Nuclear Techniques in Food and Agriculture for funding part of the genotyping. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Solano P, Bouyer J, Itard J, Cuisance D. Lefèvre PC, Chermette J, Blancou R, Uilenberg G, editors. Cyclical vectors of Trypanosomosis. Infectious and Parasitic Diseases of Livestock, Lavoisier Tec&Doc. Chap 12. 2009a. pp. 153–181. in press.
2. Budd L. DFID-funded tsetse and trypanosomiasis research and development since 1980. Chatham, UK: Department for International Development: Livestock Production Programme, Animal Health Programme/Natural Resources Systems Programme; 1999. Economic analysis.
3. Touré SM. Lutte contre Glossina palpalis gambiensis dans la région des niayes du Sénégal. Rev Elev Méd vét Pays trop. 1973;26:339–347. [PubMed]
4. Challier A, Eyraud M, Lafaye A, Laveissière C. Amélioration du rendement du piège biconique pour glossines (Diptera: Glossinidae) par l'emploi d'un cône inférieur bleu. Cah ORSTOM sér Ent Méd Parasitol. 1977;15:283–286.
5. Vreysen MJB, Saleh KM, Ali MY, Abdulla AM, Zhu ZR, et al. Glossina austeni (Diptera: Glossinidae) eradicated on the island of Unguja, Zanzibar, using the Sterile Insect Technique. J Econ Entomol. 2000;93:123–135. [PubMed]
6. Cuisance D, Itard J, Solano P, Desquesnes M, Frézil JL, et al. Editions Médicales Internationales Tec et Doc. Principales maladies infectieuses et parasitaires du bétail. Europe et Régions chaudes. Paris, France: Lavoisier; 2003. Trypanosomoses. Méthodes de lutte. pp. 139–165.
7. Kgori PM, Modo S, Torr SJ. The use of aerial spraying to eliminate tsetse from the Okavango Delta of Botswana. Acta Trop. 2006;99:184–199. [PubMed]
8. Bouyer J, Solano P, Cuisance D, Itard J, Frézil JL, et al. Control methods in Trypanosomosis. In: Lefèvre PC, Chermette J, Blancou R, Uilenberg G, editors. Infectious and Parasitic Diseases of Livestock, Chap 127. Lavoisier Tec&Doc; 2009a. in press.
9. Torr SJ, Solano P. Olfaction in tsetse hosts interactions. In: Knols B, Takken W, editors. Olfaction in vector hosts interactions. Netherlands: Wageningen University; 2009. In press.
10. Klassen W. Area-wide integrated pest management and the sterile insect technique. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile insect technique, Principles and practice in area-wide integrated pest management. Dordrecht, The Netherlands: Springer; 2005. pp. 39–68.
11. Vreysen M, Robinson AS, Hendrichs J. Area-Wide Control of Insect Pests, from research to field implementation. Dordrecht, The Netherlands: Springer; 2007.
12. Hargrove JW. Tsetse eradication: sufficiency, necessity and desirability. UK: Centre for Tropical Veterinary Medicine, University of Edinburgh; 2003.
13. Food and Agriculture Organization of the United Nations. Glossary of phytosanitary terms. ISPM, Publication Number 5, International Plant Protection Convention (IPPC) Rome, Italy: FAO; 2005.
14. Hendrichs J, Vreysen MJB, Enkerlin WR, Cayol JP. Strategic options in using sterile insects for area-wide integrated pest management. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile insect technique, Principles and practice in area-wide integrated pest management. Dordrecht, The Netherlands: Springer; 2005. pp. 563–600.
15. Solano P, de La Rocque S, de Meeûs T, Cuny G, Duvallet G, et al. Microsatellite DNA markers reveal genetic differentiation among populations of Glossina palpalis gambiensis in the agropastoral zone of Sideradougou, Burkina Faso. Insect Mol Biol. 2000;9:433–439. [PubMed]
16. Gooding RH, Krafsur ES. Tsetse genetics: contributions to biology, systematics, and control of tsetse flies. Annu Rev Entomol. 2005;50:101–123. [PMC free article] [PubMed]
17. Bouyer J, Ravel S, Dujardin JP, de Meeüs T, Vial L, et al. Population structuring of Glossina palpalis gambiensis (Diptera: Glossinidae) according to landscape fragmentation in the Mouhoun river, Burkina Faso. J Med Entomol. 2007;44:788–795. [PubMed]
18. Bouyer J, Balenghien T, Ravel S, Vial L, Sidibé I, et al. Population sizes and dispersal patterns of tsetse flies: rolling on the river? Mol Ecol. 2009b;18:2787–2797. [PubMed]
19. Solano P, Ravel S, Bouyer J, Camara M, Kagbadouno MS, et al. Population Structures of Insular and Continental Glossina palpalis gambiensis in Littoral Guinea. PLoS Negl Trop Dis. 2009b;3(3):e392. doi: 10.1371/journal.pntd.0000392. [PMC free article] [PubMed]
20. Kagbadouno M, Camara M, Bouyer J, Hervouet JP, Morifaso O, et al. Tsetse elimination: its interest and feasibility in the historical sleeping sickness focus of Loos islands, Guinea. Parasite. 2009;16:29–36. [PubMed]
21. Solano P, Duvallet G, Dumas V, Cuisance D, Cuny G. Microsatellite markers for genetic population studies in Glossina palpalis (Diptera: Glossinidae). Acta Trop. 1997;65:175–180. [PubMed]
22. Luna C, Bonizzoni MB, Cheng Q, Aksoy S, Zheng L. Microsatellite polymorphism in the tsetse fies (Diptera: Glossinidae). J Med Entomol. 2001;38:376–381. [PubMed]
23. Baker MD, Krafsur ES. Identification and properties of microsatellite markers in tsetse flies Glossina morsitans sensu lato (Diptera: Glossinidae). Mol Ecol Notes. 2001;1:234–236. [PMC free article] [PubMed]
24. Wright S. The interpretation of population structure by F-statistics with special regard to system of mating. Evolution. 1965;19:395–420.
25. Nagylaki T. Fixation indices in subdivided populations. Genetics. 1998;148:1325–1332. [PubMed]
26. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–1370. [PubMed]
27. Goudet J. FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995;86:485–486.
28. Goudet J. HIERFSTAT, a package for R to compute and test hierarchical F-statistics. Mol Ecol Notes. 2005;5:184–186.
29. R-Development-core-team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2008., ISBN 3-900051-07-0.
30. De Meeûs T, Goudet J. A step-by-step tutorial to use HierFstat to analyse populations hierarchically structured at multiple levels. Infect Genet Evol. 2007;7:731–735. [PubMed]
31. Goudet J, Raymond M, De Meeûs T, Rousset F. Testing differentiation in diploid populations. Genetics. 1996;144:1933–1940. [PubMed]
32. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–538.
33. Brookfield JFY. A simple new method for estimating null allele frequency from heterozygote deficiency. Mol Ecol. 1996;5:453–455. [PubMed]
34. De Meeus T, Teriokhin AT, GJ F. MultiTest V.1.2, a program to binomially combine independent tests and performance comparison with other related methods on proportional data. BMC Bioinformatics. 2009 in press. [PMC free article] [PubMed]
35. Kumar S, Tamura K, Nei M. MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform. 2004;5:150–163. [PubMed]
36. Cavalli-Sforza LL, Edwards AWF. Phylogenetic analysis: model and estimation procedures. Am J Hum Genet. 1967;19:233–257. [PubMed]
37. Takezaki N, Nei M. Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA. Genetics. 1996;144:3 89–399. [PubMed]
38. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Montpellier (France): Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II; 2004.
39. Vitalis R, Couvet D. ESTIM 1.0: a computer program to infer population parameters from one- and two-locus gene identity probabilities. Mol Ecol Notes. 2001;1:354–356.
40. Bartley D, Bagley M, Gall G, Bentley B. Use of linkage disequilibrium data to estimate effective size of hatchery and natural fish populations. Conserv Biol. 1992;6:365–375.
41. Hill WG. Estimation of effective population size from data on linkage disequilibrium. Genet Res. 1981;38:209–216.
42. Waples RS. Hoelzel AR, editor. Genetic methods for estimating the effective size of cetacean populations. Genetic ecology of whales and dolphins, International Whaling Commission (Special Issue No. 13) 1991. pp. 279–300.
43. Peel D, Ovenden JR, Peel SL. NeEstimator Version 1.3: software for estimating effective population size. Queensland Government, Department of Primary Industries and Fisheries; 2004.
44. Waples RS. A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci. Conserv Genet. 2006;7:167–184.
45. Weir BS. Inferences about linkage disequilibrium. Biometrics. 1979;35:235–254. [PubMed]
46. England PR, Cornuet JM, Berthier P, Tallmon DA, Luikart G. Estimating effective population size from linkage disequilibrium: severe bias in small samples. Conserv Genet. 2006;7:303–308.
47. Pudovkin AI, Zaykin DV, Hedgecock D. On the potential for estimating the effective number of breeders from heterozygote excess in progeny. Genetics. 1996;144:383–387. [PubMed]
48. Luikart G, Cornuet JM. Estimating the effective number of breeders from heterozygote excess in progeny. Genetics. 1999;151:1211–1216. [PubMed]
49. Balloux F. Heterozygote excess in small populations and the heterozygote-excess effective population size. Evolution. 2004;58:1891–1900. [PubMed]
50. Cornuet JM, Luikart G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996;144:2001–2014. [PubMed]
51. 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–503.
52. Nei M, Chesser RK. Estimation of fixation indices and gene diversities. Ann Hum Genet. 1983;47:253–259. [PubMed]
53. Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6:65–70.
54. De Meeûs T, McCoy KD, Prugnolle F, Chevillon C, Durand P, et al. Population genetics and molecular epidemiology or how to “débusquer la bête”. Infect Genet Evol. 2007;7:308–332. [PubMed]
55. Simon C, Frati F, Bechenbach A, Crespi B, Liu H, et al. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc America. 1994;87:651–701.
56. Dyer NA, Lawton SP, Ravel S, Choi KS, Lehane MJ, et al. Molecular phylogenetics of tsetse flies (Diptera: Glossinidae) based on mitochondrial (CO1, 16S, ND2) and nuclear ribosomal DNA sequences, with an emphasis on the palpalis group. Mol Phylog Evol. 2008;49:227–239. [PubMed]
57. Rozas J, Rozas R. DnaSP, DNA sequence polymorphism: an interactive program for estimating population genetics parameters from DNA sequence data. Comput Appl Biosci. 1995;11:621–625. [PubMed]
58. Nei M. Molecular Evolutionary Genetics. New York: Columbia Univ. Press; 1973.
59. Tajima F. Measurement of DNA polymorphism. In: Takahata N, Clark AG, editors. Mechanisms of Molecular Evolution. Sunderland, Massachusetts: Sinauer Associates Inc.; 1993. pp. 37–59.
60. Hudson RR, Boos DD, Kaplan ND. A statistical test for detecting geographic subdivision. Mol Biol Evol. 1992;9:138–151. [PubMed]
61. Tajima F. Evolutionary relationship of DNA sequences in finite populations. Genetics. 1983;105:437–460. [PubMed]
62. Templeton AR, Crandall KA, Sing CF. A Cladistic-Analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA-sequence data.3. Cladogram Estimation. Genetics. 1992;132:619–633. [PubMed]
63. Clement M, Posada D, Crandall KA. TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000;9:1657–1659. [PubMed]
64. Camara M, Caro-Riaño H, Ravel S, Dujardin JP, Hervouet JP, et al. Genetic and morphometric evidence for isolation of a tsetse (Diptera: Glossinidae) population (Loos islands, Guinea). J Med Entomol. 2006;43:853–860. [PubMed]
65. Rohlf F J. Rotational fit (Procrustes) methods. In: Rohlf FJ, Bookstein FL, editors. Proceedings of the Michigan Morphometrics Workshop. Ann Arbor, USA: University of Michigan Museums; 1990.
66. Rohlf FJ. Marcus LF, Corti M, Loy A, Naylor GJP, Slice D, editors. Morphometric spaces, shape components and the effects of linear transformations. Advances in Morphometrics. 1996. pp. 117–129. Proc 1993 NATO-ASI on Morphometrics, Series A, Life Sciences. Plenum Publication, New York, USA.
67. Good P. Permutation tests: a practical guide to resampling methods for testing hypotheses. New York, USA: Springer; 2000.
68. Mantel N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967;27:209–220. [PubMed]
69. Jarne P, Lagoda PJL. Microsatellites, from molecules to populations and back. Trends Ecol Evol. 1996;11:424–429. [PubMed]
70. Gooding RH, Solano P, Ravel S. X chromosome mapping experiments suggest occurrence of cryptic species in the tsetse fly, Glossina palpalis palpalis (Diptera: Glossinidae). Can J Zool. 2004;82:1902–1909.
71. Dyer NA, Furtado A, Cano J, Ferreira F, Odete Afonso M, et al. 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–3282. [PubMed]
72. Vreysen MJB, Van der Vloedt AMV. The effect of intersubspecific hybrization and gamma radiation on the reproductive biology of Glossina palpalis palpalis and G. p. gambiensis. Ann Soc belge Méd Trop. 1990;70:145–158. [PubMed]
73. Morel PC, Touré SM. Glossina palpalis gambiensis Vanderplank 1949 (Diptera) dans la région des Niayes et sur la Petite Côte (République du Sénégal). Rev Elev Méd vét Pays trop. 1967;20:571–578. [PubMed]
74. Touré SM. Notes sur quelques particularités dans l'habitat de Glossina palpalis gambiensis Vanderplank 1949 (Diptera: Glossinidae) observées au Sénégal. Rev Elev Méd vét Pavs trop. 1974a;27:81–91. [PubMed]
75. Marquez JG, Vreysen MJB, Robinson AS, Bado S, Krafsur ES. Mitochondrial diversity analysis of Glossina palpalis gambiensis from Mali and Senegal. Med Vet Entomol. 2004;18:288–295. [PMC free article] [PubMed]
76. Solano P, Ravel S, de Meeus T. How can tsetse population genetics contribute to African Trypanosomiasis Control? Trends in Parasitol. 2010 in press. [PubMed]

Articles from PLoS Neglected Tropical Diseases are provided here courtesy of Public Library of Science