David Weetman

Genetic population structure and contemporary dispersal patterns of a recent European invader, the Chinese mitten crab, Eriocheir sinensis

Genetic studies of recently established populations are challenging because the assumption of equilibrium underlying many analyses is likely to be violated. Using microsatellites, we investigated determinants of genetic structure and migration among invasive European‐Chinese mitten crab populations, applying a combination of traditional population genetic analyses and nonequilibrium Bayesian methods. Consistent with their recent history, invasive populations showed much lower levels of genetic diversity than a native Chinese population, indicative of recent bottlenecks. Population differentiation was generally low but significant and especially pronounced among recently established populations. Significant differentiation among cohorts from the same geographical location (River Thames) suggests the low effective population size and associated strong genetic drift that would be anticipated from a very recent colonization. An isolation‐by‐distance pattern appears to be driven by an underlying correlation between geographical distance and population age, suggesting that cumulative homogenizing gene flow reduces founder bottleneck‐associated genetic differentiation between longer‐established populations. This hypothesis was supported by a coalescent analysis, which supported a drift + gene flow model as more likely than a model excluding gene flow. Furthermore, admixture analysis identified several recent migrants between the UK and Continental European population clusters. Admixture proportions were significantly predicted by the volume of shipping between sites, indicating that human‐mediated transport remains a significant factor for dispersal of mitten crabs after the initial establishment of populations. Our study highlights the value of nonequilibrium methods for the study of invasive species, and also the importance of evaluating nonequilibrium explanations for isolation by distance patterns.

Multiple-Insecticide Resistance in Anopheles gambiae Mosquitoes, Southern Côte d’Ivoire

Preventing malaria used to seem as simple as killing the vector, the mosquito; however, a recent study shows that this concept is now anything but simple. The highly effective use of insecticide-treated bed nets and indoor insecticide spraying is being challenged by mosquito resistance to insecticides. In West Africa, populations of this mosquito vector are now resistant to all 4 classes of insecticide approved for this use. And no new classes of insecticide are anticipated until 2020, at the earliest. Development of newer classes of insecticide is crucial because if resistance continues unchecked, the hard-earned progress in malaria control in Africa could be quickly reversed.

Concordant Genetic Estimators of Migration Reveal Anthropogenically Enhanced Source-Sink Population Structure in the River Sculpin, Cottus gobio

River systems are vulnerable to natural and anthropogenic habitat fragmentation and will often harbor populations deviating markedly from simplified theoretical models. We investigated fine-scale population structure in the sedentary river fish Cottus gobio using microsatellites and compared migration estimates from three FST estimators, a coalescent maximum-likelihood method and Bayesian recent migration analyses. Source-sink structure was evident via asymmetry in migration and genetic diversity with smaller upstream locations emigration biased and larger downstream subpopulations immigration biased. Patterns of isolation by distance suggested that the system was largely, but not entirely, in migration-drift equilibrium, with headwater populations harboring a signal of past colonizations and in some cases also recent population bottlenecks. Up- vs. downstream asymmetry in population structure was partly attributable to the effects of flow direction, but was enhanced by weirs prohibiting compensatory upstream migration. Estimators of migration showed strong correspondence, at least in relative terms, especially if pairwise FST was used as an indirect index of relative gene flow rather than being translated to Nm. Since true parameter values are unknown in natural systems, comparisons among estimators are important, both to determine confidence in estimates of migration and to validate the performance of different methods.

Does kdr genotype predict insecticide-resistance phenotype in mosquitoes?

Several groups are developing and applying DNA-based technologies to monitor insecticide-based disease control programmes. However, several recent papers have concluded that the knockdown resistance (kdr) genotype-phenotype correlation that is observed in a wide variety of taxa might not hold in all mosquitoes. In this article, we review the evidence to support this putative breakdown and argue that the conclusion follows from unreliable data or the unparsimonious interpretation of data. We assert that the link between kdr genotype and DDT- and pyrethroid-susceptibility phenotype is clear. However, we emphasize that kdr genotype might explain only a portion of heritable variation in resistance and that diagnostic assays to test the importance of other resistance mechanisms in field populations are required.

Footprints of positive selection associated with a mutation (N1575Y) in the voltage-gated sodium channel of Anopheles gambiae

Insecticide resistance is an ideal model to study the emergence and spread of adaptative variants. In the African malaria mosquito, Anopheles gambiae, this is complemented by a strong public health rationale. In this insect, resistance to pyrethroid and DDT insecticides is strongly associated with the mutations L1014F and L1014S within the para voltage-gated sodium channel (VGSC). Across much of West Africa, 1014F frequency approaches fixation. Here, we document the emergence of a mutation, N1575Y, within the linker between domains III-IV of the VGSC. In data extending over 40 kbp of the VGSC 1575Y occurs on only a single long-range haplotype, also bearing 1014F. The 1014F-1575Y haplotype was found in both M and S molecular forms of An. gambiae in West/Central African sample sites separated by up to 2,000 km. In Burkina Faso M form, 1575Y allele frequency rose significantly from 0.053 to 0.172 between 2008 and 2010. Extended haplotype homozygosity analysis of the wild-type 1575N allele showed rapid decay of linkage disequilibrium (LD), in sharp contrast to the extended LD exhibited by 1575Y. A haplotype with long-range LD and high/increasing frequency is a classical sign of strong positive selection acting on a recent mutant. 1575Y occurs ubiquitously on a 1014F haplotypic background, suggesting that the N1575Y mutation compensates for deleterious fitness effects of 1014F and/or confers additional resistance to insecticides. Haplotypic tests of association suggest the latter: The 1014F-1575Y haplotype confers a significant additive benefit above 1014F-1575N for survival to DDT (M form P = 0.03) and permethrin (S form P = 0.003).

Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation

Adaptive introgression can provide novel genetic variation to fuel rapid evolutionary responses, though it may be counterbalanced by potential for detrimental disruption of the recipient genomic background. We examine the extent and impact of recent introgression of a strongly selected insecticide-resistance mutation (Vgsc-1014F) located within one of two exceptionally large genomic islands of divergence separating the Anopheles gambiae species pair. Here we show that transfer of the Vgsc mutation results in homogenization of the entire genomic island region (~1.5% of the genome) between species. Despite this massive disruption, introgression is clearly adaptive with a dramatic rise in frequency of Vgsc-1014F and no discernable impact on subsequent reproductive isolation between species. Our results show (1) how resilience of genomes to massive introgression can permit rapid adaptive response to anthropogenic selection and (2) that even extreme prominence of genomic islands of divergence can be an unreliable indicator of importance in speciation.

Field, Genetic, and Modeling Approaches Show Strong Positive Selection Acting upon an Insecticide Resistance Mutation in Anopheles gambiae s.s.

Alleles subject to strong, recent positive selection will be swept toward fixation together with contiguous sections of the genome. Whether the genomic signatures of such selection will be readily detectable in outbred wild populations is unclear. In this study, we employ haplotype diversity analysis to examine evidence for selective sweeps around knockdown resistance (kdr) mutations associated with resistance to dichlorodiphenyltrichloroethane and pyrethroid insecticides in the mosquito Anopheles gambiae. Both kdr mutations have significantly lower haplotype diversity than the wild-type (nonresistant) allele, with kdr L1014F showing the most pronounced footprint of selection. We complement these data with a time series of collections showing that the L1014F allele has increased in frequency from 0.05 to 0.54 in 5 years, consistent with a maximum likelihood-fitted selection coefficient of 0.16 and a dominance coefficient of 0.25. Our data show that strong, recent positive selective events, such as those caused by insecticide resistance, can be identified in wild insect populations.

Insecticide resistance monitoring of field-collected Anopheles gambiae s.l. populations from Jinja, eastern Uganda, identifies high levels of pyrethroid resistance

Insecticide resistance in the malaria vector Anopheles gambiae s.l. (Diptera: Culicidae) threatens insecticide‐based control efforts, necessitating regular monitoring. We assessed resistance in field‐collected An. gambiae s.l. from Jinja, Uganda using World Health Organization (WHO) biosassays. Only An. gambiae s.s. and An. arabiensis (≅70%) were present. Female An. gambiae exhibited extremely high pyrethroid resistance (permethrin LT50 > 2 h; deltamethrin LT50 > 5 h). Female An. arabiensis were resistant to permethrin and exhibited reduced susceptibility to deltamethrin. However, while An. gambiae were DDT resistant, An. arabiensis were fully susceptible. Both species were fully susceptible to bendiocarb and fenitrothion. Kdr 1014S has increased rapidly in the Jinja population of An. gambiae s.s. and now approaches fixation (≅95%), consistent with insecticide‐mediated selection, but is currently at a low frequency in An. arabiensis (0.07%). Kdr 1014F was also at a low frequency in An. gambiae. These frequencies preclude adequately‐powered tests for an association with phenotypic resistance. PBO synergist bioassays resulted in near complete recovery of pyrethroid susceptibility suggesting involvement of CYP450s in resistance. A small number (0.22%) of An. gambiae s.s. ×An. arabiensis hybrids were found, suggesting the possibility of introgression of resistance alleles between species. The high levels of pyrethroid resistance encountered in Jinja threaten to reduce the efficacy of vector control programmes which rely on pyrethroid‐impregnated bednets or indoor spraying of pyrethroids.

Association mapping of insecticide resistance in wild Anopheles gambiae populations: major variants identified in a low-linkage disequilbrium genome.

Background Association studies are a promising way to uncover the genetic basis of complex traits in wild populations. Data on population stratification, linkage disequilibrium and distribution of variant effect-sizes for different trait-types are required to predict study success but are lacking for most taxa. We quantified and investigated the impacts of these key variables in a large-scale association study of a strongly selected trait of medical importance: pyrethroid resistance in the African malaria vector Anopheles gambiae. Methodology/Principal Findings We genotyped ≈1500 resistance-phenotyped wild mosquitoes from Ghana and Cameroon using a 1536-SNP array enriched for candidate insecticide resistance gene SNPs. Three factors greatly impacted study power. (1) Population stratification, which was attributable to co-occurrence of molecular forms (M and S), and cryptic within-form stratification necessitating both a partitioned analysis and genomic control. (2) All SNPs of substantial effect (odds ratio, OR>2) were rare (minor allele frequency, MAF<0.05). (3) Linkage disequilibrium (LD) was very low throughout most of the genome. Nevertheless, locally high LD, consistent with a recent selective sweep, and uniformly high ORs in each subsample facilitated significant direct and indirect detection of the known insecticide target site mutation kdr L1014F (OR≈6; P<10−6), but with resistance level modified by local haplotypic background. Conclusion Primarily as a result of very low LD in wild A. Gambiae, LD-based association mapping is challenging, but is feasible at least for major effect variants, especially where LD is enhanced by selective sweeps. Such variants will be of greatest importance for predictive diagnostic screening.

Ecological Zones Rather Than Molecular Forms Predict Genetic Differentiation in the Malaria Vector Anopheles gambiae s.s. in Ghana

The malaria mosquito Anopheles gambiae s.s. is rapidly becoming a model for studies on the evolution of reproductive isolation. Debate has centered on the taxonomic status of two forms (denoted M and S) within the nominal taxon identified by point mutations in the X-linked rDNA region. Evidence is accumulating that there are significant barriers to gene flow between these forms, but that the barriers are not complete throughout the entire range of their distribution. We sampled populations from across Ghana and southern Burkina Faso, West Africa, from areas where the molecular forms occurred in both sympatry and allopatry. Neither Bayesian clustering methods nor FST-based analysis of microsatellite data found differentiation between the M and S molecular forms, but revealed strong differentiation among different ecological zones, irrespective of M/S status and with no detectable effect of geographical distance. Although no M/S hybrids were found in the samples, admixture analysis detected evidence of contemporary interform gene flow, arguably most pronounced in southern Ghana where forms occur sympatrically. Thus, in the sampled area of West Africa, lack of differentiation between M and S forms likely reflects substantial introgression, and ecological barriers appear to be of greater importance in restricting gene flow.