- Short report
- Open Access
Genetic diversity of Glossina fuscipes fuscipes along the shores of Lake Victoria in Tanzania and Kenya: implications for management
Parasites & Vectors volume 10, Article number: 268 (2017)
Tsetse flies (Diptera: Glossinidae) are sole vectors for trypanosomiasis, which affect human health and livestock productivity in Africa. Little is known about the genetic diversity of Glossina fuscipes fuscipes, which is an important species in Tanzania and Kenya. The main objective of the study was to provide baseline data to determine the genetic variability and divergence of G. f. fuscipes in the Lake Victoria basin of Tanzania and Kenya in order to guide future vector control efforts in the region.
Two hundred and seventy five G. f. fuscipes from 8 sites along the shores of Lake Victoria were screened for genetic polymorphisms at 19 microsatellite loci. Samples were collected from two sites in Kenya and six sites in Tanzania. Four of the Tanzanian sites were located in the Rorya district, on the eastern shores of Lake Victoria, while the other two sites were from Ukerewe and Bukoba districts from the southern and western Lake Victoria shores, respectively. Four genetically distinct allopatric clusters were revealed by microsatellite analysis, which sorted the sampling sites according to geography, with sites separated by as little as ~65 km belonging to distinct genetic clusters, while samples located within ~35 km from each other group in the same cluster.
Our results suggest that there is ongoing genetic admixture within sampling sites located ~35 km from each other, while sites located ~65 km apart are genetically isolated from each other. Similar patterns emerged from a parallel study on G. f. fuscipes analyzed from the Lake Victoria Uganda shores. From a control perspective these results suggest that for sites within the same genetic cluster, control efforts should be carried out in a coordinated fashion in order to avoid re-invasions. Future work should focus on better quantifying the extent and spatial patterns of the observed genetic discontinuities of the G. f. fuscipes populations along the Tanzanian shores. This will aid in their control by providing guidelines on the geographical extent of the area to be treated at the same time.
Tsetse flies (Diptera: Glossina) remain to be insects of economic and medical importance in sub-Saharan Africa. They transmit pathogenic trypanosomes that cause sleeping sickness to humans and nagana to livestock [1, 2]. The occurrence of the disease in humans and livestock has greatly limited the development of agriculture and human health in the region . Estimates by FAO show, Africa loses over 3 million cattle and other domestic livestock due to trypanosomiasis every year . The annual losses, in terms of reduced meat and milk production and in terms of the costs related to treatment and controlling the disease, have been estimated at US $1.2 billion . No vaccine has been developed for the disease to date both for humans and cattle due to the ability of trypanosome parasites to change their surface proteins by antigenic variation [6, 7]. Drugs which are used to treat cattle have been used for long time; as a result drug resistance is increasing rapidly which in the long run will seriously affect the use of these drugs for animal trypanosomiasis control. Similarly, drugs which are used to treat humans are toxic, expensive and difficult to administer at village settings and also have bad side effects [1, 8]. Therefore, tsetse control remains the most effective method of reducing trypanosome infections in animals and human in sub-Saharan countries [6,7,8,9].
Glossina fuscipes fuscipes is one of the most important tsetse species in the Palpalis group, subgenus Nermohina Robineau-Desvoidy . In Tanzania the species is widely distributed along the shores of Lake Victoria supported by vegetation growing close to the water . The species is found at the eastern margin in Uganda which extends further east along the shores of Lake Victoria in Western Kenya. The fly is also found in southern Sudan, Chad, the Central African Republic, the Democratic Republic of Congo (DRC) and Angola [12, 13]. Various methods have been employed in tsetse control in different areas in Africa, but the success of these methods varied. Some of these interventions used in the past included bush clearing (destroying vector habitat), elimination of wild animals (reservoir host of the parasites), insecticide ground spraying, live bait technology and the use of baited traps and targets. Despite the fact that tsetse fly densities were temporarily reduced, some of the methods used in the past are not used today because of environment reason and also they are against animal conservation [2, 5, 6]. Sterile Insect Technique (SIT) and Sequential Aerial Technique (SAT) have become promising methods which have been successful used to eradicate tsetse flies in some parts of Africa. SIT was successfully used to eradicate Glossina austeni in Unguja Island in Tanzania [2, 14] and SAT has been successfully applied in the Okavango Delta in northern Botswana in the eradication of Glossina morsitans centralis . The success of these control interventions depend much on the biogeographical limits of the target tsetse species in such a way that maximum benefit is gained from the natural barriers to reinvasion of the previous controlled areas . One factor that can improve the application of various control methods is knowledge on tsetse population genetics in addition to ecological information . Studying the genetic differentiation between different populations of G. f. fuscipes will lead to novel insights into the relationships between genetically distinct populations, including geographical distribution, hybridization and migration patterns. This knowledge can then be used to inform ongoing or planned vector control programs across a target geographical region to identify the most suitable areas to target control to avoid re-infestation of cleared habitats. Studies on the population structure of G. f. fuscipes have been carried out in some parts of East Africa, particularly along Lake Victoria shores in Uganda [8, 15,16,17,18], and reported high genetic structuring of the species in the area. All authors confirmed the presence of two distinct lineages of G. f. fuscipes in Uganda using mtDNA and microsatellites markers except author 18. The isolation of northern and southern G. f. fuscipes populations which were thought to occur as a result of fragmentation during extreme drought in East Africa warrants these populations to be treated separately during eradication. Likewise, population study of another tsetse riverine species G. palpalis gambiensis in Senegal also indicated the species was sufficiently isolated and currently eradication is underway [9, 19]. However, population structure of G. p. palpalis in Burkina Faso and G. p. gambiensis in Equatorial Guinea indicated high gene flow compared to other riverine tsetse species [20, 21]. Similarly, population structure study on the same species in Cameroon indicated the species had heterozygote deficit suggesting suppression to be the best option to control the tsetse species in the area .
Little is known on the population structure and gene flow of G. f. fuscipes in the Lake Victoria basin in Kenya and Tanzania. The information will be useful in the planning of effective regional control of the tsetse species in the Lake Victoria basin (the Lake Victoria basin includes parts of Uganda, Tanzania and Kenya) under Pan African Tsetse and Trypanosomiasis Eradication Campaign (PATTEC) initiative . Lake Victoria basin is one among important regions which have been identified for tsetse eradication by PATTEC. This report describes patterns of genetic differentiation among G. f. fuscipes sampled in Tanzania and Kenya and relates them to the recommendations for guiding future vector control efforts in the region.
Study sites and data collection
Glossina fuscipes fuscipes were sampled at 8 localities around Lake Victoria, with 2 sites in Kenya (KIS and MAN) and 6 sites in Tanzania. These sites include 4 locations on the border of Tanzania and Kenya in the Rorya district (MAS, TOB, RAS and KIR), one in the Bukoba district at the border between Uganda and Tanzania (BUK), and one approximately in between these two areas in the Ukerewe district (UKE; Fig. 1). The average distance among the four sampling sites in the Rorya district ranged from 2 to 10 km, the two sampling sites from Kenya (KIS and MAN) are 35 km apart, ~65 km from the Rorya district samples, and ~200 km away from the other two Tanzanian samples from the southern and western shores of Lake Victoria (BUK and UKE), which are also similarly distant from each other (187 km). Overall, 275 individual tsetse flies were collected using biconical and pyramid traps and preserved in cryotubes containing 90% ethanol. Flies from the Tanzanian sites were collected in 2011 while the ones from the two Kenyan sites were collected in 2009. All samples were collected during the dry season. All study sites experience a bimodal rainfall pattern, short rains fall between October and December and long rains from March to May. The average annual rainfall in the basin is estimated to be 1,015 mm . The vegetation of the study sites was bushy vegetation and was the same in all study sites (Fig. 2).
DNA extraction and PCR analysis
DNA extraction from tsetse legs was conducted using a Qiagen Micro Amp extraction kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. PCR and genotyping were conducted for microsatellite markers at 19 previously identified loci that are distributed across the genome . PCR conditions used were in accordance with previously described methods . Microsatellite peaks were scored for each individual sample using GeneMarker software . Table 1 reports the number of flies genetically analyzed for each sampling site.
The program Genepop 4.2  was used to test for deviations from Hardy-Weinberg equilibrium, estimate allelic richness (AR), calculate expected and observed heterozygosities (He and Ho), and evaluate levels of inbreeding (Fis). Pairwise FST values between all sampling sites were calculated in Arlequin v 3.5 . Significance was tested using 10,000 permutations. Genetic differentiation among tsetse sampling sites was analyzed using a model-based Bayesian clustering method implemented in Structure 2.3.3 . The optimal number of clusters in the data was calculated using the ad hoc statistic ΔK in STRUCUTURE HARVESTER . The STRUCTURE analysis was run for k = 1 to 8, with 10 replicates of each k value, for a million generations and a ‘burn-in’ period lasting the first 200,000 generations. Afterwards, the results from the STRUCTURE analysis were summarized across the 10 replicates using the Greedy method of the software CLUMPP , and distruct package was used to plot the results. We also used discriminant analysis of principal components (DAPC)  conducted with the R adegenet package  to determine genetic structure. In addition, we used the Bayesian Information Criterion (BIC) to determine the most likely number of clusters in our data .
To investigate patterns of migration and identify any probable migrants within our dataset, we used GeneClass2 . GeneClass2 computes the probability of an individual’s observed multilocus genotype belonging to a given population. As such, the program can detect whether an individual is a resident of the sampled population or a first generation migrant. Effective population sizes (Ne) for the clustered populations were estimated under the linkage disequilibrium model with random mating included in the program NeEstimator v2 . All alleles were taken into account regardless of frequency. Jackknife on loci calculations was implemented to generate 95% confidence bounds.
Results and discussion
In the present study, a total of 275 tsetse flies were genotyped at 19 microsatellite loci. According to inbreeding coefficient values for each population, no loci deviated significantly from Hardy-Weinberg equilibrium (Table 1, Additional file 1: Table S1). The number of alleles, a measure of genetic diversity, varied from 4.1 in Rasi Nyabero (RAS) to 5.6 in Manga (MAN). The number of alleles per locus varied from 2 in D101 to 15 in GmA06 with an average of 9.4. Mean allelic richness (AR) ranged from 2.89 to 4.94 across the sites. No significant differences were observed between observed and expected heterozygosities (HE ranged from 5.07 to 19.85 per site, and HO ranged from 3.50 to 20.53 per site). These values are similar to the ones reported for G. f. fuscipes populations around Lake Victoria shores in Uganda [6, 19] and suggest that G. f. fuscipes flies in this region harbor a wealth of genetic diversity and thus have the long term potential to adapt to ecological changes in the area.
FST values revealed significant genetic differentiation among sampling sites over 60–65 km apart, with notable exceptions being the significant FST values found between two close sampling sites in the Rorya district (KIR and RAS, < 6 km apart, Additional file 1: Table S2) and the lack of significance for FST values between samples from KIS and either RAS or MAN, located ~ 68–71 km from KIS (Additional file 1: Table S2). When the population structure was investigated using Bayesian clustering and multivariate analyses, the data was most consistent with 4 distinct population clusters (Fig. 3; Additional file 1: Figure S1, and Additional file 1: Figure S2), with FST values between clusters being all statistically significant and ranging from 0.231 to 0.341 (Additional file 1: Tables S3 and S4). Clusters 1 and 2 group flies from BUK and UK in Bukoba and Ukerewe districts, respectively. Cluster 3 comprises flies from all the Rorya district sampling sites (Rasi Nyabero, RAS; Kirongwe, KIR; Masonga, MAS; and River Tobwe, TOB). Cluster 4 includes the two Kenyan locations (Kisasi, KIS and Manga, MAN). Individual flies collected in sampling sites within these clusters assign mostly to their respective clusters with only a few individuals showing signs of genetic admixture with another cluster. The high genetic differentiation (FST = 0.341; FST = 0.308 and FST = 0.286) recorded in this study is consistent with earlier reports on the same species in other areas of East Africa and is likely a consequence of genetic drift. In the riverine and lacustrine tsetse species, such as G. f. fuscipes, wet season dispersion from dry season refugia is thought to establish demes in which genetic drift leads to differentiation [8, 16, 18]. Similar results have been reported also for other riverine tsetse species in West Africa which include G. tachinoides in Ghana [35, 36], G. p. gambiensis in Niayes-Senegal [9, 37] and G. p. palpalis in Cote d’Ivore . However different results of similar species (G. p. palpalis) in Cameroon indicated the species formed a large panmixia population suggesting suppression to be the best option for the vector control . In contrast to other riverine tsetse species, G. swynnertoni, a savannah species from northern Tanzania, have been reported to have the high level of gene flow . Geographical distance and ecological differences are factors that are involved in these genetic differences of flies .
The presence of admixed genotypes within clusters 2 and 3 (Fig. 3) suggests ongoing or recent gene flow among sampling locales. This is confirmed by the GeneClass2 analysis aimed to detect first generation migrants, which identified nine likely migrants: three between clusters 1 and 3, two between clusters 2 and 3, and four between clusters 3 and 4 (Fig. 4). However, given the lack of data for intermediate locations, it is not possible to evaluate if this mixing is due to gene flow from intermediate locations or long-range dispersal. These patterns of genetic discontinuities mirror the ones obtained for the same species along the Uganda coast of Lake Victoria, where among similarly spaced samples genetic discontinuities were identified using microsatellite loci .
Although estimates of effective population size (Ne) were not uniform among sampling sites, and given the confidence intervals for clusters 3 and 4 (Fig. 4), the Ne estimates may not be very reliable. The high genetic diversity (FST = 0.341; FST = 0.308 and FST = 0.286) recovered at some sampling sites (Table 1) points to large resident populations. Recent reports indicated high and varied G. f. fuscipes densities in the sampling sites analyzed in this study . These data indirectly seem to support the existence of local large pockets of tsetse flies at these sampling sites. Further screening of tsetse flies at these and neighboring sites, including temporal collections to capture seasonal variations in fly population sizes and densities, would shed additional light on local tsetse population dynamics, which is of critical relevance for monitoring and control.
From a control perspective these results suggest that, given the genetic homogeneity among sampling sites from the Rorya district, these sites should be controlled and monitored together to guard against reinvasion from neighboring untreated areas. A similar strategy should be adopted for the region including the two Kenyan sampling sites. However, given that the microsatellite data show that Kenya and the Rorya district samples are genetically distinct even though they are only ~65–80 km apart, control and monitoring effort may not necessarily need to include large geographical areas. Future studies should include a more uniform geographical representation of tsetse infested areas than in the current study to quantify the extent and spatial pattern of the genetic discontinuities found in this study among sampling sites located at varying geographical distances.
Bayesian Information Criterion
Primer name derived from loci (D101) for assessing genetic diversity
Discriminant Analysis of Principal components
- FIS :
Primer name derived from the loci of mutant lines of interest for assessing genetic diversity
Effective population size
Echodu R, Sistrom M, Hyseni C, Enyaru J, Okedi L, Aksoy S, Caccone A. Genetically distinct Glossina fuscipes fuscipes populations in the Lake Kyoga Region of Uganda and its relevance for human African trypanosomiasis. Bio Med Res Int. 2013;614721:2–4.
Leak SGA. Tsetse biology and ecology. Wallingford: CABI Publishing; 1999. p. 189–90.
Msangi AR, Kiwia NE, Mramba F, Malele II, Kitwika WA, Byamungu MB, et al. After Zanzibar using SIT Mafia is next. Proceedings ISCTRC. Mombasa 1999;30.
FAO. Impacts of trypanosomiasis on African agriculture. PAAT Technic Sci Series. 2000;2:4–8.
Kuzoe FAS, Schofield CJ. Strategic review of traps and targets for tsetse and African trypanosomiasis control. UNICEF/UNDP/World bank /WHO Special Programme for Research and Training in Tropical Diseases (TDR). 2004;9–15.
Aksoy S. Control of tsetse flies and trypanosomes using molecular genetics. J Vet Parasitol. 2003;115:125–45.
Mélachio TTT, Njiokou F, Ravel S, Simo G, Solano P, Meeûs T. Effect of sampling methods, effective population size and migration rate estimation in Glossina palpalis palpalis from Cameroon. Infect Genet Evol. 2015;33:150–7.
Beadell JS, Hyseni C, Abilla PP, Azabo R, Enyaru JCK, Ouma JO, et al. Phylogeography and population structure of Glossina fuscipes fuscipes in Uganda: Implications for control of tsetse. PLoS Negl Trop Dis. 2010;4:636.
Solano P, Delarocque S, Cuisance D, Geoffroy B, Demeeus T, Cuny G, Duvallet G. Intraspecific variability in natural populations of Glossina palpalis gambiensis and morphometric analyses. Med Vet Entomol. 2010;13:401–7.
Gooding RH, Krasfur ES. Tsetse genetics. Contributions to biology, systematics, and control of tsetse flies. Annu Rev Entomol. 2005;50:101123.
Manangwa O, Ouma JO, Msangi A, Mramba F, Gamba N. Distribution and population size of G. f. fuscipes (tsetse flies) along the Lake Victoria, for trypanosomiasis management in Tanzania. Livestock Res Rural Dev. 2015;27:1–2.
Albert M, Wardrop NA, Atkinson PM, Torr SJ, Welburn SC. Tsetse fly (G. f. fuscipes). distribution in the Lake Victoria Basin of Uganda. PLoS Negl Trop Diseases. 2015;9(4):2–3.
Rogers DJ, Robinson TP. Tsetse distribution. In: Maudin I, Holmes PH, Miles MA, editors. The Trypanosomes. Wallinford: CABI publishing; 2004. p. 139–79.
Vreysen MJB, Saleh KM, Ali MY, Abdulla AM, Zhu ZR, Juma KG, et al. Glossina austeni (Diptera: Glossinidae) eradicated on the island of Unguja, Zanzibar, using the sterile insect technique. J Appl Entomol. 2000;93:123–35.
Echodu R, Beadell JS, Okedi LM, Hyseni C, Aksoy S, Caccone A. Temporal stability of Glossina fuscipes fuscipes populations in Uganda. Parasit Vectors. 2011;4:2–3.
Hyseni C, Kato AB, Okedi LM, Masembe C, Ouma JO, Aksoy S, et al. The population structure of Glossina fuscipes fuscipes in the Lake Victoria basin in Uganda: implications for vector control. Parasit Vectors. 2012;5:222.
Kato BA, Hyseni C, Okedi LM, Ouma JO, Aksoy S, Caccone A, Masembe C. Mitochondrial DNA sequence divergence and diversity of Glossina fuscipes fuscipes in the Lake Victoria basin of Uganda: implications for control. Parasit Vectors. 2015;8:385.
Krafsur ES, Marquez JG, Ouma JO. Structure of some East African Glossina fuscipes fuscipes populations. Med Vet Entomol. 2008;22:222–7.
Kagbadouno M, Camara M, Bouyer J, et al. Tsetse elimination: its interest and feasibility in the historical sleeping sickness focus of loos islands, Guinea. Parasite. 2009;16(1):29–35.
Bouyer J, Sibert A, Desquesnes M, Cuisance D, De La Rocque S. A diffusion model for Glossina palpalis gambiensis in Burkina Faso. In: Vreysen MJB, Robinson AS, Hendrichs J, editors. Area-wide control of insect pests. From research to field implementation. Dordrecht: Springer; 2007. p. 221–8.
Dyer NA, Furtado A, Cano J, 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(15):3268–82.
Schofield CJ, Kabayo JP. Trypanosomiasis vector control in Africa and Latin America. Parasit Vectors. 2008;1(24):1756–3305.
Makalle AMP, Obando J, Bamutaze Y. Effects of land use practices on livelihoods in the transboundary sub-catchments of the Lake Victoria Basin. Afr J Environ Sci Technol. 2008;2(10):309.
Soft genetics. GeneMarker genotyping software. Release 1.85. Soft Genetics LLC, State College 2009.
Rousset F. GENEPOP’007 a complete re-implementation of the Genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5. A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genoptype data. J Genet. 2000;155:945–59.
Earl DA, von Holdt BM. STRUCTURE HARVESTER. A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.
Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23:1801–6.
Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94.
Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24(11):1403–5.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.
Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A. GeneClass2: A software for genetic assignment and first-generation migrant detection. J Hered. 2004;95:536–9.
Do C, Waples RS, Peel D, Macbeth GM, Tillett BJ, Ovenden JR. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (N e ) from genetic data. Mol Ecol Resour. 2014;14:209–14.
Ravel S, Rayaisse JB, Courtin F, Solano P, De Meeûs T. Genetic signature of a recent southern range shift in Glossina tachinoides in East Burkina Faso. Infect Genet Evol. 2013;18:309–14.
Kone N, Bouyer J, Ravel S, Vreysen MJB, Domagni KT, Causse S, et al. Contrasting population structures of two vectors of African trypanosomoses in Burkina Faso: consequences for control. PLoS Negl Trop Dis. 2011;5:1217–20.
Bouyer J, Ravel S, Guerrini L, Dujardin JP, Sidibé I, Vreysen MJB, et al. 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:326–9.
Ravel S, De Meeûs T, Dujardin JP, Zeze DG, Gooding RH, Dusfour I, et al. The tsetse fly Glossina palpalis palpalis is composed of several genetically differentiated small populations in the sleeping sickness focus of Bonon, Côte d’Ivoire. Infect Genet Evol. 2007;7:116–25.
Marquez JG, Malele II, Ouma JO, Krafsur ES. Glossina swynnertoni (Diptera: Glossinidae). effective population size and breeding structure estimated by mitochondrial diversity. Bull Entomol Res. 2006;96(4):353–60.
Solano P, de La Rocque S, Cuisance D, Geoffroy B, de Meeus T, et al. Intraspecific variability in natural populations of Glossina palpalis gambiensis from West Africa, revealed by genetic and morphometric analyses. Med Vet Entomol. 1999;13:401–7.
We are grateful for the logistical support which was given by Vector and Vector Borne Disease Institute which enabled the field work to be done smoothly. Assistance provided by technicians Benedect Kimbisa and Godfrey Mashenga of Vector and Vector Borne Disease Institute on data collection was highly appreciated.
This investigation received financial support from UNICEF/UNDP/World Bank/WHO Special Programme for Research and Training in Tropical Diseases (TDR) grant no. A80132 and NIH grants D43TW007391 and AI068932. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional files.
OM, JOO, AC, IIM, GN and SA conceived the study; OM collected samples from the field, and OM and FK carried out the wet lab analyses. OM, MS and FK analyzed the data, OM, JOO, AC, GN and MS interpreted the data, JOO, GN, AC, IIM, FM and SA supervised the studies. JOO, IIM, AC, GC and SA obtained funding for the studies. OM and FK wrote the initial draft of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Matrix of geographical distances among sampling sites (km). Table S2. Genetic differentiation between all population pairs. Values in bold are significant at the 0.05 level. Table S3. Genetic summary statistics for all 19 microsatellite loci. Summary statistics are shown for each of the 19 microsatellite loci: AR, allelic richness; HO, observed heterozygosity; HE, expected heterozygosity; FIS, inbreeding coefficient and its P-value. Table S4. Genetic differentiation between the 4 clusters identified by STRUCTURE. Values in bold are significant at the 0.05 level. Key: Cluster 1 = BUK, Cluster 2 = UKE, Cluster 3 = RAS, KIR, MAS, TOB, Cluster 4 = KIS, MAN. Figure S1. Delta K Log Likelihood plot for G. f. fuscipes clusters using the second order rate of change method . The ΔK plot for a given number of clusters (K) shows that the most likely number of G. f. fuscipes clusters from the samples studied is four. Figure S2. Bayesian Information Criterion (BIC) versus the number of clusters (k) generated for discriminant analysis of principal components (DAPC; Jombart et al., ) using Adegenet (Jombart, ) for all G. f. fuscipes microsatellite MLLs. A k value of four was chosen to describe the data. Figure S3. Discriminant analysis of principal components (DAPC). Assignment of individuals from 8 sampling sites to each of the 4 identified clusters. Points represent individual genotypes sampled from a sampling site and are connected by lines to the 95% confidence ellipse centroid of the respective population. Numbers refer to the clusters identified by the STRUCTURE analyses (Fig. 3) (DOC 328 kb)