Cis-acting mutation and duplication: History of molecular evolution in a P450 haplotype responsible for insecticide resistance in Culex quinquefasciatus

Whole genome sequence analysis of population structure and insecticide resistance markers in Anopheles melas from the Bijagós Archipelago, Guinea-Bissau

BACKGROUND

Anopheles melas is an understudied malaria vector with a potential role in malaria transmission on the Bijagós Archipelago of Guinea-Bissau. This study presents the first whole-genome sequencing and population genetic analysis for this species from the Bijagós. To our knowledge, this also represents the largest population genetic analysis using WGS data from non-pooled An. melas mosquitoes.

METHODS

WGS was conducted for 30 individual An. melas collected during the peak malaria transmission season in 2019 from six different islands on the Bijagós Archipelago. Bioinformatics tools were used to investigate the population structure and prevalence of insecticide resistance markers in this mosquito population.

RESULTS

Insecticide resistance mutations associated with pyrethroid resistance in Anopheles gambiae s.s. from the Bijagós were absent in the An. melas population, and no signatures of selective sweeps were identified in insecticide resistance-associated genes. Analysis of structural variants identified a large duplication encompassing the cytochrome-P450 gene cyp9k1. Phylogenetic analysis using publicly available mitochondrial genomes indicated that An. melas from the Bijagós split into two phylogenetic groups because of differentiation on the mitochondrial genome attributed to the cytochrome C oxidase subunits COX I and COX II and the NADH dehydrogenase subunits 1, 4, 4L and 5.

CONCLUSIONS

This study identified an absence of insecticide-resistant SNPs common to An. gambiae in the An. melas population, but did identify structural variation over insecticide resistance-associated genes. Furthermore, this study presents novel insights into the population structure of this malaria vector using WGS analysis. Additional studies are required to further understand the role of this vector in malaria transmission.

Genetics, genomics and mechanisms responsible for high levels of pyrethroid resistance in Musca domestica

Insecticide resistance is both economically important and evolutionarily interesting phenomenon. Identification of the mutations responsible for resistance allows for highly sensitive resistance monitoring and allows tools to study the forces (population genetics, fitness costs, etc.) that shape the evolution of resistance. Genes coding for insecticide targets have many well-characterized mutations, but the mutations responsible for enhanced detoxification have proven difficult to identify. We employed multiple strategies to identify the mutations responsible for the extraordinarily high permethrin resistance in the KS17-R strain of house fly (Musca domestica): insecticide synergist assays, linkage analysis, bulk segregant analyses (BSA), transcriptomics and long read DNA (Nanopore) sequencing. The >85,100-fold resistance in KS17-R was partially suppressed by the insecticide synergists piperonyl butoxide and S,S,S-tributylphosphorothionate, but not by diethyl maleate nor by injection. This suggests the involvement of target site insensitivity, CYP-mediated resistance, possibly hydrolase mediated resistance and potentially other unknown factors. Linkage analysis identified chromosomes 1, 2, 3 and 5 as having a role in resistance. BSA mapped resistance loci on chromosomes 3 and 5. The locus on chromosome 3 was centered on the voltage sensitive sodium channel. The locus on chromosome 5 was associated with a duplication of multiple detoxification genes. Transcriptomic analyses and long read DNA sequencing revealed overexpressed CYPs and esterases and identified a complex set of structural variants at the chromosome 5 locus.

Molecular Evolutionary Mechanisms of CYP6ER1vA-Type Variant Associated with Resistance to Neonicotinoid Insecticides in Field Populations of Nilaparvata lugens

The evolution of insecticide resistance has threatened the control of Nilaparvata lugens. Research on mechanisms behind neonicotinoid resistance in N. lugens remains incomplete. This study examined P450-mediated resistance to neonicotinoids in a resistant N. lugens strain (XA-2017-3G). The overexpression of CYP6ER1 in the XA-2017-3G strain plays a role in neonicotinoid resistance, as confirmed by RNA interference. Phenotypic analyses of CYP6ER1-mediated resistance in strains, including laboratory-susceptible, field-collected, and imidacloprid-laboratory further-selected strains, revealed that the vA-type/vL-type genotype exhibited greater resistance to neonicotinoids compared to the vA-type/vA-type genotype. The mRNA expression levels of CYP6ER1vA-type were closely correlated with the levels of neonicotinoid resistance in N. lugens strains, in which CYP6ER1vA-type overexpression is in part attributed to increased copy numbers of CYP6ER1. CYP6ER1vA-type-mediated neonicotinoid resistance was further confirmed by a CYP6ER1vA-type transgenic Drosophila melanogaster line. Taken together, our findings strongly suggest that the overexpression of CYP6ER1vA-type, which can be partially attributed to copy number variations, plays a crucial role in N. lugens resistance to neonicotinoids.

Transposon accumulation at xenobiotic gene family loci in aphids

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

Chromosome-level genome assemblies of two cotton-melon aphid Aphis gossypii biotypes unveil mechanisms of host adaption

The cotton-melon aphid Aphis gossypii is a sap-sucking insect that is considered a serious global pest. The species is distributed over a large geographical range and uses a wide variety of hosts, with some populations being specialized to attack different plant species. Here, we provide de novo chromosome-level genome assemblies of a cotton specialist population (Hap1) and a cucurbit specialist population (Hap3). We achieved this by using a combination of third-generation sequencing platforms, namely Illumina and Hi-C sequencing technologies. We were able to anchor a total of 334.89 Mb (scaffold N50 of 89.13 Mb) and 359.95 Mb (scaffold N50 of 68.88 Mb) to four chromosomes for Hap1 and Hap3, respectively. Moreover, our results showed that the X-chromosome of Hap3 (113.01 Mb) was significantly longer than that of Hap1 (100.26 Mb), with a high level of sequence conservation between the aphid species. We also report variation in the number of protein-coding genes and repeat sequences between Hap1 and Hap3. In particular, olfactory and gustatory receptor genes underwent a high level of gene duplication and expansion events in A. gossypii, including between Hap1 and Hap3. Moreover, we identified two glutathione S-transferase genes which underwent single gene duplications in Hap3, and tandem duplication and inversion events affecting the cytochrome P450 monooxygenase between Hap1 and Hap3, all of which include the CYP3 family. Our results illustrate the variance in the genomic composition of two specialized A. gossypii populations and provide a helpful resource for the study of aphid population evolution, host adaption and insecticide resistance.

High-resolution genetic mapping reveals cis-regulatory and copy number variation in loci associated with cytochrome P450-mediated detoxification in a generalist arthropod pest

Chemical control strategies are driving the evolution of pesticide resistance in pest populations. Understanding the genetic mechanisms of these evolutionary processes is of crucial importance to develop sustainable resistance management strategies. The acaricide pyflubumide is one of the most recently developed mitochondrial complex II inhibitors with a new mode of action that specifically targets spider mite pests. In this study, we characterize the molecular basis of pyflubumide resistance in a highly resistant population of the spider mite Tetranychus urticae. Classical genetic crosses indicated that pyflubumide resistance was incompletely recessive and controlled by more than one gene. To identify resistance loci, we crossed the resistant population to a highly susceptible T. urticae inbred strain and propagated resulting populations with and without pyflubumide exposure for multiple generations in an experimental evolution set-up. High-resolution genetic mapping by a bulked segregant analysis approach led to the identification of three quantitative trait loci (QTL) linked to pyflubumide resistance. Two QTLs were found on the first chromosome and centered on the cytochrome P450 CYP392A16 and a cluster of CYP392E6-8 genes. Comparative transcriptomics revealed a consistent overexpression of CYP392A16 and CYP392E8 in the experimental populations that were selected for pyflubumide resistance. We further corroborated the involvement of CYP392A16 in resistance by in vitro functional expression and metabolism studies. Collectively, these experiments uncovered that CYP392A16 N-demethylates the toxic carboxamide form of pyflubumide to a non-toxic compound. A third QTL coincided with cytochrome P450 reductase (CPR), a vital component of cytochrome P450 metabolism. We show here that the resistant population harbors three gene copies of CPR and that this copy number variation is associated with higher mRNA abundance. Together, we provide evidence for detoxification of pyflubumide by cytochrome P450s that is likely synergized by gene amplification of CPR.

Our understanding of the causal genetic variants that drive the evolution of quantitative traits, such as polygenic pesticide resistance, remains very limited. Here, we followed a high-resolution genetic mapping approach to localize the genetic variants that cause pyflubumide resistance in the two-spotted spider mite Tetranychus urticae. Three well-supported QTL were uncovered and pointed towards a major role for cytochrome P450-mediated detoxification. Cis-regulatory variation for cytochrome P450s was observed, and in vitro cytochrome P450 experiments showed that pyflubumide was metabolized into a non-toxic derivate. A third QTL centered on cytochrome P450 reductase (CPR), which is required for cytochrome P450 activity, and is amplified in pyflubumide resistant populations. Our results indicate that pyflubumide resistance is mediated by cytochrome P450 detoxification that is enhanced by gene amplification at the CPR locus.

Multiple cytochrome P450 genes: conferring high levels of permethrin resistance in mosquitoes, Culex quinquefasciatus

Insecticides, especially pyrethroids, are the most important in the insect pest control and preventing insect vector-borne human diseases. However, insect pests, including mosquitoes, have developed resistance in the insecticides that used against them. Cytochrome P450s are associated with insecticide resistance through overexpression and detoxification mechanisms in insect species. In this study, we utilized a powerful tool, the RNAi technique, to determine the roles of key P450 genes overexpressed in permethrin resistant mosquitoes that confer insecticide resistance to unravel the molecular basis of resistance mechanisms in the mosquito Culex quinquefasciatus. The results showed that knockdown of 8 key P450 genes using RNAi techniques significantly decreased resistance to permethrin in resistant mosquitoes. In silico modeling and docking analysis further revealed the potential metabolic function of overexpressed P450 genes in the development of insecticide resistance in mosquitoes. These findings not only highlighted the functional importance of these P450 genes in insecticide resistance, but also revealed that overexpression of multiple P450 genes was responsible for the high levels of insecticide resistance in a mosquito population of Culex quinquefasciatus.

From Global to Local-New Insights into Features of Pyrethroid Detoxification in Vector Mosquitoes

The threat of mosquito-borne diseases continues to be a problem for public health in subtropical and tropical regions of the world; in response, there has been increased use of adulticidal insecticides, such as pyrethroids, in human habitation areas over the last thirty years. As a result, the prevalence of pyrethroid-resistant genetic markers in natural mosquito populations has increased at an alarming rate. This review details recent advances in the understanding of specific mechanisms associated with pyrethroid resistance, with emphasis on features of insecticide detoxification and the interdependence of multiple cellular pathways. Together, these advances add important context to the understanding of the processes that are selected in resistant mosquitoes. Specifically, before pyrethroids bind to their targets on motoneurons, they must first permeate the outer cuticle and diffuse to inner tissues. Resistant mosquitoes have evolved detoxification mechanisms that rely on cytochrome P450s (CYP), esterases, carboxyesterases, and other oxidation/reduction (redox) components to effectively detoxify pyrethroids to nontoxic breakdown products that are then excreted. Enhanced resistance mechanisms have evolved to include alteration of gene copy number, transcriptional and post-transcriptional regulation of gene expression, as well as changes to cellular signaling mechanisms. Here, we outline the variety of ways in which detoxification has been selected in various mosquito populations, as well as key gene categories involved. Pathways associated with potential new genes of interest are proposed. Consideration of multiple cellular pathways could provide opportunities for development of new insecticides.

Novel CONCOMITANT mutations L932F and I936V in the Voltage-Gated Sodium Channel and Its Association With Pyrethroid Resistance in Culex quinquefasciatus (Diptera: Culicidae)

Highly residual pyrethroids such as permethrin have been used for controlling mosquitoes that transmit infectious diseases. However, the selective pressure from such insecticides may result in cross-resistance against other pyrethroids used for household insecticides. In this study, we investigated the susceptibility of Culex quinquefasciatus Say collected from Brazil and Myanmar to permethrin in addition to four types of household pyrethroids. Both strains exhibited high resistance against all pyrethroids tested, indicating cross-resistance. Furthermore, we detected the knockdown resistance (kdr) mutations L932F+I936V in the voltage-gated sodium channel gene (VGSC) in the Brazilian strain. Notably, the L932F+I936V haplotype has previously been observed in in silico data, but it should be detected not directly from living insects. In comparison, a common kdr mutation, L1014F, was detected from the Myanmar strain. Although L1014F was also detected from the Brazilian strain, the allele frequency was too low to affect resistance. Both strains harbored the resistance-associated haplotypes of the cytochrome P450 gene, CYP9M10. The Brazilian strain demonstrated comparable resistance against pyrethroids as that of the Myanmar strain even when a cytochrome P450 inhibitor, piperonyl butoxide was added to the bioassay. Our results suggested that the L932F+I936V mutations confer the Brazilian strain of Cx. Quiquefasciatus with resistance at a comparable level to that conferred by the well-recognized kdr mutation L1014F in the Myanmar strain. The identification of unexplored mutations may improve the diagnosis and understanding of resistance of this medically important species.

Chromatin Structure and Function in Mosquitoes

The principles and function of chromatin and nuclear architecture have been extensively studied in model organisms, such as Drosophila melanogaster. However, little is known about the role of these epigenetic processes in transcriptional regulation in other insects including mosquitoes, which are major disease vectors and a worldwide threat for human health. Some of these life-threatening diseases are malaria, which is caused by protozoan parasites of the genus Plasmodium and transmitted by Anopheles mosquitoes; dengue fever, which is caused by an arbovirus mainly transmitted by Aedes aegypti; and West Nile fever, which is caused by an arbovirus transmitted by Culex spp. In this contribution, we review what is known about chromatin-associated mechanisms and the 3D genome structure in various mosquito vectors, including Anopheles, Aedes, and Culex spp. We also discuss the similarities between epigenetic mechanisms in mosquitoes and the model organism Drosophila melanogaster, and advocate that the field could benefit from the cross-application of state-of-the-art functional genomic technologies that are well-developed in the fruit fly. Uncovering the mosquito regulatory genome can lead to the discovery of unique regulatory networks associated with the parasitic life-style of these insects. It is also critical to understand the molecular interactions between the vectors and the pathogens that they transmit, which could hold the key to major breakthroughs on the fight against mosquito-borne diseases. Finally, it is clear that epigenetic mechanisms controlling mosquito environmental plasticity and evolvability are also of utmost importance, particularly in the current context of globalization and climate change.

A cytochrome P450 allele confers pyrethroid resistance on a major African malaria vector, reducing insecticide-treated bednet efficacy

Metabolic resistance to insecticides such as pyrethroids in mosquito vectors threatens control of malaria in Africa. Unless it is managed, recent gains in reducing malaria transmission could be lost. To improve monitoring and assess the impact of insecticide resistance on malaria control interventions, we elucidated the molecular basis of pyrethroid resistance in the major African malaria vector, Anopheles funestus We showed that a single cytochrome P450 allele (CYP6P9a_R) in A. funestus reduced the efficacy of insecticide-treated bednets for preventing transmission of malaria in southern Africa. Expression of key insecticide resistance genes was detected in populations of this mosquito vector throughout Africa but varied according to the region. Signatures of selection and adaptive evolutionary traits including structural polymorphisms and cis-regulatory transcription factor binding sites were detected with evidence of selection due to the scale-up of insecticide-treated bednet use. A cis-regulatory polymorphism driving the overexpression of the major resistance gene CYP6P9a allowed us to design a DNA-based assay for cytochrome P450-mediated resistance to pyrethroid insecticides. Using this assay, we tracked the spread of pyrethroid resistance and found that it was almost fixed in mosquitoes from southern Africa but was absent from mosquitoes collected elsewhere in Africa. Furthermore, a field study in experimental huts in Cameroon demonstrated that mosquitoes carrying the resistance CYP6P9a_R allele survived and succeeded in blood feeding more often than did mosquitoes that lacked this allele. Our findings highlight the need to introduce a new generation of insecticide-treated bednets for malaria control that do not rely on pyrethroid insecticides.

High-throughput genotyping of a full voltage-gated sodium channel gene via genomic DNA using target capture sequencing and analytical pipeline MoNaS to discover novel insecticide resistance mutations

In insects, the voltage-gated sodium channel (VGSC) is the primary target site of pyrethroid insecticides. Various amino acid substitutions in the VGSC protein, which are selected under insecticide pressure, are known to confer insecticide resistance. In the genome, the VGSC gene consists of more than 30 exons sparsely distributed across a large genomic region, which often exceeds 100 kbp. Due to this complex genomic structure, it is often challenging to genotype full coding nucleotide sequences (CDSs) of VGSC from individual genomic DNA (gDNA). In this study, we designed biotinylated oligonucleotide probes from CDSs of VGSC of Asian tiger mosquito, Aedes albopictus. The probe set effectively concentrated (>80,000-fold) all targeted regions of gene VGSC from pooled barcoded Illumina libraries each constructed from individual A. albopictus gDNAs. The probe set also captured all orthologous VGSC CDSs, except some tiny exons, from the gDNA of other Culicinae mosquitos, A. aegypti and Culex pipiens complex, with comparable efficiency as a result of the high nucleotide-level conservation of VGSC. To improve efficiency of the downstream bioinformatic process, we developed an automated pipeline-MoNaS (Mosquito Na+ channel mutation Search)-which calls amino acid substitutions in the VGSC from NGS reads and compares those to known resistance mutations. The proposed method and our bioinformatic tool should facilitate the discovery of novel amino acid variants conferring insecticide resistance on VGSC and population genetic studies on resistance alleles (with respect to the origin, selection, and migration etc.) in both clinically and agriculturally important insect pests.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

Copy number variation (CNV) and insecticide resistance in mosquitoes: evolving knowledge or an evolving problem?

Copy number variation (CNV) in insect genomes is a rich source of potentially adaptive polymorphism which may help overcome the constraints of purifying selection on conserved genes and/or permit elevated transcription. Classic studies of amplified esterases and acetylcholinesterase duplication in Culex pipiens quantified evolutionary dynamics of CNV driven by insecticidal selection. A more complex and potentially medically impactful form of CNV is found in Anopheles gambiae, with both heterogeneous duplications and homogeneous amplifications strongly linked with insecticide resistance. Metabolic gene amplification, revealed by shotgun sequencing, appears common in Aedes aegypti, but poorly understood in other mosquito species. Many methodologies have been used to detect CNV in mosquitoes, but relatively few can detect both duplications and amplifications, and contrasting methods should be combined. Genome scans for CNV have been rare to date in mosquitoes, but offer immense potential to determine the overall role of CNV as a component of resistance mechanisms.

The Glutathione-S-Transferase, Cytochrome P450 and Carboxyl/Cholinesterase Gene Superfamilies in Predatory Mite Metaseiulus occidentalis

Pesticide-resistant populations of the predatory mite Metaseiulus (= Typhlodromus or Galendromus) occidentalis (Arthropoda: Chelicerata: Acari: Phytoseiidae) have been used in the biological control of pest mites such as phytophagous Tetranychus urticae. However, the pesticide resistance mechanisms in M. occidentalis remain largely unknown. In other arthropods, members of the glutathione-S-transferase (GST), cytochrome P450 (CYP) and carboxyl/cholinesterase (CCE) gene superfamilies are involved in the diverse biological pathways such as the metabolism of xenobiotics (e.g. pesticides) in addition to hormonal and chemosensory processes. In the current study, we report the identification and initial characterization of 123 genes in the GST, CYP and CCE superfamilies in the recently sequenced M. occidentalis genome. The gene count represents a reduction of 35% compared to T. urticae. The distribution of genes in the GST and CCE superfamilies in M. occidentalis differs significantly from those of insects and resembles that of T. urticae. Specifically, we report the presence of the Mu class GSTs, and the J’ and J” clade CCEs that, within the Arthropoda, appear unique to Acari. Interestingly, the majority of CCEs in the J’ and J” clades contain a catalytic triad, suggesting that they are catalytically active. They likely represent two Acari-specific CCE clades that may participate in detoxification of xenobiotics. The current study of genes in these superfamilies provides preliminary insights into the potential molecular components that may be involved in pesticide metabolism as well as hormonal/chemosensory processes in the agriculturally important M. occidentalis.

Testing the causality between CYP9M10 and pyrethroid resistance using the TALEN and CRISPR/Cas9 technologies

Recently-emerging genome editing technologies have enabled targeted gene knockout experiments even in non-model insect species. For studies on insecticide resistance, genome editing technologies offer some advantages over the conventional reverse genetic technique, RNA interference, for testing causal relationships between genes of detoxifying enzymes and resistance phenotypes. There were relatively abundant evidences indicating that the overexpression of a cytochrome P450 gene CYP9M10 confers strong pyrethroid resistance in larvae of the southern house mosquito Culex quinquefasciatus. However, reverse genetic verification has not yet been obtained because of the technical difficulty of microinjection into larvae. Here, we tested two genome editing technologies, transcription activator-like effector nucleases (TALEN)s and clustered regularly interspaced short palindromic repeats (CRISPR/Cas9), to disrupt CYP9M10 in a resistant strain of C. quinquefasciatus. Additionally, we developed a novel, effective approach to construct a TALE using the chemical cleavage of phosphorothioate inter-nucleotide linkages in the level 1 assembly. Both TALEN and CRISPR/Cas9 induced frame-shifting mutations in one or all copies of CYP9M10 in a pyrethroid-resistant strain. A line fixed with a completely disrupted CYP9M10 haplotype showed more than 100-fold reduction in pyrethroid resistance in the larval stage.

A single nucleotide change in a core promoter is involved in the progressive overexpression of the duplicated CYP9M10 haplotype lineage in Culex quinquefasciatus

Although the importance of cis-acting mutations on detoxification enzyme genes for insecticide resistance is widely accepted, only a few of them have been determined as concrete mutations present in genomic DNA till date. The overexpression of a cytochrome P450 gene, CYP9M10, is associated with pyrethroid resistance in the southern house mosquito Culex quinquefasciatus. The haplotypes of CYP9M10 exhibiting overexpression (resistant haplotypes) belong to one specific phylogenetic lineage that shares high nucleotide sequence homology and the same insertion of a transposable element. Among the resistant haplotypes, allelic progression involving an additional cis-acting mutation and gene duplication evolved a CYP9M10 haplotype associated with extremely high transcription and strong pyrethroid resistance. Here we show that a single nucleotide substitution G-27A, which is located near the transcription start site of CYP9M10, is involved in the progression of the duplicated haplotype lineage. The deletion of a 7-bp AT-rich sequence that includes nucleotide -27 inhibited the initiation of transcription from the original transcriptional initiation site. The mutation was suspected to reside within a core promoter, TATA-box, of CYP9M10.

Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing

The capacity of mosquitoes to resist insecticides threatens the control of diseases such as dengue and malaria. Until alternative control tools are implemented, characterizing resistance mechanisms is crucial for managing resistance in natural populations. Insecticide biodegradation by detoxification enzymes is a common resistance mechanism; however, the genomic changes underlying this mechanism have rarely been identified, precluding individual resistance genotyping. In particular, the role of copy number variations (CNVs) and polymorphisms of detoxification enzymes have never been investigated at the genome level, although they can represent robust markers of metabolic resistance. In this context, we combined target enrichment with high-throughput sequencing for conducting the first comprehensive screening of gene amplifications and polymorphisms associated with insecticide resistance in mosquitoes. More than 760 candidate genes were captured and deep sequenced in several populations of the dengue mosquito Ae. aegypti displaying distinct genetic backgrounds and contrasted resistance levels to the insecticide deltamethrin. CNV analysis identified 41 gene amplifications associated with resistance, most affecting cytochrome P450s overtranscribed in resistant populations. Polymorphism analysis detected more than 30,000 variants and strong selection footprints in specific genomic regions. Combining Bayesian and allele frequency filtering approaches identified 55 nonsynonymous variants strongly associated with resistance. Both CNVs and polymorphisms were conserved within regions but differed across continents, confirming that genomic changes underlying metabolic resistance to insecticides are not universal. By identifying novel DNA markers of insecticide resistance, this study opens the way for tracking down metabolic changes developed by mosquitoes to resist insecticides within and among populations.

Genotype to phenotype, the molecular and physiological dimensions of resistance in arthropods

The recent accumulation of molecular studies on mutations in insects, ticks and mites conferring resistance to insecticides, acaricides and biopesticides is reviewed. Resistance is traditionally classified by physiological and biochemical criteria, such as target-site insensitivity and metabolic resistance. However, mutations are discrete molecular changes that differ in their intrinsic frequency, effects on gene dosage and fitness consequences. These attributes in turn impact the population genetics of resistance and resistance management strategies, thus calling for a molecular genetic classification. Mutations in structural genes remain the most abundantly described, mostly in genes coding for target proteins. These provide the most compelling examples of parallel mutations in response to selection. Mutations causing upregulation and downregulation of genes, both in cis (in the gene itself) and in trans (in regulatory processes) remain difficult to characterize precisely. Gene duplications and gene disruption are increasingly reported. Gene disruption appears prevalent in the case of multiple, hetero-oligomeric or redundant targets.

Pyrethroid resistance in Culex pipiens mosquitoes

Mosquitoes within the Culex pipiens complex are widely distributed and important in the transmission of many human diseases. Insecticides, pyrethroids in particular, remain a mainstay for control of these important vectors. In this paper we review what is known about the levels, mechanisms and fitness costs of pyrethroid resistance in Cx. pipiens. Pyrethroid resistance in Cx. pipiens is a global problem, and resistance ratios of up to 7000-fold have been found in larvae of field collected mosquitoes. However, there is considerable variation between populations, indicating significant geographic heterogeneity of the resistance. The two major mechanisms of resistance to pyrethroids in Culex are mutations in Vssc (target site insensitivity) and overexpression of cytochrome P450(s) (increased detoxification). The most frequently reported Vssc mutation is L1014F (i.e. kdr), which has been found throughout the world. The L1014S mutation has been found in Cx. p. pallens from Japan and China, and in Cx. p. pipiens from China. The L1014C mutation has only been reported for Cx. p. pipens molestus from China and the V1016G mutation has only been reported from Saudi Arabia. Studies on the P450s of Cx. pipiens have identified several that are overexpressed (measured as transcript levels) in pyrethroid resistant strains. CYP9M10 is consistently overexpressed in pyrethroid resistant Cx. pipiens from at least seven countries, suggesting this P450 might be of global importance in resistance. Both CYP9M10-mediated pyrethroid resistance and kdr have fitness costs in the absence of insecticides under certain environmental conditions. Research needs and future directions are discussed.

Mechanisms of pyrethroid resistance in the dengue mosquito vector, Aedes aegypti: target site insensitivity, penetration, and metabolism

Aedes aegypti is the major vector of yellow and dengue fevers. After 10 generations of adult selection, an A. aegypti strain (SP) developed 1650-fold resistance to permethrin, which is one of the most widely used pyrethroid insecticides for mosquito control. SP larvae also developed 8790-fold resistance following selection of the adults. Prior to the selections, the frequencies of V1016G and F1534C mutations in domains II and III, respectively, of voltage-sensitive sodium channel (Vssc, the target site of pyrethroid insecticide) were 0.44 and 0.56, respectively. In contrast, only G1016 alleles were present after two permethrin selections, indicating that G1016 can more contribute to the insensitivity of Vssc than C1534. In vivo metabolism studies showed that the SP strain excreted permethrin metabolites more rapidly than a susceptible SMK strain. Pretreatment with piperonyl butoxide caused strong inhibition of excretion of permethrin metabolites, suggesting that cytochrome P450 monooxygenases (P450s) play an important role in resistance development. In vitro metabolism studies also indicated an association of P450s with resistance. Microarray analysis showed that multiple P450 genes were over expressed during the larval and adult stages in the SP strain. Following quantitative real time PCR, we focused on two P450 isoforms, CYP9M6 and CYP6BB2. Transcription levels of these P450s were well correlated with the rate of permethrin excretion and they were certainly capable of detoxifying permethrin to 4'-HO-permethrin. Over expression of CYP9M6 was partially due to gene amplification. There was no significant difference in the rate of permethrin reduction from cuticle between SP and SMK strains.

Gene amplification and microsatellite polymorphism underlie a recent insect host shift

Host plant shifts of herbivorous insects may be a first step toward sympatric speciation and can create new pests of agriculturally important crops; however, the molecular mechanisms that mediate this process are poorly understood. Certain races of the polyphagous aphid Myzus persicae have recently adapted to feed on tobacco (Myzus persicae nicotianae) and show a reduced sensitivity to the plant alkaloid nicotine and cross-resistance to neonicotinoids a class of synthetic insecticides widely used for control. Here we show constitutive overexpression of a cytochrome P450 (CYP6CY3) allows tobacco-adapted races of M. persicae to efficiently detoxify nicotine and has preadapted them to resist neonicotinoid insecticides. CYP6CY3, is highly overexpressed in M. persicae nicotianae clones from three continents compared with M. persicae s.s. and expression level is significantly correlated with tolerance to nicotine. CYP6CY3 is highly efficient (compared with the primary human nicotine-metabolizing P450) at metabolizing nicotine and neonicotinoids to less toxic metabolites in vitro and generation of transgenic Drosophila expressing CYP6CY3 demonstrate that it confers resistance to both compounds in vivo. Overexpression of CYP6CY3 results from the expansion of a dinucleotide microsatellite in the promoter region and a recent gene amplification, with some aphid clones carrying up to 100 copies. We conclude that the mutations leading to overexpression of CYP6CY3 were a prerequisite for the host shift of M. persicae to tobacco and that gene amplification and microsatellite polymorphism are evolutionary drivers in insect host adaptation.

Global spread and genetic variants of the two CYP9M10 haplotype forms associated with insecticide resistance in Culex quinquefasciatus Say

Insecticide resistance develops as a genetic factor (allele) conferring lower susceptibility to insecticides proliferates within a target insect population under strong positive selection. Intriguingly, a resistance allele pre-existing in a population often bears a series of further adaptive allelic variants through new mutations. This phenomenon occasionally results in replacement of the predominating resistance allele by fitter new derivatives, and consequently, development of greater resistance at the population level. The overexpression of the cytochrome P450 gene CYP9M10 is associated with pyrethroid resistance in the southern house mosquito Culex quinquefasciatus. Previously, we have found two genealogically related overexpressing CYP9M10 haplotypes, which differ in gene copy number (duplicated and non-duplicated). The duplicated haplotype was derived from the non-duplicated overproducer probably recently. In the present study, we investigated allelic series of CYP9M10 involved in three C. quinquefasciatus laboratory colonies recently collected from three different localities. Duplicated and non-duplicated overproducing haplotypes coexisted in African and Asian colonies indicating a global distribution of both haplotype lineages. The duplicated haplotypes both in the Asian and African colonies were associated with higher expression levels and stronger resistance than non-duplicated overproducing haplotypes. There were slight variation in expression level among the non-duplicated overproducing haplotypes. The nucleotide sequences in coding and upstream regions among members of this group also showed a little diversity. Non-duplicated overproducing haplotypes with relatively higher expression were genealogically closer to the duplicated haplotypes than the other non-duplicated overproducing haplotypes, suggesting multiple cis-acting mutations before duplication.

Impact of environment on mosquito response to pyrethroid insecticides: facts, evidences and prospects

By transmitting major human diseases such as malaria, dengue fever and filariasis, mosquito species represent a serious threat worldwide in terms of public health, and pose a significant economic burden for the African continent and developing tropical regions. Most vector control programmes aiming at controlling life-threatening mosquitoes rely on the use of chemical insecticides, mainly belonging to the pyrethroid class. However, resistance of mosquito populations to pyrethroids is increasing at a dramatic rate, threatening the efficacy of control programmes throughout insecticide-treated areas, where mosquito-borne diseases are still prevalent. In the absence of new insecticides and efficient alternative vector control methods, resistance management strategies are therefore critical, but these require a deep understanding of adaptive mechanisms underlying resistance. Although insecticide resistance mechanisms are intensively studied in mosquitoes, such adaptation is often considered as the unique result of the selection pressure caused by insecticides used for vector control. Indeed, additional environmental parameters, such as insecticides/pesticides usage in agriculture, the presence of anthropogenic or natural xenobiotics, and biotic interactions between vectors and other organisms, may affect both the overall mosquito responses to pyrethroids and the selection of resistance mechanisms. In this context, the present work aims at updating current knowledge on pyrethroid resistance mechanisms in mosquitoes and compiling available data, often from different research fields, on the impact of the environment on mosquito response to pyrethroids. Key environmental factors, such as the presence of urban or agricultural pollutants and biotic interactions between mosquitoes and their microbiome are discussed, and research perspectives to fill in knowledge gaps are suggested.

Directionally selected cytochrome P450 alleles are driving the spread of pyrethroid resistance in the major malaria vector Anopheles funestus

Pyrethroid insecticides are critical for malaria control in Africa. However, resistance to this insecticide class in the malaria vector Anopheles funestus is spreading rapidly across Africa, threatening the success of ongoing and future malaria control programs. The underlying resistance mechanisms driving the spread of this resistance in wild populations remain largely unknown. Here, we show that increased expression of two tandemly duplicated P450 genes, CYP6P9a and CYP6P9b, is the main mechanism driving pyrethroid resistance in Malawi and Mozambique, two southern African countries where this insecticide class forms the mainstay of malaria control. Genome-wide transcription analysis using microarray and quantitative RT-PCR consistently revealed that CYP6P9a and CYP6P9b are the two genes most highly overexpressed (>50-fold; q < 0.01) in permethrin-resistant mosquitoes. Transgenic expression of CYP6P9a and CYP6P9b in Drosophila melanogaster demonstrated that elevated expression of either of these genes confers resistance to both type I (permethrin) and type II (deltamethrin) pyrethroids. Functional characterization of recombinant CYP6P9b confirmed that this protein metabolized both type I (permethrin and bifenthrin) and type II (deltamethrin and Lambda-cyhalothrin) pyrethroids but not DDT. Variability analysis identified that a single allele of each of these genes is predominantly associated with pyrethroid resistance in field populations from both countries, which is suggestive of a single origin of this resistance that has since spread across the region. Urgent resistance management strategies should be implemented in this region to limit a further spread of this resistance and minimize its impact on the success of ongoing malaria control programs.

Resistance of Australian Helicoverpa armigera to fenvalerate is due to the chimeric P450 enzyme CYP337B3

Worldwide, increasing numbers of insects have evolved resistance to a wide range of pesticides, which hampers their control in the field and, therefore, threatens agriculture. Members of the carboxylesterase and cytochrome P450 monooxygenase superfamilies are prominent candidates to confer metabolic resistance to pyrethroid insecticides. Both carboxylesterases and P450 enzymes have been shown to be involved in pyrethroid resistance in Australian Helicoverpa armigera, the noctuid species possessing by far the most reported resistance cases worldwide. However, specific enzymes responsible for pyrethroid resistance in field populations of this species have not yet been identified. Here, we show that the resistance toward fenvalerate in an Australian strain of H. armigera is due to a unique P450 enzyme, CYP337B3, which arose from unequal crossing-over between two parental P450 genes, resulting in a chimeric enzyme. CYP337B3 is capable of metabolizing fenvalerate into 4'-hydroxyfenvalerate, which exhibits no toxic effect on susceptible larvae; enzymes from the parental P450 genes showed no detectable fenvalerate metabolism. Furthermore, a polymorphic H. armigera strain could be bred into a susceptible line possessing the parental genes CYP337B1 and CYP337B2 and a resistant line possessing only CYP337B3. The exclusive presence of CYP337B3 in resistant insects of this strain confers a 42-fold resistance to fenvalerate. Thus, in addition to previously documented genetic mechanisms of resistance, recombination can also generate selectively advantageous variants, such as this chimeric P450 enzyme with an altered substrate specificity leading to a potent resistance mechanism.

A cis-regulatory sequence driving metabolic insecticide resistance in mosquitoes: functional characterisation and signatures of selection

Although cytochrome P450 (CYP450) enzymes are frequently up-regulated in mosquitoes resistant to insecticides, no regulatory motifs driving these expression differences with relevance to wild populations have been identified. Transposable elements (TEs) are often enriched upstream of those CYP450s involved in insecticide resistance, leading to the assumption that they contribute regulatory motifs that directly underlie the resistance phenotype. A partial CuRE1 (Culex Repetitive Element 1) transposable element is found directly upstream of CYP9M10, a cytochrome P450 implicated previously in larval resistance to permethrin in the ISOP450 strain of Culex quinquefasciatus, but is absent from the equivalent genomic region of a susceptible strain. Via expression of CYP9M10 in Escherichia coli we have now demonstrated time- and NADPH-dependant permethrin metabolism, prerequisites for confirmation of a role in metabolic resistance, and through qPCR shown that CYP9M10 is >20-fold over-expressed in ISOP450 compared to a susceptible strain. In a fluorescent reporter assay the region upstream of CYP9M10 from ISOP450 drove 10× expression compared to the equivalent region (lacking CuRE1) from the susceptible strain. Close correspondence with the gene expression fold-change implicates the upstream region including CuRE1 as a cis-regulatory element involved in resistance. Only a single CuRE1 bearing allele, identical to the CuRE1 bearing allele in the resistant strain, is found throughout Sub-Saharan Africa, in contrast to the diversity encountered in non-CuRE1 alleles. This suggests a single origin and subsequent spread due to selective advantage. CuRE1 is detectable using a simple diagnostic. When applied to C. quinquefasciatus larvae from Ghana we have demonstrated a significant association with permethrin resistance in multiple field sites (mean Odds Ratio = 3.86) suggesting this marker has relevance to natural populations of vector mosquitoes. However, when CuRE1 was excised from the allele used in the reporter assay through fusion PCR, expression was unaffected, indicating that the TE has no direct role in resistance and hence that CuRE1 is acting only as a marker of an as yet unidentified regulatory motif in the association analysis. This suggests that a re-evaluation of the assumption that TEs contribute regulatory motifs involved in gene expression may be necessary.

Molecular cloning and preliminary function study of iron responsive element binding protein 1 gene from cypermethrin-resistant Culex pipiens pallens

Background

Insecticide resistance jeopardizes the control of mosquito populations and mosquito-borne disease control, which creates a major public health concern. Two-dimensional electrophoresis identified one protein segment with high sequence homology to part of Aedes aegypti iron-responsive element binding protein (IRE-BP).

Method

RT-PCR and RACE (rapid amplification of cDNA end) were used to clone a cDNA encoding full length IRE-BP 1. Real-time quantitative RT-PCR was used to evaluate the transcriptional level changes in the Cr-IRE strain Aedes aegypti compared to the susceptible strain of Cx. pipiens pallens. The expression profile of the gene was established in the mosquito life cycle. Methyl tritiated thymidine (³H-TdR) was used to observe the cypermethrin resistance changes in C6/36 cells containing the stably transfected IRE-BP 1 gene of Cx. pipiens pallens.

Results

The complete sequence of iron responsive element binding protein 1 (IRE-BP 1) has been cloned from the cypermethrin-resistant strain of Culex pipiens pallens (Cr-IRE strain). Quantitative RT-PCR analysis indicated that the IRE-BP 1 transcription level was 6.7 times higher in the Cr-IRE strain than in the susceptible strain of 4th instar larvae. The IRE-BP 1 expression was also found to be consistently higher throughout the life cycle of the Cr-IRE strain. A protein of predicted size 109.4 kDa has been detected by Western blotting in IRE-BP 1-transfected mosquito C6/36 cells. These IRE-BP 1-transfected cells also showed enhanced cypermethrin resistance compared to null-transfected or plasmid vector-transfected cells as determined by ³H-TdR incorporation.

Conclusion

IRE-BP 1 is expressed at higher levels in the Cr-IRE strain, and may confer some insecticide resistance in Cx. pipiens pallens.