- Research article
- Open Access
The mitogen-activated protein kinome from Anopheles gambiae: identification, phylogeny and functional characterization of the ERK, JNK and p38 MAP kinases
BMC Genomics volume 12, Article number: 574 (2011)
Anopheles gambiae is the primary mosquito vector of human malaria parasites in sub-Saharan Africa. To date, three innate immune signaling pathways, including the nuclear factor (NF)-kappaB-dependent Toll and immune deficient (IMD) pathways and the Janus kinase/signal transducers and activators of transcription (Jak-STAT) pathway, have been extensively characterized in An. gambiae. However, in addition to NF-kappaB-dependent signaling, three mitogen-activated protein kinase (MAPK) pathways regulated by JNK, ERK and p38 MAPK are critical mediators of innate immunity in other invertebrates and in mammals. Our understanding of the roles of the MAPK signaling cascades in anopheline innate immunity is limited, so identification of the encoded complement of these proteins, their upstream activators, and phosphorylation profiles in response to relevant immune signals was warranted.
In this study, we present the orthologs and phylogeny of 17 An. gambiae MAPKs, two of which were previously unknown and two others that were incompletely annotated. We also provide detailed temporal activation profiles for ERK, JNK, and p38 MAPK in An. gambiae cells in vitro to immune signals that are relevant to malaria parasite infection (human insulin, human transforming growth factor-beta1, hydrogen peroxide) and to bacterial lipopolysaccharide. These activation profiles and possible upstream regulatory pathways are interpreted in light of known MAPK signaling cascades.
The establishment of a MAPK "road map" based on the most advanced mosquito genome annotation can accelerate our understanding of host-pathogen interactions and broader physiology of An. gambiae and other mosquito species. Further, future efforts to develop predictive models of anopheline cell signaling responses, based on iterative construction and refinement of data-based and literature-based knowledge of the MAP kinase cascades and other networked pathways will facilitate identification of the "master signaling regulators" in biomedically important mosquito species.
Mitogen-activated protein kinases (MAPKs) are serine-threonine protein kinases that regulate a variety of cellular processes, including growth, metabolism, apoptosis, and innate immune responses [1–3]. MAPKs function in multi-tiered signaling cascades, in which an activated MAP4K phosphorylates and activates a MAP3K which, in turn, activates a downstream MAP2K, which activates a MAPK that can regulate effector proteins or transcription factors to positively or negatively regulate suites of genes [4, 5]. MAPK signaling modules provide multiple levels of regulation that confer signal amplification and specificity toward a desired outcome . A wide assortment of stimuli activate MAPKs, including inflammatory cytokines , osmotic stress , oxidative stress and redox signaling , and growth factors [9, 10].
MAPKs have been extensively studied and a wealth of information is available from many model systems, including Caenorhabditis elegans, Drosophila melanogaster and a variety of mammals [11–13]. From an evolutionary standpoint, MAPKs have diverged very little over time and several published phylogenies of MAPKs have revealed a high degree of conservation from invertebrates to vertebrates [14, 15]. Further, these analyses have contributed to our understanding of the evolution and function of the MAPKs [14, 15]. For example, a MAPK phylogeny was constructed from the encoded sequences in the genome of the human pathogenic blood fluke, Schistosoma japonicum, together with known eukaryotic MAPKs from model organisms to elucidate putative functions of previously undescribed S. japonicum MAPKs . The construction of MAPK phylogenies can, therefore, facilitate predictions of the roles of MAPKs in non-model organisms, including those of public health importance.
Malaria is a parasitic disease of great public health concern, with over 250 million new cases per year, resulting in nearly one million deaths annually . In sub-Saharan Africa, the mosquito Anopheles gambiae transmits the most deadly human malaria parasite Plasmodium falciparum. Despite highly efficient transmission, the invertebrate and vertebrate hosts of malaria parasites can mount sophisticated immune responses to infection. These responses are regulated in both hosts, in part, by MAPKs [9, 18–20].
Two prominent parasite-derived signals - glycosylphosphatidylinositols (GPIs) and hemozoin - activate MAPK signaling in both the mammalian and mosquito hosts. Mammalian JNK, ERK and p38 MAPKs transduce signals from P. falciparum glycophosphatidylinositols (PfGPIs) for inflammatory cytokine synthesis in immune cells in vitro and during parasite infection in vivo[18, 19]. Hemozoin signals principally through ERK to increase interferon-gamma-dependent production of anti-parasite nitric oxide (NO) in mammalian cells [21, 22]. In an analogous fashion, PfGPIs function as an early signal of parasite infection in An. gambiae and in Anopheles stephensi, a vector of malaria in Asia and close relative of An. gambiae. In An. stephensi, PfGPIs robustly activate MEK-ERK phosphorylation in the mosquito midgut epithelium , a site that is critical for parasite development in the insect host. As in mammalian cells, hemozoin can activate MEK-ERK signaling in the An. stephensi midgut . Further, transforming growth factor (TGF)-beta1-dependent MEK-ERK-dependent signaling can facilitate P. falciparum development at the midgut epithelium by inhibiting the expression of NO synthase  and synthesis of inflammatory levels of reactive nitrogen oxides that limit parasite development [26–28].
In contrast to our understanding of ERK signaling in An. stephensi and An. gambiae, our knowledge of the regulatory ligands and signaling pathways as well as the biological impacts of JNK and p38 MAPK signaling in these species is more limited. In particular, JNK activation appears to mediate as yet unidentified inhibitory responses to the murine parasite Plasmodium berghei in An. gambiae, although the extracellular signals and upstream regulatory proteins for JNK activation are unknown. In the case of p38 MAPK, this signaling protein appears to regulate antimicrobial responses in Aedes mosquitoes [30–32] and p38 MAPK is activated during insulin-dependent immune signaling in An. stephensi. However, the pathway elements that regulate this signaling, the biological effects of signaling, and the identity of other activating ligands in Anopheles spp., including those derived from malaria parasites , are unknown.
In this study, we used bioinformatics, existing knowledge of mammalian signaling pathways, and signaling assays in An. gambiae cells to provide a more comprehensive understanding of the likely roles of the MAPK cascades in innate immunity of An. gambiae, a critical malaria vector mosquito for which an annotated genome sequence is available . The MAPK phylogeny presented here recapitulates the conservation of the MAPKs and, for the first time, identifies the orthologous associations and evolutionary origins of the complete suite of An. gambiae MAPKs. The functional data presented herein demonstrate differing roles for An. gambiae p38 MAPK, ERK and JNK in response to stimuli that are biologically relevant for malaria parasite infection and, together with existing knowledge from mammalian models, suggest that anopheline mosquitoes utilize a unique MAPK architecture for signaling.
The initial set of 43 human MAPKs used for reference in this study was identified in the GeneCards database  (Additional File 1). Fifteen An. gambiae MAPK orthologs were identified in GeneCards and confirmed using Basic Local Alignment Search Tool (BLAST)  against An. gambiae genome sequence data. In addition, An. gambiae sequence data were analyzed using highly conserved, orthologous MAPK activation loop sequences. The activation of MAP kinases occurs by a dual phosphorylation event by the upstream kinase on a solvent-exposed activation loop that typically contains a T-X-Y motif [36, 37]. We used the encoded activation loop sequences from D. melanogaster Slipper (KTLKITDFGLAREAGTYAWMPPEVISV) and D. melanogaster Wallenda (EVVKISDFGTSREGTVAWMAPEVIRNPCSEKVDIWSY), flanked by 5-10 additional conserved amino residues, in TBLASTN queries against An. gambiae sequence data. Conserved catalytic domains for all predicted proteins were confirmed using Prosite.
For phylogenetic analyses, 108 sequences, ranging from 115-354 amino acids, of MAPK catalytic domains from An. gambiae, C. elegans, Ciona intestinalis, D. melanogaster, Homo sapiens and Saccharomyces cerevisiae were used to generate an alignment using the BLOSUM protein weight matrix and ClustalX . Orthologs of C. elegans, D. melanogaster and C. intestinalis MAPKs were confirmed using BLASTP of GeneCards, Ensemble Metazoa, Wormbase and Flybase databases. Phylogenetic analyses were conducted using Phylogeny Inference Package (PHYLIP) 3.69 . PROTDIST was used to generate a distance matrix based on the Jones-Taylor-Thornton (JTT) model of amino acid substitution. A distance-based phylogenetic tree was then inferred using the neighbor-joining (NJ) algorithm implemented in NEIGHBOR. SEQBOOT was used to prepare 1000 pseudoreplicate datasets for bootstrap analysis; these replicate datasets were then analyzed using PROTDIST with the JTT model and NEIGHBOR to create a bootstrap NJ tree. The bootstrap majority-rule consensus tree was generated using CONSENSE, visualized in FigTree v 1.3.1 http://tree.bio.ed.ac.uk/software/figtree/ and edited using Inkscape 0.48 http://inkscape.org/. The unrooted consensus tree was drawn to scale, with branch lengths representing the mean number of substitutions per site.
Mosquito cell culture and experimental treatments
The immortalized, embryo-derived An. gambiae 4a3B cell line (kindly provided by Hans-Michael Müller)  was maintained in Schneider's medium (Invitrogen, Carlsbad, CA) with 10% heat-inactivated fetal bovine serum (Invitrogen) at 28°C. For each condition, 1 × 106 4a3B cells in 2 mL medium were plated in one well of a 12-well tissue culture plate and allowed to recover overnight. Cells were then treated with the following stimuli at concentrations previously validated for signaling: 250 μM hydrogen peroxide (VWR International, Radnor, PA) , 6000 pg/ml human TGF-beta1 (R&D Systems, Minneapolis, MN) , 1.7 μM human insulin (Sigma-Aldrich, St. Louis, MO)  or 100 μg/ml lipopolysaccharide (LPS; Sigma-Aldrich)  and collected at 5 min, 15 min, 30 min, 1 h, 3 h, 6 h and 24 h after treatment for subsequent western blot analysis. To quantify and analyze ERK, JNK, and p38 MAPK transcript expression, 4a3B cells prepared as described above were treated with 1.7 μM human insulin and collected at 1 h, 3 h, 6 h and 24 h for qPCR (see below). Controls for each time point consisted of an identical aliquot of cells treated with diluent. At least three biological replicates were performed for each treatment.
Protein extraction and western blotting
Protein extracts of 4a3B cells were prepared by collecting cells in lysis buffer as previously described [9, 24]. Briefly, cell medium was removed following treatment and cells were washed with ice cold phosphate buffered saline (PBS) and lysed in the plate in 120 μl cell lysis buffer (10 mM Tris-HCl pH 7.4, 1 mM EDTA, 100 mM NaCl, 1 mM NaF, 1 mM EGTA, 2 mM Na3VO4, 20 mM Na4P2O7, 0.1% SDS, 1% Triton X-100, 0.5% sodium deoxycholate, 1 mM phenylmethylsulfonyl fluoride, 10% glycerol, 60 mg/ml aprotinin, 10 mg/ml leupeptin, 1 mg/ml pepstatin, and 1 mg/ml calyculin A). The plate was agitated for 30 min at 4°C and samples were incubated on ice for 30 min. Cell lysates were centrifuged at 14,000 × g for 10 min at 4°C to remove cellular debris and 100 μl of supernatant from each sample was mixed with 20 μl of 6 × sample buffer (125 mM Tris-HCl pH 6.8, 10% glycerol, 10% SDS, 0.006% bromophenol blue, 130 mM dithiothreitol) and heated at 95°C for 4 min.
Protein samples were electrophoretically separated on 10% polyacrylamide gels via sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to nitrocellulose membranes (BioRad, Hercules, CA). Protein loading was visually assessed by Coomassie blue staining. Membranes were blocked in 5% nonfat dry milk (NFDM) in 1 × Tris-buffered saline containing 0.1% Tween-20 (TBS-T) for 1 h at room temperature. For phosphorylated ERK detection, membranes were incubated at 4°C overnight with 1:10,000 mouse anti-phospho-ERK monoclonal antibody (Sigma-Aldrich) in 5% NFDM in TBS-T. For detection of phospho-p38, phospho-JNK and the loading control protein glyceraldehyde 3-phosphate dehydrogenase (GAPDH), membranes were incubated at 4°C overnight with 1: 1,250 rabbit anti-phospho-p38 MAPK antibody (Cayman Chemical, Ann Arbor, MI), 1:1,250 rabbit anti-phospho-JNK1/2 antibody (Biosource, Carlsbad, CA), or with 1:10,000 rabbit anti-GAPDH antibody (Abcam, San Francisco, CA) in 5% NFDM in TBS-T. Membranes were washed 3 times for 5 min in 1 × TBS-T and incubated with a 1:20,000 dilution of horse radish peroxidase (HRP)-conjugated rabbit anti-mouse IgG (Sigma-Aldrich) or with a 1:20,000 dilution of HRP-conjugated goat anti-rabbit (Fab')2 fragment (Cell Signaling Technology, Danvers, MA) at 4°C overnight. Following incubation with the secondary antibody, membranes were washed 3 times for 5 min in 1 × TBS-T. To reveal antibody-bound proteins, membranes were incubated with SuperSignal West Pico chemiluminescent reagent (Pierce, Rockford, IL) for 2-3 min. Each membrane was exposed to blue autoradiography film (ISC Bioexpress, Kaysville, UT). Phospho-MAPK levels and total GAPDH levels were quantified on scanned film using a GS-800 calibrated densitometer (BioRad, Hercules, California).
MAPK quantitative PCR (qPCR)
Primers for qPCR of An. gambiae ERK, p38, JNKa, and JNKb were designed using Primer3Plus http://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi. The primers included: ERK, 5'ATCCCGAGCACGATCACA3' (forward) and 5'CGATTTTGTGTAGCCCTTGGA3' (reverse); p38 MAPK, 5'CGGACCACATTCACCAGCTA3' (forward) and 5'CGCTAAAGTTGCGCTTCTCC3' (reverse); JNKa, 5'GCACGCAGCGATACATTAGC3' (forward) and 5'GTCCACCGAAATCCTTTCCA3' (reverse); and JNKb, 5'GACCGACTCGAACGAGCAC3' (forward) and 5'CGAACCACACGTTGATGTAGC3' (reverse). The housekeeping gene encoding ribosomal protein S7 was analyzed as a control with the primers 5'GAAGGCCTTCCAGAAGGTACAGA3' (forward) and 5'CATCGGTTTGGGCAGAATG3' (reverse).
RNA was isolated from 4a3B cells using Trizol regent (Invitrogen) and contaminating DNA was removed from the RNA samples using Turbo DNA-free™ (Ambion, Austin, TX). All RNA samples were diluted to 200 ng/μl and 6 μl (1.2 μg) RNA from each sample was reverse-transcribed using SuperScript® III (Invitrogen). Sample cDNAs were used to perform qPCR using Maxima SYBR Green/ROX qPCR Master Mix (Fermentas, Glen Burnie, MD). The qPCR cycling conditions were as follows: 2 min at 50°C, 10 min at 95°C, and 40 cycles of 15 s at 95°C and 1 min at 60°C. Three biological replicates with independent batches of 4a3B cells were analyzed for the diluent controls and treatments. No template controls were included for each reaction plate and all reactions were performed in triplicate to confirm amplification consistency (e.g., all reactions within 0.5 Ct of each other). The triplicate data were used to determine an average Ct for each reaction and the data analyzed using the 2-ΔΔCt method as described . Briefly, the average Ct value for each target was first normalized to the average Ct for ribosomal protein S7 gene to obtain ΔCt, and then the ΔCt value for each insulin-treated sample was normalized to timepoint-matched control to obtain ΔΔCt. Fold-change of expression relative to the diluent control was determined using 2-ΔΔCt.
For western blotting data, the Kolmogorov-Smirnov test revealed that non-normal distributions were present in our data sets (Graphpad Prism 5.02, San Diego, California). Overall significance values of all data sets, therefore, were determined by the Kruskal-Wallis test and Dunn's post test for pairwise comparisons (α = 0.05). For the q-PCR data, one way ANOVA and Student-Neuman-Keuls (Graphpad Prism 5.02) were used for analyses of significance (α = 0.05).
We used a combination of bioinformatic techniques to identify a total of 17 MAPKs from the An. gambiae genome (Additional File 1). Among these, AGAP002710 and AGAP006461 were newly identified MAPK orthologs. Based on homology, AGAP002710 is most appropriately defined as An. gambiae MAP3K10/11 or mixed lineage kinase (MLK) 2/3 (Additional File 1). Similarly, AGAP006461 is most appropriately defined as An. gambiae MAP3K 12/13 (Additional File 1).
AGAP012148 and AGAP009460 are annotated as p38 MAPK and JNK orthologs, respectively, in Vectorbase . However, inspection of AGAP012148 by alignment with p38 MAPK orthologs and Reverse Psi (RPS)-BLAST for conserved domains revealed that the annotated sequence was missing key N-terminal residues in the kinase catalytic domain. Specifically, RPS-BLAST of the annotated sequence predicted only 4 of 26 active site and 4 of 11 activation loop residues (not shown). Because the annotation appeared to be incomplete, we used the first 10 encoded amino acids from the N-terminus of the Ae. aegypti p38 MAPK ortholog (GenBank XP_001653239.1) to query the National Center for Biotechnology Information (NCBI) An. gambiae trace archives database (translated in all 6 reading frames) and identified the probable N-terminus of AGAP012148 (Additional File 2). Similarly, the annotation of AGAP009460 in Vectorbase  does not encode a start methionine. Approximately 150,000 base pairs of sequence upstream of the 5'-most nucleotides of AGAP009460 were translated in the three forward reading frames and analyzed by BLASTX using the first 20 amino acids of the corresponding orthologs from Ae. aegypti (AAEL008622) and from Culex quinquefasciatus (CPIJ001156, CPIJ001157). A putative start methionine with Kozak consensus and additional N-terminal amino acids extended the encoded sequence of AGAP009460 (Additional File 3).
The MAP4Ks, the MAP3Ks, the MAP2Ks are present in 1:1 orthology with gene products that have been identified in D. melanogaster (Additional File 1). Among the MAPKs, only D. melanogaster Rolled and An. gambiae ERK are in 1:1 orthology. In contrast, the D. melanogaster genome encodes two p38 MAPK orthologs (p38b, MPK2), a single JNK ortholog (Basket), an ERK5 ortholog (p38c), and an ERK7 ortholog, although the latter two are not associated with prototypical MAP kinase cascades . The An. gambiae genome, on the other hand, encodes a single p38 MAPK, two JNK orthologs (JNKa, JNKb; 57% identical), and no ERK5 and ERK7 orthologs. Based on sequence homology, An. gambiae JNKb is most closely related to D. melanogaster Basket (Additional File 1). The lack of strict 1:1 orthology among the MAPKs is perhaps not surprising because diversification among these proteins to accommodate unique aspects of host biology would have a greater likelihood of being tolerated relative to upstream changes in the signaling cascades that could result in significantly amplified non-adaptive phenotypes.
To determine whether the complement of MAPKs in the An. gambiae genome was similar to that in other mosquito genomes, we examined the available ERK, JNK, and p38 MAPK orthologs for Ae. aegypti and Cx. quinquefasciatus (Additional File 4). The genome sequences of Ae. aegypti and Cx. quinquefasciatus are currently assembled into supercontigs in Vectorbase  (AaegL1.2, September 2009; CpipJ1.2, June 2008). In this context, if identical nucleotide sequences (due perhaps to sequence assembly) are excluded, ERK, JNKa, and JNKb orthologs are in 1:1:1 orthology among the three mosquito species. In contrast, Cx. quinquefasciatus appears to possess two relatively short MAPKs that are tentatively identified as p38 MAPKs. The encoded sequences share significant identity to Ae. aegypti p38 MAPK and the putative full-length An. gambiae p38 MAPK, but are non-overlapping sequences, so it is impossible to discern whether these are indeed truly distinct gene products encoded in the Cx. quinquefasciatus genome. Because these relationships cannot be confirmed with the existing genome assemblies, the Ae. aegypti and Cx. quinquefasciatus MAPK sequences were excluded from the phylogenetic analysis.
Phylogenetic analysis of the MAPK catalytic domains (Figure 1) showed strong support for the MAP4K and MAP2K groups (100% bootstrap support). The MAP3K Raf superfamily formed a group (81% bootstrap support), but the MAP3K MEKK family did not form an exclusive group. Similarly, the two MAP3K superfamilies did not group together as an exclusive group and the sequences representing MAPK were not supported as a group in this analysis. The lack of a close relationship between the Raf and MEKK MAP3K groups was not unexpected based on previous phylogenetic analyses [44, 45]. The MEKK family proteins (MEKK1-4, ASK1-3) and the Raf superfamily proteins (MLK1-3, DLK, TAK1) are both associated with p38 and JNK signaling pathways; however, Raf proteins also figure prominently in MEK/ERK signaling (RAF, MOS) . Although not reflected in this analysis, the previously noted division of the ERK subfamily (ERK1, 2, and 5 versus ERK3, 4, 7/8) may be due to the classification of ERK3, 4 and 7/8 as MAPK-like cascade independent kinases . Frequently, sequences from the two deuterostomes (humans and C. intestinalis) and the insects (An. gambiae and D. melanogaster) grouped together. Exceptions to these groupings occurred when there were no clear homologs (e.g., HsMEK5 and HsPRKE1), which likely reflects the sequence diversification that has resulted in the many MAPKs that are found in humans relative to other species.
The identification of the complement of MAP4Ks, MAP3Ks, MAP2Ks and MAPKs provides a framework for the prediction of patterns of activation of cell signaling cascades in An. gambiae cells. To attempt to determine whether these pathways are phenotypically represented in this mosquito species, we selected three stimuli (hydrogen peroxide, human insulin, and human TGF-beta1) that we had previously implicated in ERK- and p38 MAPK-dependent An. gambiae immune signaling , as well as lipopolysaccharide (LPS), an activator of immune signaling pathways in mosquito cells [41, 46], to stimulate An. gambiae 4a3B cells over an extended timecourse from 5 min to 24 h. A representative western blot is shown in Additional File 5.
Stimulation with 250 μM hydrogen peroxide induced significant activation of p38 MAPK at 30 min relative to the control group (Figure 2A, middle panel). Signaling declined from 30 min to 3 h post-stimulation, with a return to control levels by 3 h (Figure 2A, middle panel). Despite trends toward induction within 1 h of stimulation, ERK and JNK were not significantly activated by 250 μM hydrogen peroxide in 4a3B cells (Figure 2A, left and right panels). In previous studies, stimulation of 4a3B cells with 500 μM hydrogen peroxide induced activation of p38 MAPK as well as ERK by 15 min post-treatment , a faster and more extensive MAPK induction than observed here for 250 μM hydrogen peroxide (Figure 2A, middle panel).
In accord with previous observations at 15 min post-treatment , stimulation with 1.7 μM human insulin induced significant activation of ERK and p38 MAPK at 15 min relative to the control and 5 min time point (Figure 2B, left and middle panels). In contrast to the more temporally limited activation by hydrogen peroxide, insulin-activated ERK and p38 MAPK were detectable to 3 h and 24 h post-stimulation, respectively, suggesting a greater persistence of these signaling responses. Phosphorylation of JNK was slightly, although not significantly, elevated relative to control levels within 1 h after stimulation (Figure 2B, right panel).
Stimulation with 6000 pg/ml human TGF-beta1 induced significant activation of ERK only at 15 min post-treatment relative to the 5 min time point (Figure 2C, left panel), with detectable levels of activated ERK through 1 h post-stimulation. This induction of ERK activation at 15 min, together with a lack of p38 MAPK activation, was noted previously . In contrast to slight trends toward JNK activation by hydrogen peroxide and insulin (Figures 2A and 2B, right panels), phosphorylated JNK levels never exceeded baseline control levels at all times after stimulation with human TGF-beta1 (Figure 2C, right panel).
Stimulation with 100 μg/ml LPS presented a unique pattern of MAPK activation in An. gambiae 4a3B cells. In contrast to ERK and p38 MAPK activation by hydrogen peroxide, insulin and TGF-beta1 (Figures 2A-C), LPS stimulation did not induce mean phosphorylated ERK levels above baseline (Figure 2D, left panel) and only non-significant trends toward elevated phosphorylated p38 MAPK were noted through 30 min post-stimulation (Figure 2D, middle panel). Levels of phosphorylated JNK, however, were significantly increased by LPS stimulation (Figure 2D, right panel) at 5 min and 30 min post-treatment relative to control, with slightly elevated but not significant levels of phosphorylated JNK through 3 h post-treatment that declined to baseline control levels by 6 h post-treatment.
Although the majority of signaling protein genes do not exhibit transcriptional responsiveness to environmental stimuli (reviewed in ), we questioned whether we might detect a signaling response at the MAPK transcript level in An. gambiae cells. To address this question, we used qPCR to analyze expression levels of p38 MAPK, JNKa, JNKb, and ERK expression in response to human insulin, the only stimulus for which we detected changes in phosphorylation for more than one MAPK (Figure 2). In replicated assays, none of the four genes exhibited a significant change in expression from 1 h to 24 h after treatment with 1.7 μM human insulin (Additional File 6). Together with our phosphorylation data, these data suggest that MAPK responsiveness in An. gambiae cells occurs principally at the protein level.
The protein targets of the mammalian MAPKs (e.g., STAT, AP1, p300; Figure 3) [13, 43, 48] predict a strong association of the MAPKs with immune regulation. This association is also ancient, with strong representation in invertebrates. In C. elegans, p38 MAPK- and pERK-dependent signaling pathways, as well as insulin and TGF-beta signaling pathways, are required for host defense, with the general presumptions that insulin signaling functions in pathogen surveillance, p38 MAPK-dependent signaling functions in general defense, and p38 MAPK-independent signaling pathways provide pathogen-specific responses (reviewed in ). This nematode, however, lacks orthologs of Rel and nuclear factor (NF)-kappaB, the primary immune signaling transcription factors of D. melanogaster. In the fruit fly, Toll- or Imd-dependent signaling regulates NF-kappaB and JNK pathways, while p38 MAPK-dependent signaling mediates host defense against pathogen invasion [50, 51]. The JNK and p38 MAPKs also mediate control of antimicrobial peptide expression in D. melanogaster[50, 51], but cooperativity with NF-κB appears to be restricted to JNK. In particular, Imd activation results in the activation of TAK1 upstream of JNK and I-Kappa Kinase (IKK), the upstream activator of Rel that ultimately feeds back to reduce JNK activation .
The complement of MAP kinases encoded in the An. gambiae genome would provide this organism with the significant redundancy of p38 MAPK and JNK signaling observed in model organisms, although some intriguing differences between the mosquito and other species are predicted as well. In particular, An. gambiae MAP4K1 may function similarly to its mammalian counterpart downstream of Toll and RAC1 (AGAP005445) to activate the mosquito MLK2/3 ortholog (MAP3K10/11; Additional File 1) and p38 MAPK and JNK signaling (Figure 3). The second An. gambiae MAP4K - MAP4K4 - would be predicted to activate MAP3K1 and MAP3K7. However, An. gambiae lacks an ortholog for MAP3K1, an upstream activator for MEK/ERK signaling, so MAP4K4 activation of MAP3K7 (AgTAK1) would be predicted to activate NF-kappaB signaling through IKK (AGAP009166), and p38 MAPK and JNK signaling through MAP2K3 and MAP2K7 activation, respectively.
Among the other An. gambiae MAP3K orthologs, a number are predicted to lack well-defined upstream activators or to function independently of MAP4K activation. In particular, the activators of An. gambiae ortholog of MAP3K12/13 (Additional File 1) are not well defined based on data from mammalian models (Figure 3); nevertheless, the mosquito protein would be predicted to contribute to the activation of JNK signaling (Figure 3). Among the MAP3K orthologs that can be activated independently of MAP4Ks, An. gambiae MAP3K4 would be predicted to lie downstream of TGF-beta activation of GADD45 (AGAP007651), while MAP3K5 could be activated by TGF-beta-mediated activation of DAXX (AGAP009432) or TNF/TNFR activation of TRAF2 (Figure 3). Eiger and wengen have been identified as D. melanogaster orthologs of TNF/TNFR and, in the context of innate immunity, eiger functions independently of the Toll and Imd signaling pathways to regulate pathogen-induced cell death , proliferation of extracellular pathogens  and prophenoloxidase-mediated defense responses . Although a clear ortholog of TRAF2 is not apparent in the An. gambiae genome, orthologs of eiger (AGAP006771) and wengen (AGAP000728) would suggest that signaling by this ligand and receptor could lead to MAP3K5/MAP2K3 activation of p38 MAPK.
An examination of possible pathways leading to MEK/ERK signaling based on the encoded MAP kinases reveals that, in addition to a lack of MAP3K1, An. gambiae also lacks an ortholog for MAP3K8, so the regulation of MEK/ERK signaling is likely dependent on B-RAF and MOS activation, which are the downstream targets for RAP1 (AGAP001874) and for RAS (AGAP002219, AGAP002812), a major transducer for multiple signaling inputs from growth factors and cytokines (e.g., insulin, IGF-1, and TGF-beta; Figure 3).
The data from our functional studies (Figure 2) not only confirm a subset of signaling predictions from known pathways (Figure 3), but also suggest that other pathways may be involved in An. gambiae MAPK activation. In particular, human insulin-induced activation of An. gambiae ERK was predicted (Figure 3) and this was consistent with previous observations . In contrast, activation of An. gambiae p38 MAPK by insulin, which may occur through hydrogen peroxide at a concentration that is insufficient for ERK activation (Figure 2 compare panels A and B), is not predicted to flow through MAP4K-MAP3K-MAP2K activation (Figure 3). Interestingly, several studies have reported cross talk between the insulin receptor kinase and the Janus-activated kinase (JAK) signaling pathways [56–58]. In particular, JAK2 - the major mediator of inflammatory cytokine signaling in mammalian cells - can be activated by hydrogen peroxide  and contributes to insulin-dependent activation of p38 MAPK, JNK, and ERK in a manner that is independent of insulin receptor substrate activation and activation of the phosphatidylinositol 3-kinase/Akt-dependent arm of insulin signaling . Jaramillo-Gutierrez et al.  proposed, based on the regulation of antioxidant gene expression by JNK, that this MAPK may in fact be induced by oxidative stress to control redox biology in An. gambiae. These observations - together with our data and known models of MAPK signaling (Figure 3) - predict the existence of a pathway involving insulin-induced hydrogen peroxide regulation of a JAK ortholog (AGAP008354) and downstream activation of p38 MAPK and JNK that can be tested experimentally in An. gambiae.
In contrast to the prediction of ERK-specific signaling by insulin, existing mammalian signaling data (Figure 3) predict that human TGF-beta1 should activate ERK as well as p38 MAPK in An. gambiae cells. We have extensively characterized TGF-beta1-dependent ERK signaling in the regulation of P. falciparum development in An. stephensi[9, 20] and we have also demonstrated that p38 MAPK is activated by TGF-beta1 in An. stephensi cells . Given that p38 MAPK can be activated in An. gambiae cells (Figure 2B, center panel), our understanding of mammalian MAPK signaling (Figure 3) would suggest that perhaps a deficiency in the TGFBR-GADD45-MAP3K4 cascade is responsible for the failure of TGF-beta1 to activate p38 MAPK in An. gambiae 4a3B cells. A comparative analysis of An. stephensi and An. gambiae 4a3B cells, therefore, could be used to identify critical differences in these signaling pathways.
LPS has been used as an activator of a variety of immune signaling pathways in mosquito cells [41, 46]. A putative An. gambiae RAC1-MAP4K1-MAP3K11-MAP2K7-JNK signaling cascade (Figure 3) could provide one mechanism for LPS activation of JNK in An. gambiae 4a3B cells (Figure 2D, right panel). However, we could also predict, given that this mosquito MAP3K is orthologous with MAP3K10 as well as MAP3K11 that LPS should also activate p38 MAPK signaling through MAP2K4. Because LPS does not significantly activate p38 MAPK in An. gambiae 4a3B cells (Figure 2D, middle panel), we could speculate that this An. gambiae MAP3K is functionally more analogous to MAP3K11. An examination of MAP2K4 and MAP2K3 activation downstream of An. gambiae MAP3K10/11 using inhibitors or knock down strategies could clarify the orthology and functionality of this mosquito MAP3K.
Although LPS was the only pathogen-associated molecular pattern (PAMP) used in these signaling studies, the restriction of MAPK activation by LPS to An. gambiae JNK is insightful. Akman-Anderson et al.  showed that P. falciparum hemozoin induced activation of TAK1 and ERK as well as the signaling kinases Akt/protein kinase B and atypical protein kinase C zeta/lambda in An. gambiae 4a3B cells. Activation of TAK1 by this important parasite signaling factor, together with LPS-specific activation of JNK in our studies, suggests that JNK may function specifically in pathogen detection in An. gambiae cells. Indeed, following infection with P. berghei or P. falciparum, the mosquito midgut epithelium undergoes profound cytoskeletal changes similar to those regulated by the Toll/MAP2K4/JNK signaling module in D. melanogaster[61, 62].
Additional support of the importance, and perhaps specificity, of JNK signaling for control of natural P. falciparum infection in An. gambiae is provided by two population genetics studies. We identified a single nucleotide polymorphism (SNP) in the An. gambiae MAP2K4 (MKK4) gene that was in linkage disequilibrium with a SNP in the gene encoding Toll5B that was significantly associated with P. falciparum infection . This synonymous SNP, MKK43, introduces a reduction of codon frequency greater than 3-fold, which could result in changes in downstream protein expression or function and, hence, JNK signaling . In addition to our work, Riehle et al.  mapped a locus that is strongly associated with An. gambiae resistance to P. falciparum infection in east Africa that encompasses the chromosomal locations of MAP4K4 and MAP3K12, both of which are predicted to regulate JNK signaling (Figure 3).
In sum, we have identified the complement of An. gambiae MAP kinases and used a combination of predicted signaling relationships and bioassays to suggest novel interactions and functionality that can be tested experimentally in this biomedically important mosquito species. Clearly, the MAPK signaling cascades are the major regulators of innate immunity in C. elegans, in D. melanogaster, and in mammals, but we are only beginning to appreciate the importance of these complex cascades in innate immunity of the mosquito genera that serve as vectors of globally devastating pathogens of humans. The establishment of this "road map" based on the most advanced mosquito genome annotation can, therefore, accelerate our understanding of host-pathogen interactions and broader physiology in An. gambiae and in the closely related An. stephensi and provide a reasonable architecture for similar efforts in Ae. aegypti and Cx. quinquefasciatus. Further, future efforts to develop predictive models of anopheline cell signaling responses, based on iterative construction and refinement of data-based and literature-based knowledge of the MAP kinase cascades and other networked pathways (reviewed in ), will ultimately be necessary to elucidate the "master signaling regulators" in these biomedically important insects.
Apoptosis Signal-Regulating Kinase
Basic Local Alignment Search Tool for translated query against a protein database (X)
Death-associated protein 6
Dual Leucine zipper Kinase
Extracellular signal-Regulated Kinase
Growth Arrest and DNA-Damage inducible protein
Glyceraldehyde 3-Phosphate Dehydrogenase
Horse radish peroxidase
I Kappa b Kinase
Janus Kinase/Signal Transducers and Activators of Transcription
Mitogen-Activated Protein Kinase Kinase Kinase Kinase
Mitogen-Activated Protein Kinase Kinase Kinase
Mitogen-Activated Protein Kinase Kinase
Mitogen-activated Protein Kinase
Equivalent to MAP2K
Equivalent to MAP3K
Mixed Lineage Kinase
Nonfat Dry Milk
Plasmodium falciparum glycosylphosphatidylinositols
Equivalent to MAP3K9, MLK1
Ras-related C3 botulinum toxin substrate 1
Ras-related Protein 1
TGF-beta-activated Kinase 1
Basic Local Alignment Search Tool for protein query against a translated database (T)
1 × Tris-Buffered Saline with 0.1% Tween 20
Transforming Growth Factor-beta1
Tumor Necrosis Factor Receptor
/6: TNFR-associated Factor 2/6.
Zhang W, Liu HT: MAPK signal pathways in the regulation of cell proliferation in mammalian cells. Cell Res. 2002, 12: 9-18. 10.1038/sj.cr.7290105.
Roux PP, Blenis J: ERK and p38 MAPK-activated protein kinases: a family of protein kinases with diverse biological functions. Microbiol Mol Biol Rev. 2004, 68: 320-344. 10.1128/MMBR.68.2.320-344.2004.
Huang G, Shi LZ, Chi H: Regulation of JNK and p38 MAPK in the immune system: signal integration, propagation and termination. Cytokine. 2009, 48: 161-169. 10.1016/j.cyto.2009.08.002.
Chang L, Karin M: Mammalian MAP kinase signalling cascades. Nature. 2001, 410: 37-40. 10.1038/35065000.
Krishna M, Narang H: The complexity of mitogen-activated protein kinases (MAPKs) made simple. Cell Mol Life Sci. 2008, 65: 3525-3544. 10.1007/s00018-008-8170-7.
Haddad JJ: The role of inflammatory cytokines and NF-kappaB/MAPK signaling pathways in the evolution of familial Mediterranean fever: current clinical perspectives and potential therapeutic approaches. Cell Immunol. 2009, 260: 6-13. 10.1016/j.cellimm.2009.08.003.
Lee SH, Park Y, Yoon SK, Yoon JB: Osmotic stress inhibits proteasome by p38 MAPK-dependent phosphorylation. J Biol Chem. 2010, 285: 41280-41289. 10.1074/jbc.M110.182188.
Runchel C, Matsuzawa A, Ichijo H: Mitogen-activated protein kinases in mammalian oxidative stress responses. Antioxid Redox Signal. 2011,
Surachetpong W, Singh N, Cheung KW, Luckhart S: MAPK ERK signaling regulates the TGF-beta1-dependent mosquito response to Plasmodium falciparum. PLoS Pathog. 2009, 5: e1000366-10.1371/journal.ppat.1000366.
Cuadrado A, Nebreda AR: Mechanisms and functions of p38 MAPK signalling. Biochem J. 2010, 429: 403-417. 10.1042/BJ20100323.
Sakaguchi A, Matsumoto K, Hisamoto N: Roles of MAP kinase cascades in Caenorhabditis elegans. J Biochem. 2004, 136: 7-11. 10.1093/jb/mvh097.
Ragab A, Buechling T, Gesellchen V, Spirohn K, Boettcher AL, Boutros M: Drosophila Ras/MAPK signalling regulates innate immune responses in immune and intestinal stem cells. EMBO J. 2011, 30: 1123-1136. 10.1038/emboj.2011.4.
Cargnello M, Roux PP: Activation and function of the MAPKs and their substrates, the MAPK-activated protein kinases. Microbiol Mol Biol Rev. 2011, 75: 50-83. 10.1128/MMBR.00031-10.
Kultz D: Phylogenetic and functional classification of mitogen- and stress-activated protein kinases. J Mol Evol. 1998, 46: 571-588. 10.1007/PL00006338.
Caffrey DR, O'Neill LA, Shields DC: The evolution of the MAP kinase pathways: coduplication of interacting proteins leads to new signaling cascades. J Mol Evol. 1999, 49: 567-582. 10.1007/PL00006578.
Wang L, Yang Z, Li Y, Yu F, Brindley PJ, McManus DP, Wei D, Han Z, Feng Z, Li Y, Hu W: Reconstruction and in silico analysis of the MAPK signaling pathways in the human blood fluke, Schistosoma japonicum. FEBS Lett. 2006, 580: 3677-3686. 10.1016/j.febslet.2006.05.055.
Hay SI, Guerra CA, Gething PW, Patil AP, Tatem AJ, Noor AM, Kabaria CW, Manh BH, Elyazar IR, Brooker S, Smith DL, Moyeed RA, Snow RW: A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 2009, 6: e1000048-
Zhu J, Krishnegowda G, Gowda DC: Induction of proinflammatory responses in macrophages by the glycosylphosphatidylinositols of Plasmodium falciparum: the requirement of extracellular signal-regulated kinase, p38, c-Jun N-terminal kinase and NF-kappaB pathways for the expression of proinflammatory cytokines and nitric oxide. J Biol Chem. 2005, 280 (9): 8617-8627.
Lu Z, Serghides L, Patel SN, Degousee N, Rubin BB, Krishnegowda G, Gowda DC, Karin M, Kain KC: Disruption of JNK2 decreases the cytokine response to Plasmodium falciparum glycosylphosphatidylinositol in vitro and confers protection in a cerebral malaria model. J Immunol. 2006, 177: 6344-6352.
Surachetpong W, Pakpour N, Cheung KW, Luckhart S: Reactive oxygen species-dependent cell signaling regulates the mosquito immune response to Plasmodium falciparum. Antioxid Redox Signal. 2011, 14: 943-955. 10.1089/ars.2010.3401.
Jaramillo M, Gowda DC, Radzioch D, Olivier M: Hemozoin increases IFN-gamma-inducible macrophage nitric oxide generation through extracellular signal-regulated kinase- and NF-kappa B-dependent pathways. J Immunol. 2003, 171: 4243-4253.
Shio MT, Kassa FA, Bellemare MJ, Olivier M: Innate inflammatory response to the malarial pigment hemozoin. Microbes Infect. 2010, 12: 889-899. 10.1016/j.micinf.2010.07.001.
Arrighi RB, Debierre-Grockiego F, Schwarz RT, Faye I: The immunogenic properties of protozoan glycosylphosphatidylinositols in the mosquito Anopheles gambiae. Dev Comp Immunol. 2009, 33: 216-223. 10.1016/j.dci.2008.08.009.
Lim J, Gowda DC, Krishnegowda G, Luckhart S: Induction of nitric oxide synthase in Anopheles stephensi by Plasmodium falciparum: mechanism of signaling and the role of parasite glycosylphosphatidylinositols. Infect Immun. 2005, 73: 2778-2789. 10.1128/IAI.73.5.2778-2789.2005.
Akman-Anderson L, Olivier M, Luckhart S: Induction of nitric oxide synthase and activation of signaling proteins in Anopheles mosquitoes by the malaria pigment, hemozoin. Infect Immun. 2007, 75: 4012-4019. 10.1128/IAI.00645-07.
Luckhart S, Vodovotz Y, Cui L, Rosenberg R: The mosquito Anopheles stephensi limits malaria parasite development with inducible synthesis of nitric oxide. Proc Natl Acad Sci USA. 1998, 95: 5700-5705. 10.1073/pnas.95.10.5700.
Peterson TM, Gow AJ, Luckhart S: Nitric oxide metabolites induced in Anopheles stephensi control malaria parasite infection. Free Radic Biol Med. 2007, 42: 132-142. 10.1016/j.freeradbiomed.2006.10.037.
Ramiro RS, Alpedrinha J, Carter L, Gardner A, Reece SE: Sex and death: the effects of innate immune factors on the sexual reproduction of malaria parasites. PLoS Pathog. 2011, 7: e1001309-10.1371/journal.ppat.1001309.
Jaramillo-Gutierrez G, Molina-Cruz A, Kumar S, Barillas-Mury C: The Anopheles gambiae oxidation resistance 1 (OXR1) gene regulates expression of enzymes that detoxify reactive oxygen species. PLoS One. 2010, 5: e11168-10.1371/journal.pone.0011168.
Chen-Chih Wu R, Shaio MF, Cho WL: A p38 MAP kinase regulates the expression of the Aedes aegypti defensin gene in mosquito cells. Insect Mol Biol. 2007, 16: 389-399. 10.1111/j.1365-2583.2007.00734.x.
Cancino-Rodezno A, Alexander C, Villaseñor R, Pacheco S, Porta H, Pauchet Y, Soberón M, Gill SS, Bravo A: The mitogen-activated protein kinase p38 is involved in insect defense against Cry toxins from Bacillus thuringiensis. Insect Biochem Mol Biol. 2010, 40: 58-63. 10.1016/j.ibmb.2009.12.010.
Moon AE, Walker AJ, Goodbourn S: Regulation of transcription of the Aedes albopictus cecropin A1 gene: a role for p38 mitogen-activated protein kinase. Insect Biochem Mol Biol. 2011,
Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, Salzberg SL, Loftus B, Yandell M, Majoros WH, Rusch DB, Lai Z, Kraft CL, Abril JF, Anthouard V, Arensburger P, Atkinson PW, Baden H, de Berardinis V, Baldwin D, Benes V, Biedler J, Blass C, Bolanos R, Boscus D, Barnstead M, et al: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298: 129-149. 10.1126/science.1076181.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Cobb MH, Goldsmith EJ: How MAP kinases are regulated. J Biol Chem. 1995, 270: 14843-14846. 10.1074/jbc.270.25.14843.
Zhang J, Zhang F, Ebert D, Cobb MH, Goldsmith EJ: Activity of the MAP kinase ERK2 is controlled by a flexible surface loop. Structure. 1995, 3: 299-307. 10.1016/S0969-2126(01)00160-5.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.
Felsenstein J: Mathematics vs. evolution: mathematical evolutionary theory. Science. 1989, 246: 941-942. 10.1126/science.246.4932.941.
Müller HM, Dimopoulos G, Blass C, Kafatos FC: A hemocyte-like cell line established from the malaria vector Anopheles gambiae expresses six prophenoloxidase genes. J Biol Chem. 1999, 274: 11727-11735. 10.1074/jbc.274.17.11727.
Meredith JM, Munks RJ, Grail W, Hurd H, Eggleston P, Lehane MJ: A novel association between clustered NF-kappaB and C/EBP binding sites is required for immune regulation of mosquito defensin genes. Insect Mol Biol. 2006, 15: 393-401. 10.1111/j.1365-2583.2006.00635.x.
Lawson D, Arensburger P, Atkinson P, Besansky NJ, Bruggner RV, Butler R, Campbell KS, Christophides GK, Christley S, Dialynas E, Hammond M, Hill CA, Konopinski N, Lobo NF, MacCallum RM, Madey G, Megy K, Meyer J, Redmond S, Severson DW, Stinson EO, Topalis P, Birney E, Gelbart WM, Kafatos FC, Louis C, Collins FH: VectorBase: a data resource for invertebrate vector genomics. Nucleic Acids Res. 2009, D583-D587. VectorBase http://www.vectorbase.orgAnopheles gambiae PEST annotation, AgamP3.5, 37 Database
Keshet Y, Seger R: The MAP kinase signaling cascades: a system of hundreds of components regulates a diverse array of physiological functions. Methods Molec Biol. 2010, 661: 3-38. 10.1007/978-1-60761-795-2_1.
Champion A, Picaud A, Henry Y: Reassessing the MAP3K and MAP4K relationships. Trends Plant Sci. 2004, 9: 123-129. 10.1016/j.tplants.2004.01.005.
Rao KP, Richa T, Kumar K, Raghuram B, Sinha AK: In silico analysis reveals 75 members of mitogen-activated protein kinase kinase kinase gene family in rice. DNA Res. 2010, 17: 139-153. 10.1093/dnares/dsq011.
Mizutani T, Kobayashi M, Eshita Y, Shirato K, Kimura T, Ako Y, Miyoshi H, Takasaki T, Kurane I, Kariwa H, Umemura T, Takashima I: Involvement of the JNK-like protein of the Aedes albopictus mosquito cell line, C6/36, in phagocytosis, endocytosis and infection of West Nile virus. Insect Mol Biol. 2003, 12: 491-499. 10.1046/j.1365-2583.2003.00435.x.
Graves PR, Haystead TA: A functional proteomics approach to signal transduction. Recent Prog Horm Res. 2003, 58: 1-24. 10.1210/rp.58.1.1.
KEGG (Kyoto Encyclopedia of Genes and Genomes), Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38: D355-D360. 10.1093/nar/gkp896. Pathway Database http://www.genome.jp/kegg/pathway.html Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res 2006, 34:D354-D357; Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res 2000, 28:27-30.
Shivers RP, Youngman MJ, Kim DH: Transcriptional responses to pathogens in Caenorhabditis elegans. Curr Opin Microbiol. 2008, 11: 251-256. 10.1016/j.mib.2008.05.014.
Delaney JR, Stöven S, Uvell H, Anderson KV, Engström Y, Mlodzik M: Cooperative control of Drosophila immune responses by the JNK and NF-kappaB signaling pathways. EMBO J. 2006, 25: 3068-3077. 10.1038/sj.emboj.7601182.
Chen J, Xie C, Tian L, Hong L, Wu X, Han J: Participation of the p38 pathway in Drosophila host defense against pathogenic bacteria and fungi. Proc Natl Acad Sci USA. 2010, 107: 20774-20779. 10.1073/pnas.1009223107.
Park JM, Brady H, Ruocco MG, Sun H, Williams D, Lee SJ, Kato T, Richards N, Chan K, Mercurio F, Karin M, Wasserman SA: Targeting of TAK1 by the NF-kappa B protein Relish regulates the JNK-mediated immune response in Drosophila. Genes Dev. 2004, 18: 584-594. 10.1101/gad.1168104.
Brandt SM, Dionne MS, Khush RS, Pham LN, Vigdal TJ, Schneider DS: Secreted bacterial effectors and host-produced eiger/TNF drive death in a Salmonella-infected fruit fly. PLoS Biol. 2004, 2: e418-10.1371/journal.pbio.0020418.
Schneider DS, Ayres JS, Brandt SM, Costa A, Dionne MS, Gordon MD, Mabery EM, Moule MG, Pham LN, Shirasu-Hiza MM: Drosophila eiger mutants are sensitive to extracellular pathogens. PLoS Pathog. 2007, 3: e41-10.1371/journal.ppat.0030041.
Bidla G, Dushay MS, Theopold U: Crystal cell rupture after injury in Drosophila requires the JNK pathway, small GTPases and the TNF homolog Eiger. J Cell Sci. 2007, 120: 1209-1215. 10.1242/jcs.03420.
Giorgetti-Peraldi S, Peyrade F, Baron V, Van Obberghen E: Involvement of Janus kinases in the insulin signaling pathway. Eur J Biochem. 1995, 234: 656-660. 10.1111/j.1432-1033.1995.656_b.x.
Maegawa H, Kashiwagi A, Fujita T, Ugi S, Hasegawa M, Obata T, Nishio Y, Kojima H, Hidaka H, Kikkawa R: SHPTP2 serves adapter protein linking between Janus kinase 2 and insulin receptor substrates. Biochem Biophys Res Commun. 1996, 228: 122-127. 10.1006/bbrc.1996.1626.
Gual P, Baron V, Lequoy V, Van Obberghen E: Interaction of Janus kinases JAK-1 and JAK-2 with the insulin receptor and the insulin-like growth factor-1 receptor. Endocrinology. 1998, 139: 884-893. 10.1210/en.139.3.884.
Sandberg EM, Sayeski PP: Jak2 tyrosine kinase mediates oxidative stress-induced apoptosis in vascular smooth muscle cells. J Biol Chem. 2004, 279: 34547-34552. 10.1074/jbc.M405045200.
Thirone AC, JeBailey L, Bilan PJ, Klip A: Opposite effect of JAK2 on insulin-dependent activation of mitogen-activated protein kinases and Akt in muscle cells: possible target to ameliorate insulin resistance. Diabetes. 2006, 55: 942-951. 10.2337/diabetes.55.04.06.db05-1265.
Han YS, Thompson J, Kafatos FC, Barillas-Mury C: Molecular interactions between Anopheles stephensi midgut cells and Plasmodium berghei: the time bomb theory of ookinete invasion of mosquitoes. EMBO J. 2000, 19: 6030-6040. 10.1093/emboj/19.22.6030.
Shiao SH, Whitten MM, Zachary D, Hoffmann JA, Levashina EA: Fz2 and cdc42 mediate melanization and actin polymerization but are dispensable for Plasmodium killing in the mosquito midgut. PLoS Pathog. 2006, 2: e133-10.1371/journal.ppat.0020133.
Horton AA, Lee Y, Coulibaly CA, Rashbrook VK, Cornel AJ, Lanzaro GC, Luckhart S: Identification of three single nucleotide polymorphisms in Anopheles gambiae immune signaling genes that are associated with natural Plasmodium falciparum infection. Malar J. 2010, 9: 160-10.1186/1475-2875-9-160.
Riehle MM, Markianos K, Niaré O, Xu J, Li J, Touré AM, Podiougou B, Oduol F, Diawara S, Diallo M, Coulibaly B, Ouatara A, Kruglyak L, Traoré SF, Vernick KD: Natural malaria infection in Anopheles gambiae is regulated by a single genomic control region. Science. 2006, 312: 577-579. 10.1126/science.1124153.
Vodovotz Y, Constantine G, Faeder J, Mi Q, Rubin J, Bartels J, Sarkar J, Squires RH, Okonkwo DO, Gerlach J, Zamora R, Luckhart S, Ermentrout B, An G: Translational systems approaches to the biology of inflammation and healing. Immunopharmacol Immunotoxicol. 2010, 32: 181-195. 10.3109/08923970903369867.
Acknowledgements and funding
We acknowledge support provided by Hannah Smithers, Kong Wai Cheung and Andrew Ross for technical support for cell culture experiments and Dr. Ian Korf for computational support. This work was funded by the National Institutes of Health National Institute of Allergy and Infectious Diseases R01 AI080799 and R01 AI078183 to SL.
AAH and MSP confirmed the identities of the annotated and unannotated An. gambiae MAPKs, co-wrote the initial draft, and assisted in the preparation of figures/files. AAH also edited the manuscript and revisions. BW performed the MAPK signaling and expression assays with An. gambiae cells, analyzed the data, and edited the manuscript and revisions. LC and SAN performed the phylogenetic analyses, prepared and edited the phylogram, and edited the manuscript. AA, MN and JRF created the MAPK signaling diagram based on the KEGG Pathway Database and current literature and edited the manuscript. SL conceived of this work, directed the studies, and co-wrote the initial and final drafts for publication. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1:. H. sapiens, D. melanogaster, C. elegans and C. intestinalis are listed by common names and by accession number. An. gambiae MAPKs are listed by Vectorbase AGAP numbers as well as by NCBI accession numbers and common names. Dashes indicate the absence of identified orthologs. N/A for An. gambiae MAP3K7 indicates no available NCBI accession number. (XLS 66 KB)
Additional file 2:. Alignment of An. gambiae (AGAP), Ae. aegypti (AAEL) and Cx. quinquefasciatus (CPIJ) p38 MAPKs. The Ae. aegypti and Cx. quinquefasciatus sequences were used to query the An. gambiae trace archives database to identify the probable N-terminal amino acids of An. gambiae p38 MAPK (bold, underlined). (PDF 12 KB)
Additional file 3: An. gambiae JNKa and Ae. aegypti and Cx. quinquefasciatus JNK orthologs. Alignment of the first 20 amino acids of Ae. aegypti (AAEL) and Cx. quinquefasciatus (CPIJ) JNK orthologs with upstream sequence from AGAP009460 was used to predict the start methionine and additional N-terminal amino acids (bold, underlined) of An. gambiae JNKa. (PDF 25 KB)
Additional file 4: Ae. aegypti and Cx. quinquefasciatus. Ae. aegypti and Cx. quinquefasciatus MAPK orthologs as predicted from the supercontig assemblies for genome sequences for these species. These sequences were excluded from the phylogeny and, therefore, the relationships of these orthologs with An. gambiae MAPKs cannot definitively be confirmed. (PDF 96 KB)
Additional file 5:LPS-induced MAPK phosphorylation in An. gambiae 4a3B cells. Cells were treated with 100 μg/ml LPS or an equivalent volume of PBS. Cells were collected at 5 min, 30 min, 1 h, 3 h, 6 h and 24 h after treatment. MAPK phosphorylation was examined by western blotting as described in the Methods. GAPDH levels provided an assessment of protein loading and were used to normalize corresponding phospho-MAPK levels. This figure is a representative of immunoblots from 6-9 independent experiments. (PDF 81 KB)
Additional file 6:MAPK transcript expression in insulin-stimulated An. gambiae 4a3B cells. Cells were treated with 1.7 μM human insulin or an equivalent volume of diluent and collected at 1 h, 3 h, 6 h and 24 h after treatment (n = 3 for control and treatment at each timepoint). Expression levels of ERK (A), p38 MAPK (B), JNKa (C), and JNKb (D) were analyzed by qPCR as described in the Methods. Expression of ribosomal protein S7, a housekeeping gene control, was used for normalization of treatment and control expression. Insulin-treated MAPK expression levels are shown as fold changes relative to the timepoint-matched control group levels. (PDF 35 KB)
About this article
Cite this article
Horton, A.A., Wang, B., Camp, L. et al. The mitogen-activated protein kinome from Anopheles gambiae: identification, phylogeny and functional characterization of the ERK, JNK and p38 MAP kinases. BMC Genomics 12, 574 (2011). https://doi.org/10.1186/1471-2164-12-574
- Mosquito Cell
- 4a3B Cell
- MAPK Signaling Cascade
- Immune Signaling Pathway
- Gambiae Genome