Immunogenicity and Serological Cross-Reactivity of Saliva Proteins among Different Tsetse Species

Tsetse are vectors of pathogenic trypanosomes, agents of human and animal trypanosomiasis in Africa. Components of tsetse saliva (sialome) are introduced into the mammalian host bite site during the blood feeding process and are important for tsetse’s ability to feed efficiently, but can also influence disease transmission and serve as biomarkers for host exposure. We compared the sialome components from four tsetse species in two subgenera: subgenus Morsitans: Glossina morsitans morsitans (Gmm) and Glossina pallidipes (Gpd), and subgenus Palpalis: Glossina palpalis gambiensis (Gpg) and Glossina fuscipes fuscipes (Gff), and evaluated their immunogenicity and serological cross reactivity by an immunoblot approach utilizing antibodies from experimental mice challenged with uninfected flies. The protein and immune profiles of sialome components varied with fly species in the same subgenus displaying greater similarity and cross reactivity. Sera obtained from cattle from disease endemic areas of Africa displayed an immunogenicity profile reflective of tsetse species distribution. We analyzed the sialome fractions of Gmm by LC-MS/MS, and identified TAg5, Tsal1/Tsal2, and Sgp3 as major immunogenic proteins, and the 5'-nucleotidase family as well as four members of the Adenosine Deaminase Growth Factor (ADGF) family as the major non-immunogenic proteins. Within the ADGF family, we identified four closely related proteins (TSGF-1, TSGF-2, ADGF-3 and ADGF-4), all of which are expressed in tsetse salivary glands. We describe the tsetse species-specific expression profiles and genomic localization of these proteins. Using a passive-immunity approach, we evaluated the effects of rec-TSGF (TSGF-1 and TSGF-2) polyclonal antibodies on tsetse fitness parameters. Limited exposure of tsetse to mice with circulating anti-TSGF antibodies resulted in a slight detriment to their blood feeding ability as reflected by compromised digestion, lower weight gain and less total lipid reserves although these results were not statistically significant. Long-term exposure studies of tsetse flies to antibodies corresponding to the ADGF family of proteins are warranted to evaluate the role of this conserved family in fly biology.

Insect saliva contains many proteins that are injected into the mammalian host during the blood feeding process. Saliva proteins enhance the blood feeding ability of insects, but they can also induce mammalian immune responses that inhibit successful feeding, or modulate the bite site to benefit pathogen transmission. Here we studied saliva from four different tsetse species that belong to two distant species groups. We show that the saliva protein profiles of different species groups vary. Experimental mice subjected to fly bites display varying immunological responses against the abundant saliva proteins and the antigenicity of the shared saliva proteins in different tsetse species differs. We show that one member of the ADGF family with adenosine deaminase motifs, TSGF-2, is non-immunogenic in Glossina morsitans in mice, while the same protein from Glossina fuscipes is highly immunogenic. Such species-specific immune responses could be exploited as biomarkers of host exposures in the field. We also show that short-term exposure of G. morsitans to mice passively immunized by anti-TSGF antibodies leads to slight but not statistically significant negative fitness effects. Thus, future investigations with non-antigenic saliva proteins are warranted as they can lead to novel mammalian vaccine targets to reduce tsetse populations in the field.

Mitochondrial DNA sequence divergence and diversity of Glossina fuscipes fuscipes in the Lake Victoria basin of Uganda: implications for control

BACKGROUND

Glossina fuscipes fuscipes is the main vector of African Trypanosomiasis affecting both humans and livestock in Uganda. The human disease (sleeping sickness) manifests itself in two forms: acute and chronic. The Lake Victoria basin in Uganda has the acute form and a history of tsetse re-emergence despite concerted efforts to control tsetse. The government of Uganda has targeted the basin for tsetse eradication. To provide empirical data for this initiative, we screened tsetse flies from the basin for genetic variation at the mitochondrial DNA cytochrome oxidase II (mtDNA COII) gene with the goal of investigating genetic diversity and gene flow among tsetse, tsetse demographic history; and compare these results with results from a previous study based on microsatellite loci data in the same area.

METHODS

We collected 429 Gff tsetse fly samples from 14 localities in the entire Ugandan portion of the Lake Victoria coast, covering 40,000 km(2). We performed genetic analyses on them and added data collected for 56 Gff individuals from 4 additional sampling sites in the basin. The 529 pb partial mitochondrial DNA cytochrome oxidase II (mtDNA COII) sequences totaling 485 were analysed for genetic differentiation, structuring and demographic history. The results were compared with findings from a previous study based on microsatellite loci data from the basin.

RESULTS

The differences within sampling sites explained a significant proportion of the genetic variation. We found three very closely related mtDNA population clusters, which co-occurred in multiple sites. Although Φ ST (0 - 0.592; P < 0.05) and Bayesian analyses suggest some level of weak genetic differentiation, there is no correlation between genetic divergence and geographic distance (r = 0.109, P = 0.185), and demographic tests provide evidence of locality-based demographic history.

CONCLUSION

The mtDNA data analysed here complement inferences made in a previous study based on microsatellite data. Given the differences in mutation rates, mtDNA afforded a look further back in time than microsatellites and revealed that Gff populations were more connected in the past. Microsatellite data revealed more genetic structuring than mtDNA. The differences in connectedness and structuring over time could be related to vector control efforts. Tsetse re-emergence after control interventions may be due to re-invasions from outside the treated areas, which emphasizes the need for an integrated area-wide tsetse eradication strategy for sustainable removal of the tsetse and trypanosomiasis problem from this area.

Wolbachia association with the tsetse fly, Glossina fuscipes fuscipes, reveals high levels of genetic diversity and complex evolutionary dynamics

Background

Wolbachia pipientis, a diverse group of α-proteobacteria, can alter arthropod host reproduction and confer a reproductive advantage to Wolbachia-infected females (cytoplasmic incompatibility (CI)). This advantage can alter host population genetics because Wolbachia-infected females produce more offspring with their own mitochondrial DNA (mtDNA) haplotypes than uninfected females. Thus, these host haplotypes become common or fixed (selective sweep). Although simulations suggest that for a CI-mediated sweep to occur, there must be a transient phase with repeated initial infections of multiple individual hosts by different Wolbachia strains, this has not been observed empirically. Wolbachia has been found in the tsetse fly, Glossina fuscipes fuscipes, but it is not limited to a single host haplotype, suggesting that CI did not impact its population structure. However, host population genetic differentiation could have been generated if multiple Wolbachia strains interacted in some populations. Here, we investigated Wolbachia genetic variation in G. f. fuscipes populations of known host genetic composition in Uganda. We tested for the presence of multiple Wolbachia strains using Multi-Locus Sequence Typing (MLST) and for an association between geographic region and host mtDNA haplotype using Wolbachia DNA sequence from a variable locus, groEL (heat shock protein 60).

Results

MLST demonstrated that some G. f. fuscipes carry Wolbachia strains from two lineages. GroEL revealed high levels of sequence diversity within and between individuals (Haplotype diversity = 0.945). We found Wolbachia associated with 26 host mtDNA haplotypes, an unprecedented result. We observed a geographical association of one Wolbachia lineage with southern host mtDNA haplotypes, but it was non-significant (p = 0.16). Though most Wolbachia-infected host haplotypes were those found in the contact region between host mtDNA groups, this association was non-significant (p = 0.17).

Conclusions

High Wolbachia sequence diversity and the association of Wolbachia with multiple host haplotypes suggest that different Wolbachia strains infected G. f. fuscipes multiple times independently. We suggest that these observations reflect a transient phase in Wolbachia evolution that is influenced by the long gestation and low reproductive output of tsetse. Although G. f. fuscipes is superinfected with Wolbachia, our data does not support that bidirectional CI has influenced host genetic diversity in Uganda.

The population structure of Glossina fuscipes fuscipes in the Lake Victoria basin in Uganda: implications for vector control

Background

Glossina fuscipes fuscipes is the primary vector of trypanosomiasis in humans and livestock in Uganda. The Lake Victoria basin has been targeted for tsetse eradication using a rolling carpet initiative, from west to east, with four operational blocks (3 in Uganda and 1 in Kenya), under a Pan-African Tsetse and Trypanosomiasis Eradication Campaign (PATTEC). We screened tsetse flies from the three Ugandan PATTEC blocks for genetic diversity at 15 microsatellite loci from continental and offshore populations to provide empirical data to support this initiative.

Methods

We collected tsetse samples from 11 sites across the Lake Victoria basin in Uganda. We performed genetic analyses on 409 of the collected tsetse flies and added data collected for 278 individuals in a previous study. The flies were screened across 15 microsatellite loci and the resulting data were used to assess the temporal stability of populations, to analyze patterns of genetic exchange and structuring, to estimate dispersal rates and evaluate the sex bias in dispersal, as well as to estimate demographic parameters (N_E and N_C).

Results

We found that tsetse populations in this region were stable over 4-16 generations and belong to 4 genetic clusters. Two genetic clusters (1 and 2) corresponded approximately to PATTEC blocks 1 and 2, while the other two (3 and 4) fell within PATTEC block 3. Island populations grouped into the same genetic clusters as neighboring mainland sites, suggesting presence of gene flow between these sites. There was no evidence of the stretch of water separating islands from the mainland forming a significant barrier to dispersal. Dispersal rates ranged from 2.5 km per generation in cluster 1 to 14 km per generation in clusters 3 and 4. We found evidence of male-biased dispersal. Few breeders are successfully dispersing over large distances. Effective population size estimates were low (33–310 individuals), while census size estimates ranged from 1200 (cluster 1) to 4100 (clusters 3 and 4). We present here a novel technique that adapts an existing census size estimation method to sampling without replacement, the scheme used in sampling tsetse flies.

Conclusion

Our study suggests that different control strategies should be implemented for the three PATTEC blocks and that, given the high potential for re-invasion from island sites, mainland and offshore sites in each block should be targeted at the same time.

Cryptic diversity within the major trypanosomiasis vector Glossina fuscipes revealed by molecular markers

BACKGROUND

The tsetse fly Glossina fuscipes s.l. is responsible for the transmission of approximately 90% of cases of human African trypanosomiasis (HAT) or sleeping sickness. Three G. fuscipes subspecies have been described, primarily based upon subtle differences in the morphology of their genitalia. Here we describe a study conducted across the range of this important vector to determine whether molecular evidence generated from nuclear DNA (microsatellites and gene sequence information), mitochondrial DNA and symbiont DNA support the existence of these taxa as discrete taxonomic units.

PRINCIPAL FINDINGS

The nuclear ribosomal Internal transcribed spacer 1 (ITS1) provided support for the three subspecies. However nuclear and mitochondrial sequence data did not support the monophyly of the morphological subspecies G. f. fuscipes or G. f. quanzensis. Instead, the most strongly supported monophyletic group was comprised of flies sampled from Ethiopia. Maternally inherited loci (mtDNA and symbiont) also suggested monophyly of a group from Lake Victoria basin and Tanzania, but this group was not supported by nuclear loci, suggesting different histories of these markers. Microsatellite data confirmed strong structuring across the range of G. fuscipes s.l., and was useful for deriving the interrelationship of closely related populations.

CONCLUSION/SIGNIFICANCE

We propose that the morphological classification alone is not used to classify populations of G. fuscipes for control purposes. The Ethiopian population, which is scheduled to be the target of a sterile insect release (SIT) programme, was notably discrete. From a programmatic perspective this may be both positive, given that it may reflect limited migration into the area or negative if the high levels of differentiation are also reflected in reproductive isolation between this population and the flies to be used in the release programme.

Genetic diversity and population structure of Glossina pallidipes in Uganda and western Kenya

Background

Glossina pallidipes has been implicated in the spread of sleeping sickness from southeastern Uganda into Kenya. Recent studies indicated resurgence of G. pallidipes in Lambwe Valley and southeastern Uganda after what were deemed to be effective control efforts. It is unknown whether the G. pallidipes belt in southeastern Uganda extends into western Kenya. We investigated the genetic diversity and population structure of G. pallidipes in Uganda and western Kenya.

Results

AMOVA indicated that differences among sampling sites explained a significant proportion of the genetic variation. Principal component analysis and Bayesian assignment of microsatellite genotypes identified three distinct clusters: western Uganda, southeastern Uganda/Lambwe Valley, and Nguruman in central-southern Kenya. Analyses of mtDNA confirmed the results of microsatellite analysis, except in western Uganda, where Kabunkanga and Murchison Falls populations exhibited haplotypes that differed despite homogeneous microsatellite signatures. To better understand possible causes of the contrast between mitochondrial and nuclear markers we tested for sex-biased dispersal. Mean pairwise relatedness was significantly higher in females than in males within populations, while mean genetic distance was lower and relatedness higher in males than females in between-population comparisons. Two populations sampled on the Kenya/Uganda border, exhibited the lowest levels of genetic diversity. Microsatellite alleles and mtDNA haplotypes in these two populations were a subset of those found in neighboring Lambwe Valley, suggesting that Lambwe was the source population for flies in southeastern Uganda. The relatively high genetic diversity of G. pallidipes in Lambwe Valley suggest large relict populations remained even after repeated control efforts.

Conclusion

Our research demonstrated that G. pallidipes populations in Kenya and Uganda do not form a contiguous tsetse belt. While Lambwe Valley appears to be a source population for flies colonizing southeastern Uganda, this dispersal does not extend to western Uganda. The complicated phylogeography of G. pallidipes warrants further efforts to distinguish the role of historical and modern gene flow and possible sex-biased dispersal in structuring populations.

Polyandry is a common event in wild populations of the Tsetse fly Glossina fuscipes fuscipes and may impact population reduction measures

BACKGROUND

Glossina fuscipes fuscipes is the main vector of human and animal trypanosomiasis in Africa, particularly in Uganda. Attempts to control/eradicate this species using biological methods require knowledge of its reproductive biology. An important aspect is the number of times a female mates in the wild as this influences the effective population size and may constitute a critical factor in determining the success of control methods. To date, polyandry in G.f. fuscipes has not been investigated in the laboratory or in the wild. Interest in assessing the presence of remating in Ugandan populations is driven by the fact that eradication of this species is at the planning stage in this country.

METHODOLOGY/PRINCIPAL FINDINGS

Two well established populations, Kabukanga in the West and Buvuma Island in Lake Victoria, were sampled to assess the presence and frequency of female remating. Six informative microsatellite loci were used to estimate the number of matings per female by genotyping sperm preserved in the female spermathecae. The direct count of the minimum number of males that transferred sperm to the spermathecae was compared to Maximum Likelihood and Bayesian probability estimates. The three estimates provided evidence that remating is common in the populations but the frequency is substantially different: 57% in Kabukanga and 33% in Buvuma.

CONCLUSIONS/SIGNIFICANCE

The presence of remating, with females maintaining sperm from different mates, may constitute a critical factor in cases of re-infestation of cleared areas and/or of residual populations. Remating may enhance the reproductive potential of re-invading propagules in terms of their effective population size. We suggest that population age structure may influence remating frequency. Considering the seasonal demographic changes that this fly undergoes during the dry and wet seasons, control programmes based on SIT should release large numbers of sterile males, even in residual surviving target populations, in the dry season.

Permanent genetic resources added to Molecular Ecology Resources Database 1 December 2010-31 January 2011

This article documents the addition of 238 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Alytes dickhilleni, Arapaima gigas, Austropotamobius italicus, Blumeria graminis f. sp. tritici, Cobitis lutheri, Dendroctonus ponderosae, Glossina morsitans morsitans, Haplophilus subterraneus, Kirengeshoma palmata, Lysimachia japonica, Macrolophus pygmaeus, Microtus cabrerae, Mytilus galloprovincialis, Pallisentis (Neosentis) celatus, Pulmonaria officinalis, Salminus franciscanus, Thais chocolata and Zootoca vivipara. These loci were cross-tested on the following species: Acanthina monodon, Alytes cisternasii, Alytes maurus, Alytes muletensis, Alytes obstetricans almogavarii, Alytes obstetricans boscai, Alytes obstetricans obstetricans, Alytes obstetricans pertinax, Cambarellus montezumae, Cambarellus zempoalensis, Chorus giganteus, Cobitis tetralineata, Glossina fuscipes fuscipes, Glossina pallidipes, Lysimachia japonica var. japonica, Lysimachia japonica var. minutissima, Orconectes virilis, Pacifastacus leniusculus, Procambarus clarkii, Salminus brasiliensis and Salminus hilarii.

Temporal stability of Glossina fuscipes fuscipes populations in Uganda

BACKGROUND

Glossina fuscipes, a riverine species of tsetse, is the major vector of human African trypanosomiasis (HAT) in sub-Saharan Africa. Understanding the population dynamics, and specifically the temporal stability, of G. fuscipes will be important for informing vector control activities. We evaluated genetic changes over time in seven populations of the subspecies G. f. fuscipes distributed across southeastern Uganda, including a zone of contact between two historically isolated lineages. A total of 667 tsetse flies were genotyped at 16 microsatellite loci and at one mitochondrial locus.

RESULTS

Results of an AMOVA indicated that time of sampling did not explain a significant proportion of the variance in allele frequencies observed across all samples. Estimates of differentiation between samples from a single population ranged from approximately 0 to 0.019, using Jost's DEST. Effective population size estimates using momentum-based and likelihood methods were generally large. We observed significant change in mitochondrial haplotype frequencies in just one population, located along the zone of contact. The change in haplotypes was not accompanied by changes in microsatellite frequencies, raising the possibility of asymmetric mating compatibility in this zone.

CONCLUSION

Our results suggest that populations of G. f. fuscipes were stable over the 8-12 generations studied. Future studies should aim to reconcile these data with observed seasonal fluctuations in the apparent density of tsetse.

Pyrethroid resistance in an Anopheles funestus population from Uganda

BACKGROUND

The susceptibility status of Anopheles funestus to insecticides remains largely unknown in most parts of Africa because of the difficulty in rearing field-caught mosquitoes of this malaria vector. Here we report the susceptibility status of the An. funestus population from Tororo district in Uganda and a preliminary characterisation of the putative resistance mechanisms involved.

METHODOLOGY/PRINCIPAL FINDINGS

A new forced egg laying technique used in this study significantly increased the numbers of field-caught females laying eggs and generated more than 4000 F1 adults. WHO bioassays indicated that An. funestus in Tororo is resistant to pyrethroids (62% mortality after 1 h exposure to 0.75% permethrin and 28% mortality to 0.05% deltamethrin). Suspected DDT resistance was also observed with 82% mortality. However this population is fully susceptible to bendiocarb (carbamate), malathion (organophosphate) and dieldrin with 100% mortality observed after exposure to each of these insecticides. Sequencing of a fragment of the sodium channel gene containing the 1014 codon conferring pyrethroid/DDT resistance in An. gambiae did not detect the L1014F kdr mutation but a correlation between haplotypes and resistance phenotype was observed indicating that mutations in other exons may be conferring the knockdown resistance in this species. Biochemical assays suggest that resistance in this population is mediated by metabolic resistance with elevated level of GSTs, P450s and pNPA compared to a susceptible strain of Anopheles gambiae. RT-PCR further confirmed the involvement of P450s with a 12-fold over-expression of CYP6P9b in the Tororo population compared to the fully susceptible laboratory colony FANG.

CONCLUSION

This study represents the first report of pyrethroid/DDT resistance in An. funestus from East Africa. With resistance already reported in southern and West Africa, this indicates that resistance in An. funestus may be more widespread than previously assumed and therefore this should be taken into account for the implementation and management of vector control programs in Africa.

Phylogeography and population structure of Glossina fuscipes fuscipes in Uganda: implications for control of tsetse

Background

Glossina fuscipes fuscipes, a riverine species of tsetse, is the main vector of both human and animal trypanosomiasis in Uganda. Successful implementation of vector control will require establishing an appropriate geographical scale for these activities. Population genetics can help to resolve this issue by characterizing the extent of linkage among apparently isolated groups of tsetse.

Methodology/Principal Findings

We conducted genetic analyses on mitochondrial and microsatellite data accumulated from approximately 1000 individual tsetse captured in Uganda and neighboring regions of Kenya and Sudan. Phylogeographic analyses suggested that the largest scale genetic structure in G. f. fuscipes arose from an historical event that divided two divergent mitochondrial lineages. These lineages are currently partitioned to northern and southern Uganda and co-occur only in a narrow zone of contact extending across central Uganda. Bayesian assignment tests, which provided evidence for admixture between northern and southern flies at the zone of contact and evidence for northerly gene flow across the zone of contact, indicated that this structure may be impermanent. On the other hand, microsatellite structure within the southern lineage indicated that gene flow is currently limited between populations in western and southeastern Uganda. Within regions, the average F_ST between populations separated by less than 100 km was less than ∼0.1. Significant tests of isolation by distance suggested that gene flow is ongoing between neighboring populations and that island populations are not uniformly more isolated than mainland populations.

Conclusions/Significance

Despite the presence of population structure arising from historical colonization events, our results have revealed strong signals of current gene flow within regions that should be accounted for when planning tsetse control in Uganda. Populations in southeastern Uganda appeared to receive little gene flow from populations in western or northern Uganda, supporting the feasibility of area wide control in the Lake Victoria region by the Pan African Tsetse and Trypanosomiasis Eradication Campaign.

Glossina fuscipes fuscipes is the most common species of tsetse in Uganda, where it transmits human sleeping sickness and nagana, a related disease of cattle. A consortium of African countries dedicated to controlling these diseases is poised to begin area wide control of tsetse, but a critical question remains: What is the most appropriate geographical scale for these activities? To address this question, we used population genetics to determine the extent of linkage between populations of tsetse confined to discrete patches of riverine habitat. Our results suggest that Uganda was colonized by two distinct lineages of G. f. fuscipes, which now co-occur only in a narrow band across central Uganda. Evidence for interbreeding at the zone of contact and movement of genes from the south to the north suggest that this historical genetic structure may dissolve in the future. At smaller scales, we have demonstrated that exchange of genes among neighboring populations via dispersal is at equilibrium with the differentiating force of genetic drift. Our results highlight the need for investment in vector control programs that account for the linkage observed among tsetse populations. Given its genetic isolation and its location at the far edge of G. fuscipes' range, the Lake Victoria region appears to be an appropriate target for area wide control.

Insecticide resistance and its association with target-site mutations in natural populations of Anopheles gambiae from eastern Uganda

Insecticide resistance in Anopheles gambiae threatens the success of malaria vector control programmes in sub-Saharan Africa. In order to manage insecticide resistance successfully, it is essential to assess continuously the target mosquito population. Here, we collected baseline information on the distribution and prevalence of insecticide resistance and its association with target-site mutations in eastern Uganda. Anopheles gambiae s.l. adults were raised from wild-caught larvae sampled from two ecologically distinct breeding sites and exposed to WHO discriminating concentrations of DDT, permethrin, deltamethrin, bendiocarb and malathion. Survival rates to DDT were as high as 85.4%, alongside significant resistance levels to permethrin (38.5%), reduced susceptibility to deltamethrin, but full susceptibility to bendiocarb and malathion. Using molecular diagnostics, susceptible and resistant specimens were further tested for the presence of knockdown resistance (kdr) and acetylcholinesterase 1 resistance (ace-1(R)) alleles. While ace-1(R) and kdrL1014F ('kdr west') alleles were absent, the kdr L1014S ('kdr east') allele was present in both populations. In A. gambiae s.s., L1014S was closely associated with DDT and, to a lesser degree, with permethrin resistance. Intriguingly, the association between DDT resistance and the presence of L1014S is consistent with a co-dominant effect, with heterozygous individuals showing an intermediate phenotype.

Molecular phylogenetics of tsetse flies (Diptera: Glossinidae) based on mitochondrial (COI, 16S, ND2) and nuclear ribosomal DNA sequences, with an emphasis on the palpalis group

Relationships of 13 species of the genus Glossina (tsetse flies) were inferred from mitochondrial (cytochrome oxidase 1, NADH dehydrogenase 2 and 16S) and nuclear (internal transcribed spacer 1 of rDNA) sequences. The resulting phylogeny confirms the monophyly of the morphologically defined fusca, morsitans and palpalis subgenera. Genetic distances between palpalis and morsitans subspecies suggest that their status needs revision. In particular, cytochrome oxidase 1 sequences showed large geographical differences within G. palpalis palpalis, suggesting the existence of cryptic species within this subspecies. The morphology of palpalis group female genital plates was examined, and individuals were found varying outside the ranges specified by the standard identification keys, making definitive morphological classification impossible. A diagnostic PCR to distinguish G. palpalis palpalis, G. tachinoides and G. palpalis gambiensis based on length differences of internal transcribed spacer 1 sequences is presented.