Detecting Wahlund effects together with amplification problems: Cryptic species, null alleles and short allele dominance in Glossina pallidipes populations from Tanzania

Challenges and lessons from a vector control campaign targeting Glossina palpalis palpalis in an isolated protected forest area in Abidjan, Côte d'Ivoire

Vector control (VC) is one of the strategies employed to manage African trypanosomoses. This study aimed at assessing the effectiveness of a VC campaign against Glossina palpalis palpalis using tiny targets (TTs) impregnated with insecticide in an isolated, protected forest in Abidjan, Côte d'Ivoire, while considering ecological, genetic, and operational factors. Between January 2020 and September 2022, 2,712 TTs were deployed at 684 sites, covering a total area of 1.7 km2. VC monitoring was conducted using Vavoua traps during 12 evaluation surveys, between June 2020 and March 2023. Five months after the initial TT deployment, tsetse fly density had decreased by 98.53%. Although tsetse density remained low due to TT redeployment and reinforcement, there was a significant increase a few months after the last redeployment. VC appeared to have minimal impact on the genetic structuring of G. p. palpalis. This suggested recruitment of local surviving tsetse flies all along the VC campaign due to a low probability of tsetse coming into contact with TTs, or to the evolution of behavioral or physiological resistance to control efforts. The genetic study revealed that one of the microsatellite markers used, the GPCAG locus, exhibited a selection signature possibly in response to VC. This could partly explain the challenges encountered in eliminating a seemingly isolated tsetse population thriving in a particularly favorable habitat.

Impact of long-term mass-rearing on the genetic structure of tsetse fly Glossina palpalis gambiensis colonies

Tsetse flies are the sole cyclic vectors of African trypanosomes, which cause human and animal African trypanosomiases in Africa. Tsetse fly control remains a promising option for disease management. The sterile insect technique (SIT) stands as an environmentally friendly tool to control tsetse populations. SIT requires the mass-rearing of competent sterile males to mate with wild females. However, long-term colonization might affect the genetic structure of the reared flies. This study investigated the genetic structure of four Glossina palpalis gambiensis colonies of different ages: two originating from Senegal (SEN and ICIRSEN) and two from Burkina Faso (CIR and IBD). Samples from these colonies were genotyped at ten microsatellite loci, followed by downstream population genetic analyses. The results show that the two colonies from Burkina Faso collected from close sites (∼20 km apart) over 45-year interval retained the same genetic background (FST_CIR∼IBD ≈ 0, P-value = 0.47). These flies were however, genetically different from those from the Senegal colonies (FST_CIR∼SEN ≈ 0.047; FST_IBD∼SEN ≈ 0.058, P-value = 10-4). Moreover, no significant difference was detected in the gene diversity of the CIR and IBD colonies, with HS values of 0.650 and 0.665, respectively. The inbreeding coefficient showed that all four colonies where under Hardy-Weinberg equilibrium, with FIS values of 0.026, 0.012, -0.064, and 0.001, for CIR, IBD, ICIRSEN, and SEN, respectively. Furthermore, no sign of a recent bottleneck was identified in tsetse samples from any of the four colonies. The results suggest that long-term mass-rearing of tsetse flies has no significant impact on their genetic background and diversity.

Development and characterization of microsatellite markers for population genetics of the cocoa pod borer Conopomorpha cramerella (Snellen) (Lepidoptera: Gracillaridae)

The cocoa pod borer (CPB) Conopomorpha cramerella (Snellen) (Lepidoptera: Gracillaridae) is one of the major constraints for cocoa production in South East Asia. In addition to cultural and chemical control methods, autocidal control tactics such as the Sterile Insect Technique (SIT) could be an efficient addition to the currently control strategy, however SIT implementation will depend on the population genetics of the targeted pest. The aim of the present work was to search for suitable microsatellite loci in the genome of CPB that is partially sequenced. Twelve microsatellites were initially selected and used to analyze moths collected from Indonesia, Malaysia, and the Philippines. A quality control verification process was carried out and seven microsatellites found to be suitable and efficient to distinguish differences between CPB populations from different locations. The selected microsatellites were also tested against a closely related species, i.e. the lychee fruit borer Conopomorpha sinensis (LFB) from Vietnam and eight loci were found to be suitable. The availability of these novel microsatellite loci will provide useful tools for the analysis of the population genetics and gene flow of these pests, to select suitable CPB strains to implement the SIT.

Limited impact of vector control on the population genetic structure of Glossina fuscipes fuscipes from the sleeping sickness focus of Maro, Chad

Tsetse flies (genus Glossina) transmit deadly trypanosomes to human populations and domestic animals in sub-Saharan Africa. Some foci of Human African Trypanosomiasis due to Trypanosoma brucei gambiense (g-HAT) persist in southern Chad, where a program of tsetse control was implemented against the local vector Glossina fuscipes fuscipes in 2018 in Maro. We analyzed the population genetics of G. f. fuscipes from the Maro focus before control (T0), one year (T1), and 18 months (T2) after the beginning of control efforts. Most flies captured displayed a local genetic profile (local survivors), but a few flies displayed outlier genotypes. Moreover, disturbance of isolation by distance signature (increase of genetic distance with geographic distance) and effective population size estimates, absence of any genetic signature of a bottleneck, and an increase of genetic diversity between T0 and T2 strongly suggest gene flows from various origins, and a limited impact of the vector control efforts on this tsetse population. Continuous control and surveillance of g-HAT transmission is thus recommended in Maro. Particular attention will need to be paid to the border with the Central African Republic, a country where the entomological and epidemiological status of g-HAT is unknown.

Development and characterization of microsatellite markers for the tsetse species Glossina brevipalpis and preliminary population genetics analyses

Tsetse flies, the vectors of African trypanosomes are of key medical and economic importance and one of the constraints for the development of Africa. Tsetse fly control is one of the most effective and sustainable strategies used for controlling the disease. Knowledge about population structure and level of gene flow between neighbouring populations of the target vector is of high importance to develop appropriate strategies for implementing effective management programmes. Microsatellites are commonly used to identify population structure and assess dispersal of the target populations and have been developed for several tsetse species but were lacking for Glossina brevipalpis. In this study, we screened the genome of G. brevipalpis to search for suitable microsatellite markers and nine were found to be efficient enough to distinguish between different tsetse populations. The availability of these novel microsatellite loci will help to better understand the population biology of G. brevipalpis and to assess the level of gene flow between different populations. Such information will help with the development of appropriate strategies to implement the sterile insect technique (SIT) in the framework of an area-wide integrated pest management (AW-IPM) approach to manage tsetse populations and ultimately address the trypanosomoses problem in these targeted areas.

Highly restricted dispersal in habitat-forming seaweed may impede natural recovery of disturbed populations

Cystoseira sensu lato (Class Phaeophyceae, Order Fucales, Family Sargassaceae) forests play a central role in marine Mediterranean ecosystems. Over the last decades, Cystoseira s.l. suffered from a severe loss as a result of multiple anthropogenic stressors. In particular, Gongolaria barbata has faced multiple human-induced threats, and, despite its ecological importance in structuring rocky communities and hosting a large number of species, the natural recovery of G. barbata depleted populations is uncertain. Here, we used nine microsatellite loci specifically developed for G. barbata to assess the genetic diversity of this species and its genetic connectivity among fifteen sites located in the Ionian, the Adriatic and the Black Seas. In line with strong and significant heterozygosity deficiencies across loci, likely explained by Wahlund effect, high genetic structure was observed among the three seas (ENA corrected FST = 0.355, IC = [0.283, 0.440]), with an estimated dispersal distance per generation smaller than 600 m, both in the Adriatic and Black Sea. This strong genetic structure likely results from restricted gene flow driven by geographic distances and limited dispersal abilities, along with genetic drift within isolated populations. The presence of genetically disconnected populations at small spatial scales (< 10 km) has important implications for the identification of relevant conservation and management measures for G. barbata: each population should be considered as separated evolutionary units with dedicated conservation efforts.

Negative density-dependent dispersal in tsetse (Glossina spp): red flag or red herring?

A deterministic model of the distribution of tsetse flies (Glossina spp) was used to assess the extent to which the efficacy of control operations would be affected by three different modes of density dependence in per capita adult dispersal: (i) density-independent dispersal which has been commonly adopted in previous models, (ii) positive density-dependent dispersal which has occasionally been discussed in the tsetse literature, (iii) negative density-dependent dispersal (NDDD). The last has recently been suggested, from genetic studies, to change the dispersal rate of tsetse by up to 200-fold, thereby posing a severe risk for the success of tsetse control operations. Modelling outputs showed that NDDD poses no such risk, provided the mean daily dispersal of tsetse is below about 1 km, which is greater than any rate actually recorded in the field or indicated by the genetic studies. NDDD can be problematic only if tsetse disperse at rates that appear highly unlikely, or even impossible, on energetic grounds. Under some circumstances these high rates would help rather than hinder the control officer. NDDD is not necessary to explain the results of control operations, and not sufficient to explain the results of successful control programmes.

Negative density-dependent dispersal in tsetse (Glossina spp): An artefact of inappropriate analysis

Published analysis of genetic material from field-collected tsetse (Glossina spp, primarily from the Palpalis group) has been used to predict that the distance (δ) dispersed per generation increases as effective population densities (De) decrease, displaying negative density-dependent dispersal (NDDD). Using the published data we show this result is an artefact arising primarily from errors in estimates of S, the area occupied by a subpopulation, and thereby in De. The errors arise from the assumption that S can be estimated as the area ([Formula: see text]) regarded as being covered by traps. We use modelling to show that such errors result in anomalously high correlations between [Formula: see text] and [Formula: see text] and the appearance of NDDD, with a slope of -0.5 for the regressions of log([Formula: see text]) on log([Formula: see text]), even in simulations where we specifically assume density-independent dispersal (DID). A complementary mathematical analysis confirms our findings. Modelling of field results shows, similarly, that the false signal of NDDD can be produced by varying trap deployment patterns. Errors in the estimates of δ in the published analysis were magnified because variation in estimates of S were greater than for all other variables measured, and accounted for the greatest proportion of variation in [Formula: see text]. Errors in census population estimates result from an erroneous understanding of the relationship between trap placement and expected tsetse catch, exacerbated through failure to adjust for variations in trapping intensity, trap performance, and in capture probabilities between geographical situations and between tsetse species. Claims of support in the literature for NDDD are spurious. There is no suggested explanation for how NDDD might have evolved. We reject the NDDD hypothesis and caution that the idea should not be allowed to influence policy on tsetse and trypanosomiasis control.

High dispersal capacity of Culicoides obsoletus (Diptera: Ceratopogonidae), vector of bluetongue and Schmallenberg viruses, revealed by landscape genetic analyses

BACKGROUND

In the last two decades, recurrent epizootics of bluetongue virus and Schmallenberg virus have been reported in the western Palearctic region. These viruses affect domestic cattle, sheep, goats and wild ruminants and are transmitted by native hematophagous midges of the genus Culicoides (Diptera: Ceratopogonidae). Culicoides dispersal is known to be stratified, i.e. due to a combination of dispersal processes occurring actively at short distances and passively or semi-actively at long distances, allowing individuals to jump hundreds of kilometers.

METHODS

Here, we aim to identify the environmental factors that promote or limit gene flow of Culicoides obsoletus, an abundant and widespread vector species in Europe, using an innovative framework integrating spatial, population genetics and statistical approaches. A total of 348 individuals were sampled in 46 sites in France and were genotyped using 13 newly designed microsatellite markers.

RESULTS

We found low genetic differentiation and a weak population structure for C. obsoletus across the country. Using three complementary inter-individual genetic distances, we did not detect any significant isolation by distance, but did detect significant anisotropic isolation by distance on a north-south axis. We employed a multiple regression on distance matrices approach to investigate the correlation between genetic and environmental distances. Among all the environmental factors that were tested, only cattle density seems to have an impact on C. obsoletus gene flow.

CONCLUSIONS

The high dispersal capacity of C. obsoletus over land found in the present study calls for a re-evaluation of the impact of Culicoides on virus dispersal, and highlights the urgent need to better integrate molecular, spatial and statistical information to guide vector-borne disease control.

Altitude and hillside orientation shapes the population structure of the Leishmania infantum vector Phlebotomus ariasi

Despite their role in Leishmania transmission, little is known about the organization of sand fly populations in their environment. Here, we used 11 previously described microsatellite markers to investigate the population genetic structure of Phlebotomus ariasi, the main vector of Leishmania infantum in the region of Montpellier (South of France). From May to October 2011, we captured 1,253 Ph. ariasi specimens using sticky traps in 17 sites in the North of Montpellier along a 14-km transect, and recorded the relevant environmental data (e.g., altitude and hillside). Among the selected microsatellite markers, we removed five loci because of stutter artifacts, absence of polymorphism, or non-neutral evolution. Multiple regression analyses showed the influence of altitude and hillside (51% and 15%, respectively), and the absence of influence of geographic distance on the genetic data. The observed significant isolation by elevation suggested a population structure of Ph. ariasi organized in altitudinal ecotypes with substantial rates of migration and positive assortative mating. This organization has implications on sand fly ecology and pathogen transmission. Indeed, this structure might favor the global temporal and spatial stability of sand fly populations and the spread and increase of L. infantum cases in France. Our results highlight the necessity to consider sand fly populations at small scales to study their ecology and their impact on pathogens they transmit.

Fine-scale seascape genomics of an exploited marine species, the common cockle Cerastoderma edule, using a multimodelling approach

Population dynamics of marine species that are sessile as adults are driven by oceanographic dispersal of larvae from spawning to nursery grounds. This is mediated by life-history traits such as the timing and frequency of spawning, larval behaviour and duration, and settlement success. Here, we use 1725 single nucleotide polymorphisms (SNPs) to study the fine-scale spatial genetic structure in the commercially important cockle species Cerastoderma edule and compare it to environmental variables and current-mediated larval dispersal within a modelling framework. Hydrodynamic modelling employing the NEMO Atlantic Margin Model (AMM15) was used to simulate larval transport and estimate connectivity between populations during spawning months (April-September), factoring in larval duration and interannual variability of ocean currents. Results at neutral loci reveal the existence of three separate genetic clusters (mean F ST = 0.021) within a relatively fine spatial scale in the north-west Atlantic. Environmental association analysis indicates that oceanographic currents and geographic proximity explain over 20% of the variance observed at neutral loci, while genetic variance (71%) at outlier loci was explained by sea surface temperature extremes. These results fill an important knowledge gap in the management of a commercially important and overexploited species, bringing us closer to understanding the role of larval dispersal in connecting populations at a fine geographic scale.

Developing and quality testing of microsatellite loci for four species of Glossina

Microsatellite loci still represent valuable resources for the study of the population biology of non-model organisms. Discovering or adapting new suitable microsatellite markers in species of interest still represents a useful task, especially so for non-model organisms as tsetse flies (genus Glossina), which remain a serious threat to the health of humans and animals in sub-Saharan Africa. In this paper, we present the development of new microsatellite loci for four species of Glossina: two from the Morsitans group, G. morsitans morsitans (Gmm) from Zimbabwe, G. pallidipes (Gpalli) from Tanzania; and the other two from the Palpalis group, G. fuscipes fuscipes (Gff) from Chad, and G. palpalis gambiensis (Gpg) from Guinea. We found frequent short allele dominance and null alleles. Stuttering could also be found and amended when possible. Cryptic species seemed to occur frequently in all taxa but Gff. This explains why it may be difficult finding ecumenical primers, which thus need adaptation according to each taxonomic and geographic context. Amplification problems occurred more often in published old markers, and Gmm and Gpg were the most affected (stronger heterozygote deficits). Trinucleotide markers displayed selection signature in some instances (Gmm). Combining old and new loci, for Gmm, eight loci can be safely used (with correction for null alleles); and five seem particularly promising; for Gpalli, only five to three loci worked well, depending on the clade, which means that the use of loci from other species (four morsitans loci seemed to work well), or other new primers will need to be used; for Gff, 14 loci behaved well, but with null alleles, seven of which worked very well; and for G. palpalis sl, only four loci, needing null allele and stuttering corrections seem to work well, and other loci from the literature are thus needed, including X-linked markers, five of which seem to work rather well (in females only), but new markers will probably be needed. Finally, the high proportion of X-linked markers (around 30%) was explained by the non-Y DNA quantity and chromosome structure of tsetse flies studied so far.

Strong genetic structure among populations of the tick Ixodes ricinus across its range

Ixodes ricinus is the most common and widely distributed tick species in Europe, responsible for several zoonotic diseases, including Lyme borreliosis. Population genetics of disease vectors is a useful tool for understanding the spread of pathogens and infection risks. Despite the threat to the public health due to the climate-driven distribution changes of I. ricinus, the genetic structure of tick populations, though essential for understanding epidemiology, remains unclear. Previous studies have demonstrated weak to no apparent spatial pattern of genetic differentiation between European populations. Here, we analysed the population genetic structure of 497 individuals from 28 tick populations sampled from 20 countries across Europe, the Middle-East, and northern Africa. We analysed 125 SNPs loci after quality control. We ran Bayesian and multivariate hierarchical clustering analyses to identify and describe clusters of genetically related individuals. Both clustering methods support the identification of three spatially-structured clusters. Individuals from the south and north-western parts of Eurasia form a separated cluster from northern European populations, while central European populations are a mix between the two groups. Our findings have important implications for understanding the dispersal processes that shape the spread of zoonotic diseases under anthropogenic global changes.

Phylogeography and population structure of the tsetse fly Glossina pallidipes in Kenya and the Serengeti ecosystem

Glossina pallidipes is the main vector of animal African trypanosomiasis and a potential vector of human African trypanosomiasis in eastern Africa where it poses a large economic burden and public health threat. Vector control efforts have succeeded in reducing infection rates, but recent resurgence in tsetse fly population density raises concerns that vector control programs require improved strategic planning over larger geographic and temporal scales. Detailed knowledge of population structure and dispersal patterns can provide the required information to improve planning. To this end, we investigated the phylogeography and population structure of G. pallidipes over a large spatial scale in Kenya and northern Tanzania using 11 microsatellite loci genotyped in 600 individuals. Our results indicate distinct genetic clusters east and west of the Great Rift Valley, and less distinct clustering of the northwest separate from the southwest (Serengeti ecosystem). Estimates of genetic differentiation and first-generation migration indicated high genetic connectivity within genetic clusters even across large geographic distances of more than 300 km in the east, but only occasional migration among clusters. Patterns of connectivity suggest isolation by distance across genetic breaks but not within genetic clusters, and imply a major role for river basins in facilitating gene flow in G. pallidipes. Effective population size (Ne) estimates and results from Approximate Bayesian Computation further support that there has been recent G. pallidipes population size fluctuations in the Serengeti ecosystem and the northwest during the last century, but also suggest that the full extent of differences in genetic diversity and population dynamics between the east and the west was established over evolutionary time periods (tentatively on the order of millions of years). Findings provide further support that the Serengeti ecosystem and northwestern Kenya represent independent tsetse populations. Additionally, we present evidence that three previously recognized populations (the Mbeere-Meru, Central Kenya and Coastal "fly belts") act as a single population and should be considered as a single unit in vector control.

Population genetics of Glossina palpalis palpalis in sleeping sickness foci of Côte d'Ivoire before and after vector control

Glossina palpalis palpalis remains the major vector of sleeping sickness in Côte d'Ivoire. The disease is still active at low endemic levels in Bonon and Sinfra foci in the western-central part of the country. In this study, we investigated the impact of a control campaign on G. p. palpalis population structure in Bonon and Sinfra foci in order to adapt control strategies. Genetic variation at microsatellite loci was used to examine the population structure of different G. p. palpalis cohorts before and after control campaigns. Isolation by distance was observed in our sampling sites. Before control, effective population size was high (239 individuals) with dispersal at rather short distance (731 m per generation). We found some evidence that some of the flies captured after treatment come from surrounding sites, which increased the genetic variance. One Locus, GPCAG, displayed a 1000% increase of subdivision measure after control while other loci only exhibited a substantial increase in variance of subdivision. Our data suggested a possible trap avoidance behaviour in G. p. palpalis. It is important to take into account and better understand the possible reinvasion from neighboring sites and trap avoidance for the sake of sustainability of control campaigns effects.

Parasite-mediated selection in a natural metapopulation of Daphnia magna

Parasite-mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease-related traits. However, nonadaptive processes like migration and extinction-(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life-history and disease-related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite-mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life-history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST -F_ST -like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction-(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.

Negative Density-dependent Dispersal in Tsetse Flies: A Risk for Control Campaigns?

Tsetse flies are vectors of parasites that cause diseases responsible for significant economic losses and health issues in sub-Saharan Africa, including sleeping sickness in humans and nagana in domestic animals. Efficient vector-control campaigns require good knowledge of the demographic parameters of the targeted populations. In the last decade, population genetics emerged as a convenient way to measure population densities and dispersal in tsetse flies. Here, by revealing a strong negative density-dependent dispersal in two dimensions, we suggest that control campaigns might unleash dispersal from untreated areas. If confirmed by direct measurement of dispersal before and after control campaigns, area-wide and/or sequential treatments of neighboring sites will be necessary to prevent this issue.