- Open Access
Population genetics of Glossina palpalis palpalis from central African sleeping sickness foci
Parasites & Vectorsvolume 4, Article number: 140 (2011)
Glossina palpalis palpalis (Diptera: Glossinidae) is widespread in west Africa, and is the main vector of sleeping sickness in Cameroon as well as in the Bas Congo Province of the Democratic Republic of Congo. However, little is known on the structure of its populations. We investigated G. p. palpalis population genetic structure in five sleeping sickness foci (four in Cameroon, one in Democratic Republic of Congo) using eight microsatellite DNA markers.
A strong isolation by distance explains most of the population structure observed in our sampling sites of Cameroon and DRC. The populations here are composed of panmictic subpopulations occupying fairly wide zones with a very strong isolation by distance. Effective population sizes are probably between 20 and 300 individuals and if we assume densities between 120 and 2000 individuals per km2, dispersal distance between reproducing adults and their parents extends between 60 and 300 meters.
This first investigation of population genetic structure of G. p. palpalis in Central Africa has evidenced random mating subpopulations over fairly large areas and is thus at variance with that found in West African populations of G. p. palpalis. This study brings new information on the isolation by distance at a macrogeographic scale which in turn brings useful information on how to organise regional tsetse control. Future investigations should be directed at temporal sampling to have more accurate measures of demographic parameters in order to help vector control decision.
Human African Trypanosomiasis (HAT) is a neglected tropical disease occurring in sub-Saharan Africa. After several historical cycles of epidemics followed by decreases in prevalence , WHO has recently announced the aim of elimination of HAT as a public health problem . Central Africa, in particular DRC, remains the most affected area by sleeping sickness, harbouring more than 90% of the total number of cases .
The distribution of HAT foci depends on the combined presence of the parasite, the vertebrate host, and the tsetse. The species Glossina palpalis, which is the main vector of HAT in West Africa, and which is also a vector of HAT in Central Africa, and a vector of animal trypanosomiasis in western and central Africa, is composed of two subspecies, G. p. gambiensis and G. p. palpalis. Although several studies on tsetse population genetics have been published on G. p. gambiensis, very few data are available on the population genetic structure of G. p. palpalis in central Africa [5, 6]. Investigations on population genetic structure of G. p. gambiensis have allowed observervations on genetic structuring at microgeographical scales and have allowed to measure genetic isolation between populations, which has in turn allowed control programmes to choose their control strategy (i.e. eradication or suppression) [7–9].
In the present work, we undertook a population genetic analysis of G. p. palpalis coming from different sleeping sickness foci of Cameroon and DRC using microsatellite DNA markers.
This study was undertaken in four HAT foci of Cameroon (Bafia, Bipindi, Campo and Fontem) and one HAT focus (Malanga) of the Democratic Republic of Congo. Each Cameroon focus is separated from the other by at least 100 km. No vector control activity has been undertaken in these foci so far.
Bipindi and Campo are in the South Region of Cameroon (see map in Figure 1). In both foci, previous studies identified the presence of four different tsetse species, among which G. p. palpalis was the most caught [10, 11]. Bipindi (3°2'N, 10°22'E) is an old HAT focus known since 1920 and covers several villages that are mainly located along roads . The vegetation is an equatorial forest interspersed by farmland located along the roads. This region is surrounded by hills and has a dense hydrographic network with fast running streams. It remains the most active focus of Cameroon with about 70 patients detected between 1999 and 2006. Campo (2°20'N, 9°52'E) is a hypo-endemic focus where no epidemic outbreak has been reported for several years . This focus lies along the Atlantic coast and extends along the Ntem River which constitutes the Cameroon and Equatorial Guinean border. It is an equatorial rain forest zone with a network of several rivers, swampy areas and marshes. Less than 35 cases were detected between 1999 and 2006.
The Fontem sleeping sickness focus (5°40'N, 9°55'E) is in the South-West Region of Cameroon. It has been known at a HAT focus since 1949 and has a very varied topography with hills and valleys and several fast-flowing streams. The human population, domestic animals and tsetse flies are scattered in the vegetation of the valleys and hills. Previously, the Fontem focus was among the most active foci of Cameroon . From 1998 to 2006, only 8 patients were detected among 16000 inhabitants examined (OCEAC, Unpublished data). G. p. palpalis is the only tsetse fly species known to occur in the Fontem focus [10, 15].
Bafia (4°31'N, 11°7'E) is an old HAT focus where no sleeping sickness case has been reported since 1991. It is located in the transitional zone between the forest and the savannah. G. p. palpalis is the main tsetse fly found in this locality. Nevertheless, other tsetse species such as Glossina fuscipes and Glossina fusca have been reported in this area.
Malanga (4°34'S, 14°21'E) is a HAT focus located in the Bas Congo Province of the Democratic Republic of Congo. This focus lies along the Kisantu River. The vegetation is characterized by the presence of herbaceous savannah and forest relics. G. p. palpalis is the only Glossina species found in this focus.
Entomological surveys and sampling
In Cameroon, tsetse flies were sampled in 2009 in four HAT foci, Campo in the villages (Akak, Mabiogo and Campo Beach) in March, Fontem in the villages (Bechati, Folepi and Menji) in April, Bipindi in the villages (Ebimingbang, Lambi and Memel) in July and Bafia in the village (Ombessa) in October. In each village, twelve pyramidal traps  were settled for six consecutive days in favourable tsetse fly biotopes. In the Democratic Republic of Congo, the entomological survey was carried out in August 2009 in the village Malanga (Kimpese) of the Bas Congo Province. In this village, twenty pyramidal traps  were set up for four consecutive days in tsetse fly favourable biotopes. For both Cameroon and DRC sites, tsetse flies were collected once a day. The species, sex and teneral status were identified according to routine morphological criteria .
From each G. p. palpalis individual, all legs were taken and kept into an Eppendorf tube in 95% ethanol and labelled with a code containing the trap number followed by the individual fly number. The sampling dates and the fly number were recorded in a registration book. Flies from the 12 sub-populations, except Campo Beach, Mabiogo and the two populations from DRC, came from a single trap deployed for four consecutive days. In Malanga (DRC), Campo Beach and Mabiogo, flies came from a maximum of three traps.
A total of 427 G. palpalis palpalis individuals was analysed. Flies were genotyped using 8 microsatellite loci: Gpg 55.3 , B3, B104, B110, C102 (kindly supplied by A. Robinson, Insect Pest Control Laboratory; formerly Entomology Unit, Food and Agricultural Organization of the United Nations/International Atomic Energy Agency [FAO/IAEA], Agriculture and Biotechnology Laboratories, Seibersdorf, Austria), pGp13, pGp24 , and GpCAG . From these, B104, B110, pGp13, and 55.3 are known to be located on the X chromosome [17, 18]. GpCAG and C102 have trinucleotide repeats whereas the others are dinucleotides.
DNA extraction and analysis of microsatellite loci
In each tube containing individual tsetse fly legs, 200 μl of 5% Chelex® chelating resin was added. After incubation at 56°C for 1 hour, DNA was denatured at 95°C for 30 min. The tubes were then centrifuged at 12,000 g for 2 min and frozen for later analysis.
The PCR reactions were carried out in a thermocycler (MJ Research, Cambridge, UK) as described by Ravel et al.  using 10 μl of the diluted supernatant from the extraction step in a final volume of 20 μl. After PCR amplification, allele bands were resolved on a 4300 DNA Analysis System from LI-COR (Lincoln, NE) after migration in 96-lane reloadable (3×) 6.5% denaturing polyacrylamide gels. This method allows multiplexing by the use of two infrared dyes (IRDye™), separated by 100 nm (700 and 800nm), and read by a two channel detection system that uses two separate lasers and detectors to eliminate errors due to fluorescence overlap. To determine the size of different alleles, a panel of 40 size markers was used. These markers have been previously generated for G. p. gambiensis by cloning alleles from individual tsetse flies into pGEM-T Easy Vector (Promega Corporation, Madison, WI, USA). Three clones of each allele were sequenced using the T7 primer and the Big Dye Terminator Cycle Sequencing Ready Reaction Kit (PE Applied Biosystems, Foster City, CA, USA). Sequences were analyzed on a PE Applied Biosystems 310 automatic DNA sequencer (PE Applied Biosystems) and the exact size of each cloned allele was determined. PCR products from these cloned alleles were run in the same acrylamide gel as the samples, allowing the allele size of the samples to be determined accurately. The gels were read twice by two independent readers using the LIC-OR SagaGT genotyping software.
Population parameters were assessed through Weir and Cockerham's unbiased estimators  of Wright's F-statistics  and their significance assessed through 10000 permutations with Fstat 2.9.4 (, updated from Goudet ). FIS is allele identity probability in individuals relative to allele identity between individuals from the same subsample and is thus a measure of deviation from random union of gametes within subpopulations (FIS = 0 under local panmixia). FST is allele identity between individuals from the same subsample relative to allele identity between different subsamples and is thus a measure of genetic differentiation between subpopulations (FST = 0 under random distribution of genotypes across subpopulations). Significance of FIS was assessed through randomizing alleles between individuals of the same subsamples and the statistics used was directly Weir and Cockerham's estimator as implemented in Fstat. Confidence intervals were computed with Jackknives over subsamples (for each locus) or bootstraps over loci (over all loci and subsamples) from Fstat output files .
Heterozygote deficits can be caused by Wahlund effects, null alleles, allele drop-outs or short allele dominance. Wahlund effects were first investigated through the possible role of the number of traps in a particular zone. This was investigated through a linear regression of FIS as a function of the number of traps per site. Null alleles were looked for with the software Micro-Checker v 2.2.3 . For each locus, the frequency of null alleles required to explain observed deviation from the panmictic model were computed in each subsample with van Oosterhout et al.'s method  and Brookfield's second method . These frequencies were used to compute the number of expected null homozygotes (blanks) assuming panmixia. These expected blanks numbers were summed over all subsamples for each locus and compared to the observed ones with unilateral (alternative hypothesis: there are less observed blanks than expected) exact binomial tests under R 18.104.22.168 . We also regressed these number of blanks observed at each locus and over all subsamples against mean FIS's and tested the significance of the relationship with R.
In tsetse flies, some microsatellite loci are X-linked and thus haploid in males (e.g. ). For these loci, data where coded as missing in males for heterozygote dependent analyses (FIS, null alleles, drop outs and short allele dominance), and homozygous for the allele present for other analyses (clustering, differentiation and linkage disequilibrium) following a routinely undertaken protocol [7, 8, 29, 30].
To determine the relevant unit of population structure, we used the hierarchical approach implemented in the R package HierFstat 0.04-4 . Four different hierarchical levels with four corresponding F s could be considered: The Country (Cameroon and DRC) within Total (FCT), the Village within the Country (FVC), the Site within the Village (FSV) and the Trap within the site (FTS). The significance of these different levels was tested with a G-based randomization test , the randomization unit always being lower level among the units defined by the focused level. For instance, to test the effect of the trap within site, individuals were randomized between traps of the same site. The number of randomization was set to 1000. More details on the procedures and methods implemented in HierFstat can be found elsewhere .
Linkage disequilibrium was assessed through the G-based randomization procedure per pair of locus overall subsamples, this procedure being known to be the most powerful . This was implemented in Fstat with 10000 random re-associations of alleles between loci pairs. The proportion of locus pairs that were significant was compared to the expected proportion under the null hypothesis at the 5% level of significance with a unilateral exact binomial test with alternative hypothesis "there are more than 5% significant tests in the test series" (e.g. ). In case of significance, detection of locus pairs that were responsible for this globally significant linkage was assessed using Benjamini and Hochberg's correction method . For this, the k P-values are ranked from the smallest to the largest, the highest P-value remains unchanged, the second highest P-value is multiplied by k/(k-1), the second by k(k-2) and so on until the smallest P-value that is thus multiplied by k. This procedure is thus as severe as the Bonferroni correction for the smallest P-value, and hence is very conservative (e.g. ) for a comment), but it is less stringent for the other P-values.
Sex-biased dispersal was assessed using three tests implemented in Fstat. First, Weir and Cockerham's estimate of FST, mean (mAI c ) and variance (vAI c ) of Favre et al.'s corrected assignment index AI c  were computed separately in each sex. Next, all three statistics were submitted to a permutation procedure during which the sex of each individual is randomly re-assigned in each subsample (10,000 permutations). The observed difference between male and female FST, the ratio of the largest to the smallest vAI c and the AI c -based t-statistics defined by Goudet et al.  were then compared to the resulting chance distributions. For the sex that has a higher dispersal rate, FST and mAI c are expected to be smaller and vAI c is expected to be higher than for the sex that has a lower dispersal rate (see  for more details on these tests). This choice of statistics is motivated by the work of Goudet et al.  where vAI c was shown to be the most powerful statistic when migration is low (less than 10%), while FST performs better in other circumstances. We also chose to keep mAI c because it may be more powerful in case of complex patterns of sex specific genetic structures [39, 40]. Tests were all bilateral.
Isolation by distance used the two dimensional model of Rousset . Under this model, the parameter FST/(1-FST) estimated between two subpopulations (subsamples) is a linear function of the logarithm of geographical distances Ln(G D ): FST/((1-FST) = b Ln(G D )+a and where the slope b is directly a function of demographic parameters. The product of migration rate m by subpopulation effective size N e , i.e. the number of immigrants from neighboring sites is Nm = 1/(2πb) and the product of dispersal surface σ2 by the effective density of individuals D e is D e σ2 = 1/(4πb). The significance of this regression was assessed with a Mantel test of randomization of cells of one matrix . All these isolation by distance procedures were undertaken with Genepop 4  with 1000000 iterations for the Mantel test and georeferenced coordinates in Km of traps. The software also computes bootstrap 95% confidence intervals for the slope b.
Evaluating dispersal distance between adults and their parents (σ) requires getting a proxy for effective density of adults. For this, we used three different methods for estimating mean effective population size N e . The two first methods were linkage disequilibrium based. The first is Bartley's method , from Hill  and modified by Waples  and is implemented with NeEstimator . The second is Waples and Do's method implemented by LDNe  and finally Balloux's FIS based method. This last method corresponds to heterozygote excess method from Pudovkin et al.  (see also ) corrected by Balloux . It 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 FIS in the equation N e = 1/(-2FIS)-FIS/(1+FIS)  and was only applicable in subsamples and loci with heterozygote excess, thus with very few null alleles, and probably provided overestimates in our case. These mean N e were used to estimate the effective density of tsetse flies in the different sites. A very rough estimate for the surface of a site was given by the sampling surface which was of about S = 0.15 km2 and is probably an underestimate. Hence the resulting density D e = N e /S is overestimated and dispersal distance σ = (4πbD e )-0.5 represents an underestimate.
Out of the 427 individuals constituting the 12 samples at the 8 microsatellite loci, mean genetic diversity (H s ), observed heterozygosity (H o ) and allelic richness (R s ) were greater in samples from Cameroon (H s = 0.832, H o = 0.572 and R s = 11.2 respectively) than in samples from DRC (H s = 0.723, H o = 0.488 and R s = 7.1 respectively); these differences being significant (P-values = 0.016, 0.031 and 0.013 respectively).
Agreement with genotypic proportions expected under random mating was computed within each trap. There was an important global heterozygote deficit (FIS = 0.176, P-value = 0.0001). This heterozygote deficit was highly variable across loci (Figure 2). The number of available traps in a site did not explain this heterozygote deficit (P-value = 0.91). Using the expected number of blanks (i.e. no alleles observed) computed from Micro-Checker output, we found an agreement with the hypothesis that the observed positive FIS are explained by null alleles (minimum P-value > 0.433). At least 73% of the variance of FIS is explained by the number of blank genotypes found across loci (P-value = 0.007) (Figure 3).
HierFstat analysis gave a negative (hence non-significant) FTS (no trap effect) and a very small and not significant effect for sites within villages FSV~0 (P-value = 0.823). The only significant effect was found for villages with a FVC = 0.028 (P-value = 0.001), the effect of countries displaying no additional significant effect (FCT = 0.053, P-value = 0.092). The relevant unit that we kept for the following analyses was thus the village.
Among the 28 possible pairs of loci that could be tested, four displayed significant linkage (14%). This is more than the 5% expected under the null hypothesis (exact binomial test, P-value = 0.0491). None of these tests remained significant after Benjamini and Hochberg's correction. This marginally significant linkage at the genome wide scale is thus probably coming from demographic causes (small N e ).
Sex biased dispersal
Assignment based parameters were in line with a male biased dispersal but only mAI c provided a significant test (P-value = 0.013) (Table 1). Thus the signature of sex-biased dispersal is weak.
Isolation by distance, population density and dispersal
Figure 4 shows that there is a highly significant isolation by distance with a slope b = 0.0099 with 95% confidence interval [0.006, 0.017], a number of migrants of Nm = 16 in 95% CI (9, 25) and a product D e σ2 = 8 in 95% CI (5, 13).
Estimates of effective population sizes varied greatly according to the methods used (Table 2). In particular Waples and Do's method gave inconsistent values, and the reality must lie between values given by Bartley's and Balloux's methods. Hence, effective population sizes are probably between 20 and 300 individuals leading to an estimated migration rate between m~0.05 and m~0.8. Then, assuming densities stretching over between 120 and 2000 individuals per km2, dispersal distance between reproducing adults and their parents can be estimated between σ~60 and σ~300 meters (Table 2). The fact that Balloux's method provided several defined values of N e (in fact 17 FIS values were negative and provided an N e >0) for four loci in various subsamples (nine) is also indicative that null alleles are indeed responsible for the global heterozygote deficit observed and not a Wahlund effect.
Although there have been studies dealing with palpalis group tsetse population structure [4, 8, 52], this work represents one of the first to focus on the sleeping sickness vector G. p. palpalis in Central Africa [see also ].
Population structuring was found at the geographical scale of the village, but not at the scale of the traps and above all no hidden substructure (Wahlund effect) was evidenced here, which is at variance with that reported for G. p. palpalis in Ivory Coast . In addition, in the present work, null alleles probably explain all of the FIS that was observed. As a consequence, this high frequency of null alleles may have hampered our estimates of population structure. The effective population sizes found here (between 20 and 300 as a mean for each village) are of the same order of magnitude as found for G. p. gambiensis in Guinea . It will be of interest to reassess these population sizes and densities in the future using temporal sampling to get more accuracy.
A slight male biased dispersal was observed, which does not correspond to typical Mark-Release-Recapture studies, these latter showing in general that females disperse more than males . It is however noteworthy that using microsatellite markers, very recent studies on other tsetse species report the same trend, i.e. a male sex-biased dispersal [54, 55]. Either a behaviour such as philopatry (where the females would come back to the same larviposition sites), or a sex specific local adaptation rendering immigrant females very unlikely to survive locally, may explain this.
A strong isolation by distance explains most of the population structure observed in our sampling sites of Cameroon and DRC. The populations here are composed of random mating subpopulations occupying fairly wide zones with a very strong isolation by distance that makes the probability of an allele to cross from one focus to the other (e.g: from Bipindi to Campo) very unlikely. This was reinforced by the observation of differences in number of alleles and heterozygosity between populations from Cameroon and from DRC. It has been recently reported that G. p. palpalis from Equatorial Guinea may constitute a different gene pool compared with G. p. palpalis from other countries . Here, the differences in population structure as compared to Ivory-Coast counterparts might reflect again such a taxonomic heterogeneity.
Future studies on this important vector of sleeping sickness in Central Africa ought to focus on temporal sampling to better assess population densities to help targeting efficient control strategies.
Courtin F, Jamonneau V, Duvallet G, Garcia A, Coulibaly B, Doumenge JP, Cuny G, Solano P: Sleeping sickness in West Africa (1906-2006): changes in spatial repartition and lessons from the past. Trop Med Int Health. 2008, 13: 334-344. 10.1111/j.1365-3156.2008.02007.x.
Simarro P, 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 (2): e1007-10.1371/journal.pntd.0001007.
Simarro P, Cecchi G, Paone M, Franco JR, Diarra A, Ruiz Postigo JR, Fèvre EM, Courtin F, Mattioli RC, Jannin JG: The Atlas of human African trypanosomiasis: a contribution to global mapping of neglected tropical diseases. Int J Health Geogr. 2010, 9: 57-10.1186/1476-072X-9-57.
Solano P, Ravel S, de Meeûs T: How can tsetse population genetics contribute to African Trypanosomosis control?. Trends Parasitol. 2010, 26: 255-263. 10.1016/j.pt.2010.02.006.
Dyer NA, Lawton SP, Ravel S, Choi KS, Lehane MJ, Robinson AS, Okedi LA, Hall M, Solano P, Donnelly MJ: 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. 10.1016/j.ympev.2008.07.011.
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-3282. 10.1111/j.1365-294X.2009.04265.x.
Bouyer J, Ravel S, Guerrini L, Dujardin JP, Sidibé I, Vreysen MJ, Solano P, de Meeûs T: Population structure of Glossina palpalis gambiensis (Diptera: Glossinidae) between river basins in Burkina-Faso: consequences for area-wide integrated pest management. Infect Genet Evol. 2010, 10: 321-328. 10.1016/j.meegid.2009.12.009.
Koné N, de Meeûs T, Bouyer J, Ravel S, Guerrini L, Guerrini L, N'Goran EK, Vial L: Population structuring of Glossina tachinoides (Diptera: Glossinidae) according to landscape fragmentation in the Mouhoun river basin, Burkina Faso. Med Vet Entomol. 2010, 24: 162-168. 10.1111/j.1365-2915.2010.00857.x.
Kagbadouno M, Camara M, Bouyer J, Courtin F, Morifaso O, Solano P: Progress towards the eradication of tsetse from Loos islands, Guinea. Parasites and Vectors. 2011, 4: 18-10.1186/1756-3305-4-18.
Simo G, Njiokou F, Mbida Mbida JA, Njitchouang GR, Herder S, Asonganyi T, Cuny G: Tsetse fly host preference from sleeping sickness foci in Cameroon: epidemiological implications. Infect Genet Evol. 2008, 8: 34-39. 10.1016/j.meegid.2007.09.005.
Farikou Njiokou F, Simo G, Asonganyi T, Cuny G, Geiger A: Tsetse fly blood meal modification and trypanosome identification in two sleeping sickness foci in the forest of southern Cameroon. Acta Trop. 116: 81-88.
Grébaut P, Mbida Mbida JA, Antonio Kondjio C, Njiokou F, Penchenier L, Laveissière C: Spatial and Temporal Patterns of Human African Trypanosomosis (HAT) Transmission Risk in Bipindi Focus, in Forest Zone of Southern Cameroon. Vector-Borne Zoonot. 2004, 4: 230-238.
Penchenier L, Gébaut P, Ebo'o Eyenga V, Bodo JM, Njiokou F, Binzouli JJ, Simarro P, Soula G, Herder S: Le foyer de la Trypanosomiase Humaine Africaine de Campo (Cameroun) en 1998. Aspects épidémiologiques, état de l'endémie et comparaison des CATT 1.3 et CATT Latex dans le dépistage de masse. Bull Soc Pathol Exot. 1999, 92: 185-190.
Asonganyi T, Ade SS: Review of sleeping sickness in Cameroon. J Cam Med. 1992, 3: 30-37.
Morlais I, Grébaut P, Bodo J, Djoha S, Cuny G, Herder S: Detection and identification of trypanosomes by polymerase chain reaction in wild tsetse flies in Cameroon. Acta Trop. 1998, 70: 109-117. 10.1016/S0001-706X(98)00014-X.
Gouteux J, Lancien J: Le piège pyramidal à tsé-tsé (Diptera: Glossinidae) pour la capture et la lutte. Essais comparatifs et description de nouveaux systèmes de capture. Tropenmed Parasit. 1986, 37: 61-66.
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. 10.1016/S0001-706X(97)00663-3.
Luna C, Bonizzoni M, Cheng Q, Robinson AS, Aksoy S, Zheng L: Microsatellite polymorphism in tsetse flies (Diptera: Glossinidae). J Med Entomol. 2001, 38: 376-1. 10.1603/0022-2585-38.3.376.
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. 10.1046/j.1471-8278.2001.00087.x.
Ravel S, De Meeus T, Dujardin JP, Zézé DG, Gooding RH, Sane B, Dusfour I, Cuny G, Solano P: Different genetic groups occur within Glossina palpalis palpalis in the sleeping sickness focus of Bonon, Côte d'Ivoire. Infect Genet Evol. 2007, 7: 116-125. 10.1016/j.meegid.2006.07.002.
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.
Wright S: The interpretation of population structure by F-statistics with special regard to system of mating. Evolution. 1965, 19: 395-420. 10.2307/2406450.
Goudet J: Fstat (ver. 2.9.4): A Program to Estimate and Test Population Genetics Parameters 2003. Updated from Goudet , [http://www.unil.ch/izea/softwares/fstat.html/]
Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.
De Meeûs T, McCoy KD, Prugnolle F, Chevillon C, Durand P, Hurtrez-Boussès S, Renaud F: Population genetics and molecular epidemiology or how to "débusquer la bête". Infect Genet Evol. 2007, 7: 308-332. 10.1016/j.meegid.2006.07.003.
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. 10.1111/j.1471-8286.2004.00684.x.
Brookfield JFY: A simple new method for estimating null allele frequency from heterozygote deficiency. Mol Ecol. 1996, 5: 453-455.
R-Development-core-team: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. 2011, ISBN 3-900051-07-0, [http://www.R-project.org]
Solano P, Kaba D, Ravel S, Dyer NA, Sall B, Vreysen MJ, Seck MT, Darbyshir H, Gardes L, Donnelly MJ, De Meeûs T, Bouyer J: Population genetics as a tool to select tsetse control strategies: suppression or eradication of Glossina palpalis gambiensis in the Niayes of Senegal. PLoS Negl Trop Dis. 2010, 4 (5): e692-10.1371/journal.pntd.0000692.
Solano P, Ravel S, Bouyer J, Camara M, Kagbadouno MS, Dyer N, Gardes L, Herault D, Donnelly MJ, De Meeûs T: The population structure of Glossina palpalis gambiensis from island and continental locations in coastal Guinea. PLoS Negl Trop Dis. 2009, 3: e392-10.1371/journal.pntd.0000392.
Goudet J: HIERFSTAT, a package for R to compute and test hierarchical F-statistics. Mol Ecol Notes. 2005, 5: 184-186. 10.1111/j.1471-8286.2004.00828.x.
Goudet J, Raymond M, De Meeûs T, Rousset F: Testing differentiation in diploid populations. Genetics. 1996, 144: 1933-1940.
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. 10.1016/j.meegid.2007.07.005.
De Meeûs T, Guégan JF, Teriokhin AT: MultiTest V.1.2, a program to binomially combine independent tests and performance comparison with other related methods on proportional data. BMC Bioinformatics. 2009, 10: 443-10.1186/1471-2105-10-443.
Benjamini Y, Hochberg Y: On the adaptive control of the false discovery fate in multiple testing with independent statistics. Journal of Educational and Behavioral Statistics. 2000, 25: 60-83.
Favre L, Balloux F, Goudet J, Perrin N: Female-biased dispersal in the monogamous mammal Crocidura russula: evidence from field data and microsatellite patterns. Proc R Soc London B. 1997, 264: 127-132. 10.1098/rspb.1997.0019.
Goudet J, Perrin N, Waser P: Tests for sex-biased dispersal using bi-parentally inherited genetic markers. Mol Ecol. 2002, 11: 1103-1114. 10.1046/j.1365-294X.2002.01496.x.
Prugnolle F, De Meeûs T: Inferring sex-biased dispersal from population genetic tools: a review. Heredity. 2002, 88: 161-165. 10.1038/sj.hdy.6800060.
De Meeûs T, Beati L, Delaye C, Aeschlimann A, Renaud F: Sex-biased genetic structure in the vector of Lyme disease, Ixodes ricinus. Evolution. 2002, 56: 1802-1807.
Kempf F, McCoy KD, De Meeus T: Wahlund effects and sex-biased dispersal in Ixodes ricinus, the European vector of Lyme borreliosis: new tools for old data. Infect Genet Evol. 2010, 10: 989-997. 10.1016/j.meegid.2010.06.003.
Rousset F: Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997, 145: 1219-1228.
Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27: 209-220.
Rousset F: GENEPOP ' 007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.
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. 10.1046/j.1523-1739.1992.06030365.x.
Hill WG: Estimation of effective population size from data on linkage disequilibrium. Genet Res. 1981, 38: 209-216. 10.1017/S0016672300020553.
Waples RS: Genetic methods for estimating the effective size of cetacean populations. Genetic ecology of whales and dolphins. Edited by: Hoelzel AR. 1991, International Whaling Commission, 279-300. Special Issue No. 13
Peel D, Ovenden JR, Peel SL: NeEstimator Version 1.3: software for estimating effective population size, version 1.32004. Queensland Government, Department of Primary Industries and Fisheries.
Waples RS, Do C: LDNE: a program for estimating effective population size from data on linkage disequilibrium. Mol Ecol Resour. 2008, 8: 753-756. 10.1111/j.1755-0998.2007.02061.x.
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.
Luikart G, Cornuet JM: Estimating the effective number of breeders from heterozygote excess in progeny. Genetics. 1999, 151: 1211-1216.
Balloux F: Heterozygote excess in small populations and the heterozygote-excess effective population size. Evolution. 2004, 58: 1891-1900.
Echodu R, Beadell JS, Okedi LM, Hyseni C, Aksoy S, Caccone A: Temporal stability of Glossina fuscipes fuscipes populations in Uganda. Parasites and Vectors. 2011, 4: 19-10.1186/1756-3305-4-19.
Cuisance D, Février J, Dejardin J, Filledier J: Dispersion linéaire de Glossina palpalis gambiensis et G. tachinoides dans une galerie forestière en zone soudano-guinéenne (Burkina Faso). Rev Elev Méd vét Pays Trop. 1985, 38: 153-172.
Koné N, Bouyer J, Ravel S, Vreysen MJB, Domagni KJ, Causse S, Solano P, De Meeûs T: Contrasting Population Structures of Two Vectors of African Trypanosomoses in Burkina Faso: Consequences for Control. PLoS Negl Trop Dis. 2011,
Ouma JO, Beadell JS, Hyseny C, Okedi LO, Krafsur ES, Aksoy SA, Caccone A: Genetic diversity and population structure of Glossina pallidipes in Uganda and western Kenya. Parasites and Vectors. 2011,
European Foundation Initiative for Neglected Tropical Disease (EFINTD), AIRES-SUD project n°7003, IRD, the University of Yaoundé I, the University of Dschang, and B. Noose.
The authors declare that they have no competing interests.
TDM participated to the tsetse fly sampling, the genotyping step and the manuscript drafting. GS participated to the conception, the design of the study, the tsetse fly sampling and the drifting of the manuscript. SR was involved on the genotyping step and the drifting of the manuscript. TDM performed the statistical analysis and helped to draft the manuscript. SC participated to the genotyping step and the interpretation of data. PS participated to the conception of the study, the data analysis and the drafting of the manuscript. PL participated to the conception and the tsetse fly sampling. TA participated to the conception and the design of the study. FN was involved on the conception of the study, the tsetse fly sampling and the drafting of the manuscript. All the authors read and approved the final manuscript.