Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2014 Feb 13;8(2):e2697.
doi: 10.1371/journal.pntd.0002697. eCollection 2014 Feb.

Laboratory colonisation and genetic bottlenecks in the tsetse fly Glossina pallidipes

Affiliations

Laboratory colonisation and genetic bottlenecks in the tsetse fly Glossina pallidipes

Marc Ciosi et al. PLoS Negl Trop Dis. .

Abstract

Background: The IAEA colony is the only one available for mass rearing of Glossina pallidipes, a vector of human and animal African trypanosomiasis in eastern Africa. This colony is the source for Sterile Insect Technique (SIT) programs in East Africa. The source population of this colony is unclear and its genetic diversity has not previously been evaluated and compared to field populations.

Methodology/principal findings: We examined the genetic variation within and between the IAEA colony and its potential source populations in north Zimbabwe and the Kenya/Uganda border at 9 microsatellites loci to retrace the demographic history of the IAEA colony. We performed classical population genetics analyses and also combined historical and genetic data in a quantitative analysis using Approximate Bayesian Computation (ABC). There is no evidence of introgression from the north Zimbabwean population into the IAEA colony. Moreover, the ABC analyses revealed that the foundation and establishment of the colony was associated with a genetic bottleneck that has resulted in a loss of 35.7% of alleles and 54% of expected heterozygosity compared to its source population. Also, we show that tsetse control carried out in the 1990's is likely reduced the effective population size of the Kenya/Uganda border population.

Conclusions/significance: All the analyses indicate that the area of origin of the IAEA colony is the Kenya/Uganda border and that a genetic bottleneck was associated with the foundation and establishment of the colony. Genetic diversity associated with traits that are important for SIT may potentially have been lost during this genetic bottleneck which could lead to a suboptimal competitiveness of the colony males in the field. The genetic diversity of the colony is lower than that of field populations and so, studies using colony flies should be interpreted with caution when drawing general conclusions about G. pallidipes biology.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Figure 1
Figure 1. Competing scenarios considered in the ABC analysis of IAEA colony past demography (analysis 3).
The demography of the Busia and Rukomeshi populations was determined as described in Supplementary file S2. In all scenarios, the two potential source populations merged Tanc. ago into an unsampled ancestral population and a bottleneck started in Busia in 1991 and lasted BD cont. generations. Scenario 4, 5 and 6 are respectively variations of scenarios 1, 2 and 3 in which a genetic bottleneck (of duration BD col.Bus. or BD col.Ruk.) is associated with the laboratory colonisation. In scenarios 1, 2, 4 and 5 the IAEA colony has a single population of origin. In scenarios 1 and 4, the IAEA colony was founded from the Busia population in 1975 while in scenarios 2 and 5, it was founded from the Rukomeshi population Tcol.Ruk. generations ago. In scenario 3 and 6 the IAEA colony originates from an admixture between unsampled colonies of Busia and Rukomeshi origin that were respectively founded from the Busia population in 1975 and from the Rukomeshi population Tcol.Ruk. generations ago. When admixture occurs, the admixture rate ar is the proportion from unsampled Busia colony that contributed to the admixed IAEA colony.
Figure 2
Figure 2. Genetic variations within samples.
Error bars indicate the standard deviations across loci. Na: average number of alleles per locus. AR: allelic richness. AR is based on minimum sample size (N = 23 in Friuli for locus DVV-ET1). H: mean expected heterozygosity.
Figure 3
Figure 3. Distribution of the allele frequency classes.
(A) in Busia. (B) in the IAEA colony samples from 2012 and 2013. (C) in Rukomeshi. The result of the bottleneck tests (Wilcoxon' tests on heterozygosity excess, and the mode shift tests) are indicated for each panel which correspond to a field population or to the IAEA colony. Heterozygosity excess and/or a shifted mode in the distribution of allele frequency classes indicate a recent reduction in population size.
Figure 4
Figure 4. Estimated number of clusters and population structure from the Structure analysis.
(A) Mean (±SD) natural logarithm of the likelihood of the data [LnP(X|K)] over 10 Structure replicated runs for each value of the putative number of clusters (K). (B) Estimated population structure from the Structure analysis for K = 2 and 3. Each individual is represented by a vertical line divided into K coloured segments that represent the individual's estimated membership fractions in K clusters. Black lines separate individuals from different samples. Each plot, is based on the highest-probability run (among ten) at K = 2 and 3.
Figure 5
Figure 5. Effective population sizes (Ne) estimated from the ABC analysis with scenario 4 (Figure 1).
The prior distributions are shown as black lines while the posterior distributions are shown as coloured lines using the same colour code as in Figure 1. The medians of a posterior distribution, considered as point estimate of the parameters, are indicated. Each distribution was obtained from 10000 values. (A) N Bus.: Ne of the Busia population. N Ruk.: Ne of the Rukomeshi population. N IAEA.: Ne of the IAEA colony. The prior distributions for N Bus. and N Ruk. are shown as a plain line while the prior distribution for N IAEA. is shown as a dashed line. (B) NBcol. Bus.: Ne of the colony of Busia origin (which is at the origin of the current IAEA colony) during the bottleneck associated with its establishment.
Figure 6
Figure 6. Loss of genetic diversity in the IAEA colony (2013 sample) with respect to its simulated source population at the Kenya/Uganda border in 1975.
Na; allelic diversity loss and H; expected heterozygosity loss. Both losses are statistically significant (Wilcoson's signed-rank tests over loci, p≤0.05).

References

    1. Holmes P (2013) Tsetse-transmitted trypanosomes - Their biology, disease impact and control. J Invertebr Pathol 112: S11–S14. - PubMed
    1. Vreysen MJB, Saleh KM, Lancelot R, Bouyer J (2011) Factory Tsetse Flies Must Behave Like Wild Flies: A Prerequisite for the Sterile Insect Technique. PLoS Neglect Trop D 5: e907. - PMC - PubMed
    1. Klassen W, Curtis CF (2005) History of the Sterile Insect Technique. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile Insect Technique: Springer Netherlands. pp. 3–36.
    1. Krafsur ES (2005) Role of Population Genetics in the Sterile Insect Technique. In: Dyck VA, Hendrichs J, Robinson AS, editors. Sterile Insect Technique: Springer Netherlands. pp. 389–406.
    1. Gooding RH (1990) Genetic aspects of quality control in tsetse colonies. Insect Sci Appl 11: 385–398.

Publication types