Insecticide Resistance Genes Research Articles

Complex genome evolution in Anopheles coluzzii associated with increased insecticide usage in Mali.

In certain cases, a species may have access to important genetic variation present in a related species via adaptive introgression. These novel alleles may interact with their new genetic background, resulting in unexpected phenotypes. In this study, we describe a selective sweep on standing variation on the X chromosome in the mosquito Anopheles coluzzii, a principal malaria vector in West Africa. This event may have been influenced by the recent adaptive introgression of the insecticide resistance gene known as kdr from the sister species Anopheles gambiae. Individuals carrying both kdr and a nearly fixed X‐linked haplotype, encompassing at least four genes including the P450 gene CYP9K1 and the cuticular protein CPR125, have rapidly increased in relative frequency. In parallel, a reproductively isolated insecticide‐susceptible A. gambiae population (Bamako form) has been driven to local extinction, likely due to strong selection from increased insecticide‐treated bed net usage.

In vivo functional analysis of the Drosophila melanogaster nicotinic acetylcholine receptor Dα6 using the insecticide spinosad

The vinegar fly, Drosophila melanogaster, has been used to identify and manipulate insecticide resistance genes. The advancement of genome engineering technology and the increasing availability of pest genome sequences has increased the predictive and diagnostic capacity of the Drosophila model. The Drosophila model can be extended to investigate the basic biology of the interaction between insecticides and the proteins they target. Recently we have developed an in vivo system that permits the expression and study of key insecticide targets, the nicotinic acetylcholine receptors (nAChRs), in controlled genetic backgrounds.Here this system is used to study the interaction between the insecticide spinosad and a nAChR subunit, Dα6. Reciprocal chimeric subunits were created from Dα6 and Dα7, a subunit that does not respond to spinosad. Using the in vivo system, the Dα6/Dα7 chimeric subunits were tested for their capacity to respond to spinosad. Only the subunits containing the C-terminal region of Dα6 were able to respond to spinosad, thus confirming the importance this region for spinosad binding.A new incompletely dominant, spinosad resistance mechanism that may evolve in pest species is also examined. First generated using chemical mutagenesis, the Dα6P146S mutation was recreated using the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 system, the first use of this technology to introduce a resistant mutation into a controlled genetic background. Both alleles present with the same incompletely dominant, spinosad resistance phenotype, proving the P146S replacement to be the causal mutation. The proximity of the P146S mutation to the conserved Cys-loop indicates that it may impair the gating of the receptor. The results of this study enhance the understanding of nAChR structure:function relationships.

Characterization of the transcriptome of the Asian gypsy moth Lymantria dispar identifies numerous transcripts associated with insecticide resistance

Although the Asian gypsy moth Lymantria dispar causes extensive forest damage worldwide, little is known regarding the genes involved in its development or response to insecticides. Accordingly, characterization of the transcriptome of L. dispar larvae would promote the development of toxicological methods for its control. RNA-seq analysis of L. dispar larvae messenger RNA (mRNA) generated 62,063 unigenes with N50 of 993 bp, from which 23,975 unique sequences (E-value < 10−5) were identified using a BLASTx search of the NCBI non-redundant (nr) database. Using functional classification in the Gene Ontology (GO) and Clusters of Orthologous Groups (COG) databases, 7,309 indentified sequences were categorized into 51 functional groups and 8,079 sequences were categorized into 25 functional groups, respectively. Moreover, we identified a large number of transcripts encoding known insecticide targets, or proteins involved in the metabolism of insecticides. Reads per kilobase of unigene length per million mapped reads (RPKM) analysis identified 39 high abundance transcripts, of which 27 exhibited significantly altered expression patterns across the egg, larvae, pupae, male and female adult stages. Our study provides the most comprehensive transcriptomic sequence resource for L. dispar, which will form the basis for future identification of candidate insecticide resistance genes in L. dispar.

Effect of isodillapiole on the expression of the insecticide resistance genes GSTE7 and CYP6N12 in Aedes aegypti from central Amazonia.

The yellow fever mosquito Aedes (Stegomyia) aegypti is the main vector of dengue arbovirus and other arboviruses. Dengue prevention measures for the control of A. aegypti involve mainly the use of synthetic insecticides. The constant use of insecticides has caused resistance in this mosquito. Alternative studies on plant extracts and their products have been conducted with the aim of controlling the spread of the mosquito. Dillapiole is a compound found in essential oils of the plant Piper aduncum (Piperaceae) which has been effective as a biopesticide against A. aegypti. Isodillapiole is a semisynthetic substance obtained by the isomerization of dillapiole. In the present study, isodillapiole was evaluated for its potential to induce differential expression of insecticide resistance genes (GSTE7 and CYP6N12) in 3rd instar larvae of A. aegypti. These larvae were exposed to this compound at two concentrations (20 and 40 μg/mL) for 4 h during four generations (G1, G2, G3, and G4). Quantitative RT-PCR was used to assess the expression of GSTE7 and CYP6N12 genes. GSTE7 and CYP6N12 relative expression levels were higher at 20 than at 40 μg/mL and varied among generations. The decrease in GSTE7 and CYP6N12 expression levels at the highest isodillapiole concentration suggests that larvae may have suffered from metabolic stress, revealing a potential alternative product in the control of A. aegypti.

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.

Genetic diversity of Brazilian Aedes aegypti: patterns following an eradication program.

Background Aedes aegypti is the most important vector of dengue fever in Brazil, where severe epidemics have recently taken place. Ae. aegypti in Brazil was the subject of an intense eradication program in the 1940s and 50s to control yellow fever. Brazil was the largest country declared free of this mosquito by the Pan-American Health Organization in 1958. Soon after relaxation of this program, Ae. aegypti reappeared in this country, and by the early 1980s dengue fever had been reported. The aim of this study is to analyze the present-day genetic patterns of Ae. aegypti populations in Brazil.Methodology/Principal FindingsWe studied the genetic variation in samples of 11 widely spread populations of Ae. aegypti in Brazil based on 12 well-established microsatellite loci. Our principal finding is that present-day Brazilian Ae. aegypti populations form two distinct groups, one in the northwest and one in the southeast of the country. These two groups have genetic affinities to northern South American countries and the Caribbean, respectively. This is consistent with what has been reported for other genetic markers such as mitochondrial DNA and allele frequencies at the insecticide resistance gene, kdr.Conclusions/SignificanceWe conclude that the genetic patterns in present day populations of Ae. aegypti in Brazil are more consistent with a complete eradication of the species in the recent past followed by re-colonization, rather than the alternative possibility of expansion from residual pockets of refugia. At least two colonizations are likely to have taken place, one from northern South American countries (e.g., Venezuela) that founded the northwestern group, and one from the Caribbean that founded the southeastern group. The proposed source areas were never declared free of Ae. aegypti.

De novo sequencing-based transcriptome and digital gene expression analysis reveals insecticide resistance-relevant genes in Propylaea japonica (Thunberg) (Coleoptea: Coccinellidae).

The ladybird Propylaea japonica (Thunberg) is one of most important natural enemies of aphids in China. This species is threatened by the extensive use of insecticides but genomics-based information on the molecular mechanisms underlying insecticide resistance is limited. Hence, we analyzed the transcriptome and expression profile data of P. japonica in order to gain a deeper understanding of insecticide resistance in ladybirds. We performed de novo assembly of a transcriptome using Illumina's Solexa sequencing technology and short reads. A total of 27,243,552 reads were generated. These were assembled into 81,458 contigs and 33,647 unigenes (6,862 clusters and 26,785 singletons). Of the unigenes, 23,965 (71.22%) have putative homologues in the non-redundant (nr) protein database from NCBI, using BLASTX, with a cut-off E-value of 10−5. We examined COG, GO and KEGG annotations to better understand the functions of these unigenes. Digital gene expression (DGE) libraries showed differences in gene expression profiles between two insecticide resistant strains. When compared with an insecticide susceptible profile, a total of 4,692 genes were significantly up- or down- regulated in a moderately resistant strain. Among these genes, 125 putative insecticide resistance genes were identified. To confirm the DGE results, 16 selected genes were validated using quantitative real time PCR (qRT-PCR). This study is the first to report genetic information on P. japonica and has greatly enriched the sequence data for ladybirds. The large number of gene sequences produced from the transcriptome and DGE sequencing will greatly improve our understanding of this important insect, at the molecular level, and could contribute to the in-depth research into insecticide resistance mechanisms.

Genetic structure of Bemisia tabaci Med populations from home‐range countries, inferred by nuclear and cytoplasmic markers: impact on the distribution of the insecticide resistance genes

Insecticide resistance management in Bemisia tabaci is one of the main issues facing agricultural production today. An extensive survey was undertaken in five Mediterranean countries to examine the resistance status of Med B. tabaci species in its range of geographic origin and the relationship between population genetic structure and the distribution of resistance genes. The investigation combined molecular diagnostic tests, sequence and microsatellite polymorphism studies and monitoring of endosymbionts. High frequencies of pyrethroid (L925I and T929V, VGSC gene) and organophosphate (F331W, ace1 gene) resistance mutations were found in France, Spain and Greece, but not in Morocco or Tunisia. Sequence analyses of the COI gene delineated two closely related mitochondrial groups (Q1 and Q2), which were found either sympatrically (Spain) or separately (France). Only Q1 was observed in Greece, Morocco and Tunisia. Bayesian analyses based on microsatellite loci revealed three geographically delineated genetic groups (France, Spain, Morocco/Greece/Tunisia) and high levels of genetic differentiation even between neighbouring samples. Evidence was also found for hybridisation and asymmetrical gene flow between Q1 and Q2. Med B. tabaci is more diverse and structured than reported so far. On a large geographic scale, resistance is affected by population genetic structure, whereas on a local scale, agricultural practices appear to play a major role.

Effect of Beauveria bassiana infection on detoxification enzyme transcription in pyrethroid resistant Anopheles arabiensis: a preliminary study

Fungal biopesticides are of great interest to vector control scientists as they provide a novel and environmentally friendly alternative to insecticide use. The aim of this study was to determine whether genes associated with pyrethroid resistance in Anopheles arabiensis from Sudan and South Africa are further induced following exposure to the entomopathogenic fungus, Beauveria bassiana (strain GHA). Following B. bassiana bioassays, RNA was extracted from infected mosquitoes and the transcription of four important insecticide resistance genes, CYP9L1, CYP6M2 and CYP4G16 (cytochrome P450s) and TPX4 (thioredoxin peroxidase) was investigated using quantitative real-time PCR. Beauveria bassiana strain GHA was highly infective and virulent against An. arabiensis. In terms of changes in gene transcription, overall, the fold change (FC) values for each gene in the infected strains, were lower than 1.5. The FC values of CYP9L1, CYP6M2 and TPX4, were significantly lower than the FC values of the same genes in uninfected resistant An. arabiensis. These data suggest that B. bassiana does not enhance the pyrethroid resistant phenotype on a molecular level as the two An. arabiensis strains used here, with different pyrethroid resistance mechanisms, revealed no increase in pre-existing metabolic transcripts. This supports the fact that fungal pathogens are suitable candidates for vector control, particularly with regard to the development of novel vector control strategies.

Demographic history and population structure of Anopheles pseudopunctipennis in Argentina based on the mitochondrial COI gene.

BackgroundAnopheles pseudopunctipennis is an important malaria vector in the Neotropical region and the only species involved in Plasmodium transmission in the Andean foothills. Its wide geographical distribution in America, high preference for biting humans and capacity to rest inside dwellings after feeding, are attributes contributing to its vector status. Previous reports have tried to elucidate its taxonomic status, distinguishing populations from North, Central and South America. In the present study we used a mitochondrial marker to examine the demographic history of An. pseudopunctipennis in northwestern Argentina.MethodsTwelve localities were selected across 550 km of the distribution of this species in Argentina, including two near the Bolivian border and several in South Tucumán, for sampling. A fragment of the cytochrome oxidase I (COI) gene was sequenced and haplotype relationships were analyzed by a statistical parsimony network and a Neighbor-Joining (NJ) tree. Genetic differentiation was estimated with FST. Historical demographic processes were evaluated using diversity measures, neutrality tests and mismatch distribution.ResultsForty-one haplotypes were identified, of which haplotype A was the most common and widely distributed. Neither the network nor the NJ tree showed any geographic differentiation between northern and southern populations. Haplotype diversities, Tajima’s DT and Fu & Li’s F and D neutrality tests and mismatch distribution supported a scenario of Holocene demographic expansion.ConclusionThe demographic pattern suggests that An. pseudopunctipennis has undergone a single colonization process, and the ancestral haplotype is shared by specimens from all localities, indicating mitochondrial gene flow. Genetic differentiation was minimal, observed only between one northern and one southern locality. The estimated time of the population expansion of this species was during the Holocene. These data suggest that regional vector control measures would be equally effective in both northern and southern localities sampled, but also that insecticide resistant genes may spread rapidly within this region.

Dissecting the organ specificity of insecticide resistance candidate genes in Anopheles gambiae: known and novel candidate genes.

BackgroundThe elevated expression of enzymes with insecticide metabolism activity can lead to high levels of insecticide resistance in the malaria vector, Anopheles gambiae. In this study, adult female mosquitoes from an insecticide susceptible and resistant strain were dissected into four different body parts. RNA from each of these samples was used in microarray analysis to determine the enrichment patterns of the key detoxification gene families within the mosquito and to identify additional candidate insecticide resistance genes that may have been overlooked in previous experiments on whole organisms.ResultsA general enrichment in the transcription of genes from the four major detoxification gene families (carboxylesterases, glutathione transferases, UDP glucornyltransferases and cytochrome P450s) was observed in the midgut and malpighian tubules. Yet the subset of P450 genes that have previously been implicated in insecticide resistance in An gambiae, show a surprisingly varied profile of tissue enrichment, confirmed by qPCR and, for three candidates, by immunostaining. A stringent selection process was used to define a list of 105 genes that are significantly (p ≤0.001) over expressed in body parts from the resistant versus susceptible strain. Over half of these, including all the cytochrome P450s on this list, were identified in previous whole organism comparisons between the strains, but several new candidates were detected, notably from comparisons of the transcriptomes from dissected abdomen integuments.ConclusionsThe use of RNA extracted from the whole organism to identify candidate insecticide resistance genes has a risk of missing candidates if key genes responsible for the phenotype have restricted expression within the body and/or are over expression only in certain tissues. However, as transcription of genes implicated in metabolic resistance to insecticides is not enriched in any one single organ, comparison of the transcriptome of individual dissected body parts cannot be recommended as a preferred means to identify new candidate insecticide resistant genes. Instead the rich data set on in vivo sites of transcription should be consulted when designing follow up qPCR validation steps, or for screening known candidates in field populations.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-1018) contains supplementary material, which is available to authorized users.

Pyrethroid Insecticide Resistance &amp;lt;i&amp;gt;kdr&amp;lt;/i&amp;gt; Gene in the House Fly, &amp;lt;i&amp;gt;Musca domestica&amp;lt;/i&amp;gt; (Diptera: Muscidae), in the United Arab Emirates

The house fly, Musca domestica L., is a very important insect because of its role as a vector of several human, livestock, and poultry diseases. Different groups of insecticides, including the pyrethroids, are used to control house flies. The objective of the current study was to conduct a survey of the kdr insecticide resistance gene in the house fly population in the United Arab Emirates (UAE). Flies were collected from five locations (Abu Dhabi, Dubai, Fujairah, Ajman, and Ras Al- Khaymah). Genomic DNA was extracted and polymerase chain reaction (PCR) amplification of specific allele (PASA) was performed. The maximum percentage (70%) of homozygous genotype (kdr/kdr) occurred at Ajman, followed by Dubai (59.5%), and the minimum percentage occurred at Ras Al-Khaymah (14.3%). Resistant insects of the genotype kdr/kdr were detected in all of the tested house fly populations. Ajman and Fujairah house fly populations were in Hardy-Weinberg equilibrium. The resistance kdr allele was found at a high frequency (0.54 - 0.8) at all locations except at Ras Al-Khaymah (0.21). Together, this study demonstrated that the pyrethroid insecticide resistance kdr allele was found in UAE house fly populations, strongly suggesting that a countrywide pyrethroid insecticide resistance management program needs to be implemented.

A comparative study of P450 gene expression in field and laboratory Musca domestica L. strains

The housefly is a global pest that has developed resistance to most insecticides applied for its control. Resistance has been associated with cytochrome P450 monooxygenases (P450s). The authors compare the expression of six genes possibly associated with insecticide resistance in three unselected strains: a multiresistant strain (791a), a neonicotinoid-resistant strain (766b) and a new field strain (845b). CYP4G2 was highly expressed throughout the range of strains and proved to be the one of the most interesting expression profiles of all P450s analysed. CYP6G4 was expressed up to 11-fold higher in 766b than in WHO-SRS. Significant differences between expression of P450 genes between F1 flies from 845b and established laboratory strains were shown. In general, P450 gene expression in 845b was 2-14-fold higher than in the reference strain (P < 0.0101) and 2-23-fold higher than in the multiresistant strain (P < 0.0110). The newly collected field strain 845b had significantly higher constitutive gene expression than both WHO-SRS and 791a. High constitutive expression of CYP4G2 in houseflies indicates a possible role of this gene in metabolic resistance. There is a strong indication that CYP6G4 is a major insecticide resistance gene involved in neonicotinoid resistance.

The impact of insecticide resistance on Culex pipiens immunity

Because of their role as vectors of diseases, the evolution of insecticide resistance in mosquitoes has been intensively investigated. Insecticide resistance is associated to a wide range of pleiotropic effects on several key life-history traits of mosquitoes such as longevity and behavior. However, despite its potential implications in pathogen transmission, the effects of insecticide resistance on mosquito immunity have received little, if any, attention. Here, we investigate the impact of insecticide resistance in Culex pipiens, an epidemiologically important vector of a wide array of pathogens. Using both isogenic laboratory strains and field-caught mosquitoes, we investigate the impact of two main insecticide resistance mechanisms (metabolic detoxification and target site modification) on the relative transcription of several genes involved in the immune response to pathogens, at both their constitutive and inducible levels. Our results show a discrepancy between the isogenic laboratory lines and field-collected mosquitoes: While in the isogenic strains, insecticide-resistant mosquitoes show a drastic increase in immune gene expression, no such effect appears in the field. We speculate on the different mechanisms that may underlie this discrepancy and discuss the risks of making inferences on the pleiotropic effects of insecticide-resistant genes by using laboratory-selected insecticide-resistant lines.

Identifying insecticide resistance genes in mosquito by combining AFLP genome scans and 454 pyrosequencing

AFLP-based genome scans are widely used to study the genetics of adaptation and to identify genomic regions potentially under selection. However, this approach usually fails to detect the actual genes or mutations targeted by selection owing to the difficulties of obtaining DNA sequences from AFLP fragments. Here, we combine classical AFLP outlier detection with 454 sequencing of AFLP fragments to obtain sequences from outlier loci. We applied this approach to the study of resistance to Bacillus thuringiensis israelensis (Bti) toxins in the dengue vector Aedes aegypti. A genome scan of Bti-resistant and Bti-susceptible A.aegypti laboratory strains was performed based on 432 AFLP markers. Fourteen outliers were detected using two different population genetic algorithms. Out of these, 11 were successfully sequenced. Three contained transposable elements (TEs) sequences, and the 10 outliers that could be mapped at a unique location in the reference genome were located on different supercontigs. One outlier was in the vicinity of a gene coding for an aminopeptidase potentially involved in Bti toxin-binding. Patterns of sequence variability of this gene showed significant deviation from neutrality in the resistant strain but not in the susceptible strain, even after taking into account the known demographic history of the selected strain. This gene is a promising candidate for future functional analysis.

The tomato borer Tuta absoluta : insecticide resistance and control failure.

Abstract The tomato borer Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae) evolved from a local South American concern to become a worldwide threat in tomato production less than three decades after initially spreading throughout South America, and eventually through Europe, North Africa and Middle East with an increased likelihood of invading Asia. Despite the studies performed thus far on this invasive species, management is difficult and insecticide is a major component of control. As a result, insecticide resistance has led to insecticide control failure in South America, where farmers have heavily relied on insecticide against this pest species since the 1980s. Reports from Chile, Argentina and, in particular, Brazil suggest that insecticide resistance evolves quickly in this species as a direct response to insecticide use (and overuse). Insecticide replacement for controlling the tomato borer typically leads to rapid changes in the patterns of insecticide resistance for this species. However, the resistance frequency for previously used compounds, such as pyrethroids, remains and likely results in a faster increase in resistance with selection, even in areas where this pest species was recently introduced. Recent genetic studies on the tomato borer suggest that insecticide resistance and control failure for this species in South America are also concerns elsewhere because newly introduced strains carry insecticide-resistance genes at a relatively high frequency even without local selection at the site of introduction.

Bionomics of sympatric chromosomal forms of Anopheles funestus (Diptera: Culicidae)

Anopheles funestus is one of the major vectors of malaria in Africa. Cytogenetic studies conducted on populations from West Africa have shown variable degrees of polymorphism with a genetic structure leading to the description of two chromosomal forms called "Folonzo" and "Kiribina" that exhibit limited gene flow. Because studies on allopatric populations showed bionomical heterogeneities, the present study was undertaken during three consecutive years (2006, 2007, and 2008) in an area of sympatry in Senegal, in order to assess their bionomical characteristics and compare their epidemiologic role in malaria transmission. Overall, the two forms coexisted in the study area; the Kiribina form being more abundant and exhibiting higher biting rates. Based on an enzyme-linked immunosorbent assay, the anthropophilic rates were statistically comparable and were, respectively, 30.7% and 28.6% for Kiribina and Folonzo. Plasmodium falciparum circumsporozoite rates were also comparable and were 2.7% for Kiribina and 3.1% for Folonzo. Both forms were involved in malaria transmission; Kiribina being responsible for 68% of transmission. Thus, due to the limited gene flow between the two forms, the introduction of transgene for Plasmodium resistance in one of the two forms could be a disadvantage for the implementation of control strategies based on the use of Plasmodium-refractory genetically modified individuals. Nevertheless, it could represent an advantage limiting the insurgence of insecticide resistance gene spread between forms and should be taken into account for the implementation of control strategies.

Population transcriptomics of Drosophila melanogaster females

BackgroundVariation at the level of gene expression is abundant in natural populations and is thought to contribute to the adaptive divergence of populations and species. Gene expression also differs considerably between males and females. Here we report a microarray analysis of gene expression variation among females of 16 Drosophila melanogaster strains derived from natural populations, including eight strains from the putative ancestral range in sub-Saharan Africa and eight strains from Europe. Gene expression variation among males of the same strains was reported previously.ResultsWe detected relatively low levels of expression polymorphism within populations, but much higher expression divergence between populations. A total of 569 genes showed a significant expression difference between the African and European populations at a false discovery rate of 5%. Genes with significant over-expression in Europe included the insecticide resistance gene Cyp6g1, as well as genes involved in proteolysis and olfaction. Genes with functions in carbohydrate metabolism and vision were significantly over-expressed in the African population. There was little overlap between genes expressed differently between populations in females and males.ConclusionsOur results suggest that adaptive changes in gene expression have accompanied the out-of-Africa migration of D. melanogaster. Comparison of female and male expression data indicates that the vast majority of genes differing in expression between populations do so in only one sex and suggests that most regulatory adaptation has been sex-specific.

Insecticide resistance and malaria transmission: infection rate and oocyst burden in Culex pipiens mosquitoes infected with Plasmodium relictum.

BackgroundThe control of most vectors of malaria is threatened by the spread of insecticide resistance. One factor that has been hitherto largely overlooked is the potential effects of insecticide resistance on the ability of mosquitoes to transmit malaria: are insecticide-resistant mosquitoes as good vectors of Plasmodium as susceptible ones? The drastic physiological changes that accompany the evolution of insecticide resistance may indeed alter the ability of vectors to transmit diseases, a possibility that, if confirmed, could have major epidemiological consequences.MethodsUsing a novel experimental system consisting of the avian malaria parasite (Plasmodium relictum) and its natural vector (the mosquito Culex pipiens), two of the most common mechanisms of insecticide resistance (esterase overproduction and acetylcholinesterase modification) were investigated for their effect on mosquito infection rate and parasite burden. For this purpose two types of experiments were carried out using (i) insecticide-resistant and susceptible laboratory isogenic lines of Cx. pipiens and (ii) wild Cx. pipiens collected from a population where insecticide resistant and susceptible mosquitoes coexist in sympatry.ResultsThe isogenic line and wild-caught mosquito experiments were highly consistent in showing no effect of either esterase overproduction or of acetylcholinesterase modification on either the infection rate or on the oocyst burden of mosquitoes. The only determinant of these traits was blood meal size, which was similar across the different insecticide resistant categories in both experiments.ConclusionsInsecticide resistance was found to have no effect on Plasmodium development within the mosquito. This is the first time this question has been addressed using a natural mosquito-Plasmodium combination, while taking care to standardize the genetic background against which the insecticide resistance genes operate. Infection rate and oocyst burden are but two of the factors that determine the vectorial capacity of mosquitoes. Other key determinants of parasite transmission, such as mosquito longevity and behaviour, or the parasite's incubation time, need to be investigated before concluding on whether insecticide resistance influences the ability of mosquitoes to transmit malaria.

Association Mapping of Insecticide Resistance in Wild Anopheles gambiae Populations: Major Variants Identified in a Low-Linkage Disequilbrium Genome

BackgroundAssociation 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 FindingsWe 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.ConclusionPrimarily 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.