Analysis of bacteria-challenged wild silkmoth, Antheraea mylitta (lepidoptera) transcriptome reveals potential immune genes
© Gandhe et al; licensee BioMed Central Ltd. 2006
Received: 04 April 2006
Accepted: 21 July 2006
Published: 21 July 2006
In the recent years a strong resemblance has been observed between the insect immune system and the mammalian innate immune mechanisms suggesting their common origin. Among the insects, only the dipterans (Drosophila and various mosquito species) have been widely investigated for their immune responses towards diverse pathogens. In the present study we constructed and analysed the immune transcriptome of the lepidopteran Antheraea mylitta, an economically important Indian tasar silkmoth with a view to unravel the potential immune-related genes and pathways.
An expressed sequence tag (EST) library was constructed from mRNA obtained from fat bodies of A. mylitta larvae that had been challenged by infection with Escherichia coli cells. We identified 719 unique ESTs from a total of 1412 sequences so generated. A third of the transcriptome showed similarity with previously characterized immune-related genes that included both the known and putative immune genes. Of the four putative novel defence proteins (DFPs) annotated by PSI-BLAST three showed similarity to extracellular matrix proteins from vertebrates implicated in innate immunity, while the fourth was similar to, yet distinct from, the anti-microbial protein cecropin. Finally, we analysed the expression profiles of 15 potential immune-related genes, and the majority of them were induced more prominently with E. coli compared to Micrococcus luteus. We also identified several unknown proteins, some of which could have probable immune-related functions based on the results of the ProDom analysis.
The present study has identified many potential immune-related genes in A. mylitta some of which are vertebrate homologues and others are hitherto unreported putative defence proteins. Several genes were present as members of gene families, as has also been observed in other insect species.
Insects are evolutionarily successful organisms and occupy almost all habitats in nature. An efficient immune system is one of the attributes for this evolutionary success. However, unlike mammals, the insects lack an adaptive immune system. The insect immune response is comprised of cellular and humoral components. The former involves the action of haemocytes in phagocytosis of microbes, encapsulation of large pathogens and nodule formation  whereas the latter involves activation of prophenoloxidase cascade leading to melanisation of invading microorganisms  and synthesis of a battery of anti-microbial peptides .
Insect immunity is well studied in dipterans such as fruit flies and mosquito species [4–7]. Only limited information is available on genes induced on pathogen challenge in a few lepidopteran species that include the domesticated silkmoth, Bombyx mori , Cecropia moth,Hyalophora cecropia  and tobacco hornworm,Manduca sexta  and in these too the immune response pathways employed to combat pathogen infections remain to be fully characterised.
Abundant genetic resources are now available for B. mori, with a 9 X shotgun sequence coverage of its genome and more than 100,000 ESTs in dbEST (NCBI) [11–13]. With reference to insect immunity, the ESTs have been obtained from baculovirus-infected B. mori cultured cells and pupae, but no large scale information on bacteria-induced immune genes is as yet available.
In this study, we have constructed and analysed an immune transcriptome following bacterial challenge of the Indian tasar wild silkmoth, Antheraea mylitta, an economically important lepidopteran cultivated for silk production. Prior information on immune response genes in wild silkmoths is lacking except for a few peripheral studies. Two proteins from A. mylitta – a lysozyme protein, 3-D structure of which is elucidated  and a protease inhibitor have been characterized . We chose to examine the fat body transcriptome since it is a major immune organ in insects, analogous to the mammalian liver. We generated a total of 1412 ESTs, of which 31% could be ascribed to putative immune functions. We also validated the upregulation of a selected subset of genes from the immune transcriptome by semi-quantitative RT-PCR.
Results and discussion
Classification of EST functional class categories based on similarity searches with NCBI protein database.
No. of ESTs
Total no. of ESTs
Housekeeping genes and other genes not involved in immunity
Hypothetical/unknown insect proteins
Hypothetical/unknown non-insect proteins
Weak homology in NCBI*
No homology in NCBI
The sequences, which could not be assigned any function based on homology search in NCBI, were searched for conserved domains in ProDom database. Of the 260 clusters that had no matches with known proteins in NCBI BLAST, we could assign protein domain families to 196 clusters based on the ProDom search (see Additional data file 2). The remaining sequences did not show any hits in ProDom and should be further analysed by other specialized computational tools.
Signal peptide and transmembrane domains
We further screened the unannotated proteins (no hits or hits with hypothetical proteins in NCBI database) for the presence of signal peptide and absence of transmembrane domains. This characteristic of many of the immune proteins has been utilized to screen probable immune-related genes in large-scale transcriptome studies , although the reliability of this criterion to identify the immune related genes remains to be experimentally tested. In the present study, 25 out of 260 genes tested fulfilled these criteria and could be considered as potential immune-responsive genes. However, since all the 260 clusters checked are not all full-length sequences, it is possible that we may have missed some others amongst them that represent gene products harbouring signal peptide without transmembrane domain, and hence the actual number may be higher.
Putative immune proteins
Putative anti-microbial proteins
The abundant genes in A. mylitta transcriptome.
No of ESTs
Defense protein 1
Ribosomal protein S2, Attacin
Elongation factor 1α
The data generated in the present study conformed to the trend of anti-microbial genes being present as multiple gene families as observed in many other insect species and highlighted the essentiality of these genes in the organism . Another interesting finding with respect to immune response was the presence of a protein identical to seroin. This protein is reported to express in the silk glands of B. mori and is known to protect the cocoon from microbes . A homologue of anti-fungal protein, gallerimycin, previously characterized from Galleria mellonella and shown to be induced by LPS injection , was also found.
Putative pattern recognition receptors (PRRs)
PRRs in the insects bind to and detect pathogen associated molecular patterns (PAMPs) like lipopolysaccharide, peptidoglycan, β 1–3 glucan, lipotechoic acid etc [9, 25]. The putative PRRs identified in the present immune transcriptome analysis were hemolins, peptidoglycan recognition proteins (PGRPs), gram-negative binding proteins (GNBPs), lectins and mucins. A. mylitta fat body transcriptome revealed two different proteins resembling hemolin, an immune inducible protein implicated in insect immunity  and it would be worthwhile to study the function and specificities of the different proteins of hemolin family. Three types of PGRP-like proteins were found in the A. mylitta immune repertoire in the present study as compared to D. melanogaster where 13 PGRP genes are known to be involved in the activation of the various effector pathways in immunity .
Putative proteases and protease inhibitors
Various proteases and protease inhibitors regulate the diverse immune mechanisms like melanization, phagocytosis and induction of anti-microbial peptides . A homologue of a prophenoloxidase activating protease characterized in B. mori and M. sexta  was also found in the A. mylitta immune transcriptome. Several other classes of putative proteases like cysteine proteases, serine proteases and metallo-proteases were also identified. As many as thirteen distinct serpin-like (serine protease inhibitors) and twelve potential protease inhibitors were detected in the immune transcriptome. Five different serpins have earlier been identified in M. sexta and shown to differ in the induction pattern upon immune challenge . In the light of these studies, the information on serpins and protease inhibitors from the A. mylitta immune transcriptome will prove to be invaluable in further understanding of various immune pathways in insects.
Putative defence proteins with unknown function
Several potential new members of known protein families were identified in the present study. Among the new members, one was a putative lysozyme-like protein described in the previous section. An array of putative proteinase inhibitors and proteases were also found. Many of them are new, and their study would enhance our understanding of the mechanisms of proteolytic cascades in insect innate immunity. A few immunoglobulin (Ig) like molecules were identified by ProDom search (see Additional data file 2). Ig-like molecules- hemolin  and more recently Dscam have been implicated in insect immunity [28, 29] and it would be interesting to evaluate the role of these putative Ig-like molecules in insects. Among the potential immune proteins, we describe below in more detail four putative defence proteins (DFPs), for two of which (DFP-1 and DFP-4) we have confirmed the induction upon E. coli infection by semi-quantitative RT-PCR.
DFP-1, DFP-2 and DFP-3
Proteins similar to putative defence proteins (DFPs) revealed by PSI-BLAST analysis.
Novel proteins from A. mylitta
Description and function of proteins similar to DFPs.
DFP-1, DFP-2, DFP-3
1. Stromal cell derived factor receptor 2 homologue
Homo sapiens [GenBank:NM_001013660.1]
(catecholamine catabolism) reeler domain
Bos taurus [GenBank:Q9GLX9]
H. sapiens [GenBank:Q9HCB6]
Extracellular matrix, cell adhesion protein that promotes the attachment of spinal cord and sensory neuron cells
3. Reelin precursor
H. sapiens [GenBank:P78509]
Extracellular matrix serine protease. Enzymatic activity is important for the modulation of cell adhesion
Cecropin D precursor
L. obliqua [GenBank:AAV91462]
The immune response in insects is dynamic and different effector genes are likely expressed at different time points during infection, contributing to the ability of the insects to ward off infections in spite of the absence of adaptive immunity. The current transcriptome represents genes likely expressed upon E. coli infection in the A. mylitta fat body at 24 hrs post infection. Unexpectedly, the Imd pathway components that are implicated in the activation of various effector pathways upon gram-negative bacterial infection in Drosophila were not present in our transcriptome. The 24-hour post infection period may have been non-optimal for expression of some genes, and since the EST library was not normalized the less abundant transcripts may have gone undetected. Alternatively, it is possible that other pathways are involved in the immune induction in the moths.
The present study has increased the repertoire of lepidopteran-specific putative immune response genes by several hundred-fold. This will be a valuable resource for lepidopteran-specific immune studies in particular and insect immune studies in general.
Insects, bacterial inoculation and tissue collection
A. mylitta, 5th instar, day 3 larvae were procured from Regional Research Station, Warangal, Andhra Pradesh. Log phase E. coli cells (DH5α), washed and resuspended in saline (0.3 M NaCl, 0.005 M KCl), were injected into the haemocoel of the larvae as described earlier . At 24 hours post infection (hpi), larvae were dissected to isolate fat body, and the tissue was flash frozen in liquid N2 and then stored at -70°C till further use.
cDNA library construction and generation of ESTs
Total RNA was extracted from the fat body using Trizol reagent (Invitrogen). The complementary DNA synthesis was carried out using Stratagene ZAP-cDNA® synthesis kit following manufacturer's instructions. Directional cDNA library was constructed by cloning of cDNA fragments into pBluescript II SK (+) vector and electroporation into E. coli strain DH10B. Insert-containing plasmid clones were sequenced with RV-M primer (5'GAG CGG ATA ACA ATT TCA CAC AGG 3') with the aid of MegaBACE3000 sequencer.
Raw sequences obtained from sequence chromatograms were processed using several programs. A cut-off Phred Quality Value of ≥15 was assigned to extract quality sequences from chromatograms. The quality sequences were screened for the presence of vector sequences using 'Cross Match' program . Then masked vector sequences were automatically removed by in-house developed trimming tool. Sequences shorter than 50 bases were removed. The resulting high-quality sequences were assembled into sequence contigs with the TGICL program , which initially makes clusters using MegaBLAST and thereafter makes an assembly using CAP3 for each cluster generated in the first step. A cluster is defined as a unique sequence obtained either by multiple alignment of many sequences that are > 95% similar or derived from a single sequence. A cluster containing ≥2 ESTs is termed a contig and that containing only one sequence, a singleton.
The unique putative gene sequences obtained by clustering and assembly were annotated by running BLAST  against non-redundant (nr) protein database of NCBI. Further, BLAST output was parsed to classify the putative gene transcripts into different functional classes.
Analysis of unknown proteins
Proteins that did not show any significant hits in NCBI nr database or showed similarity to unknown or hypothetical proteins were characterized by additional computational tools.
a. Domain search
For the sequences showing high similarity to hypothetical and/or unknown proteins, and those showing weak similarity to nr protein database, domain search was performed using ProDom . Putative function was assigned based on the type of domains found.
b. Signal peptide and transmembrane domain analysis
Presence of transmembrane domains and signal peptide analysis was done on the transcripts not showing any significant hits in NCBI database. The signal peptide analysis was done by SignalP software  and trans-membrane domain analysis was done with TMHMM program .
c. Functional annotation with PSI-BLAST
The functional annotation of four novel immune upregulated transcripts (DFPs) was done using PSI-BLAST .
A. mylitta larvae were differentially challenged- a) Unchallenged, b) Saline-injected, c) M. luteus and d) E. coli. Log phase E. coli and M. luteus bacteria, washed and resuspended in insect saline (0.3 M NaCl, 0.005 M KCl) were injected (30 μl) into a set of A. mylitta larvae. One set each of saline-injected and uninjected larvae were kept as a control. Four tissues, namely fat body, epidermis, mid gut, and silk gland were dissected out and flash frozen in liquid nitrogen. Total RNA was isolated using Trizol reagent. To remove genomic DNA contamination, total RNA was treated with RNAse free DNAse (NEB) as prescribed by the manufacturers. cDNA was synthesized using MMLV reverse transcriptase (Invitrogen) and oligo dT primers from 1 μg of total RNA. The primers were designed for the selected ESTs by Primer3 software .
Semi-quantitative RT-PCR was carried out for all the four differentially challenged tissues using an Eppendorf master cycler under the following conditions- 94°C, 2 min- initial denaturation, 27 cycles (94°C – 30 s, 58°C- 30s, 72°C-2 mins) and a final elongation at 72°C for 10 mins. Actin cDNA was amplified as an endogenous control. PCR reaction components included: 1X buffer, 100μM dNTPs, 1.5 mM MgCl2, 0.5 units Taq polymerase (MBI), 0.5 μM primers. Primer sequences are enlisted in additional data (see Additional data file 5).
Obtaining full-length cDNA by 5' RACE
We had obtained full-length coding sequences of the DFP-1, 2 and 3 through EST sequencing. To acquire full-length DFP-4 cDNA, we carried out 5' RACE PCR using the 5' RACE kit (Clontech). The 5' ends were amplified by using an adaptor primer and a reverse gene specific primer. PCR was performed for 25 cycles in an Eppendorf master cycler. A 300 bp band was isolated, sequenced and confirmed to be the 5' DFP-4 sequences.
We acknowledge Ms Riti Mohapatra for her assistance in the bioinformatics analysis. We are thankful to Mr. Riccardo M. Bennett-Lovsey and Mr. Alex D. Herbert, graduate students, Imperial College, London, for their suggestions in annotating the unknown proteins. We also acknowledge Mr. Venkateswara Rao, Department of Sericulture, Andhra Pradesh, India for supplying A. mylitta larvae. JN acknowledges financial support from Department of Biotechnology, India. AG and AKP are recipients of the fellowship from University Grants Commission (UGC) and Council of Scientific and Industrial Research (CSIR), India respectively.
- Kanost MR, Jiang H, Yu XQ: Innate immune responses of a lepidopteran insect, Manduca sexta. Immunol Rev. 2004, 198: 97-105. 10.1111/j.0105-2896.2004.0121.x.PubMedView ArticleGoogle Scholar
- Cerenius L, Soderhall K: The prophenoloxidase-activating system in invertebrates. Immunol Rev. 2004, 198: 116-126. 10.1111/j.0105-2896.2004.00116.x.PubMedView ArticleGoogle Scholar
- Steiner H: Peptidoglycan recognition proteins: on and off switches for innate immunity. Immunol Rev. 2004, 198: 83-96. 10.1111/j.0105-2896.2004.0120.x.PubMedView ArticleGoogle Scholar
- Vodovar N, Vinals M, Liehl P, Basset A, Degrouard J, Spellman P, Boccard F, Lemaitre B: Drosophila host defense after oral infection by an entomopathogenic Pseudomonas species. Proc Natl Acad Sci U S A. 2005, 102 (32): 11414-11419. 10.1073/pnas.0502240102.PubMedPubMed CentralView ArticleGoogle Scholar
- Johansson KC, Metzendorf C, Soderhall K: Microarray analysis of immune challenged Drosophila hemocytes. Exp Cell Res. 2005, 305 (1): 145-155. 10.1016/j.yexcr.2004.12.018.PubMedView ArticleGoogle Scholar
- Xu X, Dong Y, Abraham EG, Kocan A, Srinivasan P, Ghosh AK, Sinden RE, Ribeiro JM, Jacobs-Lorena M, Kafatos FC, Dimopoulos G: Transcriptome analysis of Anopheles stephensi-Plasmodium berghei interactions. Mol Biochem Parasitol. 2005, 142 (1): 76-87. 10.1016/j.molbiopara.2005.02.013.PubMedView ArticleGoogle Scholar
- Aguilar R, Jedlicka AE, Mintz M, Mahairaki V, Scott AL, Dimopoulos G: Global gene expression analysis of Anopheles gambiae responses to microbial challenge. Insect Biochem Mol Biol. 2005, 35 (7): 709-719. 10.1016/j.ibmb.2005.02.019.PubMedView ArticleGoogle Scholar
- Okano K, Shimada T, Mita K, Maeda S: Comparative expressed-sequence-tag analysis of differential gene expression profiles in BmNPV-infected BmN cells. Virology. 2001, 282 (2): 348-356. 10.1006/viro.2001.0823.PubMedView ArticleGoogle Scholar
- Hultmark D, Engstrom A, Andersson K, Steiner H, Bennich H, Boman HG: Insect immunity. Attacins, a family of antibacterial proteins from Hyalophora cecropia. Embo J. 1983, 2 (4): 571-576.PubMedPubMed CentralGoogle Scholar
- Zhu Y, Johnson TJ, Myers AA, Kanost MR: Identification by subtractive suppression hybridization of bacteria-induced genes expressed in Manduca sexta fat body. Insect Biochem Mol Biol. 2003, 33 (5): 541-559. 10.1016/S0965-1748(03)00028-6.PubMedView ArticleGoogle Scholar
- Wang J, Xia Q, He X, Dai M, Ruan J, Chen J, Yu G, Yuan H, Hu Y, Li R, Feng T, Ye C, Lu C, Wang J, Li S, Wong GK, Yang H, Wang J, Xiang Z, Zhou Z, Yu J: SilkDB: a knowledgebase for silkworm biology and genomics. Nucleic Acids Res. 2005, 33 (Database issue): D399-402. 10.1093/nar/gki116.PubMedPubMed CentralView ArticleGoogle Scholar
- Mita K, Kasahara M, Sasaki S, Nagayasu Y, Yamada T, Kanamori H, Namiki N, Kitagawa M, Yamashita H, Yasukochi Y, Kadono-Okuda K, Yamamoto K, Ajimura M, Ravikumar G, Shimomura M, Nagamura Y, Shin IT, Abe H, Shimada T, Morishita S, Sasaki T: The genome sequence of silkworm, Bombyx mori. DNA Res. 2004, 11 (1): 27-35. 10.1093/dnares/11.1.27.PubMedView ArticleGoogle Scholar
- Mita K, Morimyo M, Okano K, Koike Y, Nohata J, Kawasaki H, Kadono-Okuda K, Yamamoto K, Suzuki MG, Shimada T, Goldsmith MR, Maeda S: The construction of an EST database for Bombyx mori and its application. Proc Natl Acad Sci U S A. 2003, 100 (24): 14121-14126. 10.1073/pnas.2234984100.PubMedPubMed CentralView ArticleGoogle Scholar
- Jain D, Nair DT, Swaminathan GJ, Abraham EG, Nagaraju J, Salunke DM: Structure of the induced antibacterial protein from tasar silkworm, Antheraea mylitta. Implications to molecular evolution. J Biol Chem. 2001, 276 (44): 41377-41382. 10.1074/jbc.M104674200.PubMedView ArticleGoogle Scholar
- Shrivastava B, Ghosh AK: Protein purification, cDNA cloning and characterization of a protease inhibitor from the Indian tasar silkworm, Antheraea mylitta. Insect Biochem Mol Biol. 2003, 33 (10): 1025-1033. 10.1016/S0965-1748(03)00117-6.PubMedView ArticleGoogle Scholar
- 719 clusters of A. mylitta. [http://188.8.131.52:9999/PHP/SILKSAT/uniqueseqs/]Google Scholar
- Roxstrom-Lindquist K, Terenius O, Faye I: Parasite-specific immune response in adult Drosophila melanogaster: a genomic study. EMBO Rep. 2004, 5 (2): 207-212. 10.1038/sj.embor.7400073.PubMedPubMed CentralView ArticleGoogle Scholar
- Irving P, Troxler L, Heuer TS, Belvin M, Kopczynski C, Reichhart JM, Hoffmann JA, Hetru C: A genome-wide analysis of immune responses in Drosophila. Proc Natl Acad Sci U S A. 2001, 98 (26): 15119-15124. 10.1073/pnas.261573998.PubMedPubMed CentralView ArticleGoogle Scholar
- Gudmundsson GH, Lidholm DA, Asling B, Gan R, Boman HG: The cecropin locus. Cloning and expression of a gene cluster encoding three antibacterial peptides in Hyalophora cecropia. J Biol Chem. 1991, 266 (18): 11510-11517.PubMedGoogle Scholar
- Ochiai M, Ashida M: A pattern recognition protein for peptidoglycan. Cloning the cDNA and the gene of the silkworm, Bombyx mori. J Biol Chem. 1999, 274 (17): 11854-11858. 10.1074/jbc.274.17.11854.PubMedView ArticleGoogle Scholar
- Hultmark D, Steiner H, Rasmuson T, Boman HG: Insect immunity. Purification and properties of three inducible bactericidal proteins from hemolymph of immunized pupae of Hyalophora cecropia. Eur J Biochem. 1980, 106 (1): 7-16.PubMedView ArticleGoogle Scholar
- Cheng T, Zhao P, Liu C, Xu P, Gao Z, Xia Q, Xiang Z: Structures, regulatory regions, and inductive expression patterns of antimicrobial peptide genes in the silkworm Bombyx mori. Genomics. 2006, 87 (3): 356-365. 10.1016/j.ygeno.2005.11.018.PubMedView ArticleGoogle Scholar
- Nirmala X, Mita K, Vanisree V, Zurovec M, Sehnal F: Identification of four small molecular mass proteins in the silk of Bombyx mori. Insect Mol Biol. 2001, 10 (5): 437-445. 10.1046/j.0962-1075.2001.00282.x.PubMedView ArticleGoogle Scholar
- Schuhmann B, Seitz V, Vilcinskas A, Podsiadlowski L: Cloning and expression of gallerimycin, an antifungal peptide expressed in immune response of greater wax moth larvae, Galleria mellonella. Arch Insect Biochem Physiol. 2003, 53 (3): 125-133. 10.1002/arch.10091.PubMedView ArticleGoogle Scholar
- Hoffmann JA: The immune response of Drosophila. Nature. 2003, 426 (6962): 33-38. 10.1038/nature02021.PubMedView ArticleGoogle Scholar
- Sun SC, Lindstrom I, Boman HG, Faye I, Schmidt O: Hemolin: an insect-immune protein belonging to the immunoglobulin superfamily. Science. 1990, 250 (4988): 1729-1732.PubMedView ArticleGoogle Scholar
- Dziarski R: Peptidoglycan recognition proteins (PGRPs). Mol Immunol. 2004, 40 (12): 877-886. 10.1016/j.molimm.2003.10.011.PubMedView ArticleGoogle Scholar
- Watson FL, Puttmann-Holgado R, Thomas F, Lamar DL, Hughes M, Kondo M, Rebel VI, Schmucker D: Extensive diversity of Ig-superfamily proteins in the immune system of insects. Science. 2005, 309 (5742): 1874-1878. 10.1126/science.1116887.PubMedView ArticleGoogle Scholar
- Dong Y, Taylor HE, Dimopoulos G: AgDscam, a Hypervariable Immunoglobulin Domain-Containing Receptor of theAnopheles gambiae Innate Immune System. PLoS Biol. 2006, 4 (7): e229-10.1371/journal.pbio.0040229.PubMedPubMed CentralView ArticleGoogle Scholar
- Veiga AB, Ribeiro JM, Guimaraes JA, Francischetti IM: A catalog for the transcripts from the venomous structures of the caterpillar Lonomia obliqua: identification of the proteins potentially involved in the coagulation disorder and hemorrhagic syndrome. Gene. 2005, 355: 11-27. 10.1016/j.gene.2005.05.002.PubMedPubMed CentralView ArticleGoogle Scholar
- Shin SW, Park SS, Park DS, Kim MG, Kim SC, Brey PT, Park HY: Isolation and characterization of immune-related genes from the fall webworm, Hyphantria cunea, using PCR-based differential display and subtractive cloning. Insect Biochem Mol Biol. 1998, 28 (11): 827-837. 10.1016/S0965-1748(98)00077-0.PubMedView ArticleGoogle Scholar
- Bao Y, Mega K, Yamano Y, Morishima I: cDNA cloning and expression of bacteria-induced Hdd11 gene from eri-silkworm, Samia cynthia ricini. Comp Biochem Physiol C Toxicol Pharmacol. 2003, 136 (4): 337-342. 10.1016/j.cca.2003.10.007.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- He YW, Li H, Zhang J, Hsu CL, Lin E, Zhang N, Guo J, Forbush KA, Bevan MJ: The extracellular matrix protein mindin is a pattern-recognition molecule for microbial pathogens. Nat Immunol. 2004, 5 (1): 88-97. 10.1038/ni1021.PubMedView ArticleGoogle Scholar
- Zhang XF, Wang JF, Matczak E, Proper JA, Groopman JE: Janus kinase 2 is involved in stromal cell-derived factor-1alpha-induced tyrosine phosphorylation of focal adhesion proteins and migration of hematopoietic progenitor cells. Blood. 2001, 97 (11): 3342-3348. 10.1182/blood.V97.11.3342.PubMedView ArticleGoogle Scholar
- De Gregorio E, Spellman PT, Rubin GM, Lemaitre B: Genome-wide analysis of the Drosophila immune response by using oligonucleotide microarrays. Proc Natl Acad Sci U S A. 2001, 98 (22): 12590-12595. 10.1073/pnas.221458698.PubMedPubMed CentralView ArticleGoogle Scholar
- Hedengren-Olcott M, Olcott MC, Mooney DT, Ekengren S, Geller BL, Taylor BJ: Differential activation of the NF-kappaB-like factors Relish and Dif in Drosophila melanogaster by fungi and Gram-positive bacteria. J Biol Chem. 2004, 279 (20): 21121-21127. 10.1074/jbc.M313856200.PubMedView ArticleGoogle Scholar
- Abraham EG, Nagaraju J, Salunke D, Gupta HM, Datta RK: Purification and partial characterization of an induced antibacterial protein in the silkworm, Bombyx mori. J Invertebr Pathol. 1995, 65 (1): 17-24. 10.1006/jipa.1995.1003.PubMedView ArticleGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.PubMedView ArticleGoogle Scholar
- The TGICL program . [www.tigr.org]Google Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1006/jmbi.1990.9999.PubMedView ArticleGoogle Scholar
- Bru C, Courcelle E, Carrere S, Beausse Y, Dalmar S, Kahn D: The ProDom database of protein domain families: more emphasis on 3D. Nucleic Acids Res. 2005, 33 (Database issue): D212-5. 10.1093/nar/gki034.PubMedPubMed CentralView ArticleGoogle Scholar
- SignalP software. [http://www.cbs.dtu.dk/services/SignalP/]Google Scholar
- TMHMM program. http://www.cbs.dtu.dk/services/TMHMM-2.0/Google Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.