Translator Disclaimer
1 May 2010 Differential Gene Expression from Midguts of Refractory and Susceptible Lines of the Mosquito, Aedes aegypti, Infected with Dengue-2 Virus
Author Affiliations +

Suppressive subtractive hybridization was used to evaluate the differential expression of midgut genes of feral populations of Aedes aegypti (Diptera: Culicidae) from Colombia that are naturally refractory or susceptible to Dengue-2 virus infection. A total of 165 differentially expressed sequence tags (ESTs) were identified in the subtracted libraries. The analysis showed a higher number of differentially expressed genes in the susceptible Ae. aegypti individuals than the refractory mosquitoes. The functional annotation of ESTs revealed a broad response in the susceptible library that included immune molecules, metabolic molecules and transcription factors. In the refractory strain, there was the presence of a trypsin inhibitor gene, which could play a role in the infection. These results serve as a template for more detailed studies aiming to characterize the genetic components of refractoriness, which in turn can be used to devise new approaches to combat transmission of dengue fever.


Mosquitoes are the vectors of pathogens that cause important human diseases including malaria, filariasis, dengue, yellow fever among others (Paul et al. 2006). Aedes aegypti (Diptera: Culicidae) is the major vector of dengue viruses that cause more human mortality and morbidity than any other arthropod-transmitted viral disease (Alphey et al. 2002). An estimated 50 to 100 million cases of dengue fever occur annually, and 2.5 billion people are at risk of infection (Black et al. 2002; Guzman and Kouri 2002; Gubler 2004; Guha-Sapir and Schimmer 2005). There is no vaccine available, and there are no drugs to cure dengue fever. Dengue control is based on surveillance measures and mosquito control using insecticides and larval habitat reduction strategies (Spiegel et al. 2005). However, despite these efforts, the number of cases of dengue fever and dengue hemorrhagic fever continue to rise each year, and, therefore, alternative control avenues are being investigated.

Some of these efforts have focused on the genetic manipulation of insect vectors (Beaty 2000; Aksoy et al. 2001; Alphey et al. 2002) to modulate characteristics such as vector competence, the intrinsic ability of a vector to transmit a pathogen (Woodring et al. 1996). Recent advances in molecular biology and the availability of genomic databases have enabled the development of new strategies for the control of vector-borne diseases. Manipulation of vector competence requires extensive knowledge on the molecular aspects of vector-parasite interactions. In this context, transgenic techniques have been used to introduce and achieve expression of foreign, antipathogenic genes in insect vectors (Aksoy et al. 2001; Dotson et al. 2003; Abraham et al. 2005; Riehle and Jacobs-Lorena 2005). In Ae. aegypti, an engineered construct of Sindbis virus has been used to express insect immune peptides (Cheng et al. 2001) or heterologous virus sequences to induce an RNA interference-like response to the target virus (Adelman et al. 2001).

Population genetics studies of vector competence in Ae. aegypti have demonstrated a high variation of this characteristic among different populations (Tabachnick 1982; Apostol et al. 1996; Paupy et al. 2000; Vazeille-Falcoz et al. 2001; Garcia-Franco et al. 2002; Gorrochotegui-Escalante et al. 2002). Susceptible and refractory strains obtained using isofemale selection demonstrated an association of vector competence with genetic components that can be affected by environmental changes (Wallis et al. 1985; Miller and Mitchell 1991). Quantitative genetic studies have revealed that at least two genes or sets of genes control vector competence (Bosio et al. 2000; Black et al. 2002). These studies have allowed associating a genetic component with vector competence, but the role of specific receptors or immune response-related genes that modulate arbovirus infection and replication, and the factors that determine resistance or susceptibility to arboviruses such as dengue, are still unknown.

Insects have developed precise mechanisms to protect themselves against bacterial, fungal, and parasitic infections. This immune response is innate and, depending on the type and size of the pathogen, may involve a combination of phagocytosis (Pearson et al. 1995; Kocks et al. 2005; Watson et al. 2005; Lemaitre and Hoffmann 2007), encapsulation and melanization (Karlsson et al. 2004; Bidla et al. 2005; Paskewitz et al. 2006), and production of lethal antimicrobial peptides (Lowenberger 2001; Imler and Bulet 2005). Immune reactions begin with the recognition of cell-surface molecules of pathogens, known as pathogen-associated molecular patterns, by specific insect receptors (pattern recognition receptors) (Michel and Kafatos 2005). This interaction between pathogen-associated molecular patterns and pattern recognition receptors selectively activates either of two intracellular signaling pathways. In Aedes, as in Drosophila, Gram + bacterial and fungal infections induce the Toll pathway, which results in the translocation of the NF-KB transcription factor, and Gram — bacteria triggers the Imd pathway, which results in the nuclear translocation of Relish (a NF-κ;B-type transcription factor) and the induction of antimicrobial peptides such as cecropin and defensin (Lowenberger 2001; Bartholomay et al. 2004). Although these immune pathways are conserved among mosquito species, there are differences in the molecules involved (Hoffmann and Reichhart 2002; Shin et al. 2003; Meister et al. 2005). Despite the knowledge of antimicrobial responses, the immune response in mosquitoes against viruses such as dengue has not been thoroughly defined (Sanders et al. 2005). Other insect-virus models such as Drosophila/Drosophila virus C suggest that, in addition to Toll and Imd, a third pathway, Jak/Stat, functions as a part of an antiviral innate immune response (Dostert et al. 2005; Zambon et al. 2005). In vertebrates, apoptosis is the first response to viral infections and later stimulates the adaptive immune system. This mechanism has also been described in some insect models infected with baculovirus as an antiviral defense strategy, but the role of apoptosis as an immune response has not been well-characterized in insects that transmit human parasites (Clarke and Clem 2003; Cooper et al. 2007a), despite the fact that apoptotic-like activity has been associated with Plasmodium infection in Anopheles gambiae (Al-Olayan et al. 2002). Recent studies in Ae. aegypti have identified and characterized two initiator caspases associated with apoptosis pathways, suggesting that this immune response might function as one of the mechanisms that insect vectors use to regulate the establishment and replication of intracellular parasites such as viruses (Cooper et al. 2007a, 2007b).

Although Ae. aegypti is the main vector of dengue virus, there are populations or fractions of populations that do not permit virus development, presumably because they have biological barriers that impair the establishment and dissemination process (Black et al. 2002). A high variability in vector competence among local populations of Ae. aegypti in Cali, Colombia ranging between 19% to 60% was identified (Ocampo and Wesson 2004), indicating the presence of naturally susceptible and refractory mosquitoes to Dengue-2 virus with different infection barriers. The midgut infection barrier is one of the initial mechanisms that viruses must overcome to establish a successful infection and is one that might be genetically altered to render mosquitoes resistant to arboviruses. Therefore, it is the focus of this study.

The differential expression of midgut genes between susceptible and refractory Ae. aegypti after exposure to Dengue-2 virus is reported here. The functional annotation of pathogen-specific vector-expressed sequence tags (ESTs) that could play a role in determining or contributing to vector competence in Ae. aegypti are also described.

Materials and Methods

Study rationale

In an attempt to identify factors that determine the susceptibility of Ae. aegypti to dengue virus, midgut gene expression was evaluated in Ae. aegypti individuals that were susceptible or refractory to Dengue-2 virus 48 h post infection.

Mosquito strains

Ae. aegypti were collected in different localities from the city of Cali, Colombia and colonized at the Centro Internacional de Entrenamiento e Investigaciones Médicas insectary at 26 ± 2° C with 80% relative humidity and a 12:12 lightdark photoperiod. Ae. aegypti (Rockefeller strain) provided by the Centers for Disease Control and Prevention (Puerto Rico) were maintained in the insectary and were used as positive and negative controls in all PCR assays. The variability of susceptibility of different local populations of mosquitoes to Dengue-2 was described previously (Ocampo and Wesson 2004), indicating the presence of susceptible and refractory mosquitoes in Cali. To increase the number of refractory mosquitoes, an isofemale selection was carried out. For this selection, females were allowed to feed on an infectious blood meal and were placed individually in oviposition cages. After 14 days of incubation, the phenotype of the mother, with respect to the biological barriers, was identified (midgut infection barrier, midgut escape barrier, and susceptibility) as described by Bennett et al. (2005). The eggs of each female were collected. These eggs were hatched and the emerging adults were sorted based on the phenotype of the mother. These descendents from susceptible and midgut infection barrier females were infected with a Dengue-2 infectious blood meal and midgut tissues were collected 48 h later.

Virus maintenance and mosquito infection

Dengue-2 virus New Guinea C strain, freshly grown in C6/36HT (Aedes albopictus larvae cells) was used in oral challenges. Infected cells were incubated for 14 days at 32° C in L 15 medium supplemented with 2% heatinactivated fetal bovine serum, 1% penicillin/streptomycin, and 1% L-glutamine (Higgs and Beaty 1996). Virus and cells were harvested and collected in a 15-ml conical centrifuge tube. Aliquots of the infected cell suspension and the mixture of blood and virus before and after the infection process were titred using the methodology described by Bennett et al. (2002). Titres in the cell suspensions ranged from 108 to 108.5 TCID50/ml in all oral challenges. Oral infections were done in a BSL2+ (biosafety laboratory) insectary with eight protection barriers. Artificial blood feeding was carried out using a membrane feeder. Infected blood was prepared by mixing defibrinated rabbit blood and Dengue-2 virus suspension (1:1 v/v) (Higgs and Beaty 1996). Adult females, six to seven days after eclosion, were deprived of sucrose and water for 24 h prior to blood feeding. Mosquitoes were allowed 1 to 1.5 hours to feed ad libitum. Fully engorged mosquitoes were separated and kept in a separate cage with access to a 10% sugar solution.

Tissue dissection and RNA isolation

Forty-eight hours after infection, midguts from bloodied mosquitoes were dissected on a chilled table and thoroughly rinsed in cold DEPC-PBS to remove traces of the blood meal. Tissues were stored individually in RNA later (Qiagen, at -20° C for subsequent RNA isolation. The RNA later solution was removed by pipetting. Total RNA extraction from individual midguts was performed using RNeasy Mini Kit (Qiagen) according to the manufacturer's instructions. Total RNA was quantified using a NanoDrop Spectrophotometer ND-1000 (NanoDrop Technologies,

Detection of infection

An established nested reverse transcriptase PCR protocol (Lanciotti et al. 1992) was standardized using three groups of Ae. aegypti Rockefeller strain to determine the sensitivity of reverse transcriptase PCR to detect the virus in individual midguts. These groups were: mosquitoes inoculated with Dengue-2 virus (positive controls), infected-bloodfed mosquitoes, and non-bloodfed (naïve) mosquitoes. Midguts were dissected, and RNA was extracted as described above.

In the reverse transcriptase PCR reactions, 50 ng of total RNA were reverse transcribed in a 20 µl reaction mixture containing 1X firststrand buffer (50 mM Tris-HCl (pH 8.3), 75 mM KCl, 3 mM MgCl2), 5 mM DTT, 500 µM of dNTPs mix, 50 pmol of primer D2 (5′TTGCACCAACAGTCAATGTCTTCAGGT TC-3′) and 50 units of Superscript II Reverse Transcriptase (Invitrogen, Reverse transcription was conducted at 42° C for 60 min and 95° C for 5 min. The resulting cDNA was used in a 50 ìl PCR reaction containing 1X PCR buffer (50 mM KCl, 10 mM Tris-HCl (pH 9.0), 0.1% Triton® X-100), 1.5 mM MgCl2, 125 µM of each dNTP, 50 pmol of primers D1 (5′TCAATATGCTGAAACGCGCGAGAAACC G-3′) and D2, and 0.05 U of Taq DNA polymerase (Invitrogen). PCR was performed with the following parameters: 95° C for 1 min; 30 cycles of 94° C for 45 s, 58° C for 45 s, and 72° C for 1 min; and a final extension at 72° C for 7 min. A second-round PCR was run with a 1:100 dilution from the first PCR reaction. PCR was performed under the same conditions used for the primary PCR with the following modifications: primer D2 was replaced with the Dengue-2 virus-specific primer TS2 (5′-CGCCACAAGGGCCATGA ACAG-3′, 50 pmol) and 35 amplification cycles were used. PCR products were resolved by 2% agarose gel electrophoresis with a 100bp DNA ladder (Invitrogen) stained with ethidium bromide and visualized under UV light.

Subtractive library construction

According to the PCR result, positive (infected) and negative (non-infected) midgut RNA samples for each phenotype were pooled separately. A total of 60 midguts were pooled for each phenotype to obtain sufficient RNA to generate the suppressive subtractive hybridization (SSH) libraries. All RNA pools were precipitated and treated with DNAse (Qiagen).

Total RNA from each pool was used to generate cDNA using the SMART PCR cDNA Synthesis Kit (Clontech, according to the manufacturer's recommendations. This procedure generated a sufficient quantity of high-quality cDNA from small quantities of RNA for subtractive library procedures.

Libraries were built using PCR-Select cDNA Subtraction kit (Clontech) according to the manufacturer's specifications. SSH is a PCR-based technique that facilitates the detection of differentially expressed sequences in two samples by allowing exponential amplification of differentially expressed genes and suppressing the amplification of sequences common to both samples. This technique has been used previously to identify differentially expressed genes in Rhodnius prolixus in response to pathogens and parasites (Ursic-Bedoya and Lowenberger 2007).

Three subtractive libraries were constructed: a library of differentially expressed genes in mosquito midguts after injection of Escherichia coli (control library), and two cDNA libraries from the midguts of Dengue-2 virus-susceptible and virus-refractory mosquitoes after the ingestion of a Dengues-infected blood meal. A predictable E. coli control library was built to confirm that the small amount of RNA available for the dengue-susceptible and dengue-refractory was sufficient to build SSH libraries.

The products of the subtracted procedure were ligated into pGemT Easy plasmid vector (Promega, and transformed by heat shock into E. coli JM109 ultra-competent cells (Promega) as previously described by Ursic-Bedoya and Lowenberger (2007). Putative transformant colonies were grown overnight in 96-well plates with 100 µl of LB medium and 0.1 µl of ampicillin (100 µg/µl). For forward and reverse libraries, a total of 384 colonies (four plates) were selected from each library for differential screening. For the bacteria-induced library, 192 colonies (2 plates) were selected.

Subtraction efficiency analysis and differential screening

The subtraction efficiency of the SSH process in all libraries was measured using PCR to amplify, before and after subtraction, a housekeeping gene that should be present in both libraries and an induced gene that should be present in only the enriched library. The Beta-actin sequence from Ae. aegypti with the primers actinF637LVP: 5′- ATTAAGGAGAAGCTGTGCTACGTC and actinR942LVP: 5′-CATACGATCAGCA TTACCTGGG was used. The PCR program was 94° C for 1 min, followed by 33 cycles of 94° C for 20 s, 60° C for 20 s, 68° C for 30 s and a final extension of 68° C for 2 min. To measure a differentially expressed gene, the Ae. aegypti lysozyme was used as described by Ursic-Bedoya and Lowenberger (2007).

Both midgut subtracted libraries were screened for differentially expressed ESTs using the PCR-select differential screening kit (Clontech) following the manufacturer's instructions. One hundred and fifty nanograms from the forward and reverse subtracted libraries were used to create a 32P-labeled probe by random priming. Forward and reverse subtracted probes were hybridized in individual tubes with Hybond+ DNA membranes (Amersham Biosciences, containing individually spotted EST clones (Ursic-Bedoya and Lowenberger 2007). 32P-labeled probes and target EST membranes were hybridized at 65° C for 2.5 h in a rotatory oven using Rapid—Hyb buffer (Amersham Biosciences). Following hybridization, the membranes were washed with low stringency (2X SSC, 0.5% SDS; 3 times, 20 min each) and high stringency (0.2X SSC, 0.5% SDS; 3 times, 20 min each) buffers at 65° C to eliminate non-specific binding due to excess probe. Membranes were exposed to a Kodak BioMax MS film (Eastman Kodak, overnight at room temperature. Selected colonies (strong signal with the forward and low signal with the reverse subtracted probe) were sent to BC Genome Sciences Centre (Vancouver, BC) for plasmid purification and sequencing.

Sequence analysis

Sequence homology searches were carried out using NCBI's BLAST-X ( against nr databases with default parameters. The best annotated matches were retained. Sequences with no significant matches in NCBI's databases were translated in all possible reading frames and were analyzed using INTERPRO SCAN to identify conserved protein domains so that putative function could be assigned. Additionally, sequences with no significant match in the NCBI program were analyzed against the VECTORBASE database ( Homologies were considered statistically significant if they generated an Expect value (E) < 0.1. The EST sequences reported in this paper were submitted to the NCBI dbEST and assigned accession numbers 56768811 to 567689975 (GenBank FG107129 to FG107293).


Library of ESTs found only in Dengue-2 susceptible midguts

This library was created using susceptible insects as the tester and refractory insects as the driver (forward library) in order to identify genes differentially expressed in the susceptible population. In order to increase selection of susceptibility-related genes, recombinant colonies were differentially screened by hybridization with forward (susceptible) and reverse (refractory) probes. Differentially expressed and over-expressed clones were selected. Of the 384 clones screened, only 125 were confirmed to be up-regulated by differential screening as described above.

All 125 clones were sequenced, from which 22 clones (17.6%) did not have similarities with other sequences in the databases (data not shown). A similarity search identified 57 putative genes from 103 clones that matched with annotated sequences in databases (Tables 1, 2). All identified genes were clustered in functional groups according to their putative function as cytoskeleton, nucleic acid binding, metabolism, transcription factors, immunity, ion binding and transport, receptors, mitochondrial, signaling and digestion genes. Out of 57 ESTs, 9 sequences coded for ribosomal genes (normally repressed in the suppressive subtractive hybridization), and 15 clones corresponded to hypothetical proteins. Only 11 of the putative genes had more than one copy, and 4 of them were highly repetitive (more than 3 copies) as DNA binding, Cytoskeleton, Cytochrome P450 and calcium ion binding genes that are potentially related with intracellular infections, cellular distress, and immune responses.

Library of ESTs found only in Dengue-2 refractory midguts

A total of 384 clones were spotted on membranes, but only 40 were confirmed to be up-regulated in the refractory tissues compared with the susceptible library after hybridization with the forward and reverse probes. Of those 40 clones, 5 had no significant match to other genes in the databases (data not shown). Bioinformatic analyses showed 23 (65.7%) of the EST sequences corresponded to different putative genes (Tables 3, 4). Among these transcripts, 9 sequences were hypothetical proteins. Four clones (1%) had more than one copy, but, in contrast with the susceptible library, they were not highly repetitive. Subtraction efficiency analysis by PCR showed better quality in this subtraction than in the susceptible library since only one ribosomal gene was detected.

Contrary to the susceptible subtracted library, genes related with cellular stress or immune responses were not detected (Figure 1). An interesting finding was the presence of a trypsin inhibitor gene that was differentially expressed in this library (Table 3). This protein could affect dengue virus infection; the inhibition or knockdown of specific trypsin molecules has been reported to reduce (Molina-Cruz et al. 2005) or increase (Brackney et al. 2008) the infectivity of dengue virus in Ae. aegypti.

Table 1.

EST identified in the susceptible midgut subtracted library using BLAST database.




Figure 1.

Differential gene expression in susceptible and refractory subtracted libraries. High quality figures are available online.


Table 2.

EST identified in the susceptible midgut subtracted library using VECTORBASE database.







Insect innate immune responses have been studied against bacteria, parasites, and fungi, but antiviral responses have not been well-characterized (Sanders et al. 2005). In eukaryotic organisms, antiviral innate immune mechanisms involve mechanical barriers (Paskewitz and Christensen 1996; Schmid-Hempel 2001; Schmid-Hempel 2005), gene silencing (RNAi and miRNA) (Ausubel 2005; Fritz et al. 2006; Wang et al. 2006), production of humoral and effector mechanisms (Cherry and Silverman 2006; Seth et al. 2006; Zhong et al. 2006) and apoptosis (Everett and McFadden 1999; Irusta et al. 2003).

It is known that the mosquito immune response mechanisms are similar to those of Drosophila melanogaster. However, D. melanogaster is not a vector of arboviruses, and therefore, certain immune mechanisms could be specific to mosquitoes such as Ae. aegypti. Although there are conserved genes among D. melanogaster, An. gambiae and Ae. aegypti, immune genes are the most divergent group even among closely related species. An estimated 285 genes related to immune response, apoptosis and oxidative stress were identified in D. melanogaster, while 338 such genes were identified in An. gambiae, and 353 in Ae. aegypti (Nene et al. 2007). Phylogenetic studies indicate that genes related to pathogen recognition and signaling intracellular pathways are conserved in the three species. Genes encoding effector molecules such as antimicrobial peptides, however, may be more diverse or species specific (Nene et al. 2007; Waterhouse et al. 2007).

Table 3.

EST identified in the refractory midgut subtracted library using BLAST database.


Table 4.

EST identified in the refractory midgut subtracted library using VECTORBASE database.


Viral infection studies in D. melanogaster have revealed the role of RNA interference, Toll, Imd, Jak-Stat pathways and apoptosis mechanisms as antiviral responses. The diverse kinds of responses seen in these experiments depend on the virus life cycle and the method used in the experimental infections (oral or intrathoracic inoculation) (Keene et al. 2004; Lemaitre and Hoffmann 2007). This study focused on identifying immune response related genes in the mosquito midgut, the primary barrier that the virus must overcome in order to infect the vector. It has been estimated that only 40–50% of the initial viral load is capable of crossing the midgut barrier (Paskewitz and Christensen 1996).

In this study, differential expression analysis showed that susceptible Ae. aegypti express a higher number of metabolic and immune response-related genes than the refractory mosquitoes. Immune genes were primarily associated with both inhibition and execution mechanisms of apoptosis. The identified molecules associated with this process were: inhibitor of apoptosis, ubitiquin ligase complex, Cytochrome c, cytoskeleton genes and proteins with calcium ion binding domain. The over-expression of these molecules suggests that cell stress and apoptosis could be playing a key role during the infection process. Other studies have hypothesized that mosquitoes might use apoptosis to eliminate intracellular parasites such as viruses (Cooper et al. 2007a, 2007b), and that viruses might prevent it by activating the expression of inhibitors of apoptosis. The presence of anti and pro-apoptotic molecules suggests that the virus manipulates the cellular machinery to allow its replication and dissemination, as has been demonstrated in cell cultures (O'Brien 1998). To this point, the differential apoptotic responses in refractory and susceptible Ae. aegypti strains after ingesting Dengue-2 were not measured, but experiments are ongoing.

In contrast, an upregulation of immune-related genes was not observed in refractory mosquitoes possibly because: a) the virus does not enter the midgut cells, b) the virus cannot establish the infection, c) there are early events of apoptosis that eliminate virus-infected cells, d) the technique used or availability of all genes in the databases did not allow the identification of other potential refractory mechanisms, or e) the mechanism is not mediated at the transcriptional level, but may be mediated by previously produced zymogens. In this context, an interesting finding in the refractory strain was the presence of a trypsin inhibitor gene. Trypsin is a digestive enzyme that has been implicated in the dengue virus infectious process (Molina-Cruz et al. 2005, Brackney et al. 2008). It is possible the trypsin inhibitor gene found in this library could affect dengue infection in the mosquito. If it is supported, this finding may suggest that refractoriness may not be due to an active immune response based on well-known and well-characterized immune processes (i.e. antimicrobial peptide expression, phagocytosis activation by TEP, phenoloxidase, melanization) but may be a function of the virus not being able to exit midgut epithelial cells, if it was able to enter initially.

The SSH technique and its ability to identify differentially expressed genes in the midguts of Dengue-susceptible and -refractory individuals were validated, but there are few similar studies with which to compare the results. The E. coli-injected control library identified some immune peptides that have been demonstrated broadly in other studies as cecropin, serine proteases, and conserved protein related to cell death (data not shown).

The limitations of the SSH technique were recognized in terms of its sensitivity, since it detects only highly over-expressed genes; however, this technique was selected for its accessibility and as a primary step in identifying potential differences between the susceptible and refractory strains. Additionally, this technique allowed us to work with the small amount of RNA that was obtained by pooling individual midguts that were previously tested.

As described above, many of the ESTs generated in this study have no known match in the databases, and they will continue to be submitted to the growing number of databases as more motifs and genomes are sequenced. The lack of match may indicate a true lack of comparable sequences in the databases, or may indicate that the ESTs map to 3′ and 5′ untranslated regions. The ESTs were examined, but none of the classic motifs found in 3′ untranslated regions were found.

Some of the more interesting proteins to which the ESTs map, and which could play a key role in the susceptibility or refractoriness to Dengue-2 virus in Ae. aegypti, are discussed below. Further studies are underway to evaluate these molecules in more detail.

Inhibitor of apoptosis

Insect inhibitor-of-apoptosis proteins contain two baculoviral inhibitor-of-apoptosis repeat domains and a Zinc RING domain. Inhibitors of apoptosis impede activation of initiator and executioner caspases preventing either their dimerization or their binding to the active catalytic site of these enzymes (Huh et al. 2007; Leu et al. 2007). Some inhibitors of apoptosis have been identified and characterized in insects, but their significance during arbovirus infection in mosquitoes has not been completely elucidated (Blitvich et al. 2002; Li et al. 2007; Bryant et al. 2008). In an An. gambiae functional genomics study, 6 inhibitors of apoptosis were differentially expressed during Plasmodium berghei infection in midgut epithelial cells (Vlachou et al. 2005). Likewise, there was an upregulation of apoptosis related-molecules in Ae. aegypti infected with Sindbis virus, and, among these, one inhibitor of apoptosis was over-expressed (Sanders et al. 2005). What is most interesting in this study is, for the first time, the identification of apoptosis as an antiviral response in a natural Dengue/Ae. aegypti model using wild mosquito populations.

Apoptosis as an immune response mechanism in vertebrates has been widely described. However, in invertebrates, this process has not been clearly characterized. The results of this study are consistent with other studies, suggesting this cell death process is one of the mechanisms that insect vectors use to regulate intracellular parasites such as viruses (Cooper et al. 2007a, 2007b).

Ubiquitin ligase complex

This enzyme complex participates in protein degradation by the proteasome in a number of key biological processes, including cell cycle progression and signal transduction (Maniatis 1999). Ubiquitin-dependent proteolysis controls the abundance of many regulatory proteins and caspase activation (Wojcik 2002; Arama et al. 2007). Several studies have shown a crosstalk between the apoptotic pathways and the ubiquitin- proteasome system (Orlowski 1999; Schreader et al. 2003; Arama et al. 2007). During the cell death process, the ubiquitin ligase complex promotes caspase activation via ubiquitination and degradation of caspase inhibitors. Inhibitors of apoptosis, with ubiquitin protease ligase (E3) activity in their RING finger domain, undergo auto-ubiquitination and degradation by proteasome (Grimm and Osborne 1999; Hu and Yang 2003). The ubiquitin-mediated pathway also regulates NF-KB factors in activation. In the D. melanogaster Toll pathway, upregulation of ubiquitin ligase levels leads to the degradation of Cactus, allowing the nuclear translocation of Dorsal (Spencer et al. 1999).

Cytochromes (P450 and c)

Cytochromes are proteins involved in several cellular functions such as oxidative stress, respiration, apoptosis and xenobiotic metabolism (Scott and Kasai 2004; Arama et al. 2006). In mammal cells, release of Cytochrome c and other proapoptotic molecules induce caspase activation and cell death via the mitochondrial apoptosis pathway (Hengartner 2000; Wang 2001). In insects, the role of mitochondria and Cytochrome c in apoptosis has been contradictory and not completely characterized (Abdelwahid et al. 2007; Goyal et al. 2007). Some reports have suggested that Cytochrome c release is not a necessary step to trigger apoptosis in some D. melanogaster cells (Dorstyn et al. 2004; Abdelwahid et al. 2007). In contrast, Liu (2007) showed that Cytochrome c is required for caspase activation in Baculovirus-induced apoptosis in Spodoptera litura cells.

Trypsin inhibitor

Trypsin, a major midgut proteolytic enzyme, is essential for blood digestion in Ae. aegypti (Noriega and Wells 1999) and its presence has been reported to enhance arbovirus infectivity, replication, and dissemination (Ludwig et al. 1991; Xu et al. 1997, Molina-Cruz et al. 2005). The addition of soybean trypsin inhibitor has been reported either to increase midgut infection rates (Brackney et al. 2008) or to decrease Dengue-2 infectivity and dissemination (Molina-Cruz et al. 2005).

Feeding different trypsin inhibitors blocked early trypsin activity but did not reduce late trypsin expression, and RNAi knockdown of early trypsin did not affect late trypsin expression (Lu et al. 2006). The interaction between the expression of early and late trypsin is not clear. RNAi silencing of chymotrypsin, early and late trypsin had no effect on Dengue-2 infectivity whereas RNAi knockdown of a third trypsin, 5G1, reduced trypsin activity and increased dengue infectivity in the midgut (Brackney et al. 2008). These studies suggest that some midgut serine proteases, acting through digestion or direct activity on viral proteins, may affect Dengue-2 infectivity of Ae. aegypti. The presence of the EST with high homology to a trypsin inhibitor (Table 3) in the refractory strain suggests that inhibition of trypsin activity as a digestive enzyme or in cleaving viral proteins could contribute to the refractoriness of the wild Ae. aegypti population used in this study. However, further studies are required to determine which trypsins are affected by this inhibitor and subsequently their specific roles in limiting or enhancing Dengue-2 infectivity.

The data presented here have identified differences in gene expression between feral populations of Ae. aegypti that are naturally susceptible or refractory to Dengue-2 virus. There was an over-expression of numerous molecules and the involvement of diverse biological processes showed the complexity of viral infection and immune responses against the virus. The functional characterizations of the apoptosis-related genes have begun to be evaluated in order to elucidate their role in the susceptible or refractory phenotypes. In addition, more investigations need to be done in order to evaluate whether known immune pathways (Toll and Imd) are activated after dengue virus infection.


The authors are grateful to Paola Caicedo and Andrea Arévalo, young investigators at CIDEEVI who participated in the initial vector competence studies. We also thank Dawn Cooper and Mathias Schuetz, who contributed, respectively, to the primers design and radioactive screening of libraries in Canada. This project was supported by WHO-TDR (A30848) and COLCIENCIAS (Contract: 22290416474). Other partial support was from the CHIR grant (69558), the Canada Research Chairs program, and a Michael Smith Scholar award to CL.




expressed sequence tag;



suppressive subtractive hybridization



E Abdelwahid , T Yokokura , et al. 2007. Mitochondrial disruption in Drosophila apoptosis. Developmental Cell 12(5): 793– 806. Google Scholar


EG Abraham , M Donnelly-Doman , et al. 2005. Driving midgut-specific expression and secretion of a foreign protein in transgenic mosquitoes with AgAper1 regulatory elements. Insect Molecular Biology 14(3): 271–9. Google Scholar


ZN Adelman , CD Blair , et al. 2001. Sindbis virus-induced silencing of dengue viruses in mosquitoes. Insect Molecular Biology 10(3): 265–73. Google Scholar


S Aksoy , I Maudlin , et al. 2001. Prospects for control of African trypanosomiasis by tsetse vector manipulation. Trends in Parasitology 17(1): 29–35. Google Scholar


EM Al-Olayan , GT Williams , et al. 2002. Apoptosis in the malaria protozoan, Plasmodium berghei: a possible mechanism for limiting intensity of infection in the mosquito. International Journal for Parasitology 32(9): 1133–43. Google Scholar


L Alphey , CB Beard 2002. Malaria control with genetically manipulated insect vectors. Science 298(5591): 119–121. Google Scholar


BL Apostol , WC Black 1996. Population genetics with RAPD-PCR markers: the breeding structure of Aedes aegypti in Puerto Rico. Heredity 76(Pt 4): 325–34. Google Scholar


E Arama , M Bader , et al. 2007. A Ubiquitin Ligase Complex Regulates Caspase Activation During Sperm Differentiation in Drosophila. PLoS Biology 5(10): e251. Google Scholar


E Arama , M Bader 2006. The two Drosophila cytochrome C proteins can function in both respiration and caspase activation. Embo Journal 25(1): 232–43. Google Scholar


FM Ausubel , 2005. Are innate immune signaling pathways in plants and animals conserved? Nature Immunology 6(10): 973–9. Google Scholar


LC Bartholomay , JF Fuchs 2004. Reassessing the role of defensin in the innate immune response of the mosquito, Aedes aegypti. Insect Molecular Biology 13(2): 125– 32. Google Scholar


BJ Beaty , 2000. Genetic manipulation of vectors: A potential novel approach for control of vector-borne diseases. Proceedings of the National Academy of Sciences of the USA 97(19): 10295–10297. Google Scholar


KE Bennett , KE Olson 2002. Variation in vector competence for dengue 2 virus among 24 collections of Aedes aegypti from Mexico and the United States. American Journal of Tropical Medicine & Hygiene 67(1): 85–92. Google Scholar


KE Bennett BJ Beaty , et al. 2005. Selection of D2S3, an Aedes aegypti (Diptera: Culicidae) Strain with High Oral Susceptibility to Dengue 2 Virus and D2MEB, a Strain with a Midgut Barrier to Dengue 2 Escape. Journal of Medical Entomology 42(2): 110–119. Google Scholar


G Bidla , M Lindgren 2005. Hemolymph coagulation and phenoloxidase in Drosophila larvae. Developmental & Comparative Immunology 29(8): 669–79. Google Scholar


WC Black , KE Bennett , et al. 2002. Flavivirus susceptibility in Aedes aegypti. Archive of Medical Research 33 (4): 379–88. Google Scholar


BJ Blitvich , CD Blair 2002. Developmental- and tissue-specific expression of an inhibitor of apoptosis protein 1 homologue from Aedes triseriatus mosquitoes. Insect Molecular Biology 11(5): 431–42. Google Scholar


CF Bosio , RE Fulton 2000. Quantitative trait loci that control vector competence for dengue-2 virus in the mosquito Aedes aegypti. Genetics 156(2): 687–98. Google Scholar


DE Brackney , BD Foy , KE Olson . 2008. The Effects of Midgut Serine Proteases on Dengue Virus Type 2 Infectivity of Aedes aegypti. American Journal of Tropical Medicine & Hygiene 79(2): 267–274. Google Scholar


B Bryant , CD Blair 2008. Annotation and expression profiling of apoptosis-related genes in the yellow fever mosquito, Aedes aegypti. Insect Biochemistry and Molecular Biology 38(3): 331–45. Google Scholar


LL Cheng , LC Bartholomay 2001. Characterization of an endogenous gene expressed in Aedes aegypti using an orally infectious recombinant Sindbis virus. Journal of Insect Science 1: 10. Available online:  http://www.insectscience.Org/1.10  Google Scholar


S Cherry , N Silverman . 2006. Host-pathogen interactions in Drosophila: new tricks from an old friend. Nature Immunology 7(9): 911–7. Google Scholar


TE Clarke , RJ Clem . 2003. Insect defenses against virus infection: the role of apoptosis. International Reviews of Immunology 22(5–6): 401–24. Google Scholar


DM Cooper , F Pio 2007a. Characterization of Aedes Dredd: a novel initiator caspase from the yellow fever mosquito, Aedes aegypti. Insect biochemistry and molecular biology 37(6): 559–69. Google Scholar


DM Cooper , EP Thi 2007b. Aedes Drone: a novel ecdysone-inducible caspase in the yellow fever mosquito, Aedes aegypti. Insect molecular biology 16(5): 563–72. Google Scholar


L Dorstyn , K Mills 2004. The two cytochrome c species, DC3 and DC4, are not required for caspase activation and apoptosis in Drosophila cells. The Journal of Cell Biology 167(3): 405–10. Google Scholar


C Dostert , E Jouanguy 2005. The JakSTAT signaling pathway is required but not sufficient for the antiviral response of Drosophila. Nature Immunology 6(9): 946–53. Google Scholar


EM Dotson , B Plikaytis 2003. Transformation of Rhodococcus rhodnii, a symbiont of the Chagas disease vector Rhodnius prolixus, with integrative elements of the Ll mycobacteriophage. Infection Genetics and Evolution 3(2): 103–9. Google Scholar


H Everett , G McFadden . 1999. Apoptosis: an innate immune response to virus infection. Trends in Microbiology 7(4): 160–5. Google Scholar


JH Fritz , SE Girardin , et al. 2006. Innate immune defense through RNA interference. Science Signaling: STKE 2006(339): pe27. Google Scholar


F Garcia-Franco , ML Munoz 2002. Large genetic distances among Aedes aegypti populations along the South Pacific coast of Mexico. American Journal of Tropical Medicine & Hygiene 66(5): 594–8. Google Scholar


N Gorrochotegui-Escalante , C Gomez-Machorro 2002. Breeding structure of Aedes aegypti populations in Mexico varies by region. American Journal of Tropical Medicine & Hygiene 66(2): 213–22. Google Scholar


G Goyal , B Fell 2007. Role of mitochondrial remodeling in programmed cell death in Drosophila melanogaster. Developmental Cell 12(5): 807–16. Google Scholar


LM Grimm , BA Osborne 1999. Apoptosis and the proteasome. Results and Problems in Cell Differentiation 23: 209–28. Google Scholar


DJ Gubler , 2004. The changing epidemiology of yellow fever and dengue, 1900 to 2003: full circle? Comparative Immunology, Microbiology & Infectious Diseases 27(5): 319–30. Google Scholar


Guha-SapirD , B Schimmer . 2005. Dengue fever: new paradigms for a changing epidemiology. Emerging Themes in Epidemiology 2(1): 1. Google Scholar


MG Guzman , G Kouri . 2002. Dengue: an update. The Lancet Infectious Diseases 2(1): 33–42. Google Scholar


MO Hengartner . 2000. The biochemistry of apoptosis. Nature 407(6805): 770–6. Google Scholar


S Higgs , BJ Beaty 1996. The biology of Disease Vectors. 1st edition. University Press of Colorado. Google Scholar


JA Hoffmann , M Reichhart . 2002. Drosophila innate immunity: an evolutionary perspective. Nature Immunology 3(2): 121–6. Google Scholar


S Hu , X Yang . 2003. Cellular inhibitor of apoptosis 1 and 2 are ubiquitin ligases for the apoptosis inducer Smac/DIABLO. Journal of Biological Chemistry 278(12): 10055–60. Google Scholar


JR Huh , I Foe 2007. The Drosophila inhibitor of apoptosis (IAP) DIAP2 is dispensable for cell survival, required for the innate immune response to gram-negative bacterial infection, and can be negatively regulated by the reaper/hid/grim family of IAP-binding apoptosis inducers. Journal of Biological Chemistry 282(3): 2056–68. Google Scholar


JL Imler , P Bulet . 2005. Antimicrobial peptides in Drosophila: structures, activities and gene regulation. Chemical Immunology and Allergy 86: 1–21. Google Scholar


PM Irusta , YB Chen 2003. Viral modulators of cell death provide new links to old pathways. Current Opinion in Cell Biology 15(6): 700–5. Google Scholar


C Karlsson , AM Korayem 2004. Proteomic analysis of the Drosophila larval hemolymph clot. Journal of Biological Chemistry 279(50): 52033–41. Google Scholar


KM Keene , BD Foy 2004. RNA interference acts as a natural antiviral response to O'nyong-nyong virus (Alphavirus; Togaviridae) infection of Anopheles gambiae. Proceedings of the National Academy of Sciences of the United States of America 101(49): 17240–5. Google Scholar


C Kocks , JH Cho 2005. Eater, a transmembrane protein mediating phagocytosis of bacterial pathogens in Drosophila. Cell 123(2): 335–46. Google Scholar


RS Lanciotti , CH Calisher 1992. Rapid detection and typing of dengue viruses from clinical samples by using reverse transcriptase-polymerase chain reaction. Journal of Clinical Microbiology 30(3): 545– 51. Google Scholar


B Lemaitre , J Hoffmann . 2007. The host defense of Drosophila melanogaster. Annual Reviews Immunology 25: 697–743. Google Scholar


JH Leu , YC Kuo 2008. Molecular cloning and characterization of an inhibitor of apoptosis protein (IAP) from the tiger shrimp, Penaeus monodon. Developmental & Comparative Immunology 32(2): 121–33. Google Scholar


Q Li , H Li 2007. The Aedes albopictus inhibitor of apoptosis 1 gene protects vertebrate cells from bluetongue virusinduced apoptosis. Insect Molecular Biology 16(1): 93–105. Google Scholar


L Liu , J Peng 2007. Influence of cytochrome c on apoptosis induced by Anagrapha (Syngrapha) falcifera multiple nuclear polyhedrosis virus (AfMNPV) in insect Spodoptera litura cells. Cell Biology international 31(9): 996–1001. Google Scholar


CA Lowenberger , 2001. Form, function and phylogenetic relationships of mosquito immune peptides. Advances in experimental medicine and biology 484: 113–29. Google Scholar


SJ Lu , JE Pennington 2006. Reevaluation of the role of early trypsin activity in the transcriptional activation of the late trypsin gene in the mosquito Aedes aegypti. Insect biochemistry and molecular biology 36(4): 336–43. Google Scholar


GV Ludwig , BA Israel 1991. Role of La Crosse virus glycoproteins in attachment of virus to host cells. Virology 181(2): 564–71. Google Scholar


T Maniatis , ( 1999). A ubiquitin ligase complex essential for the NF-kappaB, Wnt/Wingless, and Hedgehog signaling pathways. Genes & Development 13(5): 505– 10. Google Scholar


S Meister , SM Kanzok 2005. Immune signaling pathways regulating bacterial and malaria parasite infection of the mosquito Anopheles gambiae. Proceedings of the National Academy of Sciences of the United States of America 102(32): 11420–5. Google Scholar


K Michel , FC Kafatos . 2005. Mosquito immunity against Plasmodium. Insect Biochemestry and Molecular Biology 35: 677– 89. Google Scholar


BR Miller , CJ Mitchell . 1991. Genetic selection of a flavivirus-refractory strain of the yellow fever mosquito Aedes aegypti. American Journal of Tropical Medicine & Hygiene 45(4): 399–407. Google Scholar


A Molina-Cruz , L Gupta , et al. 2005. Effect of mosquito midgut trypsin activity on Dengue-2 virus infection and dissemination in Aedes aegypti. American Journal of Tropical Medicine & Hygiene 72 (5): 631–37. Google Scholar


V Nene , JR Wortman 2007. Genome sequence of Aedes aegypti, a major arbovirus vector. Science 316(5832): 1718–23. Google Scholar


F Noriega , M Wells . 1999. A molecular view of trypsin synthesis in the midgut of Aedes aegypti. Journal of Insect Physiology 45: 613– 20. Google Scholar


V O'Brien , 1998. Viruses and apoptosis. Journal of General Virology 79 (8): 1833–45. Google Scholar


CB Ocampo , DM Wesson 2004. Population dynamics of Aedes aegypti from a dengue hyperendemic urban setting in Colombia. American Journal of Tropical Medicine & Hygiene 71(4): 506–13. Google Scholar


RZ Orlowski , 1999. The role of the ubiquitinproteasome pathway in apoptosis. Cell Death & Differentiation 6(4): 303–13. Google Scholar


S Paskewitz , B. M. Christensen 1996. Immune responses of vectors. In BJ Beaty , WC Marquardt editors. The Biology of Disease Vectors, 632. University Press of Colorado Google Scholar


SM Paskewitz , O Andreev 2006. Gene silencing of serine proteases affects melanization of Sephadex beads in Anopheles gambiae. Insect Biochemistry and Molecular Biology 36(9): 701–11. Google Scholar


A Paul , LC Harrington 2006. Evaluation of novel insecticides for control of dengue vector Aedes aegypti (Diptera: Culicidae). Journal of Medical Entomology 43(1): 55–60. Google Scholar


C Paupy , M Vazeille-Falcoz 2000. Aedes aegypti in Tahiti and Moorea (French Polynesia): isoenzyme differentiation in the mosquito population according to human population density. American Journal of Tropical Medicine & Hygiene 62(2): 217–24. Google Scholar


A Pearson , A Lux 1995. Expression cloning of dSR-CI, a class C macrophagespecific scavenger receptor from Drosophila melanogaster. Proceedings of the National Academy of Sciences of the USA 92(9): 4056– 60. Google Scholar


MA Riehle , M Jacobs-Lorena . 2005. Using bacteria to express and display anti-parasite molecules in mosquitoes: current and future strategies. Insect Biochemistry and Molecular Biology 35(7): 699–707. Google Scholar


HR Sanders , BD Foy 2005. Sindbis virus induces transport processes and alters expression of innate immunity pathway genes in the midgut of the disease vector, Aedes aegypti. Insect Biochemistry and Molecular Biology 35(11): 1293–307. Google Scholar


P Schmid-Hempel , 2001. On the evolutionary ecology of host-parasite interactions: addressing the question with regard to bumblebees and their parasites. Naturwissenschaften 88(4): 147–58. Google Scholar


P Schmid-Hempel , 2005. Evolutionary ecology of insect immune defenses. Annual Review of Entomology 50: 529–51. Google Scholar


BA Schreader , Y Wang 2003. Drosophila morgue and the intersection between protein ubiquitination and programmed cell death. Apoptosis 8(2): 129– 39. Google Scholar


JG Scott , S Kasai . 2004. Evolutionary plasticity of monooxygenase-mediated resistance. Pesticide Biochemistry and Physiology 78: 171–178. Google Scholar


RB Seth , L Sun 2006. Antiviral innate immunity pathways. Cell Research 16(2): 141–7. Google Scholar


SW Shin , V Kokoza 2003. Relishmediated immune deficiency in the transgenic mosquito Aedes aegypti. Proceedings of the National Academy of Sciences of the USA 100(5): 2616–21. Google Scholar


E Spencer , J Jiang 1999. Signal-induced ubiquitination of Iκ;Bα by the F-box protein Slimb/β-TrCP. Genes & Development 13(3): 284–94. Google Scholar


J Spiegel , S Bennett 2005. Barriers and Bridges to Prevention and Control of Dengue: The Need for a Social-Ecological Approach. EcoHealth 2: 273–290. Google Scholar


WJ Tabachnick , 1982. Geographic and temporal patterns of genetic variation of Aedes aegypti in New Orleans. American Journal of Tropical Medicine & Hygiene 31(4): 849–53. Google Scholar


RJ Ursic-Bedoya , CA Lowenberger . 2007. Rhodnius prolixus: identification of immunerelated genes up-regulated in response to pathogens and parasites using suppressive subtractive hybridization. Developmental & Comparative Immunology 31(2): 109–20. Google Scholar


M Vazeille-Falcoz , L Mousson 2001. Population genetic structure and competence as a vector for dengue type 2 virus of Aedes aegypti and Aedes albopictus from Madagascar. American Journal of Tropical Medicine & Hygiene 65(5): 491–7. Google Scholar


D Vlachou , T Schlegelmilch 2005. Functional genomic analysis of midgut epithelial responses in Anopheles during Plasmodium invasion. Current Biology 15(13): 1185–95. Google Scholar


GP Wallis , TH Aitken 1985. Selection for susceptibility and refractoriness of Aedes aegypti to oral infection with yellow fever virus. American Journal of Tropical Medicine & Hygiene 34(6): 1225–31. Google Scholar


X Wang , 2001. The expanding role of mitochondria in apoptosis. Genes & Development 15(22): 2922–33. Google Scholar


XH Wang , R Aliyari 2006. RNA interference directs innate immunity against viruses in adult Drosophila. Science 312(5772): 452–4. Google Scholar


RM Waterhouse , EV Kriventseva 2007. Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes. Science 316(5832): 1738–43. Google Scholar


FL Watson , R Puttmann-Holgado 2005. Extensive diversity of Ig-superfamily proteins in the immune system of insects. Science 309(5742): 1874–8. Google Scholar


C Wojcik , 2002. Regulation of apoptosis by the ubiquitin and proteasome pathway. Journal of Cellular and Molecular Medicine 6(1): 25–48. Google Scholar


JL Woodring , S Higgs 1996. Natural cycles of vector-borne pathogens. In: BJ B. Beaty, WC Marquardt . The biology of disease vectors. 51–72.University Press of Colorado. Google Scholar


G Xu , W Wilson , et al. 1997. VP7: an attachment protein of bluetongue virus for cellular receptors in Culicoides variipennis. The Journal of General Virology 78 (7): 1617–23. Google Scholar


RA Zambon , M Nandakumar 2005. The Toll pathway is important for an antiviral response in Drosophila. Proceedings of the National Academy of Sciences of the USA 102(20): 7257–62. Google Scholar


B Zhong , P Tien 2006. Innate immune responses: crosstalk of signaling and regulation of gene transcription. Virology 352(1): 14–21. Google Scholar
This is an open access paper. We use the Creative Commons Attribution 3.0 license that permits unrestricted use, provided that the paper is properly attributed.
Olga L. Barón, Raul J. Ursic-Bedoya, Carl A. Lowenberger, and Clara B. Ocampo "Differential Gene Expression from Midguts of Refractory and Susceptible Lines of the Mosquito, Aedes aegypti, Infected with Dengue-2 Virus," Journal of Insect Science 10(41), 1-23, (1 May 2010).
Received: 28 April 2008; Accepted: 1 June 2008; Published: 1 May 2010

Back to Top