Population genetics of Glossina palpalis palpalis from central African sleeping sickness foci
© Tanekou et al; licensee BioMed Central Ltd. 2011
Received: 31 May 2011
Accepted: 18 July 2011
Published: 18 July 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.
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 188.8.131.52 . 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).
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
Results of sex-biased dispersal tests on G. palpalis palpalis from Cameroon and DRC
Isolation by distance, population density and dispersal
Estimation of effective population sizes with three different methods.
LD Waples and Do
1442 [0, 4428]
321 [102, 539]
18 [7, 29]
D e (individuals/km2)
0.01 [0.006, 0.02]
0.05 [0.03, 0.08]
0.89 [0.51, 1]
σ (in meters)
29 [22, 36]
61 [46, 77]
258 [196, 322]
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.
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.
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle 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-3282. 10.1111/j.1365-294X.2009.04265.x.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Asonganyi T, Ade SS: Review of sleeping sickness in Cameroon. J Cam Med. 1992, 3: 30-37.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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/]Google Scholar
- Goudet J: FSTAT (Version 1.2): A computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Brookfield JFY: A simple new method for estimating null allele frequency from heterozygote deficiency. Mol Ecol. 1996, 5: 453-455.View ArticlePubMedGoogle Scholar
- 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]Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Goudet J, Raymond M, De Meeûs T, Rousset F: Testing differentiation in diploid populations. Genetics. 1996, 144: 1933-1940.PubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Rousset F: Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997, 145: 1219-1228.PubMed CentralPubMedGoogle Scholar
- Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27: 209-220.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Hill WG: Estimation of effective population size from data on linkage disequilibrium. Genet Res. 1981, 38: 209-216. 10.1017/S0016672300020553.View ArticleGoogle Scholar
- 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. 13Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- Luikart G, Cornuet JM: Estimating the effective number of breeders from heterozygote excess in progeny. Genetics. 1999, 151: 1211-1216.PubMed CentralPubMedGoogle Scholar
- Balloux F: Heterozygote excess in small populations and the heterozygote-excess effective population size. Evolution. 2004, 58: 1891-1900.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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,Google Scholar
- 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,Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.