Population genetic structure of Anopheles gambiae and Anopheles arabiensis in Niger

Tracing the origin of the early wet-season Anopheles coluzzii in the Sahel

In arid environments, the source of the malaria mosquito populations that re‐establish soon after first rains remains a puzzle and alternative explanations have been proposed. Using genetic data, we evaluated whether the early rainy season (RS) population of Anopheles coluzzii is descended from the preceding late RS generation at the same locality, consistent with dry season (DS) dormancy (aestivation), or from migrants from distant locations. Distinct predictions derived from these two hypotheses were assessed, based on variation in 738 SNPs in eleven A. coluzzii samples, including seven samples spanning 2 years in a Sahelian village. As predicted by the “local origin under aestivation hypothesis,” temporal samples from the late RS and those collected after the first rain of the following RS were clustered together, while larger genetic distances were found among samples spanning the RS. Likewise, multilocus genotype composition of samples from the end of the RS was similar across samples until the following RS, unlike samples that spanned the RS. Consistent with reproductive arrest during the DS, no genetic drift was detected between samples taken over that period, despite encompassing extreme population minima, whereas it was detected between samples spanning the RS. Accordingly, the variance in allele frequency increased with time over the RS, but not over the DS. However, not all the results agreed with aestivation. Large genetic distances separated samples taken a year apart, and during the first year, within‐sample genetic diversity declined and increased back during the late RS, suggesting a bottleneck followed by migration. The decline of genetic diversity followed by a mass distribution of insecticide‐treated nets was accompanied by a reduced mosquito density and a rise in the mutation conferring resistance to pyrethroids, indicating a bottleneck due to insecticidal selection. Overall, our results support aestivation in A. coluzzii during the DS that is accompanied by long‐distance migration in the late RS.

Islands and stepping-stones: comparative population structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and implications for the spread of insecticide resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise FST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average FST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.

Diversity, differentiation, and linkage disequilibrium: prospects for association mapping in the malaria vector Anopheles arabiensis

Association mapping is a widely applied method for elucidating the genetic basis of phenotypic traits. However, factors such as linkage disequilibrium and levels of genetic diversity influence the power and resolution of this approach. Moreover, the presence of population subdivision among samples can result in spurious associations if not accounted for. As such, it is useful to have a detailed understanding of these factors before conducting association mapping experiments. Here we conducted whole-genome sequencing on 24 specimens of the malaria mosquito vector, Anopheles arabiensis, to further understanding of patterns of genetic diversity, population subdivision and linkage disequilibrium in this species. We found high levels of genetic diversity within the An. arabiensis genome, with ~800,000 high-confidence, single- nucleotide polymorphisms detected. However, levels of nucleotide diversity varied significantly both within and between chromosomes. We observed lower diversity on the X chromosome, within some inversions, and near centromeres. Population structure was absent at the local scale (Kilombero Valley, Tanzania) but detected between distant populations (Cameroon vs. Tanzania) where differentiation was largely restricted to certain autosomal chromosomal inversions such as 2Rb. Overall, linkage disequilibrium within An. arabiensis decayed very rapidly (within 200 bp) across all chromosomes. However, elevated linkage disequilibrium was observed within some inversions, suggesting that recombination is reduced in those regions. The overall low levels of linkage disequilibrium suggests that association studies in this taxon will be very challenging for all but variants of large effect, and will require large sample sizes.

Longitudinal follow-up of malaria transmission dynamics in two villages in a Sahelian area of Niger during a nationwide insecticide-treated bednet distribution programme

Malaria transmission was monitored in two villages in the Sahel zone of Niger over 4 years. During this period, a nationwide vector control programme was carried out in which insecticide-treated bednets were distributed free to mothers of children aged <5 years. Anopheles gambiae and Anopheles arabiensis (Diptera: Culicidae) were found to be the major malaria vectors. The dynamics of An. gambiae s.l. did not vary dramatically over the study period although the proportion of female mosquitoes found resting indoors decreased in both villages and, in one village, the parity rate and sporozoite index were significantly reduced after bednet distribution. By contrast with An. gambiae, the dynamics of Anopheles funestus altered greatly after the bednet distribution period, when adult density, endophagous rate and sporozoite rates decreased dramatically. Our observations highlight the importance of quantifying and monitoring the dynamics and infections of malaria vectors during large-scale vector control interventions.