Why are there so few resistance-associated mutations in insecticide target genes?

Adipose Tissue Exosome circ_sxc Mediates the Modulatory of Adiposomes on Brain Aging by Inhibiting Brain dme-miR-87-3p

Aging of the brain usually leads to the decline of neurological processes and is a major risk factor for various neurodegenerative diseases, including sleep disturbances and cognitive decline. Adipose tissue exosomes, as adipocyte-derived vesicles, may mediate the regulatory processes of adipose tissue on other organs, including the brain; however, the regulatory mechanisms remain unclear. We analyzed the sleep-wake behavior of young (10 days) and old (40 days) Drosophila and found that older Drosophila showed increased sleep fragmentation, which is similar to mammalian aging characteristics. To investigate the cross-tissue regulatory mechanisms of adiposity on brain aging, we extracted 10-day and 40-day Drosophila adipose tissue exosomes and identified circRNAs with age-dependent expression differences by RNA-seq and differential analysis. Furthermore, by combining data from 3 datasets of the GEO database (GSE130158, GSE24992, and GSE184559), circ_sxc that was significantly downregulated with age was finally screened out. Moreover, dme-miR-87-3p, a conserved target of circ_sxc, accumulates in the brain with age and exhibits inhibitory effects in predicted binding relationships with neuroreceptor ligand genes. In summary, the current study showed that the Drosophila brain could obtain circ_sxc by uptake of adipose tissue exosomes which crossed the blood-brain barrier. And circ_sxc suppressed brain miR-87-3p expression through sponge adsorption, which in turn regulated the expression of neurological receptor ligand proteins (5-HT1B, GABA-B-R1, Rdl, Rh7, qvr, NaCP60E) and ensured brain neuronal synaptic signaling normal function of synaptic signaling. However, with aging, this regulatory mechanism is dysregulated by the downregulation of the adipose exosome circ_sxc, which contributes to the brain exhibiting sleep disturbances and other "aging" features.

Molecular diversity of genes related to biological rhythms (period and timeless) and insecticide resistance (Na V and ace-1) in Anopheles darlingi

BACKGROUND

Malaria is a public health concern in the Amazonian Region, where Anopheles darlingi is the main vector of Plasmodium spp. Several studies hypothesised the existence of cryptic species in An. darlingi, considering variations in behaviour, morphological and genetic aspects. Determining their overall genetic background for vector competence, insecticide resistance, and other elements is essential to better guide strategies for malaria control.

OBJECTIVES

This study aimed to evaluate the molecular diversity in genes related to behaviour and insecticide resistance, estimating genetic differentiation in An. darlingi populations from Amazonian localities in Brazil and Pacific Colombian region.

METHODS

We amplified, cloned and sequenced fragments of genes related to behaviour: timeless (tim) and period (per), and to insecticide resistance: voltage-gated sodium channel (Na V ) and acetylcholinesterase (ace-1) from 516 An. darlingi DNA samples from Manaus, Unini River, Jaú River and Porto Velho - Brazil, and Chocó - Colombia. We discriminated single nucleotide polymorphisms (SNPs), determined haplotypes and evaluate the phylogenetic relationship among the populations.

FINDINGS

The genes per, tim and ace-1 were more polymorphic than Na V . The classical kdr and ace-1 R mutations were not observed. Phylogenetic analyses suggested a significant differentiation between An. darlingi populations from Brazil and Colombia, except for the Na V gene. There was a geographic differentiation within Brazilian populations considering per and ace-1.

CONCLUSIONS

Our results add genetic data to the discussion about polymorphisms at population levels in An. darlingi. The search for insecticide resistance-related mechanisms should be extended to more populations, especially from localities with a vector control failure scenario.

Predicting the insecticide-driven mutations in a crop pest insect: Evidence for multiple polymorphisms of acetylcholinesterase gene with potential relevance for resistance to chemicals

The silverleaf whitefly Bemisia tabaci (Gennadius, 1889) (Homoptera: Aleyrodidae) is a serious invasive herbivorous insect pest worldwide. The excessive use of pesticides has progressively selected B. tabaci specimens, reducing the effectiveness of the treatments, and ultimately ending in the selection of pesticide-resistant strains. The management of this crop pest has thus become challenging owing to the level of resistance to all major classes of recommended insecticides. Here, we used in silico techniques for detecting sequence polymorphisms in ace1 gene from naturally occurring B. tabaci variants, and monitor the presence and frequency of the detected putative mutations from 30 populations of the silverleaf whitefly from Egypt and Pakistan. We found several point mutations in ace1-type acetylcholinesterase (ace1) in the studied B. tabaci variants naturally occurring in the field. By comparing ace1 sequence data from an organophosphate-susceptible and an organophosphate-resistant strains of B. tabaci to ace1 sequence data retrieved from GenBank for that species and to nucleotide polymorphisms from other arthropods, we identified novel mutations that could potentially influence insecticide resistance. Homology modeling and molecular docking analyses were performed to determine if the mutation-induced changes in form 1 acetylcholinesterase (AChE1) structure could confer resistance to carbamate and organophosphate insecticides. Mutations had small effects on binding energy (ΔG_b) interactions between mutant AChE1 and insecticides; they altered the conformation of the peripheral anionic site of AChE1, and modified the enzyme surface, and these changes have potential effects on the target-site sensitivity. Altogether, the results from this study provide information on genic variants of B. tabaci ace1 for future monitoring insecticide resistance development and report a potential case of environmentally driven gene variations.

Mortality of Orius insidiosus by contact with spinosad in the laboratory as well as in the field and a perspective of these as controllers of Frankliniella occidentalis

Orius insidiosus, known as the pirate bug, is widely distributed throughout the Americas. It is employed for the biological control of Frankliniella occidentalis in organic berry crops in Mexico. In conventional crops, spinosad is the main control method for this pest. The LD₅₀ of spinosad on O. insidiosus was determined. In addition, we monitored the population density of F. occidentalis in blackberry crops under two types of management (biochemical+mass trapping, and biological control). The LD₅₀ was 225.65 ppm 3.8 times greater than the 60 ppm dose commonly used in blackberry crops. Both types of control are efficient; however, spinosad is less effective and should be combined with other environmentally friendly strategies. The possibility of combining chromatic traps+spinosad application and chromatic traps+strategic release of O. insidiosus to effectively control thrips without compromising fruit quality is discussed.

Distinct roles of two RDL GABA receptors in fipronil action in the diamondback moth (Plutella xylostella)

The phenylpyrazole insecticide fipronil blocks resistance to dieldrin (RDL) γ-aminobutyric acid (GABA) receptors in insects, thereby impairing inhibitory neurotransmission. Some insect species, such as the diamondback moth (Plutella xylostella), possess more than one Rdl gene. The involvement of multiple Rdls in fipronil toxicity and resistance remains largely unknown. In this study, we investigated the roles of two Rdl genes, PxRdl1 and PxRdl2, in P. xylostella fipronil action. In Xenopus oocytes, PxRDL2 receptors were 40 times less sensitive to fipronil than PxRDL1. PxRDL2 receptors were also less sensitive to GABA compared with PxRDL1. Knockout of the fipronil-sensitive PxRdl1 reduced the fipronil potency 10-fold, whereas knockout of the fipronil-resistant PxRdl2 enhanced the fipronil potency 4.4-fold. Furthermore, in two fipronil-resistant diamondback moth field populations, PxRdl2 expression was elevated 3.7- and 4.1-fold compared with a susceptible strain, whereas PxRdl1 expression was comparable among the resistant and susceptible strains. Collectively, our results indicate antagonistic effects of PxRDL1 and PxRDL2 on fipronil action in vivo and suggest that enhanced expression of fipronil-resistant PxRdl2 is potentially a new mechanism of fipronil resistance in insects.

Reduced neuronal sensitivity and susceptibility of the fall armyworm, Spodoptera frugiperda, to pyrethroids in the absence of known knockdown mutations

Neurophysiological recordings were employed to quantify neuronal sensitivity to neurotoxic insecticides and assessed toxicity across field and laboratory fall armyworm (FAW) populations. Topical toxicity resistance ratios (RR) in field-collected FAW was 767-fold compared to laboratory strains and, importantly, a 1750-fold reduction in potency was observed for λ-cyhalothrin in neurophysiological assays. Field collected FAW were found to have a RR of 12 to chlorpyrifos when compared to the susceptible strain and was 8-fold less sensitive in neurophysiological assays. Surprisingly, there were no point mutations identified in the voltage-gated sodium channel known to cause pyrethroid resistance. For acetylcholinesterase, FAW had more than 80% of their nucleotide sequences consistent with A201 and F290 of the susceptible strains although 60% of the tested population was heterozygous for the G227A mutation. These data indicate that point mutations did not contribute to the high level of pyrethroid resistance and nerve insensitivity in this population of field collected FAW. Additionally, these data suggest the kdr phenotype only explains a portion of the heritable variation in FAW resistance and indicates kdr is not the only predictor of high pyrethroid resistance. Phenotypic assays, such as toxicity bioassays or neurophysiological recordings, using field-collected populations are necessary to reliably predict resistant phenotypes and product failures.

Kdr genotyping in Aedes aegypti from Brazil on a nation-wide scale from 2017 to 2018

Insecticide resistance is currently a threat to the control of Aedes agypti, the main vector of arboviruses in urban centers. Mutations in the voltage gated sodium channel (NaV), known as kdr (knockdown resistance), constitute an important selection mechanism for resistance against pyrethroids. In the present study, we investigated the kdr distribution for the Val1016Ile and Phe1534Cys alterations in Ae. aegypti from 123 Brazilian municipalities, based on SNP genotyping assays in over 5,500 mosquitoes. The alleles NaVS (1016Val+ + 1534Phe+), NaVR1 (1016Val+ + 1534Cyskdr) and NaVR2 (1016Ilekdr + 1534Cyskdr) were consistently observed, whereas kdr alleles have rapidly spread and increased in frequency. NaVS was the less frequent allele, mostly found in Northeastern populations. The highest allelic frequencies were observed for NaVR1, especially in the North, which was fixed in one Amazonian population. The double kdr NaVR2 was more prevalent in the Central-west and South-eastern populations. We introduce the 'kdr index', which revealed significant spatial patterns highlighting two to three distinct Brazilian regions. The 410L kdr mutation was additionally evaluated in 25 localities, evidencing that it generally occurs in the NaVR2 allele. This nationwide screening of a genetic mechanism for insecticide resistance is an important indication on how pyrethroid resistance in Ae. aegypti is evolving in Brazil.

Tracing temporal and geographic distribution of resistance to pyrethroids in the arboviral vector Aedes albopictus

BACKGROUND

The arboviral vector Aedes albopictus became established on all continents except Antarctica in the past 50 years. A consequence of its rapid global invasion is the transmission of diseases previously confined to the tropics and subtropics occurring in temperate regions of the world, including the re-emergence of chikungunya and dengue in Europe. Application of pyrethroids is among the most widely-used interventions for vector control, especially in the presence of an arboviral outbreak. Studies are emerging that reveal phenotypic resistance and monitor mutations at the target site, the para sodium channel gene, primarily on a local scale.

METHODS

A total of 512 Ae. albopictus mosquitoes from twelve geographic sites, including those from the native home range and invaded areas, were sampled between 2011 and 2018, and were analyzed at five codons of the para sodium channel gene with mutations predictive of resistance phenotype. Additionally, to test for the origin of unique kdr mutations in Mexico, we analyzed the genetic connectivity of southern Mexico mosquitoes with mosquitoes from home range, the Reunion Island, America and Europe.

RESULTS

We detected mutations at all tested positions of the para sodium channel gene, with heterozygotes predominating and rare instance of double mutants. We observed an increase in the distribution and frequency of F1534C/L/S mutations in the ancestral China population and populations in the Mediterranean Greece, the appearance of the V1016G/I mutations as early as 2011 in Italy and mutations at position 410 and 989 in Mexico. The analyses of the distribution pattern of kdr alleles and haplotype network analyses showed evidence for multiple origins of all kdr mutations.

CONCLUSIONS

Here we provide the most-up-to-date survey on the geographic and temporal distribution of pyrethroid-predictive mutations in Ae. albopictus by combining kdr genotyping on current and historical samples with published data. While we confirm low levels of pyrethroid resistance in most analyzed samples, we find increasing frequencies of F1534C/S and V1016G in China and Greece or Italy, respectively. The observed patterns of kdr allele distribution support the hypothesis that on site emergence of resistance has contributed more than spread of resistance through mosquito migration/invasions to the current widespread of kdr alleles, emphasizing the importance of local surveillance programs and resistance management.

Evolution of kdr haplotypes in worldwide populations of Aedes aegypti: Independent origins of the F1534C kdr mutation

Aedes aegypti is the primary vector of dengue, chikungunya, Zika, and urban yellow fever. Insecticides are often the most effective tools to rapidly decrease the density of vector populations, especially during arbovirus disease outbreaks. However, the intense use of insecticides, particularly pyrethroids, has selected for resistant mosquito populations worldwide. Mutations in the voltage gated sodium channel (Na_V) are among the principal mechanisms of resistance to pyrethroids and DDT, also known as “knockdown resistance,” kdr. Here we report studies on the origin and dispersion of kdr haplotypes in samples of Ae. aegypti from its worldwide distribution. We amplified the IIS6 and IIIS6 Na_V segments from pools of Ae. aegypti populations from 15 countries, in South and North America, Africa, Asia, Pacific, and Australia. The amplicons were barcoded and sequenced using NGS Ion Torrent. Output data were filtered and analyzed using the bioinformatic pipeline Seekdeep to determine frequencies of the IIS6 and IIIS6 haplotypes per population. Phylogenetic relationships among the haplotypes were used to infer whether the kdr mutations have a single or multiple origin. We found 26 and 18 haplotypes, respectively for the IIS6 and IIIS6 segments, among which were the known kdr mutations 989P, 1011M, 1016I and 1016G (IIS6), 1520I, and 1534C (IIIS6). The highest diversity of haplotypes was found in African samples. Kdr mutations 1011M and 1016I were found only in American and African populations, 989P + 1016G and 1520I + 1534C in Asia, while 1534C was present in samples from all continents, except Australia. Based primarily on the intron sequence, IIS6 haplotypes were subdivided into two well-defined clades (A and B). Subsequent phasing of the IIS6 + IIIS6 haplotypes indicates two distinct origins for the 1534C kdr mutation. These results provide evidence of kdr mutations arising de novo at specific locations within the Ae. aegypti geographic distribution. In addition, our results suggest that the 1534C kdr mutation had at least two independent origins. We can thus conclude that insecticide selection pressure with DDT and more recently with pyrethroids is selecting for independent convergent mutations in Na_V.

Insecticide resistance is a global threat for the control of Aedes aegypti, the mosquito vector of aboviruses such as dengue, chikungunya and Zika. Mutations in the voltage gated sodium channel (Na_V), known as kdr, are one of the principal mechanisms related to resistance to pyrethroids, the class of insecticide most employed worldwide inside and around residences. We investigate whether the same kdr mutations found in Ae. aegypti populations from distinct regions of the world have a common origin and subsequently dispersed or if they emerged in unrelated populations at distinct moments. By evaluating the sequences of two fragments of the Na_V gene, obtained from DNA collections of Ae. aegypti from several countries, we found at least two independent origins for the F1534C kdr mutation in American, African and Asian populations. There was no evidence for multiple origins of the common kdr mutations V1016I and P989S + V1016G, which were exclusive to American and Asian populations. Our results increase our knowledge of insecticide resistance evolution in one of the main arboviral mosquito vectors of major global diseases.

Properties, toxicity and current applications of the biolarvicide spinosad

Characterized as a highly valuable bioactive natural product, spinosad is a pesticide with a complex chemical structure, composed of spinosyn A and D, molecules synthesized by the actinomycete Saccharopolyspora spinosa. The larvicidal activity of spinosad was postulated to be a promising approach to combat crop pests and control species responsible to transmit mosquito-borne illness, including Aedes aegypti. Although initially deemed as relatively safe for non-target organisms and highly effective against insects and crop pests, recent studies focused on the toxicity profile detected the occurrence of side effects in different living species. Thus, the present review was undertaken to describe the properties and characteristics of spinosad. In addition to indicating potential adverse effects on living organisms, alternative uses of the biopesticide as a mixture with different compounds are provided.

New ways to acquire resistance: imperfect convergence in insect adaptations to a potent plant toxin

Evolution of insensitivity to the toxic effects of cardiac glycosides has become a model in the study of convergent evolution, as five taxonomic orders of insects use the same few similar amino acid substitutions in the otherwise highly conserved Na,K-ATPase α. We show here that insensitivity in pyrgomorphid grasshoppers evolved along a slightly divergent path. As in other lineages, duplication of the Na,K-ATPase α gene paved the way for subfunctionalization: one copy maintains the ancestral, sensitive state, while the other copy is resistant. Nonetheless, in contrast with all other investigated insects, the grasshoppers' resistant copy shows length variation by two amino acids in the first extracellular loop, the main part of the cardiac glycoside-binding pocket. RT-qPCR analyses confirmed that this copy is predominantly expressed in tissues exposed to the toxins, while the ancestral copy predominates in the nervous tissue. Functional tests with genetically engineered Drosophila Na,K-ATPases bearing the first extracellular loop of the pyrgomorphid genes showed the derived form to be highly resistant, while the ancestral state is sensitive. Thus, we report convergence in gene duplication and in the gene targets for toxin insensitivity; however, the means to the phenotypic end have been novel in pyrgomorphid grasshoppers.

Readapting to DCV Infection without Wolbachia: Frequency Changes of Drosophila Antiviral Alleles Can Replace Endosymbiont Protection

There is now ample evidence that endosymbionts can contribute to host adaptation to environmental challenges. However, how endosymbiont presence affects the adaptive trajectory and outcome of the host is yet largely unexplored. In Drosophila, Wolbachia confers protection to RNA virus infection, an effect that differs between Wolbachia strains and can be targeted by selection. Adaptation to RNA virus infections is mediated by both Wolbachia and the host, raising the question of whether adaptive genetic changes in the host vary with the presence/absence of the endosymbiont. Here, we address this question using a polymorphic D. melanogaster population previously adapted to DCV infection for 35 generations in the presence of Wolbachia, from which we removed the endosymbiont and followed survival over the subsequent 20 generations of infection. After an initial severe drop, survival frequencies upon DCV selection increased significantly, as seen before in the presence of Wolbachia. Whole-genome sequencing, revealed that the major genes involved in the first selection experiment, pastrel and Ubc-E2H, continued to be selected in Wolbachia-free D. melanogaster, with the frequencies of protective alleles being closer to fixation in the absence of Wolbachia. Our results suggest that heterogeneity in Wolbachia infection status may be sufficient to maintain polymorphisms even in the absence of costs.

Parallel genetic adaptation across environments differing in mode of growth or resource availability

Evolution experiments have demonstrated high levels of genetic parallelism between populations evolving in identical environments. However, natural populations evolve in complex environments that can vary in many ways, likely sharing some characteristics but not others. Here, we ask whether shared selection pressures drive parallel evolution across distinct environments. We addressed this question in experimentally evolved populations founded from a clone of the bacterium Burkholderia cenocepacia. These populations evolved for 90 days (approximately 600 generations) under all combinations of high or low carbon availability and selection for either planktonic or biofilm modes of growth. Populations that evolved in environments with shared selection pressures (either level of carbon availability or mode of growth) were more genetically similar to each other than populations from environments that shared neither characteristic. However, not all shared selection pressures led to parallel evolution. Genetic parallelism between low-carbon biofilm and low-carbon planktonic populations was very low despite shared selection for growth under low-carbon conditions, suggesting that evolution in low-carbon environments may generate stronger trade-offs between biofilm and planktonic modes of growth. For all environments, a population's fitness in a particular environment was positively correlated with the genetic similarity between that population and the populations that evolved in that particular environment. Although genetic similarity was low between low-carbon environments, overall, evolution in similar environments led to higher levels of genetic parallelism and that genetic parallelism, in turn, was correlated with fitness in a particular environment.

Levels of Resistance to Pyrethroid among Distinct kdr Alleles in Aedes aegypti Laboratory Lines and Frequency of kdr Alleles in 27 Natural Populations from Rio de Janeiro, Brazil

Background

Several mutations in voltage gated sodium channel (Na_V) have been identified in Aedes aegypti populations worldwide. However, only few are related to knockdown resistance to pyrethroids, most of which with variations in the 1016 and 1534 Na_V sites. In Brazil, at least two Na_V alleles are known: Na_VR1, with a substitution in the 1534 (1016 Val⁺ + 1534 Ile^kdr) and Na_VR2, with substitutions in both 1016 and sites (1016Ile^kdr + 1534Cys^kdr). There is also the duplication in the Na_V gene, with one copy carrying the substitution Ile1011Met, although its effects on pyrethroid resistance remain to be clarified. Our goals in this study were (1) to determine the role of each kdr Na_V allele and the duplication on pyrethroid resistance and (2) to screen the frequency of the kdr alleles in 27 several natural Ae. aegypti populations from the metropolitan region of Rio de Janeiro.

Methods

Pyrethroid resistance was evaluated by a knockdown time (KdT) assay, an adaptation of the WHO test tubes with paper impregnated with deltamethrin. We used laboratory-selected Ae. aegypti lineages: R1R1 and R2R2 (homozygous for the kdr Na_VR1 and Na_VR2 alleles, respectively), Dup (with duplication in the Na_V gene), Rockefeller (the susceptibility reference control), and F1 hybrids among them. Genotyping of both 1016 and 1534 Na_V sites was performed in 811 Ae. aegypti sampled from 27 localities from Rio de Janeiro (17), Niterói (6) and Nova Iguaçu (4) cities, Rio de Janeiro State, Brazil, with a TaqMan real time PCR approach.

Results

The laboratory lineages R1R1, R2R2, and R1R2 were the only ones that needed more than 60 minutes to knock down all the insects exposed to the pyrethroid, being the KdT R2R2 > R1R2 > R1R1, corroborating the recessive nature of the kdr mutations. Frequency of kdr alleles Na_VR1 and Na_VR2 in field-caught mosquitoes varied from 0 to 52% and 43 to 86%, respectively, evidencing high levels of “resistant genotypes” (R1R1, R1R2, and R2R2), which together summed 60 to 100% in Ae. aegypti populations from Rio de Janeiro.

Conclusions

The Na_VR1 and Na_VR2 kdr alleles confer resistance to the pyrethroid deltamethrin in homozygotes and R1R2 heterozygotes, being the R2R2 most resistant genotype. The allele containing duplication in the Na_V gene, with a mutation in the 1011 site, did not confer resistance under the tested conditions. The frequencies of the “resistant genotypes” are elevated in Ae. aegypti natural populations from Rio de Janeiro.

Toward a Predictive Framework for Convergent Evolution: Integrating Natural History, Genetic Mechanisms, and Consequences for the Diversity of Life

A charm of biology as a scientific discipline is the diversity of life. Although this diversity can make laws of biology challenging to discover, several repeated patterns and general principles govern evolutionary diversification. Convergent evolution, the independent evolution of similar phenotypes, has been at the heart of one approach to understand generality in the evolutionary process. Yet understanding when and why organismal traits and strategies repeatedly evolve has been a central challenge. These issues were the focus of the American Society of Naturalists Vice Presidential Symposium in 2016 and are the subject of this collection of articles. Although naturalists have long made inferences about convergent evolution and its importance, there has been confusion in the interpretation of the pattern of convergence. Does convergence primarily indicate adaptation or constraint? How often should convergence be expected? Are there general principles that would allow us to predict where and when and by what mechanisms convergent evolution should occur? What role does natural history play in advancing our understanding of general evolutionary principles? In this introductory article, I address these questions, review several generalizations about convergent evolution that have emerged over the past 15 years, and present a framework for advancing the study and interpretation of convergence. Perhaps the most important emerging conclusion is that the genetic mechanisms of convergent evolution are phylogenetically conserved; that is, more closely related species tend to share the same genetic basis of traits, even when independently evolved. Finally, I highlight how the articles in this special issue further develop concepts, methodologies, and case studies at the frontier of our understanding of the causes and consequences of convergent evolution.

Synergistic and compensatory effects of two point mutations conferring target-site resistance to fipronil in the insect GABA receptor RDL

Insecticide resistance can arise from a variety of mechanisms, including changes to the target site, but is often associated with substantial fitness costs to insects. Here we describe two resistance-associated target-site mutations that have synergistic and compensatory effects that combine to produce high and persistent levels of resistance to fipronil, an insecticide targeting on γ-aminobytyric acid (GABA) receptors. In Nilaparvata lugens, a major pest of rice crops in many parts of Asia, we have identified a single point mutation (A302S) in the GABA receptor RDL that has been identified previously in other species and which confers low levels of resistance to fipronil (23-fold) in N. lugans. In addition, we have identified a second resistance-associated RDL mutation (R300Q) that, in combination with A302S, is associated with much higher levels of resistance (237-fold). The R300Q mutation has not been detected in the absence of A302S in either laboratory-selected or field populations, presumably due to the high fitness cost associated with this mutation. Significantly, it appears that the A302S mutation is able to compensate for deleterious effects of R300Q mutation on fitness cost. These findings identify a novel resistance mechanism and may have important implications for the spread of insecticide resistance.

A CRISPR/Cas9 mediated point mutation in the alpha 6 subunit of the nicotinic acetylcholine receptor confers resistance to spinosad in Drosophila melanogaster

Spinosad, a widely used and economically important insecticide, targets the nicotinic acetylcholine receptor (nAChRs) of the insect nervous system. Several studies have associated loss of function mutations in the insect nAChR α6 subunit with resistance to spinosad, and in the process identified this particular subunit as the specific target site. More recently a single non-synonymous point mutation, that does not result in loss of function, was identified in spinosad resistant strains of three insect species that results in an amino acid substitution (G275E) of the nAChR α6 subunit. The causal role of this mutation has been called into question as, to date, functional evidence proving its involvement in resistance has been limited to the study of vertebrate receptors. Here we use the CRISPR/Cas9 gene editing platform to introduce the G275E mutation into the nAChR α6 subunit of Drosophila melanogaster. Reverse transcriptase-PCR and sequencing confirmed the presence of the mutation in Dα6 transcripts of mutant flies and verified that it does not disrupt the normal splicing of the two exons in close vicinity to the mutation site. A marked decrease in sensitivity to spinosad (66-fold) was observed in flies with the mutation compared to flies of the same genetic background minus the mutation, clearly demonstrating the functional role of this amino acid substitution in resistance to spinosad. Although the resistance levels observed are 4.7-fold lower than exhibited by a fly strain with a null mutation of Dα6, they are nevertheless predicated to be sufficient to result in resistance to spinosad at recommended field rates. Reciprocal crossings with susceptible fly strains followed by spinosad bioassays revealed G275E is inherited as an incompletely recessive trait, thus resembling the mode of inheritance described for this mutation in the western flower thrips, Frankliniella occidentalis. This study both resolves a debate on the functional significance of a target-site mutation and provides an example of how recent advances in genome editing can be harnessed to study insecticide resistance.

Causes of molecular convergence and parallelism in protein evolution

To what extent is the convergent evolution of protein function attributable to convergent or parallel changes at the amino acid level? The mutations that contribute to adaptive protein evolution may represent a biased subset of all possible beneficial mutations owing to mutation bias and/or variation in the magnitude of deleterious pleiotropy. A key finding is that the fitness effects of amino acid mutations are often conditional on genetic background. This context dependence (epistasis) can reduce the probability of convergence and parallelism because it reduces the number of possible mutations that are unconditionally acceptable in divergent genetic backgrounds. Here, I review factors that influence the probability of replicated evolution at the molecular level.

The Molecular Evolution of Xenobiotic Metabolism and Resistance in Chelicerate Mites

Chelicerate mites diverged from other arthropod lineages more than 400 million years ago and subsequently developed specific and remarkable xenobiotic adaptations. The study of the two-spotted spider mite, Tetranychus urticae, for which a high-quality Sanger-sequenced genome was first available, revealed expansions and radiations in all major detoxification gene families, including P450 monooxygenases, carboxyl/cholinesterases, glutathione-S-transferases, and ATP-binding cassette transporters. Novel gene families that are not well studied in other arthropods, such as major facilitator family transporters and lipocalins, also reflect the evolution of xenobiotic adaptation. The acquisition of genes by horizontal gene transfer provided new routes to handle toxins, for example, the β-cyanoalanine synthase enzyme that metabolizes cyanide. The availability of genomic resources for other mite species has allowed researchers to study the lineage specificity of these gene family expansions and the distinct evolution of genes involved in xenobiotic metabolism in mites. Genome-based tools have been crucial in supporting the idiosyncrasies of mite detoxification and will further support the expanding field of mite-plant interactions.

Identification of a point mutation in the ace1 gene of Therioaphis trifolli maculata and detection of insecticide resistance by a diagnostic PCR-RFLP assay

Aphids are important agricultural pests worldwide. Their control is largely based on chemical insecticides. One species that shows important invasive abilities and host-plant-related differences is Therioaphis trifolii (Monell) (Hemiptera: Aphididae). T. trifolii maculata, also known as spotted alfalfa aphid (SAA), can be very injurious to alfalfa crops in certain regions, such as in Saudi Arabia for effective control it is essential to diagnose and monitor the resistance mechanisms in the SAA populations. In the present study, we analysed acetylcholinesterase (ace) target site insensitivity mechanisms. A 650 bp length DNA containing the putative acetylcholinesterase (ace1) precursor was obtained and compared with other Hemipteran species. The sequences of many individual aphids collected from alfalfa crops in Saudi Arabia were analysed for the presence of resistance mutations: no resistance mutations were found at the resistance mutation loci 302; however, the presence of a serine-phenylalanine substitution (S431F) was identified in one individual. The S431F substitution, has been shown to confer significant levels of both organophosphate and carbamate resistance in other aphid species, and is now found for the first time in T. trifolii. We subsequently developed a simple polymerase chain reaction-restriction fragment length polymorphism assays for the S431F mutation, using a TaqI restriction site destroyed by the S431F mutation. The novel diagnostic assay may support the implementation of Insecticide Resistance Management strategies, for the control of SAA in alfalfa crops in the Kingdom of Saudi Arabia, and other countries worldwide.

Geographical distribution and frequencies of organophosphate-resistant Ace alleles and morphometric variations in olive fruit fly populations

BACKGROUND

In the Mediterranean basin, organophosphate (OP) insecticides have been used intensively to control olive fly populations. Acetylcholinesterase (Ace) is the molecular target of OP insecticides, and three resistance-associated mutations that confer different levels of OP insensitivity have been identified. In this study, genotypes of olive fly Ace were determined in field-collected populations from broad geographical areas in Turkey. In addition, the levels of asymmetry of wing and leg characters were compared in these populations.

RESULTS

Our study revealed the existence of a genetically smooth stratification pattern in OP resistance allele distribution in the olive fly populations of Turkey. In contrast to earlier findings, the frequency of Δ3Q was found to be lower in the Aegean region, where the populations have been subjected to high selection pressure. Results based on the morphological differences among the samples revealed a similar pattern for both sides and did not demonstrate a clear separation.

CONCLUSION

The frequencies and geographic range of resistance alleles indicate that they were selected in the Aegean coast of Turkey and then spread westward towards Europe. One possible explanation for the absence of morphological asymmetry in olive fly samples might be the presence of modifier allele(s) that compensate for the increase in asymmetry.

Involvement of metabolic resistance and F1534C kdr mutation in the pyrethroid resistance mechanisms of Aedes aegypti in India

Pesticide resistance poses a serious problem for worldwide mosquito control programs. Resistance to insecticides can be caused by an increased metabolic detoxification of the insecticide and/or by target site insensitivity. In the present study, we estimated the tolerance of Indian Aedes aegypti populations using adult bioassays that revealed high resistance levels of the field populations to permethrin (RR-6, 5.8 and 5.1 folds) compared to our susceptible population. Enzymatic assays revealed increased activities of glutathione S-transferase and carboxylesterase enzymes in the field populations comparatively to the susceptible population. PBO synergist assays did not confirm that cytochrome P450 monooxygenase metabolic detoxification acted as a major cause of resistance. Hence the role of target site resistance was therefore investigated. A single substitution Phe1534Cys in the voltage gated sodium channel was found in domain III, segment 6 (III-S6) of the resistance populations (allele frequency=0.59, 0.51 and 0.47) suggesting its potential role in permethrin resistance in A. aegypti.

Cloning of Two Acetylcholinesterase Genes and Analysis of Point Mutations Putatively Associated with Triazophos Resistance in Chilo auricilius (Lepidoptera: Pyralidae)

Acetylcholinesterase (AChE) is the target of organophosphate (OP) and carbamate insecticides. Mutations in the AChE gene (ace) leading to decreased insecticide susceptibility is the main resistance mechanism in insects. In this study, two Chilo auricilius acetylcholinesterase genes, designated as Caace1 and Caace2, were cloned using RT-PCR and RACE. Caace1 cDNA is 2534 bp, with ORF of 2082 bp, and it encodes an acetylcholinesterase 1 (CaAChE1) protein comprising a calculated 693 amino acid (aa) residues. Caace2 cDNA contains 2280 bp, with a full-length ORF of 1917 bp, encoding acetylcholinesterase 2 (CaAChE2) comprising a calculated 638 aa residues. At the aa level, CaAChE1 displays the highest similarity (97%) with the Chilo suppressalis AChE1, and CaAChE2 shows the highest similarity with the C. suppressalis AChE2 (99%). From the restriction fragment length polymorphism (RFLP) PCR (RFLP-PCR) analysis, one mutation in Caace1, similar to the ace1 mutation associated with triazophos resistance in C. suppressalis, was detected. Detailed examination of field populations of C. auricilius indicated this resistance mutation in C. auricilius is still quite infrequent. Based on the assay of AChE activity and RFLP-PCR testing, an individual that contains resistance mutation has lower AChE activities, while the individual that does not contain the resistance mutation has higher AChE activities. This study provides a basis for future investigations into the mechanism of OP resistance in C. auricilius, as well as a guidance for C. auricilius control with reasonable choice of pesticides.

Mutations in the voltage-gated sodium channel gene of anophelines and their association with resistance to pyrethroids - a review

Constant and extensive use of chemical insecticides has created a selection pressure and favored resistance development in many insect species worldwide. One of the most important pyrethroid resistance mechanisms is classified as target site insensitivity, due to conformational changes in the target site that impair a proper binding of the insecticide molecule. The voltage-gated sodium channel (NaV) is the target of pyrethroids and DDT insecticides, used to control insects of medical, agricultural and veterinary importance, such as anophelines. It has been reported that the presence of a few non-silent point mutations in the NaV gene are associated with pyrethroid resistance, termed as 'kdr' (knockdown resistance) for preventing the knockdown effect of these insecticides. The presence of these mutations, as well as their effects, has been thoroughly studied in Anopheles mosquitoes. So far, kdr mutations have already been detected in at least 13 species (Anopheles gambiae, Anopheles arabiensis, Anopheles sinensis, Anopheles stephensi, Anopheles subpictus, Anopheles sacharovi, Anopheles culicifacies, Anopheles sundaicus, Anopheles aconitus, Anopheles vagus, Anopheles paraliae, Anopheles peditaeniatus and Anopheles albimanus) from populations of African, Asian and, more recently, American continents. Seven mutational variants (L1014F, L1014S, L1014C, L1014W, N1013S, N1575Y and V1010L) were described, with the highest prevalence of L1014F, which occurs at the 1014 site in NaV IIS6 domain. The increase of frequency and distribution of kdr mutations clearly shows the importance of this mechanism in the process of pyrethroid resistance. In this sense, several species-specific and highly sensitive methods have been designed in order to genotype individual mosquitoes for kdr in large scale, which may serve as important tolls for monitoring the dynamics of pyrethroid resistance in natural populations. We also briefly discuss investigations concerning the course of Plasmodium infection in kdr individuals. Considering the limitation of insecticides available for employment in public health campaigns and the absence of a vaccine able to brake the life cycle of the malaria parasites, the use of pyrethroids is likely to remain as the main strategy against mosquitoes by either indoor residual spraying (IR) and insecticide treated nets (ITN). Therefore, monitoring insecticide resistance programs is a crucial need in malaria endemic countries.

Parallel evolution of tetrodotoxin resistance in three voltage-gated sodium channel genes in the garter snake Thamnophis sirtalis

Members of a gene family expressed in a single species often experience common selection pressures. Consequently, the molecular basis of complex adaptations may be expected to involve parallel evolutionary changes in multiple paralogs. Here, we use bacterial artificial chromosome library scans to investigate the evolution of the voltage-gated sodium channel (Na_v) family in the garter snake Thamnophis sirtalis, a predator of highly toxic Taricha newts. Newts possess tetrodotoxin (TTX), which blocks Na_v’s, arresting action potentials in nerves and muscle. Some Thamnophis populations have evolved resistance to extremely high levels of TTX. Previous work has identified amino acid sites in the skeletal muscle sodium channel Na_v1.4 that confer resistance to TTX and vary across populations. We identify parallel evolution of TTX resistance in two additional Na_v paralogs, Na_v1.6 and 1.7, which are known to be expressed in the peripheral nervous system and should thus be exposed to ingested TTX. Each paralog contains at least one TTX-resistant substitution identical to a substitution previously identified in Na_v1.4. These sites are fixed across populations, suggesting that the resistant peripheral nerves antedate resistant muscle. In contrast, three sodium channels expressed solely in the central nervous system (Na_v1.1–1.3) showed no evidence of TTX resistance, consistent with protection from toxins by the blood–brain barrier. We also report the exon–intron structure of six Na_v paralogs, the first such analysis for snake genes. Our results demonstrate that the molecular basis of adaptation may be both repeatable across members of a gene family and predictable based on functional considerations.

The evolution of insecticide resistance in the peach potato aphid, Myzus persicae

The peach potato aphid, Myzus persicae is a globally distributed crop pest with a host range of over 400 species including many economically important crop plants. The intensive use of insecticides to control this species over many years has led to populations that are now resistant to several classes of insecticide. Work spanning over 40 years has shown that M. persicae has a remarkable ability to evolve mechanisms that avoid or overcome the toxic effect of insecticides with at least seven independent mechanisms of resistance described in this species to date. The array of novel resistance mechanisms, including several 'first examples', that have evolved in this species represents an important case study for the evolution of insecticide resistance and also rapid adaptive change in insects more generally. In this review we summarise the biochemical and molecular mechanisms underlying resistance in M. persicae and the insights study of this topic has provided on how resistance evolves, the selectivity of insecticides, and the link between resistance and host plant adaptation.

A de novo transcriptome of European pollen beetle populations and its analysis, with special reference to insecticide action and resistance

The pollen beetle Meligethes aeneus is the most important coleopteran pest in European oilseed rape cultivation, annually infesting millions of hectares and responsible for substantial yield losses if not kept under economic damage thresholds. This species is primarily controlled with insecticides but has recently developed high levels of resistance to the pyrethroid class. The aim of the present study was to provide a transcriptomic resource to investigate mechanisms of resistance. cDNA was sequenced on both Roche (Indianapolis, IN, USA) and Illumina (LGC Genomics, Berlin, Germany) platforms, resulting in a total of ∼53 m reads which assembled into 43 396 expressed sequence tags (ESTs). Manual annotation revealed good coverage of genes encoding insecticide target sites and detoxification enzymes. A total of 77 nonredundant cytochrome P450 genes were identified. Mapping of Illumina RNAseq sequences (from susceptible and pyrethroid-resistant strains) against the reference transcriptome identified a cytochrome P450 (CYP6BQ23) as highly overexpressed in pyrethroid resistance strains. Single-nucleotide polymorphism analysis confirmed the presence of a target-site resistance mutation (L1014F) in the voltage-gated sodium channel of one resistant strain. Our results provide new insights into the important genes associated with pyrethroid resistance in M. aeneus. Furthermore, a comprehensive EST resource is provided for future studies on insecticide modes of action and resistance mechanisms in pollen beetle.

Molecular population genetics of the NADPH cytochrome P450 reductase (CPR) gene in Anopheles minimus

Development of insecticide resistance (IR) in mosquito vectors is a primary huddle to malaria control program. Since IR has genetic basis, and genes constantly evolve with response to environment for adaptation to organisms, it is important to know evolutionary pattern of genes conferring IR in malaria vectors. The mosquito Anopheles minimus is a major malaria vector of the Southeast (SE) Asia and India and is susceptible to all insecticides, and thus of interest to know if natural selection has shaped variations in the gene conferring IR. If not, the DNA fragment of such a gene could be used to infer population structure and demography of this species of malaria vector. We have therefore sequenced a ~569 bp DNA segment of the NADPH cytochrome P450 reductase (CPR) gene (widely known to confer IR) in 123 individuals of An. minimus collected in 10 different locations (eight Indian, one Thai and one Vietnamese). Two Indian population samples were completely mono-morphic in the CPR gene. In general, low genetic diversity was found with no evidence of natural selection in this gene. The data were therefore analyzed to infer population structure and demography of this species. The 10 populations could be genetically differentiated into four different groups; the samples from Thailand and Vietnam contained high nucleotide diversity. All the 10 populations conform to demographic equilibrium model with signature of past population expansion in four populations. The results in general indicate that the An. minimus mosquitoes sampled in the two SE Asian localities contain several genetic characteristics of being parts of the ancestral population.

Resistance is not Futile: It Shapes Insecticide Discovery

Conventional chemical control compounds used for the management of insect pests have been much maligned, but still serve a critical role in protecting people and agricultural products from insect pests, as well as conserving biodiversity by eradicating invasive species. Although biological control can be an effective option for area-wide management of established pests, chemical control methods are important for use in integrated pest management (IPM) programs, as well as in export treatments, eradicating recently arrived invasive species, and minimizing population explosions of vectors of human disease. Cogitated research and development programs have continued the innovation of insecticides, with a particular focus on combating insecticide resistance. Recent developments in the fields of human health, protecting the global food supply, and biosecurity will be highlighted.

Distribution and dissemination of the Val1016Ile and Phe1534Cys Kdr mutations in Aedes aegypti Brazilian natural populations

BACKGROUND

The chemical control of the mosquito Aedes aegypti, the major vector of dengue, is being seriously threatened due to the development of pyrethroid resistance. Substitutions in the 1016 and 1534 sites of the voltage gated sodium channel (AaNaV), commonly known as kdr mutations, confer the mosquito with knockdown resistance. Our aim was to evaluate the allelic composition of natural populations of Brazilian Ae. aegypti at both kdr sites.

METHODS

The AaNaV IIIS6 region was cloned and sequenced from three Brazilian populations. Additionally, individual mosquitoes from 30 populations throughout the country were genotyped for 1016 and 1534 sites, based in allele-specific PCR. For individual genotypes both sites were considered as a single locus.

RESULTS

The 350 bp sequence spanning the IIIS6 region of the AaNaV gene revealed the occurrence of the kdr mutation Phe1534Cys in Brazil. Concerning the individual genotyping, beyond the susceptible wild-type (NaVS), two kdr alleles were identified: substitutions restricted to the 1534 position (NaVR1) or simultaneous substitutions in both 1016 and 1534 sites (NaVR2). A clear regional distribution pattern of these alleles was observed. The NaVR1kdr allele occurred in all localities, while NaVR2 was more frequent in the Central and Southeastern localities. Locations that were sampled multiple times in the course of a decade revealed an increase in frequency of the kdr mutations, mainly the double mutant allele NaVR2. Recent samples also indicate that NaVR2 is spreading towards the Northern region.

CONCLUSIONS

We have found that in addition to the previously reported Val1016Ile kdr mutation, the Phe1534Cys mutation also occurs in Brazil. Allelic composition at both sites was important to elucidate the actual distribution of kdr mutations throughout the country. Studies to determine gene flow and the fitness costs of these kdr alleles are underway and will be important to better understand the dynamics of Ae. aegypti pyrethroid resistance.

Multiple insecticide resistances in the disease vector Culex p. quinquefasciatus from Western Indian Ocean

Several mosquito-borne diseases affect the Western Indian Ocean islands. Culex pipiens quinquefasciatus is one of these vectors and transmits filariasis, Rift Valley and West Nile viruses and the Japanese encephalitis. To limit the impact of these diseases on public health, considerable vector control efforts have been implemented since the 50s, mainly through the use of neurotoxic insecticides belonging to Organochlorines (OC), Organophosphates (OP) and pyrethroids (PYR) families. However, mosquito control failures have been reported on site, and they were probably due to the selection of resistant individuals in response to insecticide exposure. In this study, we used different approaches to establish a first regional assessment of the levels and mechanisms of resistance to various insecticides. Bioassays were used to evaluate resistance to various insecticides, enzyme activity was measured to assess the presence of metabolic resistances through elevated detoxification, and molecular identification of known resistance alleles was investigated to determine the frequency of target-site mutations. These complementary approaches showed that resistance to the most used insecticides families (OC, OP and PYR) is widespread at a regional scale. However, the distribution of the different resistance genes is quite heterogeneous among the islands, some being found at high frequencies everywhere, others being frequent in some islands and absent in others. Moreover, two resistance alleles displayed clinal distributions in Mayotte and La Réunion, probably as a result of a heterogeneous selection due to local treatment practices. These widespread and diverse resistance mechanisms reduce the capacity of resistance management through classical strategies (e.g. insecticide rotation). In case of a disease outbreak, it could undermine the efforts of the vector control services, as only few compounds could be used. It thus becomes urgent to find alternatives to control populations of Cx. p. quinquefasciatus in the Indian Ocean.

Gene duplication in the major insecticide target site, Rdl, in Drosophila melanogaster

The Resistance to Dieldrin gene, Rdl, encodes a GABA-gated chloride channel subunit that is targeted by cyclodiene and phenylpyrazole insecticides. The gene was first characterized in Drosophila melanogaster by genetic mapping of resistance to the cyclodiene dieldrin. The 4,000-fold resistance observed was due to a single amino acid replacement, Ala(301) to Ser. The equivalent change was subsequently identified in Rdl orthologs of a large range of resistant insect species. Here, we report identification of a duplication at the Rdl locus in D. melanogaster. The 113-kb duplication contains one WT copy of Rdl and a second copy with two point mutations: an Ala(301) to Ser resistance mutation and Met(360) to Ile replacement. Individuals with this duplication exhibit intermediate dieldrin resistance compared with single copy Ser(301) homozygotes, reduced temperature sensitivity, and altered RNA editing associated with the resistant allele. Ectopic recombination between Roo transposable elements is involved in generating this genomic rearrangement. The duplication phenotypes were confirmed by construction of a transgenic, artificial duplication integrating the 55.7-kb Rdl locus with a Ser(301) change into an Ala(301) background. Gene duplications can contribute significantly to the evolution of insecticide resistance, most commonly by increasing the amount of gene product produced. Here however, duplication of the Rdl target site creates permanent heterozygosity, providing unique potential for adaptive mutations to accrue in one copy, without abolishing the endogenous role of an essential gene.

Gene flow across a hybrid zone maintained by a weak heterogametic incompatibility and positive selection of incompatible alleles

Hybridization between incipient species is more likely to produce sterile or inviable F(1) offspring in the heterogametic (XY or ZW) sex than in the homogametic (XX or ZZ) sex, a phenomenon known as Haldane's rule. Population dynamics associated with Haldane's rule may play an important role in early speciation of sexually reproducing organisms. The dynamics of the hybrid zone maintained by incomplete hybrid inferiority (sterility/inviability) in the heterogametic sex (a 'weak' Haldane's rule) caused by a Bateson-Dobzhansky-Muller incompatibility was modelled. The influences and interplays of the strengths of incompatibility, dispersal, density-dependent regulation (DDR) and local adaptation of incompatible alleles in a scenario of short-range dispersal (the stepping-stone model) were examined. It was found that a partial heterogametic hybrid incompatibility could efficiently impede gene flow and maintain characteristic clinal noncoincidence and discordance of alleles. Density-dependent regulation appears to be an important factor affecting hybrid zone dynamics: it can effectively skew the effects of the partial incompatibility and dispersal as measured by effective dispersal, clinal structures and density depression. Unexpectedly, local adaptation of incompatible alleles in the parental populations, which would be critical for the establishment of the incompatibility, exerts little effect on hybrid zone dynamics. These results strongly support the plausibility of the adaptive origin of hybrid incompatibility and ecological speciation: an adaptive mutation, if it confers a marginal fitness advantage in the local population and happens to cause epistatic inferiority in hybrids, could efficiently drive further genetic divergence that may result in the gene becoming an evolutionary hotspot.

Expression profile of genes during resistance reversal in a temephos selected strain of the dengue vector, Aedes aegypti

Background

The mosquito Aedes aegypti is one of the most important disease vectors because it transmits two major arboviruses, dengue and yellow fever, which cause significant global morbidity and mortality. Chemical insecticides form the cornerstone of vector control. The organophosphate temephos a larvicide recommended by WHO for controlling Ae. aegypti, however, resistance to this compound has been reported in many countries, including Brazil.

Methodology/Principal Findings

The aim of this study was to identify genes implicated in metabolic resistance in an Ae. aegypti temephos resistant strain, named RecR, through microarray analysis. We utilized a custom ‘Ae. aegypti detox chip’ and validated microarray data through RT-PCR comparing susceptible and resistant individuals. In addition, we analyzed gene expression in 4^th instar larvae from a reversed susceptible strain (RecRev), exposed and unexposed to temephos. The results obtained revealed a set of 13 and 6 genes significantly over expressed in resistant adult mosquitoes and larvae, respectively. One of these genes, the cytochrome P450 CYP6N12, was up-regulated in both stages. RT-PCR confirmed the microarray results and, additionally, showed no difference in gene expression between temephos exposed and unexposed RecRev mosquitoes. This suggested that the differences in the transcript profiles among the strains are heritable due to a selection process and are not caused by immediate insecticide exposure. Reversal of temephos resistance was demonstrated and, importantly, there was a positive correlation between a decrease in the resistance ratio and an accompanying decrease in the expression levels of previously over expressed genes. Some of the genes identified here have also been implicated in metabolic resistance in other mosquito species and insecticide resistant populations of Ae. aegypti.

Conclusions/Significance

The identification of gene expression signatures associated to insecticide resistance and their suppression could greatly aid the development of improved strategies of vector control.

Existence of the rdl mutant alleles among the anopheles malaria vector in Indonesia

BACKGROUND

The gamma-aminobutyric acid (GABA) receptor-chloride channel complex is known to be the target site of dieldrin, a cyclodiene insecticide. GABA-receptors, with a naturally occurring amino acid substitution, A302S/G in the putative ion-channel lining region, confer resistance to cyclodiene insecticides that includes aldrin, chlordane, dieldrin, heptachlor, endrin and endosulphan.

METHODS

A total of 154 mosquito samples from 10 provinces of malaria-endemic areas across Indonesia (Aceh, North Sumatra, Bangka Belitung, Lampung, Central Java, East Nusa Tenggara, West Nusa Tenggara, West Sulawesi, Molucca and North Molucca) were obtained and identified by species, using morphological characteristic. The DNA was individually extracted using chelex-ion exchanger and the DNA obtained was used for analyses using sequencing method.

RESULTS

Molecular analysis indicated 11% of the total 154 Anopheles samples examined, carried Rdl mutant alleles. All of the alleles were found in homozygous form. Rdl 302S allele was observed in Anopheles vagus (from Central Java, Lampung, and West Nusa Tenggara), Anopheles aconitus (from Central Java), Anopheles barbirostris (from Central Java and Lampung), Anopheles sundaicus (from North Sumatra and Lampung), Anopheles nigerrimus (from North Sumatra), whereas the 302 G allele was only found in Anopheles farauti from Molucca.

CONCLUSION

The existence of the Rdl mutant allele indicates that, either insecticide pressure on the Anopheles population in these areas might still be ongoing (though not directly associated with the malaria control programme) or that the mutant form of the Rdl allele is relatively stable in the absence of insecticide. Nonetheless, the finding suggests that integrated pest management is warranted in malaria-endemic areas where insecticides are widely used for other purposes.

The rph1 gene is a common contributor to the evolution of phosphine resistance in independent field isolates of Rhyzopertha dominica

Phosphine is the only economically viable fumigant for routine control of insect pests of stored food products, but its continued use is now threatened by the world-wide emergence of high-level resistance in key pest species. Phosphine has a unique mode of action relative to well-characterised contact pesticides. Similarly, the selective pressures that lead to resistance against field sprays differ dramatically from those encountered during fumigation. The consequences of these differences have not been investigated adequately. We determine the genetic basis of phosphine resistance in Rhyzopertha dominica strains collected from New South Wales and South Australia and compare this with resistance in a previously characterised strain from Queensland. The resistance levels range from 225 and 100 times the baseline response of a sensitive reference strain. Moreover, molecular and phenotypic data indicate that high-level resistance was derived independently in each of the three widely separated geographical regions. Despite the independent origins, resistance was due to two interacting genes in each instance. Furthermore, complementation analysis reveals that all three strains contain an incompletely recessive resistance allele of the autosomal rph1 resistance gene. This is particularly noteworthy as a resistance allele at rph1 was previously proposed to be a necessary first step in the evolution of high-level resistance. Despite the capacity of phosphine to disrupt a wide range of enzymes and biological processes, it is remarkable that the initial step in the selection of resistance is so similar in isolated outbreaks.

Structural and regulatory evolution of cellular electrophysiological systems

Cellular electrophysiological systems, like developmental systems, appear to evolve primarily by means of regulatory evolution. It is suggested that electrophysiological systems share two key features with developmental systems that account for this dependence on regulatory evolution. For both systems, structural evolution has the potential to create significant problems of pleiotropy and both systems are predominantly computational in nature. It is concluded that the relative balance of physical and computational tasks that a biological system has to perform, combined with the probability that these tasks may have to change significantly during the course of evolution, will be major factors in determining the relative mix of regulatory and structural evolution that is observed for a given system. Physiological systems that directly interface with the environment will almost always perform some low-level physical task. In the majority of cases this will require evolution of protein function in order for the tasks themselves to evolve. For complex physiological systems a large fraction of their function will be devoted to high-level control functions that are predominantly computational in nature. In most cases regulatory evolution will be sufficient in order for these computational tasks to evolve.

Adaptive divergence in experimental populations of Pseudomonas fluorescens. IV. Genetic constraints guide evolutionary trajectories in a parallel adaptive radiation

The capacity for phenotypic evolution is dependent upon complex webs of functional interactions that connect genotype and phenotype. Wrinkly spreader (WS) genotypes arise repeatedly during the course of a model Pseudomonas adaptive radiation. Previous work showed that the evolution of WS variation was explained in part by spontaneous mutations in wspF, a component of the Wsp-signaling module, but also drew attention to the existence of unknown mutational causes. Here, we identify two new mutational pathways (Aws and Mws) that allow realization of the WS phenotype: in common with the Wsp module these pathways contain a di-guanylate cyclase-encoding gene subject to negative regulation. Together, mutations in the Wsp, Aws, and Mws regulatory modules account for the spectrum of WS phenotype-generating mutations found among a collection of 26 spontaneously arising WS genotypes obtained from independent adaptive radiations. Despite a large number of potential mutational pathways, the repeated discovery of mutations in a small number of loci (parallel evolution) prompted the construction of an ancestral genotype devoid of known (Wsp, Aws, and Mws) regulatory modules to see whether the types derived from this genotype could converge upon the WS phenotype via a novel route. Such types-with equivalent fitness effects-did emerge, although they took significantly longer to do so. Together our data provide an explanation for why WS evolution follows a limited number of mutational pathways and show how genetic architecture can bias the molecular variation presented to selection.

Is genetic evolution predictable?

Ever since the integration of Mendelian genetics into evolutionary biology in the early 20th century, evolutionary geneticists have for the most part treated genes and mutations as generic entities. However, recent observations indicate that all genes are not equal in the eyes of evolution. Evolutionarily relevant mutations tend to accumulate in hotspot genes and at specific positions within genes. Genetic evolution is constrained by gene function, the structure of genetic networks, and population biology. The genetic basis of evolution may be predictable to some extent, and further understanding of this predictability requires incorporation of the specific functions and characteristics of genes into evolutionary theory.

Acetylcholinesterase genes in the basal Hexapod Orchesella villosa

Acetylcholinesterase (AChE) is a key enzyme of the cholinergic nerve system. Of the two forms found in insects, the predominant one is active in the synapses and is the target of organophosphate and carbamate insecticides, while the role of the second is currently unknown. Two acetylcholinesterase cDNAs from the basal hexapod Orchesella villosa have been characterized and compared with others reported form insects. One form conforms well to the typical structure, while the other is characterized by an unusual 3' region. No amino acid mutation could be directly associated with known resistance mutations in other insect species or to a clear signal of selection in the distribution of alleles, although the action of some population process is suggested.

The dominant cold-sensitive Out-cold mutants of Drosophila melanogaster have novel missense mutations in the voltage-gated sodium channel gene paralytic

Here we report the molecular characterization of Out-cold (Ocd) mutants of Drosophila melanogaster, which produce a dominant, X-linked, cold-sensitive paralytic phenotype. From its initial 1.5-Mb cytological location within 13F1-16A2, P-element and SNP mapping reduced the Ocd critical region to <100 kb and to six candidate genes: hangover, CG9947, CG4420, eIF2a, Rbp2, and paralytic (para). Complementation testing with para null mutations strongly suggests Ocd and para are allelic, as does gene rescue of Ocd semilethality with a wild-type para transgene. Pesticide resistance and electrophysiological phenotypes of Ocd mutants support this conclusion. The para gene encodes a voltage-gated sodium channel. Sequencing the Ocd lines revealed mutations within highly conserved regions of the para coding sequence, in the transmembrane segment S6 of domain III (I1545M and T1551I), and in the linker between domains III and IV (G1571R), the location of the channel inactivation gate. The G1571R mutation is of particular interest as mutations of the orthologous residue (G1306) in the human skeletal muscle sodium channel gene SCN4A are associated with cases of periodic paralysis and myotonia, including the human cold-sensitive disorder paramyotonia congenita. The mechanisms by which sodium channel mutations cause cold sensitivity are not well understood. Therefore, in the absence of suitable vertebrate models, Ocd provides a system in which genetic, molecular, physiological, and behavioral tools can be exploited to determine mechanisms underlying sodium channel periodic paralyses.

A small deletion in the olive fly acetylcholinesterase gene associated with high levels of organophosphate resistance

Organophosphate resistance in the olive fly was previously shown to associate with two point mutations in the ace gene. The frequency of these mutations was monitored in Bactrocera oleae individuals of increasing resistance. In spite of the difference in resistance among the individuals, there was no correlation between mutation frequencies and resistance level, indicating that other factors may contribute to this variation. The search for additional mutations in the ace gene of highly resistant insects revealed a small deletion at the carboxyl terminal of the protein (termed Delta3Q). Significant correlation was shown between the mutation frequency and resistance level in natural populations. In addition, remaining activity of acetylcholinesterase enzyme (AChE) after dimethoate inhibition was higher in genotypes carrying the mutation. These results strongly suggest a role of Delta3Q in high levels of organophosphate (OP) resistance. Interestingly, the carboxyl terminal of AChE is normally cleaved and substituted by a glycosylphosphatidylinositol (GPI) anchor. We hypothesize that Delta3Q may improve GPI anchoring, thus increasing the amount of AChE that reaches the synaptic cleft. In this way, despite the presence of insecticide, enough enzyme would remain in the cleft for its normal role of acetylcholine hydrolysis, allowing the insect to survive. This provides a previously un-described mechanism of resistance.

Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family

Approximately 75% of vertebrate proteins belong to protein families encoded by multiple evolutionarily related genes, a pattern that emerged as a result of gene and genome duplications over the course of vertebrate evolution. In families of genes with similar or related functions, adaptation to a strong selective agent should involve multiple adaptive changes across the entire gene family. However, we know of no evolutionary studies that have explicitly addressed this point. Here, we show how 4 taxonomically diverse species of pufferfishes (Tetraodontidae) each evolved resistance to the guanidinium toxins tetrodotoxin (TTX) and saxitoxin (STX) via parallel amino acid replacements across all 8 sodium channels present in teleost fish genomes. This resulted in diverse suites of coexisting sodium channel types that all confer varying degrees of toxin resistance, yet show remarkable convergence among genes and phylogenetically diverse species. Using site-directed mutagenesis and expression of a vertebrate sodium channel, we also demonstrate that resistance to TTX/STX is enhanced up to 15-fold by single, frequently observed replacements at 2 sites that have not previously been implicated in toxin binding but show similar or identical replacements in pufferfishes and in distantly related vertebrate and nonvertebrate animals. This study presents an example of natural selection acting upon a complete gene family, repeatedly arriving at a diverse but limited number of adaptive changes within the same genome. To be maximally informative, we suggest that future studies of molecular adaptation should consider all functionally similar paralogs of the affected gene family.

Molecular mechanisms of metabolic resistance to synthetic and natural xenobiotics

Xenobiotic resistance in insects has evolved predominantly by increasing the metabolic capability of detoxificative systems and/or reducing xenobiotic target site sensitivity. In contrast to the limited range of nucleotide changes that lead to target site insensitivity, many molecular mechanisms lead to enhancements in xenobiotic metabolism. The genomic changes that lead to amplification, overexpression, and coding sequence variation in the three major groups of genes encoding metabolic enzymes, i.e., cytochrome P450 monooxygenases (P450s), esterases, and glutathione-S-transferases (GSTs), are the focus of this review. A substantial number of the adaptive genomic changes associated with insecticide resistance that have been characterized to date are transposon mediated. Several lines of evidence suggest that P450 genes involved in insecticide resistance, and perhaps insecticide detoxification genes in general, may share an evolutionary association with genes involved in allelochemical metabolism. Differences in the selective regime imposed by allelochemicals and insecticides may account for the relative importance of regulatory or structural mutations in conferring resistance.