Transcription profiling of resistance to Bti toxins in the mosquito Aedes aegypti using next-generation sequencing

Comparative transcriptional analysis between susceptible and resistant populations of Aedes (Stegomyia) aegypti (Linnaeus, 1762) after malathion exposure

Aedes aegypti is an important vector of arboviruses, including dengue, chikungunya and Zika. The application of synthetic insecticides is a frequently used strategy to control this insect. Malathion is an organophosphate insecticide that was widely used in Brazil in the 1980s and 1990s to control the adult form of A. aegypti. In situations where resistance to currently used insecticides is detected, the use of malathion may be resumed as a control measure. Many studies have confirmed resistance to malathion, however, comparative studies of differential gene expression of the entire transcriptome of resistant and susceptible insects are scarce. Therefore, understanding the molecular basis of resistance to this insecticide in this species is extremely important. In this paper, we present the first transcriptomic description of susceptible and resistant strains of A. aegypti challenged with malathion. Guided transcriptome assembly resulted in 39,904 transcripts, where 2133 differentially expressed transcripts were detected, and three were validated by RT-qPCR. Enrichment analysis for these identified transcripts resulted in 13 significant pathways (padj < 0.05), 8 associated with down-regulated and 5 with up-regulated transcripts in treated resistant insects. It was possible to divide the transcripts according to the mechanism of action into three main groups: (i) genes involved in detoxification metabolic pathways; (ii) genes of proteins located in the membrane/extracellular region; and (iii) genes related to DNA integration/function. These results are important in advancing knowledge of genes related to resistance mechanisms in this insect, enabling the development of effective technologies and strategies for managing insecticide resistance.

Discovery and structure-activity relationship of Morita-Baylis-Hillman adducts as larvicides against dengue mosquito vector, Aedes aegypti (Diptera: Culicidae)

Neglected tropical diseases (NTDs) have become a significant public health problem worldwide, notably the life-threatening dengue hemorrhagic fever borne by the Aedes aegypti mosquito. Thus, mosquito vector control measures remain essential in public health vector surveillance and control to combat Aedes-borne infections. Therefore, a series of MBH adducts were synthesized and assessed towards the fourth instar mosquito larvae, Aedes aegypti, along with the preliminary structure-activity relationship (SAR). Noteworthy, this compound class might be synthetized by an efficient eco-friendly synthesismethod and a rapid route for the synthesis of commercial larvicide through a single synthetic step. The bioassays showed that this compound class is a promising larvicide to control Aedes aegypti mosquito larvae, mainly 3g, with an LC₅₀ of 41.35 µg/mL, which was higher than evaluated positive controls. Nevertheless, it is a viable larvicidalhit candidate for further hit-to-leadproperties optimization of its biphenyl backbone scaffold with enhanced insecticidalbioactivity. Moreover, scanning electron microscopy analysis suggested a disruption of the osmoregulatory/ionoregulatory functions by the complete deterioration of the terminal exoskeleton hindgut and anal papillae. Therefore, this new study shows the larvicidal efficacy of the tested compounds against the Aedes aegypti larvae.

Aedes aegypti Strain Subjected to Long-Term Exposure to Bacillus thuringiensis svar. israelensis Larvicides Displays an Altered Transcriptional Response to Zika Virus Infection

Bacillus thuringiensis svar. israelensis (Bti) larvicides are effective in controlling Aedes aegypti; however, the effects of long-term exposure need to be properly evaluated. We established an Ae. aegypti strain that has been treated with Bti for 30 generations (RecBti) and is still susceptible to Bti, but females exhibited increased susceptibility to Zika virus (ZIKV). This study compared the RecBti strain to a reference strain regarding: first, the relative transcription of selected immune genes in ZIKV-challenged females (F30) with increased susceptibility detected in a previous study; then, the whole transcriptomic profile using unchallenged females (F35). Among the genes compared by RT-qPCR in the ZIKV-infected and uninfected females from RecBti (F30) and the reference strain, hop, domeless, relish 1, defensin A, cecropin D, and gambicin showed a trend of repression in RecBti infected females. The transcriptome of RecBti (F35) unchallenged females, compared with a reference strain by RNA-seq, showed a similar profile and only 59 differentially expressed genes were found among 9202 genes analyzed. Our dataset showed that the long-term Bti exposure of the RecBti strain was associated with an alteration of the expression of genes potentially involved in the response to ZIKV infection in challenged females, which is an important feature found under this condition.

The Challenges of Microbial Control of Mosquito-Borne Diseases Due to the Gut Microbiome

Mosquitoes are one of the deadliest animals on earth because of their ability to transmit a wide range of human pathogens. Traditional mosquito control methods use chemical insecticides, but with dwindling long-term effectiveness and negative effects on the environment, microbial forms of control have become common alternatives. The insecticide Bacillus thuringiensis subspecies israelensis (Bti) is the most popular of these alternatives, although it can also have direct effects on lowering environmental biodiversity and indirect effects on food-web relationships in the ecosystems where it is deployed. In addition, microbial control agents that impede pathogen development or transmission from mosquito to human are under investigation, including Wolbachia and Asaia, but unexpected interactions with mosquito gut bacteria can hinder their effectiveness. Improved characterization of mosquito gut bacterial communities is needed to determine the taxa that interfere with microbial controls and their effectiveness in wild populations. This mini-review briefly discusses relationships between mosquito gut bacteria and microbial forms of control, and the challenges in ensuring their success.

Mosquito-Borne Diseases Emergence/Resurgence and How to Effectively Control It Biologically

Deadly pathogens and parasites are transmitted by vectors and the mosquito is considered the most threatening vector in public health, transmitting these pathogens to humans and animals. We are currently witnessing the emergence/resurgence in new regions/populations of the most important mosquito-borne diseases, such as arboviruses and malaria. This resurgence may be the consequence of numerous complex parameters, but the major cause remains the mismanagement of insecticide use and the emergence of resistance. Biological control programmes have rendered promising results but several highly effective techniques, such as genetic manipulation, remain insufficiently considered as a control mechanism. Currently, new strategies based on attractive toxic sugar baits and new agents, such as Wolbachia and Asaia, are being intensively studied for potential use as alternatives to chemicals. Research into new insecticides, Insect Growth Regulators, and repellent compounds is pressing, and the improvement of biological strategies may provide key solutions to prevent outbreaks, decrease the danger to at-risk populations, and mitigate resistance.

Reduced cadherin expression associated with resistance to Bt toxin Cry1Ac in pink bollworm

BACKGROUND

Better understanding of the molecular basis of resistance is needed to improve management of pest resistance to transgenic crops that produce insecticidal proteins from Bacillus thuringiensis (Bt). Here we analyzed resistance of the pink bollworm (Pectinophora gossypiella) to Bt toxin Cry1Ac, which is used widely in transgenic Bt cotton. Field-evolved practical resistance of pink bollworm to Cry1Ac is widespread in India, but not in China or the United States. Previous work with laboratory- and field-selected pink bollworm indicated that resistance to Cry1Ac is caused by changes in the amino acid sequence of a midgut cadherin protein (PgCad1) that binds Cry1Ac in susceptible larvae.

RESULTS

Relative to a susceptible strain, the laboratory-selected APHIS-R strain had 530-fold resistance to Cry1Ac with autosomal recessive inheritance. Unlike previous results, resistance in this strain was not consistently associated with insertions or deletions in the expected amino acid sequence of PgCad1. However, this resistance was associated with 79- to 190-fold reduced transcription of the PgCad1 gene and markedly lower abundance of PgCad1 protein.

CONCLUSION

The ability of pink bollworm and other major pests to evolve resistance to Bt toxins via both qualitative and quantitative changes in receptor proteins demonstrates their remarkable adaptability and presents challenges for monitoring and managing resistance to Bt crops. © 2019 Society of Chemical Industry.

Transcriptome and Proteome Alternation With Resistance to Bacillus thuringiensis Cry1Ah Toxin in Ostrinia furnacalis

Background: Asian corn borer (ACB), Ostrinia furnacalis can develop resistance to transgenic Bacillus thuringiensis (Bt) maize expressing Cry1Ah-toxin. However, the mechanisms that regulate the resistance of ACB to Cry1Ah-toxin are unknown.

Objective: In order to understand the molecular basis of the Cry1Ah-toxin resistance in ACB, “omics” analyses were performed to examine the difference between Cry1Ah-resistant (ACB-AhR) and susceptible (ACB-BtS) strains of ACB at both transcriptional and translational levels.

Results: A total of 7,007 differentially expressed genes (DEGs) and 182 differentially expressed proteins (DEPs) were identified between ACB-AhR and ACB-BtS and 90 genes had simultaneous transcription and translation profiles. Down-regulated genes associated with Cry1Ah resistance included aminopeptidase N, ABCC3, DIMBOA-induced cytochrome P450, alkaline phosphatase, glutathione S-transferase, cadherin-like protein, and V-ATPase. Whereas, anti-stress genes, such as heat shock protein 70 and carboxylesterase were up-regulated in ACB-AhR, displaying that a higher proportion of genes/proteins related to resistance was down-regulated compared to up-regulated. The Kyoto encyclopedia of genes and genomes (KEGG) analysis mapped 578 and 29 DEGs and DEPs, to 27 and 10 pathways, respectively (P < 0.05). Furthermore, real-time quantitative (qRT-PCR) results based on relative expression levels of randomly selected genes confirmed the “omics” response.

Conclusion: Despite the previous studies, this is the first combination of a study using RNA-Seq and iTRAQ approaches on Cry1Ah-toxin binding, which led to the identification of longer length of unigenes in ACB. The DEGs and DEPs results are valuable for further clarifying Cry1Ah-mediated resistance.

Technological Microbiology: Development and Applications

Over thousands of years, modernization could be predicted for the use of microorganisms in the production of foods and beverages. However, the current accelerated pace of new food production is due to the rapid incorporation of biotechnological techniques that allow the rapid identification of new molecules and microorganisms or even the genetic improvement of known species. At no other time in history have microorganisms been so present in areas such as agriculture and medicine, except as recognized villains. Currently, however, beneficial microorganisms such as plant growth promoters and phytopathogen controllers are required by various agricultural crops, and many species are being used as biofactories of important pharmacological molecules. The use of biofactories does not end there: microorganisms have been explored for the synthesis of diverse chemicals, fuel molecules, and industrial polymers, and strains environmentally important due to their biodecomposing or biosorption capacity have gained interest in research laboratories and in industrial activities. We call this new microbiology Technological Microbiology, and we believe that complex techniques, such as heterologous expression and metabolic engineering, can be increasingly incorporated into this applied science, allowing the generation of new and improved products and services.

Transcriptome Sequencing Reveals Large-Scale Changes in Axenic Aedes aegypti Larvae

Mosquitoes host communities of microbes in their digestive tract that consist primarily of bacteria. We previously reported that Aedes aegypti larvae colonized by a native community of bacteria and gnotobiotic larvae colonized by only Escherichia coli develop very similarly into adults, whereas axenic larvae never molt and die as first instars. In this study, we extended these findings by first comparing the growth and abundance of bacteria in conventional, gnotobiotic, and axenic larvae during the first instar. Results showed that conventional and gnotobiotic larvae exhibited no differences in growth, timing of molting, or number of bacteria in their digestive tract. Axenic larvae in contrast grew minimally and never achieved the critical size associated with molting by conventional and gnotobiotic larvae. In the second part of the study we compared patterns of gene expression in conventional, gnotobiotic and axenic larvae by conducting an RNAseq analysis of gut and nongut tissues (carcass) at 22 h post-hatching. Approximately 12% of Ae. aegypti transcripts were differentially expressed in axenic versus conventional or gnotobiotic larvae. However, this profile consisted primarily of transcripts in seven categories that included the down-regulation of select peptidases in the gut and up-regulation of several genes in the gut and carcass with roles in amino acid transport, hormonal signaling, and metabolism. Overall, our results indicate that axenic larvae exhibit alterations in gene expression consistent with defects in acquisition and assimilation of nutrients required for growth.

Transcriptional profiling analysis of Spodoptera litura larvae challenged with Vip3Aa toxin and possible involvement of trypsin in the toxin activation

Vip proteins, a new group of insecticidal toxins produced by Bacillus thuringiensis, are effective against specific pests including Spodoptera litura. Here, we report construction of a transcriptome database of S. litura by de novo assembly along with detection of the transcriptional response of S. litura larvae to Vip3Aa toxin. In total, 56,498 unigenes with an N50 value of 1,853 bp were obtained. Results of transcriptome abundance showed that Vip3Aa toxin provoked a wide transcriptional response of the S. litura midgut. The differentially expressed genes were enriched for immunity-related, metabolic-related and Bt-related genes. Twenty-nine immunity-related genes, 102 metabolic-related genes and 62 Bt-related genes with differential expression were found. On the basis of transcriptional profiling analysis, we focus on the functional validation of trypsin which potentially participated in the activation of Vip3Aa protoxin. Zymogram analysis indicated that the presence of many proteases, including trypsin, in S. litura larvae midgut. Results of enzymolysis in vitro of Vip3Aa by trypsin, and bioassay and histopathology of the trypsin-digested Vip3Aa toxin showed that trypsin was possibly involved in the Vip3Aa activation. This study provides a transcriptome foundation for the identification and functional validation of the differentially expressed genes in an agricultural important pest, S. litura.

The cuticle proteins: a putative role for deltamethrin resistance in Culex pipiens pallens

Insecticide resistance has been a major public health challenge. It is impendent to study the mechanism on insecticide resistance. In our previous study, 14 differentially accumulated insect cuticle proteins (ICPs) based on insecticide resistance proteomes and transcriptomes were found in the deltamethrin-resistant (DR) and -susceptible (DS) strains of Culex pipiens pallens. To investigate if these ICPs are associated with deltamethrin resistance, different transcriptional levels of the 14 ICPs were detected in the DS and DR strains from laboratory and field populations by using quantitative real-time polymerase chain reaction (qRT-PCR). The expression levels of the 14 ICPs were also measured after short-term exposure of the DS strain to deltamethrin. The full-length complementary DNA (cDNA) of CpCPLCG5 gene, which encodes one of the 14 ICPs, was cloned from Cx. pipiens pallens. Homology analysis and phylogenetic analysis were carried out with some other insects. Furthermore, small interfering RNA (siRNA) was used to knockdown the expression level of CpCPLCG5 gene for characterizing its contribution to deltamethrin resistance. The results showed that the expression level of CpCPLCG5 gene was higher in DR strain than in DS strain both in laboratory and field populations while the other 13 ICPs were downregulated. The full-length cDNA of CpCPLCG5 gene was 732 bp, with the ORF of 390 bp and deduced 129 amino acids (GenBank/KF723314,2013). Knockdown of CpCPLCG5 gene increased the susceptibility of the DR strain while the expression level of the other 13 ICPs elevated. Our findings indicate that the cuticle proteins are associated with deltamethrin resistance in Cx. pipiens pallens.

The genetic architecture of a complex trait: Resistance to multiple toxins produced by Bacillus thuringiensis israelensis in the dengue and yellow fever vector, the mosquito Aedes aegypti

The bacterial insecticide Bacillus thuringiensis subsp. israelensis (Bti) is an increasingly popular alternative to chemical insecticides for controlling mosquito populations. Because Bti toxicity relies on the action of four main toxins, resistance to Bti is very likely a complex phenotype involving several genes simultaneously. Dissecting the underlying genetic basis thus requires associating a quantitative measure of resistance to genetic variation at many loci in a segregating population. Here, we undertake this task using the dengue and yellow fever vector, the mosquito Aedes aegypti, as a study model. We conducted QTL (Quantitative Trait Locus) and admixture mapping analyses on two controlled crosses and on an artificial admixed population, respectively, all obtained from resistant and susceptible lab strains. We detected 16 QTL regions, among which four QTLs were revealed by different analysis methods. These four robust QTLs explained altogether 29.2% and 62.2% of the total phenotypic variance in the two QTL crosses, respectively. They also all showed a dominant mode of action. In addition, we found six loci showing statistical association with Bti resistance in the admixed population. Five of the supercontigs highlighted in this study contained candidate genes as suggested by their function, or by prior evidence from expression and/or outlier analyses. These genomic regions are thus good starting points for fine mapping of resistance to Bti or functional analyses aiming at identifying the underlying genes and mutations. Moreover, for the purpose of this work, we built the first Ae. aegypti genetic map based on markers associated with genes expressed in larvae. This genetic map harbors 229 SNP markers mapped across the three chromosomes for a total length of 311.9cM. It brought to light several assembly discrepancies with the reference genome, suggesting a high level of genome plasticity in Ae. aegypti.

Thiosemicarbazones as Aedes aegypti larvicidal

A set of aryl- and phenoxymethyl-(thio)semicarbazones were synthetized, characterized and biologically evaluated against the larvae of Aedes aegypti (A. aegypti), the vector responsible for diseases like Dengue and Yellow Fever. (Q)SAR studies were useful for predicting the activities of the compounds not included to create the QSAR model as well as to predict the features of a new compound with improved activity. Docking studies corroborated experimental evidence of AeSCP-2 as a potential target able to explain the larvicidal properties of its compounds. The trend observed between the in silico Docking scores and the in vitro pLC50 (equals -log LC50, at molar concentration) data indicated that the highest larvicidal compounds, or the compounds with the highest values for pLC50, are usually those with the higher docking scores (i.e., greater in silico affinity for the AeSCP-2 target). Determination of cytotoxicity for these compounds in mammal cells demonstrated that the top larvicide compounds are non-toxic.

Identification of proteins associated with pyrethroid resistance by iTRAQ-based quantitative proteomic analysis in Culex pipiens pallens

Background

Mosquito control based on chemical insecticides is considered as an important element in the current global strategies for the control of mosquito-borne diseases. Unfortunately, the development of pyrethroid resistance in important vector mosquito species jeopardizes the effectiveness of insecticide-based mosquito control. To date, the mechanisms of pyrethroid resistance are still unclear. Recent advances in proteomic techniques can facilitate to identify pyrethroid resistance-associated proteins at a large-scale for improving our understanding of resistance mechanisms, and more importantly, for seeking some genetic markers used for monitoring and predicting the development of resistance.

Methods

We performed a quantitative proteomic analysis between a deltamethrin-susceptible strain and a deltamethrin-resistant strain of laboratory population of Culex pipiens pallens using isobaric tags for relative and absolute quantitation (iTRAQ) analysis. Gene Ontology (GO) analysis was used to find the relative processes that these differentially expressed proteins were involved in. One differentially expressed protein was chosen to confirm by Western blot in the laboratory and field populations of Cx. pipiens pallens.

Results

We identified 30 differentially expressed proteins assigned into 10 different categories, including oxidoreductase activity, transporter activity, catalytic activity, structural constituent of cuticle and hypothetical proteins. GO analysis revealed that 25 proteins were sub-categorized into 35 hierarchically-structured GO classifications. Western blot results showed that CYP6AA9 as one of the up-regulated proteins was confirmed to be overexpressed in the deltamethrin-resistant strains compared with the deltamethrin-susceptible strains both in the laboratory and field populations.

Conclusions

This is the first study to use modern proteomic tools for identifying pyrethroid resistance-related proteins in Cx. pipiens. The present study brought to light many proteins that were not previously thought to be associated with pyrethroid resistance, which further expands our understanding of pyrethroid resistance mechanisms. CYP6AA9 was overexpressed in the deltamethrin-resistant strains, indicating that CYP6AA9 may be involved in pyrethroid resistance and may be used as a potential genetic marker to monitor and predict the pyrethroid resistance level of field populations.

Electronic supplementary material

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

Differential protein expression in the susceptible and resistant Myzus persicae (Sulzer) to imidacloprid

Myzus persicae, a serious economic agricultural pest, has developed resistance to imidacloprid (IMI), which was widely used to control this aphid worldwide. To gain a better understanding of the mechanisms of IMI resistance in M. persicae, we carried out a comparative proteomic analysis. Total proteins of the IMI-susceptible and resistant strains were extracted and separated by two-dimensional gel electrophoresis. More than 1300 protein spots were reproducibly detected, including 14 that were more abundant and 14 less abundant. Mass spectrometry analysis and database searching helped us to identify 25 differentially abundant proteins. The identified proteins were categorized into several functional groups including signal transduction, RNA processing, protein processing, transport processing, stress response, metabolisms, and cytoskeleton structure, etc. This study is the first analysis of differentially expressed proteins in IMI-susceptible and resistant M. Persicae, and gives new insights into the mechanisms of IMI resistance in M. persicae.

Comprehensive analysis of gene expression profiles of the beet armyworm Spodoptera exigua larvae challenged with Bacillus thuringiensis Vip3Aa toxin

Host-pathogen interactions result in complex relationship, many aspects of which are not completely understood. Vip proteins, which are Bacillus thuringensis (Bt) insecticidal toxins produced during the vegetative stage, are selectively effective against specific insect pests. This new group of Bt proteins represents an interesting alternative to the classical Bt Cry toxins because current data suggests that they do not share the same mode of action. We have designed and developed a genome-wide microarray for the beet armyworm Spodoptera exigua, a serious lepidopteran pest of many agricultural crops, and used it to better understand how lepidopteran larvae respond to the treatment with the insecticidal protein Vip3Aa. With this approach, the goal of our study was to evaluate the changes in gene expression levels caused by treatment with sublethal doses of Vip3Aa (causing 99% growth inhibition) at 8 and 24 h after feeding. Results indicated that the toxin provoked a wide transcriptional response, with 19% of the microarray unigenes responding significantly to treatment. The number of up- and down-regulated unigenes was very similar. The number of genes whose expression was regulated at 8 h was similar to the number of genes whose expression was regulated after 24 h of treatment. The up-regulated sequences were enriched for genes involved in innate immune response and in pathogen response such as antimicrobial peptides (AMPs) and repat genes. The down-regulated sequences were mainly unigenes with homology to genes involved in metabolism. Genes related to the mode of action of Bt Cry proteins were found, in general, to be slightly overexpressed. The present study is the first genome-wide analysis of the response of lepidopteran insects to Vip3Aa intoxication. An insight into the molecular mechanisms and components related to Vip intoxication will allow designing of more effective management strategies for pest control.

Gustatory receptor expression in the labella and tarsi of Aedes aegypti

The yellow-fever mosquito, Aedes aegypti, infects a growing number of people every year with dengue, yellow fever and chikungunya viruses. Contact chemoreception in mosquitoes influences a number of behaviors including host-selection, oviposition and feeding. While these behaviors are in many instances well documented, the molecular mechanisms mediating them are not well understood. Here we report the results of sequencing total messenger RNA in the labella and tarsi of both male and female Ae. aegypti to reveal Gustatory Receptor (GR) gene expression profiles in these major gustatory appendages. Gene expression levels in each tissue were verified by RT-qPCR. We discuss potential functions for the GRs revealed here by considering homologous GRs in other insects. Specific GRs provide molecular targets for modification of gustatory-mediated behaviors in this important disease vector.

Insecticide-driven patterns of genetic variation in the dengue vector Aedes aegypti in Martinique Island

Effective vector control is currently challenged worldwide by the evolution of resistance to all classes of chemical insecticides in mosquitoes. In Martinique, populations of the dengue vector Aedes aegypti have been intensively treated with temephos and deltamethrin insecticides over the last fifty years, resulting in heterogeneous levels of resistance across the island. Resistance spreading depends on standing genetic variation, selection intensity and gene flow among populations. To determine gene flow intensity, we first investigated neutral patterns of genetic variability in sixteen populations representative of the many environments found in Martinique and experiencing various levels of insecticide pressure, using 6 microsatellites. Allelic richness was lower in populations resistant to deltamethrin, and consanguinity was higher in populations resistant to temephos, consistent with a negative effect of insecticide pressure on neutral genetic diversity. The global genetic differentiation was low, suggesting high gene flow among populations, but significant structure was found, with a pattern of isolation-by-distance at the global scale. Then, we investigated adaptive patterns of divergence in six out of the 16 populations using 319 single nucleotide polymorphisms (SNPs). Five SNP outliers displaying levels of genetic differentiation out of neutral expectations were detected, including the kdr-V1016I mutation in the voltage-gated sodium channel gene. Association tests revealed a total of seven SNPs associated with deltamethrin resistance. Six other SNPs were associated with temephos resistance, including two non-synonymous substitutions in an alkaline phosphatase and in a sulfotransferase respectively. Altogether, both neutral and adaptive patterns of genetic variation in mosquito populations appear to be largely driven by insecticide pressure in Martinique.

The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector

Mosquitoes are vectors of a number of important human and animal diseases. The development of novel vector control strategies requires a thorough understanding of mosquito biology. To facilitate this, we used RNA-seq to identify novel genes and provide the first high-resolution view of the transcriptome throughout development and in response to blood feeding in a mosquito vector of human disease, Aedes aegypti, the primary vector for Dengue and yellow fever. We characterized mRNA expression at 34 distinct time points throughout Aedes development, including adult somatic and germline tissues, by using polyA+ RNA-seq. We identify a total of 14,238 novel new transcribed regions corresponding to 12,597 new loci, as well as many novel transcript isoforms of previously annotated genes. Altogether these results increase the annotated fraction of the transcribed genome into long polyA+ RNAs by more than twofold. We also identified a number of patterns of shared gene expression, as well as genes and/or exons expressed sex-specifically or sex-differentially. Expression profiles of small RNAs in ovaries, early embryos, testes, and adult male and female somatic tissues also were determined, resulting in the identification of 38 new Aedes-specific miRNAs, and ~291,000 small RNA new transcribed regions, many of which are likely to be endogenous small-interfering RNAs and Piwi-interacting RNAs. Genes of potential interest for transgene-based vector control strategies also are highlighted. Our data have been incorporated into a user-friendly genome browser located at www.Aedes.caltech.edu, with relevant links to Vectorbase (www.vectorbase.org)

Investigating the genetics of Bti resistance using mRNA tag sequencing: application on laboratory strains and natural populations of the dengue vector Aedes aegypti

Mosquito control is often the main method used to reduce mosquito-transmitted diseases. In order to investigate the genetic basis of resistance to the bio-insecticide Bacillus thuringiensis subsp. israelensis (Bti), we used information on polymorphism obtained from cDNA tag sequences from pooled larvae of laboratory Bti-resistant and susceptible Aedes aegypti mosquito strains to identify and analyse 1520 single nucleotide polymorphisms (SNPs). Of the 372 SNPs tested, 99.2% were validated using DNA Illumina GoldenGate® array, with a strong correlation between the allelic frequencies inferred from the pooled and individual data (r = 0.85). A total of 11 genomic regions and five candidate genes were detected using a genome scan approach. One of these candidate genes showed significant departures from neutrality in the resistant strain at sequence level. Six natural populations from Martinique Island were sequenced for the 372 tested SNPs with a high transferability (87%), and association mapping analyses detected 14 loci associated with Bti resistance, including one located in a putative receptor for Cry11 toxins. Three of these loci were also significantly differentiated between the laboratory strains, suggesting that most of the genes associated with resistance might differ between the two environments. It also suggests that common selected regions might harbour key genes for Bti resistance.

Molecular mechanisms associated with increased tolerance to the neonicotinoid insecticide imidacloprid in the dengue vector Aedes aegypti

Mosquitoes are vectors of several major human diseases and their control is mainly based on the use of chemical insecticides. Resistance of mosquitoes to organochlorines, organophosphates, carbamates and pyrethroids led to a regain of interest for the use of neonicotinoid insecticides in vector control. The present study investigated the molecular basis of neonicotinoid resistance in the mosquito Aedes aegypti. A strain susceptible to insecticides was selected at the larval stage with imidacloprid. After eight generations of selection, larvae of the selected strain (Imida-R) showed a 5.4-fold increased tolerance to imidacloprid while adult tolerance level remained low. Imida-R larvae showed significant cross-tolerance to other neonicotinoids but not to pyrethroids, organophosphates and carbamates. Transcriptome profiling identified 344 and 108 genes differentially transcribed in larvae and adults of the Imida-R strain compared to the parental strain. Most of these genes encode detoxification enzymes, cuticle proteins, hexamerins as well as other proteins involved in cell metabolism. Among detoxification enzymes, cytochrome P450 monooxygenases (CYPs) and glucosyl/glucuronosyl transferases (UDPGTs) were over-represented. Bioassays with enzyme inhibitors and biochemical assays confirmed the contribution of P450s with an increased capacity of the Imida-R microsomes to metabolize imidacloprid in presence of NADPH. Comparison of substrate recognition sites and imidacloprid docking models of six CYP6s over-transcribed in the Imida-R strain together with Bemisia tabaci CYP6CM1vQ and Drosophila melanogaster CYP6G1, both able to metabolize imidacloprid, suggested that CYP6BB2 and CYP6N12 are good candidates for imidacloprid metabolism in Ae. aegypti. The present study revealed that imidacloprid tolerance in mosquitoes can arise after few generations of selection at the larval stage but does not lead to a significant tolerance of adults. As in other insects, P450-mediated insecticide metabolism appears to play a major role in imidacloprid tolerance in mosquitoes.

Larval midgut modifications associated with Bti resistance in the yellow fever mosquito using proteomic and transcriptomic approaches

BACKGROUND

Bacillus thuringiensis var. israelensis (Bti) is a natural larval mosquito pathogen producing pore-forming toxins targeting the midgut of Diptera larvae. It is used worldwide for mosquito control. Resistance mechanisms of an Aedes aegypti laboratory strain selected for 30 generations with field-collected leaf litter containing Bti toxins were investigated in larval midguts at two levels: 1. gene transcription using DNA microarray and RT-qPCR and 2. differential expression of brush border membrane proteins using DIGE (Differential In Gel Electrophoresis).

RESULTS

Several Bti Cry toxin receptors including alkaline phosphatases and N-aminopeptidases and toxin-binding V-ATPases exhibited altered expression levels in the resistant strain. The under-expression of putative Bti-receptors is consistent with Bt-resistance mechanisms previously described in Lepidoptera. Four soluble metalloproteinases were found under-transcribed together with a drastic decrease of metalloproteinases activity in the resistant strain, suggesting a role in resistance by decreasing the amount of activated Cry toxins in the larval midgut.

CONCLUSIONS

By combining transcriptomic and proteomic approaches, we detected expression changes at nearly each step of the ingestion-to-infection process, providing a short list of genes and proteins potentially involved in Bti-resistance whose implication needs to be validated. Collectively, these results open the way to further functional analyses to better characterize Bti-resistance mechanisms in mosquitoes.

Transcriptome profiling of the intoxication response of Tenebrio molitor larvae to Bacillus thuringiensis Cry3Aa protoxin

Bacillus thuringiensis (Bt) crystal (Cry) proteins are effective against a select number of insect pests, but improvements are needed to increase efficacy and decrease time to mortality for coleopteran pests. To gain insight into the Bt intoxication process in Coleoptera, we performed RNA-Seq on cDNA generated from the guts of Tenebrio molitor larvae that consumed either a control diet or a diet containing Cry3Aa protoxin. Approximately 134,090 and 124,287 sequence reads from the control and Cry3Aa-treated groups were assembled into 1,318 and 1,140 contigs, respectively. Enrichment analyses indicated that functions associated with mitochondrial respiration, signalling, maintenance of cell structure, membrane integrity, protein recycling/synthesis, and glycosyl hydrolases were significantly increased in Cry3Aa-treated larvae, whereas functions associated with many metabolic processes were reduced, especially glycolysis, tricarboxylic acid cycle, and fatty acid synthesis. Microarray analysis was used to evaluate temporal changes in gene expression after 6, 12 or 24 h of Cry3Aa exposure. Overall, microarray analysis indicated that transcripts related to allergens, chitin-binding proteins, glycosyl hydrolases, and tubulins were induced, and those related to immunity and metabolism were repressed in Cry3Aa-intoxicated larvae. The 24 h microarray data validated most of the RNA-Seq data. Of the three intoxication intervals, larvae demonstrated more differential expression of transcripts after 12 h exposure to Cry3Aa. Gene expression examined by three different methods in control vs. Cry3Aa-treated larvae at the 24 h time point indicated that transcripts encoding proteins with chitin-binding domain 3 were the most differentially expressed in Cry3Aa-intoxicated larvae. Overall, the data suggest that T. molitor larvae mount a complex response to Cry3Aa during the initial 24 h of intoxication. Data from this study represent the largest genetic sequence dataset for T. molitor to date. Furthermore, the methods in this study are useful for comparative analyses in organisms lacking a sequenced genome.