Genome-wide sequencing and an open reading frame analysis of dichlorodiphenyltrichloroethane (DDT) susceptible (91-C) and resistant (91-R) Drosophila melanogaster laboratory populations

Phenotypic Plasticity of Common Wasps in an Industrially Polluted Environment in Southwestern Finland

Simple Summary

Social insects are ecologically and economically important as ecosystem engineers, agricultural pest predators, pollinators, and seed dispersers. Many of the vespid wasps are social insects. Our study species, Common wasp Vespula vulgaris, is native to Finland and classified as invasive in some other parts of the world. The Common wasp have conspicuous yellow and black pigmentation. Their functions and activities in the environment expose the species to environmental pollutants and this study assessed the effect of heavy-metals on common wasps collected from the vicinity of a metal smelter in southwestern Finland. The samples collected were analyzed using various methods such as color morph categorization, electron microscopy, metal analysis, and energy dispersive X-ray analysis (EDX). The methods were used to understand the effects of metal pollution on the species and the adaptive response. Our results indicated phenotypic variation between common wasp samples across the pollution gradient and an adaptive melanin encapsulation process.

Abstract

Insects vary in the degree of their adaptability to environmental contamination. Determining the responses with phenotypic plasticity in ecologically important species in polluted environments will ease further conservation and control actions. Here, we investigated morphological characteristics such as body size, body mass, and color of the common wasp Vespula vulgaris in an industrially polluted environment, considering different levels of metal pollution, and we studied the localization of contaminants in the guts of wasps. We revealed some differences in morphological characteristics and melanization of wasps collected in habitats with high, moderate, and low levels of pollution. The results indicated that V. vulgaris from highly polluted environments had reduced melanin pigmentation on the face but increased melanin pigmentation on the 2nd tergite of the abdomen. In addition, with transmission electron microscopy (TEM) and energy dispersive X-ray analysis (EDX), we found metal particles from the midgut of wasps originating from the polluted environment. Most of the particles were encapsulated with melanin pigment. This finding confirmed that in wasps, ingested metal particles are accumulated in guts and covered by melanin layers. Our data suggest that wasps can tolerate metal contamination but respond phenotypically with modification of their size, coloration, and probably with the directions of the melanin investments (immunity or coloration). Thus, in industrially polluted areas, wasps might probably survive by engaging phenotypic plasticity with no significant or visible impact on the population.

Dietary antioxidant vitamin C influences the evolutionary path of insecticide resistance in Drosophila melanogaster

Herbivorous insects encounter a variety of toxic environmental substances ranging from ingested plant defensive compounds to human-introduced insecticidal agents. Dietary antioxidants are known to reduce the negative physiological impacts of toxins in mammalian systems through amelioration of reactive oxygen-related cellular damage. The analogous impacts to insects caused by multigenerational exposure to pesticides and the effects on adaptive responses within insect populations, however, are currently unknown. To address these research gaps, we used Drosophila as a model system to explore adaptive phenotypic responses to acute dichlorodiphenyltrichloroethane (DDT) exposure in the presence of the dietary antioxidant vitamin C and to examine the structural genomic consequences of this exposure. DDT resistance increased significantly among four replicates exposed to a low concentration of DDT for 10 generations. In contrast, dietary intake of vitamin C significantly reduced DDT resistance after mutigenerational exposure to the same concentration of DDT. As to the genomic consequences, no significant differences were predicted in overall nucleotide substitution rates across the genome between any of the treatments. Despite this, replicates exposed to a low concentration of DDT without vitamin C showed the highest number of synonymous and non-synonymous variants (3196 in total), followed by the DDT plus vitamin C (1174 in total), and vitamin C alone (728 in total) treatments. This study demonstrates the potential role of diet (specifically, antioxidant intake) on adaptive genome responses, and thus on the evolution of pesticide resistance within insect populations.

A review of DDT resistance as it pertains to the 91-C and 91-R strains in Drosophila melanogaster

While insecticide resistance presents a challenge for those intent on controlling insect populations, these challenges have also generated a set of tools that can be used to ask fundamental biological questions about that resistance. Numerous species of insects have evolved resistance to multiple classes of insecticides. Each one of these species and their respective resistant populations represent a potential tool for understanding the molecular basis of the evolution of resistance. However, in-laboratory maintenance of resistant insect populations (and their comparative susceptible populations) suitable for asking the needed set of questions around the molecular consequences of long-term pesticide exposure requires a significant, in places prohibitive, level of resources. Drosophila melanogaster (hereafter referred to as Drosophila) is a model insect system with populations easily selected with pesticides and readily maintainable over decades. Even within Drosophila, however, few populations exist where long-term pesticide selection has occurred along with contrasting non-selected population. As such, the Drosophila 91-C and 91-R populations, which exhibit insecticide resistance to DDT (91-R), compared to a non-selection population (91-C), represent a unique resource for the study of high level DDT resistance. Moreover, with the availability of "omics" technologies over the past several decades, this paired population has emerged as a useful tool for understanding both the molecular basis of pesticide resistance and the molecular consequences of long-term pesticide exposure. In this review, we summarize the studies with these aforementioned populations over the past several decades, addressing what has been learned from these efforts.

Cytochrome P450s Cyp4p1 and Cyp4p2 associated with the DDT tolerance in the Drosophila melanogaster strain 91-R

Cytochrome P450s are part of a super-gene family that has undergone gene duplication, divergence, over-expression and, in some cases, loss of function. One such case is the 91-R and 91-C strains of common origin, in Drosophila melanogaster, whereby 91-R (DDT resistant strain) overexpresses Cyp4p1 and Cyp4p2 and both genes are lost in 91-C (DDT susceptible strain). In this study, we used a comparative approach to demonstrate that transcription of Cyp4p1 and Cyp4p2 were constitutively up-regulated in the Drosophila melanogaster strain 91-R as compared to another DDT susceptible strain Canton-S which does not have a loss of function of these genes. Furthermore, significantly increased expression of Cyp4p1 and Cyp4p2 was induced in 91-R in response to sublethal DDT exposure, however, such induction did not occur in the DDT treated Canton-S. Additionally, fixed nucleotide variation within putative transcription factor binding sites of Cyp4p1 and Cyp4p2 promoters were observed between 91-R and Canton-S, however, their impact on transcription remains to be determined. Two GAL4/UAS transgenic strains with integrated heat shock-inducible Cyp4p1- or Cyp4p2-RNAi constructs within wild-type genetic backgrounds were developed. Following heat shock induction of Cyp4p1 and Cyp4p2 knockdown, these transgenic lines showed increased DDT mortality as compared to their corresponding non-heat shock controls. These results provide a functional link of Cyp4p1 and Cyp4p2 in conferring tolerance to DDT exposure.

Comparative CYP-omic analysis between the DDT-susceptible and -resistant Drosophila melanogaster strains 91-C and 91-R

BACKGROUND

Cytochrome P450 monooxygenases (P450s) are involved in the biosynthesis of endogenous intracellular compounds and the metabolism of xenobiotics, including chemical insecticides. We investigated the structural and expression level variance across all P450 genes with respect to the evolution of insecticide resistance under multigenerational dichlorodiphenyltrichloroethane (DDT) selection.

RESULTS

RNA-sequencing (RNA-seq) and reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) indicated that the transcript levels of seven P450 genes were significantly up-regulated and three P450 genes were down-regulated in the DDT-resistant strain 91-R, as compared to the control strain 91-C. The overexpression of Cyp6g1 was associated with the presence of an Accord and an HMS-Beagle element insertion in the 5' upstream region in conjunction with copy number variation in the 91-R strain, but not in the 91-C strain. A total of 122 (50.2%) fixed nonsynonymous (amino acid-changing) mutations were found between 91-C and 91-R, and 20 (8.2%) resulted in amino acid changes within functional domains. Three P450 proteins were truncated as a result of premature stop codons and fixed between strains.

CONCLUSION

Our results demonstrate that a combination of changes in P450 protein-coding regions and transcript levels are possibly associated with DDT resistance, and thereby suggest that selection for variant function may occur within this gene family in response to chronic DDT exposure. © 2018 Society of Chemical Industry.

Variation in Mitochondria-Derived Transcript Levels Associated With DDT Resistance in the 91-R Strain of Drosophila melanogaster (Diptera: Drosophilidae)

The organochloride insecticide dichlorodiphenyltrichloroethane (DDT) and its metabolites can increase cellular levels of reactive oxygen species (ROS), cause mitochondrial dysfunction, and induce apoptosis. The highly DDT-resistant Drosophila melanogaster Meigen 1830 (Drosophila) strain, 91-R, and its susceptible control, 91-C, were used to investigate functional and structural changes among mitochondrial-derived pathways. Resequencing of mitochondrial genomes (mitogenomes) detected no structural differences between 91-R and 91-C, whereas RNA-seq suggested the differential expression of 221 mitochondrial-associated genes. Reverse transcriptase-quantitative PCR validation of 33 candidates confirmed that transcripts for six genes (Cyp12d1-p, Cyp12a4, cyt-c-d, COX5BL, COX7AL, CG17140) were significantly upregulated and two genes (Dif, Rel) were significantly downregulated in 91-R. Among the upregulated genes, four genes are duplicated within the reference genome (cyt-c-d, CG17140, COX5BL, and COX7AL). The predicted functions of the differentially expressed genes, or known functions of closely related genes, suggest that 91-R utilizes existing ROS regulation pathways of the mitochondria to combat increased ROS levels from exposure to DDT. This study represents, to our knowledge, the initial investigation of mitochondrial genome sequence variants and functional adaptations in responses to intense DDT selection and provides insights into potential adaptations of ROS management associated with DDT selection in Drosophila.

Changes in Neuronal Signaling and Cell Stress Response Pathways are Associated with a Multigenic Response of Drosophila melanogaster to DDT Selection

The adaptation of insect populations to insecticidal control is a continual threat to human health and sustainable agricultural practices, but many complex genomic mechanisms involved in this adaption remain poorly understood. This study applied a systems approach to investigate the interconnections between structural and functional variance in response to dichlorodiphenyltrichloroethane (DDT) within the Drosophila melanogaster strain 91-R. Directional selection in 6 selective sweeps coincided with constitutive gene expression differences in DDT resistant flies, including the most highly upregulated transcript, Unc-115 b, which plays a central role in axon guidance, and the most highly downregulated transcript, the angiopoietin-like CG31832, which is involved in directing vascular branching and dendrite outgrowth but likely may be under trans-regulatory control. Direct functions and protein–protein interactions mediated by differentially expressed transcripts control changes in cell migration, signal transduction, and gene regulatory cascades that impact the nervous system. Although changes to cellular stress response pathways involve 8 different cytochrome P450s, stress response, and apoptosis is controlled by a multifacetted regulatory mechanism. These data demonstrate that DDT selection in 91-R may have resulted in genome-wide adaptations that impacts genetic and signal transduction pathways that converge to modify stress response, cell survival, and neurological functions. This study implicates the involvement of a multigenic mechanism in the adaptation to a chemical insecticide, which impact interconnected regulatory cascades. We propose that DDT selection within 91-R might act systemically, wherein pathway interactions function to reinforce the epistatic effects of individual adaptive changes on an additive or nonadditive basis.

Differentially expressed microRNAs associated with changes of transcript levels in detoxification pathways and DDT-resistance in the Drosophila melanogaster strain 91-R

Dichloro-diphenyl-trichloroethane (DDT) resistance among arthropod species is a model for understanding the molecular adaptations in response to insecticide exposures. Previous studies reported that DDT resistance may involve one or multiple detoxification genes, such as cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs), esterases, and ATP binding cassette (ABC) transporters, or changes in the voltage-sensitive sodium channel. However, the possible involvement of microRNAs (miRNAs) in the post-transcriptional regulation of genes associated with DDT resistance in the Drosophila melanogaster strain 91-R remains poorly understood. In this study, the majority of the resulting miRNAs discovered in small RNA libraries from 91-R and the susceptible control strain, 91-C, ranged from 16-25 nt, and contained 163 precursors and 256 mature forms of previously-known miRNAs along with 17 putative novel miRNAs. Quantitative analyses predicted the differential expression of ten miRNAs between 91-R and 91-C, and, based on Gene Ontology and pathway analysis, these ten miRNAs putatively target transcripts encoding proteins involved in detoxification mechanisms. RT-qPCR validated an inverse correlation between levels of differentially-expressed miRNAs and their putatively targeted transcripts, which implies a role of these miRNAs in the differential regulation of detoxification pathways in 91-R compared to 91-C. This study provides evidence associating the differential expression of miRNAs in response to multigenerational DDT selection in Drosophila melanogaster and provides important clues for understanding the possible roles of miRNAs in mediating insecticide resistance traits.

Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel

Insecticide resistance is considered a classic model of microevolution, where a strong selective agent is applied to a large natural population, resulting in a change in frequency of alleles that confer resistance. While many insecticide resistance variants have been characterized at the gene level, they are typically single genes of large effect identified in highly resistant pest species. In contrast, multiple variants have been implicated in DDT resistance in Drosophila melanogaster; however, only the Cyp6g1 locus has previously been shown to be relevant to field populations. Here we use genome-wide association studies (GWAS) to identify DDT-associated polygenes and use selective sweep analyses to assess their adaptive significance. We identify and verify two candidate DDT resistance loci. A largely uncharacterized gene, CG10737, has a function in muscles that ameliorates the effects of DDT, while a putative detoxifying P450, Cyp6w1, shows compelling evidence of positive selection.

Selective sweep analysis in the genomes of the 91-R and 91-C Drosophila melanogaster strains reveals few of the 'usual suspects' in dichlorodiphenyltrichloroethane (DDT) resistance

Adaptation of insect phenotypes for survival after exposure to xenobiotics can result from selection at multiple loci with additive genetic effects. To the authors' knowledge, no selective sweep analysis has been performed to identify such loci in highly dichlorodiphenyltrichloroethane (DDT) resistant insects. Here we compared a highly DDT resistant phenotype in the Drosophila melanogaster (Drosophila) 91-R strain to the DDT susceptible 91-C strain, both of common origin. Whole genome re-sequencing data from pools of individuals was generated separately for 91-R and 91-C, and mapped to the reference Drosophila genome assembly (v. 5.72). Thirteen major and three minor effect chromosome intervals with reduced nucleotide diversity (π) were identified only in the 91-R population. Estimates of Tajima's D (D) showed corresponding evidence of directional selection in these same genome regions of 91-R, however, no similar reductions in π or D estimates were detected in 91-C. An overabundance of non-synonymous proteins coding to synonymous changes were identified in putative open reading frames associated with 91-R. Except for NinaC and Cyp4g1, none of the identified genes were the 'usual suspects' previously observed to be associated with DDT resistance. Additionally, up-regulated ATP-binding cassette transporters have been previously associated with DDT resistance; however, here we identified a structurally altered MDR49 candidate resistance gene. The remaining fourteen genes have not previously been shown to be associated with DDT resistance. These results suggest hitherto unknown mechanisms of DDT resistance, most of which have been overlooked in previous transcriptional studies, with some genes having orthologs in mammals.