Molecular evolution of glutathione S-transferases in the genus Drosophila

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Identification and Characterization of Glutathione S-transferase Genes in Spodoptera frugiperda (Lepidoptera: Noctuidae) under Insecticides Stress

Insect glutathione S-transferases (GSTs) serve critical roles in insecticides and other forms of xenobiotic chemical detoxification. The fall armyworm, Spodoptera frugiperda (J. E. Smith), is a major agricultural pest in several countries, especially Egypt. This is the first study to identify and characterize GST genes in S. frugiperda under insecticidal stress. The present work evaluated the toxicity of emamectin benzoate (EBZ) and chlorantraniliprole (CHP) against the third-instar larvae of S. frugiperda using the leaf disk method. The LC₅₀ values of EBZ and CHP were 0.029 and 1.250 mg/L after 24 h of exposure. Moreover, we identified 31 GST genes, including 28 cytosolic and 3 microsomal SfGSTs from a transcriptome analysis and the genome data of S. frugiperda. Depending on the phylogenetic analysis, sfGSTs were divided into six classes (delta, epsilon, omega, sigma, theta, and microsomal). Furthermore, we investigated the mRNA levels of 28 GST genes using qRT-PCR under EBZ and CHP stress in the third-instar larvae of S. frugiperda. Interestingly, SfGSTe10 and SfGSTe13 stood out with the highest expression after the EBZ and CHP treatments. Finally, a molecular docking model was constructed between EBZ and CHP using the most upregulated genes (SfGSTe10 and SfGSTe13) and the least upregulated genes (SfGSTs1 and SfGSTe2) of S. frugiperda larvae. The molecular docking study showed EBZ and CHP have a high binding affinity with SfGSTe10, with docking energy values of -24.41 and -26.72 kcal/mol, respectively, and sfGSTe13, with docking energy values of -26.85 and -26.78 kcal/mol, respectively. Our findings are important for understanding the role of GSTs in S. frugiperda regarding detoxification processes for EBZ and CHP.

De novo transcriptomes of six calanoid copepods (Crustacea): a resource for the discovery of novel genes

This study presents eight new high-quality de novo transcriptomes from six co-occurring species of calanoid copepods, the first published for Neocalanus plumchrus, N. cristatus, Eucalanus bungii and Metridia pacifica and additional ones for N. flemingeri and Calanus marshallae. They are ecologically-important members of sub-arctic North Pacific marine zooplankton communities. 'Omics data for this diverse and numerous taxonomic group are sparse and difficult to obtain. Total RNA from single individuals was used to construct gene libraries that were sequenced on an Illumina Next-Seq platform. Quality filtered reads were assembled with Trinity software and validated using multiple criteria. The study's primary purpose is to provide a resource for gene expression studies. The integrated database can be used for quantitative inter- and intra-species comparisons of gene expression patterns across biological processes. An example of an additional use is provided for discovering novel and evolutionarily-significant proteins within the Calanoida. A workflow was designed to find and characterize unannotated transcripts with homologies across de novo assemblies that have also been shown to be eco-responsive.

Genome-wide analyses of Glutathione S-transferase gene family and expression profiling under deltamethrin exposure in non-biting midge Propsilocerus akamusi

Glutathione S-transferases (GSTs) are major enzymes in detoxification phase II, and have been functioned in resistance to various insecticides or oxidative stress. Herein, we selected the non-biting midge, Propsilocerus akamusi, widespread in Asian aquatic ecosystems, to uncover the gene location, structure, and phylogenetics relationship of GSTs at genome scale first time. Thirty-three cytosolic and four microsomal GST genes were identified and located on the four chromosomes. The cytosolic GSTs involved in the eight subclasses and five GST genes were unclassified. The expansion of GST genes in P. akamusi experienced duplication events on the delta, theta, xi, iota, and unclassified subclasses. The RNA-Seq analyses and RT-qPCR validation showed that the expression of PaGSTt2 gene is significantly elevated, with deltamethrin concentration increasing. The tertiary structure of PaGSTt2 enzyme was reconstructed, which was different from the other theta gene in the active site. In addition, the GST genes of six chironomids were first described based on the assembled genomes to explore the difference of those in the adaptation to kinds of environments. The GST frame for P. akmusi and its expression profiles provide valuable resources to understand their role in insecticide resistance of this species, as well as those of other biting midges.

Alpha-class glutathione S-transferases involved in the detoxification of aflatoxin B1 in ducklings

The objective of this study was to identify the key glutathione S-transferase (GST) isozymes involved in the detoxification of Aflatoxin B₁ (AFB₁) in ducks' primary hepatocytes. The full-length cDNA encoding the 10 GST isozymes (GST, GST3, GSTM3, MGST1, MGST2, MGST3, GSTK1, GSTT1, GSTO1 and GSTZ1) were isolated/synthesized from ducks' liver and cloned into the pcDNA3.1(+) vector. The results showed that pcDNA3.1(+)-GSTs plasmids were successfully transferred into the ducks' primary hepatocytes and the mRNA of the 10 GST isozymes were overexpressed by 1.9-3274.7 times. Compared to the control, 75 μg/L (IC₃₀) or 150 μg/L (IC₅₀) AFB₁ treatment reduced the cell viability by 30.0-50.0% and increased the LDH activity by 19.8-58.2% in the ducks' primary hepatocytes. Notably, the AFB₁-induced changes in cell viability and LDH activity were mitigated by overexpression of GST and GST3. Compared to the cells treated with AFB₁, exo-AFB₁-8,9-epoxide (AFBO)-GSH, as the major detoxified product of AFB₁, was increased in the cells overexpression of GST and GST3. Moreover, the sequences, phylogenetic and domain analysis revealed that the GST and GST3 were orthologous to Meleagris gallopavo GSTA3 and GSTA4. In conclusion, this study found that the ducks' GST and GST3 were orthologous to Meleagris gallopavo GSTA3 and GSTA4, which were involved in the detoxification of AFB₁ in ducks' primary hepatocytes.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families - genes directly mediating interactions with plant chemical defenses - underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila .

Chronic exposure to the star polycation (SPc) nanocarrier in the larval stage adversely impairs life history traits in Drosophila melanogaster

BACKGROUND

Nanomaterials are widely used as pesticide adjuvants to increase pesticide efficiency and minimize environmental pollution. But it is increasingly recognized that nanocarrier is a double-edged sword, as nanoparticles are emerging as new environmental pollutants. This study aimed to determine the biotoxicity of a widely applied star polycation (SPc) nanocarrier using Drosophila melanogaster, the fruit fly, as an in vivo model.

RESULTS

The lethal concentration 50 (LC₅₀) value of SPc was identified as 2.14 g/L toward third-instar larvae and 26.33 g/L for adults. Chronic exposure to a sub lethal concentration of SPc (1 g/L) in the larval stage showed long-lasting adverse effects on key life history traits. Exposure to SPc at larval stage adversely impacted the lifespan, fertility, climbing ability as well as stresses resistance of emerged adults. RNA-sequencing analysis found that SPc resulted in aberrant expression of genes involved in metabolism, innate immunity, stress response and hormone production in the larvae. Orally administrated SPc nanoparticles were mainly accumulated in intestine cells, while systemic responses were observed.

CONCLUSIONS

These findings indicate that SPc nanoparticles are hazardous to fruit flies at multiple levels, which could help us to develop guidelines for further large-scale application.

A near-chromosome level genome assembly of the European hoverfly, Sphaerophoria rueppellii (Diptera: Syrphidae), provides comparative insights into insecticide resistance-related gene family evolution

BACKGROUND

Sphaerophoria rueppellii, a European species of hoverfly, is a highly effective beneficial predator of hemipteran crop pests including aphids, thrips and coleopteran/lepidopteran larvae in integrated pest management (IPM) programmes. It is also a key pollinator of a wide variety of important agricultural crops. No genomic information is currently available for S. rueppellii. Without genomic information for such beneficial predator species, we are unable to perform comparative analyses of insecticide target-sites and genes encoding metabolic enzymes potentially responsible for insecticide resistance, between crop pests and their predators. These metabolic mechanisms include several gene families - cytochrome P450 monooxygenases (P450s), ATP binding cassette transporters (ABCs), glutathione-S-transferases (GSTs), UDP-glycosyltransferases (UGTs) and carboxyl/choline esterases (CCEs).

METHODS AND FINDINGS

In this study, a high-quality near-chromosome level de novo genome assembly (as well as a mitochondrial genome assembly) for S. rueppellii has been generated using a hybrid approach with PacBio long-read and Illumina short-read data, followed by super scaffolding using Hi-C data. The final assembly achieved a scaffold N50 of 87Mb, a total genome size of 537.6Mb and a level of completeness of 96% using a set of 1,658 core insect genes present as full-length genes. The assembly was annotated with 14,249 protein-coding genes. Comparative analysis revealed gene expansions of CYP6Zx P450s, epsilon-class GSTs, dietary CCEs and multiple UGT families (UGT37/302/308/430/431). Conversely, ABCs, delta-class GSTs and non-CYP6Zx P450s showed limited expansion. Differences were seen in the distributions of resistance-associated gene families across subfamilies between S. rueppellii and some hemipteran crop pests. Additionally, S. rueppellii had larger numbers of detoxification genes than other pollinator species.

CONCLUSION AND SIGNIFICANCE

This assembly is the first published genome for a predatory member of the Syrphidae family and will serve as a useful resource for further research into selectivity and potential tolerance of insecticides by beneficial predators. Furthermore, the expansion of some gene families often linked to insecticide resistance and selectivity may be an indicator of the capacity of this predator to detoxify IPM selective insecticides. These findings could be exploited by targeted insecticide screens and functional studies to increase effectiveness of IPM strategies, which aim to increase crop yields by sustainably and effectively controlling pests without impacting beneficial predator populations.

The transcription factor Xrp1 is required for PERK-mediated antioxidant gene induction in Drosophila

PERK is an endoplasmic reticulum (ER) transmembrane sensor that phosphorylates eIF2α to initiate the Unfolded Protein Response (UPR). eIF2α phosphorylation promotes stress-responsive gene expression most notably through the transcription factor ATF4 that contains a regulatory 5’ leader. Possible PERK effectors other than ATF4 remain poorly understood. Here, we report that the bZIP transcription factor Xrp1 is required for ATF4-independent PERK signaling. Cell-type-specific gene expression profiling in Drosophila indicated that delta-family glutathione-S-transferases (gstD) are prominently induced by the UPR-activating transgene Rh1^G69D. Perk was necessary and sufficient for such gstD induction, but ATF4 was not required. Instead, Perk and other regulators of eIF2α phosphorylation regulated Xrp1 protein levels to induce gstDs. The Xrp1 5’ leader has a conserved upstream Open Reading Frame (uORF) analogous to those that regulate ATF4 translation. The gstD-GFP reporter induction required putative Xrp1 binding sites. These results indicate that antioxidant genes are highly induced by a previously unrecognized UPR signaling axis consisting of PERK and Xrp1.

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding.

Transcriptome analysis and comparison reveal divergence between the Mediterranean and the greenhouse whiteflies

Both the Mediterranean (MED) species of the Bemisia tabaci whitefly complex and the greenhouse whitefly (Trialeurodes vaporariorum, TV) are important agricultural pests. The two species of whiteflies differ in many aspects such as morphology, geographical distribution, host plant range, plant virus transmission, and resistance to insecticides. However, the molecular basis underlying their differences remains largely unknown. In this study, we analyzed the genetic divergences between the transcriptomes of MED and TV. In total, 2,944 pairs of orthologous genes were identified. The average identity of amino acid sequences between the two species is 93.6%. The average nonsynonymous (Ka) and synonymous (Ks) substitution rates and the ratio of Ka/Ks of the orthologous genes are 0.0389, 2.23 and 0.0204, respectively. The low average Ka/Ks ratio indicates that orthologous genes tend to be under strong purified selection. The most divergent gene classes are related to the metabolisms of xenobiotics, cofactors, vitamins and amino acids, and this divergence may underlie the different biological characteristics between the two species of whiteflies. Genes of differential expression between the two species are enriched in carbohydrate metabolism and regulation of autophagy. These findings provide molecular clues to uncover the biological and molecular differences between the two species of whiteflies.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Behavioral and Transcriptional Response to Selection for Olfactory Behavior in Drosophila

The detection, discrimination, and behavioral responses to chemical cues in the environment can have marked effects on organismal survival and reproduction, eliciting attractive or aversive behavior. To gain insight into mechanisms mediating this hedonic valence, we applied thirty generations of divergent artificial selection for Drosophila melanogaster olfactory behavior. We independently selected for positive and negative behavioral responses to two ecologically relevant chemical compounds: 2,3-butanedione and cyclohexanone. We also tested the correlated responses to selection by testing behavioral responses to other odorants and life history traits. Measurements of behavioral responses of the selected lines and unselected controls to additional odorants showed that the mechanisms underlying responses to these odorants are, in some cases, differentially affected by selection regime and generalization of the response to other odorants was only detected in the 2,3-butanedione selection lines. Food consumption and lifespan varied with selection regime and, at times, sex. An analysis of gene expression of both selection regimes identified multiple differentially expressed genes. New genes and genes previously identified in mediating olfactory behavior were identified. In particular, we found functional enrichment of several gene ontology terms, including cell-cell adhesion and sulfur compound metabolic process, the latter including genes belonging to the glutathione S-transferase family. These findings highlight a potential role for glutathione S-transferases in the evolution of hedonic valence to ecologically relevant volatile compounds and set the stage for a detailed investigation into mechanisms by which these genes mediate attraction and aversion.

Molecular Evolution of the Glutathione S-Transferase Family in the Bemisia tabaci Species Complex

The glutathione S-transferase (GST) family plays an important role in the adaptation of herbivorous insects to new host plants and other environmental constrains. The family codes for enzymes that neutralize reactive oxygen species and phytotoxins through the conjugation of reduced glutathione. Here, we studied the molecular evolution of the GST family in Bemisia tabaci, a complex of >35 sibling species, differing in their geographic and host ranges. We tested if some enzymes evolved different functionality, by comparing their sequences in six species, representing five of the six major genetic clades in the complex. Comparisons of the nonsynonymous to synonymous substitution ratios detected positive selection events in 11 codons of 5 cytosolic GSTs. Ten of them are located in the periphery of the GST dimer, suggesting a putative involvement in interactions with other proteins. Modeling the tertiary structure of orthologous enzymes, identified additional 19 mutations in 9 GSTs, likely affecting the enzymes' functionality. Most of the mutation events were found in the environmentally responsive classes Delta and Sigma, indicating a slightly different delta/sigma tool box in each species. At a broader genomic perspective, our analyses indicated a significant expansion of the Delta GST class in B. tabaci and a general association between the diet breadth of hemipteran species and their total number of GST genes. We raise the possibility that at least some of the identified changes improve the fitness of the B. tabaci species carrying them, leading to their better adaptation to specific environments.

Contraction of the ROS Scavenging Enzyme Glutathione S-Transferase Gene Family in Cetaceans

Cetaceans are a group of marine mammals whose ancestors were adaptated for life on land. Life in an aquatic environment poses many challenges for air-breathing mammals. Diving marine mammals have adapted to rapid reoxygenation and reactive oxygen species (ROS)-mediated reperfusion injury. Here, we considered the evolution of the glutathione transferase (GST) gene family which has important roles in the detoxification of endogenously-derived ROS and environmental pollutants. We characterized the cytosolic GST gene family in 21 mammalian species; cetaceans, sirenians, pinnipeds, and their terrestrial relatives. All seven GST classes were identified, showing that GSTs are ubiquitous in mammals. Some GST genes are the product of lineage-specific duplications and losses, in line with a birth-and-death evolutionary model. We detected sites with signatures of positive selection that possibly influence GST structure and function, suggesting that adaptive evolution of GST genes is important for defending mammals from various types of noxious environmental compounds. We also found evidence for loss of alpha and mu GST subclass genes in cetacean lineages. Notably, cetaceans have retained a homolog of at least one of the genes GSTA1, GSTA4, and GSTM1; GSTs that are present in both the cytosol and mitochondria. The observed variation in number and selection pressure on GST genes suggest that the gene family structure is dynamic within cetaceans.

Host use and host shifts in Drosophila

Over a thousand Drosophila species have radiated onto a wide range of feeding and breeding sites. These radiations involve adaptations for locating, accepting, and growing in hosts with highly differing characteristics. In a number of species, owing to the availability of sequenced genomes, particular steps in host specialization and genes that control them, are being identified. Many cases of specialization involve the ability to detoxify some component of the host. Examples include Drosophila sechellia and the octanoic acid in Morinda citrifolia, alpha-amanitin in mycophagous drosophilids, and the alkaloids in cactophilic species. Owing to the known ecologies of many species for which genomes exist, the Drosophila model system provides an unprecedented opportunity to simultaneously examine the genes underlying HOST LOCATION, HOST ACCEPTANCE and HOST USE, the types of selection acting upon them and any coevolutionary interactions among the genes underlying these steps.

Rapid birth-death evolution and positive selection in detoxification-type glutathione S-transferases in mammals

Glutathione S-Transferases (GSTs) are phase II detoxification enzymes that may have evolved in response to changes of environmental substrates. GST genes formed a multigene family and in mammals, there are six classes known as Alpha, Mu, Omega, Pi, Theta, and Zeta. Recent studies in phase I detoxification system specifically the cytochrome P450s provided a general explanation on why genes from a common origin such as those in a multigene family have both phylogenetically stable and unstable genes. Genes that participate in core functions of organisms such as development and physiology are stable whereas genes that play a role in detoxification are unstable and evolve in a process known as birth-death evolution, which is characterised by frequent gene gains and losses. The generality of the birth-death model at explaining the evolution of detoxification enzymes beyond the phase I enzyme has not been comprehensively explored. This work utilized 383 Gst genes and 300 pseudogenes across 22 mammalian species to study gene gains and losses. GSTs vary greatly in their phylogenetic stability despite their overall sequence similarity. Stable Gst genes from Omega and Zeta classes do not show fluctuation in gene numbers from human to opossum. These genes play a role in biosynthesis related functions. Unstable genes that include Alpha, Mu, Pi and Theta undergo frequent gene gain and loss in a process known as birth-death evolution. Gene members of these four classes are well known for their roles in detoxification. Our positive selection screen identified five positively selected sites in mouse GSTA3. Previous studies showed two of these sites (108H and 208E) were biochemically tested as important residues that conferred catalytic activity against the toxic aflatoxin B1-8,9-epoxide. The functional significance against aflatoxin of the remaining three positively selected sites warrant further investigation.

Complex evolution of the GSTM gene family involves sharing of GSTM1 deletion polymorphism in humans and chimpanzees

Background

The common deletion of the glutathione S-transferase Mu 1 (GSTM1) gene in humans has been shown to be involved in xenobiotic metabolism and associated with bladder cancer. However, the evolution of this deletion has not been investigated.

Results

In this study, we conducted comparative analyses of primate genomes. We demonstrated that the GSTM gene family has evolved through multiple structural variations, involving gene duplications, losses, large inversions and gene conversions. We further showed experimentally that the GSTM1 was polymorphically deleted in both humans and also in chimpanzees, through independent deletion events. To generalize our results, we searched for genic deletions that are polymorphic in both humans and chimpanzees. Consequently, we found only two such deletions among the thousands that we have searched, one of them being the GSTM1 deletion and the other surprisingly being another metabolizing gene, the UGT2B17.

Conclusions

Overall, our results support the emerging notion that metabolizing gene families, such as the GSTM, NAT, UGT and CYP, have been evolving rapidly through gene duplication and deletion events in primates, leading to complex structural variation within and among species with unknown evolutionary consequences.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4676-z) contains supplementary material, which is available to authorized users.

Dual Roles of Glutathione in Ecdysone Biosynthesis and Antioxidant Function During Larval Development in Drosophila

Ecdysteroids, including the biologically active hormone 20-hydroxyecdysone (20E), play essential roles in controlling many developmental and physiological events in insects. Ecdysteroid biosynthesis is achieved by a series of specialized enzymes encoded by the Halloween genes. Recently, a new class of Halloween gene, noppera-bo (nobo), encoding a glutathione S-transferase (GST) in dipteran and lepidopteran species, has been identified and characterized. GSTs are well known to conjugate substrates with the reduced form of glutathione (GSH), a bioactive tripeptide composed of glutamate, cysteine, and glycine. We hypothesized that GSH itself is required for ecdysteroid biosynthesis. However, the role of GSH in steroid hormone biosynthesis has not been examined in any organisms. Here, we report phenotypic analysis of a complete loss-of-function mutant in the γ-glutamylcysteine synthetase catalytic subunit (Gclc) gene in the fruit fly Drosophila melanogaster. Gclc encodes the evolutionarily conserved catalytic component of the enzyme that conjugates glutamate and cysteine in the GSH biosynthesis pathway. Complete Gclc loss-of-function leads to drastic GSH deficiency in the larval body fluid. Gclc mutant animals show a larval-arrest phenotype. Ecdysteroid titer in Gclc mutant larvae decreases, and the larval-arrest phenotype is rescued by oral administration of 20E or cholesterol. Moreover, Gclc mutant animals exhibit abnormal lipid deposition in the prothoracic gland, a steroidogenic organ during larval development. All of these phenotypes are reminiscent to nobo loss-of-function animals. On the other hand, Gclc mutant larvae also exhibit a significant reduction in antioxidant capacity. Consistent with this phenotype, Gclc mutant larvae are more sensitive to oxidative stress response as compared to wild-type. Nevertheless, the ecdysteroid biosynthesis defect in Gclc mutant animals is not associated with loss of antioxidant function. Our data raise the unexpected hypothesis that a primary role of GSH in early D. melanogaster larval development is ecdysteroid biosynthesis, independent from the antioxidant role of GSH.

Senecio brasiliensis impairs eclosion rate and induces apoptotic cell death in larvae of Drosophila melanogaster

Senecio brasilienis (Spreng) Less., is a species native from Brazil, popularly known as "Maria mole", and known to induce hepatotoxicity due to its high content of Pyrrolizidine alkaloids. Despite its toxicity, this plant is widely used in Brazilian folk medicine. Considering the antagonizing effects described for S. brasiliensis, we describe here molecular markers involved in the toxicity of hydroalcoholic extract from leaves of S. brasiliensis (HESB) in Drosophila melanogaster. Phytochemical analysis of HESB revealed the presence of phenolic acids and flavonoids. A significant antioxidant potential against ABTS⁺ and DPPH radical was found in parallel. Ingestion of extract did not alter the survival and locomotor activity of adult flies. However when ingested along the larval developmental phase, the eclosion rate of flies was interrupted at higher concentration of extract. To comprehend this phenomenon several analysis were conducted in larvae. HESB stimulated activity of antioxidant enzymes SOD and GST, and increased GSH/GSSG ratio and ROS production. Additionally, HESB caused a significant decrease of cell viability. The mRNA expression of Nrf2, TrxR, CAT, Drice and Dilp6 were also significantly up-regulated. HESB caused significant decrease on the phosphorylation of MAPKs and AKT. In parallel, PARP cleavage and caspases 3/7 activity were stimulated. In addition, glucose, glycogen and triglycerides levels were decreased. Taken together our study depicts a disruption in the eclosion of D. melanogaster possibly attributed to the inhibition of kinases implied in developmental process, energetic demand and induction of apoptotic cell death process.

Arabidopsis glucosinolates trigger a contrasting transcriptomic response in a generalist and a specialist herbivore

Phytophagous insects have to deal with toxic defense compounds from their host plants. Although it is known that insects have evolved genes and mechanisms to detoxify plant allochemicals, how specialist and generalist precisely respond to specific secondary metabolites at the molecular level is less understood. Here we studied the larval performance and transcriptome of the generalist moth Heliothis virescens and the specialist butterfly Pieris brassicae feeding on Arabidopsis thaliana genotypes with different glucosinolate (GS) levels. H. virescens larvae gained significantly more weight on the GS-deficient mutant quadGS compared to wild-type (Col-0) plants. On the contrary, P. brassicae was unaffected by the presence of GS and performed equally well on both genotypes. Strikingly, there was a considerable differential gene expression in H. virescens larvae feeding on Col-0 compared to quadGS. In contrast, compared to H. virescens, P. brassicae displayed a much-reduced transcriptional activation when fed on both plant genotypes. Transcripts coding for putative detoxification enzymes were significantly upregulated in H. virescens, along with digestive enzymes and transposable elements. These data provide an unprecedented view on transcriptional changes that are specifically activated by GS and illustrate differential molecular responses that are linked to adaptation to diet in lepidopteran herbivores.

Transcriptome modulation during host shift is driven by secondary metabolites in desert Drosophila

High-throughput transcriptome studies are breaking new ground to investigate the responses that organisms deploy in alternative environments. Nevertheless, much remains to be understood about the genetic basis of host plant adaptation. Here, we investigate genome-wide expression in the fly Drosophila buzzatii raised in different conditions. This species uses decaying tissues of cactus of the genus Opuntia as primary rearing substrate and secondarily, the necrotic tissues of the columnar cactus Trichocereus terscheckii. The latter constitutes a harmful host, rich in mescaline and other related phenylethylamine alkaloids. We assessed the transcriptomic responses of larvae reared in Opuntia sulphurea and T. terscheckii, with and without the addition of alkaloids extracted from the latter. Whole-genome expression profiles were massively modulated by the rearing environment, mainly by the presence of T. terscheckii alkaloids. Differentially expressed genes were mainly related to detoxification, oxidation-reduction and stress response; however, we also found genes involved in development and neurobiological processes. In conclusion, our study contributes new data onto the role of transcriptional plasticity in response to alternative rearing environments.

The Biochemical Adaptations of Spotted Wing Drosophila (Diptera: Drosophilidae) to Fresh Fruits Reduced Fructose Concentrations and Glutathione-S Transferase Activities

Spotted wing drosophila, Drosophila suzukii Matsumura, is an invasive and economically damaging pest in Europe and North America. The females have a serrated ovipositor that enables them to infest almost all ripening small fruits. To understand the physiological and metabolic basis of spotted wing drosophila food preferences for healthy ripening fruits, we investigated the biological and biochemical characteristics of spotted wing drosophila and compared them with those of Drosophila melanogaster Meigen. We found that the susceptibility to oxidative stressors was significantly increased in spotted wing drosophila compared with those of D. melanogaster. In addition, we found that spotted wing drosophila had significantly reduced glutathione-S transferase (GST) activity and gene numbers. Furthermore, fructose concentrations found in spotted wing drosophila were significantly lower than those of D. melanogaster. Our data strongly suggest that the altered food preferences of spotted wing drosophila may stem from evolutionary adaptations to fresh foods accompanied by alterations in carbohydrate metabolism and GST activities.

Status of pyrethroid resistance and mechanisms in Brazilian populations of Tuta absoluta

The tomato leafminer, Tuta absoluta, is a major pest of tomato crops worldwide. This study surveyed the resistance of T. absoluta populations from four regions in Brazil to pyrethroid insecticides, the frequencies of L1014F, T929I and M918T Na channel mutations, and the role of detoxification metabolism in the resistance. Resistance ratios varied from 1- to 11-times among populations and insecticides, but control failure likelihood assays showed that all pyrethroids assessed exhibited no efficacy at all (and thus, 98-100% control failure likelihood) against all T. absoluta populations. The activity of glutathione S-transferase and cytochrome P450-mediated N-demethylation in biochemical assays was significantly correlated with the level of resistance to deltamethrin and permethrin suggesting that these enzymes may play a role in resistance. TaqMan assays were used to screen for the presence of knockdown resistance (kdr) mutations and revealed that the L1014F kdr mutation was fixed in all populations and associated with two super-kdr mutations, M918T and particularly T929I, at high frequency. Altogether, results suggest that control failures are because of mutations in the domain II of the sodium channel, as a prevailing mechanism of resistance to pyrethroids in populations of T. absoluta in Brazil. But, enhanced cytochrome P450-dependent monooxygenases and GST activities also play an important role in the resistance of some populations, which reinforce that pyrethroids must not be used overall to control T. absoluta.

The secret lives of Drosophila flies

Flies of the genus Drosophila, and particularly those of the species Drosophila melanogaster, are best known as laboratory organisms. As with all model organisms, they were domesticated for empirical studies, but they also continue to exist as wild populations.

Decades of research on these flies in the laboratory have produced astounding and important insights into basic biological processes, but we have only scratched the surface of what they have to offer as research organisms. An outstanding challenge now is to build on this knowledge and explore how natural history has shaped D. melanogaster in order to advance our understanding of biology more generally.

DOI:http://dx.doi.org/10.7554/eLife.06793.001

Characterization and expression profiling of glutathione S-transferases in the diamondback moth, Plutella xylostella (L.)

Background

Glutathione S-transferases (GSTs) are multifunctional detoxification enzymes that play important roles in insects. The completion of several insect genome projects has enabled the identification and characterization of GST genes over recent years. This study presents a genome-wide investigation of the diamondback moth (DBM), Plutella xylostella, a species in which the GSTs are of special importance because this pest is highly resistant to many insecticides.

Results

A total of 22 putative cytosolic GSTs were identified from a published P. xylostella genome and grouped into 6 subclasses (with two unclassified). Delta, Epsilon and Omega GSTs were numerically superior with 5 genes for each of the subclasses. The resulting phylogenetic tree showed that the P. xylostella GSTs were all clustered into Lepidoptera-specific branches. Intron sites and phases as well as GSH binding sites were strongly conserved within each of the subclasses in the GSTs of P. xylostella. Transcriptome-, RNA-seq- and qRT-PCR-based analyses showed that the GST genes were developmental stage- and strain-specifically expressed. Most of the highly expressed genes in insecticide resistant strains were also predominantly expressed in the Malpighian tubules, midgut or epidermis.

Conclusions

To date, this is the most comprehensive study on genome-wide identification, characterization and expression profiling of the GST family in P. xylostella. The diversified features and expression patterns of the GSTs are inferred to be associated with the capacity of this species to develop resistance to a wide range of pesticides and biological toxins. Our findings provide a base for functional research on specific GST genes, a better understanding of the evolution of insecticide resistance, and strategies for more sustainable management of the pest.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1343-5) contains supplementary material, which is available to authorized users.

Aniline exposure associated with up-regulated transcriptional responses of three glutathione S-transferase Delta genes in Drosophila melanogaster

Complex transcriptional profile of glutathione S-transferase Delta cluster genes occurred in the developmental process of the fruit fly Drosophila melanogaster. The purpose of this project was to quantify the expression levels of Gst Delta class genes altered by aniline exposure and to understand the relationship between aniline dosages and the variation of Gst Delta genes expressed in D. melanogaster. Using RT-PCR expression assays, the expression patterns of the transcript mRNAs of the glutathione S-transferase Delta genes were revealed and their expression levels were measured at eggs, larvae, pupae and adults. The adult stage was selected for further dose-response assays. After analysis, the results indicated that three Gst Delta genes (Gst D2, Gst D5 and Gst D6) were found to show a peak of up-regulated transcriptional response at 6-8h of exposure of aniline. Furthermore, the dose-response relationship of their induction levels within the dose regiments (from 1.2 to 2.0 μl/tube) had been measured. The expression patterns and annotations of these genes were discussed in the context.

Glutathione transferases in the bioactivation of azathioprine

The prodrug azathioprine is primarily used for maintaining remission in inflammatory bowel disease, but approximately 30% of the patients suffer adverse side effects. The prodrug is activated by glutathione conjugation and release of 6-mercaptopurine, a reaction most efficiently catalyzed by glutathione transferase (GST) A2-2. Among five genotypes of GST A2-2, the variant A2*E has threefold-fourfold higher catalytic efficiency with azathioprine, suggesting that the expression of A2*E could boost 6-mercaptopurine release and adverse side effects in treated patients. Structure-activity studies of the GST A2-2 variants and homologous alpha class GSTs were made to delineate the determinants of high catalytic efficiency compared to other alpha class GSTs. Engineered chimeras identified GST peptide segments of importance, and replacing the corresponding regions in low-activity GSTs by these short segments produced chimeras with higher azathioprine activity. By contrast, H-site mutagenesis led to decreased azathioprine activity when active-site positions 208 and 213 in these favored segments were mutagenized. Alternative substitutions indicated that hydrophobic residues were favored. A pertinent question is whether variant A2*E represents the highest azathioprine activity achievable within the GST structural framework. This issue was addressed by mutagenesis of H-site residues assumed to interact with the substrate based on molecular modeling. The mutants with notably enhanced activities had small or polar residues in the mutated positions. The most active mutant L107G/L108D/F222H displayed a 70-fold enhanced catalytic efficiency with azathioprine. The determination of its structure by X-ray crystallography showed an expanded H-site, suggesting improved accommodation of the transition state for catalysis.

Micro-plasticity of genomes as illustrated by the evolution of glutathione transferases in 12 Drosophila species

Glutathione transferases (GST) are an ancient superfamily comprising a large number of paralogous proteins in a single organism. This multiplicity of GSTs has allowed the copies to diverge for neofunctionalization with proposed roles ranging from detoxication and oxidative stress response to involvement in signal transduction cascades. We performed a comparative genomic analysis using FlyBase annotations and Drosophila melanogaster GST sequences as templates to further annotate the GST orthologs in the 12 Drosophila sequenced genomes. We found that GST genes in the Drosophila subgenera have undergone repeated local duplications followed by transposition, inversion, and micro-rearrangements of these copies. The colinearity and orientations of the orthologous GST genes appear to be unique in many of the species which suggests that genomic rearrangement events have occurred multiple times during speciation. The high micro-plasticity of the genomes appears to have a functional contribution utilized for evolution of this gene family.

Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogaster

The metabolism of volatile signal molecules by odorant degrading enzymes (ODEs) is crucial to the ongoing sensitivity and specificity of chemoreception in various insects, and a few specific esterases, cytochrome P450s, glutathione S-transferases (GSTs) and UDP-glycosyltransferases (UGTs) have previously been implicated in this process. Significant progress has been made in characterizing ODEs in Lepidoptera but very little is known about them in Diptera, including in Drosophila melanogaster, a major insect model. We have therefore carried out a transcriptomic analysis of the antennae of D. melanogaster in order to identify candidate ODEs. Virgin male and female and mated female antennal transcriptomes were determined by RNAseq. As with the Lepidoptera, we found that many esterases, cytochrome P450 enzymes, GSTs and UGTs are expressed in D. melanogaster antennae. As olfactory genes generally show selective expression in the antennae, a comparison to previously published transcriptomes for other tissues has been performed, showing preferential expression in the antennae for one esterase, JHEdup, one cytochrome P450, CYP308a1, and one GST, GSTE4. These largely uncharacterized enzymes are now prime candidates for ODE functions. JHEdup was expressed heterologously and found to have high catalytic activity against a chemically diverse group of known ester odorants for this species. This is a finding consistent with an ODE although it might suggest a general role in clearing several odorants rather than a specific role in clearing a particular odorant. Our findings do not preclude the possibility of odorant degrading functions for other antennally expressed esterases, P450s, GSTs and UGTs but, if so, they suggest that these enzymes also have additional functions in other tissues.

Whole-genome expression analysis in the third instar larval midgut of Drosophila melanogaster

Survival of insects on a substrate containing toxic substances such as plant secondary metabolites or insecticides is dependent on the metabolism or excretion of those xenobiotics. The primary sites of xenobiotic metabolism are the midgut, Malpighian tubules, and fat body. In general, gene expression in these organs is reported for the entire tissue by online databases, but several studies have shown that gene expression within the midgut is compartmentalized. Here, RNA sequencing is used to investigate whole-genome expression in subsections of third instar larval midguts of Drosophila melanogaster. The data support functional diversification in subsections of the midgut. Analysis of the expression of gene families that are implicated in the metabolism of xenobiotics suggests that metabolism may not be uniform along the midgut. These data provide a starting point for investigating gene expression and xenobiotic metabolism and other functions of the larval midgut.

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.

Evolution in an ancient detoxification pathway is coupled with a transition to herbivory in the drosophilidae

Chemically defended plant tissues present formidable barriers to herbivores. Although mechanisms to resist plant defenses have been identified in ancient herbivorous lineages, adaptations to overcome plant defenses during transitions to herbivory remain relatively unexplored. The fly genus Scaptomyza is nested within the genus Drosophila and includes species that feed on the living tissue of mustard plants (Brassicaceae), yet this lineage is derived from microbe-feeding ancestors. We found that mustard-feeding Scaptomyza species and microbe-feeding Drosophila melanogaster detoxify mustard oils, the primary chemical defenses in the Brassicaceae, using the widely conserved mercapturic acid pathway. This detoxification strategy differs from other specialist herbivores of mustard plants, which possess derived mechanisms to obviate mustard oil formation. To investigate whether mustard feeding is coupled with evolution in the mercapturic acid pathway, we profiled functional and molecular evolutionary changes in the enzyme glutathione S-transferase D1 (GSTD1), which catalyzes the first step of the mercapturic acid pathway and is induced by mustard defense products in Scaptomyza. GSTD1 acquired elevated activity against mustard oils in one mustard-feeding Scaptomyza species in which GstD1 was duplicated. Structural analysis and mutagenesis revealed that substitutions at conserved residues within and near the substrate-binding cleft account for most of this increase in activity against mustard oils. Functional evolution of GSTD1 was coupled with signatures of episodic positive selection in GstD1 after the evolution of herbivory. Overall, we found that preexisting functions of generalized detoxification systems, and their refinement by natural selection, could play a central role in the evolution of herbivory.

Analysis of insecticide resistance-related genes of the Carmine spider mite Tetranychus cinnabarinus based on a de novo assembled transcriptome

The carmine spider mite (CSM), Tetranychus cinnabarinus, is an important pest mite in agriculture, because it can develop insecticide resistance easily. To gain valuable gene information and molecular basis for the future insecticide resistance study of CSM, the first transcriptome analysis of CSM was conducted. A total of 45,016 contigs and 25,519 unigenes were generated from the de novo transcriptome assembly, and 15,167 unigenes were annotated via BLAST querying against current databases, including nr, SwissProt, the Clusters of Orthologous Groups (COGs), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). Aligning the transcript to Tetranychus urticae genome, the 19255 (75.45%) of the transcripts had significant (e-value <10-5) matches to T. urticae DNA genome, 19111 sequences matched to T. urticae proteome with an average protein length coverage of 42.55%. Core Eukaryotic Genes Mapping Approach (CEGMA) analysis identified 435 core eukaryotic genes (CEGs) in the CSM dataset corresponding to 95% coverage. Ten gene categories that relate to insecticide resistance in arthropod were generated from CSM transcriptome, including 53 P450-, 22 GSTs-, 23 CarEs-, 1 AChE-, 7 GluCls-, 9 nAChRs-, 8 GABA receptor-, 1 sodium channel-, 6 ATPase- and 12 Cyt b genes. We developed significant molecular resources for T. cinnabarinus putatively involved in insecticide resistance. The transcriptome assembly analysis will significantly facilitate our study on the mechanism of adapting environmental stress (including insecticide) in CSM at the molecular level, and will be very important for developing new control strategies against this pest mite.

The molecular evolution of cytochrome P450 genes within and between drosophila species

We map 114 gene gains and 74 gene losses in the P450 gene family across the phylogeny of 12 Drosophila species by examining the congruence of gene trees and species trees. Although the number of P450 genes varies from 74 to 94 in the species examined, we infer that there were at least 77 P450 genes in the ancestral Drosophila genome. One of the most striking observations in the data set is the elevated loss of P450 genes in the Drosophila sechellia lineage. The gain and loss events are not evenly distributed among the P450 genes-with 30 genes showing no gene gains or losses whereas others show as many as 20 copy number changes among the species examined. The P450 gene clades showing the fewest number of gene gain and loss events tend to be those evolving with the most purifying selection acting on the protein sequences, although there are exceptions, such as the rapid rate of amino acid replacement observed in the single copy phantom (Cyp306a1) gene. Within D. melanogaster, we observe gene copy number polymorphism in ten P450 genes including multiple cases of interparalog chimeras. Nonallelic homologous recombination (NAHR) has been associated with deleterious mutations in humans, but here we provide a second possible example of an NAHR event in insect P450s being adaptive. Specifically, we find that a polymorphic Cyp12a4/Cyp12a5 chimera correlates with resistance to an insecticide. Although we observe such interparalog exchange in our within-species data sets, we have little evidence of it between species, raising the possibility that such events may occur more frequently than appreciated but are masked by subsequent sequence change.

Evolutionary changes in gene expression, coding sequence and copy-number at the Cyp6g1 locus contribute to resistance to multiple insecticides in Drosophila

Widespread use of insecticides has led to insecticide resistance in many populations of insects. In some populations, resistance has evolved to multiple pesticides. In Drosophila melanogaster, resistance to multiple classes of insecticide is due to the overexpression of a single cytochrome P450 gene, Cyp6g1. Overexpression of Cyp6g1 appears to have evolved in parallel in Drosophila simulans, a sibling species of D. melanogaster, where it is also associated with insecticide resistance. However, it is not known whether the ability of the CYP6G1 enzyme to provide resistance to multiple insecticides evolved recently in D. melanogaster or if this function is present in all Drosophila species. Here we show that duplication of the Cyp6g1 gene occurred at least four times during the evolution of different Drosophila species, and the ability of CYP6G1 to confer resistance to multiple insecticides exists in D. melanogaster and D. simulans but not in Drosophila willistoni or Drosophila virilis. In D. virilis, which has multiple copies of Cyp6g1, one copy confers resistance to DDT and another to nitenpyram, suggesting that the divergence of protein sequence between copies subsequent to the duplication affected the activity of the enzyme. All orthologs tested conferred resistance to one or more insecticides, suggesting that CYP6G1 had the capacity to provide resistance to anthropogenic chemicals before they existed. Finally, we show that expression of Cyp6g1 in the Malpighian tubules, which contributes to DDT resistance in D. melanogaster, is specific to the D. melanogaster–D. simulans lineage. Our results suggest that a combination of gene duplication, regulatory changes and protein coding changes has taken place at the Cyp6g1 locus during evolution and this locus may play a role in providing resistance to different environmental toxins in different Drosophila species.

Transcriptomic survey of the midgut of Anthonomus grandis (Coleoptera: Curculionidae)

Anthonomus grandis Boheman is a key pest in cotton crops in the New World. Its larval stage develops within the flower bud using it as food and as protection against its predators. This behavior limits the effectiveness of its control using conventional insecticide applications and biocontrol techniques. In spite of its importance, little is known about its genome sequence and, more important, its specific expression in key organs like the midgut. Total mRNA isolated from larval midguts was used for pyrosequencing. Sequence reads were assembled and annotated to generate a unigene data set. In total, 400,000 reads from A. grandis midgut with an average length of 237 bp were assembled and combined into 20,915 contigs. The assembled reads fell into 6,621 genes models. BlastX search using the NCBI-NR database showed that 3,006 unigenes had significant matches to known sequences. Gene Ontology (GO) mapping analysis evidenced that A. grandis is able to transcripts coding for proteins involved in catalytic processing of macromolecules that allows its adaptation to very different feeding source scenarios. Furthermore, transcripts encoding for proteins involved in detoxification mechanisms such as p450 genes, glutathione-S-transferase, and carboxylesterases are also expressed. This is the first report of a transcriptomic study in A. grandis and the largest set of sequence data reported for this species. These data are valuable resources to expand the knowledge of this insect group and could be used in the design of new control strategies based in molecular information.

Mining genes involved in insecticide resistance of Liposcelis bostrychophila Badonnel by transcriptome and expression profile analysis

BACKGROUND

Recent studies indicate that infestations of psocids pose a new risk for global food security. Among the psocids species, Liposcelis bostrychophila Badonnel has gained recognition in importance because of its parthenogenic reproduction, rapid adaptation, and increased worldwide distribution. To date, the molecular data available for L. bostrychophila is largely limited to genes identified through homology. Also, no transcriptome data relevant to psocids infection is available.

METHODOLOGY AND PRINCIPAL FINDINGS

In this study, we generated de novo assembly of L. bostrychophila transcriptome performed through the short read sequencing technology (Illumina). In a single run, we obtained more than 51 million sequencing reads that were assembled into 60,012 unigenes (mean size = 711 bp) by Trinity. The transcriptome sequences from different developmental stages of L. bostrychophila including egg, nymph and adult were annotated with non-redundant (Nr) protein database, gene ontology (GO), cluster of orthologous groups of proteins (COG), and KEGG orthology (KO). The analysis revealed three major enzyme families involved in insecticide metabolism as differentially expressed in the L. bostrychophila transcriptome. A total of 49 P450-, 31 GST- and 21 CES-specific genes representing the three enzyme families were identified. Besides, 16 transcripts were identified to contain target site sequences of resistance genes. Furthermore, we profiled gene expression patterns upon insecticide (malathion and deltamethrin) exposure using the tag-based digital gene expression (DGE) method.

CONCLUSION

The L. bostrychophila transcriptome and DGE data provide gene expression data that would further our understanding of molecular mechanisms in psocids. In particular, the findings of this investigation will facilitate identification of genes involved in insecticide resistance and designing of new compounds for control of psocids.

Variation of gene expression associated with colonisation of an anthropized environment: comparison between African and European populations of Drosophila simulans

The comparison of transcriptome profiles among populations is a powerful tool for investigating the role of gene expression change in adaptation to new environments. In this study, we use massively parallel sequencing of 3' cDNAs obtained from large samples of adult males, to compare a population of Drosophila simulans from a natural reserve within its ancestral range (eastern Africa) with a derived population collected in the strongly anthropized Rhône valley (France). The goal was to scan for adaptation linked to the invasion of new environments by the species. Among 15,090 genes retained for the analysis, 794 were found to be differentially expressed between the two populations. We observed an increase in expression of reproduction-related genes in eastern Africa, and an even stronger increase in expression of Cytochrome P450, Glutathione transferase and Glucuronosyl transferase genes in the derived population. These three gene families are involved in detoxification processes, which suggests that pesticides are a major environmental pressure for the species in this area. The survey of the Cyp6g1 upstream region revealed the insertion of a transposable element, Juan, in the regulatory sequence that is almost fixed in the Rhône Valley, but barely present in Mayotte. This shows that Cyp6g1 has undergone parallel evolution in derived populations of D. simulans as previously shown for D. melanogaster. The increasing amount of data produced by comparative population genomics and transcriptomics should permit the identification of additional genes associated with functional divergence among those differentially expressed.

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human-mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)-based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST -based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross-validated by four FST -based methods, and 8.9% (34 of 383) were cross-validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human-mediated environmental changes.

Identification, genomic organization, and oxidative stress response of a sigma class glutathione S-transferase gene (AccGSTS1) in the honey bee, Apis cerana cerana

Glutathione S-transferases (GSTs) are members of a multifunctional antioxidant enzyme superfamily that play pivotal roles in both detoxification and protection against oxidative damage caused by reactive oxygen species. In this study, a complementary DNA (cDNA) encoding a sigma class GST was identified in the Chinese honey bee, Apis cerana cerana (AccGSTS1). AccGSTS1 was constitutively expressed in all tissues of adult worker bees, including the brain, fat body, epidermis, muscle, and midgut, with particularly robust transcription in the fat body. Relative messenger RNA expression levels of AccGSTS1 at different developmental stages varied, with the highest levels of expression observed in adults. The potential function of AccGSTS1 in cellular defenses against abiotic stresses (cold, heat, UV, H2O2, HgCl2, and insecticides) was investigated. AccGSTS1 was significantly upregulated in response to all of the treatment conditions examined, although the induction levels were varied. Recombinant AccGSTS1 protein showed characteristic glutathione-conjugating catalytic activity toward 1-chloro-2,4-dinitrobenzene. Functional assays revealed that AccGSTS1 could remove H2O2, thereby protecting DNA from oxidative damage. Escherichia coli overexpressing AccGSTS1 showed long-term resistance under conditions of oxidative stress. Together, these results suggest that AccGSTS1 is a crucial antioxidant enzyme involved in cellular antioxidant defenses and honey bee survival.

Transcriptomic and proteomic responses of sweetpotato whitefly, Bemisia tabaci, to thiamethoxam

BACKGROUND

The sweetpotato whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae), is one of the most widely distributed agricultural pests. Although it has developed resistance to many registered insecticides including the neonicotinoid insecticide thiamethoxam, the mechanisms that regulate the resistance are poorly understood. To understand the molecular basis of thiamethoxam resistance, "omics" analyses were carried out to examine differences between resistant and susceptible B. tabaci at both transcriptional and translational levels.

RESULTS

A total of 1,338 mRNAs and 52 proteins were differentially expressed between resistant and susceptible B. tabaci. Among them, 11 transcripts had concurrent transcription and translation profiles. KEGG analysis mapped 318 and 35 differentially expressed genes and proteins, respectively, to 160 and 59 pathways (p<0.05). Thiamethoxam treatment activated metabolic pathways (e.g., drug metabolism), in which 118 transcripts were putatively linked to insecticide resistance, including up-regulated glutathione-S-transferase, UDP glucuronosyltransferase, glucosyl/glucuronosyl transferase, and cytochrome P450. Gene Ontology analysis placed these genes and proteins into protein complex, metabolic process, cellular process, signaling, and response to stimulus categories. Quantitative real-time PCR analysis validated "omics" response, and suggested a highly overexpressed P450, CYP6CX1, as a candidate molecular basis for the mechanistic study of thiamethoxam resistance in whiteflies. Finally, enzymatic activity assays showed elevated detoxification activities in the resistant B. tabaci.

CONCLUSIONS

This study demonstrates the applicability of high-throughput omics tools for identifying molecular candidates related to thiamethoxam resistance in an agricultural important insect pest. In addition, transcriptomic and proteomic analyses provide a solid foundation for future functional investigations into the complex molecular mechanisms governing the neonicotinoid resistance in whiteflies.

Alpha-class glutathione S-transferases in wild turkeys (Meleagris gallopavo): characterization and role in resistance to the carcinogenic mycotoxin aflatoxin B1

Domestic turkeys (Meleagris gallopavo) are one of the most susceptible animals known to the toxic effects of the mycotoxin aflatoxin B₁ (AFB₁), a potent human hepatocarcinogen, and universal maize contaminant. We have demonstrated that such susceptibility is associated with the inability of hepatic glutathione S-transferases (GSTs) to detoxify the reactive electrophilic metabolite exo-AFB₁-8,9-epoxide (AFBO). Unlike their domestic counterparts, wild turkeys, which are relatively AFB₁-resistant, possess hepatic GST-mediated AFBO conjugating activity_. Here, we characterized the molecular and functional properties of hepatic alpha-class GSTs (GSTAs) from wild and domestic turkeys to shed light on the differences in resistance between these closely related strains. Six alpha-class GST genes (GSTA) amplified from wild turkeys (Eastern and Rio Grande subspecies), heritage breed turkeys (Royal Palm) and modern domestic (Nicholas strain) turkeys were sequenced, and catalytic activities of heterologously-expressed recombinant enzymes determined. Alpha-class identity was affirmed by conserved GST domains and four signature motifs. All GSTAs contained single nucleotide polymorphisms (SNPs) in their coding regions: GSTA1.1 (5 SNPs), GSTA1.2 (7), GSTA1.3 (3), GSTA2 (3), GSTA3 (1) and GSTA4 (2). E. coli-expressed GSTAs possessed varying activities toward GST substrates 1-chloro-2,4-dinitrobenzene (CDNB), 1,2-dichloro-4-nitrobenzene (DCNB), ethacrynic acid (ECA), cumene hydroperoxide (CHP). As predicted by their relative resistance, livers from domestic turkeys lacked detectable GST-mediated AFBO detoxification activity, whereas those from wild and heritage birds possessed this critical activity, suggesting that intensive breeding and selection resulted in loss of AFB₁-protective alleles during domestication. Our observation that recombinant tGSTAs detoxify AFBO, whereas their hepatic forms do not, implies that the hepatic forms of these enzymes are down-regulated, silenced, or otherwise modified by one or more mechanisms. These data may inform of possible molecular mechanisms of resistance to AFB₁, and may also have the benefit of identifying genetic markers which could be used to enhance AFB₁ resistance in modern domestic strains.

Using Drosophila melanogaster to validate metabolism-based insecticide resistance from insect pests

Identifying molecular mechanisms of insecticide resistance is important for preserving insecticide efficacy, developing new insecticides and implementing insect control. The metabolic detoxification of insecticides is a widespread resistance mechanism. Enzymes with the potential to detoxify insecticides are commonly encoded by members of the large cytochrome P450, glutathione S-transferase and carboxylesterase gene families, all rapidly evolving in insects. Here, we demonstrate that the model insect Drosophila melanogaster is useful for functionally validating the role of metabolic enzymes in conferring metabolism-based insecticide resistance. Alleles of three well-characterized genes from different pest insects were expressed in transgenic D. melanogaster : a carboxylesterase gene (αE7) from the Australian sheep blowfly Lucilia cuprina, a glutathione S-transferase gene (GstE2) from the mosquito Anopheles gambiae and a cytochrome P450 gene (Cyp6cm1) from the whitefly Bemisia tabaci. For all genes, expression in D. melanogaster resulted in insecticide resistance phenotypes mirroring those observed in resistant populations of the pest species. Using D. melanogaster to assess the potential for novel metabolic resistance mechanisms to evolve in pest species is discussed.