Chemical and biological insecticides select distinct gene expression patterns in Aedes aegypti mosquito

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.

Chromobacterium Biopesticide Exposure Does Not Select for Resistance in Aedes Mosquitoes

Developing effective tools to control mosquito populations is essential for reducing the incidence of diseases like malaria and dengue. Biopesticides of microbial origin are a rich, underexplored source of mosquitocidal compounds. We previously developed a biopesticide from the bacterium Chromobacterium sp. Panama that rapidly kills vector mosquito larvae, including Aedes aegypti and Anopheles gambiae. Here, we demonstrate that two independent Ae. aegypti colonies exposed to a sublethal dose of that biopesticide over consecutive generations persistently exhibited high mortality and developmental delays, indicating that resistance did not develop during the study period. Critically, the descendants of biopesticide-exposed mosquitoes experienced decreased longevity and did not display increased susceptibility to dengue virus or decreased susceptibility to common chemical insecticides. Through RNA sequencing, we observed no link between biopesticide exposure and the increased activity of xenobiotic metabolism and detoxification genes typically associated with insecticide resistance. These findings indicate that the Chromobacterium biopesticide is an exciting, emerging mosquito control tool. IMPORTANCE Vector control is an essential part of mitigating diseases caused by pathogens that mosquitoes spread. Modern vector control is highly reliant on using synthetic insecticides to eliminate mosquito populations before they can cause disease. However, many of these populations have become resistant to commonly used insecticides. There is a strong need to explore alternative vector control strategies that aim to mitigate disease burden. Biopesticides, insecticides of biological origin, can have unique mosquitocidal activities, meaning they can effectively kill mosquitoes that are already resistant to other insecticides. We previously developed a highly effective mosquito biopesticide from the bacterium Chromobacterium sp. Csp_P. Here, we investigate whether exposure to a sublethal dose of this Csp_P biopesticide over 9 to 10 generations causes resistance to arise in Aedes aegypti mosquitoes. We find no evidence of resistance at the physiological or molecular levels, confirming that the Csp_P biopesticide is a highly promising new tool for controlling mosquito populations.

Control of Aedes mosquito populations using recombinant microalgae expressing short hairpin RNAs and their effect on plankton

New biocontrol strategies are urgently needed to combat vector-borne infectious diseases. This study presents a low-cost method to produce a potential mosquito insecticide that utilizes the microalgae released into suburban water sources to control mosquito populations. Chlorella microalgae are ubiquitous in local waters, which were chosen as the host for genetic transfection. This species facilitated the recombinant algae to adapt to the prevailing environmental conditions with rapid growth and high relative abundance. The procedure involved microalgae RNAi-based insecticides developed using short hairpin RNAs targeting the Aedes aegypti chitin synthase A (chsa) gene in Chlorella. These insecticides effectively silenced the chsa gene, inhibiting Aedes metamorphosis in the laboratory and simulated-field trials. This study explored the impact of recombinant microalgae on the phytoplankton and zooplankton in suburban waters. High-throughput sequencing revealed that rapid reproduction of recombinant Chlorella indirectly caused the disappearance of some phytoplankton and reduced the protozoan species. This study demonstrated that a recombinant microalgae-based insecticide could effectively reduce the population of Aedes mosquitoes in the laboratory and simulated field trials. However, the impact of this technology on the environment and ecology requires further investigation.

Effect of marker-free transgenic Chlamydomonas on the control of Aedes mosquito population and on plankton

BACKGROUND

More than half of the world's population suffers from epidemic diseases that are spread by mosquitoes. The primary strategy used to stop the spread of mosquito-borne diseases is vector control. Interference RNA (RNAi) is a powerful tool for controlling insect populations and may be less susceptible to insect resistance than other strategies. However, public concerns have been raised because of the transfer of antibiotic resistance marker genes to environmental microorganisms after integration into the recipient genome, thus allowing the pathogen to acquire resistance. Therefore, in the present study, we modified the 3-hydroxykynurenine transaminase (3hkt) and hormone receptor 3 (hr3) RNAi vectors to remove antibiotic resistance marker genes and retain the expression cassette of the inverse repeat sequence of the 3hkt/hr3 target gene. This recombinant microalgal marker-free RNAi insecticide was subsequently added to the suburban water in a simulated-field trial to test its ability to control mosquito population.

METHODS

The expression cassette of the 3hkt/hr3 inverted repeat sequence and a DNA fragment of the argininosuccinate lyase gene without the ampicillin resistance gene were obtained using restriction enzyme digestion and recovery. After the cotransformation of Chlamydomonas, the recombinant algae was then employed to feed Aedes albopictus larvae. Ten and 300 larvae were used in small- and large-scale laboratory Ae.albopictus feeding trials, respectively. Simulated field trials were conducted using Meishe River water that was complemented with recombinant Chlamydomonas. Moreover, the impact of recombinant microalgae on phytoplankton and zooplankton in the released water was explored via high-throughput sequencing.

RESULTS

The marker-free RNAi-recombinant Chlamydomonas effectively silenced the 3hkt/hr3 target gene, resulting in the inhibition of Ae. albopictus development and also in the high rate of Ae. albopictus larvae mortality in the laboratory and simulated field trials. In addition, the results confirmed that the effect of recombinant Chlamydomonas on plankton in the released water was similar to that of the nontransgenic Chlamydomonas, which could reduce the abundance and species of plankton.

CONCLUSIONS

The marker-free RNAi-recombinant Chlamydomonas are highly lethal to the Ae. albopictus mosquito, and their effect on plankton in released water is similar to that of the nontransgenic algal strains, which reduces the abundance and species of plankton. Thus, marker-free recombinant Chlamydomonas can be used for mosquito biorational control and mosquito-borne disease prevention.

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.

Expansive and Diverse Phenotypic Landscape of Field Aedes aegypti (Diptera: Culicidae) Larvae with Differential Susceptibility to Temephos: Beyond Metabolic Detoxification

Arboviruses including dengue, Zika, and chikungunya are amongst the most significant public health concerns worldwide. Arbovirus control relies on the use of insecticides to control the vector mosquito Aedes aegypti (Linnaeus), the success of which is threatened by widespread insecticide resistance. The work presented here profiled the gene expression of Ae. aegypti larvae from field populations of Ae. aegypti with differential susceptibility to temephos originating from two Colombian urban locations, Bello and Cúcuta, previously reported to have distinctive disease incidence, socioeconomics, and climate. We demonstrated that an exclusive field-to-lab (Ae. aegypti strain New Orleans) comparison generates an over estimation of differential gene expression (DGE) and that the inclusion of a geographically relevant field control yields a more discrete, and likely, more specific set of genes. The composition of the obtained DGE profiles is varied, with commonly reported resistance associated genes including detoxifying enzymes having only a small representation. We identify cuticle biosynthesis, ion exchange homeostasis, an extensive number of long noncoding RNAs, and chromatin modelling among the differentially expressed genes in field resistant Ae. aegypti larvae. It was also shown that temephos resistant larvae undertake further gene expression responses when temporarily exposed to temephos. The results from the sampling triangulation approach here contribute a discrete DGE profiling with reduced noise that permitted the observation of a greater gene diversity, increasing the number of potential targets for the control of insecticide resistant mosquitoes and widening our knowledge base on the complex phenotypic network of the Ae. aegypti response to insecticides.

Aedes aegypti continuously exposed to Bacillus thuringiensis svar. israelensis does not exhibit changes in life traits but displays increased susceptibility for Zika virus

Background

Aedes aegypti can transmit arboviruses worldwide, and Bacillus thuringiensis svar. israelensis (Bti)-based larvicides represent an effective tool for controlling this species. The safety of Bti and lack of resistance have been widely reported; however, little is known regarding the impact of the extensive use of these larvicides on the life traits of mosquitoes. Therefore, this study investigated biological parameters, including susceptibility to arbovirus, of an Ae. aegypti strain (RecBti) subjected to 29 generations of exposure to Bti compared with the RecL reference strain.

Methods

The biological parameters of individuals reared under controlled conditions were compared. Also, the viral susceptibility of females not exposed to Bti during their larval stage was analysed by oral infection and followed until 14 or 21 days post-infection (dpi).

Results

RecBti individuals did not display alterations in the traits that were assessed (fecundity, fertility, pupal weight, developmental time, emergence rate, sex ratio and haematophagic capacity) compared to RecL individuals. Females from both strains were susceptible to dengue serotype 2 (DENV-2) and Zika virus (ZIKV). However, RecBti females showed significantly higher rates of ZIKV infection compared with RecL females at 7 (90% versus 68%, Chi-square: χ² = 7.27, df = 1, P = 0.006) and 14 dpi (100% versus 87%, Chi-square: χ² = 7.69, df = 1, P = 0.005) and for dissemination at 7 dpi (83.3% versus 36%, Fisher’s exact test: P < 0.0001, OR = 0.11, 95% CI 0.03–0.32). Quantification of DENV-2 and ZIKV viral particles produced statistically similar results for females from both strains.

Conclusions

Prolonged exposure of Ae. aegypti larvae to Bti did not alter most of the evaluated biological parameters, except that RecBti females exhibited a higher vector susceptibility for ZIKV. This finding is related to a background of Bti exposure for several generations but not to a previous exposure of the tested females during the larval stage. This study highlights mosquito responses that could be associated with the chronic exposure to Bti in addition to the primary larvicidal effect elicited by this control agent.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04880-6.

Permethrin Resistance in Aedes aegypti Affects Aspects of Vectorial Capacity

Aedes aegypti, as one of the vectors transmitting several arboviruses, is the main target in mosquito control programs. Permethrin is used to control mosquitoes and Aedes aegypti get exposed due to its overuse and are now resistant. The increasing percentage of permethrin resistant Aedes aegypti has become an important issue around the world and the potential influence on vectorial capacity needs to be studied. Here we selected a permethrin resistant (p-s) Aedes aegypti population from a wild Florida population and confirmed the resistance ratio to its parental population. We used allele-specific PCR genotyping of the V1016I and F1534C sites in the sodium channel gene to map mutations responsible for the resistance. Two important factors, survival rate and vector competence, that impact vectorial capacity were checked. Results indicated the p-s population had 20 times more resistance to permethrin based on LD50 compared to the parental population. In the genotyping study, the p-s population had more homozygous mutations in both mutant sites of the sodium channel gene. The p-s adults survived longer and had a higher dissemination rate for dengue virus than the parental population. These results suggest that highly permethrin resistant Aedes aegypti populations might affect the vectorial capacity, moreover, resistance increased the survival time and vector competence, which should be of concern in areas where permethrin is applied.

Comparative analysis of diet-associated responses in two rice planthopper species

Background

Host adaptation is the primary determinant of insect diversification. However, knowledge of different host ranges in closely related species remains scarce. The brown planthopper (Nilaparvata lugens, BPH) and the small brown planthopper (Laodelphax striatellus, SBPH) are the most destructive insect pests within the family Delphacidae. These two species differ in their host range (SBPH can well colonize rice and wheat plants, whereas BPH survives on only rice plants), but the underlying mechanism of this difference remains unknown. High-throughput sequencing provides a powerful approach for analyzing the association between changes in gene expression and the physiological responses of insects. Therefore, gut transcriptomes were performed to elucidate the genes associated with host adaptation in planthoppers. The comparative analysis of planthopper responses to different diets will improve our knowledge of host adaptation regarding herbivorous insects.

Results

In the present study, we analyzed the change in gene expression of SBPHs that were transferred from rice plants to wheat plants over the short term (rSBPH vs tSBPH) or were colonized on wheat plants over the long term (rSBPH vs wSBPH). The results showed that the majority of differentially expressed genes in SBPH showed similar changes in expression for short-term transfer and long-term colonization. Based on a comparative analysis of BPH and SBPH after transfer, the genes associated with sugar transporters and heat-shock proteins showed similar variation. However, most of the genes were differentially regulated between the two species. The detoxification-related genes were upregulated in SBPH after transfer from the rice plants to the wheat plants, but these genes were downregulated in BPH under the same conditions. In contrast, ribosomal-related genes were downregulated in SBPH after transfer, but these genes were upregulated in BPH under the same conditions.

Conclusions

The results of this study provide evidence that host plants played a dominant role in shaping gene expression and that the low fitness of BPH on wheat plants might be determined within 24 h after transfer. This study deepens our understanding of different host ranges for the two planthopper species, which may provide a potential strategy for pest management.

Effects of predator cues and pesticide resistance on the toxicity of a (bio)pesticide mixture

BACKGROUND

Populations of target species are typically exposed to pesticide mixtures and natural stressors such as predator cues, and are increasingly developing resistance to single pesticides. Nevertheless, we have poor knowledge whether natural stressors and the presence of pesticide resistance shape mixture toxicity. We tested the single and combined effects of the pesticide chlorpyrifos and the biopesticide Bacillus thuringiensis israelensis (Bti) on the survival of the Southern house mosquito (Culex quinquefasciatus, Say) and whether these effects were magnified by synthetic predator cues of Notonecta water bugs and differed between a chlorpyrifos-resistant (Ace-1R) and non-resistant (S-Lab) strain.

RESULTS

Single exposure to Bti caused mortality in both strains (S-Lab ∼27%, Ace-1R ∼41%) and single exposure to chlorpyrifos caused only mortality in the S-Lab strain (∼33%), while predator cues did not induce mortality. The chlorpyrifos-resistant strain was 1.5-fold more sensitive to Bti, indicating a cost of resistance. The interaction types between chlorpyrifos and Bti (additive), between chlorpyrifos and predator cues (additive), and between Bti and predator cues (synergistic) were consistent in both strains. Despite predator cues making Bti approximately 8% more lethal, they did not change the additive interaction between Bti and chlorpyrifos in their mixture in either strain.

CONCLUSION

These results indicate that the resistance against chlorpyrifos was not partly lifted when chlorpyrifos exposure was combined with Bti and predator cues. Identifying the interaction type within pesticide mixtures and how this depends on natural stressors is important to select control strategies that give a disadvantage to resistant individuals compared to non-resistant individuals. © 2019 Society of Chemical Industry.

A marker of glutathione S-transferase-mediated resistance to insecticides is associated with higher Plasmodium infection in the African malaria vector Anopheles funestus

Metabolic resistance to insecticides is threatening malaria control in Africa. However, the extent to which it impacts malaria transmission remains unclear. Here, we investigated the association between a marker of glutathione S-transferase mediated metabolic resistance and Plasmodium infection in field population of Anopheles funestus s.s. in comparison to the A296S-RDL target site mutation. The 119F-GSTe2 resistant allele was present in southern (Obout) (56%) and central (Mibellon) (25%) regions of Cameroon whereas the 296S-RDL resistant allele was detected at 98.5% and 15% respectively. The whole mosquito Plasmodium and sporozoite infection rates were 57% and 14.8% respectively in Obout (n = 508) and 19.7% and 5% in Mibellon (n = 360). No association was found between L119F-GSTe2 genotypes and whole mosquito infection status. However, when analyzing oocyst and sporozoite infection rates separately, the resistant homozygote 119F/F genotype was significantly more associated with Plasmodium infection in Obout than both heterozygote (OR = 2.5; P = 0.012) and homozygote susceptible (L/L119) genotypes (OR = 2.10; P = 0.013). In contrast, homozygote RDL susceptible mosquitoes (A/A296) were associated more frequently with Plasmodium infection than other genotypes (OR = 4; P = 0.03). No additive interaction was found between L119F and A296S. Sequencing of the GSTe2 gene showed no association between the polymorphism of this gene and Plasmodium infection. Glutathione S-transferase metabolic resistance is potentially increasing the vectorial capacity of resistant An. funestus mosquitoes. This could result in a possible exacerbation of malaria transmission in areas of high GSTe2-based metabolic resistance to insecticides.

RNAi-based bioinsecticide for Aedes mosquito control

Zika virus infection and dengue and chikungunya fevers are emerging viral diseases that have become public health threats. Their aetiologic agents are transmitted by the bite of genus Aedes mosquitoes. Without effective therapies or vaccines, vector control is the main strategy for preventing the spread of these diseases. Increased insecticide resistance calls for biorational actions focused on control of the target vector population. The chitin required for larval survival structures is a good target for biorational control. Chitin synthases A and B (CHS) are enzymes in the chitin synthesis pathway. Double-stranded RNA (dsRNA)-mediated gene silencing (RNAi) achieves specific knockdown of target proteins. Our goal in this work, a new proposed RNAi-based bioinsecticide, was developed as a potential strategy for mosquito population control. DsRNA molecules that target five different regions in the CHSA and B transcript sequences were produced in vitro and in vivo through expression in E. coli HT115 and tested by direct addition to larval breeding water. Mature and immature larvae treated with dsRNA targeting CHS catalytic sites showed significantly decreased viability associated with a reduction in CHS transcript levels. The few larval and adult survivors displayed an altered morphology and chitin content. In association with diflubenzuron, this bioinsecticide exhibited insecticidal adjuvant properties.

Chitinous Structures as Potential Targets for Insect Pest Control

Chitinous structures are physiologically fundamental in insects. They form the insect exoskeleton, play important roles in physiological systems and provide physical, chemical and biological protections in insects. As critically important structures in insects, chitinous structures are attractive target sites for the development of new insect-pest-control strategies. Chitinous structures in insects are complex and their formation and maintenance are dynamically regulated with the growth and development of insects. In the past few decades, studies on insect chitinous structures have shed lights on the physiological functions, compositions, structural formation, and regulation of the chitinous structures. Current understanding of the chitinous structures has indicated opportunities for exploring new target sites for insect control. Mechanisms to disrupt chitinous structures in insects have been studied and strategies for the potential development of new means of insect control by targeting chitinous structures have been proposed and are practically to be explored.

Long-term exposure of Aedes aegypti to Bacillus thuringiensis svar. israelensis did not involve altered susceptibility to this microbial larvicide or to other control agents

Background

Bacillus thuringiensis svar. israelensis (Bti) is an effective and safe biolarvicide to control Aedes aegypti. Its mode of action based on four protoxins disfavors resistance; however, control in endemic areas that display high mosquito infestation throughout the year requires continuous larvicide applications, which imposes a strong selection pressure. Therefore, this study aimed to investigate the effects of an intensive Bti exposure on an Ae. aegypti strain (RecBti), regarding its susceptibility to Bti and two of its protoxins tested individually, to other control agents temephos and diflubenzuron, and its profile of detoxifying enzymes.

Methods

The RecBti strain was established using a large egg sample (10,000) from Recife city (Brazil) and more than 290,000 larvae were subjected to Bti throughout 30 generations. Larvae susceptibility to larvicides and the activity of detoxifying enzymes were determined by bioassays and catalytic assays, respectively. The Rockefeller strain was the reference used for these evaluations.

Results

Bti exposure yielded an average of 74% mortality at each generation. Larvae assessed in seven time points throughout the 30 generations were susceptible to Bti crystal (resistance ratio RR ≤ 2.8) and to its individual toxins Cry11Aa and Cry4Ba (RR ≤ 4.1). Early signs of altered susceptibility to Cry11Aa were detected in the last evaluations, suggesting that this toxin was a marker of the selection pressure imposed. RecBti larvae were also susceptible (RR ≤ 1.6) to the other control agents, temephos and diflubenzuron. The activity of the detoxifying enzymes α- and β-esterases, glutathione-S-transferases and mixed-function oxidases was classified as unaltered in larvae from two generations (F₁₉ and F₂₅), except for a β-esterases increase in F₂₅.

Conclusions

Prolonged exposure of Ae. aegypti larvae to Bti did not evolve into resistance to the crystal, and no cross-resistance with temephos and diflubenzuron were recorded, which supports their sustainable use with Bti for integrated control practices. The unaltered activity of most detoxifying enzymes suggests that they might not play a major role in the metabolism of Bti toxins, therefore resistance by this mechanism is unlikely to occur. This study also highlights the need to establish suitable criteria to classify the status of larval susceptibility/resistance.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3246-1) contains supplementary material, which is available to authorized users.

Larval Exposure to the Bacterial Insecticide Bti Enhances Dengue Virus Susceptibility of Adult Aedes aegypti Mosquitoes

Understanding the interactions between pathogens sharing the same host can be complicated for holometabolous animals when larval and adult stages are exposed to distinct pathogens. In medically important insect vectors, the effect of pathogen exposure at the larval stage may influence susceptibility to human pathogens at the adult stage. We addressed this hypothesis in the mosquito Aedes aegypti, a major vector of arthropod-borne viruses (arboviruses), such as the dengue virus (DENV) and the chikungunya virus (CHIKV). We experimentally assessed the consequences of sub-lethal exposure to the bacterial pathogen Bacillus thuringiensis subsp. israelensis (Bti), during larval development, on arbovirus susceptibility at the adult stage in three Ae. aegypti strains that differ in their genetic resistance to Bti. We found that larval exposure to Bti significantly increased DENV susceptibility, but not CHIKV susceptibility, in the Bti-resistant strains. However, there was no major difference in the baseline arbovirus susceptibility between the Bti-resistant strains and their Bti-susceptible parental strain. Although the generality of our results remains to be tested with additional arbovirus strains, this study supports the idea that the outcome of an infection by a pathogen depends on other pathogens sharing the same host even when they do not affect the same life stage of the host. Our findings may also have implications for Bti as a mosquito biocontrol agent, indicating that the sub-optimal Bti efficacy may have counter-productive effects by increasing vector competence, at least for some combinations of arbovirus and mosquito strains.

Transcriptomic Analysis of Aedes aegypti in Response to Mosquitocidal Bacillus thuringiensis LLP29 Toxin

Globally, Aedes aegypti is one of the most dangerous mosquitoes that plays a crucial role as a vector for human diseases, such as yellow fever, dengue, and chikungunya. To identify (1) transcriptomic basis of midgut (2) key genes that are involved in the toxicity process by a comparative transcriptomic analysis between the control and Bacillus thuringiensis (Bt) toxin (LLP29 proteins)-treated groups. Next-generation sequencing technology was used to sequence the midgut transcriptome of A. aegypti. A total of 17130 unigenes, including 574 new unigenes, were identified containing 16358 (95.49%) unigenes that were functionally annotated. According to differentially expressed gene (DEG) analysis, 557 DEGs were annotated, including 226 upregulated and 231 downregulated unigenes in the Bt toxin-treated group. A total of 442 DEGs were functionally annotated; among these, 33 were specific to multidrug resistance, 6 were immune-system-related (Lectin, Defensin, Lysozyme), 28 were related to putative proteases, 7 were lipase-related, 8 were related to phosphatases, and 30 were related to other transporters. In addition, the relative expression of 28 DEGs was further confirmed through quantitative real time polymerase chain reaction. The results provide a transcriptomic basis for the identification and functional authentication of DEGs in A. aegypti.

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 cellular responses in midgut tissue of Aedes aegypti larvae following intoxication with Cry11Aa toxin from Bacillus thuringiensis

Background

Although much is known about the mechanism of action of Bacillus thuringiensis Cry toxins, the target tissue cellular responses to toxin activity is less understood. Previous transcriptomic studies indicated that significant changes in gene expression occurred during intoxication. However, most of these studies were done in organisms without a sequenced and annotated reference genome. A reference genome and transcriptome is available for the mosquito Aedes aegypti, and its importance as a disease vector has positioned its biological control as a primary health concern. Through RNA sequencing we sought to determine the transcriptional changes observed during intoxication by Cry11Aa in A. aegypti and to analyze possible defense and recovery mechanisms engaged after toxin ingestion.

Results

In this work the changes in the transcriptome of 4^th instar A. aegypti larvae exposed to Cry11Aa toxin for 0, 3, 6, 9, and 12 h were analyzed. A total of 1060 differentially expressed genes after toxin ingestion were identified with two bioconductoR packages: DESeq2 and EdgeR. The most important transcriptional changes were observed after 9 or 12 h of toxin exposure. GO enrichment analysis of molecular function and biological process were performed as well as Interpro protein functional domains and pBLAST analyses. Up regulated processes include vesicular trafficking, small GTPase signaling, MAPK pathways, and lipid metabolism. In contrast, down regulated functions are related to transmembrane transport, detoxification mechanisms, cell proliferation and metabolism enzymes. Validation with RT-qPCR showed large agreement with Cry11Aa intoxication since these changes were not observed with untreated larvae or larvae treated with non-toxic Cry11Aa mutants, indicating that a fully functional pore forming Cry toxin is required for the observed transcriptional responses.

Conclusions

This study presents the first transcriptome of Cry intoxication response in a fully sequenced insect, and reveals possible conserved cellular processes that enable larvae to contend with Cry intoxication in the disease vector A. aegypti. We found some similarities of the mosquito responses to Cry11Aa toxin with previously observed responses to other Cry toxins in different insect orders and in nematodes suggesting a conserved response to pore forming toxins. Surprisingly some of these responses also correlate with transcriptional changes observed in Bti-resistant and Cry11Aa-resistant mosquito larvae.

Electronic supplementary material

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