Keith Steen

Mapping a Quantitative Trait Locus (QTL) conferring pyrethroid resistance in the African malaria vector Anopheles funestus

BackgroundPyrethroid resistance in Anopheles funestus populations has led to an increase in malaria transmission in southern Africa. Resistance has been attributed to elevated activities of cytochrome P450s but the molecular basis underlying this metabolic resistance is unknown. Microsatellite and SNP markers were used to construct a linkage map and to detect a quantitative trait locus (QTL) associated with pyrethroid resistance in the FUMOZ-R strain of An. funestus from Mozambique.ResultsBy genotyping 349 F2 individuals from 11 independent families, a single major QTL, rp1, at the telomeric end of chromosome 2R was identified. The rp1 QTL appears to present a major effect since it accounts for more than 60% of the variance in susceptibility to permethrin. This QTL has a strong additive genetic effect with respect to susceptibility. Candidate genes associated with pyrethroid resistance in other species were physically mapped to An. funestus polytene chromosomes. This showed that rp1 is genetically linked to a cluster of CYP6 cytochrome P450 genes located on division 9 of chromosome 2R and confirmed earlier reports that pyrethroid resistance in this strain is not associated with target site mutations (knockdown resistance).ConclusionWe hypothesize that one or more of these CYP6 P450s clustered on chromosome 2R confers pyrethroid resistance in the FUMOZ-R strain of An. funestus.

Differential expression of the detoxification genes in the different life stages of the malaria vector Anopheles gambiae.

The diverse habitats and diets encountered during the life cycle of an Anopheles mosquito have necessitated the development of extensive families of detoxification enzymes. Expansion of the three detoxification enzyme families (cytochrome P450s, carboxylesterases and glutathione transfereases), has occurred in mosquitoes compared with Drosophila, however, very little is known regarding the developmental expression of theses genes. Using a custom made microarray we determined the expression profile of the detoxification genes in adults, larvae and pupae of the malaria vector A. gambiae. The expression of approximately one quarter of these genes was developmentally regulated. The expression profile of each of these genes and the information this data provides on putative functions of the mosquito detoxification enzymes is discussed.

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.

High, clustered, nucleotide diversity in the genome of Anopheles gambiae revealed through pooled-template sequencing: implications for high-throughput genotyping protocols

BackgroundAssociation mapping approaches are dependent upon discovery and validation of single nucleotide polymorphisms (SNPs). To further association studies in Anopheles gambiae we conducted a major resequencing programme, primarily targeting regions within or close to candidate genes for insecticide resistance.ResultsUsing two pools of mosquito template DNA we sequenced over 300 kbp across 660 distinct amplicons of the An. gambiae genome. Comparison of SNPs identified from pooled templates with those from individual sequences revealed a very low false positive rate. False negative rates were much higher and mostly resulted from SNPs with a low minor allele frequency. Pooled-template sequencing also provided good estimates of SNP allele frequencies. Allele frequency estimation success, along with false positive and negative call rates, improved significantly when using a qualitative measure of SNP call quality. We identified a total of 7062 polymorphic features comprising 6995 SNPs and 67 indels, with, on average, a SNP every 34 bp; a high rate of polymorphism that is comparable to other studies of mosquitoes. SNPs were significantly more frequent in members of the cytochrome p450 mono-oxygenases and carboxy/cholinesterase gene-families than in glutathione-S-transferases, other detoxification genes, and control genomic regions. Polymorphic sites showed a significantly clustered distribution, but the degree of SNP clustering (independent of SNP frequency) did not vary among gene families, suggesting that clustering of polymorphisms is a general property of the An. gambiae genome.ConclusionThe high frequency and clustering of SNPs has important ramifications for the design of high-throughput genotyping assays based on allele specific primer extension or probe hybridisation. We illustrate these issues in the context of the design of Illumina GoldenGate assays.

Contemporary evolution of resistance at the major insecticide target site gene Ace-1 by mutation and copy number variation in the malaria mosquito Anopheles gambiae

Functionally constrained genes are ideal insecticide targets because disruption is often fatal, and resistance mutations are typically costly. Synaptic acetylcholinesterase (AChE) is an essential neurotransmission enzyme targeted by insecticides used increasingly in malaria control. In Anopheles and Culex mosquitoes, a glycine–serine substitution at codon 119 of the Ace‐1 gene confers both resistance and fitness costs, especially for 119S/S homozygotes. G119S in Anopheles gambiae from Accra (Ghana) is strongly associated with resistance, and, despite expectations of cost, resistant 119S alleles are increasing significantly in frequency. Sequencing of Accra females detected only a single Ace‐1 119S haplotype, whereas 119G diversity was high overall but very low at non‐synonymous sites, evidence of strong purifying selection driven by functional constraint. Flanking microsatellites showed reduced diversity, elevated linkage disequilibrium and high differentiation of 119S, relative to 119G homozygotes across up to two megabases of the genome. Yet these signals of selection were inconsistent and sometimes weak tens of kilobases from Ace‐1. This unexpected finding is attributable to apparently ubiquitous amplification of 119S alleles as part of a large copy number variant (CNV) far exceeding the size of the Ace‐1 gene, whereas 119G alleles were unduplicated. Ace‐1 CNV was detectable in archived samples collected when the 119S allele was rare in Ghana. Multicopy amplification of resistant alleles has not been observed previously and is likely to underpin the recent increase in 119S frequency. The large CNV compromised localization of the strong selective sweep around Ace‐1, emphasizing the need to integrate CNV analysis into genome scans for selection.

An Integrated Genetic and Physical Map for the Malaria Vector Anopheles funestus

We have constructed a genetic map of the major African malaria vector, Anopheles funestus, using genetic markers segregating in F2 progeny from crosses between two strains colonized from different field sites. Genotyping was performed on 174 progeny from three families using 33 microsatellite markers, a single RFLP, and 15 single nucleotide polymorphism (SNP) loci. Four linkage groups were resolved and these were anchored to chromosomes X and 2 and chromosomal arms 3R and 3L by comparison with a physical map of this species. Five markers were linked to the X chromosome, 16 markers to chromosome 2, and 10 and 11 markers to chromosomal arms 3R and 3L, respectively. This significantly increases the number of chromosomally defined genetic markers for this species and will facilitate the identification of genes controlling epidemiologically important traits such as resistance to insecticides or vector competence.

Reduced susceptibility to DDT in field populations of Anopheles quadriannulatus and Anopheles arabiensis in Malawi: evidence for larval selection

Abstract Bioassays for insecticide resistance in adult mosquitoes were conducted on samples of Anopheles gambiae Giles s.l. (Diptera: Culicidae) species collected as larvae from breeding sites in the lower Shire Valley, Malawi. The results indicate full susceptibility to permethrin, deltamethrin and malathion, but reduced susceptibility to DDT in one sample from Thom (LT50 of 8.39 min for females and 25.09 min for males). Polymerase chain reaction‐based species identification of the mosquitoes assayed revealed a mixture of Anopheles arabiensis Patton and Anopheles quadriannulatus (Theobold). The LT50 did not differ significantly between species. Genotyping of the L1014F and L1014S kdr alleles showed all mosquito specimens to be homozygous wild type; thus the reduced susceptibility detected is not attributable to target site insensitivity and instead is likely to be metabolic in nature. Anopheles quadriannulatus is characteristically zoophagic and exophilic. Indeed, of 82 Anopheles collected through knockdown collections within dwellings, only one was An. quadriannulatus and the rest were An. arabiensis. They are unlikely, therefore, to have been exposed to selection pressure arising from insecticide‐treated net usage or to DDT indoor residual spraying. Therefore, it is suggested that this example of reduced susceptibility to DDT in An. quadriannulatus reflects selection in the larval stages.

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 basis of pyrethroid resistance in a population of Anopheles arabiensis, the primary malaria vector in Lower Moshi, north-eastern Tanzania

BackgroundPyrethroid resistance has been slower to emerge in Anopheles arabiensis than in An. gambiae s.s and An. funestus and, consequently, studies are only just beginning to unravel the genes involved. Permethrin resistance in An. arabiensis in Lower Moshi, Tanzania has been linked to elevated levels of both P450 monooxygenases and β-esterases. We have conducted a gene expression study to identify specific genes linked with metabolic resistance in the Lower Moshi An. arabiensis population.MethodsMicroarray experiments employing an An. gambiae whole genome expression chip were performed on An. arabiensis, using interwoven loop designs. Permethrin-exposed survivors were compared to three separate unexposed mosquitoes from the same or a nearby population. A subsection of detoxification genes were chosen for subsequent quantitative real-time PCR (qRT-PCR).ResultsMicroarray analysis revealed significant over expression of 87 probes and under expression of 85 probes (in pairwise comparisons between permethrin survivors and unexposed sympatric and allopatric samples from Dar es Salaam (controls). For qRT-PCR we targeted over expressed ABC transporter genes (ABC ‘2060’), a glutathione-S-transferase, P450s and esterases. Design of efficient, specific primers was successful for ABC ‘2060’and two P450s (CYP6P3, CYP6M2). For the CYP4G16 gene, we used the primers that were previously used in a microarray study of An. arabiensis from Zanzibar islands. Over expression of CYP4G16 and ABC ‘2060’ was detected though with contrasting patterns in pairwise comparisons between survivors and controls. CYP4G16 was only up regulated in survivors, whereas ABC ‘2060’ was similar in survivors and controls but over expressed in Lower Moshi samples compared to the Dar es Salaam samples. Increased transcription of CYP4G16 and ABC ‘2060’ are linked directly and indirectly respectively, with permethrin resistance in Lower Moshi An. arabiensis.ConclusionsIncreased transcription of a P450 (CYP4G16) and an ABC transporter (ABC 2060) are linked directly and indirectly respectively, with permethrin resistance in Lower Moshi An. arabiensis. Our study provides replication of CYP4G16 as a candidate gene for pyrethroid resistance in An. arabiensis, although its role may not be in detoxification, and requires further investigation.

Contemporary gene flow between wild An. gambiae s.s. and An. arabiensis.

BackgroundIn areas where the morphologically indistinguishable malaria mosquitoes Anopheles gambiae Giles and An. arabiensis Patton are sympatric, hybrids are detected occasionally via species-diagnostic molecular assays. An. gambiae and An. arabiensis exhibit both pre- and post-reproductive mating barriers, with swarms largely species-specific and male F1 (first-generation) hybrids sterile. Consequently advanced-stage hybrids (back-crosses to parental species), which would represent a route for potentially-adaptive introgression, are expected to be very rare in natural populations. Yet the use of one or two physically linked single-locus diagnostic assays renders them indistinguishable from F1 hybrids and levels of interspecific gene flow are unknown.MethodsWe used data from over 350 polymorphic autosomal SNPs to investigate post F1 gene flow via patterns of genomic admixture between An. gambiae and An. arabiensis from eastern Uganda. Simulations were used to investigate the statistical power to detect hybrids with different levels of crossing and to identify the hybrid category significantly admixed genotypes could represent.ResultsA range of admixture proportions were detected for 11 field-collected hybrids identified via single-locus species-diagnostic PCRs. Comparison of admixture data with simulations indicated that at least seven of these hybrids were advanced generation crosses, with backcrosses to each species identified. In addition, of 36 individuals typing as An. gambiae or An. arabiensis that exhibited outlying admixture proportions, ten were identified as significantly mixed backcrosses, and at least four of these were second or third generation crosses.ConclusionsOur results show that hybrids detected using standard diagnostics will often be hybrid generations beyond F1, and that in our study area around 5% (95% confidence intervals 3%-9%) of apparently ‘pure’ species samples may also be backcrosses. This is likely an underestimate because of rapidly-declining detection power beyond the first two backcross generations. Post-F1 gene flow occurs at a far from inconsequential rate between An. gambiae and An. arabiensis, and, especially for traits under strong selection, could readily lead to adaptive introgression of genetic variants relevant for vector control.