Skip to main content
  • Research article
  • Open access
  • Published:

Identification of the Weevil immune genes and their expression in the bacteriome tissue

Abstract

Background

Persistent infections with mutualistic intracellular bacteria (endosymbionts) are well represented in insects and are considered to be a driving force in evolution. However, while pathogenic relationships have been well studied over the last decades very little is known about the recognition of the endosymbionts by the host immune system and the mechanism that limits their infection to the bacteria-bearing host tissue (the bacteriome).

Results

To study bacteriome immune specificity, we first identified immune-relevant genes of the weevil Sitophilus zeamais by using suppressive subtractive hybridization (SSH) and then analyzed their full-length coding sequences obtained by RACE-PCR experiments. We then measured immune gene expression in the bacteriome, and in the aposymbiotic larvae following S. zeamais primary endosymbiont (SZPE) injection into the hemolymph, in order to consider the questions of bacteriome immune specificity and the insect humoral response to symbionts. We show that larval challenge with the endosymbiont results in a significant induction of antibacterial peptide genes, providing evidence that, outside the bacteriome, SZPE are recognized as microbial intruders by the host. In the bacteriome, gene expression analysis shows the overexpression of one antibacterial peptide from the coleoptericin family and, intriguingly, homologs to genes described as immune modulators (that is, PGRP-LB, Tollip) were also shown to be highly expressed in the bacteriome.

Conclusion

The current data provide the first description of immune gene expression in the insect bacteriome. Compared with the insect humoral response to SZPE, the bacteriome expresses few genes among those investigated in this work. This local immune gene expression may help to maintain the endosymbiont in the bacteriome and prevent its invasion into insect tissues. Further investigations of the coleoptericin, the PGRP and the Tollip genes should elucidate the role of the host immune system in the maintenance and regulation of endosymbiosis.

Background

Chronic bacterial infections are widespread in nature and exhibit a large range of interactions with their host, from mutualism to parasitism. In insects, symbiotic intracellular bacteria (endosymbionts) are deeply integrated into host cell biology and development as they are transmitted maternally through hundreds of host generations, and early on in the insect embryogenesis they invade specialized host cells called bacteriocytes that sometimes form a specific organ, the bacteriome [13]. Physiological and molecular investigations have provided evidence that endosymbionts supply their host's diet with limiting nutrients, thereby improving their invasive power and their ability to settle on nutritionally poor sources and habitats, such as blood (Glossina, Rhodnius), plant sap (aphids, psyllids, whiteflies, mealybugs), and cereal grains (Sitophilus) [2, 4].

However, while the physiological and the evolutionary aspects of insect endosymbiosis have been thoroughly investigated over the past two decades, very little is known about the molecular mechanisms that permit the establishment, then the maintenance and the regulation of such beneficial interactions. One striking question concerns the interaction between the bacteria and the host innate immune system, an area which has been relatively well investigated in pathogenic relationships compared with mutualistic associations that have been recently approached in only a few systems [58].

To combat infection, insects rely on multiple innate defense reactions. Insect immunity includes physical barriers, together with local and systemic immune responses involving both cellular and humoral pathways (reviewed in [9]). Activation of the humoral pathway consists of microbe-associated molecular pattern recognition by pattern recognition receptors, such as peptidoglycan recognition proteins (PGRPs), and the activation of intracellular signaling pathways, such as the Toll and the Immune deficiency (Imd) pathways. These pathways activate, in particular, the production and the secretion of a panel of antimicrobial peptides (AMPs) in response to Gram-positive and Gram-negative bacteria. In addition to AMP production, the insect humoral response also involves a proteolytic cascade leading to prophenoloxidase (PPO) activation and subsequent melanin synthesis at the site of cuticular injury. This reaction, called melanization, plays a key role in wound healing, encapsulation, sequestration of microbes and production of toxic intermediates [10]. While the systemic response is by far the best characterized pathway among immune reactions, the local immune response (also known as epithelial immunity) was only recently shown to significantly contribute to protection against invading microorganisms in the alimentary tract and tracheae. In the gut, for example, there is an inducible local production of AMPs. This response is triggered upon natural infection by Gram-negative bacteria and is mediated by the Imd pathway [11]. The gut epithelium also expresses amidase PGRP. It has been shown that these PGRPs, which scavenge peptidoglycan released by commensal bacteria, reduce gut immune reactivity and avoid a state of permanent immune activation in this tissue [12, 13].

In a previous work on the mutualistic interaction between the weevil Sitophilus zeamais and a γ-proteobacterium called S. zeamais primary endosymbiont (SZPE) (see [14] for a review on the model), we discovered the overexpression of a member of the PGRP gene family in the bacteriome tissue of the host [15, 16]. We showed that wPGRP is induced after a bacterial challenge and that wPGRP gene expression depends on bacterial growth. Moreover, we have shown that the wPGRP gene was induced concomitantly with an endosymbiont release from the bacteriocytes during nymphosis, demonstrating that the PGRP gene family is involved in host-symbiont interaction [16].

This study was dedicated to enlarging the panel of insect immune genes involved in the host-symbiont interaction. We first identified the host immune-relevant genes by suppressive subtractive hybridization (SSH), and analyzed in silico their full-length coding sequences completed with RACE-PCR. We then measured their steady-state levels in aposymbiotic larvae challenged with SZPE and in the bacteriocyte cells of naturally infected larvae. We show that experimental infection of larvae with SZPE results in a significant induction of AMPs and that only few immune gene transcripts, including homologs to an antibacterial peptide and immunomodulators, are accumulated in the bacteriome tissue. These data reveal that endosymbionts are perceived as intruders while being present in hemolymph and that there is a local immune gene expression in the bacteriome. This study provides the first indication of how insects may maintain endosymbionts within the bacteriome and prevent their invasion into insect tissues.

Results

Identification of immune-relevant genes in Sitophilus zeamais

As the Sitophilus genome has not been sequenced, we applied SSH technology to cDNA from E. coli-infected larvae and cDNA from naive larvae to identify the immune genes of interest to this work. To obtain genes expressed at different phases of the immune response, three RNA samples were extracted 3, 6 and 12 hours after Escherichia coli infection and mixed prior to cDNA synthesis. We sequenced 485 expressed sequence tags (ESTs) from the subtracted library. Following quality analysis, trimming and chimeric sequence digestion (see Materials and Methods), we assembled 475 sequences into 273 putative transcripts, consisting of 62 contigs and 211 singletons (Table 1).

Table 1 General characteristics of Sitophilus zeamais ESTs from suppressive subtractive hybridization between E. coli-infected and naive larvae

To gain insight into the function of EST products, we compared the 273 resulting sequences with the UniProt databases and classified them using the Gene Ontology (GO) scheme. Moreover, we constructed a database where all the ESTs and a complete list of BlastX matches, together with their functional classification, can be found. Of these 273 sequences, 83% successfully matched against the UniProt database and 63% were functionally classified (Table 1). Concerning the immune genes, around 12% of the ESTs have similarities with transcripts that are known to encode proteins involved in an immune function. According to sequence identity, some ESTs showed only a weak similarity to antibacterial peptides, such as diptericin A, cecropin A1, sarcotoxin II-1 and luxuriosin. Nevertheless, taking into account their high redundancy (12, 4, 15, and 4 copies, respectively), we have included them in this study. In addition, sequence analysis of ESTs with similarities to diptericin, tenecin and luxuriocin has uncovered various isoforms of these peptides (see Additional file 1). However, due to the high DNA and protein sequence identities (76% to 97% and 91% to 100%, respectively) we have analyzed only one isoform of each gene family in this work.

Analysis of full-length cDNA of immune genes obtained by RACE-PCR

To improve sequence similarities, we applied the RACE-PCR technology and have sequenced the full cDNA sequences of ESTs (one EST for each cluster) with similarities to antibacterial peptides, lysozymes, PGRPs and Tollip (Figure 1). We also analyzed the cDNA sequences in silico to identify conserved protein domains and to predict the cellular localization of the peptides and proteins.

Figure 1
figure 1

Schematic representations of putative Sitophilus zeamais proteins with similarity to antibacterial peptides, lysozymes, PGRP and Tollip. For each gene, the whole cDNA sequence was obtained from the corresponding EST by RACE-PCR and was then confirmed by whole cDNA amplification and sequencing. For each cDNA, ORF was predicted using the MacMolly software package. The top scale shows the length of the various domains of the proteins (aa, amino acid). Black regions indicate the predicted signal peptide (TargetP) and gray regions the putative propeptide domain according to the conserved R-x-(K/R)-R motif. mcs: minimal cleavage site corresponding to an R-x-x-R motif with an additional arginine in P6 position, which may enhance cleavage. Regions with similarity to conserved domains detected by InterProScan are indicated as hatched regions with the associated E-value. Coleoptericin, IPR009382; Defensin, IPR001542; Kunitz, IPR002223; Lysozyme, IPR000974; Destabilase, IPR008597; PGRP, IPR002502; C2, IPR000008; CUE, IPR003892. The accession numbers of the complete coding sequences are: INF-18, EU282111; INF-145, EU282117; INF-42, EU282112; INF-163, EU282118; INF-165, EU282113; INF-479, EU282119; INF-217, EU282115; INF-282, EU282114; INF-152, EU282120; INF-441, EU282121; INF-9, EU282122; INF-359, EU282116.

Antimicrobial peptides

AMPs typically contain fewer than 100 amino acid residues, divided in most cases into three domains: a signal peptide, a propeptide (R-x-(K/R)-R motifs), and the mature peptide. The first two domains are proteolytically cleaved to release the mature peptide active against bacteria [1719].

A signal peptide and an R-x-(K/R)-R motif have been predicted for all the putative antibacterial peptides of S. zeamais, except for INF-163 and INF-479, for which only a signal peptide is reliably predicted (Figure 1). According to protein similarity, we confirmed the identification of two AMPs (INF-18 and INF-145, 40% identity) belonging to the coleoptericin family, a defensin with sequence identity of 66% (32/48) to tenecin (INF-217), and a homolog of luxuriosin, an AMP characterized by a Kunitz domain [20], with a sequence identity of 42% (45/106). For ESTs INF-42, INF-163 and INF-165, the full-length transcripts still show weak similarities to AMPs and no characteristic domain has been identified for the putative peptides, which cannot be considered as antibacterial peptides without further analysis (Figure 1).

Lysozymes

Lysozymes are widespread enzymes characterized by their ability to cleave bacterial peptidoglycans. Among lysozymes, several types have been described according to the genomic structure and phylogenetic data: the c (chicken), g (goose), i (invertebrate), phage, bacteria and plant types [21]. Here, according to sequence homology, we have identified two lysozymes: an i-type (INF-152), with a destabilase-related domain identified in an i-type lysozyme from the leech [22], and a c-type lysozyme (INF-282) with a characteristic lysozyme domain. Both proteins have a predicted signal peptide.

PGRP

Two PGRP transcripts were identified in the subtracted library. One EST (INF-9) corresponds to the wPGRP gene previously identified [15, 16] while the EST INF-441 is an additional weevil PGRP gene. These two genes will be referenced as wPGRP1 and wPGRP2, respectively. RACE-PCR product analysis has shown that the two wPGRP proteins have 30% identity and have both conserved the residues necessary for amidase activity [23]. According to in silico analysis, wPGRP1 may be an intracellular amidase while wPGRP2 may be secreted in hemolymph.

Tollip (Toll interacting protein)

Among the ESTs of the subtracted library, we found a sequence (INF-359) with homology to Tollip, a regulator of the Toll-like receptors pathway described in mammals [2426]. Analysis of the RACE-PCR product has confirmed homology to the mammal genes since the predicted protein of S. zeamais shares 46% identity (127/275) with the Tollip protein from Mus musculus.

Gene expression analysis after sterile or septic injury of weevil aposymbiotic larvae

To determine genes induced by wounding (pricking effect) and genes specifically induced by bacterial challenge, using qRT-PCR we have compared the transcript levels in aposymbiotic larvae after sterile and septic injury, and in naive aposymbiotic larvae taken as a control (Table 2).

Table 2 Immune-related ESTs and comparison of gene expression between naive larvae, mock-infected larvae and larvae challenged with Escherichia coli.

According to the Kruskal-Wallis rank sum test, most of the genes with sequence similarity to AMPs are highly induced (30 to 300-fold) after E. coli infection, including genes without any significant antibacterial domain (Table 2 and Figure 1). Some peptide genes, such as INF-18, INF-42 and INF-217, are also slightly induced (10-fold) after a sterile injury. INF-479, the homolog of luxuriosin, is the only peptide gene induced after sterile injury but not in response to E. coli challenge. Finally, both lysozyme genes (INF-152 and INF-282) are induced after sterile injury, independently of the bacterial infection.

The qRT-PCR data show that both wPGRP genes are upregulated. However, while wPGRP1 (INF-9) is weakly induced by injury (2.3-fold) and more strongly induced after E. coli infection (6.7-fold), wPGRP2 (INF-441) is induced by sterile injury only (7.8 to 11.6-fold).

As regulators we quantified the expression of genes with similarity to proteases and protease inhibitors, in addition to the Tollip gene. The Tollip gene (INF-359) and two genes with similarity to proteinases (INF-459 and INF-515) were shown to be constitutively expressed. All the other genes were induced after sterile injury, except the cysteine-rich venom-like proteinase homolog (INF-91) that was induced after infection with E. coli. However, no data are available concerning the function of this protein identified in the salivary gland of Aedes albopictus.

In addition to humoral immune response genes we quantified two cytoskeleton genes, as the participation of actin cytoskeleton regulation proteins in innate immunity has been established by functional genomic analysis of phagocytosis [27]. However, no variation in the transcript levels has been observed, either for the actin gene or for the homolog of profilin, an actin polymerization regulator (INF-13).

Gene expression analysis in aposymbiotic larvae challenged with SZPE

To examine immune response to SZPE while present in the larval hemolymph, we quantified gene transcript levels in aposymbiotic larvae following an injection of SZPE. Knowing that SZPE fails to divide in vitro, we injected approximately 1 × 105 viable or heat-killed SZPE (Figure 2), the amount of E. coli required to induce immune response (data not shown).

Figure 2
figure 2

Analysis of immune gene expression in aposymbiotic larvae challenged with Sitophilus zeamais primary endosymbiont. Transcript of genes induced by an E. coli challenge were quantitated by qRT-PCR in untreated aposymbiotic larvae (control) and in larvae six hours after a 69 μl injection of either sterile water (sterile) or SZPE cells (heat-killed or viable) obtained from 50 bacteriomes dissected from symbiotic fourth-instar larvae. ESTs were classified according to sequence similarity as in Table 2. Each bar represents the mean of three independent measurements with standard error.

No significant difference was shown between water and symbiont injections for wPGRP1 (INF-9) whereas all the genes with sequence similarity to AMPs were significantly induced by the symbiont. Moreover, heat-killed E. coli was shown to induce a weaker immune response than the untreated E. coli (data not shown), possibly because bacterial growth and division in larvae may result in an increase of bacterial density and then in a stronger immune response [16]. Unlike E. coli, no difference was found between larvae injected with viable or heat-killed SZPE, which is consistent with an absence of symbiont proliferation in the hemocoel.

Gene expression analysis in the larval bacteriome

As the bacteriome fails to develop in the absence of endosymbionts, the immune gene expression profile was examined in the bacteriome and compared with the mean transcript levels from whole aposymbiotic larval tissues. As shown in Figure 3, most of the genes overexpressed in larvae infected with E. coli (Table 2) or SZPE (Figure 2) were slightly (or not at all) expressed in the bacteriome, with the exception of two genes: the wPGRP1 gene (INF-9; 39-fold) and one of the two coleoptericin genes (INF-18; 10-fold). Interestingly, this experiment also revealed a high expression level of the Tollip gene (INF-359; 5-fold) in the bacteriome tissue. Moreover, cytoskeleton-encoding genes were shown to be significantly underexpressed in the bacteriome. Actin transcripts were 800-fold less represented in the bacteriome when compared with whole aposymbiotic larvae. Although many cases of host cytoskeleton manipulation by bacteria have been described [28, 29], we cannot exclude the possibility that the difference in the levels of actin transcription between the bacteriome and the larvae is due to the high muscle actin content in the larvae. Similarly, the abundance (or absence) of certain gene products in major tissues (for example lysozyme abundance in gut and salivary gland) can create an apparent under (or over-) expression in the bacteriome, and for such genes these results should be considered with caution. However, it is noteworthy that comparisons between the bacteriocyte and the whole body expression profiles have provided previously meaningful results [15, 30].

Figure 3
figure 3

Analysis of gene expression in the bacteriome tissue. As described in Materials and Methods, transcripts of candidate genes were quantitated by qRT-PCR in whole aposymbiotic fourth-instar larvae (control) and in bacteriomes dissected from symbiotic fourth-instar larvae. ESTs were classified according to sequence similarity as in Table 2. Each bar represents the mean of three independent measurements with standard error. The asterisk represents a significant (p < 0.05) difference between the bacteriome and the control, and a significantly high expression in a bacteriome is indicated by an arrow. ND, non-determined.

Discussion

There is compelling evidence that intracellular bacterial infections can establish chronic non-septic relationships with host cells for many generations. Some bacteria (for example, Chlamydophila and Salmonella) invade and proliferate rapidly in their host, causing diseases, while others (for example, Mycobacterium tuberculosis) remain latent and may reactivate months or years after initial exposure to cause chronic disease [31]. In the insect world, some intracellular bacteria can generate a long-term relationship within the host without causing any disease or invading any host tissues except the germ cells and the bacteriocytes. The overexpression of the wPGRP1 gene in the weevil bacteriome suggests a close interaction between the host immune system and the endosymbionts.

To determine how the insect immune system recognizes the bacteria and prevents bacterial invasion of insect tissues, we first identified, by SSH technology, the immune genes induced in insects challenged with E. coli. Bioinformatic analysis scored a relatively high proportion of ESTs with similarities to immune-related genes, including several antibacterial peptides, lysozymes, PGRPs and proteins of the PPO pathway (Table 2), which highlights the efficiency of the subtractive approach. Further in silico analysis on full-length coding sequences has confirmed the identification of an i-type and a c-type lysozyme gene (INF-152, INF-282) and at least four antibacterial peptide genes: two members of the coleoptericin family (INF-18, INF-145), a defensin (INF-217) and a homolog of luxuriosin (INF-479), an antimicrobial peptide from Acalolepta luxuriosa [20]. We also validated a homolog of Tollip, a gene so far only studied in mammals that seems to be absent from the Drosophila genome, and a PGRP gene (wPGRP2, INF-441) additional to the one previously described (wPGRP1 gene, INF-9). The expression profile of all the putative immune genes was then characterized by larvae challenge experiments combined with RT-PCR transcript quantification (Table 2). Thus we have demonstrated that all the genes studied were inducible, except for homologs of Tollip (INF-359) and hemolymph proteinase 17 (INF-459), which appear to be constitutive.

To elucidate the particular bacteriome immune features, we first examined the systemic host response following aposymbiotic larvae injection with SZPE (Figure 2) and we then measured immune gene expression in the bacteriome of naturally infected larvae (Figure 3). Experimental infection of the larvae led to a systemic immune response with the up-regulation of all the genes with similarities to AMPs, providing evidence that SZPE is recognized as a microbial intruder while present in the hemolymph. It is noteworthy that in contrast to SZPE, Spiroplasma does not induce host transcripts encoding AMPs in Drosophila while present in hemolymph. However, the lack of bacterial cell wall structure of Spiroplasma is probably the best explanation for the apparent absence of host humoral response to this bacterium [32].

In the weevil symbiotic larvae, in addition to the wPGRP1 gene, one antibacterial peptide induced by bacterial challenge (that is, the coleoptericin gene) is overexpressed in the bacteriome. These data confirm the existence of an immune response in the bacteriome tissue and provide the first elements to explain how host-endosymbiont association can persist with regard to the host immune system and bacterial proliferation.

Since all AMP genes are up-regulated following a hemolymph challenge with SZPE, the unique overexpression of the coleoptericin gene in the bacteriome suggests a constitutive expression in this tissue rather than an induction by SZPE, unless coleoptericin gene regulation involves a separate pathway from those that trigger synthesis of other AMPs. We can only speculate about this question, as comparison between symbiont-full and symbiont-free bacteriocytes cannot be conducted since no symbiont-free bacteriome has been detected in Sitophilus aposymbiotic insects so far. Taking into consideration that some AMPs are expressed constitutively in cells potentially in contact with the environment [11, 33, 34], and that the coleoptericin gene encodes a signal peptide, these data suggest that the coleoptericin gene is expressed constitutively and secreted outside the bacteriocyte. Coleoptericin is currently under investigation to determine whether this peptide may prevent host tissue invasion by the endosymbiont, and/or may protect the endosymbiont population against infection with environmental bacteria, as the bacteriome tissue is intimately attached to the insect foregut. Moreover, another interesting issue would be to investigate whether coleptericin expression in the bacteriome can protect the host from pathogens. This would represent an additional example of 'symbiont-mediated protection' [35].

Although the absence of a host humoral immune response to an intracellular symbiont and a constitutive expression of coleoptericin in the bacteriome are consistent with current knowledge of the insect immune response, the overexpression of gene homologs to PGRP-LB (that is, wPGRP1) and Tollip is quite intriguing. PGRP-LB gene inhibition in Drosophila was indeed shown to result in a significant induction of antibacterial peptides [13], while Tollip has been implicated as a negative regulator of the mammalian immune response [24, 25]. Moreover, these two genes are overexpressed in the gut epithelia and have been proposed to play a role in the control of immune reactivity of the host to the presence of bacteria in the gut. Therefore, the overexpression of these two genes in the bacteriome may be consistent with an immune modulation that inhibits the production of AMPs in this organ, as it was described in the gut epithelia. Together with the overexpression of an AMP (that is, coleoptericin), these results uncover striking similarities between the bacteriome immune profile and the local immune response in gut epithelium that is in permanent contact with commensal and mutualistic bacteria [9, 11, 13, 33, 36, 37]. A gene silencing study in weevil is currently under development to validate coleoptericin, wPGRP1 and Tollip gene function and their contribution to the maintenance of symbiosis.

It is noteworthy that this work was performed on a recently established symbiosis where, contrary to ancient symbiotic associations known for a drastically reduced endosymbiont genome size such as in the aphid/Buchnera association [38], SZPE exhibits similar features to those of free-living Gram-negative bacteria [3943]. Hence, whether and how the host immune response evolves in parallel to symbiont genome reduction constitutes a key aspect in the understanding of host-symbiont interaction in the course of evolution. Recently, a transcriptomic study was performed on the aphid bacteriocyte to identify a gene of interest in intracellular symbiosis [30]. Intriguingly, no homolog to any known immune genes was shown to be expressed in the aphid bacteriocyte except for i-type lysozymes, whose function in bacteriocytes remains unknown. Since both the weevil i-type and c-type lysozymes were shown to be weakly expressed in the bacteriocytes, these data suggest that the bacteriocyte immune profile may have evolved with bacterial coevolution with the host, unless the immune response of weevils (holometabolous insects) has diverged significantly from that of aphids (heterometabolous insects).

Whatever the antimicrobial effectors (that is, antibacterial peptides, lysozyme), their permanent expression in bacteriocytes could represent a common mechanism to restrict the localization of mutualistic symbionts, which is probably necessary to optimize the host-symbiont interaction. Without this confinement, intracellular symbiotic bacteria may invade the whole organism without inducing any antibacterial peptide synthesis, as has been shown in Drosophila infected by Wolbachia, a parasitic symbiont widely distributed in host tissues [44].

Conclusion

This work provides the first immune gene expression profile in the insect bacteriome and reveals the overexpression of at least three homologs to immune genes: a member of the PGRP gene family involved in bacteria recognition, an antibacterial peptide involved in bacterial clearance and a gene involved in immunomodulation in mammals. This immune profile uncovers some striking similarities between the bacteriome and the gut epithelium, which is in constant contact with environmental bacteria. Taking into account that the endosymbiont is recognized as an intruder in the host hemolymph, these findings also indicate that the host immune system may prevent endosymbiont invasion into insect tissues in a manner similar to the gut local immune response that helps to confine the microbiota to the gut, avoiding a permanent systemic response to the commensal bacteria.

Methods

Insect rearing

Insects from a SZPE-monosymbiotic strain (S. zeamais Lagoa) and from the corresponding aposymbiotic strain were reared and collected as described in Anselme et al. [16].

Bacterial challenge

Fourth-instar aposymbiotic larvae were challenged by pricking with sterile sharpened needles (mock infection) or with needles previously dipped into a pellet from E. coli (TOP10, Invitrogen) overnight cultures, and kept in a moist atmosphere at 27.5°C for 3, 6, and 12 hours. Living larvae were stored at -80°C for RNA preparation. Unchallenged larvae (naive larvae) were treated in parallel as controls.

For the study of the host immune response to SZPE, larvae were injected with 69 nl of either sterile water or bacterial solution containing approximately 1 × 105 viable or heat-killed (5 min at 95°C) bacterial cells. SZPE solution was freshly prepared from bacteriomes of fourth-instar larvae. Fifty bacteriomes were dissected in buffer A (25 mM KCl, 10 mM MgCl2, 250 mM sucrose, 35 mM Tris-HCl, pH 7.5), transferred in a Dounce Teflon Homogenizer and gently crushed in buffer A. After removal of cellular debris by low-speed centrifugation (400 g, 10 min), bacteria were pelleted (10000 g, 5 min) and resuspended to a concentration of approximately 1.45 × 106 bacteria/μl in buffer A.

Total RNA isolation

Total RNA from infected and naive larvae was extracted using TRIzol® Reagent (Invitrogen), treated with RNase-Free DNase (Promega) and purified through RNeasy mini kit columns (Qiagen) as described in the manufacturer's procedures. After purification, the RNA concentration of each sample was measured by the Nanodrop® spectrophotometer and total RNA quality was checked by electrophoresis.

Suppression subtractive hybridization

We applied SSH technology by using a PCR-selected cDNA subtraction kit (Clontech laboratories). For the synthesis of cDNA from E. coli-infected larvae used as a tester, three RNA samples were extracted 3, 6 and 12 hours after infection and mixed prior to cDNA synthesis.

Briefly, cDNA synthesis was carried out on a 1 μg aliquot from a mix of 1 μg of each infected RNA sample for the tester and from naive larvae RNA for the driver (SMART PCR cDNA synthesis kit, Clontech laboratories). After phenol-chloroform purification, the infected larval cDNA was digested with RsaI, ligated to two adaptors separately (1 and 2R), then hybridized to RsaI-digested naive larvae cDNA. After hybridization, the subtracted cDNAs were amplified by PCR according to the manufacturer's instructions (Clontech, cDNA subtraction kit). The amplified products were directly cloned into a pCR 2.1-TOPO plasmid (Invitrogen).

Clone sequencing, Blast homology searching and GO assignment

After transformation by electroporation, around 500 colonies were recovered from Luria Broth agar plates. They were subjected to plasmid extraction (NucleoSpin® Plasmid Kit, Macherey-Nagel) and sequenced in the Institut Cavanilles de Biodiversitat i Biologia Evolutiva (València, Spain). The 485 available sequences were trimmed using SeqClean http://compbio.dfci.harvard.edu/tgi/software/ to remove flanking vector sequences, adaptors (match with at least 98% identity over at least 11 base pairs) and poly(A/T) tails, and remaining sequences shorter than 60 base pairs were discarded. Because subtracted library construction includes a step of RsaI digestion and adaptor ligation where chimeric cDNA could be formed by cDNA ligation, we checked for RsaI sites in the ESTs and discarded the chimeras from our dataset.

Remaining sequences were clustered using the TGICL assembly program [45] and consensus sequences were conserved for the subtracted library. All sequences (consensus and singletons) were then compared against UniProt using BlastX. For GO assignment, we retained the first hit with at least one GO annotation among the top five hits. GO annotation results were then classified using WEGO [46].

Sequencing of immune-relevant full-length cDNA

The complete sequences of the transcripts of interest were obtained by 3' and 5' RACE, performed with the SMART RACE cDNA Amplification Kit including the Advantage II PCR kit (Clontech Laboratories). For each gene, the nucleotide sequences of the 3' and/or 5' primers (GSP1 and GSP2) were designed on the corresponding EST (see Additional file 2). The first-strand cDNA used for 5' and 3'RACE were produced by using 1 μg of the infected RNA mix prepared for the SSH method, and using the primers provided in the kit. Amplification of the RACE products was carried out according to the manufacturer's instructions. PCR fragments were gel-purified with the Nucleotrap Gel Extraction Kit (Clontech Laboratories) and inserted into the plasmid vector pCR2.1-TOPO (Invitrogen). The sequences were generated by Genome Express Company (Grenoble, France) with the M13 and the M13rev vector primers.

The full-length sequence of the transcript was predicted using the MacMolly software package according to the RACE PCR product and according to the presence of both in-frame stop codons upstream of the prospective methionine start site and of a poly(A) tail following the prospective stop codon. The predicted sequence was then confirmed by PCR amplification and sequencing. The putative proteins were compared with protein sequences in Swiss-Prot-Trembl to confirm the first homology, and analyzed using the InterProScan http://www.ebi.ac.uk/InterProScan that combines different protein signature recognition methods.

Real-Time RT-PCR transcript quantification

Real-time RT-PCR transcript quantification was performed with a LightCycler® instrument using the LightCycler Fast Start DNA Master SYBR Green I kit (Roche Diagnostics), as described in Anselme et al. [16]. RNA extractions on three independent biological samples were carried out for each condition (naive larvae, mock-infected larvae and infected larvae) 6 hours after treatment and on samples of 100 bacteriomes dissected from naive symbiotic larvae. After sample purification, reverse transcription into the first strand cDNA was made with the First-Strand Synthesis System for RT-PCR kit (Invitrogen) using oligo(dT) primers. PCR primers were designed on ESTs or cDNA sequences, when necessary, to amplify fragments of less than 300 bp (see Additional file 2).

The PCR reactions were carried out in LightCycler capillaries in a final volume of 20 μl containing 2 μl of cDNA samples (diluted 5-fold), 3.5 mM MgCl2, 0.5 μM of each primer and 2 μl of LC-Fast Start Reaction Master SYBR Green I. After 8 min at 95°C, the cycling conditions were as follows: 45 cycles at 95°C for 10 s, primer annealing temperature for 20 s and then 72°C for 30 s. For product identification, a melting curve was constructed at the end of each PCR by heating for 5 s at 95°C and for 15 s at a temperature corresponding to 10°C more than the primer annealing temperature, and then increasing the temperature up to 95°C with increment rates of 0.1°C/s. Reactions were stopped by cooling at 4°C.

For the individual samples, the crossing point (Cp) and, according to the standard curve, the concentrations of the transcripts were determined. The melting curves of each sample were analyzed and the concentrations of samples presenting primer-dimer formation were considered as 'non-determined' (ND). Data were normalized using the ratio of the target cDNA concentration to that of the house keeping gene, glyceraldehyde 3-phosphate dehydrogenase (gapdh). The expression of this gene is not significantly influenced by the treatments and it is similar to the expression of the ribosomal protein L29 gene (data not shown). Normalized data were analyzed using nonparametric tests.

The sequences of the EST and the full-length cDNA reported in this paper have been deposited in the GenBank database (Accession nos. from EY122775 to EY123248 and from EU282111 to EU282122, respectively). A more detailed analysis on the ESTs (for example, Accession nos., ESTs assembly, Blast results) is available at http://mandragore.univ-lyon1.fr/sitozea (login: sitophilus, password: zeamais).

Authors' note

During the preparation of this manuscript, the experimental identification of genes that are induced in response to septic injury in the red flour beetle Tribolium castaneum using the SSH method has been published [47]. By comparing EST sequence to Tribolium sequence databases, the authors have identified 75 immune-inducible genes in T. castaneum potentially involved in immune defense, signaling, and other immunity-linked cellular processes including homologs of, for example, Toll, PGRP-SC, lysozyme, and multiple isoforms of defensins and thaumatin-like peptides. They have also performed a qRT-PCR to determine transcriptional regulation of selected genes (that is, defensin, thaumatin and stress-related genes) in response to either septic versus sterile injury or heat shock.

While this study explored the Tribolium immune response considering all cellular processes, we focused our study on Sitophilus immune pathways and effectors. It should be noted that, in contrast to this study, we have identified members of different antibacterial peptide families (for example, coleoptericin, defensin). We also identified some homologs of thaumatin-like peptides (for example, INF-475, INF-332, INF-CL57Contig1), but were unable to quantify them without specific amplification.

References

  1. Buchner P: Endosymbiosis of Animals with Plant Microorganisms. 1965, New York, London, Sydney: Interscience Publishers, a division of Wiley & sons, Inc

    Google Scholar 

  2. Heddi A, Grenier AM, Khatchadourian C, Charles H, Nardon P: Four intracellular genomes direct weevil biology: Nuclear, mitochondrial, principal endosymbionts, and Wolbachia. Proc Natl Acad Sci U S A. 1999, 96: 6814-6819. 10.1073/pnas.96.12.6814.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Braendle C, Miura T, Bickel R, Shingleton AW, Kambhampati S, Stern DL: Developmental origin and evolution of bacteriocytes in the aphid-Buchnera symbiosis. PLoS Biol. 2003, 1 (1): e21-10.1371/journal.pbio.0000021.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Douglas AE: Host benefit and the evolution of specialization in symbiosis. Heredity. 1998, 81: 599-603. 10.1038/sj.hdy.6884550.

    Article  Google Scholar 

  5. Rakoff-Nahoum S, Paglino J, Eslami-Varzaneh F, Edberg S, Medzhitov R: Recognition of commensal microflora by toll-like receptors is required for intestinal homeostasis. Cell. 2004, 118: 229-241. 10.1016/j.cell.2004.07.002.

    Article  CAS  PubMed  Google Scholar 

  6. Goodson MS, Kojadinovic M, Troll JV, Scheetz TE, Casavant TL, Soares MB, McFall-Ngai MJ: Identifying components of the NF-kappaB pathway in the beneficial Euprymna scolopes-Vibrio fischeri light organ symbiosis. Appl Environ Microbiol. 2005, 71 (11): 6934-6946. 10.1128/AEM.71.11.6934-6946.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Mazmanian SK, Liu CH, Tzianabos AO, Kasper DL: An immunomodulatory molecule of symbiotic bacteria directs maturation of host immune system. Cell. 2005, 122: 107-118. 10.1016/j.cell.2005.05.007.

    Article  CAS  PubMed  Google Scholar 

  8. Cash HL, Whitham CV, Behrendt CL, Hooper LV: Symbiotic bacteria direct expression of an intestinal bactericidal lectin. Science. 2006, 313: 1126-1030. 10.1126/science.1127119.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Lemaitre B, Hoffmann J: The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007, 25: 697-743. 10.1146/annurev.immunol.25.022106.141615.

    Article  CAS  PubMed  Google Scholar 

  10. Christensen BM, Li J, Chen CC, Nappi AJ: Melanization immune responses in mosquito vectors. Trends Parasitol. 2005, 21: 192-199. 10.1016/j.pt.2005.02.007.

    Article  CAS  PubMed  Google Scholar 

  11. Tzou P, Ohresser S, Ferrandon D, Capovilla M, Reichhart JM, Lemaitre B, Hoffmann JA, Imler JL: Tissue-specific inducible expression of antimicrobial peptide genes in Drosophila surface epithelia. Immunity. 2000, 13: 737-748. 10.1016/S1074-7613(00)00072-8.

    Article  CAS  PubMed  Google Scholar 

  12. Bischoff V, Vignal C, Duvic B, Boneca IG, Hoffmann JA, Royet J: Downregulation of the Drosophila immune response by peptidoglycan-recognition proteins SC1 and SC2. PLoS Pathog. 2006, 2 ( 2): e14-10.1371/journal.ppat.0020014.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Zaidman-Remy A, Herve M, Poidevin M, Pili-Floury S, Kim MS, Blanot D, Oh BH, Ueda R, Mengin-Lecreulx D, Lemaitre B: The Drosophila Amidase PGRP-LB Modulates the Immune Response to Bacterial Infection. Immunity. 2006, 24: 463-473. 10.1016/j.immuni.2006.02.012.

    Article  CAS  PubMed  Google Scholar 

  14. Heddi A, Nardon P: Sitophilus oryzae L.: A model for intracellular symbiosis in the Dryophthoridae weevils (Coleoptera). Symbiosis. 2005, 39: 1-11.

    Google Scholar 

  15. Heddi A, Vallier A, Anselme C, Xin H, Rahbé Y, Wäckers F: Molecular and cellular profiles of insect bacteriocytes: mutualism and harm at the initial evolutionary step of symbiogenesis. Cell Microbiol. 2005, 7: 293-305.

    Article  CAS  PubMed  Google Scholar 

  16. Anselme C, Vallier A, Balmand S, Fauvarque MO, Heddi A: Host PGRP gene expression and bacterial release in endosymbiosis of the weevil Sitophilus zeamais. Appl Environ Microbiol. 2006, 72: 6766-6772. 10.1128/AEM.00942-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Hosaka M, Nagahama M, Kim WS, Watanabe T, Hatsuzawa K, Ikemizu J, Murakami K, Nakayama K: Arg-X-Lys/Arg-Arg motif as a signal for precursor cleavage catalyzed by furin within the constitutive secretory pathway. J Biol Chem. 1991, 266: 12127-12130.

    CAS  PubMed  Google Scholar 

  18. Lazzaro BP, Clark AG: Molecular population genetics of inducible antibacterial peptide genes in Drosophila melanogaster. Mol Biol Evol. 2003, 20: 914-923. 10.1093/molbev/msg109.

    Article  CAS  PubMed  Google Scholar 

  19. Bulet P, Stocklin R: Insect antimicrobial peptides: structures, properties and gene regulation. Protein Pept Lett. 2005, 12: 3-11. 10.2174/0929866053406011.

    Article  CAS  PubMed  Google Scholar 

  20. Ueda K, Saito A, Imamura M, Miura N, Atsumi S, Tabunoki H, Watanabe A, Kitami M, Sato R: Purification and cDNA cloning of luxuriosin, a novel antibacterial peptide with Kunitz domain from the longicorn beetle, Acalolepta luxuriosa. Biochim Biophys Acta. 2005, 1722: 36-42.

    Article  CAS  PubMed  Google Scholar 

  21. Jolles P, Jolles J: What's new in lysozyme research? Always a model system, today as yesterday. Mol Cell Biochem. 1984, 63: 165-189. 10.1007/BF00285225.

    Article  CAS  PubMed  Google Scholar 

  22. Zavalova LL, Baskova IP, Lukyanov SA, Sass AV, Snezhkov EV, Akopov SB, Artamonova II, Archipova VS, Nesmeyanov VA, Kozlov DG, et al: Destabilase from the medicinal leech is a representative of a novel family of lysozymes. Biochim Biophys Acta. 2000, 1478: 69-77.

    Article  CAS  PubMed  Google Scholar 

  23. Kim MS, Byun M, Oh BH: Crystal structure of peptidoglycan recognition protein LB from Drosophila melanogaster. Nat Immunol. 2003, 4: 787-793. 10.1038/ni952.

    Article  CAS  PubMed  Google Scholar 

  24. Burns K, Clatworthy J, Martin L, Martinon F, Plumpton C, Maschera B, Lewis A, Ray K, Tschopp J, Volpe F: Tollip, a new component of the IL-1RI pathway, links IRAK to the IL-1 receptor. Nat Cell Biol. 2000, 2: 346-351. 10.1038/35014038.

    Article  CAS  PubMed  Google Scholar 

  25. Zhang G, Ghosh S: Negative regulation of Toll-like receptor-mediated signaling by Tollip. J Biol Chem. 2002, 277: 7059-7065. 10.1074/jbc.M109537200.

    Article  CAS  PubMed  Google Scholar 

  26. Didierlaurent A, Brissoni B, Velin D, Aebi N, Tardivel A, Kaslin E, Sirard JC, Angelov G, Tschopp J, Burns K: Tollip regulates proinflammatory responses to interleukin-1 and lipopolysaccharide. Mol Cell Biol. 2006, 26 (3): 735-742. 10.1128/MCB.26.3.735-742.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Ramet M, Manfruelli P, Pearson A, Mathey-Prevot B, Ezekowitz RA: Functional genomic analysis of phagocytosis and identification of a Drosophila receptor for E. coli. Nature. 2002, 416: 644-648. 10.1038/nature735.

    Article  CAS  PubMed  Google Scholar 

  28. Barbieri JT, Riese MJ, Aktories K: Bacterial toxins that modify the actin cytoskeleton. Annu Rev Cell Dev Biol. 2002, 18: 315-344. 10.1146/annurev.cellbio.18.012502.134748.

    Article  CAS  PubMed  Google Scholar 

  29. Kimbell JR, McFall-Ngai MJ: Symbiont-induced changes in host actin during the onset of a beneficial animal-bacterial association. Appl Environ Microbiol. 2004, 70 (3): 1434-1441. 10.1128/AEM.70.3.1434-1441.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Nakabachi A, Shigenobu S, Sakazume N, Shiraki T, Hayashizaki Y, Carninci P, Ishikawa H, Kudo T, Fukatsu T: Transcriptome analysis of the aphid bacteriocyte, the symbiotic host cell that harbors an endocellular mutualistic bacterium, Buchnera. Proc Natl Acad Sci U S A. 2005, 102 (15): 5477-5482. 10.1073/pnas.0409034102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Scanga CA, Mohan VP, Yu K, Joseph H, Tanaka K, Chan J, Flynn JL: Depletion of CD4(+) T cells causes reactivation of murine persistent tuberculosis despite continued expression of interferon gamma and nitric oxide synthase 2. J Exp Med. 2000, 192 (3): 347-358. 10.1084/jem.192.3.347.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Hurst GD, Anbutsu H, Kutsukake M, Fukatsu T: Hidden from the host: Spiroplasma bacteria infecting Drosophila do not cause an immune response, but are suppressed by ectopic immune activation. Insect Mol Biol. 2003, 12: 93-97. 10.1046/j.1365-2583.2003.00380.x.

    Article  CAS  PubMed  Google Scholar 

  33. Schröder J-M: Epithelial peptide antibiotics. Biochem Pharmacol. 1999, 57: 121-134. 10.1016/S0006-2952(98)00226-3.

    Article  PubMed  Google Scholar 

  34. Ryu JH, Nam KB, Oh CT, Nam HJ, Kim SH, Yoon JH, Seong JK, Yoo MA, Jang IH, Brey PT, et al: The homeobox gene Caudal regulates constitutive local expression of antimicrobial peptide genes in Drosophila epithelia. Mol Cell Biol. 2004, 24 (1): 172-185. 10.1128/MCB.24.1.172-185.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Haines LR: Symbiont-mediated protection. Proc Biol Sci. 2008, 275 (1633): 353-361. 10.1098/rspb.2007.1211.

    Article  Google Scholar 

  36. Werner T, Liu G, Kang D, Ekengren S, Steiner H, Hultmark D: A family of peptidoglycan recognition proteins in the fruit fly Drosophila melanogaster. Proc Natl Acad Sci U S A. 2000, 97: 13772-13777. 10.1073/pnas.97.25.13772.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Melmed G, Thomas LS, Lee N, Tesfay SY, Lukasek K, Michelsen KS, Zhou Y, Hu B, Arditi M, Abreu MT: Human intestinal epithelial cells are broadly unresponsive to Toll-like receptor 2-dependent bacterial ligands: implications for host-microbial interactions in the gut. J Immunol. 2003, 170: 1406-1415.

    Article  CAS  PubMed  Google Scholar 

  38. Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H: Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature. 2000, 407: 81-86. 10.1038/35024074.

    Article  CAS  PubMed  Google Scholar 

  39. Charles H, Condemine G, Nardon C, Nardon P: Genome size characterization of the principal endocellular symbiotic bacteria of the weevil Sitophilus oryzae, using pulsed field gel electrophoresis. Insect Biochem Mol Biol. 1997, 27: 345-350. 10.1016/S0965-1748(97)00005-2.

    Article  CAS  Google Scholar 

  40. Heddi A, Charles H, Khatchadourian C, Bonnot G, Nardon P: Molecular characterization of the principal symbiotic bacteria of the weevil Sitophilus oryzae: A peculiar G – C content of an endocytobiotic DNA. J Mol Evol. 1998, 47: 52-61. 10.1007/PL00006362.

    Article  CAS  PubMed  Google Scholar 

  41. Charles H, Heddi A, Rahbé Y: A putative insect intracellular endosymbiont stem clade, within the Enterobacteriaceae, inferred from phylogenetic analysis based on a heterogeneous model of DNA evolution. C R Acad Sci III. 2001, 324: 489-494.

    Article  CAS  PubMed  Google Scholar 

  42. Dale C, Plague GR, Wang B, Ochman H, Moran NA: Type III secretion systems and the evolution of mutualistic endosymbiosis. Proc Natl Acad Sci U S A. 2002, 99 (19): 12397-12402. 10.1073/pnas.182213299.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Lefèvre C, Charles H, Vallier A, Delobel B, Farrell B, Heddi A: Endosymbiont phylogenesis in the Dryophthoridae weevils: Evidence for bacterial replacement. Mol Biol Evol. 2004, 21: 965-973. 10.1093/molbev/msh063.

    Article  PubMed  Google Scholar 

  44. Bourtzis K, Pettigrew MM, O'Neill SL: Wolbachia neither induces nor suppresses transcripts encoding antimicrobial peptides. Insect Mol Biol. 2000, 9: 635-639. 10.1046/j.1365-2583.2000.00224.x.

    Article  CAS  PubMed  Google Scholar 

  45. Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.

    Article  CAS  PubMed  Google Scholar 

  46. Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (Web Server issue): W293-297. 10.1093/nar/gkl031.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Altincicek B, Knorr E, Vilcinskas A: Beetle immunity: Identification of immune-inducible genes from the model insect Tribolium castaneum. Dev Comp Immunol. 2008, 32: 585-595. 10.1016/j.dci.2007.09.005.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We gratefully acknowledge M. Veyre and S. Guyon for their technical assistance, Y. Rahbé for helpful discussions and V. James for English corrections. The authors would also like to thank the anonymous reviewers for their constructive criticisms. This work was supported by the Institut National de la Recherche Agronomique (INRA), the Institut National des Sciences Appliquées (INSA) and The French ANR-06-BLAN-0316 (EndoSymArt).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Abdelaziz Heddi.

Additional information

Authors' contributions

CA designed and performed experiments, analyzed data (statistics and bioinformatics), wrote the paper and participated in bioinformatic analysis of the subtracted library; VPB participated in molecular studies and sequenced the subtracted library. AV carried out dissections and real-time RT-PCR, and CVM participated in molecular studies and sequence analysis. DC carried out bioinformatics analysis and the URL construction for the subtracted library; AL and AM sequenced the subtracted library and provided critical comments on the manuscript. AH conceived the study, coordinated the work and helped to draft and write the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Anselme, C., Pérez-Brocal, V., Vallier, A. et al. Identification of the Weevil immune genes and their expression in the bacteriome tissue . BMC Biol 6, 43 (2008). https://0-doi-org.brum.beds.ac.uk/10.1186/1741-7007-6-43

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1741-7007-6-43

Keywords