- Research article
- Open Access
Dissecting the organ specificity of insecticide resistance candidate genes in Anopheles gambiae: known and novel candidate genes
BMC Genomics volume 15, Article number: 1018 (2014)
The 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.
A 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.
The 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.
Insecticides play a vital role in reducing malaria transmission in Africa. An escalation in the use of two very effective tools, indoor residual spraying with insecticides and insecticide treated bednets, has led to impressive reductions in malaria with child death rates halved and more than 3.3 million lives saved since 2000 . Inevitably, as insecticide use has intensified, malaria vectors have developed resistance to these chemicals [2–4]. With just four classes of insecticides available for public health and only the pyrethroids approved for bednet treatment, this poses a major challenge to sustaining and extending recent achievements in malaria reduction.
Numerous studies have attempted to identify the genes responsible for insecticide resistance in the major malaria vectors. One of the most potent mechanisms identified to date is increased activity of enzymes that detoxify insecticides [4–6]. Four enzyme families are known to be associated with insecticide metabolism (carboxylesterases (CCEs), glutathione transferases (GSTs), UDP glucornyltransferases (UGTs) and cytochrome P450s (P450s)) and a number of individual enzymes, most notably from the cytochrome P450 family, have been implicated in conferring resistance to one or more insecticide classes [2, 7–9]. More recently, the importance of interactions between different enzymes and transporters in the insecticide detoxification pathway has been recognised [10, 11]. To dissect these pathways further, and to distinguish members of these large gene families with housekeeping functions from those more likely to have detoxification roles, further information on their sites of expression is required. The first objective of this study was to characterise expression patterns of the key gene families associated with insecticide resistance across the major organs linked to xenobiotic detoxification in insects, with a particular focus on the P450s. Although transcriptomes of many of the key tissues in Anopheles gambiae have already been described [12–15], this study used material from a highly insecticide resistant strain and from a susceptible strain to identify genes whose tissue-specific enrichment might be linked to the resistance phenotype.
To date, all comparisons of the transcriptome between insecticide resistant and susceptible malaria vectors have compared gene expression in the whole organism. This approach has the potential to miss candidates. If, for example, expression of a gene is restricted to an organ that contributes only a small proportion of mRNA to the total RNA pool, or differential expression occurs in only one tissue, even large differences in expression between a resistant and susceptible population may not be detectable . Thus the second objective of the study was to compare gene expression in key body parts between an insecticide resistant and susceptible strain of mosquito to identify candidates not immediately apparent in whole organism microarray studies. Adult mosquitoes were dissected into body parts that could be readily separated with minimal risk of contamination, and are suspected to be involved in metabolic resistance.
Results and discussion
RNA was extracted from three dissected ‘body parts’: the malpighian tubules, the midgut and the abdomen integument (containing the fat body, but also epidermal, neuronal, muscle and oenocyte cells) with the remaining undissected body parts forming a fourth sample group. Each biological replicate consisted of 15–20 adult female mosquitoes from the major African malaria vector, An. gambiae. Dissections were performed on an insecticide susceptible strain (N’Gousso) originally from Cameroon  and the Tiassalé strain from Cote d’Ivoire, which is highly resistant to all four classes of insecticide [2, 3].
Transcription in these body parts was compared in two ways (i) each body part against the corresponding whole organism for both strains and (ii) resistant against the corresponding susceptible body parts (Additional file 1: Figure S1).
Gene enrichment in individual body part versus whole organism
Transcripts showing enriched transcription in each of the four body parts were determined using a multiple test correction significance cut-off of p ≤0.05 for both the resistant and susceptible strains. As expected, a clear positive correlation can be observed for local transcription between the resistant and susceptible strains in each body part. This is illustrated in Figure 1, where the vast majority of probes follow a y = x trend.
In all body parts, similar numbers of transcripts showed enriched or depleted transcription compared to the whole organism, with a range of 1.4% to 4.4% of the total probes (Additional file 2: Table S1). The magnitude of change in transcription of individual transcripts in the midgut and remaining body parts is relatively low with no transcripts exceeding log2 fold change of 2, compared to the abdomen integument and malpighian tubules where 7 and 42 genes respectively are above log2 3-fold enriched (Figure 1). The full lists of transcripts enriched in each body parts are listed in Additional file 3: Table S2.
The An. gambiae genome contains 211 genes encoding four major detoxification families (111 cytochrome P450s, 31 GSTs, 43 CCES (only putatively catalytically active enzymes included) and 26 UGTs), together comprising 1.5% of the probes on the Agilent array . Each of the dissected body parts have an over representation of detoxification transcripts with more members of these gene families overtranscribed relative to the whole organism than found in the ‘remaining’ (undissected body part). Additional file 2: Table S1 shows that 8.5% of the total number of detoxification transcripts are enriched in the malpighian tubules of the resistant strain, 13.7% in the midgut and 7.6% in the abdomen integument. In contrast, the ‘remaining’ (undissected body part) show a depletion of detoxification transcripts in the overtranscribed subset when compared to the whole organism. These data reinforce the importance of the selected body parts in xenobiotic detoxification.
The cytochrome P450 family has been most strongly linked with insecticide resistance in Anopheles mosquitoes, with several enzymes capable of detoxifying insecticides from more than one class [2, 8]. Identifying the primary sites of transcription of this enzyme family will aid prediction of function  and help identify the key organs largely responsible for insecticide detoxification in resistant mosquitoes. Within the P450 gene family, body part enrichment shows some relationship with the gene tree clustering (Additional file 4: Figure S2) with the CYP9J and 6P families largely enriched in the midgut, the CYP4Gs, 6Ys and 325Cs enriched in the abdomen integument and the CYP6Z family enriched in the malpighian tubules. The diversity in enrichment patterns within this gene family led us to look specifically at thirteen cytochrome P450s that had been implicated in insecticide resistance to see if these were enriched in a particular body part (Figure 2). The criteria for inclusion of these P450s as candidate insecticide resistant genes was that they had been found to be significantly over expressed in pyrethroid and/or DDT resistant An. gambiae populations in more than one independent study. Four of these transcripts were significantly (p < 0.05) and highly (Log2 Fold Change >1.5) enriched in the malpighian tubules (CYP6M3, CYP6Z1, CYP6Z2 and CYP6Z3) compared to whole organism and one was enriched in the midgut and two in the abdomen integument (CYP4H24 and CYP4G16/CYP4G17 respectively). The remaining 5 showed no significant tissue enrichment.
The enrichment of a subset of seven cytochrome P450s in particular body parts was confirmed by qPCR, although for two of these midgut enriched transcripts, a much greater over expression was observed by qPCR than for microarray (Additional file 5: Figure S3) (A level of discordance between qPCR and microarray is frequently observed  (but it in the current study, the direction of change was consistent between the two methodologies for all seven genes). Antibodies were available for three cytochrome P450s, which were used to confirm the major sites of transcription within the abdomen integument of the Tiassalé resistant strain. In agreement with microarray transcription data, CYP6Z1 and CYP6Z2 were detected in the malpighian tubules of resistant mosquitoes, whilst CYP4G17, identified as enriched in the abdomen integument by both microarray and qPCR, was found only in the oenocytes, a cellular layer located under the abdomen cuticle (Figure 3). A CYP4G from Drosophila melanogaster, CYP4G1, has also been shown to be highly enriched in oenocyte cells; this Drosophila enzyme catalyses a key step in the formation of cuticular hydrocarbons . The potential role of over expression of 4G17, and its paralogue 4G16, in altering the cuticular structure in insecticide resistant mosquitoes is currently under investigation.
A new list of candidate insecticide resistance genes from the body part specific arrays
The two strains used in this study originate from sites separated by approximately 2,500 km, so substantial variation in gene transcription between them is to be expected, regardless of the difference in their insecticide resistance profile. Thus caution must be applied when correlating gene transcription levels with the resistance phenotype. However, with a goal of identifying further candidate insecticide resistance genes from the direct comparisons between dissected body parts from the resistant and susceptible strain for further functional validation, we applied a stringent selection process to derive a new gene list. This involved selecting only genes detected as enriched in body parts from Tiassalé versus N’Gousso via both the limma and GaGa methods with an adjusted p value ≤0.001 and setting a cut off of 1.4-fold differential expression (see Methods section). This identified a list of 134 transcripts, representing 105 genes transcribed at higher levels in the Tiassalé strain and 16 genes with higher transcription in the susceptible N’Gousso strain (Additional file 6: Table S3). Eleven of these transcripts were selected for qPCR validation yielding a positive correlation with the array (Pearson’s correlation (r = 0.870)) (Figure 4).
All of the 9 cytochrome P450s (CYP6P2, P3 and P4, CYP6Z2, Z3 CYP4H17, CYP4K2, CYP4C35 and CYP4G16) on this candidate list were also detected in the direct whole comparisons of the whole organism transcriptomes between these strains [3 and C. Strode, unpublished data]. Indeed the majority (89) of the 134 transcripts on our candidate list derived from the comparisons of dissected body parts were also detected in comparisons of the entire transcriptome between the two strains. Nevertheless, 22 (16.4%) would have been missed using whole organism arrays and 23 (17.2%) of the genes are regulated in the opposite direction between the whole organism and dissected body part comparisons. Of the 22 transcripts not detected in strain comparisons at the whole organism level, the majority (13) were detected from the abdomen integument (Table 1). Abdomen integuments from other resistant strains of An. gambiae are being dissected to search for further supporting evidence for a role of these transcripts in conferring resistance prior to follow up functional analysis.
Microarrays are widely used to identify insecticide resistance mechanisms in mosquito populations [3, 9, 21, 22]. However, all of these studies have used whole-organism RNA as a template, which may dilute or distort the final level of transcript expression detectable. By dissecting some of the major body parts involved in xenobiotic detoxification from two strains of An. gambiae, differing in their resistance phenotype we have been able to simultaneously identify the local expression profiles of known insecticide resistance candidates and compare transcription between the two strains within individual body parts.
This rich data set will be useful for establishing pathways of detoxification as genes catalysing the three classic classes of drug metabolism (oxidation, conjugation and excretion)  would be expected to be co-regulated. Pyrethroid mimetic activity–based probes, used to detect pyrethroid metabolising enzymes in the rat liver, identified a potential network of drug metabolising enzymes from multiple families involved in pyrethroid metabolism . Applying the same approach to insecticide resistant mosquitoes and using the data on local transcription from the current data set, will help unravel the pathways of insecticide metabolism selected for by intensive use of pyrethroids.
Although potential new insecticide resistance candidates have emerged from this study, it is encouraging that the majority of candidate insecticide resistant transcripts identified from direct comparison of the transcriptomes of dissected body part were also detected in the whole organisms comparisons. No single body part emerged as the key site of overtranscription of putative insecticide resistance genes in this study and it is therefore recommended that, unless resources enable a more comprehensive study design involving multiple dissected tissues, transcriptional approaches to identify candidate insecticide resistance transcripts continue to use the whole body transcriptome. Nevertheless this data set on local sites of transcription should be consulted when designing follow up qPCR validation steps, or for screening known candidates in field populations.
Mosquito rearing conditions
The An. gambiae used in these experiments were all reared under standard insectary conditions at 27°C and 70-80% humidity under a 16:8 hour photoperiod. The N’Gousso strain is originally from Cameroon and is susceptible to all classes of insecticide . N’Guosso is the M molecular form of A gambiae, recently re-classified as a separate species, Ano coluzzi. In contrast, the Tiassalé strain from Côte D’Ivoire is resistant to all classes of insecticide [2, 3]. This strain was colonised from the field site in 2012 and is a mixture of the M and S molecular forms. At the time of the study, the LD50 for the Tiassalé strain was 68 and 81 fold higher than the corresponding value for the N’Gousso strain for permethrin and deltamethrin respectively. Further details of the resistance profile of this strain are contained within references by Edi et al. [2, 3].
RNA was extracted from three dissected body parts: the malpighian tubules, the midgut and the abdomen integument (containing the fat body but also epidermal, neuronal, muscle and oenocyte cells) with the remaining undissected body parts forming a fourth sample group. Mosquitoes were collected between the hours of 8AM and 2 pm and dissected immediately on a CO2 block. Post dissection, each body part was added to extraction buffer from the PicoPure RNA extraction kit, heated for 30 minutes at 42°C and frozen at -80°C as per manufacturers instructions. Each biological replicate for each strain consisted of RNA, extracted using PicoPure RNA Isolation kit (Arcturus), from 12 3–5 day old non-blood fed, presumed mated females. The quantity and quality of the RNA was assessed using a nanodrop spectrophotometer (Nanodrop Technologies UK) and Bioanalyser (Agilent) respectively. Four biological replicates were prepared for each body parts per strain. RNA from the four dissections was pooled according to the proportion of RNA extracted from each body part to reconstitute the ‘whole organism’ (7%, 6%, 24% and 63% RNA from abdomen integument, malpighian tubules, midgut, and remaining material respectively). The use of a reconstituted reference sample minimised potential sources of bias that could have arisen from circadian changes in gene expression and changes in the proportion of the M or S molecular form in the different biological replicates. 100 ng of RNA was amplified and labelled with Cy3 and Cy5, using the Two colour low input Quick Amp labelling kit (Agilent) following the manufacturers instructions. Samples were then purified (Qiagen) with the cRNA yield and quality assessed using the nano-drop and Bioanalyser respectively. RNA from each Tiassalé body part was competitively hybridised with the respective N’Gousso body part, as well as each body part from the resistant and susceptible strain being compared to the re-constituted whole organism (Additional file 1: Figure S1). Dye swaps were performed on two out of four technical replicates for each array, to correct for dye bias.
Labelled cRNAs were hybridised to the whole genome 8×15k An gambiae array (ArrayExpress accession number A-MEXP-2211). Microarray hybridisation, washing and scanning were performed according to previously described protocols .
The resulting data were analysed using R. Within-array normalisation was carried out by Loess, and between array normalisation by Aquantile. Signals were corrected for dye automatically. The limma package  was used to fit linear models to the normalised data. In the case of complete loop designs a design matrix was used to infer the appropriate contrast matrices for each array. All parameters used were default. A bespoke pipeline using the GaGa package  was used to fit gamma-gamma models of variation to normalised corrected signals, in order to assign probes to one of two patterns of expression X equals Y or X does not equal Y, where X represents the resistant population arrays and Y the susceptible arrays. These data were subsequently used to assess enrichment in each expression pattern, through GO term analysis using the TopGO package . A standard FDR adjusted p value cut off of p ≤0.05 was applied to all data describing localisation of detoxification candidates. A second stringent selection method was used to reduce the probe list based on previously published methodology , requiring that the following criteria were met: adjusted p-value ≤0.001, raw fluorescence intensity > median, and Tiassalé vs. N’Gousso ±0.485 Log2 fold change between the strains. All candidates selected also demonstrated a positive GaGa analysis fold change, indicative of higher transcript localisation in the resistant tissue to the susceptible, thereby utilising all available array data.
RNA (4 μg) from each biological replicate was reverse transcribed using Oligo dT (Invitrogen) and Superscript III (Invitrogen) according to manufacturers instructions. Quantitative real-time PCR was performed using SYBR Green Supermix III (Applied Biosystems) using an MX3005 and the associated MxPro software (Agilent). Primer Blast (NCBI)  was used to design primer pairs (Additional file 7: Table S4). Where possible, primers were designed to span an exon junction but this was not possible for six of the P450 genes (CYP325A1, CYP6P3, CYP4G17, CYP6Z3, CYP12F2 and CYP6Z2) due to the high degree of polymorphisms in their DNA sequence. Each 20 μl reaction contained 10 μl SYBR Green Supermix, 0.3 μM of each primer and 1 μl of 1:10 diluted cDNA. Standard curves were produced using whole N’Gousso cDNA, in 1, 1:5, 1:25, 1:125 dilutions, (48.3 ng/μl to 0.386 ng/μl). qPCR was performed with the following conditions: 3 minutes at 95°C, with 40 cycles of 10 seconds at 95°C and 10 seconds at 60°C. All amplification efficiencies of designed primers were within acceptable range (90-120%), following MIQE guidelines .
Preparation of antibodies
Fragments encoding unique peptides for CYP6Z1 and CYP4G17 were cloned into the pET 16b vector. Upon expression, the resulting His-tagged peptide was purified to homogeneity by Ni-NTA affinity chromatography and used to raise rabbit polyclonal antibodies. The CYP6Z1 peptide sequence was: VALRDLNNPDSFINNIRTAGVFLCPGLLKFTGINSLSPPMKKFTTEVISSHLHQRETGQVTRKDFIQMLTDLRRKAGSSGEETLTDA and the CYP4G17 peptide: KRQLKIHLRLDPLFNLTGVKKEQERLLQIIHGLTRKVVREKKQLYERQMAEGKMPSPSLTEIIGKEEKPGEGQLGGSPAFISQ. The antibody for CYP6Z2 was a gift from Dr Mark Paine (LSTM, UK). Rabbit antibodies to CYP6Z2 were prepared to the C terminal peptide sequence MRIDHRK by Moravian Biotechnology, Brno, Czech Republic.
Immunofluorescence and microscopy
Female mosquitoes (3–5 days old) were fixed in cold solution of 4% PFA (methanol free, Thermo scientific) in phosphate-buffered saline (PBS) for 4 h and then were cryo-protected in 30% sucrose/PBS at 4°C for 12 h. Finally, mosquitoes were immobilized in O.C.T. (Tissue-Tek, SAKURA) and stored at -80°C until use.
Immunofluorescent analysis, followed by confocal microscopy, was performed to longitudinal sections of frozen pre-fixed mosquito specimens. More detailed, 10 mm sections, obtained in Leitz kryostat 1720 digital, were washed (3 × 5 min) with 0,05% Tween in PBS and blocked for 3 h in blocking solution (1% Fetal Bovine Serum, biosera, in 0,05% Triton/PBS). Then, the sections were stained with rabbit primary antibodies in 1/500 dilution, followed by goat anti-rabbit (Alexa Fluor 488, Molecular Probes) (1/1000) that gave the green colour. Also To-PRO 3-Iodide (Molecular Probes), which stains DNA specifically (red colour), was used, after RNAse A treatment. As controls, pre-immune serums (in 1/500 dilutions) and anti-rabbit (Alexa Fluor 488, 1/1000) were tested, in parallel with a-P450’s to check specificity of each primary antibody. Finally images were obtained on Leica TCS-NT Laser Scanning microscope using the 40-objective.
The data sets supporting the results of this article are included within the article and its additional files. All microarray datasets are MIAME compliant  and deposited in ArrayExpress and VectorBase (Accession numbers E-MTAB-2808).
False discovery rate
World Health Organisation: World Malaria Report. 2013, Geneva, http://www.who.int/malaria/publications/world_malaria_report_2013/en/,
Edi CV, Djogbénou L, Jenkins AM, Regna K, Muskavitch MAT, Poupardin R, Jones CM, Essandoh J, Kétoh GK, Paine MJI: CYP6 P450 enzymes and ACE-1 duplication produce extreme and multiple insecticide resistance in the malaria mosquito anopheles gambiae. PLoS Genet. 2014, 10: e1004236-10.1371/journal.pgen.1004236.
Edi CV, Benjamin KG, Jones CM, Weetman D, Ranson H: Multiple-insecticide resistance in anopheles gambiae mosquitoes, Southern Côte d’Ivoire. Emerg Infect Dis. 2012, 18: 1508-1511. 10.3201/eid1809.120262.
Ranson H, N’Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V: Pyrethroid resistance in African anopheline mosquitoes: what are the implications for malaria control?. Trends Parasitol. 2011, 27: 91-98. 10.1016/j.pt.2010.08.004.
Liu N, Xu Q, Zhu F, Zhang L: Pyrethroid resistance in mosquitoes. Insect Sci. 2006, 13: 159-166. 10.1111/j.1744-7917.2006.00078.x.
Ranson H, Claudianos C, Ortelli F, Abgrall C, Hemingway J, Sharakhova MV, Unger MF, Collins FH, Feyereisen R: Evolution of supergene families associated with insecticide resistance. Science. 2002, 298 (80-): 179-181.
Chiu T-L, Wen Z, Rupasinghe SG, Schuler MA: Comparative molecular modeling of Anopheles gambiae CYP6Z1, a mosquito P450 capable of metabolizing DDT. Proc Natl Acad Sci. 2008, 105: 8855-8860. 10.1073/pnas.0709249105.
Mitchell SN, Stevenson BJ, Müller P, Wilding CS, Egyir-Yawson A, Field SG, Hemingway J, Paine MJI, Ranson H, Donnelly MJ: Identification and validation of a gene causing cross-resistance between insecticide classes in anopheles gambiae from Ghana. Proc Natl Acad Sci. 2012, 109 (16): 6147-6152. 10.1073/pnas.1203452109.
Müller P, Warr E, Stevenson BJ, Pignatelli PM, Morgan JC, Steven A, Yawson AE, Mitchell SN, Ranson H, Hemingway J, Paine MJI, Donnelly MJ: Field-caught permethrin-resistant anopheles gambiae overexpress CYP6P3, a P450 that metabolises pyrethroids. PLoS Genet. 2008, 4: e1000286-10.1371/journal.pgen.1000286.
Ismail HM, O’Neill PM, Hong DW, Finn RD, Henderson CJ, Wright AT, Cravatt BF, Hemingway J, Paine MJI: Pyrethroid activity-based probes for profiling cytochrome P450 activities associated with insecticide interactions. Proc Natl Acad Sci. 2013, 10 (49): 19766-19771.
Chandor-Proust A, Bibby J, Regent-Kloeckner M, Roux J, Guittard-Crilat E, Poupardin R, Asam Riaz M, Paine M, Dauphin-Villemant C, Reynaud S: The central role of mosquito cytochrome P450 CYP6Zs in insecticide detoxification revealed by functional expression and structural modelling. Biochem J. 2013, 455 (1): 75-85. 10.1042/BJ20130577.
Baker DA, Nolan T, Fischer B, Pinder A, Crisanti A, Russell S: A comprehensive gene expression atlas of sex-and tissue-specificity in the malaria vector. Anopheles gambiae. BMC Genomics. 2011, 12: 296-10.1186/1471-2164-12-296.
Koutsos AC, Blass C, Meister S, Schmidt S, MacCallum RM, Soares MB, Collins FH, Benes V, Zdobnov E, Kafatos FC, Christophides GK: Life cycle transcriptome of the malaria mosquito Anopheles gambiae and comparison with the fruitfly Drosophila melanogaster. Proc Natl Acad Sci. 2007, 104 (27): 11304-11309. 10.1073/pnas.0703988104.
Marinotti O, Calvo E, Nguyen QK, Dissanayake S, Ribeiro JMC, James AA: Genome-wide analysis of gene expression in adult Anopheles gambiae. Insect Mol Biol. 2006, 15: 1-12. 10.1111/j.1365-2583.2006.00610.x.
Neira Oviedo M, Vanekeris L, Corena-Mcleod MDP, Linser PJ: A microarray-based analysis of transcriptional compartmentalization in the alimentary canal of Anopheles gambiae (Diptera: Culicidae) larvae. Insect Mol Biol. 2008, 17: 61-72. 10.1111/j.1365-2583.2008.00779.x.
Johnson BR, Atallah J, Plachetzki DC: The importance of tissue specificity for RNA-seq: highlighting the errors of composite structure extractions. BMC Genomics. 2013, 14: 586-10.1186/1471-2164-14-586.
Harris C, Lambrechts L, Rousset F, Abate L, Nsango SE, Fontenille D, Morlais I, Cohuet A: Polymorphisms in Anopheles gambiae immune genes associated with natural resistance to Plasmodium falciparum. PLoS Pathog. 2010, 6: e1001112-10.1371/journal.ppat.1001112.
Chung H, Sztal T, Pasricha S, Sridhar M, Batterham P, Daborn PJ: Characterization of Drosophila melanogaster cytochrome P450 genes. Proc Natl Acad Sci. 2009, 106: 5731-5736. 10.1073/pnas.0812141106.
Morey JS, Ryan JC, Van Dolah FM: Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR. Biol Proced Online. 2006, 8: 175-193. 10.1251/bpo126.
Qiu Y, Tittiger C, Wicker-Thomas C, Le Goff G, Young S, Wajnberg E, Fricaux T, Taquet N, Blomquist GJ, Feyereisen R: An insect-specific P450 oxidative decarbonylase for cuticular hydrocarbon biosynthesis. Proc Natl Acad Sci. 2012, 109: 14858-14863. 10.1073/pnas.1208650109.
Jones CM, Toé HK, Sanou A, Namountougou M, Hughes A, Diabaté A, Dabiré R, Simard F, Ranson H: Additional Selection for Insecticide Resistance in Urban Malaria Vectors: DDT Resistance in Anopheles arabiensis from Bobo-Dioulasso, Burkina Faso. PLoS One. 2012, 7: e45995-10.1371/journal.pone.0045995.
Tene BF, Poupardin R, Costantini C, Awono-Ambene P, Wondji CS, Ranson H, Antonio-Nkondjio C: Resistance to DDT in an urban setting: common mechanisms implicated in both M and S forms of anopheles gambiae in the City of Yaoundé Cameroon. PLoS One. 2013, 8: e61408-10.1371/journal.pone.0061408.
Gibson GG, Skett P: Introduction to Drug Metabolism. 2001, Cheltenha, UK: Nelson Thornes
Coetzee M, Hunt RH, Wilkerson R, Della Torre A, Coulibaly MB, Besansky NJ: Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013, 3619: 246-274.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: 3-
Rossell D: GaGa: a parsimonious and flexible model for differential expression analysis. Ann Appl Stat. 2009, 3: 1035-1051. 10.1214/09-AOAS244.
Alexa A, Rahnenfuhrer J: topGO: Enrichment Analysis for Gene Ontology. R Package Version 2.18.0. 2010, http://www.bioconductor.org/packages/release/bioc/html/topGO.html,
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL: Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012, 13: 134-10.1186/1471-2105-13-134.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC: Minimum information about a microarray experiment (MIAME)—toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.
The research leading to these results was supported by the European Union Seventh Framework Programme FP7 (2007–2013) under grant agreement no 265660 AvecNet. We thank Dr Clare Strode, LSTM, for sharing unpublished microarray data and Dr Gareth Lycett, LSTM, for critical review of the manuscript.
The authors declare that they have no competing interests.
V.I, C.M.J. P.P and V.B. performed the experiments. H.R and J.V, designed the study. V.I, S.W and J.M. analysed the data. V.I and H.R drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Schematic of design of microarray experiment. a. Sample vs reference design, for each of the susceptible lab population N’Gousso and the resistant population Tiassalé b. Resistant (T) vs Susceptible (N) design, for each individual tissue. (PNG 332 KB)
Additional file 2: Table S1: Overview of probes over or under transcribed in each body part for both the resistant and susceptible strains when compared to the whole organism. The local transcription of probes in each of the two mosquito strains is expressed as total number of probes over or under transcribed in a particular body part, compared to the reconstituted whole organism, as a percentage of the total probes on the array. For genes that were represented by multiple probes, an average of all the probes was used. (XLSX 56 KB)
Additional file 3: Table S2: Genes significantly enriched or depleted in one or more body part compared to the whole organism. AGAP identifiers, description and the body part in which the gene is enriched with log2 fold change (relative to reconstituted whole transcriptome) indicated for each strain. Only probes outside the 95% intervals are listed. Blue indicates cytochrome p450s, green GSTs, purple COEs, orange UDPs and teal carboxylesterases. (XLS 253 KB)
Additional file 4: Table S3: Probe list from stringent analysis of direct transcriptome comparisons of dissected body parts from susceptible and resistant strains. AGAP identifier, description, body part that the probe fits the stringent selection criteria, whole organism array log2 fold change, resistant vs susceptible array log2 fold change, GaGa log2 fold change and log2 body part qPCR results for given tissue. qPCR validation has been performed on several candidates. Cells are coloured with a gradient dependent upon the directionality of the fold change, down regulated transcripts are indicated in red and up regulated in green. Sheets for both up regulated and down regulated genes are present. (PNG 151 KB)
Additional file 5: Figure S2: Local expression of all Cytochrome p450s, following a phylogenetic dendrogam. Full length protein sequence alignment and neighbour joining tree as computed on MEGA5 decorated with local log2 transcription of Tiassalé cytochrome p450s. (PNG 62 KB)
Additional file 6: Figure S3: qPCR validation of body part enrichment of cytochrome P450s. Difference in transcription levels in Tiassalé RNA from individual dissected body parts, compared to the reconstituted whole, were measured by qPCR for six P450 genes and compared to the data obtained from the microarray. Data represent log2 fold change with associated standard error bars. (XLSX 53 KB)
Additional file 7: Table S4: qPCR primer list. Forward and reverse primers used for all qPCR reactions. All primer products are between 80 and 150 base pairs and follow MIQE guidelines. (XLSX 46 KB)
About this article
Cite this article
Ingham, V.A., Jones, C.M., Pignatelli, P. et al. Dissecting the organ specificity of insecticide resistance candidate genes in Anopheles gambiae: known and novel candidate genes. BMC Genomics 15, 1018 (2014). https://doi.org/10.1186/1471-2164-15-1018
- Insecticide resistance