Copy Number Variation and Transposable Elements Feature in Recent, Ongoing Adaptation at the Cyp6g1 Locus

The increased transcription of the Cyp6g1 gene of Drosophila melanogaster, and consequent resistance to insecticides such as DDT, is a widely cited example of adaptation mediated by cis-regulatory change. A fragment of an Accord transposable element inserted upstream of the Cyp6g1 gene is causally associated with resistance and has spread to high frequencies in populations around the world since the 1940s. Here we report the existence of a natural allelic series at this locus of D. melanogaster, involving copy number variation of Cyp6g1, and two additional transposable element insertions (a P and an HMS-Beagle). We provide evidence that this genetic variation underpins phenotypic variation, as the more derived the allele, the greater the level of DDT resistance. Tracking the spatial and temporal patterns of allele frequency changes indicates that the multiple steps of the allelic series are adaptive. Further, a DDT association study shows that the most resistant allele, Cyp6g1-[BP], is greatly enriched in the top 5% of the phenotypic distribution and accounts for ∼16% of the underlying phenotypic variation in resistance to DDT. In contrast, copy number variation for another candidate resistance gene, Cyp12d1, is not associated with resistance. Thus the Cyp6g1 locus is a major contributor to DDT resistance in field populations, and evolution at this locus features multiple adaptive steps occurring in rapid succession

Sexual conflict maintains variation at an insecticide resistance locus

Background: The maintenance of genetic variation through sexually antagonistic selection is controversial, partly because specific sexually-antagonistic alleles have not been identified. The Drosophila DDT resistance allele (DDT-R) is an exception. This allele increases female fitness, but simultaneously decreases male fitness, and it has been suggested that this sexual antagonism could explain why polymorphism was maintained at the locus prior to DDT use. We tested this possibility using a genetic model and then used evolving fly populations to test model predictions. Results: Theory predicted that sexual antagonism is able to maintain genetic variation at this locus, hence explaining why DDT-R did not fix prior to DDT use despite increasing female fitness, and experimentally evolving fly populations verified theoretical predictions. Conclusions: This demonstrates that sexually antagonistic selection can maintain genetic variation and explains the DDT-R frequencies observed in nature

Use of Mutagenesis, Genetic Mapping and Next Generation Transcriptomics to Investigate Insecticide Resistance Mechanisms

Insecticide resistance is a worldwide problem with major impact on agriculture and human health. Understanding the underlying molecular mechanisms is crucial for the management of the phenomenon; however, this information often comes late with respect to the implementation of efficient counter-measures, particularly in the case of metabolism-based resistance mechanisms. We employed a genome-wide insertional mutagenesis screen to Drosophila melanogaster, using a Minos-based construct, and retrieved a line (MiT[w−]3R2) resistant to the neonicotinoid insecticide Imidacloprid. Biochemical and bioassay data indicated that resistance was due to increased P450 detoxification. Deep sequencing transcriptomic analysis revealed substantial over- and under-representation of 357 transcripts in the resistant line, including statistically significant changes in mixed function oxidases, peptidases and cuticular proteins. Three P450 genes (Cyp4p2, Cyp6a2 and Cyp6g1) located on the 2R chromosome, are highly up-regulated in mutant flies compared to susceptible Drosophila. One of them (Cyp6g1) has been already described as a major factor for Imidacloprid resistance, which validated the approach. Elevated expression of the Cyp4p2 was not previously documented in Drosophila lines resistant to neonicotinoids. In silico analysis using the Drosophila reference genome failed to detect transcription binding factors or microRNAs associated with the over-expressed Cyp genes. The resistant line did not contain a Minos insertion in its chromosomes, suggesting a hit-and-run event, i.e. an insertion of the transposable element, followed by an excision which caused the mutation. Genetic mapping placed the resistance locus to the right arm of the second chromosome, within a ∼1 Mb region, where the highly up-regulated Cyp6g1 gene is located. The nature of the unknown mutation that causes resistance is discussed on the basis of these results

Pleiotropic Effects of DDT Resistance on Male Size and Behaviour

Understanding the evolution and spread of insecticide resistance requires knowing the relative fitness of resistant organisms. In the absence of insecticides, resistance is predicted to be costly. The Drosophila melanogaster DDT resistance allele (DDT-R) is associated with a male mating cost. This could be because resistant males are generally smaller, but DDT-R may also alter courtship behaviours. Here we tested for body size and courtship effects of DDT-R on mating success in competitive and non-competitive mating trials respectively. We also assessed relative aggression in resistant and susceptible males because aggression can also influence mating success. While the effect of DDT-R on male size partly contributed to reduced mating success, resistant males also had lower rates of courtship and were less aggressive than susceptible males. These differences contribute to the observed DDT-R mating costs. Additionally, these pleiotropic effects of DDT-R are consistent with the history and spread of resistance alleles in nature

Fungicide-Driven Evolution and Molecular Basis of Multidrug Resistance in Field Populations of the Grey Mould Fungus Botrytis cinerea

The grey mould fungus Botrytis cinerea causes losses of commercially important fruits, vegetables and ornamentals worldwide. Fungicide treatments are effective for disease control, but bear the risk of resistance development. The major resistance mechanism in fungi is target protein modification resulting in reduced drug binding. Multiple drug resistance (MDR) caused by increased efflux activity is common in human pathogenic microbes, but rarely described for plant pathogens. Annual monitoring for fungicide resistance in field isolates from fungicide-treated vineyards in France and Germany revealed a rapidly increasing appearance of B. cinerea field populations with three distinct MDR phenotypes. All MDR strains showed increased fungicide efflux activity and overexpression of efflux transporter genes. Similar to clinical MDR isolates of Candida yeasts that are due to transcription factor mutations, all MDR1 strains were shown to harbor activating mutations in a transcription factor (Mrr1) that controls the gene encoding ABC transporter AtrB. MDR2 strains had undergone a unique rearrangement in the promoter region of the major facilitator superfamily transporter gene mfsM2, induced by insertion of a retrotransposon-derived sequence. MDR2 strains carrying the same rearranged mfsM2 allele have probably migrated from French to German wine-growing regions. The roles of atrB, mrr1 and mfsM2 were proven by the phenotypes of knock-out and overexpression mutants. As confirmed by sexual crosses, combinations of mrr1 and mfsM2 mutations lead to MDR3 strains with higher broad-spectrum resistance. An MDR3 strain was shown in field experiments to be selected against sensitive strains by fungicide treatments. Our data document for the first time the rising prevalence, spread and molecular basis of MDR populations in a major plant pathogen in agricultural environments. These populations will increase the risk of grey mould rot and hamper the effectiveness of current strategies for fungicide resistance management

Transient exposure to low levels of insecticide affects metabolic networks of honeybee larvae

The survival of a species depends on its capacity to adjust to changing environmental conditions, and new stressors. Such new, anthropogenic stressors include the neonicotinoid class of crop-protecting agents, which have been implicated in the population declines of pollinating insects, including honeybees (Apis mellifera). The low-dose effects of these compounds on larval development and physiological responses have remained largely unknown. Over a period of 15 days, we provided syrup tainted with low levels (2 µg/L−1) of the neonicotinoid insecticide imidacloprid to beehives located in the field. We measured transcript levels by RNA sequencing and established lipid profiles using liquid chromatography coupled with mass spectrometry from worker-bee larvae of imidacloprid-exposed (IE) and unexposed, control (C) hives. Within a catalogue of 300 differentially expressed transcripts in larvae from IE hives, we detect significant enrichment of genes functioning in lipid-carbohydrate-mitochondrial metabolic networks. Myc-involved transcriptional response to exposure of this neonicotinoid is indicated by overrepresentation of E-box elements in the promoter regions of genes with altered expression. RNA levels for a cluster of genes encoding detoxifying P450 enzymes are elevated, with coordinated downregulation of genes in glycolytic and sugar-metabolising pathways. Expression of the environmentally responsive Hsp90 gene is also reduced, suggesting diminished buffering and stability of the developmental program. The multifaceted, physiological response described here may be of importance to our general understanding of pollinator health. Muscles, for instance, work at high glycolytic rates and flight performance could be impacted should low levels of this evolutionarily novel stressor likewise induce downregulation of energy metabolising genes in adult pollinators

Genomic Footprints of Selective Sweeps from Metabolic Resistance to Pyrethroids in African Malaria Vectors Are Driven by Scale up of Insecticide-Based Vector Control

Insecticide resistance in mosquito populations threatens recent successes in malaria prevention. Elucidating patterns of genetic structure in malaria vectors to predict the speed and direction of the spread of resistance is essential to get ahead of the `resistance curve' and to avert a public health catastrophe. Here, applying a combination of microsatellite analysis, whole genome sequencing and targeted sequencing of a resistance locus, we elucidated the continent-wide population structure of a major African malaria vector, Anopheles funestus. We identified a major selective sweep in a genomic region controlling cytochrome P450-based metabolic resistance conferring high resistance to pyrethroids. This selective sweep occurred since 2002, likely as a direct consequence of scaled up vector control as revealed by whole genome and fine-scale sequencing of pre- and post-intervention populations. Fine-scaled analysis of the pyrethroid resistance locus revealed that a resistanceassociated allele of the cytochrome P450 monooxygenase CYP6P9a has swept through southern Africa to near fixation, in contrast to high polymorphism levels before interventions, conferring high levels of pyrethroid resistance linked to control failure. Population structure analysis revealed a barrier to gene flow between southern Africa and other areas, which may prevent or slow the spread of the southern mechanism of pyrethroid resistance to other regions. By identifying a genetic signature of pyrethroid-based interventions, we have demonstrated the intense selective pressure that control interventions exert on mosquito populations. If this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised

Population transcriptomics of Drosophila melanogaster females

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

Multiple Origins of Knockdown Resistance Mutations in the Afrotropical Mosquito Vector Anopheles gambiae

How often insecticide resistance mutations arise in natural insect populations is a fundamental question for understanding the evolution of resistance and also for modeling its spread. Moreover, the development of resistance is regarded as a favored model to study the molecular evolution of adaptive traits. In the malaria vector Anopheles gambiae two point mutations (L1014F and L1014S) in the voltage-gated sodium channel gene, that confer knockdown resistance (kdr) to DDT and pyrethroid insecticides, have been described. In order to determine whether resistance alleles result from single or multiple mutation events, genotyping of the kdr locus and partial sequencing of the upstream intron-1 was performed on a total of 288 A. gambiae S-form collected from 28 localities in 15 countries. Knockdown resistance alleles were found to be widespread in West Africa with co-occurrence of both 1014S and 1014F in West-Central localities. Differences in intron-1 haplotype composition suggest that kdr alleles may have arisen from at least four independent mutation events. Neutrality tests provided evidence for a selective sweep acting on this genomic region, particularly in West Africa. The frequency and distribution of these kdr haplotypes varied geographically, being influenced by an interplay between different mutational occurrences, gene flow and local selection. This has important practical implications for the management and sustainability of malaria vector control programs