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

Comparative in vivo gene expression of the closely related bacteria Photorhabdus temperata and Xenorhabdus koppenhoeferi upon infection of the same insect host, Rhizotrogus majalis

Abstract

Background

Photorhabdus and Xenorhabdus are Gram-negative, phylogenetically related, enterobacteria, forming mutualism with the entomopathogenic nematodes Heterorhabditis and Steinernema, respectively. The mutualistic bacteria living in the intestines of the nematode infective juveniles are pathogenic to the insect upon release by the nematodes into the insect hemolymph. Such a switch needs activation of genes that promote bacterial virulence. We studied in vivo gene expression in Photorhabdus temperata and Xenorhabdus koppenhoeferi upon infection of the white grub Rhizotrogus majalis using selective capture of transcribed sequences technique.

Results

A total of 40 genes in P. temperata and 39 in X. koppenhoeferi were found to be upregulated in R. majalis hemolymph at 24 h post infection. Genomic presence or upregulation of these genes specific in either one of the bacterium was confirmed by the assay of comparative hybridization, and the changes of randomly selected genes were further validated by quantitative real-time PCR. The identified genes could be broadly divided into seven functional groups including cell surface structure, regulation, virulence and secretion, stress response, intracellular metabolism, nutrient scavenging, and unknown. The two bacteria shared more genes in stress response category than any other functional group. More than 60% of the identified genes were uniquely induced in either bacterium suggesting vastly different molecular mechanisms of pathogenicity to the same insect host. In P. temperata lysR gene encoding transcriptional activator was induced, while genes yijC and rseA encoding transcriptional repressors were induced in X. koppenhoeferi. Lipopolysaccharide synthesis gene lpsE was induced in X. koppenhoeferi but not in P. temperata. Except tcaC and hemolysin related genes, other virulence genes were different between the two bacteria. Genes involved in TCA cycle were induced in P. temperata whereas those involved in glyoxylate pathway were induced in X. koppenhoeferi, suggesting differences in metabolism between the two bacteria in the same insect host. Upregulation of genes encoding different types of nutrient uptake systems further emphasized the differences in nutritional requirements of the two bacteria in the same insect host. Photorhabdus temperata displayed upregulation of genes encoding siderophore-dependent iron uptake system, but X. koppenhoeferi upregulated genes encoding siderophore-independent ion uptake system. Photorhabdus temperata induced genes for amino acid acquisition but X. koppenhoeferi upregulated malF gene, encoding a maltose uptake system. Further analyses identified possible mechanistic associations between the identified gene products in metabolic pathways, providing an interactive model of pathogenesis for each bacterium species.

Conclusion

This study identifies set of genes induced in P. temperata and X. koppenhoeferi upon infection of R. majalis, and highlights differences in molecular features used by these two closely related bacteria to promote their pathogenicity in the same insect host.

Background

Entomopathogenic Gram-negative enterobacteria Photorhabdus and Xenorhabdus form symbioses with the entomopathogenic nematodes Heterorhabditis and Steinernema, respectively [1]. These bacteria not only have similar biology but are also phylogenetically related based on 16s rDNA sequence identities [2]. They naturally colonize intestines of the nematode infective juveniles which invade susceptible insects to release the bacteria into the hemolymph. The bacteria multiply in the hemolymph, killing the insect host within 24-48 h and converting the cadaver into a food source suitable for nematode growth and reproduction. After 1-3 rounds of nematode reproduction, the bacteria recolonize the emerging infective juveniles ensuring their transmission to a new host [1].

Available evidence suggests that Photorhabdus and Xenorhabdus encode specific factors to engage in a pathogenic relationship with the insect host [3]. The published genome sequence of Photorhabdus luminescens TT01 strain indicates that virulence genes are encoded within a number of pathogenicity islands located on the bacterial chromosome [4, 5]. Besides producing toxins to cause insect death, Photorhabdus and Xenorhabdus have to first evade the insect's immune response to establish a successful infection. The two bacteria differ in mechanisms by which they evade host immune responses. For example, in Photorhabdus, mutational inactivation of phoP gene results in increased sensitivity to insect immune response and decreased virulence towards insects [6, 7], while in X. nematophila, phoPQ mutants are more susceptible to immune response but are fully virulent [3]. P. luminescens produces a signaling molecule AI-2 to resist reactive oxygen species [8] and phenylpropanoid chemical ST to inhibit the activity of antimicrobial enzyme PO and formation of melanotic nodules [9], but the strategy used by X. nematophila appears to be that of suppression of transcripts involved in the insect immune response [10–12]. In addition, P. luminescens encodes a type III secretion system and one of the effectors, LopT, suppresses phagocytosis and reduces nodulation by haemocytes [13, 14]. However, the genomes of Xenorhabdus bovienii and X. nematophila do not show homologues of LopT or a dedicated type III secretion system [3], and Xenorhabdus likely uses flagellar regulatory cascade to control cellular immunity [15].

Despite the identification of the above regulatory and virulence factors from Photorhabdus and Xenorhabdus, the full profiles of genes expressed during insect infection by these two closely related bacteria are unknown. Münch et al. [16] identified 29 promoters induced by P. luminescens in the presence of isolated Galleria mellonella homogenate using a differential fluorescence induction approach. However, treatment with insect homogenate might neglect physicochemical parameters as inducers which could result in different sets of genes being upregulated during natural infection of live insects. We used selective capture of transcribed sequences (SCOTS) technique to study in vivo gene expression in Photorhabdus temperata (associated with the nematode Heterorhabditis bacteriophora GPS11) and Xenorhabdus koppenhoeferi (associated with the nematode Steinernema scarabaei AMK001) during infection of the same insect host Rhizotrogus majalis. They were selected according to their distinct pathogenicity in insects of the family of Scarabaeidae - S. scarabaei is more virulent to R. majalis than H. bacteriophora. However, unlike their nematode partners, both bacterial species are highly virulent to Rhizotrogus majalis. In addition, although X. koppenhoeferi grows slower than P. temperata both in vitro and in vivo, the LD50 and LT50 value do not differ between these two bacteria [17]. We hypothesize that the two bacteria species will use both common and distinct molecular mechanisms during infection of the same insect host.

Results

Isolation of bacterial transcripts induced in Rhizotrogus majalis

To isolate bacterial genes induced upon infection, total RNA was isolated from live grubs at 24 h post injection with P. temperata or X. koppenhoeferi, as well as from 48 h log-phase bacterial cultures grown in Brain Heart Infusion (BHI) broth that were used to inject R. majalis. Random-primed bacterial cDNAs were normalized by hybridization to biotinylated bacterial genomic DNA that had been blocked beforehand using bacterial ribosomal RNA operon, resulting in sampling of bacterial mRNA transcripts apart from its ribosomal and insect transcripts. The normalized bacterial cDNAs representing total mRNA transcripts produced by bacteria grown in vitro or within the infected insects were named as normalized in vitro and in vivo cDNA libraries, respectively. To isolate bacterial mRNA transcripts preferentially induced during infection of the insect compared to the culture, normalized in vivo cDNAs were enriched by hybridization to biotinylated bacterial genomic DNA that had been pre-hybridized with rRNA operon and in vitro normalized cDNAs. The enriched cDNAs representing in vivo-induced genes, which were either lower in abundance or absent in 48 h in vitro bacterial cultures were then cloned into a TA cloning vector to construct in vivo enriched cDNA libraries. A total of 384 clones (192 for each bacterium) derived from enriched cDNA libraries were randomly picked for screening. Out of these, the clones (150 for P. temperata and 140 for X. koppenhoeferi) that showed stronger signal upon hybridization with in vivo-derived cDNAs compared to in vitro-derived cDNAs (Fig. 1) were sequenced and analyzed using the non-redundant algorithms of BLAST in NCBI and database of Xenorhabdus Genomes. These sequenced clones were therefore considered to represent reliable transcripts specifically expressed in P. temperata or X. koppenhoeferi after 24 h infection of R. majalis.

Figure 1
figure 1

Schematic presentation of the Selective Capture of Transcribed Sequences (SCOTS) technique followed by Southern blot analysis of SCOTS identified sequences. In panel A, normalized bacterial cDNAs were obtained directly from bacteria grown in vitro in the Brain Heart Infusion broth or in vivo in the infected insect Rhizotrogus majalis. In panel B, cDNAs corresponding to genes preferentially expressed in R. majalis relative to the broth were enriched by differential cDNA hybridization. The enriched cDNAs were transformed into a cloning vector to build the cDNA library. Cloned inserts were amplified by PCR, equally transferred to two nylon membranes, and probed with digoxigenin labeled normalized in vivo cDNAs (left) or normalized in vitro cDNAs (right) as described in Methods. The dots at the same position in the two arrays were loaded with the same amplicon of each individual clone from the enriched cDNA library, and the concentration of probes was standardized to be same.

Analysis of bacterial genes induced in Rhizotrogus majalis

Most identified transcripts were either similar to sequences from P. luminescens or X. nematophila, suggesting that bacterial transcripts were successfully enriched from the infected insects. For both bacteria, each isolated transcript was identified at least twice from the sequenced clones. The coverage of in vivo enriched libraries was then evaluated using Analytic Rarefaction according to the redundancy of each identical transcript in the sequenced clones. As indicated by the rarefaction curves (Fig. 2), saturation was achieved for the sequenced clones, indicating that most of representative genes in the in vivo enriched cDNA libraries were identified. Totally, we isolated 40 different sequences from P. temperata (Table 1) and 39 from X. koppenhoeferi (Table 2). Based on homology searches, the genes corresponding to distinct sequences could be divided into seven functional groups: cell surface structure (Class I), intracellular metabolism (Class II), nutrient scavenging (Class III), regulation (Class IV), stress response (Class V), virulence (Class VI) and uncharacterized genes (Class VII). Genes involved in metabolism constitute about 25% to 33% of the total number of in vivo-induced genes (Fig. 3), showing large numbers of metabolic changes during insect infection. Semi-quantitative reverse transcriptase PCR performed on 10 randomly selected identified genes (five from P. temperata: cysK, fliA, lysR, pchC and pmt1 and five from X. koppenhoeferi: aceE, dacC, cobJ, malF, and yijC) showed that these genes were detectable from the in vivo-derived cDNA but not from in vitro-derived cDNA. In addition, 18S rRNA gene, a eukaryotic housekeeping gene, was detected in the in vivo cDNA population before normalization, but not after normalization, suggesting that bacterial cDNAs were successfully purified from the insect cDNAs after SCOTS normalization. Also gyrase A gene, a prokaryotic housekeeping gene, was not detected from the enriched in vivo cDNAs, but was detected in all other cDNA populations including in vitro and in vivo cDNAs before and after normalization, indicating that only differentially expressed genes were captured after enrichment.

Figure 2
figure 2

Rarefaction analysis curves demonstrating coverage of cDNA libraries for genes upregulated in bacteria Photorhabdus temperata or Xenorhabdus koppenhoeferi during infection of the insect Rhizotrogus majalis.

Figure 3
figure 3

Distribution of SCOTS isolated in vivo -induced Photorhabdus temperata and Xenorhabdus koppenhoeferi genes among different functional classes: cell structure (I), intracellular metabolism (II), nutrient scavenging (III), regulation (IV), adaptation to environmental stresses (V), virulence and secretion (VI), and unknown function. The percentage of genes for P. temperata or X. koppenhoeferi was presented as ratio of the number of genes in each functional class to the total number of genes from the respective bacterium, and the percentage of genes for common to both was calculated as ratio of the number of common genes in each functional class to the sum of genes from both bacteria.

Table 1 SCOTS identified Photorhabdus temperata genes induced upon infection of the white grub Rhizotrogus majalis.
Table 2 SCOTS identified Xenorhabdus koppenhoeferi genes induced upon infection of the white grub Rhizotrogus majalis.

Identification of Photorhabdus temperata or Xenorhabdus koppenhoeferi specific genes

To identify genes that were distinctive in P. temperata or X. koppenhoeferi genome or exclusively induced in one or the other bacterium, comparative hybridization was performed based upon the identified genes. Twenty nine of the 40 distinct P. temperata transcripts have similar sequences in X. koppenhoeferi genome, thus 11 were P. temperata specific; 30 of the 39 distinct X. koppenhoeferi transcripts share sequence similarity with P. temperata genome, thus 9 were specific to X. koppenhoeferi (Tables 1 and 2). Although these 59 transcripts (29 + 30) corresponding to genes present in the genome of both bacteria, 21 from P. temperata and 22 from X. koppenhoeferi were identified to be exclusively induced in each of them. Only a small number of transcripts were induced in both bacteria, and most of them corresponded to genes required for stress response. Further, the level of changes in the transcripts exclusively induced in one or the other bacterium during infection of R. majalis was determined by quantitative real-time PCR (qRT-PCR) on following representative genes: cysK, def, dnaB, lysR, pchC, pckA, phoP, ptst1, sctL, selD, tolB and virB from P. temperata, and aceK, clpP, cobJ, dacC, dld, malF, res, rtxC, tilS, virH, xhlA, and xkst4 from X. koppenhoeferi. qRT-PCR results displayed consistence with the SCOTS assay. Most tested genes displayed 6-16 fold induction in qRT-PCR assays, and a few genes including sctL, phoP, rtxC and clpP exhibited about 3-fold induction (Fig. 4).

Figure 4
figure 4

Quantitative real time PCR results showing fold changes in the expression of selected Photorhabdus temperata and Xenorhabdus koppenhoeferi genes identified by SCOTS upon infection of the insect host Rhizotrogus majalis.

Analysis of biological associations of in vivo-induced genes

To understand the function and interactions of the identified gene products, interactive networks were built to provide an overview of how gene products relate to each other by leveraging databases of published literature. Using PathwayStudio (Ariadne, Rockville, MD, USA) program, we were able to distill all published information about biological relationships of bacterial homologs of the identified genes. Molecular interactions linking each identified gene product in terms of binding, direct molecular regulation, effects on expression, transportation, and molecule synthesis were elucidated in PathwayStudio. Based on this analysis, the biological significance of each gene was inferred. Linkage between proteins shows how upregulated genes may influence other genes or gene products. For example, gene rpoB, induced in both bacteria, interacts with many gene products, suggesting its overall importance during insect infection (see Additional files 1 and 2). Other identified genes with multiple connections to different gene products were different between P. temperata and X. koppenhoeferi, indicating distinct biological processes used by the two bacteria for insect infection. Moreover, linkages with small molecules indicate how normal functions in gene regulation networks are influenced by small molecules that may be available in the insect body. For example, regulation of genes yncB and purH in P. temperata are influenced by aromatic compounds and formate, respectively (see Additional file 1). In X. koppenhoeferi, gene trxB is influenced by arabinose and selenium, and gene znuA is influenced by etoposide (see Additional file 2). Compared to the biological interactions (see Additional files 1 and 2), molecular pathways or networks in which identified genes link to each other were further illustrated (see Additional file 3). In total, 20 genes from P. temperata and 18 from X. koppenhoeferi (see Additional file 3) were found to be part of the dense protein interaction networks, indicating that these genes may act in concert to achieve successful infection. Some products encoded by induced genes were found to interact with as many as eight other gene products in the network. For example, P. temperata pckA and X. koppenhoeferi aceE both influence 8 other genes in the network (see Additional file 3). As a whole, the protein networks of the identified genes were different for the two bacteria except that NAD synthetase gene nadE was found to be common to both bacteria.

Discussion

Our gene expression data reveal similarities and differences in molecular mechanisms of pathogenicity by the two bacteria with apparently similar biology. As depicted in the conceptual molecular model (Fig. 5), apart from a relatively small number of common strategies used by P. temperata and X. koppenhoeferi to cause infection of R. majalis, most upregulated genes were different between the two bacteria, suggesting vastly different pathways to bring about infection of the same insect host.

Figure 5
figure 5

Conceptual molecular model illustrating comparative contributions of the SCOTS identified genes in bacteria Photorhabdus temperata and Xenorhabdus koppenhoeferi during infection of the insect Rhizotrogus majalis. The set of common and different genes revealed similar and different molecular mechanisms of pathogenicity in P. temperata and X. koppenhoeferi during infection of R. majalis. The common gene products were defined by either sequence similarity or similar functions. The detailed possible functions for individual gene products are listed in Tables 1 and 2.

Genes involved in cell structure modification

As the physical contact between bacteria and their host is accomplished by the outer surface, the bacterial cell envelope is crucial for communication and interaction with host cells during infection. Genes involved in the synthesis of outer membrane proteins (skp in P. temperata and ompF in X. koppenhoeferi) and flagella (fliA in P. temperata and fliM in X. koppenhoeferi) were induced in both bacteria, indicating the importance of surface proteins in this bacteria-insect interface. Outer membrane proteins may contribute to the modulation of bacterial cell surface properties during interaction with the host, thus playing a role in evading host immune response [18]. Flagellum-mediated motility may provide an advantage for infection by attachment to the host surface [19, 20], whereas it has been reported in X. nematophila that the expression of virulence factors is not controlled by the regulation of flagella motility [21]. Gene dacC encoding a peptidoglycan synthesis enzyme D-alanyl-D-alanine carboxypeptidase was induced in X. koppenhoeferi but not in P. temperata, implying that peptidoglycan undergoes structural alterations in X. koppenhoeferi that probably add to its fitness when residing in the host. As a well-defined virulence factor playing a major role in the development of septicemia [22], lipopolysaccharide biosynthesis gene lpsE was induced in X. koppenhoeferi, but not in P. temperata, which may explain in part why relatively slower growing X. koppenhoeferi did not differ from P. temperata in virulence to R. majalis in our previous studies [17].

Genes involved in regulatory networks

Facing stimuli, different regulatory genes were upregulated in P. temperata and X. koppenhoeferi, reflecting their unique ability to respond to environmental challenges. phoP and cpxR genes identified from P. temperata, encode response regulators in two-component regulatory systems PhoPQ and CpxRA, respectively [23, 24]. PhoP and CpxR are dual transcriptional regulators that are activated in response to environmental stimuli and then act to modulate activity of other genes [23, 24]. Responding to low Mg2+ levels, phoP expression is auto-regulated by PhoP and PhoQ proteins and phoQ is constitutive [25]. Auto-regulation of phoP likely regulates expression of other genes that mediate various cellular functions such as LPS modification, cell structure [26], and type III secretion systems [27]. phoP has been previously identified to be involved in virulence as mutational inactivation of this gene rendered P. luminescens avirulent to Spodoptera littoralis larvae [6]. Similarly, the regulatory protein CpxR can be activated in response to signals associated with growth and metabolic pathways [28]. By sensing cell envelope stress CpxR appears to control expression of other genes involved in envelope stress response, secretion, motility and multidrug resistance [29]. Sometimes, CpxR acts to modulate other transcriptional factors [30]. Although there is no evidence showing links between phoP/cpxR and other isolated genes, overexpression of these two response regulator genes may possibly relate to genes involved in stress responses (gshB, gor, msrA and surA), cell structure (skp, dacC and ompF) and type III secretion systems (sctL) identified in this study. Further, in P. temperata, lysR gene encoding a transcriptional activator protein was found to be upregulated during infection of R. majalis, while two genes yijC and rseA encoding transcriptional repressors were identified in X. koppenhoeferi. The pathway network analysis indicates that the product of yijC may function as a repressor to potentially control the expression of the fabB gene, which in turn, modulates the physical properties of the membrane by altering the level of unsaturated fatty acid production. Besides transcriptional regulatory genes, gene tilS, encoding a putative cell cycle protein [4], was induced in X. koppenhoeferi. Upregulation of distinct regulatory genes in these two bacteria indicates the different stress response mechanisms between the bacterial species in the same host. However, none of these regulatory factors were identified in P. luminescens exposed to G. mellonella homogenate in a previous study [16], suggesting that regulation of these factors either depends on host (R. majalis not G. mellonella) or bacteria themselves.

Genes involved in virulence and secretion

Identification of genes encoding hemolysin related proteins in P. temperata (fhaC) and X. koppenhoeferi (xhlA) may suggest their role in dealing with the insect cellular immune response and iron scavenge [31, 32]. Subsequently, bacterial virulence factors actively contribute to a successful infection by colonization of and toxicity towards the insect host. Among the virulence genes identified in this study, tcaC was identified to be induced in both P. temperata and X. koppenhoeferi during infection of R. majalis. The tca gene family encodes four toxin complexes TcaA, TcaB, TcaC and TcaD [33], and has been found in several bacterial pathogens including P. luminescens, Yersinia and Xenorhabdus [34]. Interestingly, Tca toxin encoding genes were also found to be induced in a related species P. luminescens treated with the homogenate of G. mellonella [16]. This leads us to speculate that TcaC serves as either an activator or a chaperone in secretion of toxin complex from the cell as proposed by ffrench-Constant et al. [35] and Waterfield et al. [36]. Induction of gene pplA, encoding a phospholipase, in P. temperata is expected to contribute to its faster growth rate [37] compared to X. koppenhoeferi in R. majalis. Further, genes encoding components of different types of secretion machineries were isolated in this study. Upregulation of gene sctL, encoding components of type III secretion machinery, possibly coordinates protein secretion system in P. temperata for translocating effector molecules across the bacterial cell envelope into host cells [38–40]. Unlike Photorhabdus, the genomes of Xenorhabdus bovienii and X. nematophila do not encode homologues of a dedicated type III secretion system [3], indicating that Xenorhabdus likely uses different effectors or secretion systems. Notably, in X. koppenhoeferi, virH with sequence similarity to a gene encoding putative toxin secretion transporter in Vibrio cholera, was upregulated during infection of R. majalis. In addition, several other genes encoding proteins with putative virulence functions identified in this study may represent novel toxin candidates involved in insect infection. These include pmt1, toxC and virB from P. temperata, and virH from X. koppenhoeferi. Differences in expression of virulence and secretion related genes imply the differential virulence mechanisms employed by the two bacteria.

Genes involved in adaptation to host-induced stresses

Compared to all other functional classes, more stress response genes were identified to be shared by P. temperata and X. koppenhoeferi, suggesting common themes in adaptation to host-imposed stresses. Upregulation of both surA and uspB in both bacteria indicates their importance in bacterial adaptation to global stress during insect infection. Gene surA is important for bacterial survival during infection [41, 42] to avoid accumulation of non-native protein conformations [42–44], and uspB gene product is usually accumulated in the stationary phase and probably important for ethanol stress [45]. Similarly, induction of a chaperone protein encoding gene dnaK in both bacteria may picture common stress responses between P. temperata and X. koppenhoeferi in a number of cellular processes including rescue of misfolded proteins and control of the activity of folded regulatory proteins [46]. The glutathione synthetase gene gshB which mediates condensation of gamma-glutamylcysteine and glycine to form glutathione [47, 48] was induced in P. temperata, but glutathione reductase gene gor which reduces glutathione disulfide to the sulfhydryl form of glutathione [47, 48] was induced in X. koppenhoeferi. Glutathione plays a major role in protection against oxidative stress and in detoxification of hazardous chemicals or heavy metals [49]. Thus, although gshB and gor corresponding products are involved in different glutathione biosynthesis pathways, expression of these two genes may suggests possible role of detoxification mechanism in P. temperata and X. koppenhoeferi during infection of R. majalis. Further, as Photorhabdus and Xenorhabdus bacteria can switch their phenotypes upon insect infection, we speculate that upregulation of transposase genes in both bacteria may contribute to genetic variability for better adaptation and survival in a particular niche [50]. Apart from the common mechanisms handling stresses, gene msrA encoding a peptide methionine sulfoxide reductase was unique to X. koppenhoeferi. As an important antioxidant enzyme, peptide methionine sulfoxide reductase mediates the repair of proteins damaged by sulfoxidation of methionine residues [51, 52]. Products of msrA are found to interact with several other gene products (see Additional file 2), also suggesting its importance during insect infection.

Genes involved in intracellular metabolism

Induction of different metabolism genes in Photorhabdus and Xenorhabdus during insect infection highlights the differences in metabolic changes between the two bacteria. Three genes, aceE, aceK and dld, whose products catalyze subsequent steps of the glyoxylate pathway [53], were induced in X. koppenhoeferi. In contrast, two genes kdgK and pckA involved in energy metabolism in the TCA cycle were upregulated in P. temperata [54, 55]. Genes related to the TCA cycle have been also found to be induced in other bacteria upon host infection, for example in P. luminescens in the presence of G. mellonella homogenate [16] and in Vibrio cholerae upon infection of infant mice [56, 57]. Further, as a cofactor in the metabolism of fatty acids and leucine, biotin plays a role in TCA cycle by assisting in metabolic reactions and helping to transfer carbon dioxide [58]. Therefore, upregulation of bioB encoding a biontin synthetase protein further suggests that TCA cycle is involved in P. temperata. Induction of different amino acid biosynthetic genes may suggest the differences in amino acid requirements in the two bacteria during infection. Induction of cysK gene in P. temperata suggests that cysteine may be required when growing in the insect, but lysine may be needed for X. koppenhoeferi as lysC gene was induced. Isolation of P. temperata def, encoding an amino acid tRNA synthesis enzyme N-formylmethionylaminoacyl-tRNA deformylase, may reflect the importance of specific amino acids involved in the translation process during growth in the insect host. Indeed, the pathway analysis reveals that def gene product interacts with CysS, a cysteinyl-tRNA synthetase (see Additional file 1), also suggesting the requirement for cysteine in P. temperata during infection. Genes involved in protein synthesis and degradation can also be distinguished between the two bacteria. A protease gene, clpP, involved in protein degradation, was induced in X. koppenhoeferi upon infection. In fact proteases have been identified as virulence factors in some other bacterial pathogens [59]. The signal peptidase II encoded by lspA gene induced in P. temperata belongs to a class of aspartyl proteases. Since the product of lspA has been suggested as an important virulence factor operating via maturation of several lipoproteins [60], we postulate that the product of lspA may play a pathological role in P. temperata by cleaving some core proteins and thus influencing bacterial life cycle in the insect host.

Genes involved in nutrient scavenging

Induction of different types of nutrient uptake systems in Photorhabdus and Xenorhabdus further emphasized the possible differences in nutritional requirements of the two bacteria in the insect host. Ion scavenging systems induced in P. temperata and X. koppenhoeferi are completely different. In P. temperata, pchC, a gene involved in the biogenesis of pyochelin type of siderophore [61], was upregulated indicating the importance of iron, whereas a siderophore-independent zinc ion transport system encoded by gene znuA [62] was induced in X. koppenhoeferi indicating the importance of zinc. Also supported by the pathway analysis, product of glnD identified from P. temperata likely controls the production of siderophore (see Additional file 1). In addition to ion scavenging systems, artM, encoding a putative ABC amino acid transport system [63] was identified in P. temperata, but not in X. koppenhoeferi. Besides transporting L-arginine across the inner membrane [63], the pathway network reveals that the product of artM participates in cystine transport (see Additional file 1), suggesting that this gene may also mediate cystine amino acid acquisition. Upregulation of artM in P. temperata also suggests that amino acids may be available for uptake in the insect. Similarly, it has been discussed previously that histidine and phosphatidylethanolamine, the possible components presented in the insect host hemolymph, may be important food sources used by P. luminescens for growth within the insect host [64]. As cysteine synthesis gene cysK was also identified from P. temperata as discussed above, it may indicate that the level of cysteine either synthesized by P. temperata or acquired from the insect host is not sufficient during infection. Interestingly, X. koppenhoeferi seems to require sugar as malF gene which is involved in uptake for maltose and lactose was upregulated. Further, P. temperata but not X. koppenhoeferi induced tolB gene that encodes a Tol protein upon insect infection. As a periplasmic protein, TolB is an extra member in the Tol system. In E. coli, Tol system stabilizes the outer membrane structure and mutations in Tol encoding genes result in hypersensitivity to deleterious agents [65]. Despite unknown mechanistic details of Tol system, upregulation of P. temperata tolB gene suggests that it may be involved in importing biomolecules by somehow supplying energy or transport as noted by Nikaido [66] in other bacteria.

Genes with unknown function

Ten genes encoding products with unknown functions were identified to be upregulated in both bacteria during infection of R. majalis. These include two P. temperata and three X. koppenhoeferi genes which show sequence similarity to genes encoding hypothetical proteins with unknown functions, and two P. temperata and three X. koppenhoeferi sequences that exhibit no similarity to any gene or gene product in current databases. These unknown genes may represent novel pathogenicity factors used by bacteria during insect infection and thus require further investigation.

Conclusion

This research indicates that selective capture of transcribed sequences (SCOTS) is a powerful technique to study in vivo gene expression in bacterial pathogens of insects. Results reveal that phylogenetically related bacteria with similar biology differ dramatically in molecular mechanisms of pathogenicity in the same host. Photorhabdus and Xenorhabdus bacteria appear to share more genes in host-induced stress protection category, including tolerance against universal stresses, glutathione-mediated protection against oxidative stress, and protection against nucleic acid damage, compared to genes in all other categories. While genes encoding TcaC toxin and hemolysin related proteins were identified to be induced in both P. temperata and X. koppenhoeferi, the two bacteria seem to use different regulatory genes, secretion systems and virulence factors to achieve infection of the same insect host R. majalis. Further, P. temperata and X. koppenhoeferi displayed different metabolic adaptations in the same insect host as reflected by the high proportion of identified genes involved in metabolism. Overall, this study provides a better understanding of the evolution of bacterial pathogenesis and provides molecular targets for genetic manipulation of these economically important nematode symbiotic insect pathogens.

Methods

Bacteria, insects, and culture conditions

The symbiotic bacteria P. temperata strain GPS11 [17] and X. koppenhoeferi strain USNJ01 [67] were isolated from nematodes H. bacteriophora strain GPS11 [68] and S. scarabaei strain AMK001 [69], respectively. The H. bacteriophora GPS11 were obtained from our liquid nitrogen frozen stock, and of S. scarabaei AMK001 were obtained from Dr. Albrecht M. Koppenhöfer (Rutgers University, New Brunswick, New Jersey). The bacteria were cultured in BHI broth at 25°C unless otherwise stated. The second instar larvae of the white grub R. majalis were collected from Sunleaf Nursery (Madison, Ohio). The field collected insect larvae were kept at room temperature for 10 days, and only healthy, actively moving larvae were used in all experiments.

General techniques

Bacterial genomic DNA was prepared using standard method for Gram negative bacteria [70]. Biotinylation of bacterial genomic DNA was obtained with EZ-Link Psoralen-PEO-Biotin (Pierce) according to the manufacturer's instructions. The total RNA was isolated using TRIzol reagent (Invitrogen) according to the manufacturer's instructions. RNA samples were treated with RNase-free DNase I (Ambion, Austin, TX) according to the manufacturer's guidelines, and were checked by spectrophotometer and gel electrophoresis. Total RNA was converted to first-strand cDNA using Superscript II reverse transcriptase (Invitrogen RT-PCR kit) according to the manufacturer's instructions. First strand cDNA was made double-strand with Klenow fragment (NEB, Beverly, MA) as described by Froussard [71].

Isolation of bacterial transcripts induced in Rhizotrogus majalis

SCOTS technique was used to identify in vivo gene expression of P. temperata and X. koppenhoeferi during infection of R. majalis. For the control, three 5 μg total RNA samples isolated from 48 h bacterial cultures were used to make double-strand cDNAs using the random primer SCOT09 (5-ATCCACCTATCCCAGTAGGAG NNNNNNNNN). The synthesized cDNAs were PCR amplified using the primer SCOT0 (5-ATCCACCTATCCCAGTAGGAG) for 30 cycles. For the treatment, 10 μl suspension of 48 h P. temperata or X. koppenhoeferi cultures containing 1 × 104 cells was injected into the hemolymph of R. majalis from the base of the foreleg and incubated at 25°C. We evaluated bacterial colonization and found that bacteria of P. temperata and X. koppenhoeferi were 5 × 105 and 2 × 105 colony-forming units (CFU) per insect, respectively, at 24 h post injection, suggesting development of septicemia. Total RNA was isolated from the whole insects at 24 h post injection when the insects were still alive (The insects died within 24 - 48 h). Three 5 μg RNA samples obtained from three infected insects for each bacterium were used to made double-strand by the random primer SCOT189 (5-GACAGATTCGCACTTAACCCT NNNNNNNNN). Double-strand cDNAs were then PCR amplified for 30 cycles by primer SCOT18 (5-GACAGATTCGCACTTAACCCT). Bacterial cDNA normalization was done as described below. In this study, the ribosomal operons (rRNA) of P. temperata and X. koppenhoeferi were amplified and used to block the abundant rRNA sequence in order to effectively capture the cDNA molecules representing mRNA transcripts. The rRNA operon was added to biotinylated genomic DNA at a ratio of 10:1. The genomic DNA-rRNA mixture was sonicated to a size range of 1 to 5 kb. The sonicated, biotinylated genomic DNA-rRNA mixture containing 6 μg rRNA and 0.6 μg genomic DNA was denatured and hybridized for 30 min at 65°C. PCR amplified cDNAs (6 μg) were denatured and added to the genomic DNA-rRNA prehybridized mixture, and hybridized at 65°C for 24 h. Streptavidin magnesphere paramagnetic particles (Invitrogen) were used to capture the bacterial cDNAs that hybridized to biotinylated genomic DNA according to the manufacturer's instructions. Captured cDNAs were then eluted, precipitated, and amplified by PCR with the primer SCOT0 for additional two successive rounds of SCOTS. After three rounds of SCOTS, the three amplified cDNA samples for each bacterium were pooled, and the pooled cDNA mixtures were considered to be normalized in vitro or in vivo cDNA library.

To enrich the normalized in vivo cDNAs, the normalized in vivo cDNAs from each bacterium were hybridized to biotinylated genomic DNA that has been prehybridized with both rRNA and normalized in vitro cDNAs. After hybridization, the bacterial cDNAs were captured and PCR amplified for next round of enrichment. Finally, the amplified enriched bacterial cDNAs were cloned into an original TA cloning vector (Invitrogen) to construct the enriched in vivo cDNA library.

Southern blot screening of SCOTS enriched transcripts

Inserts from the enriched in vivo cDNA library were screened to confirm to be upregulated by southern blot hybridization which has been used as a confirmatory test in most other SCOTS-based studies [72–76]. Individual clones from each subtractive library were randomly picked and amplified by PCR using M13 primers. Ten micro liter of each amplified product was mixed with 70 μl 20 × SSC (1 × SSC is 0.15 M NaCl plus 0.015 M sodium citrate), and then 80 μl sample mixture was transferred to each well of dot-blot containing the positively charged nylon membrane at the low vacuum. The nylon membrane with samples was denatured with denature buffer (3 M NaCl, 0.4 M NaOH) for 10 min at room temperature, and neutralized in 1 × PBS buffer (0.1 M NaCl, 7 mM Na2HPO4, and 3 mM NaH2PO4, pH 6.8) for 10 min at room temperature. The membrane was dried by baking at 80°C for 2 h, and then the baked membrane was rinsed with 2 × SSC and soaked in 2 × SSC for 5 min. The membrane was transferred into a hybridization bottle for hybridization using Dig easy hyb granules (Roche) according to the manufacturer's instruction. Normalized, insect-derived and broth-derived cDNA pools were labeled using the PCR DIG Probe Synthesis kit (Roche) for use as probes. The probes were denatured and added to the hybridization bottles containing the membrane and the hybridization buffer. Hybridization continued at 65°C for approximately 24 h. The membrane was washed briefly with 2 × SSC at room temperature and then twice with 1 × SSC - 0.1% SDS for 15 min each time; these washes were done at 65°C. A final brief rinse with 0.1 × SSC at room temperature completed the washing process. The membrane was incubated at room temperature with 4 ml 1 × SSC with 8% dry milk for 30 m. Four milliliter dilution of anti-digoxigenin-HRP conjugate (1: 800) was added, and incubation continues 1 h at room temperature. The membrane was briefly washed as describe above, and successful hybridization was detected by Amersham ECL Plus western blotting detection reagents (GE healthcare Bio-Sciences Corp, Piscataway, NJ) using chemiluminescent detection. The individual clones that only hybridized to the probe made from normalized in vivo cDNAs were chosen for sequence analysis.

Analysis of bacterial genes induced in Rhizotrogus majalis

The samples were sequenced at the Biotechnology Center, Madison, WI, USA. These sequences were edited and assembled using Editseq (DNASTAR) and Contig (Vector NTI). Similar sequences were identified using BLAST algorithms (blastx and tblastx) in GenBank of National Center for Biotechnology Information http://www.ncbi.nlm.nih.gov and database of Xenorhabdus genomes (http://xenorhabdus.danforthcenter.org. The functions of identified sequences were assigned using Gene Ontology http://www.geneontology.org. The redundancy of each gene was counted cumulatively and the library coverage was calculated using Analytic Rarefaction (http://www.uga.edu/strata/software/Software.html [77] which has been used for estimation of total transcript diversity in the cDNA library [78–81].

Validation of SCOTS procedure

Semi-quantitative reverse transcriptase PCR (RT-PCR) was performed as described previously [82] to verify results obtained by SCOTS analysis. To confirm that the identified genes were differentially induced during infection of the insect hemolymph compared to the culture, RNA was prepared from either infected insect or broth cultured log-phase bacterial cells and DNAse treated. RT-PCR was carried out using 5 μg of RNA in the presence of 5 pmol of gene specific reverse primers and the primer targeting 16S rRNA gene for an internal positive control together in the same reaction. The cDNA template was then diluted 1:100 with sterile ddH2O for use in subsequent PCR. Twenty-five cycles of PCR were performed using primer pairs for specific genes and 16S rRNA gene in separate vials, but from the same master mix, using 5 μl of diluted cDNA as template. In addition, a eukaryotic housekeeping gene 18S rRNA was used as a control to evaluate if bacterial cDNAs were purified apart from the insect cDNAs after SCOTS normalization. For this purpose, presence of 18S rRNA gene in the in vivo cDNA populations before and after SCOTS normalization was measured by PCR amplification using 50 ng cDNA samples and primers 18SF (5-GGAATTGACGGAAGGGCACCA) and 18SR (5-CCAGACAAATCGCTCCACCAAC). Also, a prokaryotic housekeeping gene gyrase A (gyrA) was used as another control to validate that only differentially expressed genes were captured after enrichment. The presence of gyrA in cDNA populations of enriched in vivo cDNAs, and in vitro and in vivo cDNAs before and after SCOTS normalization was evaluated by PCR using primers gyrAF (5-ACGCGACGGTGTACCGGCTT) and gyrAR (5-GCCAGAGAAATCACCCCGGTC). All experiments were performed at least three times.

Comparative hybridization and quantitative real-time PCR

Comparative hybridization was performed to assess exclusivity of the identified genes to P. temperata or X. koppenhoeferi. This approach has been validated for cross-species comparisons of gene expression in many studies [83–86]. To determine the genomic presence of the identified genes, genes identified from one bacterial species were individually screened by cross-hybridization to the sonicated biotinylated genomic DNA of the other, followed by PCR amplification of the eluted hybrids using the defined primers SCOT18. The genes with sequence present in both bacterial genomes were singled out to further evaluate their induction specificity to P. temperata or X. koppenhoeferi. The selected genes induced in one bacterial species were individually screened by cross-hybridization to a heterologous probe made by pre-hybridization of sonicated biotinylated genomic DNA with the respective rRNA operon and normalized in vitro cDNAs of the other, and subsequently amplified using the defined primer SCOT18. Quantitative real-time PCR (qRT-PCR) was further used to validate and quantify the change of genes identified exclusively in one or the other bacterium, excluding genes without sequence similarity between the two bacteria but with similar functions. These representative genes were chosen for analysis based on the availability of designing qRT-PCR primers using identified sequences. Total RNA was prepared from bacteria grown in intro or in vivo as described above. Real-time PCR was performed in an IQ5 machine (Bio-Rad) using QuantiTect SybrGreen PCR Kit (Qiagen) according to the manufacturer's instructions. The specific gene primers for real-time PCR were listed in the Additional file 4. All reactions were run in triplicate with three independent cDNA samples and a no cDNA template negative control. According to a previous study on Photorhabdus bacteria [87], 16S rRNA gene was chosen as an internal control for normalization. The calculated threshold cycle (Ct) was normalized to the Ct of the 16s rRNA gene amplified from the corresponding sample, and the relative change in expression of a gene in vivo as compared to in vitro was calculated as previously described [88].

Pathway analysis

To depict biological interactions influenced by identified genes, we built information rich predicted interaction networks using the ResnetCore bacterial molecular database of the PathwayStudio 5.0 (Ariadne, Rockville, MD, USA). PathwayStudio is a program for visualization and analysis of biological pathways and gene regulation networks. This program comes with the ResNet database of more than 500,000 functional relationships and the MedScan tool for automatic extraction of information from the scientific literature. First, the direct interaction with proteins and small molecules for each identified gene was visualized using PathwayStudio. Then, the interaction network of identified genes was built based on relationships between gene products or small molecules in form of binding, direct molecular regulation, effects on expression, molecule transport, protein modification, and molecule synthesis.

Abbreviations

IJs:

infective juveniles

SCOTS:

selective capture of transcribed sequences.

References

  1. Boemare N, Akhurst P: The genera Photorhabdus and Xenorhabdus. Prokaryotes. 2006, 6: 451-494. full_text.

    Google Scholar 

  2. Forst S, Dowds B, Boemare N, Stackebrandt E: Xenorhabdus and Photorhabdus spp.: bugs that kill bugs. Annu Rev Microbiol. 1997, 51: 47-72. 10.1146/annurev.micro.51.1.47.

    Article  CAS  PubMed  Google Scholar 

  3. Goodrich-Blair H, Clarke DJ: Mutualism and pathogenesis in Xenorhabdus and Photorhabdus: two roads to the same destination. Mol Microbiol. 2007, 64: 260-268. 10.1111/j.1365-2958.2007.05671.x.

    Article  CAS  PubMed  Google Scholar 

  4. Duchaud E, Rusniok C, Frangeul L, Buchrieser C, Givaudan A, Taourit S, Bocs S, Boursaux-Eude C, Chandler M, Charles JF, Dassa E, Derose R, Derzelle S, Freyssinet G, Gaudriault S, Medigue C, Lanois A, Powell K, Siguier P, Vincent R, Wingate V, Zouine M, Glaser P, Boemare N, Danchin A, Kunst F: The genome sequence of the entomopathogenic bacterium Photorhabdus luminescens. Nat Biotechnol. 2003, 21: 1307-1313. 10.1038/nbt886.

    Article  CAS  PubMed  Google Scholar 

  5. Waterfield NR, Daborn PJ, ffrench-Constant RH: Genomic islands in Photorhabdus. Trends Microbiol. 2002, 10: 541-545. 10.1016/S0966-842X(02)02463-0.

    Article  CAS  PubMed  Google Scholar 

  6. Derzelle S, Turlin E, Duchaud E, Pages S, Kunst F, Givaudan A, Danchin A: The PhoP-PhoQ two-component regulatory system of Photorhabdus luminescens is essential for virulence in insects. J Bacteriol. 2004, 186: 1270-1279. 10.1128/JB.186.5.1270-1279.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Bennett HP, Clarke DJ: The pbgPE operon in Photorhabdus luminescens is required for pathogenicity and symbiosis. J Bacteriol. 2005, 187: 77-84. 10.1128/JB.187.1.77-84.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Krin E, Chakroun N, Turlin E, Givaudan A, Gaboriau F, Bonne I, Rousselle JC, Frangeul L, Lacroix C, Hullo MF, Marisa L, Danchin A, Derzelle S: Pleiotropic role of quorum-sensing autoinducer 2 in Photorhabdus luminescens. Appl Environ Microbiol. 2006, 72: 6439-6451. 10.1128/AEM.00398-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Eleftherianos I, Boundy S, Joyce SA, Aslam S, Marshall JW, Cox RJ, Simpson TJ, Clarke DJ, ffrench-Constant RH, Reynolds SE: An antibiotic produced by an insect-pathogenic bacterium suppresses host defenses through phenoloxidase inhibition. PNAS. 2007, 104: 2419-2424. 10.1073/pnas.0610525104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Ji D, Kim Y: An entomopathogenic bacterium, Xenorhabdus nematophila, inhibits the expression of an antibacterial peptide, cecropin, of the beet armyworm, Spodoptera exigua. J Insect Physiol. 2004, 50: 489-496. 10.1016/j.jinsphys.2004.03.005.

    Article  CAS  PubMed  Google Scholar 

  11. Park Y, Herbert EE, Cowles CE, Cowles KN, Menard ML, Orchard SS, Goodrich-Blair H: Clonal variation in Xenorhabdus nematophila virulence and suppression of Manduca sexta immunity. Cell Microbiol. 2007, 9: 645-656. 10.1111/j.1462-5822.2006.00815.x.

    Article  CAS  PubMed  Google Scholar 

  12. Cowles KN, Cowles CE, Richards GR, Martens EC, Goodrich-Blair H: The global regulator Lrp contributes to mutualism, pathogenesis and phenotypic variation in the bacterium Xenorhabdus nematophila. Cell Microbiol. 2007, 9: 1311-1323. 10.1111/j.1462-5822.2006.00873.x.

    Article  CAS  PubMed  Google Scholar 

  13. Brugirard-Ricaud K, Givaudan A, Parkhill J, Boemare N, Kunst F, Zumbihl R, Duchaud E: Variation in the effectors of the type III secretion system among Photorhabdus species as revealed by genomic analysis. J Bacteriol. 2004, 186: 4376-4381. 10.1128/JB.186.13.4376-4381.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Brugirard-Ricaud K, Duchaud E, Givaudan A, Girard PA, Kunst F, Boemare N, Brehelin M, Zumbihl R: Site-specific antiphagocytic function of the Photorhabdus luminescens type III secretion system during insect colonization. Cell Microbiol. 2005, 7: 363-3671. 10.1111/j.1462-5822.2004.00466.x.

    Article  CAS  PubMed  Google Scholar 

  15. Park D, Forst S: Co-regulation of motility, exoenzyme and antibiotic production by the EnvZ-OmpR-FlhDC-FliA pathway in Xenorhabdus nematophila. Mol Microbiol. 2006, 61: 1397-1412. 10.1111/j.1365-2958.2006.05320.x.

    Article  CAS  PubMed  Google Scholar 

  16. Münch A, Stingl L, Jung K, Heermann R: Photorhabdus luminescens genes induced upon insect infection. BMC Genomics. 2008, 9: 229-

    Article  PubMed Central  PubMed  Google Scholar 

  17. An R, Grewal PS: Differences in the virulence of Heterorhabditis bacteriophora and Steinernema scarabaei to three white grub species: The relative contribution of the nematodes and their symbiotic bacteria. Biol Contr. 2007, 43: 310-316. 10.1016/j.biocontrol.2007.07.004.

    Article  Google Scholar 

  18. Rediers H, Rainey PB, Vanderleyden J, Mot RD: Unraveling the secret lives of bacteria: use of in vivo expression technology and differential fluorescence induction promoter traps as tools for exploring niche-specific gene expression. Microbiol Mol Biol Rev. 2005, 69: 217-261. 10.1128/MMBR.69.2.217-261.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Moens S, Vanderleyden J: Functions of bacterial flagella. Crit Rev Microbiol. 1996, 22: 67-100. 10.3109/10408419609106456.

    Article  CAS  PubMed  Google Scholar 

  20. Dalton HM, March PE: Molecular genetics of bacterial attachment and biofouling. Curr Opin Biotechnol. 1998, 9: 252-255. 10.1016/S0958-1669(98)80055-4.

    Article  CAS  PubMed  Google Scholar 

  21. Lanois A, Jubelin G, Givaudan A: FliZ, a flagellar regulator, is at the crossroads between motility, haemolysin expression and virulence in the insect pathogenic bacterium Xenorhabdus. Mol Microbiol. 2008, 68: 516-533. 10.1111/j.1365-2958.2008.06168.x.

    Article  CAS  PubMed  Google Scholar 

  22. Mayeux PR: Pathobiology of lipopolysaccharide. J Toxicol Environ Health. 1997, 51: 415-435. 10.1080/00984109708984034.

    Article  CAS  PubMed  Google Scholar 

  23. De WP, Kwon O, Lin EC: The CpxRA signal transduction system of Escherichia coli: growth-related autoactivation and control of unanticipated target operons. J Bacteriol. 1999, 181: 6772-6778.

    Google Scholar 

  24. Kato A, Tanabe H, Utsumi R: Molecular characterization of the PhoP-PhoQ two-component system in Escherichia coli K-12: identification of extracellular Mg2+-responsive promoters. J Bacteriol. 1999, 181: 5516-5520.

    PubMed Central  CAS  PubMed  Google Scholar 

  25. Soncini FC, Vescovi EG, Groisman EA: Transcriptional autoregulation of the Salmonella typhimurium PhoPQ operon. J Bacteriol. 1995, 177: 4364-4371.

    PubMed Central  CAS  PubMed  Google Scholar 

  26. Groisman EA: The pleiotropic two-component regulatory system PhoP-PhoQ. J Bacteriol. 2001, 183: 1835-1842. 10.1128/JB.183.6.1835-1842.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Francis MS, Wolf-Watz H, Forsberg A: Regulation of type III secretion systems. Curr Opin Microbiol. 2002, 5: 166-172. 10.1016/S1369-5274(02)00301-6.

    Article  CAS  PubMed  Google Scholar 

  28. Wolfe AJ, Parikh N, Lima BP, Zemaitaitis B: Signal integration by the two-component signal transduction response regulator CpxR. J Bacteriol. 2008, 190: 2314-2322. 10.1128/JB.01906-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Raivio TL: Envelope stress responses and Gram-negative bacterial pathogenesis. Mol Microbiol. 2005, 56: 1119-1128. 10.1111/j.1365-2958.2005.04625.x.

    Article  CAS  PubMed  Google Scholar 

  30. Dorel D, Lejeune P, Rodrigue A: The Cpx system of Escherichia coli, a strategic signaling pathway for confronting adverse conditions and for settling biofilm communities?. Res Microbiol. 2006, 157: 306-314. 10.1016/j.resmic.2005.12.003.

    Article  CAS  PubMed  Google Scholar 

  31. Brillard J, Ribeiro C, Boemare N, Brehelin M, Givaudan A: Two distinct hemolytic activities in Xenorhabdus nematophila are active against immunocompetent insect cells. Appl Environ Microbiol. 2001, 67: 2515-2525. 10.1128/AEM.67.6.2515-2525.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Brillard J, Duchaud E, Boemare N, Kunst F, Givaudan A: The PhlA hemolysin from the entomopathogenic bacterium Photorhabdus luminescens belongs to the two-partner secretion family of hemolysins. J Bacteriol. 2002, 184: 3871-3878. 10.1128/JB.184.14.3871-3878.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Bowen D, Rocheleau TA, Blackburn M, Andreev O, Golubeva E, Bhartia R, ffrench-Constant RH: Insecticidal toxins from the bacterium Photorhabdus luminescens. Science. 1998, 280: 2129-2132. 10.1126/science.280.5372.2129.

    Article  CAS  PubMed  Google Scholar 

  34. ffrench-Constant RH, Bowen D: Photorhabdus toxins: novel biological insecticides. Curr Opin Microbiol. 1999, 2: 284-288. 10.1016/S1369-5274(99)80049-6.

    Article  CAS  PubMed  Google Scholar 

  35. ffrench-Constant R, Waterfield N, Daborn P, Joyce S, Bennett H, Au C, Dowling A, Boundy S, Reynolds S, Clarke D: Photorhabdus: towards a functional genomic analysis of a symbiont and pathogen. FEMS Microbiol Rev. 2003, 26: 433-456. 10.1111/j.1574-6976.2003.tb00625.x.

    Article  CAS  PubMed  Google Scholar 

  36. Waterfield NR, Bowen DJ, Fetherston JD, Perry RD, ffrench-Constant RH: The tc genes of Photorhabdus: a growing family. Trends Microbiol. 2001, 9: 185-191. 10.1016/S0966-842X(01)01978-3.

    Article  CAS  PubMed  Google Scholar 

  37. Istivan TS, Coloe PJ: Phospholipase A in Gram-negative bacteria and its role in pathogenesis. Microbiology-Sgm. 2006, 152: 1263-1274. 10.1099/mic.0.28609-0.

    Article  CAS  Google Scholar 

  38. Alfano JR, Collmer A: Type III secretion system effector proteins: Double agents in bacterial disease and plant defense. Ann Rev Phytopathol. 2004, 42: 385-414. 10.1146/annurev.phyto.42.040103.110731.

    Article  CAS  Google Scholar 

  39. Buttner D, Bonas U: Common infection strategies of plant and animal pathogenic bacteria. Curr Opin Plant Biol. 2003, 6: 312-319. 10.1016/S1369-5266(03)00064-5.

    Article  PubMed  Google Scholar 

  40. Collmer A, Lindeberg M, Petnicki-Ocwieja T, Schneider DJ, Alfano JR: Genomic mining type III secretion system effectors in Pseudomonas syringae yields new picks for all TTSS prospectors. Trends Microbiol. 2002, 10: 462-469. 10.1016/S0966-842X(02)02451-4.

    Article  CAS  PubMed  Google Scholar 

  41. Leuzzi R, Serino L, Scarselli M, Savino S, Fontana MR, Monaci E, Taddei A, Fischer G, Rappuoli R, Pizza M: Ng-MIP, a surface-exposed lipoprotein of Neisseria gonorrhoeae, has a peptidyl-prolyl cis/trans isomerase (PPIase) activity and is involved in persistence in macrophages. Mol Microbiol. 2005, 58: 669-681. 10.1111/j.1365-2958.2005.04859.x.

    Article  CAS  PubMed  Google Scholar 

  42. Ren P, Rossettini A, Chaturvedi V, Hanes SD: The Ess1 prolyl isomerase is dispensable for growth but required for virulence in Cryptococcus neoformans. Microbiology-Sgm. 2005, 151: 1593-1605. 10.1099/mic.0.27786-0.

    Article  CAS  Google Scholar 

  43. Fischer G, Tradler T, Zarnt T: The mode of action of peptidyl prolyl cis/trans isomerases in vivo: binding vs. catalysis. FEBS Lett. 1998, 426: 17-20. 10.1016/S0014-5793(98)00242-7.

    Article  CAS  PubMed  Google Scholar 

  44. Schiene C, Fischer G: Enzymes that catalyse the restructuring of proteins. Curr Opin Struct Biol. 2000, 10: 40-45. 10.1016/S0959-440X(99)00046-9.

    Article  CAS  PubMed  Google Scholar 

  45. Farewell A, Kvint K, Nystrom T: uspB, a new sigma S-regulated gene in Escherichia coli which is required for stationary-phase resistance to ethanol. J Bacteriol. 1998, 180: 6140-6147.

    PubMed Central  CAS  PubMed  Google Scholar 

  46. Bukau B, Horwich AL: The Hsp70 and Hsp60 chaperone machines. Cell. 1998, 92: 351-366. 10.1016/S0092-8674(00)80928-9.

    Article  CAS  PubMed  Google Scholar 

  47. Meister A: Glutathione metabolism and its selective modification. J Biol Chem. 1988, 263: 17205-17208.

    CAS  PubMed  Google Scholar 

  48. Mannervik B: The enzymes of glutathione metabolism - an overview. Biochem Soc Trans. 1987, 15: 717-718.

    Article  CAS  PubMed  Google Scholar 

  49. Hayes JD, Flanagan JU, Jowsey IR: Glutathione transferases. Ann Rev Pharmacol Toxicol. 2005, 45: 51-88. 10.1146/annurev.pharmtox.45.120403.095857.

    Article  CAS  Google Scholar 

  50. Morschhauser J, Kohler G, Ziebuhr W, Blum-Oehler G, Dobrindt U, Hacker J: Evolution of microbial pathogens. Philos Trans R Soc Lond B Biol Sci. 2000, 355: 695-704. 10.1098/rstb.2000.0609.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Alamuri P, Maier RJ: Methionine sulphoxide reductase is an important antioxidant enzyme in the gastric pathogen Helicobacter pylori. Mol Microbiol. 2004, 53: 1397-1406. 10.1111/j.1365-2958.2004.04190.x.

    Article  CAS  PubMed  Google Scholar 

  52. Weissbach H, Etienne F, Hoshi T, Heinemann SH, Lowther WT, Matthews B, St John G, Nathan C, Brot N: Peptide methionine sulfoxide reductase: Structure, mechanism of action, and biological function. Arch Biochem Biophys. 2002, 397: 172-178. 10.1006/abbi.2001.2664.

    Article  CAS  PubMed  Google Scholar 

  53. Ikeda T, Laporte DC: Isocitrate dehydrogenase kinase phosphatase - Acek alleles that express kinase but not phosphatase-activity. J Bacteriol. 1991, 173: 1801-1806.

    PubMed Central  CAS  PubMed  Google Scholar 

  54. Cynkin MA, Ashwell G: Uronic acid metabolism in bacteria. IV. Purification and properties of 2-keto-3-deoxy-D-gluconokinase in Escherichia coli. J Biol Chem. 1960, 235: 1576-1579.

    CAS  PubMed  Google Scholar 

  55. Delbaere LT, Sudom AM, Prasad L, Leduc Y, Goldie H: Structure/function studies of phosphoryl transfer by phosphoenolpyruvate carboxykinase. Biochim Biophys Acta. 2004, 1697: 271-278.

    Article  CAS  PubMed  Google Scholar 

  56. Camilli A, Mekalanos JJ: Use of recombinase gene fusions to identify Vibrio cholerae genes induced during infection. Mol Microbiol. 1995, 18: 671-683. 10.1111/j.1365-2958.1995.mmi_18040671.x.

    Article  CAS  PubMed  Google Scholar 

  57. Osorio CG, Crawford JA, Michalski J, Martinez-Wilson H, Kaper JB, Camilli A: Second-generation recombination-based in vivo expression technology for large-scale screening for Vibrio cholerae genes induced during infection of the mouse small intestine. Infect Immun. 2005, 73: 972-980. 10.1128/IAI.73.2.972-980.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Sanyal I, Cohen G, Flint DH: Biotin synthase: purification, characterization as a [2Fe-2S] cluster protein, and in vitro activity of the Escherichia coli bioB gene product. Biochemistry. 1994, 33: 3625-3631. 10.1021/bi00178a020.

    Article  CAS  PubMed  Google Scholar 

  59. Travis J, Potempa J, Maeda H: Are bacterial proteinases pathogenic factors. Trends Microbiol. 1995, 3: 405-407. 10.1016/S0966-842X(00)88988-X.

    Article  CAS  PubMed  Google Scholar 

  60. Geukens N, De Buck E, Meyen E, Maes L, Vranckx L, Van Mellaert L, Anné J, Lammertyn E: The type II signal peptidase of Legionella pneumophila. Res Microbiol. 2006, 157: 836-841. 10.1016/j.resmic.2006.06.003.

    Article  CAS  PubMed  Google Scholar 

  61. Eppinger M, Rosovitz MJ, Fricke WF, Rasko DA, Kokorina G, Fayolle C, Lindler LE, Carniel E, Ravel J: The complete genome sequence of Yersinia pseudotuberculosis IP31758, the causative agent of Far East scarlet-like fever. PloS Genet. 2007, 3: e142-10.1371/journal.pgen.0030142.

    Article  PubMed Central  PubMed  Google Scholar 

  62. Patzer SI, Hantke K: The ZnuABC high-affinity zinc uptake system and its regulator Zur in Escherichia coli. Mol Microbiol. 1998, 28: 1199-1210. 10.1046/j.1365-2958.1998.00883.x.

    Article  CAS  PubMed  Google Scholar 

  63. Linton KJ, Higgins CF: The Escherichia coli ATP-binding cassette (ABC) proteins. Mol Microbiol. 1998, 28: 5-13. 10.1046/j.1365-2958.1998.00764.x.

    Article  CAS  PubMed  Google Scholar 

  64. Heermann R, Fuchs TM: Comparative analysis of the Photorhabdus luminescens and the Yersinia enterocolitica genomes: uncovering candidate genes involved in insect pathogenicity. BMC Genomics. 2008, 9: 40-10.1186/1471-2164-9-40.

    Article  PubMed Central  PubMed  Google Scholar 

  65. Lazzaroni JC, Germon P, Ray MC, Vianney A: The Tol proteins of Escherichia coli and their involvement in the uptake of biomolecules and outer membrane stability. FEMS Microbiol Lett. 1999, 177: 191-197. 10.1111/j.1574-6968.1999.tb13731.x.

    Article  CAS  PubMed  Google Scholar 

  66. Nikaido H: Molecular basis of bacterial outer membrane permeability revisited. Microbiol Mol Biol Rev. 2003, 67: 593-656. 10.1128/MMBR.67.4.593-656.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  67. Tailliez P, Pages S, Ginibre N, Boemare N: New insight into diversity in the genus Xenorhabdus, including the description of ten novel species. Int J Syst Evol Microbiol. 2006, 56: 2805-2818. 10.1099/ijs.0.64287-0.

    Article  CAS  PubMed  Google Scholar 

  68. Grewal PS, Grewal SK, Malik VS, Klein MG: Differences in susceptibility of introduced and native white grub species to entomopathogenic nematodes from various geographic localities. Biol Contr. 2002, 24: 230-237. 10.1016/S1049-9644(02)00025-7.

    Article  Google Scholar 

  69. Stock SP, Koppenhöfer AM: Steinernema scarabaei n. sp. (Rhabditida: Steinernematidae), a natural pathogen of scarab beetle larvae (Coleoptera: Scarabaeidae) from New Jersey, USA. Nematology. 2003, 5: 191-204. 10.1163/156854103767139680.

    Article  Google Scholar 

  70. Sambrook J, Russell DW, Russell D: Plasmids and their usefulness in molecular cloning. Molecular cloning: A laboratory manual. 2000, New York: Cold Spring Harbor Laboratory Press, 7-13. 3

    Google Scholar 

  71. Froussard P: A random-PCR method (rPCR) to construct whole cDNA library from low amounts of RNA. Nucl Acids Res. 1992, 20: 2900-10.1093/nar/20.11.2900.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Baltes N, Gerlach GF: Identification of genes transcribed by Actinobacillus pleuropneumoniae in necrotic porcine lung tissue by using selective capture of transcribed sequences. Infect Immun. 2004, 72: 6711-6716. 10.1128/IAI.72.11.6711-6716.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  73. Daigle F, Graham JE, Curtiss R: Identification of Salmonella typhi genes expressed within macrophages by selective capture of transcribed sequences (SCOTS). Mol Microbiol. 2001, 41: 1211-1222. 10.1046/j.1365-2958.2001.02593.x.

    Article  CAS  PubMed  Google Scholar 

  74. Graham JE, Clark-Curtiss JE: Identification of Mycobacterium tuberculosis RNAs synthesized in response to phagocytosis by human macrophages by selective capture of transcribed sequences (SCOTS). Proc Natl Acad Sci USA. 1999, 96: 11554-11559. 10.1073/pnas.96.20.11554.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  75. Haydel SE, Benjamin WH, Dunlap NE, Clark-Curtiss JE: Expression, autoregulation, and DNA binding properties of the Mycobacterium tuberculosis TrcR response regulator. J Bacteriol. 2002, 184: 2192-2203. 10.1128/JB.184.8.2192-2203.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Hou JY, Graham JE, Clark-Curtiss JE: Mycobacterium avium genes expressed during growth in human macrophages detected by selective capture of transcribed sequences (SCOTS). Infect Immun. 2002, 70: 3714-3726. 10.1128/IAI.70.7.3714-3726.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  77. Heck KL, Belle Gv, Simberloff D: Explicit calculation of the rarefaction diversity measurement and the determination of sufficient sample size. Ecology. 1975, 56: 1459-1461. 10.2307/1934716.

    Article  Google Scholar 

  78. Wang YL, Morse D: Rampant polyuridylylation of plastid gene transcripts in the dinoflagellate Lingulodinium. Nucleic Acids Res. 2006, 34: 613-619. 10.1093/nar/gkj438.

    Article  PubMed Central  PubMed  Google Scholar 

  79. Suga K, Mark Welch D, Tanaka Y, Sakakura Y, Hagiwara A: Analysis of expressed sequence tags of the cyclically parthenogenetic rotifer Brachionus plicatilis. PLoS ONE. 2007, 2: e671-10.1371/journal.pone.0000671.

    Article  PubMed Central  PubMed  Google Scholar 

  80. Zhu XC, Tu ZJ, Coussens PM, Kapur V, Janagama H, Naser S, Sreevatsan S: Transcriptional analysis of diverse strains Mycobacterium avium subspecies paratuberculosis in primary bovine monocyte derived macrophages. Microb Infect. 2008, 10: 1274-1282. 10.1016/j.micinf.2008.07.025.

    Article  CAS  Google Scholar 

  81. Frias-Lopez J, Shi Y, Tyson GW, Coleman ML, Schuster SC, Chisholm SW, DeLong EF: Microbial community gene expression in ocean surface waters. PNAS. 2008, 105: 3805-3810. 10.1073/pnas.0708897105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  82. Baltes N, Buettner FF, Gerlach GF: Selective capture of transcribed sequences (SCOTS) of Actinobacillus pleuropneumoniae in the chronic stage of disease reveals an HlyX-regulated autotransporter protein. Vet Microbiol. 2007, 123: 110-121. 10.1016/j.vetmic.2007.03.026.

    Article  CAS  PubMed  Google Scholar 

  83. Mazur BJ, Rice D, Haselkorn R: Identification of blue-green-algal nitrogen-fixation genes by using heterologous DNA hybridization probes. Proceedings of the National Academy of Sciences of the United States of America-Biological Sciences. 1980, 77: 186-190. 10.1073/pnas.77.1.186.

    Article  CAS  Google Scholar 

  84. Renn SCP, Aubin-Horth N, Hofmann HA: Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray. BMC Genomics. 2004, 5: 42-10.1186/1471-2164-5-42.

    Article  PubMed Central  PubMed  Google Scholar 

  85. Bar-Or C, Bar-Eyal M, Gal TZ, Kapulnik Y, Czosnek H, Koltai H: Derivation of species-specific hybridization-like knowledge out of cross-species hybridization results. BMC Genomics. 2006, 7: 110-10.1186/1471-2164-7-110.

    Article  PubMed Central  PubMed  Google Scholar 

  86. Buckley BA: Comparative environmental genomics in non-model species: using heterologous hybridization to DNA-based microarrays. J Exp Biol. 2007, 210: 1602-1606. 10.1242/jeb.002402.

    Article  CAS  PubMed  Google Scholar 

  87. Daborn PJ, Waterfield N, Blight MA, ffrench-Constant RH: Measuring virulence factor expression by the pathogenic bacterium Photorhabdus luminescens in culture and during insect infection. J Bacteriol. 2001, 183: 5834-5839. 10.1128/JB.183.20.5834-5839.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  88. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. Ganpati Jagdale for technical assistance. This project was supported by an interdisciplinary grant from the Ohio Agricultural Research and Development Center, Wooster, Ohio.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Parwinder S Grewal.

Additional information

Authors' contributions

RA performed all the experiments, analyzed the data, wrote the first draft, and participated in the design of the study. SS helped in the design of the study and provided comments on the manuscript. PSG obtained the funding for the research, designed the study, and helped in writing of the paper. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2008_2317_MOESM1_ESM.tiff

Additional file 1: Direct linkages of proteins and small molecules to the SCOTS identified Photorhabdus temperata genes induced upon infection of Rhizotrogus majalis. Linkages of molecules to identified genes were built in the PathwayStudio program. Gray and green ovals indicate small molecules. The gene products are represented by red or blue ovals, where blue ovals indicate genes identified in this study and red ovals represent genes in the database of the PathwayStudio program. The relationships are indicated by lines as follows: Binding - violet links with violet circles, MolTransport - gray arrows with green rectangles, MolSynthesis - blue arrows with blue rectangles, ProtModification - brown arrows, Regulation - gray links with gray rectangles, PromoterBinding - green arrows with green circles, and Expression - blue arrows with blue rectangles. Arrows with "+" indicate positive regulation and with "-" indicate negative regulation. (TIFF 12 MB)

12864_2008_2317_MOESM2_ESM.tiff

Additional file 2: Direct linkages of proteins and small molecules to the SCOTS identified Xenorhabdus koppenhoeferi genes induced upon infection of Rhizotrogus majalis. Linkages of molecules to identified genes were built in the PathwayStudio program. Gray and green ovals indicate small molecules. The gene products are represented by red or blue ovals, where blue ovals indicate genes identified in this study and red ovals represent genes in the database of the PathwayStudio program. The relationships are indicated by lines as follows: Binding - violet links with violet circles, MolTransport - gray arrows with green rectangles, MolSynthesis - blue arrows with blue rectangles, ProtModification - brown arrows, Regulation - gray links with gray rectangles, PromoterBinding - green arrows with green circles, and Expression - blue arrows with blue rectangles. Arrows with "+" indicate positive regulation and with "-" indicate negative regulation. (TIFF 13 MB)

12864_2008_2317_MOESM3_ESM.tiff

Additional file 3: The network interactions between gene products of Photorhabdus temperata (panel A) and Xenorhabdus koppenhoeferi (panel B) built by leveraging databases of published literature. The gene expression data obtained in this study were analyzed and visualized by PathwayStudio program. The gene products are represented by ovals, where blue ovals represent products encoded by genes induced in P. temperata (panel A) and X. koppenhoeferi (panel B) upon infection of the insect Rhizotrogus majalis and red ovals represent genes in the database of the PathwayStudio program. The product encoded by gene nadE was highlighted with a green circle as this gene was common to both bacteria and interacts with many other gene products in both networks. The relationships between gene products are indicated by lines as follows: Binding - violet links with violet circles, MolTransport - gray arrows with green rectangles, MolSynthesis - blue arrows with blue rectangles, ProtModification - brown arrows, Regulation - gray links with gray rectangles, PromoterBinding - green arrows with green circles, and Expression - blue arrows with blue rectangles. Arrows with "+" indicate positive regulation and with "-" indicate negative regulation. (TIFF 8 MB)

12864_2008_2317_MOESM4_ESM.pdf

Additional file 4: Supplementary Table 1. Oligonucleotide sequences used for quantitative real-time PCR analyses. (PDF 8 KB)

Authors’ original submitted files for images

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

An, R., Sreevatsan, S. & Grewal, P.S. Comparative in vivo gene expression of the closely related bacteria Photorhabdus temperata and Xenorhabdus koppenhoeferi upon infection of the same insect host, Rhizotrogus majalis. BMC Genomics 10, 433 (2009). https://doi.org/10.1186/1471-2164-10-433

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-433

Keywords