- Research article
- Open Access
Comparative analysis of response to selection with three insecticides in the dengue mosquito Aedes aegyptiusing mRNA sequencing
© David et al.; licensee BioMed Central Ltd. 2014
- Received: 7 November 2013
- Accepted: 21 February 2014
- Published: 5 March 2014
Mosquito control programmes using chemical insecticides are increasingly threatened by the development of resistance. Such resistance can be the consequence of changes in proteins targeted by insecticides (target site mediated resistance), increased insecticide biodegradation (metabolic resistance), altered transport, sequestration or other mechanisms. As opposed to target site resistance, other mechanisms are far from being fully understood. Indeed, insecticide selection often affects a large number of genes and various biological processes can hypothetically confer resistance. In this context, the aim of the present study was to use RNA sequencing (RNA-seq) for comparing transcription level and polymorphism variations associated with adaptation to chemical insecticides in the mosquito Aedes aegypti. Biological materials consisted of a parental susceptible strain together with three child strains selected across multiple generations with three insecticides from different classes: the pyrethroid permethrin, the neonicotinoid imidacloprid and the carbamate propoxur.
After ten generations, insecticide-selected strains showed elevated resistance levels to the insecticides used for selection. RNA-seq data allowed detecting over 13,000 transcripts, of which 413 were differentially transcribed in insecticide-selected strains as compared to the susceptible strain. Among them, a significant enrichment of transcripts encoding cuticle proteins, transporters and enzymes was observed. Polymorphism analysis revealed over 2500 SNPs showing > 50% allele frequency variations in insecticide-selected strains as compared to the susceptible strain, affecting over 1000 transcripts. Comparing gene transcription and polymorphism patterns revealed marked differences among strains. While imidacloprid selection was linked to the over transcription of many genes, permethrin selection was rather linked to polymorphism variations. Focusing on detoxification enzymes revealed that permethrin selection strongly affected the polymorphism of several transcripts encoding cytochrome P450 monooxygenases likely involved in insecticide biodegradation.
The present study confirmed the power of RNA-seq for identifying concomitantly quantitative and qualitative transcriptome changes associated with insecticide resistance in mosquitoes. Our results suggest that transcriptome modifications can be selected rapidly by insecticides and affect multiple biological functions. Previously neglected by molecular screenings, polymorphism variations of detoxification enzymes may play an important role in the adaptive response of mosquitoes to insecticides.
- RNA sequencing
- Insecticide resistance
- Detoxification enzymes
- Cytochrome P450 monooxygenase
Mosquitoes are vectors of several human diseases representing a major burden for public health worldwide . Half of the world’s population is exposed to malaria while dengue fever represents a burden in over 100 countries with 2.5 billion people at risk [2, 3]. Since the 1950s, chemical insecticides have been massively used for controlling mosquito populations but their efficacy is now threatened by resistance mechanisms developed by insects. In absence of efficient alternatives, characterizing molecular mechanisms underlying resistance is a key step for improving resistance management strategies.
Resistance to insecticides can be the consequence of different mechanisms, such as a mutation of the proteins targeted by the insecticide (target-site insensitivity), a lower penetration of the insecticide, its sequestration, or its biodegradation (metabolic resistance) [4, 5]. Target-site insensitivity and metabolic resistance are known as the two main resistance mechanisms in mosquitoes [4, 6]. Mutations causing target-site insensitivity are well-characterized in mosquitoes and molecular tests detecting these mutations are available for species of public health importance [7–12]. Metabolic resistance has been reported worldwide and usually involves detoxification enzymes such as cytochrome P450 monooxygenases (P450s or CYPs for genes), carboxy/cholinesterases (CCEs), glutathione S-transferases (GSTs) and UDP glucosyl-transferases (UGTs) [4, 6, 13, 14]. However, due to the large number of mosquito genes encoding detoxification enzymes [15–17] pinpointing those responsible for resistance remains challenging . Metabolic resistance has been mostly associated with an increased level of detoxification enzymes in resistant populations and multiple candidate genes have been identified by microarray screenings [6, 18–21]. In contrast, polymorphism variations potentially affecting the functionality of detoxification enzymes have been hardly studied in mosquitoes despite evidences suggesting that this phenomenon may play a role in insecticide resistance [14, 22]. Recently, polymorphism of a P450 gene has been associated with pyrethroid resistance in the mosquito culex pipiens and a reduction of sequence diversity in two P450 genes conferring resistance to pyrethroids has been observed in Anopheles funestus. This suggests that a deep analysis of the polymorphism associated with resistance can improve our understanding of mechanisms developed by mosquitoes to resist insecticides. Today, such knowledge gap can be overcome by high throughput sequencing approaches such as mRNA sequencing (RNA-seq), which can generates concomitantly gene expression and polymorphism data over the whole transcriptome from a single experiment [25–27].
In this context, the aim of the present study was to use RNA-seq for investigating transcription level and polymorphism variations associated with adaptation to three insecticides from distinct chemical families in the mosquito Aedes aegypti. A susceptible strain was selected with the pyrethroid permethrin, the neonicotinoid imidacloprid or the carbamate propoxur, to produce three resistant strains. After ten generations of selection, the constitutive resistance level of each resistant strain was measured and the transcriptome of each strain was deep sequenced. After mapping cDNA reads to the genome, gene expression and polymorphism variations linked to insecticide selection were identified and compared across strains. Results are discussed in regards of known and new putative adaptive mechanisms conferring insecticide resistance in mosquitoes.
Insecticide resistance levels
Resistance levels after insecticide selection
2.29 - 2.77
6.21 - 14.43
2.24 - 6.30
208.74 - 276.67
1236.24 - 1599.87
4.47 - 7.66
392.16 - 496.34
743.47 - 890.23
1.50 - 2.27
Sequencing, read mapping and genome re-annotation
More than 269 million 75 bp cDNA reads were sequenced across all samples (Additional file 1: Table S1). Each mosquito strain was represented by two cDNA libraries with an average of 33.6 million reads per library. More than 80% of reads were successfully mapped on the Ae. aegypti genome. Filtering on sequence quality and mapping score retained more than 174 million reads (~ 65%). For each strain, RPKMs (Reads Per Kilobase of exon model per Million reads) obtained from the two cDNA library replicates were well correlated (Additional file 2: Figure S1), indicating moderate variations between technical replicates. In consequence, reads from technical replicates were pooled for further analyses.
A transcription signal was detected for 85% of known Ae. aegypti genes. Comparing transcript coverage between strains revealed similar distributions with RPKMs spanning more than 6 Logs and a median transcript coverage of 4.4 RPKM (~300 reads/Kb) (Additional file 3: Figure S2). Comparing genome annotation with the distribution and coverage of reads suggested incorrect gene boundary annotations for more than 3000 transcripts and identified more than 500 isolated novel transcription events (NTEs) based on their transcription signal, structure and high distance to known transcripts (Additional file 4: Table S2 and Additional file 5: Table S3). Distribution of mapped reads across the whole Ae. aegypti genome can be accessed at http://vectorbase.org using the ‘configure this page’ and ‘RNAseq alignments’ options of the genome browser (tracks ‘Bora-Bora control’, ‘Perm-R’, ‘Imida-R’ and ‘Propo-R’).
Differential transcription in insecticide-selected strains
Differential transcription analysis overview
AC test P value
AC test P value and FC >3
New putative transcripts
New putative transcripts
Polymorphism variations in insecticide-selected strains
Differential polymorphism analysis overview
Gene ontology terms enrichment analyses
Gene Ontology (GO) term enrichment analyses were performed independently on transcripts significantly over- or under transcribed or affected by differential SNPs in their coding region (Additional file 9: Figure S4). These analyses confirmed the distinct response of each strain to insecticide selection. Among GO terms representing genes over transcribed in insecticide-selected strains, terms related to cuticle proteins were strongly enriched in the Imida-R strain. Terms related to oxygen transport or to hexamerins were enriched in both the Perm-R and Imida-R strains, together with those related to iron homeostasis and iron binding. Among GO terms represented by under transcribed genes, multiple terms related to immune response were enriched in insecticide-selected strains. When considering transcripts affected by differential SNPs, no GO term was enriched in the Imida-R strain while only ‘dehydrogenase activity’ was enriched in the Propo-R strain. In contrast, the Perm-R strain revealed an enrichment of all terms related to cytochrome P450s such as ‘monooxygenases activity', ‘electron carrier activity', ‘tetrapyrrole binding', ‘heme binding’ and ‘iron ion binding’. Terms related to hexamerins were also enriched in Perm-R together with terms related to various enzyme families.
Focus on transcripts potentially involved in insecticide detoxification
Characterizing resistance mechanisms is essential for improving resistance management. Although target site modifications play a major role in resistance [4, 5, 28, 29], other mechanisms such as insecticide biodegradation, altered transport, sequestration and modification of the insect cuticle also account for a significant part of resistance [6, 18, 30, 31]. However, the intricacy of these mechanisms makes challenging the identification of candidate genes for functional validation . Gene expression microarrays are mostly used for identifying genes differentially transcribed in resistant populations but suffer from technical biases [21, 32, 33]. In contrast RNA-seq generates transcription data with a higher resolution, better dynamic range and lower technical variation . In addition, RNA-seq produces polymorphism data and useful information to re-annotate gene models [25, 26]. Despite recent studies pointing out the role of polymorphism variations in insecticide resistance [22–24], such aspect has never been investigated at the transcriptome level in mosquitoes. Indeed, the only study using RNA-seq to investigate insecticide resistance in mosquitoes did not consider polymorphism variations . In this regard, the present study represents the first attempt to use RNA-seq for examining concomitantly quantitative and qualitative transcriptome changes associated with resistance to different insecticides in mosquitoes.
Insecticide selection and resistance levels
After 10 generations of selection, bioassays revealed a constitutive increased resistance of each selected strain to its respective insecticide compared to the susceptible parental strain. Although resistance levels were low as compared to what can be observed in natura, they were significant regarding the few generations of selection, the full susceptibility of the parental strain and the absence of target-site resistance alleles. The rapid rise of resistance suggests that alleles conferring a better fitness in presence of insecticides are already present in susceptible populations and can be promptly selected under a constant selection pressure. Although laboratory selection does not fully mimic adaptive processes occurring in natura (e.g. lower population sizes, different environment, no introduction of resistant alleles by migration …), selecting all resistant strains from a single susceptible strain allowed minimizing variations related to different genetic backgrounds. Using a fully susceptible parental strain also allowed focusing on non-target site resistance mechanisms (no target-site mutations detected in parental or selected strains).
RNA-seq as a tool for studying insecticide resistance in mosquitoes
A total of 270 million cDNA reads were sequenced across all samples and 80% of them were successfully mapped to Ae. aegypti reference genome. Such mapping efficiency appeared acceptable considering polymorphism variations occurring between the reference genome (Liverpool strain) and the parental strain used in our study (Bora-Bora strain). The high correlation of expression data obtained from cDNA library replicates confirmed the robustness of cDNA library construction and sequencing procedures. By applying high-stringency sequence quality and transcript coverage filtering, high fidelity transcription data were recovered for more than 13000 transcripts. Such detection level was comparable to those obtained with DNA microarrays at the same life stage [35, 36]. More than 2500 differential SNPs linked to insecticide selection were identified. Although our experimental design did not control for stochastic effects (genetic drift) and the presence of false positives is likely, transcripts affected by these differential SNPs represent strong candidates for further functional validation studies. Little overlap was found between transcripts differentially transcribed in selected strains and those affected by differential SNPs. This was expected as RNA-seq data are restricted to transcripts and did not cover regulatory regions often located outside transcript boundaries. Challenging reads distribution with genome annotation identified > 3000 transcripts incorrectly annotated with most of them showing wrong UTR boundaries or modification of exon/intron structure. In addition, more than 500 lonely genomic regions showing high transcription signals and realistic exon/intron structures were identified. Further analyses are now required for assigning them to new exons of known transcripts, novel transcripts, pseudogenes or non-coding RNAs.
Transcriptome changes associated with insecticide selection
Overall, our study revealed diverse response patterns depending on the insecticide used for selection. The Imida-R strain selected with the neonicotinoid imidacloprid showed the highest resistance level together with the strongest differential transcription response with numerous genes being over transcribed. In contrast, the Perm-R strain selected with the pyrethroid permethrin revealed a moderate differential transcription response but a lot of polymorphism variations following insecticide selection. The Propo-R strain selected with the carbamate propoxur showed an intermediate pattern. These dissimilar quantitative and qualitative transcriptome changes may reflect different adaptive strategies driven by costs and benefits associated with resistance mechanisms to each insecticide.
Insecticide resistance and immunity
Multiple transcripts related to immune response were affected by insecticide selection with several of them under transcribed in all resistant strains (cecropins, defensins, lectins, …) and others affected by differential SNPs (defensins, clip-domain proteases, spatzle proteins…). Mosquito humoral response is involved in their capacity to host and transmit viruses  or parasites . As mentioned in Alout et al. , our results support an association between insecticide resistance and the capacity of mosquitoes to host and transmit pathogens, which may affect the control of mosquito-borne diseases .
Altered insecticide transport and sequestration
Multiple transcripts encoding cuticular proteins were over transcribed in the Imida-R strain while most of them were under transcribed in the Perm-R strain and not affected in the Propo-R strain. Such strong effect in the Imida-R strain was not associated to changes in the polymorphism of these transcripts, suggesting that insecticide selection is affecting their regulation rather than their conformation. Cuticle plays a crucial role in protecting insects from their environment. The vast majority of chemical insecticides are active by contact and changes in cuticle thickness or conformation have been suggested to contribute to resistance in mosquitoes [31, 41, 42]. Our data support previous studies suggesting a significant role of cuticle in the adaptation to neonicotinoid insecticides [36, 43, 44]. The role of cuticle in response to imidacloprid selection was supported by the specific over transcription of the multi-copper oxidase AAEL007415 in the Imida-R strain as its An. gambiae orthologue is involved in cuticle and egg shell tanning [45, 46]. Numerous protein families are involved in cuticle biosynthesis and homeostasis including enzymes, transporters and transcription factors and further studies are now required for pinpointing those controlling cuticular resistance in mosquitoes.
Several kinases were specifically over-transcribed in the Imida-R strain. Kinases are involved in multiple regulatory mechanisms and their involvement in the response of insects to insecticides is likely. Indeed, recent studies showed that the phosphorylation state of acetylcholine nicotinic receptors can modulate the efficacy of neonicotinoid insecticides [47–49].
Five transcripts encoding hexamerins were strongly over transcribed in the Perm-R strain and in a lesser extend in other resistant strains. These transcripts, located on different supercontigs, were also affected by differential SNPs suggesting a selection imprint on these proteins. Insect hexamerins may be involved in hormone transport, energy and amino acid storage, cuticle biosynthesis and immune defense [50, 51]. Hexamerins of the lepidopteran Heliothis zea have been shown to bind insecticides, suggesting a direct role in resistance by sequestration or altered transport . Deciphering if hexamerins are impacted by insecticide selection because of their ability to bind insecticides, their interaction with cuticle homeostasis or their role in resources re-allocation associated to fitness costs remains unclear.
Finally, ATP-binding cassette transporters (ABC transporters) can play a role in adaptation to xenobiotics  and have been associated to insecticide resistance in mosquitoes [30, 54]. In our study, the response of ABC transporters to insecticide selection was marginal. Indeed, only one ABC transporter (AAEL008624) was found differentially transcribed in response to insecticide selection and this gene was down regulated in all resistant strains. Four others were affected by differential SNPs but their allele frequency variations in resistance strains were low.
Detoxification enzymes play a major role in insecticide resistance [4, 6, 13, 14, 18, 55]. As expected, these enzymes were well represented in our data set but showed distinct patterns depending on the nature of the insecticide used for selection.
Response to selection with the neonicotinoid imidacloprid (Imida-R strain) was characterized by the over transcription of multiple P450s, oxidases, transferases and one alcohol dehydrogenase, supporting the involvement of multiple enzymes in imidacloprid biodegradation pathways. Among them, the P450 CYP6BB2 was recently pointed out as a solid candidate for imidacloprid metabolism based on gene expression data and substrate binding predictions . The aldehyde oxidase AAEL002683 was over transcribed and affected by polymorphism variations in the Imida-R strain. A recent study confirmed that aldehyde oxidases can contribute to neonicotinoid metabolism through nitro-reduction .
Several P450s over transcribed in the Imida-R strain were also found over transcribed in response to propoxur selection. Cross resistance between these two insecticides was identified  and the potency of particular P450s to confer resistance to multiple insecticides has been shown. For example, An. gambiae CYP6Z1 metabolizes the organochlorine DDT and the carbamate carbaryl  while CYP6M2 metabolizes both DDT and pyrethroids [57, 58].
Metabolic resistance of mosquitoes to pyrethroids has been mainly associated with an over expression of P450s able to metabolize them. Among the multiple candidates identified by microarray screenings, An. gambiae CYP6M2 and CYP6P3, An. funestus CYP6P9, Ae. aegypti CYP9J32, CYP9J24, CYP9J28, and Cx quinquefasciatus CYP6M10 and CYP4H24 have been validated as pyrethroid metabolizers [review in 18]. CYP6Zs and CYP6Ms have also been associated with pyrethroid resistance by QTL in An. funestus where resistance mainly relies on metabolic mechanisms . Recently, the central role of mosquito CYP6Zs in pyrethroid degradation pathway was revealed . In the present study few P450s were found over transcribed in the Perm-R strain. Among them, CYP6BB2 was also found strongly over transcribed in pyrethroid resistant populations from Cuba and Cayman islands . Unexpectedly, known Ae. aegypti pyrethroid metabolizers or paralogs of those validated from other mosquito species were not found over transcribed in the Perm-R strain. However, a strong selection imprint was detected in the Perm-R strain in a P450 cluster in supercontig 1.371containing strong candidates CYP6Ms, CYP6Ns and CYP6Zs, suggesting that the selection of particular P450 variants can contribute to pyrethroid resistance. In mammals, it is well known that P450 variants can display different substrate specificity [60–64]. To date, such variations have been neglected in mosquitoes with only few studies pointing out P450 alleles associated with insecticide resistance [22, 23]. Finally, the strong selection imprint observed in these P450s might explain their unexpected under transcription in the Perm-R strain. Indeed, the apparent under expression of transcripts affected by differential SNPs might be the consequence of a mapping bias due to a higher divergence from the reference genome or an enrichment in low-expressed alleles [65, 66]. Although further analyses are required for investigating the role of allele-specific expression in resistance, our data supports the selection of particular detoxification enzyme variants by insecticides.
The present study primarily aimed at assessing the usefulness of RNA-seq to investigate insecticide resistance mechanisms in mosquitoes. Results confirmed that this technique produces high-quality gene expression data together with solid polymorphism data. Distinct responses to selection with insecticides from different chemical families were observed with a balance between gene expression and polymorphism variations. Polymorphism variations of P450 enzymes were strongly linked to pyrethroid selection. Although additional analyses are required to validate variants linked to resistance, such finding highlights the necessity to consider both gene expression and polymorphism variations for identifying candidate genes potentially involved in insecticide resistance. As sequencing costs are decreasing and new sequencing strategies and bioinformatics pipelines are developed, obtaining gene expression and polymorphism data from the same samples using high-throughput sequencing should now be considered as a valuable alternative to microarrays.
Mosquito selection with insecticides and bioassays
The mosquito Ae. aegypti, was used in the present study. Mosquitoes were reared in standard insectary conditions (26°C, 14 h/10 h light/dark, 80% relative humidity) in tap water (larvae) and net cages (adults). Larvae and adults were fed with hay pellets and papers impregnated with honey respectively. Blood feeding of adult females was performed on mice. Mice were maintained in an animal house agreed by French Ministry of animal welfare (n° B 38 421 10 001) and used in accordance to EU laws and the relevant ethic committee recommendations (ComEth Grenoble - C2EA - 12). The laboratory strain Bora-Bora, originating from French Polynesia and fully susceptible to insecticides, was used as a parental strain to select three independent strains with the pyrethroid insecticide permethrin (Perm-R strain), the neonicotinoid insecticide imidacloprid (Imida-R strain) and the carbamate insecticide propoxur (Propo-R strain). Both pyrethroid and carbamate insecticides are heavily used against mosquitoes and resistance to these insecticides has been reported worldwide. Neonicotinoids are marginaly used against mosquitoes but represent one possible alternative when resistance to other insecticide threatens the efficacy of mosquito control . Selection was performed by exposing early 4th-stage larvae for 24 h to a lethal dose of each insecticide. For each strain, the dose of insecticide was adjusted at each generation (4 to 6 μg/L permethrin, 500 to 1000 μg/L imidacloprid and 500 to 800 μg/L propoxur) in order to reach 60-80% larval mortality after 24 h exposure. Surviving larvae were transferred in clean tap water, fed with standard larval food and allowed to emerge. Adults were allowed to reproduce for 4 days and blood-fed to obtain eggs for the next generation. In order to limit bottleneck effects, each generation was seeded with more than 6000 individuals. Selection process was carried out in parallel for all strains during ten generations. During this process, the susceptible parental strain was maintained in similar condition without insecticide selection. Bioassays and molecular work were performed on early 4th-stage larvae of the 11th generation (G11 larvae) bred in standard conditions and not exposed to insecticides.
To assess the constitutive resistance level of each selected strain, larval bioassays were conducted with permethrin, imidacloprid and propoxur comparatively to the susceptible parental strain. Five doses of each insecticide and four replicates of 25 larvae per dose were used. Doses of permethrin (1.5 to 6.5 μg/L), imidacloprid (150 to 2200 μg/L) and propoxur (100 to 2000 μg/L) were chosen in order to cover the whole mortality range after 24 h exposure. Lethal concentrations corresponding to 50% mortality (LC50) and their 95% confident intervals (CI95%) were then calculated with a probit approach for each strain using XL-Stat (Addinsoft, Paris, France). Resistance ratios (RR50 based on LC50 values) were calculated by comparison to the susceptible parental strain.
RNA extraction and cDNA libraries preparation
For each strain, total RNA was extracted from G11 larvae obtained from three independent egg batches (three biological replicates per strain). Each biological replicate consisted of 180 larvae reared in 200 mL tap water in standardized insectary conditions. For each biological replicate, total RNA was extracted from 60 4th-stage larvae using the RNAqueous-4PCR kit (Applied Biosystems/Ambion, Austin, TX, USA), according to the manufacturer’s instructions. Total RNA quality and quantity were assessed with a Nanodrop ND1000 (Thermo Scientific, USA) and a 2100 Bioanalyzer (Agilent, USA). After extraction, total RNAs from each biological replicate were pooled in equal quantity to obtain a total RNA mixture representative of 180 individuals. Total RNA pools from each strain were then used for preparing cDNA libraries using the mRNA-Seq Sample Prep Kit (Part 1004898 Rev D, Illumina, USA). Two replicates of cDNA libraries were prepared for each strain as follows. Briefly, mRNAs were purified using poly-T beads and chemically fragmented. These fragmented mRNAs were reverse-transcribed using Superscript II (Invitrogen) at 42°C for 50 min. Double-stranded cDNAs were then synthesized and mRNAs removed using DNA polymerase I and RNase H at 16°C for 2.5 hours. Double-stranded cDNAs were purified using QIAquick PCR Purification Kit (QIAGEN, Germany) and processed for end-repair and 3′ adenylation using Klenow polymerase. Sequencing adaptors were then ligated using DNA ligase. Adapter-ligated cDNA libraries were then purified on 2% agarose gel based on a size range of 200 ± 25 bp. Adapter-ligated cDNA libraries were then enriched by 15 PCR cycles with adaptor-specific primers using Phusion DNA polymerase (Finnzymes Oy). Enriched cDNA libraries were then purified using QIAquick PCR Purification Kit (QIAGEN, Germany) before quality control analysis on an Agilent 2100 Bioanalyzer.
Sequencing, read mapping and genome re-annotation
Each double-stranded cDNA library was sequenced as single reads of 75 bp in a distinct flow cell lane with a Genome Analyzer II (Illumina) at the National Sequencing Center (Genoscope, Evry, France). Reads were then mapped to the Ae.aegypti genome sequence (Aaeg L1.2 gene set, http://vectorbase.org) using the Tophat algorithm with default parameters (http://tophat.cbcb.umd.edu, release 1.0.14) , leading to an average mapping rate of 80%. Bam files were then loaded into Genespring NGS version 12.5 (Agilent Technologies) for further analysis. Reads were then filtered based on the sequence quality (mean read quality ≥ 30 and with < 10 Ns) and mapping quality (alignment score ≥ 98), and non-primary multiply mapped reads were removed. The remaining reads were used for all subsequent analyses. Read coverage was computed and genome annotation was challenged based on read distribution and coverage. Default feature detection parameters were used (min exon length percentile 10; min intro length percentile 10; max intron length percentile 90; min exon RPKM percentile 50; min gene RPKM percentile 50; min gene length percentile 10; min exon RPKM with respect to host gene percentage 75; min number of reads in exon 10). Modified exon-intron structures and new putative transcripts were proposed based on read coverage, read splicing events and distance with respect to existing genes and transcripts. Reads per Kilobase exon Model per Million sequenced reads (RPKMs) were calculated for all known and putative new transcripts. Transcript RPKM values were used for assessing the variability between library replicates for each strain.
Differential transcription analysis
In order to avoid estimating transcription ratios from low coverage transcripts, only the 13105 transcripts showing at least 0.5 RPKM per condition (~ 30 reads/kb) were considered for differential transcription analysis. Differential transcription of each transcript was then tested between insecticide-selected strains and the susceptible parental strain using an Audic-Claverie test (AC test) based on read counts  with a Benjamini and Hochberg’s multiple testing correction . Transcripts showing an adjusted P value < 10-15 and a fold change > 3 in either direction were considered as differentially transcribed between insecticide-selected strains and the parental susceptible strain. RNA-seq transcription data from the Imida-R strain were compared with transcription data obtained from the same biological samples with the Aedes aegypti 15 K DNA microarray as previously described . Only the 326 transcripts showing significant adjusted P values in both techniques were considered.
Polymorphism detection and differential analysis
Detection of polymorphisms was performed based on the 174 million reads passing quality filters with the following parameters: confidence score threshold = 100, coverage > 20 reads, base quality cut off = 5, ignore locations within or next to homopolymer stretches > 10 nucleotides. Among all detected polymorphism variations, only SNP substitutions were considered for differential polymorphism analyses. SNP allele frequencies were then computed between each insecticide-selected strain and the susceptible strain. Allele frequencies were considered as differential between an insecticide-selected strain and the susceptible strain (hereafter named as differential SNPs) if the following conditions were fulfilled: i) Total read coverage at SNP position between both strains ≥ 50, ii) Strand bias at SNP position ≤ 50 for both strains and iii) Allele frequency difference between both strains > 50% in either direction. Potential genic effects of SNPs were computed by comparing SNP with reference genome annotation. Genic effects were defined as 5′UTR, Synonymous, Non-synonymous, Intronic, 3′UTR and Intergenic (i.e. close but not within gene boundaries).
GO term enrichment analyses
Go term enrichment analyses were conducted on i) transcripts significantly over transcribed in insecticide-selected strains compared to the susceptible strain, ii) transcripts significantly under transcribed in insecticide-selected strains compared to the susceptible strain and iii) transcripts affected by differential SNPs in their coding sequence as compared to the susceptible strain. GO term enrichment analyses were performed separately on each transcript list versus all detected transcripts. Enrichment statistics were computed based on hypergeometric distribution and adjusted with Benjamini-Yekutieli’s multiple testing correction  which takes into account the dependency among GO terms. GO terms showing an adjusted P value ≤ 0.05 were considered as significantly enriched in insecticide-selected strains compared to the susceptible strain.
All clustering analyses were performed with TM4 MEV version 4.3.02  using Euclidean distance and complete linkage algorithm with optimization of gene and condition trees. Transcripts significantly differentially transcribed (see above for thresholds) in any insecticide-selected strain compared to the parental susceptible strain were clustered based on their TR (log2 ratios). Annotated transcripts represented in the main clusters were then assigned to nine biological categories based on their annotation (including GO terms). Categories were defined as follow: ‘Detoxification', ‘Kinases/Phosphatases', ‘Other enzymes', ‘Cuticle', ‘Immunity', ‘Hormones and neurotransmitters signaling', ‘Transcription factors', ‘Intra or extracellular trafficking/chaperonins’ and ‘Structure’. Enrichment of these categories compared to all detected transcripts was then computed using a one-side Fisher’s exact test followed by Benjamini and Hochberg multiple testing correction . Categories showing a corrected P value < 0.05 were considered significantly enriched. A similar clustering analysis was performed with differential SNPs (see above for detailed criteria) falling in coding regions. Clustering was performed as described above using allele frequency variations (-100% to +100%) in each strain as compared to the susceptible strain. Differential SNPs were assigned to nine different biological categories (see above) in regard of the nature of the transcript affected. For each cluster, biological category enrichment was computed versus all detected transcripts as described above.
Clustering analyses were then focused on transcripts encoding genes potentially involved in insecticide detoxification pathways (detoxification sensus lato). Clustering of both transcription level variations and differential SNPs were performed as described above. Separate analyses were conducted for cytochrome P450 monoooxygenases (CYPs) and other detoxification transcripts.
RNA-seq sequence data have been deposited at EBI/SRA/ArrayExpress under the accession number E-MTAB-1635. Final gene expression data (RPKMs, Fold changes and adjusted P values) are available at http://vectorbase.org. Distribution of mapped reads across mosquito genome can visualized using the ‘configure this page’ and ‘RNAseq alignments’ options of the vectorbase Ae. aegypti genome browser (tracks ‘Bora-Bora control’, ‘Perm-R’, ‘Imida-R’ and ‘Propo-R’).
This work was funded by a research grant from the French Agence Nationale de la Recherche (ANR Mosquito-Env 07SEST014 project). We also acknowledge funding from the federative structure Environnemental and Systems Biology (BEeSy) of Grenoble – Alpes University. We thank Dr. F. Boyer, and Dr. L. Zinger for their helps with bioinformatics and statistics. We are grateful to T. Gaude and S. Veyrenc for their technical assistance in mosquito rearing and sample preparation. We thank Prof. C. Gauthier (PRABI) and Dr. C. Oger (DTAMB/PRABI) for technical help and/or discussion in the project setup. We also acknowledge the UMR 5558 “Pôle informatique” and the PRABI-AMSB to provide access to their mutualised computational cluster. We thank Prof. V. Marchais and Prof. U. Schlattner for useful comments on the manuscript. We thank Vectorbase (Dr. D. Lawson and Dr. R. MacCallum) for their help with data deposition. We appreciated efforts of the Genespring NGS technical support team for their help with data analysis.
- Lounibos LP: Invasions by insect vectors of human disease. Annu Rev Entomol. 2002, 47: 233-266. 10.1146/annurev.ento.47.091201.145206.PubMedView ArticleGoogle Scholar
- WHO: Dengue and dengue hemorrhagic fever. Fact sheet N 117. 2009, http://www.who.int/mediacentre/factsheets/fs117/en/,Google Scholar
- WHO: World Malaria Report. World Health Organization technical report. 2012, Geneva, Switzerland: WHO press, 195-http://www.who.int/malaria/publications/world_malaria_report_2012/report/en/ ISBN 978 92 4 156453 3,Google Scholar
- Hemingway J, Hawkes NJ, McCarroll L, Ranson H: The molecular basis of insecticide resistance in mosquitoes. Insect Biochem Mol Biol. 2004, 34 (7): 653-665. 10.1016/j.ibmb.2004.03.018.PubMedView ArticleGoogle Scholar
- Hemingway J, Ranson H: Insecticide resistance in insect vectors of human disease. Annu Rev Entomol. 2000, 45: 371-391. 10.1146/annurev.ento.45.1.371.PubMedView ArticleGoogle Scholar
- Nkya TE, Akhouayri I, Kisinza W, David JP: Impact of environment on mosquito response to pyrethroid insecticides: facts, evidences and prospects. Insect Biochem Mol Biol. 2013, 43: 407-416. 10.1016/j.ibmb.2012.10.006.PubMedView ArticleGoogle Scholar
- Du W, Awolola TS, Howell P, Koekemoer LL, Brooke BD, Benedict MQ, Coetzee M, Zheng L: Independent mutations in the Rdl locus confer dieldrin resistance to Anopheles gambiae and An. arabiensis. Insect Mol Biol. 2005, 14 (2): 179-183. 10.1111/j.1365-2583.2005.00544.x.PubMedView ArticleGoogle Scholar
- Kasai S, Ng LC, Lam-Phua SG, Tang CS, Itokawa K, Komagata O, Kobayashi M, Tomita T: First detection of a putative knockdown resistance gene in major mosquito vector, Aedes albopictus. Jpn J Infect Dis. 2011, 64 (3): 217-221.PubMedGoogle Scholar
- Martinez-Torres D, Chandre F, Williamson MS, Darriet F, Berge JB, Devonshire AL, Guillet P, Pasteur N, Pauron D: Molecular characterization of pyrethroid knockdown resistance (kdr) in the major malaria vector Anopheles gambiae s.s. Insect Mol Biol. 1998, 7 (2): 179-184. 10.1046/j.1365-2583.1998.72062.x.PubMedView ArticleGoogle Scholar
- Martinez-Torres D, Chevillon C, Brun-Barale A, Bergé JB, Pasteur N, Pauron D: Voltage-dependent Na + channels in pyrethroid-resistant Culex pipiens mosquitoes. Pestic Sci. 1999, 55: 1012-1020. 10.1002/(SICI)1096-9063(199910)55:10<1012::AID-PS39>3.0.CO;2-5.View ArticleGoogle Scholar
- Ranson H, Jensen B, Vulule JM, Wang X, Hemingway J, Collins FH: Identification of a point mutation in the voltage-gated sodium channel gene of Kenyan Anopheles gambiae associated with resistance to DDT and pyrethroids. Insect Mol Biol. 2000, 9 (5): 491-497. 10.1046/j.1365-2583.2000.00209.x.PubMedView ArticleGoogle Scholar
- Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, Raymond M: The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004, 13 (1): 1-7. 10.1111/j.1365-2583.2004.00452.x.PubMedView ArticleGoogle Scholar
- Feyereisen R: Insect cytochrome P450. Comprehensive Molecular Insect Science. Edited by: Gilbert LI, Iatrou K, Gill S. 2005, Amsterdam, Netherlands: Elsevier, 1-77.View ArticleGoogle Scholar
- Li X, Schuler MA, Berenbaum MR: Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics. Annu Rev Entomol. 2007, 52: 231-253. 10.1146/annurev.ento.51.110104.151104.PubMedView ArticleGoogle Scholar
- Arensburger P, Megy K, Waterhouse RM, Abrudan J, Amedeo P, Antelo B, Bartholomay L, Bidwell S, Caler E, Camara F, Campbell CL, Campbell KS, Casola C, Castro MT, Chandramouliswaran I, Chapman SB, Christley S, Costas J, Eisenstadt E, Feschotte C, Fraser-Liggett C, Guigo R, Haas B, Hammond M, Hansson BS, Hemingway J, Hill SR, Howarth C, Ignell R, Kennedy RC, et al: Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics. Science. 2010, 330 (6000): 86-88. 10.1126/science.1191864.PubMed CentralPubMedView ArticleGoogle Scholar
- Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, Salzberg SL, Loftus B, Yandell M, Majoros WH, Rusch DB, Lai Z, Kraft CL, Abril JF, Anthouard V, Arensburger P, Atkinson PW, Baden H, de Berardinis V, Baldwin D, Benes V, Biedler J, Blass C, Bolanos R, Boscus D, Barnstead M, et al: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298 (5591): 129-149. 10.1126/science.1076181.PubMedView ArticleGoogle Scholar
- Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, Loftus B, Xi Z, Megy K, Grabherr M, Ren Q, Zdobnov EM, Lobo NF, Campbell KS, Brown SE, Bonaldo MF, Zhu J, Sinkins SP, Hogenkamp DG, Amedeo P, Arensburger P, Atkinson PW, Bidwell S, Biedler J, Birney E, Bruggner RV, Costas J, Coy MR, Crabtree J, Crawford M, et al: Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007, 316 (5832): 1718-1723. 10.1126/science.1138878.PubMedView ArticleGoogle Scholar
- David JP, Ismail HM, Chandor-Proust A, Paine MJ: Role of cytochrome P450s in insecticide resistance: impact on the control of mosquito-borne diseases and use of insecticides on Earth. Philos Trans R Soc Lond B Biol Sci. 2013, 368: 20120429-10.1098/rstb.2012.0429.PubMed CentralPubMedView ArticleGoogle Scholar
- David JP, Strode C, Vontas J, Nikou D, Vaughan A, Pignatelli PM, Louis C, Hemingway J, Ranson H: The Anopheles gambiae detoxification chip: a highly specific microarray to study metabolic-based insecticide resistance in malaria vectors. Proc Natl Acad Sci U S A. 2005, 102 (11): 4080-4084. 10.1073/pnas.0409348102.PubMed CentralPubMedView ArticleGoogle Scholar
- Strode C, Wondji CS, David JP, Hawkes NJ, Lumjuan N, Nelson DR, Drane DR, Karunaratne SH, Hemingway J, Black WC, Ranson H: Genomic analysis of detoxification genes in the mosquito Aedes aegypti. Insect Biochem Mol Biol. 2008, 38 (1): 113-123. 10.1016/j.ibmb.2007.09.007.PubMedView ArticleGoogle Scholar
- Vontas J, Ranson H, Alphey L: Transcriptomics and disease vector control. BMC Biol. 2010, 8: 52-10.1186/1741-7007-8-52.PubMed CentralPubMedView ArticleGoogle Scholar
- Chiu TL, Wen Z, Rupasinghe SG, Schuler MA: Comparative molecular modeling of Anopheles gambiae CYP6Z1, a mosquito P450 capable of metabolizing DDT. Proc Natl Acad Sci USA. 2008, 105 (26): 8855-8860. 10.1073/pnas.0709249105.PubMed CentralPubMedView ArticleGoogle Scholar
- Hardstone MC, Komagata O, Kasai S, Tomita T, Scott JG: Use of isogenic strains indicates CYP9M10 is linked to permethrin resistance in Culex pipiens quinquefasciatus. Insect Mol Biol. 2010, 19 (6): 717-726. 10.1111/j.1365-2583.2010.01030.x.PubMedView ArticleGoogle Scholar
- Riveron JM, Irving H, Ndula M, Barnes KG, Ibrahim SS, Paine MJ, Wondji CS: Directionally selected cytochrome P450 alleles are driving the spread of pyrethroid resistance in the major malaria vector Anopheles funestus. Proc Natl Acad Sci USA. 2013, 110 (1): 252-257. 10.1073/pnas.1216705110.PubMed CentralPubMedView ArticleGoogle Scholar
- Ozsolak F, Milos PM: RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011, 12 (2): 87-98. 10.1038/nrg2934.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63. 10.1038/nrg2484.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilhelm BT, Landry JR: RNA-Seq-quantitative measurement of expression through massively parallel RNA-sequencing. Methods. 2009, 48 (3): 249-257. 10.1016/j.ymeth.2009.03.016.PubMedView ArticleGoogle Scholar
- Donnelly MJ, Corbel V, Weetman D, Wilding CS, Williamson MS, Black WC: Does kdr genotype predict insecticide-resistance phenotype in mosquitoes?. Trends Parasitol. 2009, 25 (5): 213-219. 10.1016/j.pt.2009.02.007.PubMedView ArticleGoogle Scholar
- Labbe P, Lenormand T, Raymond M: On the worldwide spread of an insecticide resistance gene: a role for local selection. J Evol Biol. 2005, 18 (6): 1471-1484. 10.1111/j.1420-9101.2005.00938.x.PubMedView ArticleGoogle Scholar
- Bariami V, Jones CM, Poupardin R, Vontas J, Ranson H: Gene amplification, ABC transporters and cytochrome P450s: unraveling the molecular basis of pyrethroid resistance in the dengue vector, Aedes aegypti. PLoS Negl Trop Dis. 2012, 6 (6): e1692-10.1371/journal.pntd.0001692.PubMed CentralPubMedView ArticleGoogle Scholar
- Wood O, Hanrahan S, Coetzee M, Koekemoer L, Brooke B: Cuticle thickening associated with pyrethroid resistance in the major malaria vector Anopheles funestus. Parasites Vectors. 2010, 3: 67-10.1186/1756-3305-3-67.PubMed CentralPubMedView ArticleGoogle Scholar
- Martin-Magniette ML, Aubert J, Cabannes E, Daudin JJ: Evaluation of the gene-specific dye bias in cDNA microarray experiments. Bioinformatics. 2005, 21 (9): 1995-2000. 10.1093/bioinformatics/bti302.PubMedView ArticleGoogle Scholar
- Spiess AN, Mueller N, Ivell R: Amplified RNA degradation in T7-amplification methods results in biased microarray hybridizations. BMC Genomics. 2003, 4 (1): 44-10.1186/1471-2164-4-44.PubMed CentralPubMedView ArticleGoogle Scholar
- Bonizzoni M, Afrane Y, Dunn WA, Atieli FK, Zhou G, Zhong D, Li J, Githeko A, Yan G: Comparative transcriptome analyses of deltamethrin-resistant and -susceptible Anopheles gambiae mosquitoes from Kenya by RNA-Seq. PLoS One. 2012, 7 (9): e44607-10.1371/journal.pone.0044607.PubMed CentralPubMedView ArticleGoogle Scholar
- Poupardin R, Riaz MA, Jones CM, Chandor-Proust A, Reynaud S, David JP: Do pollutants affect insecticide-driven gene selection in mosquitoes? Experimental evidence from transcriptomics. Aquat Toxicol. 2012, 114–115: 49-57.PubMedView ArticleGoogle Scholar
- Riaz MA, Chandor-Proust A, Dauphin-Villemant C, Poupardin R, Jones CM, Strode C, Regent-Kloeckner M, David JP, Reynaud S: Molecular mechanisms associated with increased tolerance to the neonicotinoid insecticide imidacloprid in the dengue vector Aedes aegypti. Aquat Toxicol. 2013, 126: 326-337.PubMedView ArticleGoogle Scholar
- Luplertlop N, Surasombatpattana P, Patramool S, Dumas E, Wasinpiyamongkol L, Saune L, Hamel R, Bernard E, Sereno D, Thomas F, Piquemal D, Yssel H, Briant L, Missé D: Induction of a peptide with activity against a broad spectrum of pathogens in the Aedes aegypti salivary gland, following Infection with Dengue Virus. PLoS Pathog. 2011, 7 (1): e1001252-10.1371/journal.ppat.1001252.PubMed CentralPubMedView ArticleGoogle Scholar
- Dimopoulos G: Insect immunity and its implication in mosquito-malaria interactions. Cell Microbiol. 2003, 5 (1): 3-14. 10.1046/j.1462-5822.2003.00252.x.PubMedView ArticleGoogle Scholar
- Alout H, Ndam NT, Sandeu MM, Djegbe I, Chandre F, Dabire RK, Djogbenou LS, Corbel V, Cohuet A: Insecticide resistance alleles affect vector competence of Anopheles gambiae s.s. for Plasmodium falciparum field isolates. PLoS One. 2013, 8 (5): e63849-10.1371/journal.pone.0063849.PubMed CentralPubMedView ArticleGoogle Scholar
- Rivero A, Vezilier J, Weill M, Read AF, Gandon S: Insecticide control of vector-borne diseases: when is insecticide resistance a problem?. PLoS Pathog. 2010, 6 (8): e1001000-10.1371/journal.ppat.1001000.PubMed CentralPubMedView ArticleGoogle Scholar
- Djouaka RF, Bakare AA, Coulibaly ON, Akogbeto MC, Ranson H, Hemingway J, Strode C: Expression of the cytochrome P450s, CYP6P3 and CYP6M2 are significantly elevated in multiple pyrethroid resistant populations of Anopheles gambiae s.s. from Southern Benin and Nigeria. BMC Genomics. 2008, 9: 538-10.1186/1471-2164-9-538.PubMed CentralPubMedView ArticleGoogle Scholar
- Vontas J, David JP, Nikou D, Hemingway J, Christophides GK, Louis C, Ranson H: Transcriptional analysis of insecticide resistance in Anopheles stephensi using cross-species microarray hybridization. Insect Mol Biol. 2007, 16 (3): 315-324. 10.1111/j.1365-2583.2007.00728.x.PubMedView ArticleGoogle Scholar
- David JP, Coissac E, Melodelima C, Poupardin R, Riaz MA, Chandor-Proust A, Reynaud S: Transcriptome response to pollutants and insecticides in the dengue vector Aedes aegypti using next-generation sequencing technology. BMC Genomics. 2010, 11: 216-10.1186/1471-2164-11-216.PubMed CentralPubMedView ArticleGoogle Scholar
- Puinean AM, Foster SP, Oliphant L, Denholm I, Field LM, Millar NS, Williamson MS, Bass C: Amplification of a cytochrome P450 gene is associated with resistance to neonicotinoid insecticides in the aphid Myzus persicae. PLoS Genet. 2010, 6 (6): e1000999-10.1371/journal.pgen.1000999.PubMed CentralPubMedView ArticleGoogle Scholar
- Dittmer NT, Kanost MR: Insect multicopper oxidases: diversity, properties, and physiological roles. Insect Biochem Mol Biol. 2010, 40 (3): 179-188. 10.1016/j.ibmb.2010.02.006.PubMedView ArticleGoogle Scholar
- Gorman MJ, Dittmer NT, Marshall JL, Kanost MR: Characterization of the multicopper oxidase gene family in Anopheles gambiae. Insect Biochem Mol Biol. 2008, 38 (9): 817-824. 10.1016/j.ibmb.2008.07.001.PubMed CentralPubMedView ArticleGoogle Scholar
- Courjaret R, Lapied B: Complex intracellular messenger pathways regulate one type of neuronal alpha-bungarotoxin-resistant nicotinic acetylcholine receptors expressed in insect neurosecretory cells (dorsal unpaired median neurons). Mol Pharmacol. 2001, 60 (1): 80-91.PubMedGoogle Scholar
- Thany SH, Lenaers G, Raymond-Delpech V, Sattelle DB, Lapied B: Exploring the pharmacological properties of insect nicotinic acetylcholine receptors. Trends Pharmacol Sci. 2007, 28 (1): 14-22. 10.1016/j.tips.2006.11.006.PubMedView ArticleGoogle Scholar
- Bodereau-Dubois B, List O, Calas-List D, Marques O, Communal PY, Thany SH, Lapied B: Transmembrane potential polarization, calcium influx, and receptor conformational state modulate the sensitivity of the imidacloprid-insensitive neuronal insect nicotinic acetylcholine receptor to neonicotinoid insecticides. J Pharmacol Exp Ther. 2012, 341 (2): 326-339. 10.1124/jpet.111.188060.PubMedView ArticleGoogle Scholar
- Burmester T, Scheller K: Ligands and receptors: common theme in insect storage protein transport. Naturwissenschaften. 1999, 86 (10): 468-474. 10.1007/s001140050656.PubMedView ArticleGoogle Scholar
- Jinwal UK, Zakharkin SO, Litvinova OV, Jain S, Benes H: Sex-, stage- and tissue-specific regulation by a mosquito hexamerin promoter. Insect Mol Biol. 2006, 15 (3): 301-311. 10.1111/j.1365-2583.2006.00644.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Haunerland NH, Bowers WS: Binding of insecticides to lipophorin and arylphorin, 2 hemolymph-proteins of Heliothis zea. Arch Insect Biochem Physiol. 1986, 3 (1): 87-96. 10.1002/arch.940030110.View ArticleGoogle Scholar
- Buss DS, Callaghan A: Interaction of pesticides with p-glycoprotein and other ABC proteins: a survey of the possible importance to insecticide, herbicide and fungicide resistance. Pest Biochem Physiol. 2008, 90 (3): 141-153. 10.1016/j.pestbp.2007.12.001.View ArticleGoogle Scholar
- Jones CM, Toe HK, Sanou A, Namountougou M, Hughes A, Diabate A, Dabire 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 (9): e45995-10.1371/journal.pone.0045995.PubMed CentralPubMedView ArticleGoogle Scholar
- Chandor-Proust A, Bibby J, Regent-Kloeckner M, Roux J, Guittard-Crilat E, Poupardin R, Riaz MA, Paine M, Dauphin-Villemant C, Reynaud S, David JP: 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Swenson TL, Casida JE: Aldehyde oxidase importance in vivo in xenobiotic metabolism: imidacloprid nitroreduction in mice. Toxicol Sci. 2013, 133 (1): 22-28. 10.1093/toxsci/kft066.PubMedView ArticleGoogle Scholar
- Mitchell SN, Stevenson BJ, Muller P, Wilding CS, Egyir-Yawson A, Field SG, Hemingway J, Paine MJ, 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 U S A. 2012, 109 (16): 6147-6152. 10.1073/pnas.1203452109.PubMed CentralPubMedView ArticleGoogle Scholar
- Stevenson BJ, Bibby J, Pignatelli P, Muangnoicharoen S, O’Neill PM, Lian LY, Muller P, Nikou D, Steven A, Hemingway J, Sutcliffe MJ, Paine MJ: Cytochrome P450 6 M2 from the malaria vector Anopheles gambiae metabolizes pyrethroids: sequential metabolism of deltamethrin revealed. Insect Biochem Mol Biol. 2011, 41 (7): 492-502. 10.1016/j.ibmb.2011.02.003.PubMedView ArticleGoogle Scholar
- Irving H, Riveron JM, Ibrahim SS, Lobo NF, Wondji CS: Positional cloning of rp2 QTL associates the P450 genes CYP6Z1, CYP6Z3 and CYP6M7 with pyrethroid resistance in the malaria vector Anopheles funestus. Heredity. 2012, 109 (6): 383-392. 10.1038/hdy.2012.53.PubMed CentralPubMedView ArticleGoogle Scholar
- Dai D, Tang J, Rose R, Hodgson E, Bienstock RJ, Mohrenweiser HW, Goldstein JA: Identification of variants of CYP3A4 and characterization of their abilities to metabolize testosterone and chlorpyrifos. J Pharmacol Exp Ther. 2001, 299 (3): 825-831.PubMedGoogle Scholar
- Dai D, Zeldin DC, Blaisdell JA, Chanas B, Coulter SJ, Ghanayem BI, Goldstein JA: Polymorphisms in human CYP2C8 decrease metabolism of the anticancer drug paclitaxel and arachidonic acid. Pharmacogenetics. 2001, 11 (7): 597-607. 10.1097/00008571-200110000-00006.PubMedView ArticleGoogle Scholar
- Eiselt R, Domanski TL, Zibat A, Mueller R, Presecan-Siedel E, Hustert E, Zanger UM, Brockmoller J, Klenk HP, Meyer UA, Khan KK, He YA, Halpert JR, Wojnowski L: Identification and functional characterization of eight CYP3A4 protein variants. Pharmacogenetics. 2001, 11 (5): 447-458. 10.1097/00008571-200107000-00008.PubMedView ArticleGoogle Scholar
- Li DN, Seidel A, Pritchard MP, Wolf CR, Friedberg T: Polymorphisms in P450 CYP1B1 affect the conversion of estradiol to the potentially carcinogenic metabolite 4-hydroxyestradiol. Pharmacogenetics. 2000, 10 (4): 343-353. 10.1097/00008571-200006000-00008.PubMedView ArticleGoogle Scholar
- Sata F, Sapone A, Elizondo G, Stocker P, Miller VP, Zheng W, Raunio H, Crespi CL, Gonzalez FJ: CYP3A4 allelic variants with amino acid substitutions in exons 7 and 12: evidence for an allelic variant with altered catalytic activity. Clin Pharmacol Ther. 2000, 67 (1): 48-56. 10.1067/mcp.2000.104391.PubMedView ArticleGoogle Scholar
- Bell CG, Beck S: Advances in the identification and analysis of allele-specific expression. Genome Med. 2009, 1 (5): 56-10.1186/gm56.PubMed CentralPubMedView ArticleGoogle Scholar
- Gaur U, Li K, Mei S, Liu G: Research progress in allele-specific expression and its regulatory mechanisms. J Appl Genet. 2013, 54 (3): 271-283. 10.1007/s13353-013-0148-y.PubMedView ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralPubMedView ArticleGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.PubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J Royal Stat Soc B. 1995, 57: 289-300.Google Scholar
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Annu Stat. 2001, 29 (4): 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, Li J, Thiagarajan M, White JA, Quackenbush J: TM4 microarray software suite. Methods Enzymol. 2006, 411: 134-193.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.