Characterization of the inheritance of field-evolved resistance to diamides in the fall armyworm (Spodoptera frugiperda) (Lepidoptera: Noctuidae) population from Puerto Rico

The fall armyworm (Spodoptera frugiperda) is one of the most destructive pests of corn. New infestations have been reported in the East Hemisphere, reaching India, China, Malaysia, and Australia, causing severe destruction to corn and other crops. In Puerto Rico, practical resistance to different mode of action compounds has been reported in cornfields. In this study, we characterized the inheritance of resistance to chlorantraniliprole and flubendiamide and identified the possible cross-resistance to cyantraniliprole and cyclaniliprole. The Puerto Rican (PR) strain showed high levels of resistance to flubendiamide (RR50 = 2,762-fold) and chlorantraniliprole (RR50 = 96-fold). The inheritance of resistance showed an autosomal inheritance for chlorantraniliprole and an X-linked inheritance for flubendiamide. The trend of the dominance of resistance demonstrated an incompletely recessive trait for H1 (♂ SUS × ♀ PR) × and an incompletely dominant trait for H2 (♀ SUS × ♂ PR) × for flubendiamide and chlorantraniliprole. The PR strain showed no significant presence of detoxification enzymes (using synergists: PBO, DEF, DEM, and VER) to chlorantraniliprole; however, for flubendiamide the SR = 2.7 (DEM), SR = 3.2 (DEF) and SR = 7.6 (VER) indicated the role of esterases, glutathione S- transferases and ABC transporters in the metabolism of flubendiamide. The PR strain showed high and low cross-resistance to cyantraniliprole (74-fold) and cyclaniliprole (11-fold), respectively. Incomplete recessiveness might lead to the survival of heterozygous individuals when the decay of diamide residue occurs in plant tissues. These results highlight the importance of adopting diverse pest management strategies, including insecticide rotating to manage FAW populations in Puerto Rico and other continents.

Identification of differentially expressed miRNAs associated with diamide detoxification pathways in Spodoptera frugiperda

The fall armyworm (FAW) Spodoptera frugiperda is a severe economic pest of multiple crops globally. Control of this pest is often achieved using insecticides; however, over time, S. frugiperda has developed resistance to new mode of action compounds, including diamides. Previous studies have indicated diamide resistance is a complex developmental process involving multiple detoxification genes. Still, the mechanism underlying the possible involvement of microRNAs in post-transcriptional regulation of resistance has not yet been elucidated. In this study, a global screen of microRNAs (miRNAs) revealed 109 known and 63 novel miRNAs. Nine miRNAs (four known and five novel) were differentially expressed between insecticide-resistant and -susceptible strains. Gene Ontology analysis predicted putative target transcripts of the differentially expressed miRNAs encoding significant genes belonging to detoxification pathways. Additionally, miRNAs are involved in response to diamide exposure, indicating they are probably associated with the detoxification pathway. Thus, this study provides comprehensive evidence for the link between repressed miRNA expression and induced target transcripts that possibly mediate diamide resistance through post-transcriptional regulation. These findings highlight important clues for further research to unravel the roles and mechanisms of miRNAs in conferring diamide resistance.

The insulin signaling pathway in Drosophila melanogaster: A nexus revealing an "Achilles' heel" in DDT resistance

Insecticide resistance is an ongoing challenge in agriculture and disease vector control. Here, we demonstrate a novel strategy to attenuate resistance. We used genomics tools to target fundamental energy-associated pathways and identified a potential "Achilles' heel" for resistance, a resistance-associated protein that, upon inhibition, results in a substantial loss in the resistance phenotype. Specifically, we compared the gene expression profiles and structural variations of the insulin/insulin-like growth factor signaling (IIS) pathway genes in DDT-susceptible (91-C) and -resistant (91-R) Drosophila melanogaster (Drosophila) strains. A total of eight and seven IIS transcripts were up- and down-regulated, respectively, in 91-R compared to 91-C. A total of 114 nonsynonymous mutations were observed between 91-C and 91-R, of which 51.8% were fixed. Among the differentially expressed transcripts, phosphoenolpyruvate carboxykinase (PEPCK), down-regulated in 91-R, encoded the greatest number of amino acid changes, prompting us to perform PEPCK inhibitor-pesticide exposure bioassays. The inhibitor of PEPCK, hydrazine sulfate, resulted in a 161- to 218-fold decrease in the DDT resistance phenotype (91-R) and more than a 4- to 5-fold increase in susceptibility in 91-C. A second target protein, Glycogen synthase kinase 3β (GSK3β-PO), had one amino acid difference between 91-C and 91-R, and the corresponding transcript was also down-regulated in 91-R. A GSK3β-PO inhibitor, lithium chloride, likewise reduced the resistance but to a lesser extent than did hydrazine sulfate for PEPCK. We demonstrate the potential role of IIS genes in DDT resistance and the potential discovery of an "Achilles' heel" against pesticide resistance in this pathway.

Dietary antioxidants impact DDT resistance in Drosophila melanogaster

Insects experience a diversity of subtoxic and/or toxic xenobiotics through exposure to pesticides and, in the case of herbivorous insects, through plant defensive compounds in their diets. Many insects are also concurrently exposed to antioxidants in their diets. The impact of dietary antioxidants on the toxicity of xenobiotics in insects is not well understood, in part due to the challenge of developing appropriate systems in which doses and exposure times (of both the antioxidants and the xenobiotics) can be controlled and outcomes can be easily measured. However, in Drosophila melanogaster, a well-established insect model system, both dietary factors and pesticide exposure can be easily controlled. Additionally, the mode of action and xenobiotic metabolism of dichlorodiphenyltrichloroethane (DDT), a highly persistent neurotoxic organochlorine insecticide that is detected widely in the environment, have been well studied in DDT-susceptible and -resistant strains. Using a glass-vial bioassay system with blue diet as the food source, seven compounds with known antioxidant effects (ascorbic acid, β-carotene, glutathione, α-lipoic acid, melatonin, minocycline, and serotonin) were orally tested for their impact on DDT toxicity across three strains of D. melanogaster: one highly susceptible to DDT (Canton-S), one mildly susceptible (91-C), and one highly resistant (91-R). Three of the antioxidants (serotonin, ascorbic acid, and β-carotene) significantly impacted the toxicity of DDT in one or more strains. Fly strain and gender, antioxidant type, and antioxidant dose all affected the relative toxicity of DDT. Our work demonstrates that dietary antioxidants can potentially alter the toxicity of a xenobiotic in an insect population.

Identification of transcriptional responsive genes to acetic acid, ethanol, and 2-phenylethanol exposure in Drosophila melanogaster

The fruit fly, Drosophila melanogaster, is predominantly found in overripe, rotten, fermenting, or decaying fruits and is constantly exposed to chemical stressors such as acetic acid, ethanol, and 2-phenylethanol. D. melanogaster has been employed as a model system for studying the molecular bases of various types of chemical-induced tolerance. Expression profiling using Illumina sequencing has been performed for identifying changes in gene expression that may be associated with evolutionary adaptation to exposure of acetic acid, ethanol, and 2-phenylethanol. We identified a total of 457 differentially expressed genes that may affect sensitivity or tolerance to three chemicals in the chemical treatment group as opposed to the control group. Gene-set enrichment analysis revealed that the genes involved in metabolism, multicellular organism reproduction, olfaction, regulation of signal transduction, and stress tolerance were over-represented in response to chemical exposure. Furthermore, we also detected a coordinated upregulation of genes in the Toll- and Imd-signaling pathways after the chemical exposure. Quantitative reverse transcription PCR analysis revealed that the expression levels of nine genes within the set of genes identified by RNA sequencing were up- or downregulated owing to chemical exposure. Taken together, our data suggest that such differentially expressed genes are coordinately affected by chemical exposure. Transcriptional analyses after exposure of D. melanogaster with three chemicals provide unique insights into subsequent functional studies on the mechanisms underlying the evolutionary adaptation of insect species to environmental chemical stressors.

Efficacy of a cry1Ab Gene for Control of Maruca vitrata (Lepidoptera: Crambidae) in Cowpea (Fabales: Fabaceae)

Cowpea [Vigna unguiculata (L) Walp.] is an important staple legume in the diet of many households in sub-Saharan Africa. Its production, however, is negatively impacted by many insect pests including bean pod borer, Maruca vitrata F., which can cause 20-80% yield loss. Several genetically engineered cowpea events that contain a cry1Ab gene from Bacillus thuringiensis (Bt) for resistance against M. vitrata were evaluated in Nigeria, Burkina Faso, and Ghana (West Africa), where cowpea is commonly grown. As part of the regulatory safety package, these efficacy data were developed and evaluated by in-country scientists. The Bt-cowpea lines were planted in confined field trials under Insect-proof netting and artificially infested with up to 500 M. vitrata larvae per plant during bud formation and flowering periods. Bt-cowpea lines provided nearly complete pod and seed protection and in most cases resulted in significantly increased seed yield over non-Bt control lines. An integrated pest management strategy that includes use of Bt-cowpea augmented with minimal insecticide treatment for protection against other insects is recommended to control pod borer to enhance cowpea production. The insect resistance management plan is based on the high-dose refuge strategy where non-Bt-cowpea and natural refuges are expected to provide M. vitrata susceptible to Cry1Ab protein. In addition, there will be a limited release of this product until a two-toxin cowpea pyramid is released. Other than South African genetically engineered crops, Bt-cowpea is the first genetically engineered food crop developed by the public sector and approved for release in sub-Saharan Africa.

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.

Impacts of Sub-lethal DDT Exposures on microRNA and Putative Target Transcript Expression in DDT Resistant and Susceptible Drosophila melanogaster Strains

Ten constitutively differentially expressed miRNAs were previously described between DDT-resistant 91-R and -susceptible control Drosophila melanogaster strains, and among their predicted target genes were those associated with metabolic DDT resistance mechanisms. The present study evaluated the inducibility of miRNA expression and putative downstream regulation of cytochrome P450s in response to DDT exposure in a time-dependent manner in 91-R and the susceptible Canton-S strain. Specifically, RT-qPCR analysis showed that DDT exposures led to the significant down-regulation (repression) of miR-310-3p, miR-311-3p, miR-312-3p, miR-313-3p, and miR-92a-3p levels in Canton-S. This is contrasted with the lack of significant changes in 91-R at most time-points following DDT exposure. The levels of expression among miRNAs exhibited opposite expression patterns compared to their corresponding putative target cytochrome P450s at the same time points after DDT exposure. Collectively, results from this study suggest that miR-310-3p, miR-311-3p, miR-312-3p, miR-313-3p, and miR-92a-3p might have a potential role in the control of DDT detoxification through the post-transcriptional regulation of target cytochrome P450s in Canton-S. Conversely, the lack of significant changes of these same miRNAs in 91-R following DDT-exposure suggests a possible adaptive mutation that removes repressive control mechanisms. These data are important for the understanding impact of adaptive changes in miRNA expression on post-transcriptional regulatory mechanism involved in the evolution of DDT resistance in 91-R.

Effects of Volatiles from Clavigralla tomentosicollis Stål. (Hemiptera: Coreidae) Adults on the Host Location Behavior of the Egg Parasitoid Gryon fulviventre (Crawford) (Hymenoptera: Scelionidae)

The egg parasitoid Gryon fulviventre is a potential biological control agent of Clavigralla tomentosicollis, a coreid pod-sucking pest of Vigna unguiculata. The host location behavior of naive parasitoid females was studied using a four-armed olfactometer. Two strains of G. fulviventre parasitoids from Burkina Faso and Benin were exposed to odors provided by healthy and infested pods as well as C. tomentosicollis females and males. The time spent in each odor zone was recorded to determine the preference of parasitoid females. Results show that odors from healthy pods, infested pods, and pest females did not attract the parasitoid. However, a significantly attractive response of both strains of G. fulviventre was recorded in the presence of volatiles from males of C. tomentosicollis. Moreover, experiments testing G. fulviventre females' behavior when simultaneously exposed to volatiles from cowpea pods (healthy and infested) and increasing numbers of C. tomentosicollis males revealed a significantly higher attraction of parasitoid females of both strains by volatiles from ten males of C. tomentosicollis. The results suggest that the males of the insect pest emit a pheromone used as kairomone by parasitoids to locate their host. The conditions determining this attractiveness at field level and its impact on host-searching efficiency are discussed.

Global Solutions for International Development Partnerships: Beyond Insider/Outsider Binaries

In recent years, information and communication technologies for development (ICT4D) approaches have facilitated international development work, but still more effective ICT4D deployments are needed. This article examines how one ICT4D initiative, Scientific Animations without Borders (SAWBO), works with partners in Africa not only to transcend problematic insider / outsider binaries that impact solution delivery but also to implement inter-organizational collaborations on research and mission-critical knowledge-transfer goals that effectively reach the widest diversity of target populations.

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.

Identification and interaction of multiple genes resulting in DDT resistance in the 91-R strain of Drosophila melanogaster by RNAi approaches

4,4'-dichlorodiphenyltrichloroethane (DDT) has been re-recommended by the World Health Organization for malaria mosquito control. Previous DDT use has resulted in resistance, and with continued use resistance will likely increase in terms of level and extent. Drosophila melanogaster is a model dipteran with a well annotated genome allowing both forward and reverse genetic manipulation, numerous studies done on insecticide resistance mechanisms, and is related to malaria mosquitoes allowing for extrapolation. The 91-R strain of D. melanogaster is highly resistant to DDT (>1500-fold) and recently, reduced penetration, increased detoxification, and direct excretion have been identified as resistance mechanisms. Their interactions, however, remain unclear. Use of Gal4/UAS-RNAi transgenic lines of D. melanogaster allowed for the targeted knockdown of genes putatively involved in DDT resistance and has identified the role of several cuticular proteins (Cyp4g1 and Lcp1), cytochrome P450 monooxygenases (Cyp6g1 and Cyp12d1), and ATP binding cassette transporters (Mdr50, Mdr65, and Mrp1) involved in decreased sensitivity to DDT. These above findings have been further validated in 91-R flies using a nanoparticle-enhanced RNAi strategy, directly implication these genes in DDT resistance in 91-R flies.

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.

Variability Assessment of Aromatic Rice Germplasm by Pheno-Genomic traits and Population Structure Analysis

While the pleasant scent of aromatic rice is making it more popular, with demand for aromatic rice expected to rise in future, varieties of this have low yield potential. Genetic diversity and population structure of aromatic germplasm provide valuable information for yield improvement which has potential market value and farm profit. Here, we show diversity and population structure of 113 rice germplasm based on phenotypic and genotypic traits. Phenotypic traits showed that considerable variation existed across the germplasm. Based on Shannon-Weaver index, the most variable phenotypic trait was lemma-palea color. Detecting 140 alleles, 11 were unique and suitable as a germplasm diagnostic tool. Phylogenetic cluster analysis using genotypic traits classified germplasm into three major groups. Moreover, model-based population structure analysis divided all germplasm into three groups, confirmed by principal component and neighbors joining tree analyses. An analysis of molecular variance (AMOVA) and pairwise FST test showed significant differentiation among all population pairs, ranging from 0.023 to 0.068, suggesting that all three groups differed. Significant correlation coefficient was detected between phenotypic and genotypic traits which could be valuable to select further improvement of germplasm. Findings from this study have the potential for future use in aromatic rice molecular breeding programs.

Monitoring of the Cowpea Bruchid, Callosobruchus maculatus (Coleoptera: Bruchidae), Feeding Activity in Cowpea Seeds: Advances in Sensing Technologies Reveals New Insights

Cowpea provides a significant source of protein for over 200 million people in Sub-Saharan Africa. The cowpea bruchid, Callosobruchus maculatus (F) (Coleoptera: Bruchidae), is a major pest of cowpea as the larval stage attacks stored cowpea grains, causing postharvest loss. Cowpea bruchid larvae spend all their time feeding within the cowpea seed. Past research findings, published over 25 yr ago, have shown that feeding activity of several bruchids within a cowpea seed emit mechanical vibrations within the frequency range 5-75 kHz. This work led to the development of monitoring technologies that are both important for basic research and practical application. Here, we use newer and significantly improved technologies to re-explore the nature of the vibration signals produced by an individual C. maculatus, when it feeds in cowpea seeds. Utilizing broadband frequency sensing, individual fourth-instar bruchid larvae feeding activities (vibration events) were recorded to identify specific key emission frequencies. Verification of recorded events and association to actual feeding activities was achieved through mass measurements over 24 h for a series of replicates. The measurements identified variable peak event emission frequencies across the replicate sample set ranging in frequency from 16.4 to 26.5 kHz. A positive correlation between the number of events recorded and the measured mass loss of the cowpea seed was observed. The procedure and verification reported in this work provide an improved basis for laboratory-based monitoring of single larval feeding. From the rich dataset captured, additional analysis can be carried out to identify new key variables of hidden bruchid larval activity.

Phylogenomics from Whole Genome Sequences Using aTRAM

Novel sequencing technologies are rapidly expanding the size of data sets that can be applied to phylogenetic studies. Currently the most commonly used phylogenomic approaches involve some form of genome reduction. While these approaches make assembling phylogenomic data sets more economical for organisms with large genomes, they reduce the genomic coverage and thereby the long-term utility of the data. Currently, for organisms with moderate to small genomes ($<$1000 Mbp) it is feasible to sequence the entire genome at modest coverage ($10-30\times$). Computational challenges for handling these large data sets can be alleviated by assembling targeted reads, rather than assembling the entire genome, to produce a phylogenomic data matrix. Here we demonstrate the use of automated Target Restricted Assembly Method (aTRAM) to assemble 1107 single-copy ortholog genes from whole genome sequencing of sucking lice (Anoplura) and out-groups. We developed a pipeline to extract exon sequences from the aTRAM assemblies by annotating them with respect to the original target protein. We aligned these protein sequences with the inferred amino acids and then performed phylogenetic analyses on both the concatenated matrix of genes and on each gene separately in a coalescent analysis. Finally, we tested the limits of successful assembly in aTRAM by assembling 100 genes from close- to distantly related taxa at high to low levels of coverage.Both the concatenated analysis and the coalescent-based analysis produced the same tree topology, which was consistent with previously published results and resolved weakly supported nodes. These results demonstrate that this approach is successful at developing phylogenomic data sets from raw genome sequencing reads. Further, we found that with coverages above $5-10\times$, aTRAM was successful at assembling 80-90% of the contigs for both close and distantly related taxa. As sequencing costs continue to decline, we expect full genome sequencing will become more feasible for a wider array of organisms, and aTRAM will enable mining of these genomic data sets for an extensive variety of applications, including phylogenomics. [aTRAM; gene assembly; genome sequencing; phylogenomics.].

Lethal and sub-lethal effects of select macrocyclic lactones insecticides on forager worker honey bees under laboratory experimental conditions

Selective insecticide application is one important strategy for more precisely targeting harmful insects while avoiding or mitigating collateral damage to beneficial insects like honey bees. Recently, macrocyclic lactone-class insecticides have been introduced into the market place as selective bio-insecticides for controlling many arthropod pests, but how to target this selectivity only to harmful insects has yet to be achieved. In this study, the authors investigated the acute toxicity of fourmacrocyclic lactone insecticides (commercialized as abamectin, emamectin benzoate, spinetoram, and spinosad) both topically and through feeding studies with adult forager honey bees. Results indicated emamectin benzoate as topically 133.3, 750.0, and 38.3-fold and orally 3.3, 7.6, and 31.7-fold more toxic, respectively than abamectin, spinetoram and spinosad. Using Hazard Quotients for estimates of field toxicity, abamectin was measured as the safest insecticide both topically and orally for honey bees. Moreover, a significant reduction of sugar solution consumption by treatment group honey bees for orally applied emamectin benzoate and spinetoram suggests that these insecticides may have repellent properties.

Molecular Ecological Basis of Grasshopper (Oedaleus asiaticus) Phenotypic Plasticity under Environmental Selection

While ecological adaptation in insects can be reflected by plasticity of phenotype, determining the causes and molecular mechanisms for phenotypic plasticity (PP) remains a crucial and still difficult question in ecology, especially where control of insect pests is involved. Oedaleus asiaticus is one of the most dominant pests in the Inner Mongolia steppe and represents an excellent system to study phenotypic plasticity. To better understand ecological factors affecting grasshopper phenotypic plasticity and its molecular control, we conducted a full transcriptional screening of O. asiaticus grasshoppers reared in four different grassland patches in Inner Mongolia. Grasshoppers showed different degrees of PP associated with unique gene expressions and different habitat plant community compositions. Grasshopper performance variables were susceptible to habitat environment conditions and closely associated with plant architectures. Intriguingly, eco-transcriptome analysis revealed five potential candidate genes playing important roles in grasshopper performance, with gene expression closely relating to PP and plant community factors. By linking the grasshopper performances to gene profiles and ecological factors using canonical regression, we first demonstrated the eco-transcriptomic architecture (ETA) of grasshopper phenotypic traits (ETAGPTs). ETAGPTs revealed plant food type, plant density, coverage, and height were the main ecological factors influencing PP, while insect cuticle protein (ICP), negative elongation factor A (NELFA), and lactase-phlorizin hydrolase (LCT) were the key genes associated with PP. Our study gives a clear picture of gene-environment interaction in the formation and maintenance of PP and enriches our understanding of the transcriptional events underlying molecular control of rapid phenotypic plasticity associated with environmental variability. The findings of this study may also provide new targets for pest control and highlight the significance of ecological management practice on grassland conservation.

The complete mitochondrial genome of Anoplocnemis curvipes F. (Coreinea, Coreidae, Heteroptera), a pest of fresh cowpea pods

The complete 16,345-bp mitochondrial genome of the agriculturally destructive pod sucking pest, the giant coreid bug, Anoplocnemis curvipes (Hemiptera: Coreidae), was assembled from paired-end Illumina HiSeq 2500 reads. The A. curvipes mitochondrial genome consists of 13 protein coding genes (PCGs), 22 tRNAs, 2 rRNAs and a control region in the order and orientation typical among insects. PCG initiation codons (ATG, ATC, ATT and ATA) with termination codon (TAA) are used with the exception of TAG stop codons by Cytb and ND3. All tRNA genes fold into predicted cloverleaf secondary structures having requisite triplets on the anticodon loop, apart from tRNA-Ser1 (AGN) whose dihydrouridine (DHU) arm forms a simple loop. The phylogenetic analysis of hemipteran mitogenomes clusters to the family level and supports the monophyly of the five superfamilies in Pentatomomorpha of Hemiptera. The Coreoidea and Pyrrhocoroidea are sister groups, while Coreidae and Alydidae are sister groups to Rhopalidae. These analyses provide insight to mitogenomics and evolutionary relationships among pentatomoid insects.

Comparison of the proliferation and excretion of Bartonella quintana between body and head lice following oral challenge

Human body and head lice are highly related haematophagous ectoparasites but only the body louse has been shown to transmit Bartonella quintana, the causative agent of trench fever. The mechanisms by which body lice became a vector for B. quintana, however, are poorly understood. Following oral challenge, green fluorescent protein-expressing B. quintana proliferated over 9 days postchallenge with the number of bacteria being significantly higher in whole body vs. head lice. The numbers of B. quintana detected in faeces from infected lice, however, were approximately the same in both lice. Nevertheless, the viability of B. quintana was significantly higher in body louse faeces. Comparison of immune responses in alimentary tract tissues revealed that basal transcription levels of peptidoglycan recognition protein and defensins were lower in body lice and the transcription of defensin 1 was up-regulated by oral challenge with wild-type B. quintana in head but not in body lice. In addition, the level of cytotoxic reactive oxygen species generated by epithelial cells was significantly lower in body lice. Although speculative at this time, the reduced immune response is consistent with the higher vector competence seen in body vs. head lice in terms of B. quintana infection.

Plant Defense Inhibitors Affect the Structures of Midgut Cells in Drosophila melanogaster and Callosobruchus maculatus

Plants produce proteins such as protease inhibitors and lectins as defenses against herbivorous insects and pathogens. However, no systematic studies have explored the structural responses in the midguts of insects when challenged with plant defensive proteins and lectins across different species. In this study, we fed two kinds of protease inhibitors and lectins to the fruit fly Drosophila melanogaster and alpha-amylase inhibitors and lectins to the cowpea bruchid Callosobruchus maculatus. We assessed the changes in midgut cell structures by comparing them with such structures in insects receiving normal diets or subjected to food deprivation. Using light and transmission electron microscopy in both species, we observed structural changes in the midgut peritrophic matrix as well as shortened microvilli on the surfaces of midgut epithelial cells in D. melanogaster. Dietary inhibitors and lectins caused similar lesions in the epithelial cells but not much change in the peritrophic matrix in both species. We also noted structural damages in the Drosophila midgut after six hours of starvation and changes were still present after 12 hours. Our study provided the first evidence of key structural changes of midguts using a comparative approach between a dipteran and a coleopteran. Our particular observation and discussion on plant-insect interaction and dietary stress are relevant for future mode of action studies of plant defensive protein in insect physiology.

Alternative Splice in Alternative Lice

Genomic and transcriptomics analyses have revealed human head and body lice to be almost genetically identical; although con-specific, they nevertheless occupy distinct ecological niches and have differing feeding patterns. Most importantly, while head lice are not known to be vector competent, body lice can transmit three serious bacterial diseases; epidemictyphus, trench fever, and relapsing fever. In order to gain insights into the molecular bases for these differences, we analyzed alternative splicing (AS) using next-generation sequencing data for one strain of head lice and one strain of body lice. We identified a total of 3,598 AS events which were head or body lice specific. Exon skipping AS events were overrepresented among both head and body lice, whereas intron retention events were underrepresented in both. However, both the enrichment of exon skipping and the underrepresentation of intron retention are significantly stronger in body lice compared with head lice. Genes containing body louse-specific AS events were found to be significantly enriched for functions associated with development of the nervous system, salivary gland, trachea, and ovarian follicle cells, as well as regulation of transcription. In contrast, no functional categories were overrepresented among genes with head louse-specific AS events. Together, our results constitute the first evidence for transcript pool differences in head and body lice, providing insights into molecular adaptations that enabled human lice to adapt to clothing, and representing a powerful illustration of the pivotal role AS can play in functional adaptation.

RNAi validation of resistance genes and their interactions in the highly DDT-resistant 91-R strain of Drosophila melanogaster

4,4'-dichlorodiphenyltrichloroethane (DDT) has been re-recommended by the World Health Organization for malaria mosquito control. Previous DDT use has resulted in resistance, and with continued use resistance will increase in terms of level and extent. Drosophila melanogaster is a model dipteran that has many available genetic tools, numerous studies done on insecticide resistance mechanisms, and is related to malaria mosquitoes allowing for extrapolation. The 91-R strain of D. melanogaster is highly resistant to DDT (>1500-fold), however, there is no mechanistic scheme that accounts for this level of resistance. Recently, reduced penetration, increased detoxification, and direct excretion have been identified as resistance mechanisms in the 91-R strain. Their interactions, however, remain unclear. Use of UAS-RNAi transgenic lines of D. melanogaster allowed for the targeted knockdown of genes putatively involved in DDT resistance and has validated the role of several cuticular proteins (Cyp4g1 and Lcp1), cytochrome P450 monooxygenases (Cyp6g1 and Cyp12d1), and ATP binding cassette transporters (Mdr50, Mdr65, and Mrp1) involved in DDT resistance. Further, increased sensitivity to DDT in the 91-R strain after intra-abdominal dsRNA injection for Mdr50, Mdr65, and Mrp1 was determined by a DDT contact bioassay, directly implicating these genes in DDT efflux and resistance.

Utilization of the human louse genome to study insecticide resistance and innate immune response

Since sequencing the human body louse genome, substantial advances have occurred in the utilization of the information gathered from louse genomes and transcriptomes. Comparatively, the body louse genome contains far fewer genes involved in environmental response, such as xenobiotic detoxification and innate immune response. Additionally, the body louse maintains a primary bacterial endosymbiont, Candidatus Riesia pediculicola, and a number of bacterial pathogens that it vectors, which have genomes that are also reduced in size. Thus, human louse genomes offer unique information and tools for use in advancing our understanding of coevolution among vectors, endosymbionts and pathogens. In this review, we summarize the current literature on the extent of pediculicide resistance, the availability of new pediculicides and information establishing this organism as an efficient model to study how xenobiotic metabolism, which is involved in insecticide resistance, is induced and how insects modify their innate immune response upon bacterial challenge resulting in enhanced vector competence.

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.

A glycine insertion in the estrogen-related receptor (ERR) is associated with enhanced expression of three cytochrome P450 genes in transgenic Drosophila melanogaster

Insecticide-resistant Drosophila melanogaster strains represent a resource for the discovery of the underlying molecular mechanisms of cytochrome P450 constitutive over-expression, even if some of these P450s are not directly involved in the resistance phenotype. For example, in select 4,4'-dichlorodiphenyltrichloroethane (DDT) resistant strains the glucocorticoid receptor-like (GR-like) potential transcription factor binding motifs (TFBMs) have previously been shown to be associated with constitutively differentially-expressed cytochrome P450s, Cyp12d1, Cyp6g2 and Cyp9c1. However, insects are not known to have glucocorticoids. The only ortholog to the mammalian glucocorticoid receptor (GR) in D. melanogaster is an estrogen-related receptor (ERR) gene, which has two predicted alternative splice isoforms (ERRa and ERRb). Sequencing of ERRa and ERRb in select DDT susceptible and resistant D. melanogaster strains has revealed a glycine (G) codon insertion which was only observed in the ligand binding domain of ERR from the resistant strains tested (ERR-G). Transgenic flies, expressing the ERRa-G allele, constitutively over-expressed Cyp12d1, Cyp6g2 and Cyp9c1. Only Cyp12d1 and Cyp6g2 were over-expressed in the ERRb-G transgenic flies. Phylogenetic studies show that the G-insertion appeared to be located in a less conserved domain in ERR and this insertion is found in multiple species across the Sophophora subgenera.

Transcriptomic response of cowpea bruchids to N-acetylglucosamine-specific lectins

Griffonia simplicifolia lectin II (GSII) and wheat germ agglutinin (WGA) are N-acetylglucosamine-binding lectins. Previous studies demonstrated that they have anti-insect activity, a property potentially useful in pest control. To gain some insight into the insect response to dietary lectins, we performed transcriptomic analysis using the cowpea bruchid (Callosobruchus maculatus) midgut microarray platform we built. Compared to the nonnutritional cellulose treatment, dietary lectins induced more profound changes in gene expression. Ingestion of relatively high doses of lectins for 24 h resulted in alteration of gene expression involved in sugar and lipid metabolism, transport, development, defense, and stress tolerance. Metabolic genes were largely downregulated. Moreover, we observed disorganized microvilli resulting from ingestion of WGA. This morphological change is consistent with the lectin-induced changes in genes related to midgut epithelial cell repair. In addition, suboptimal nutrient conditions may serve as a stress signal to trigger senescence processes, leading to growth arrest and developmental delay.

Socially responsive effects of brain oxidative metabolism on aggression

Despite ongoing high energetic demands, brains do not always use glucose and oxygen in a ratio that produces maximal ATP through oxidative phosphorylation. In some cases glucose consumption exceeds oxygen use despite adequate oxygen availability, a phenomenon known as aerobic glycolysis. Although metabolic plasticity seems essential for normal cognition, studying its functional significance has been challenging because few experimental systems link brain metabolic patterns to distinct behavioral states. Our recent transcriptomic analysis established a correlation between aggression and decreased whole-brain oxidative phosphorylation activity in the honey bee (Apis mellifera), suggesting that brain metabolic plasticity may modulate this naturally occurring behavior. Here we demonstrate that the relationship between brain metabolism and aggression is causal, conserved over evolutionary time, cell type-specific, and modulated by the social environment. Pharmacologically treating honey bees to inhibit complexes I or V in the oxidative phosphorylation pathway resulted in increased aggression. In addition, transgenic RNAi lines and genetic manipulation to knock down gene expression in complex I in fruit fly (Drosophila melanogaster) neurons resulted in increased aggression, but knockdown in glia had no effect. Finally, honey bee colony-level social manipulations that decrease individual aggression attenuated the effects of oxidative phosphorylation inhibition on aggression, demonstrating a specific effect of the social environment on brain function. Because decreased neuronal oxidative phosphorylation is usually associated with brain disease, these findings provide a powerful context for understanding brain metabolic plasticity and naturally occurring behavioral plasticity.

Nrf2/Maf-binding-site-containing functional Cyp6a2 allele is associated with DDT resistance in Drosophila melanogaster

BACKGROUND

Increased insecticide detoxification mediated by cytochrome P450s is a common mechanism of insecticide resistance. Although Cyp6a2 has been observed to be overexpressed in many 4,4'-dichlorodiphenyltrichloroethane (DDT)-resistant strains of Drosophila melanogaster, how Cyp6a2 is regulated and whether its overproduction confers DDT resistance remain elusive.

RESULTS

Molecular analysis identified five Cyp6a2 alleles (Cyp6a2(Canton) (-S-1) , Cyp6a2(Canton) (-S-2) , Cyp6a2(91-C) , Cyp6a2(91-R) and Cyp6a2(Wisconsin) (-) (WD) ) from four D. melanogaster strains, notably differing in the presence or absence of an intact Nrf2/Maf (a transcription factor) binding site in the 5'-promoter core region, a 'G1410' frameshift deletion mutation in the heme-binding region and a long terminal repeat (LTR) of transposable element 17.6 in the 3'-untranslated region (UTR). Linkage analysis confirmed that DDT resistance was genetically linked to a Nrf2/Maf-binding-site-containing, LTR-lacking functional allele of Cyp6a2 (Cyp6a2(91-R) ). The qRT-PCR results showed that overexpression of functional Cyp6a2 was consistently associated with DDT resistance. Luciferase reporter gene assays revealed that an intact Nrf2/Maf binding site in the 5'-promoter core region enhanced the constitutive transcription of Cyp6a2.

CONCLUSION

The results suggest that the Nrf2/Maf binding-site-containing functional Cyp6a2 allele is associated with DDT resistance in the D. melanogaster strains under study.

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

The Drosophila melanogaster 91-R and 91-C strains are of common origin, however, 91-R has been intensely selected for dichlorodiphenyltrichloroethane (DDT) resistance over six decades while 91-C has been maintained as the non-selected control strain. These fly strains represent a unique genetic resource to understand the accumulation and fixation of mutations under laboratory conditions over decades of pesticide selection. Considerable research has been done to investigate the differential expression of genes associated with the highly DDT resistant strain 91-R, however, with the advent of whole genome sequencing we can now begin to develop an in depth understanding of the genomic changes associated with this intense decades-long xenobiotic selection pressure. Here we present the first whole genome sequencing analysis of the 91-R and 91-C fly strains to identify genome-wide structural changes within the open reading frames. Between-strain changes in allele frequencies revealed a higher percent of new alleles going to fixation for the 91-R strain, as compared to 91-C (P<0.0001). These results suggest that resistance to DDT in the 91-R laboratory strain could potentially be due primarily to new mutations, as well as being polygenic rather than the result of a few major mutations, two hypotheses that remain to be tested.

Differential gene expression in laboratory strains of human head and body lice when challenged with Bartonella quintana, a pathogenic bacterium

Human head and body lice are obligatory hematophagous ectoparasites that belong to a single species, Pediculus humanus. Only body lice, however, are vectors of the infectious Gram-negative bacterium Bartonella quintana. Because of their near identical genomes, yet differential vector competence, head and body lice provide a unique model system to study the gain or loss of vector competence. Using our in vitro louse-rearing system, we infected head and body lice with blood containing B. quintana in order to detect both differences in the proliferation of B. quintana and transcriptional differences of immune-related genes in the lice. B. quintana proliferated rapidly in body lice at 6 days post-infection, but plateaued in head lice at 4 days post-infection. RNAseq and quantitative real-time PCR validation analyses determined gene expression differences. Eight immunoresponse genes were observed to be significantly different with many associated with the Toll pathway: Fibrinogen-like protein, Spaetzle, Defensin 1, Serpin, Scavenger receptor A and Apolipoporhrin 2. Our findings support the hypothesis that body lice, unlike head lice, fight infection from B. quintana only at the later stages of its proliferation.

Genetic differentiation among Maruca vitrata F. (Lepidoptera: Crambidae) populations on cultivated cowpea and wild host plants: implications for insect resistance management and biological control strategies

Maruca vitrata Fabricius (Lepidoptera: Crambidae) is a polyphagous insect pest that feeds on a variety of leguminous plants in the tropics and subtropics. The contribution of host-associated genetic variation on population structure was investigated using analysis of mitochondrial cytochrome oxidase 1 (cox1) sequence and microsatellite marker data from M. vitrata collected from cultivated cowpea (Vigna unguiculata L. Walp.), and alternative host plants Pueraria phaseoloides (Roxb.) Benth. var. javanica (Benth.) Baker, Loncocarpus sericeus (Poir), and Tephrosia candida (Roxb.). Analyses of microsatellite data revealed a significant global FST estimate of 0.05 (P≤0.001). The program STRUCTURE estimated 2 genotypic clusters (co-ancestries) on the four host plants across 3 geographic locations, but little geographic variation was predicted among genotypes from different geographic locations using analysis of molecular variance (AMOVA; among group variation -0.68%) or F-statistics (FSTLoc = -0.01; P = 0.62). These results were corroborated by mitochondrial haplotype data (φSTLoc = 0.05; P = 0.92). In contrast, genotypes obtained from different host plants showed low but significant levels of genetic variation (FSTHost = 0.04; P = 0.01), which accounted for 4.08% of the total genetic variation, but was not congruent with mitochondrial haplotype analyses (φSTHost = 0.06; P = 0.27). Variation among host plants at a location and host plants among locations showed no consistent evidence for M. vitrata population subdivision. These results suggest that host plants do not significantly influence the genetic structure of M. vitrata, and this has implications for biocontrol agent releases as well as insecticide resistance management (IRM) for M. vitrata in West Africa.

Differential effects of RNAi treatments on field populations of the western corn rootworm

RNA interference (RNAi) mediated crop protection against insect pests is a technology that is greatly anticipated by the academic and industrial pest control communities. Prior to commercialization, factors influencing the potential for evolution of insect resistance to RNAi should be evaluated. While mutations in genes encoding the RNAi machinery or the sequences targeted for interference may serve as a prominent mechanism of resistance evolution, differential effects of RNAi on target pests may also facilitate such evolution. However, to date, little is known about how variation of field insect populations could influence the effectiveness of RNAi treatments. To approach this question, we evaluated the effects of RNAi treatments on adults of three western corn rootworm (WCR; Diabrotica virgifera virgifera LeConte) populations exhibiting different levels of gut cysteine protease activity, tolerance of soybean herbivory, and immune gene expression; two populations were collected from crop rotation-resistant (RR) problem areas and one from a location where RR was not observed (wild type; WT). Our results demonstrated that RNAi targeting DvRS5 (a highly expressed cysteine protease gene) reduced gut cysteine protease activity in all three WCR populations. However, the proportion of the cysteine protease activity that was inhibited varied across populations. When WCR adults were treated with double-stranded RNA of an immune gene att1, different changes in survival among WT and RR populations on soybean diets occurred. Notably, for both genes, the sequences targeted for RNAi were the same across all populations examined. These findings indicate that the effectiveness of RNAi treatments could vary among field populations depending on their physiological and genetic backgrounds and that the consistency of an RNAi trait's effectiveness on phenotypically different populations should be considered or tested prior to wide deployment. Also, genes that are potentially subjected to differential selection in the field should be avoided for RNAi-based pest control.

Rates of genomic divergence in humans, chimpanzees and their lice

The rate of DNA mutation and divergence is highly variable across the tree of life. However, the reasons underlying this variation are not well understood. Comparing the rates of genetic changes between hosts and parasite lineages that diverged at the same time is one way to begin to understand differences in genetic mutation and substitution rates. Such studies have indicated that the rate of genetic divergence in parasites is often faster than that of their hosts when comparing single genes. However, the variation in this relative rate of molecular evolution across different genes in the genome is unknown. We compared the rate of DNA sequence divergence between humans, chimpanzees and their ectoparasitic lice for 1534 protein-coding genes across their genomes. The rate of DNA substitution in these orthologous genes was on average 14 times faster for lice than for humans and chimpanzees. In addition, these rates were positively correlated across genes. Because this correlation only occurred for substitutions that changed the amino acid, this pattern is probably produced by similar functional constraints across the same genes in humans, chimpanzees and their ectoparasites.

Insect resistance management for stored product pests: a case study of cowpea weevil (Coleoptera: Bruchidae)

The cowpea weevil, Callosobruchus maculatus F. (Coleoptera: Bruchidae), can cause up to 100% yield loss of stored cowpea seeds in a few months in West Africa. Genes expressing toxins delaying insect maturation (MDTs) are available for genetic engineering. A simulation model was used to investigate the possible use of MDTs for managing C. maculatus. Specifically, we studied the effect of transgenic cowpea expressing an MDT, an insecticide, or both, on the evolution of resistance by C. maculatus at constant temperature. Transgenic cowpea expressing only a nonlethal MDT causing 50-100% maturation delay did not control C. maculatus well. Mortality caused by a maturation delay improved the efficacy of transgenic cowpea expressing only a lethal MDT, but significantly reduced the durability of transgenic cowpea Transgenic cowpea expressing only a lethal MDT causing 50% maturation delay and 90% mortality controlled C. maculatus better than one expressing only a nonlethal MDT, but its durability was only 2 yr. We concluded that transgenic cowpea expressing only an MDT has little value for managing C. maculatus. The resistance by C. maculatus to transgenic cowpea expressing only an insecticide rapidly evolved. Stacking a gene expressing a nonlethal MDT and a gene expressing an insecticide in transgenic cowpea did not significantly improve the durability of an insecticide, but stacking a gene expressing a lethal MDT and a gene expressing an insecticide in transgenic cowpea significantly improved the durability of an insecticide and an MDT. We also discussed this approach within the idea of using transgenic RNAi in pest control strategies.

Development of reference transcriptomes for the major field insect pests of cowpea: a toolbox for insect pest management approaches in west Africa

Cowpea is a widely cultivated and major nutritional source of protein for many people that live in West Africa. Annual yields and longevity of grain storage is greatly reduced by feeding damage caused by a complex of insect pests that include the pod sucking bugs, Anoplocnemis curvipes Fabricius (Hemiptera: Coreidae) and Clavigralla tomentosicollis Stål (Hemiptera: Coreidae); as well as phloem-feeding cowpea aphids, Aphis craccivora Koch (Hemiptera: Aphididae) and flower thrips, Megalurothrips sjostedti Trybom (Thysanoptera: Thripidae). Efforts to control these pests remain a challenge and there is a need to understand the structure and movement of these pest populations in order to facilitate the development of integrated pest management strategies (IPM). Molecular tools have the potential to help facilitate a better understanding of pest populations. Towards this goal, we used 454 pyrosequencing technology to generate 319,126, 176,262, 320,722 and 227,882 raw reads from A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. The reads were de novo assembled into 11,687, 7,647, 10,652 and 7,348 transcripts for A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. Functional annotation of the resulting transcripts identified genes putatively involved in insecticide resistance, pathogen defense and immunity. Additionally, sequences that matched the primary aphid endosymbiont, Buchnera aphidicola, were identified among A. craccivora transcripts. Furthermore, 742, 97, 607 and 180 single nucleotide polymorphisms (SNPs) were respectively predicted among A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti transcripts, and will likely be valuable tools for future molecular genetic marker development. These results demonstrate that Roche 454-based transcriptome sequencing could be useful for the development of genomic resources for cowpea pest insects in West Africa.

Genome-wide analysis of genes associated with moderate and high DDT resistance in Drosophila melanogaster

BACKGROUND

Moderate to high DDT resistance in generally associated with overexpression of multiple genes and therefore has been considered to be polygenic. However, very little information is available about the molecular mechanisms that insect populations employ when evolving increased levels of resistance. The presence of common regulatory motifs among resistance-associated genes may help to explain how and why certain suites of genes are preferentially represented in genomic-scale analyses.

RESULTS

A set of commonly differentially expressed genes associated with DDT resistance in the fruit fly was identified on the basis of genome-wide microarray analysis followed by qRT-PCR verification. More genes were observed to be overtranscribed in the highly resistant strain (91-R) than in the moderately resistant strain (Wisconsin) and susceptible strain (Canton-S). Furthermore, possible transcription factor binding sites that occurred in coexpressed resistance-associated genes were discovered by computational motif discovery methods.

CONCLUSION

A glucocorticoid receptor (GR)-like putative transcription factor binding motif (TFBM) was observed to be associated with genes commonly differentially transcribed in both the 91-R and Wisconsin lines of DDT-resistant Drosophila.

Modeling evolution of resistance by Maruca vitrata (Lepidoptera: Crambidae) to transgenic insecticidal cowpea in Africa

We created a detailed model of the Maruca vitrata (F.) and cowpea [Vigna unguiculata (L.) Walp] system to study the possible evolution of resistance by the insect to transgenic insecticidal cowpea, which is under development. We focused on population dynamics and genetics in a region of west Africa. We simulated single-toxin and pyramided (two-toxin) cowpea and emphasized conservative, worst-case scenarios in our analysis. The results indicate that as long as a pyramided, transgenic cowpea can be developed, seed saving by farmers and reliance on natural refuge are not major problems for resistance management. Furthermore, it is possible that one or both toxins in the pyramid may not need to be high dose for evolution to be delayed significantly (>20 yr or 80 generations for resistance to become a concern if transgenic cowpea is deployed in areas where M. vitrata is endemic). If efforts are made to deploy transgenic cowpea only into the regions where M. vitrata is not endemic, then there is little to no concern with resistance emerging in the M. vitrata population.

The spatial genetic differentiation of the legume pod borer, Maruca vitrata F. (Lepidoptera: Crambidae) populations in West Africa

The legume pod borer, Maruca vitrata, is an endemic insect pest that causes significant yield loss to the cowpea crop in West Africa. The application of population genetic tools is important in the management of insect pests but such data on M. vitrata is lacking. We applied a set of six microsatellite markers to assess the population structure of M. vitrata collected at five sites from Burkina Faso, Niger and Nigeria. Observed polymorphisms ranged from one (marker 3393) to eight (marker 32008) alleles per locus. Observed and expected heterozygosities ranged from 0.0 to 0.8 and 0.0 to 0.6, respectively. Three of the loci in samples from Nigeria and Burkina Faso deviated significantly from Hardy-Weinberg Equilibrium (HWE), whereas no loci deviated significantly in samples from Niger. Analysis of molecular variance (AMOVA) indicated that 67.3% level of the genetic variation was within individuals compared to 17.3% among populations. A global estimate of F ST=0.1 (ENA corrected F ST=0.1) was significant (P⩽0.05) and corroborated by pairwise F ST values that were significant among all possible comparisons. A significant correlation was predicted between genetic divergence and geographic distance between subpopulations (R2=0.6, P=0.04), and cluster analysis by the program STRUCTURE predicted that co-ancestry of genotypes were indicative of three distinct populations. The spatial genetic variance among M. vitrata in West Africa may be due to limited gene flow, south-north seasonal movement pattern or other reproductive barriers. This information is important for the cultural, chemical and biological control strategies for managing M. vitrata.

Evolutionary toxicogenomics: diversification of the Cyp12d1 and Cyp12d3 genes in Drosophila species

Gene duplication and divergence are overwhelmingly considered to be the primary mechanisms by which cytochrome P450 monooxygenases (P450s) have radiated into a large and diverse gene superfamily. To address how environmental stress drives the fixation and diversification of gene duplications, we have analyzed Cyp12d1 and Cyp12d3, a pair of duplicated genes found in the sequenced Drosophila genomes of the melanogaster group. The paralog Cyp12d3, which is not found in Drosophila melanogaster, is basal to the melanogaster group, after it split from the obscura group (ca. 50 mya), and has a significant signature of positive selection in two species (D. sechellia and D. ananassae). Examination of the Cyp12d1 region in D. melanogaster wildtype and isoline populations revealed variation both in copy number and sequence, including splice-site variations, which certainly alter gene function. Further investigations of several strains have identified three cases in which differences in the Cyp12d1 gene region are associated with the differences in transcript abundance and transcriptional responses to the environmental stresses that have not been seen for other detoxificative loci. Together, these data highlight the value of using both macro- and microevolutionary approaches in studying the duplication and divergence events associated with detoxification genes and lay important groundwork for future studies in the field of evolutionary toxicogenomics, which uses the principles of phylogenetic analysis to predict possible enzymatic functions.

Comparison of the transcriptional profiles of head and body lice

Head and body lice are both blood-feeding parasites of humans although only the body louse is a potent disease vector. In spite of numerous morphological and life history differences, head and body lice have recently been hypothesized to be ecotypes of the same species. We took a comparative genomics approach to measure nucleotide diversity by comparing expressed sequence tag data sets from head and body lice. A total of 10 771 body louse and 10 770 head louse transcripts were predicted from a combined assembly of Roche 454 and Illumina sequenced cDNAs from whole body tissues collected at all life stages and during pesticide exposure and bacterial infection treatments. Illumina reads mapped to the 10 775 draft body louse gene models from the whole genome assembly predicted nine presence/absence differences, but PCR confirmation resulted in a single gene difference. Read per million base pair estimates indicated that 14 genes showed significant differential expression between head and body lice under our treatment conditions. One novel microRNA was predicted in both lice species and 99% of the 544 transcripts from Candidatus riesia indicate that they share the same endosymbiont. Overall, few differences exist, which supports the hypothesis that these two organisms are ecotypes of the same species.

Simplify, simplify: Lifestyle and compact genome of the body louse provide a unique functional genomics opportunity

The body louse, with its recently sequenced genome, is now primed to serve as a powerful model organism for addressing fundamental questions relating to how insects interact with their environment. One characteristic of the body louse that facilitates this research is the size of its genome-the smallest insect genome sequenced to date. This diminutive genome must nonetheless control an organism that senses and responds to its environment, reacting to threats of corporal and genomic integrity. Additionally, the body louse transmits several important human diseases compared to its very close relative, the head louse, which does not. Therefore, these two organisms comprise an excellent model system for studying molecular mechanisms associated with vector competence. To understand more fully the development of vector/pathogen interactions, we have developed an in vitro bioassay system and determined that the body louse genome appears to contain the genes necessary for RNAi. The body louse will therefore be useful for determining the set of conditions permissive to the evolution of vector competence.

Transcriptome sequencing, and rapid development and application of SNP markers for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae)

The legume pod borer, Maruca vitrata (Lepidoptera: Crambidae), is an insect pest species of crops grown by subsistence farmers in tropical regions of Africa. We present the de novo assembly of 3729 contigs from 454- and Sanger-derived sequencing reads for midgut, salivary, and whole adult tissues of this non-model species. Functional annotation predicted that 1320 M. vitrata protein coding genes are present, of which 631 have orthologs within the Bombyx mori gene model. A homology-based analysis assigned M. vitrata genes into a group of paralogs, but these were subsequently partitioned into putative orthologs following phylogenetic analyses. Following sequence quality filtering, a total of 1542 putative single nucleotide polymorphisms (SNPs) were predicted within M. vitrata contig assemblies. Seventy one of 1078 designed molecular genetic markers were used to screen M. vitrata samples from five collection sites in West Africa. Population substructure may be present with significant implications in the insect resistance management recommendations pertaining to the release of biological control agents or transgenic cowpea that express Bacillus thuringiensis crystal toxins. Mutation data derived from transcriptome sequencing is an expeditious and economical source for genetic markers that allow evaluation of ecological differentiation.

Systems-scale analysis reveals pathways involved in cellular response to methamphetamine

Background

Methamphetamine (METH), an abused illicit drug, disrupts many cellular processes, including energy metabolism, spermatogenesis, and maintenance of oxidative status. However, many components of the molecular underpinnings of METH toxicity have yet to be established. Network analyses of integrated proteomic, transcriptomic and metabolomic data are particularly well suited for identifying cellular responses to toxins, such as METH, which might otherwise be obscured by the numerous and dynamic changes that are induced.

Methodology/Results

We used network analyses of proteomic and transcriptomic data to evaluate pathways in Drosophila melanogaster that are affected by acute METH toxicity. METH exposure caused changes in the expression of genes involved with energy metabolism, suggesting a Warburg-like effect (aerobic glycolysis), which is normally associated with cancerous cells. Therefore, we tested the hypothesis that carbohydrate metabolism plays an important role in METH toxicity. In agreement with our hypothesis, we observed that increased dietary sugars partially alleviated the toxic effects of METH. Our systems analysis also showed that METH impacted genes and proteins known to be associated with muscular homeostasis/contraction, maintenance of oxidative status, oxidative phosphorylation, spermatogenesis, iron and calcium homeostasis. Our results also provide numerous candidate genes for the METH-induced dysfunction of spermatogenesis, which have not been previously characterized at the molecular level.

Conclusion

Our results support our overall hypothesis that METH causes a toxic syndrome that is characterized by the altered carbohydrate metabolism, dysregulation of calcium and iron homeostasis, increased oxidative stress, and disruption of mitochondrial functions.

Mitochondrial genome sequence and expression profiling for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae)

We report the assembly of the 14,054 bp near complete sequencing of the mitochondrial genome of the legume pod borer (LPB), Maruca vitrata (Lepidoptera: Crambidae), which we subsequently used to estimate divergence and relationships within the lepidopteran lineage. The arrangement and orientation of the 13 protein-coding, 2 rRNA, and 19 tRNA genes sequenced was typical of insect mitochondrial DNA sequences described to date. The sequence contained a high A+T content of 80.1% and a bias for the use of codons with A or T nucleotides in the 3rd position. Transcript mapping with midgut and salivary gland ESTs for mitochondrial genome annotation showed that translation from protein-coding genes initiates and terminates at standard mitochondrial codons, except for the coxI gene, which may start from an arginine CGA codon. The genomic copy of coxII terminates at a T nucleotide, and a proposed polyadenylation mechanism for completion of the TAA stop codon was confirmed by comparisons to EST data. EST contig data further showed that mature M. vitrata mitochondrial transcripts are monocistronic, except for bicistronic transcripts for overlapping genes nd4/nd4L and nd6/cytb, and a tricistronic transcript for atp8/atp6/coxIII. This processing of polycistronic mitochondrial transcripts adheres to the tRNA punctuated cleavage mechanism, whereby mature transcripts are cleaved only at intervening tRNA gene sequences. In contrast, the tricistronic atp8/atp6/coxIII in Drosophila is present as separate atp8/atp6 and coxIII transcripts despite the lack of an intervening tRNA. Our results indicate that mitochondrial processing mechanisms vary between arthropod species, and that it is crucial to use transcriptional information to obtain full annotation of mitochondrial genomes.