Systems analysis of multiple regulator perturbations allows discovery of virulence factors in Salmonella
- Hyunjin Yoon†1,
- Charles Ansong†2,
- Jason E McDermott2,
- Marina Gritsenko2,
- Richard D Smith2,
- Fred Heffron1 and
- Joshua N Adkins2Email author
© Yoon et al; licensee BioMed Central Ltd. 2011
Received: 24 March 2011
Accepted: 28 June 2011
Published: 28 June 2011
Systemic bacterial infections are highly regulated and complex processes that are orchestrated by numerous virulence factors. Genes that are coordinately controlled by the set of regulators required for systemic infection are potentially required for pathogenicity.
In this study we present a systems biology approach in which sample-matched multi-omic measurements of fourteen virulence-essential regulator mutants were coupled with computational network analysis to efficiently identify Salmonella virulence factors. Immunoblot experiments verified network-predicted virulence factors and a subset was determined to be secreted into the host cytoplasm, suggesting that they are virulence factors directly interacting with host cellular components. Two of these, SrfN and PagK2, were required for full mouse virulence and were shown to be translocated independent of either of the type III secretion systems in Salmonella or the type III injectisome-related flagellar mechanism.
Integrating multi-omic datasets from Salmonella mutants lacking virulence regulators not only identified novel virulence factors but also defined a new class of translocated effectors involved in pathogenesis. The success of this strategy at discovery of known and novel virulence factors suggests that the approach may have applicability for other bacterial pathogens.
The interactions between intracellular pathogen and host can be complex involving sophisticated offensive and defensive strategies by both organisms. Developing a systems level understanding of the virulence program of a pathogen, both in terms of the regulatory pathways and the virulence-related proteins that execute this program is important to effectively combat persistent and adapting pathogens [1–3]. Combining high-throughput characterization of proteins and gene transcripts under multiple different conditions relevant to virulence provides a wealth of information that can be mined to provide useful leads for further investigation or used as the basis of predictive models.
Salmonella enterica serovar Typhimurium (STM) is a facultative intracellular bacterial pathogen with a broad host range capable of infecting birds, reptiles, mice, humans and other mammals. In humans, it is a leading causative agent of gastroenteritis with significant impacts on childhood mortality in the developing world  and among HIV positive patients in Sub-Saharan Africa [4, 5]. In susceptible mice, STM causes a lethal systemic infection. Because its symptoms resemble human typhoid fever caused by the S. enterica serovar Typhi, which only infects man, STM-mediated systemic infection in mice represents an established model system to investigate the pathogenesis and immunology of typhoid fever in humans [6, 7]. STM is an intracellular pathogen that can replicate in a variety of cell types, but is most frequently found in monocytes and neutrophils following infection where it is located within a specialized host-membrane bound body, the Salmonella-containing vacuole (SCV) [8, 9]. STM is able to avoid lysosomal fusion with the SCV and thus evade destruction [10, 11].
During systemic infection a sophisticated regulatory network processes and integrates a variety of hostile environmental cues including acidic pH , antimicrobial peptides  and reactive oxygen species  within macrophages triggering induction of specific subsets of genes to adapt to the growth environment, evade the innate immune response, and prevent lysosomal fusion. The coordinated action by this regulatory network involves tight regulation of numerous virulence-related factors including the Salmonella Pathogenicity Island-2 (SPI-2) type III secretion system (T3SS)  and secreted proteins. SPI-2, located at 31 minutes on the chromosome, encodes the structural components of SPI-2 T3SS as well as regulators and secretion effectors. SPI-2 T3SS delivers a variety of effector proteins into host cells to manipulate host cellular activities. Functions of SPI-2 effector proteins are diverse, including blocking phago-lysosome fusion, subverting inflammatory response, and modulating motility of infected cells .
The identification of genes that are required for STM virulence using high-throughput "omics" technologies is an area of considerable experimental investigation [16–20]. The general approach employed in most of these studies involves profiling STM mRNA or protein abundances following infection of host cells or under infection-mimicking conditions versus standard laboratory conditions, with pathogen genes exhibiting higher mRNA or protein abundance under infectious conditions identified as potential virulence factors. However a challenge with this approach is deconvoluting virulence-relevant and -irrelevant effects. Omics methods applied to identify virulence factors frequently identify gene products related to cellular stresses and basal metabolism, which are common responses to growth in a hostile host environment but may not be directly involved in virulence. In order to weed out virulence-irrelevant responses, expression profiles of isogenic avirulent strains lacking essential virulence regulators [21, 22] were compared in parallel. When considered individually, any strain missing a virulence regulator may show alterations in a complex regulatory network that includes both virulence factors and genes necessary to respond to specific environmental conditions encountered by free living bacteria as well. Therefore, evaluating commonly converging effects of multiple virulence regulators and in the context of the global regulatory network during infection or, at the very least, under conditions that mimic the host environment represents an attractive novel strategy for the identification of proteins that are required for STM virulence.
In this study, we exploited multiple omics datasets (transcriptomic and proteomic) across multiple regulatory deletions under selected environmental conditions to determine a co-expression network that contains information about regulatory interactions for the infectious environment, which necessitated integration of the omics-derived data utilizing computer-aided network analysis facilitated by the context likelihood of relatedness (CLR) algorithm . Using this approach we identified previously unrecognized virulence factors as well as well-characterized virulence factors that were co-regulated by essential virulence regulators in a coordinated manner. Excitingly, the novel virulence proteins discovered through this approach were categorized into a new class of virulence effectors, which were translocated into the host cytoplasm independently of Salmonella type III secretion mechanisms.
Rationale for strains used and study design
To date, more than 300 genes have been annotated as regulatory genes in Salmonella and of those tested (83) we showed that 14 regulators including SpvR, FruR, IHF, PhoP/PhoQ, SsrA/SsrB, SlyA, Hnr, RpoE, SmpB, CsrA, RpoS, CRP, OmpR/EnvZ, and Hfq are required for virulence regulation during systemic infection in an acute mouse infection model . Salmonella strains lacking these regulators were not only attenuated for mouse virulence, but also exhibited dysregulation of many known Salmonella virulence factors, including SPI-2 genes. These regulators were linked to each other through positive or negative feedback regulation and activated SPI-2 genes in a coordinated manner . Thus we reasoned that genes that are coordinately controlled by this set of 14 virulence regulators and co-regulated with known virulence factors are likely to be important for Salmonella pathogenesis. We employed large-scale sample-matched multi-omic measurements of Salmonella strains lacking these regulators, coupled with computational network analysis to efficiently identify such genes, which represent conservatively selected putative virulence factors.
Sample-matched multi-omic profiling to infer putative virulence factors
Recent studies have discovered that a large number of Salmonella genes, estimated to be at least 20% of all genes, are post-transcriptionally regulated by a single regulatory mechanism [25, 26]. Thus mRNA and protein measurements from each regulatory mutant and wild-type strain were used as complementary methods to characterize the global/genome-wide effects for each regulator. Because some of the Salmonella regulatory mutant strains (e.g. ΔhimD, ΔrpoE, and Δhfq) survived so poorly within macrophages that preparing mRNA or protein from intracellular bacteria was not readily practical, we have used in vitro growth conditions that duplicate many aspects intracellular conditions, as previously described [24, 26, 27]. Salmonella strains were grown independently in triplicate in four in vitro conditions: Luria-Bertani (LB) logarithmic and LB stationary phases, and two acidic minimal medium (AMM) conditions, AMM1 and AMM2 [24, 26]. LB broth, given its high osmolarity and nutrient-rich condition, partially reproduces the small intestine lumen environment, while the AMM conditions, providing a low pH, low magnesium, and nutrient-deficient condition, partially mimic the intracellular milieu within the Salmonella-containing vacuole (SCV) [28, 29].
Changes in gene expression and protein abundance between the parent wild-type strain and each isogenic mutant were evaluated across the four growth conditions. Transcriptomic analysis quantified 4517 genes across 59 profiling datasets ; mass spectrometry-based label free proteomics combined with the AMT tag approach , quantified a total of 1349 proteins (Additional file 1, Table S1) across 63 profiling datasets. The large volume of data generated (over 300,000 data points from the combined omics analyses) necessitated a computational approach to integrate the data into a single systems-level view that allows efficient and systematic identification of biomolecules important for Salmonella pathogenesis.
Computational analysis of an integrated multi-omics network reveals putative virulence factors
To identify proteins important for Salmonella pathogenesis (i.e. virulence factors) in an efficient and systematic fashion we employed the context likelihood of relatedness (CLR) algorithm . The CLR algorithm can infer regulatory networks from transcriptomic data for experimentally verified regulatory interactions ; and the topology of a gene in the inferred network is correlated with phenotypic properties of the gene, for example, its role in virulence . We have previously inferred and validated a regulatory network from the transcriptomic profiles of the described regulatory mutants using CLR . In the current study we applied CLR to the proteomics data to infer protein association networks. These networks contain information about the general regulatory structure of Salmonella at the protein abundance level, but are not explicit regulatory networks. This resulted in a network of proteins and general associations between them that might indicate regulation (as in inferred transcriptional networks ), complex formation, or pathway membership.
Validation of computational predictions
SrfN, PagC, PagD, and STM2585A are translocated into the host cytosol
Virulence properties of SrfN and PagK homologues
Exploring the translcoation pathway for SrfN, PagJ, PagK1, and PagK2
In this study, we have investigated the global regulatory network required for Salmonella virulence under specific in vitro (infection-mimicking) conditions by employing regulatory protein perturbations and high throughput sample-matched omics measurements coupled with computational network analysis. This systems-level analysis of the Salmonella virulence program inferred 168 proteins which were clustered close to SPI-2 virulence proteins in the regulatory network and thus were likely to be involved in pathogenicity (Additional file 4, Table S2). A set of these predicted virulence candidate proteins were tested to verify the network prediction, including SrfN, PagD, PagC, STM1548, PdgL, STM1633, PagK2 and STM3595. These eight proteins were expressed more highly within macrophages than under LB log phase condition (Figure 2 and Additional file 7, Figure S3), which is characteristic of the SPI-2-encoded virulence proteins and suggests that their functions may be required for intracellular growth. Interestingly, SrfN, PagC, PagD, PagK2 (and its close homologues PagK1 and PagJ) were translocated to the macrophage cytoplasm (Figure 3), suggesting that these six proteins are secreted virulence effectors interacting with host cellular components to promote bacterial proliferation. In fact, SrfN and PagK homologues were required for Salmonella systemic infection in mice (Figure 4).
We showed that the integration of transcriptomics and proteomics data at the network level could provide enhanced predictions of components important for virulence. This is a novel way of integrating these two disparate data types, and should be applicable to other systems. Three of the proteins chosen for investigation had very high betweenness measures in the network (SrfN, PdgL, and STM2585A/PagK2; see Additional file 4, Table S2). It is interesting that the two translocated proteins from this later group, SrfN and PagK2, were also shown to have an effect on virulence in mice. Though a limited validation, this supports our hypothesis that the topology of the network can be informative about the importance of genes/proteins. PdgL was not translocated, but possesses a PhoP box as other PhoP-regulated virulence proteins (e.g., MgtC, PagK, VirK, PipD) responding to the intracellular environment [47, 48] and its inactivation rendered Salmonella hypervirulent in mice infection, suggesting its crucial role concerned with Salmonella intracellular fitness , a possibility we are currently examining. Our approach, combining phenotype-specific regulatory mutants with computational analysis of resultant multi-omics data, provides a novel and useful method for prediction of virulence-related genes.
Osborne et al.  recently reported that srfN, an ancestral PhoP-regulated gene, acquired an SsrB-regulatory module during Salmonella's evolution to a pathogenic bacterium. We also verified SsrB- and PhoP-dependent transcription of srfN measuring mRNA levels during growth in AMM1 medium (Additional file 12, Figure S8). However, they did not observe SrfN translocation to the host cytoplasm . The most probable explanation for this apparent discrepancy is due to the different cell types used, i.e., epithelial cells used in that study versus macrophages used in our study. Secretion of effector proteins may be specific to the intracellular environment, which has been reported for SseJ . SseJ, a SPI-2 effector was not translocated into the cytosol of HeLa cells at least at 2 h post-infection . Macrophage-biased expression was also reported for a number of SsrB-regulated SPI-2 genes .
PagK was first identified as a PhoP-activated virulence gene, where TnphoA insertion decreased Salmonella virulence significantly in mice . However, a deletion strain lacking pagK exhibited a wild-type phenotype in virulence . Although Gunn et al. identified a homologue, pagJ, in a 1.6 kb duplicated DNA region; they could not reproduce virulence attenuation with deletions of these homologue genes individually or in combination, suggesting that the original transposon insertion may have affected expression of additional genes . This virulence trait is explained by the results described here because the presence of a third related gene was not detected in that study and all three must be deleted to see the largest virulence defect. PagK (or PagK1) has a high amino acid identity with PagK2 (95%) as well as PagJ (83%). Furthermore, all of these three homologues were translocated into the host cytoplasm and negatively regulated in a similar way by SPI-2 type III secretion, suggesting a close correlation and shared properties among three PagK-homologous proteins. Supporting the possibility of a partially redundant function among the three homologues, the attenuated survival of bacteria lacking all three genes was able to be complemented only in part by PagK2 expression in trans (Figure 4).
Comparing mRNA levels in ΔssrA/ssrB and ΔphoP/phoQ strains, transcription of srfN and pag genes (pagC, pagD, and pagK1/pagK2/pagJ) was strongly dependent on PhoP, as noted previously [27, 38], and to a lesser extent on the SPI-2-encoded two-component regulator SsrA/SsrB, suggesting a positive role of the SPI-2 system in their regulation (Additional file 11, Figure S7). However, the absence of SPI-2 T3SS increased the levels of SrfN and Pag proteins inside macrophages, even though bacterial growth was restrained due to the lack of SPI-2 T3SS (Figure 5). Considering the result that the lack of SPI-2 T3SS did not increase mRNA levels of srfN and pagK1/pagK2/pagJ inside macrophages (Additional file 10, Figure S6), the increase in SrfN and Pag proteins must be a consequence of post-transcriptional regulation. The negative regulation of the translation of srfN and pagK homologues by SPI-2 T3SS raises the possibility that an activator such as a small non-coding RNA differentially regulates their translation in response to the lack of SPI-2 T3SS secretion following macrophage infection. The unusually long 5'-untranslated region of srfN (654-bp upstream of the translation start site; ) may form a complex hairpin stem-loop secondary structure that would block the ribosomal binding site and thereby interfere with translation initiation. We observed that deleting N-terminal amino acids in SrfN abolished the negative regulation of SPI-2 T3SS on srfN translation and, furthermore, Hfq overexpression increased srfN translation, but not its transcription (unpublished data). These observations support the hypothesis that a small RNA, as yet unidentified, might bind to the srfN mRNA 5'-region covering the N-terminus and alleviate an inhibitory secondary structure at the translation initiation site, sensing the absence of SPI-2 type III secretion, as is being investigated further.
The T3SS is a typical mechanism for pathogenic bacteria to deliver virulence factors into extracellular environments. However, SrfN and PagK homologues were translocated into the host cytoplasm independently from any T3SSs tested in this study. As a vehicle system to deliver bacterial components to the extracellular milieu, outer membrane vesicles (OMV) have been studied over the past few decades [51–53]. Deatherage et al. recently proposed a mechanism for OMV biogenesis wherein envelope protein interconnections modulate OMV release . OMV, observed in a variety of Gram-negative pathogens, are composed of outer membrane proteins, periplasmic proteins, lipopolysaccharide (LPS), and phospholipids and transfer bacterial DNAs and virulence factors to adjacent bacterial or host cells [51, 52, 55]. PagC was recently found to be secreted extracellularly via OMV as a major component of the outer membrane . The result that PagC was regulated and translocated in a similar manner with SrfN and PagK homologues in our study suggests the possibility of OMV-mediated transfer of SrfN and Pag proteins. In fact, PagK homologues were distributed as punctate compartments apart from intracellular bacteria inside macrophages in microscopic observation, supporting the possibility .
The approach of integrating multiple omics datasets (transcriptomic and proteomic) across multiple regulatory perturbations described in this study generated a robust association network in Salmonella highlighting regulatory interactions and provided a precise picture of putative virulence factors. The success of this strategy to discover known and novel virulence factors suggests that the approach may have applicability for other bacterial pathogens, and thereby help to innovate the host-pathogen research paradigm .
Salmonella strains and growth conditions
All deletion and tagged strains were constructed using Salmonella enterica serovar Typhimuirum 14028s as the parent strain. The phase λ Red recombination system was employed to delete or tag genes of interest as described previously [24, 59]. Linearized PCR products containing an antibiotic resistance cassette (kan or cat) and 40-nt sequences, homologous to the target sites, at both termini were introduced into recipient cells to replace genes of interest. FLP recombinase produced from pCP20 eliminated the antibiotic resistance gene via site-specific recombination  and resulted in in-frame and non-polar deletions of the target genes. Construction of mutant strains lacking regulators was described in detail previously . Translational fusions of double HA, CyaA' and Bla were constructed in a similar manner using pKD13-2HA , pMini-Tn5-cycler  and pMini-Tn5-BLAM  as PCR template plasmids. A DNA fragment encoding 2HA, CyaA' or Bla as well as Kan was introduced prior to the stop codon sequence of target gene. All primers used in strains construction are listed in Additional file 13, Table S4. Construction of chromosomal sipA::cyaA' and sseJ::cyaA' was described previously .
pSrfN and pPagK2 expressing SrfN and PagK2 respectively were constructed by cloning DNA fragments containing 900 bp upstream sequences from each start codon and coding sequences of srfN or pagK2 on pWKS30 via EcoRI and XbaI. Primers used in the PCR amplification of srfN and pagK2 are shown in Additional file 14, Table S5.
To visualize Salmonella using red fluorescent protein (Tomato) in microscopic analysis, bacteria were transformed with pWKS30-Tomato .
Bacterial cells were grown in Luria-Bertani (LB) medium or acidic minimal medium (AMM) as described previously .
Bacteria were grown under four different conditions: log phase in LB medium, stationary phase in LB medium, MgM dilution (AMM1), and MgM medium shock (AMM2) as described previously . Following RNAprotect (Qiagen) treatment to prevent RNA degradation, cells were lysed and processed to isolate total RNA as described in the instructions of Qiagen RNeasy midi kit (Qiagen) in combination with DNase I (Qiagen) treatment. Isolated total RNA was analysed by RNA 6000 nano assay (Agilent 2100 bioanalyzer, Agilent Technologies). Three biological replicates of RNAs from each strain under each growth condition were pooled with the same amount of RNA from each sample due to the number of microarray necessary for this experiment. A Δcrp strain had a growth defect in an acidic minimal medium and thereby was excluded from transcriptional profiling for AMM1 condition. Overall, a total of 59 RNA pools from 15 strains in 4 conditions were subjected to microarray analysis according to the procedures described . Strains tested included fourteen mutant strains lacking virulence regulators and one wild-type strain (S. Typhimurium 14028s). Total RNAs from 15 strains were reverse transcribed with random hexamers and Superscript II (Invitrogen) and the subsequent reverse transcripts were labeled by Cy3-linked dCTP (Amersham Biosciences). As the reference panel in all microarray analyses, genomic DNA was isolated from S. Typhimurium14028s grown in LB broth overnight using GeneElute bacterial genomic DNA kit (Sigma) and then labeled with Cy5-linked dCTP (Amersham Biosciences). A hybridization mixture containing equivalent amounts of Cy3-labled probes and Cy5-labeled probes was applied on the Salmonella array slide and incubated in a hybridization chamber (Corning) at 42°C overnight. The Salmonella-specific microarray used in this study represented PCR-amplified sequences from the annotated open reading frames (ORFs) of S. Typhimurium LT2 and was further supplemented with annotated ORFs from the serovar Typhi CT18 strain that were > 10% divergent from those of serovar Typhimurium . Each array slide has identical triplicate arrays for statistical analysis. Slides treated with Cy3/Cy5-labled probes were scanned using ScanArray Express (Packard Bioscience, BioChip Technologies) and the fluorescent signal intensities were quantified using QuantArray software (Packard Bioscience, BioChip Technologies) and exported to Excel files. The raw data was processed in WebArray  for statistical normalization. The results were displayed in a log2 scale (log2[signal intensity in a mutant/signal intensity in a wild-type]). Links to raw proteomics and transcriptomics data are available at our project website http://www.Sysbep.org.
Proteomic analysis: capillary liquid chromatography-mass spectrometry (LC-MS) analysis
Soluble Protein Preparation
Cell pellets were resuspended in 100 mM NH4HCO3, pH 8.4 buffer and lysed by using 0.1 mM zirconia/silica beads in a 2.0-ml Cryovial with vigorous vortexing for a total of 3 min with cooling steps. The supernatant and subsequent washes were transferred from the beads into new Cryovials. The beads were repeatedly washed until the supernatant was clear. After protein concentration was determined for the samples, urea and thiourea were added to final concentrations of 7 and 2 M, respectively. Following addition of DTT (5 mM), the samples were incubated at 60°C for 30 min. The samples were then diluted 10-fold with buffer, and CaCl2 was added (1 mM) followed by trypsin in a 1:50 trypsin:protein ratio. The samples were digested for 3 h at 37°C and subsequently cleaned using a C18 solid phase extraction (SPE) column (Supelco). Each 1-ml 100- or 50-mg SPE column was conditioned with MeOH and rinsed with 0.1% TFA in water. Samples were introduced to the columns and then washed with 95:5 H2O:ACN that contained 0.1% TFA. Excess liquid was removed from the columns under vacuum, and the samples were eluted with 80:20 ACN:H2O that contained 0.1% TFA and concentrated in a SpeedVac (Thermo-Savant) to a final volume of ~100 μl. A BCA protein assay was performed to determine peptide concentrations prior to analysis.
Insoluble Protein Preparation
Cell pellets were treated and lysed as described above for soluble protein preparations. The lysate was centrifuged at 1,300 × g at 4°C for 2 min, and the supernatant was transferred to polycarbonate ultracentrifuge tubes (Beckman) and centrifuged at 4°C at 356,000 × g for 10 min. Pellets were resuspended in 50 mM NH4HCO3, pH 7.8, and ultracentrifuged under the same conditions as used in the previous step. A BCA protein assay was performed on the pellets resuspended in water, and the samples were ultracentrifuged once again (as described above) before discarding the supernatant. Pellets were resuspended in ~200 μL of a solubilization solution (7 M urea, 2 M thiourea, 1% CHAPS in 50 mM ammonium bicarbonate, pH 7.8), and DTT was added to a final concentration of 9.7 mM. Samples were incubated at 60°C for 30 min, then diluted, and digested in the same manner as described for the global and soluble protein preparation. Samples were cleaned by using an appropriately sized strong cation exchange (SCX) SPE column (Supelco). Each 1-ml 100-mg column was conditioned with MeOH, rinsed in varying sequences and amounts with 10 mM ammonium formate in 25% ACN, pH 3.0; 500 mM ammonium formate in 25% ACN; and Nanopure water. Samples were acidified to a pH of ≤ 4.0 by adding 20% formic acid followed by centrifugation at 16,000 × g for 5 min. Samples were then introduced to the columns and washed with 10 mM ammonium formate in 25% ACN, pH 3.0. Excess liquid was removed from the columns under vacuum. The samples were eluted with 80:15:5 MeOH:H2O:NH4OH and concentrated to ~100 μl using a SpeedVac. Final peptide concentrations were determined using a BCA protein assay.
Capillary Liquid Chromatography (LC)-Mass Spectrometry (MS) Analysis
LC-MS spectra were analyzed using the accurate mass and elution time (AMT) tag approach . Briefly, the theoretical mass and the observed normalized elution time (NET) of each peptide identified by LC-MS/MS was used to construct a reference database of AMT tags, which served as two-dimensional markers for identifying peptides in subsequent high resolution and high mass accuracy LC-MS analyses. A reference database of AMT tags for S. Typhimurium has been generated through exhaustive SCX fractionation and LC-MS/MS analysis as described in our previous publications  and the LC-MS/MS data from this experiment, obviating the need for further time consuming SCX LC-MS/MS analyses. This approach to proteomics research was enabled by a number of published and unpublished in-group developed tools, which are available for download at http://omics.pnl.gov[62–66]. Prior to the samples being analyzed they were subjected to a blocking and randomization treatment to minimize the effects of systematic biases and ensure the even distribution of known and unknown confounding factors across the entire experimental dataset. Peptides from each of the soluble and insoluble protein preparations were separated by an automated in-house designed reverse-phase capillary HPLC system as described elsewhere . Eluate from the HPLC was directly electrosprayed into an LTQ-Orbitrap mass spectrometer (LTQ-Orbitrap, Thermo Fisher Scientific) using electrospray ionization (ESI) with emitters described previously  and the ESI interface modified with an electrodynamic ion funnel . Three biological replicates for each sample were pooled and analyzed on the ESI interface modified LTQ-Orbitrap mass spectrometer. Relevant information such as the elution time from the capillary LC column, the abundance of the signal (integrated area under the peptide peak), and the monoisotopic mass (determined from charge state and the high accuracy m/z measurement) of each feature observed in the LTQ-Orbitrap was used to match the peptide identifications contained within the AMT tag database. A matching feature was required to be within 1.6 ppm of the peptide mass and the normalized elution time (NET) within 2%. To assess false discovery rates each feature observed in the LTQ-Orbitrap was used to match the peptide identifications contained within a decoy AMT tag database. An FDR of < 1% was calculated. The abundances of these identified peptides are quantified using the area under the peptide peak and were used to infer the protein composition and relative protein abundances as described  and implemented in the software program DAnTE . Links to raw proteomics and transcriptomics data are available at our project website http://www.Sysbep.org.
Network analysis of omics data
The context likelihood of relatedness (CLR) algorithm  was utilized to systematically infer association relationships between proteins or genes. For the proteomics network we used the LC-MS-derived abundance profiles (as described above) from all the deletions under each of the four conditions with missing values filled with values corresponding to one-half the minimal observed abundance value as input to CLR. CLR calculates the mutual information between all pairs of proteins and then calculates a normalized relationship score for all resulting edges. We used 5 bins for calculation of mutual information in CLR and a spline degree of 3 for curve fitting of the binned data. For the transcriptomic network we used the transcriptional data described above as input into CLR.
To expand the integrated network to include genes not observed by proteomics, we used an approach similar to that used for topological analysis of the integrated networks (Additional file 2, Methods S1 and Additional file 5, Table S3). Each network was first filtered individually to retain only highly confident edges. The proteomics network was filtered to include only edges with Z scores above 6.0 and the transcriptomics network was filtered to include edges with Z scores above 8.0. A more stringent threshold was used with the transcriptomics to limit the influence of the much larger dataset and retain a focus on the proteomics data. Other thresholds were also applied and these resulted in very similar clustering. These networks also retained the enrichment of bottlenecks in virulence essential genes, observed in the integrated network above (data not shown). Clustering was accomplished by considering the final expanded integrated network as a similarity matrix, with a value of 0 for those gene/protein pairs that did not have an edge in the final network. The similarity matrix was used as input to the hclust function in R using the mcquitty agglomeration method .
Macrophage infection and cAMP assay
Salmonella cells were grown in Luria-Bertani (LB) broth and then opsonized with 1% mouse serum (Innovative Research) for 20 min prior to infection. RAW264.7 cells, murine macrophage-like cells (ATCC TIP-71), were infected with Salmonella at an input multiplicity of infection (MOI) of 100 by centrifuging the bacteria onto the macrophage monolayers at 1,000 × g for 5 min. The conditioned Dulbecco's Modified Eagle's Medium (DMEM) was replaced with DMEM containing gentamycin (Gibco) at 100 μg/ml to remove extracellular bacteria following 30 min incubation at 37°C with 5% CO2. After gentamycin (100 μg/ml) treatment for 1 hour, cells were washed with PBS and overlaid with DMEM containing 20 μg/ml gentamycin for the remainder of the experiments. To examine translocation of Salmonella proteins into the cytosol, macrophage cells were washed with PBS and then lysed with 0.1 M HCl as directed in the manufacturer's instructions for the cyclic AMP EIA kit (Assay Designs). Macrophage debris containing bacteria was removed by centrifugation and the cytosolic fraction was used to determine cAMP levels in the cytoplasm.
To analyze proteins from Salmonella grown in laboratory conditions, cell lysate from 5 × 107 colony-forming units (CFU) was loaded on SDS-PAGE gel. For Western blotting on intracellular Salmonella cells, infected macrophage cells were disrupted with a lysis buffer (50 mM HEPES, 1 mM EDTA, 1 mM EGTA, 1% Triton X-100, 100 mM PMSF, protease inhibitor cocktail (Roche), 50 mM NaF, DNase I (Qiagen; DNase Set), and 2 mM sodium orthovanadate) on ice for 10 min and centrifuged at 10,000 × g for 10 min to pellet intracellular bacteria. The pellet was directly resuspended in Laemmli sample buffer, boiled for 5 min, and then loaded on SDS-PAGE gel. Primary antibodies used are anti-HA.11 antibody (Covance), ant-CyaA (3D1) antibody (Santa Cruz Biotechnology), and anti-DnaK antibody (Assay Designs). As a secondary antibody, anti-mouse IgG conjugated with peroxidase (Sigma) was used in all immunoblot experiments.
CCF4-AM cleavage detection
To detect translocation of Bla fusions, macrophages were seeded in Lab-Tek II chamber coverglass slides (Nunc) and infected with Salmonella Bla fusion strains as described above. CCF4-AM (Invitrogen) solution was loaded on macrophages monolayers for 2 hours per manufacturer's instructions and the cleavage was analyzed using an Applied Precision (Issaquah) DeltaVision image restoration system with emission filter sets for green (528 nm) and blue (457 nm) fluorescence by CCF4-AM and Tomato fluorescent protein (617 nm). CCF4-AM cleavage assays were performed three times and the images shown are representatives of three independent assays.
Mouse infection experiments
In order to estimate competitive index for each test strain, Salmonella mutant/wild-type strains and the reference strain MA6054  were cultivated in LB medium overnight. Salmonella Typhimurium MA6054 is a genetically manipulated strain expressing β-galactosidase in the presence of L-arabinose. A mixture of cells containing equivalent numbers from each test strain and the reference strain was inoculated into five female 129SvJ mice of 5 weeks old at 10,000 CFU/mouse. Mice were dissected to isolate spleens at 7 days post-infection. The spleen homogenate and the inoculum were plated on LB agar medium containing 40 μg/ml X-gal (5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside) and 1 mM arabinose. The competitive index (CI) was calculated as [% of test strain recovered/% of reference strain recovered]/[% of test strain inoculated/% of reference strain inoculated].
We would like to thank George Niemann, Roslyn Brown, Liang Shi, and Penny Colton for helpful discussions. For additional technical assistance we would like to thank Ron Moore and Matt Monroe for proteomic analysis and Aurelie Snyder for immunofluorescence microscopy. This work was supported in part by the National Institute of Allergy and Infectious Diseases NIH/DHHS through interagency agreements Y1-AI-4894-01 and Y1-AI-8401-01 (project website http://www.SysBEP.org with links to raw proteomics and transcriptomics data). This work used instrumentation and capabilities developed under support from the NIH National Center for Research Resources (Grant RR 018522) and the U. S. Department of Energy Office of Biological and Environmental Research (DOE/BER). Significant portions of this work were performed using EMSL, a DOE/BER national scientific user facility located at Pacific Northwest National Laboratory. The Pacific Northwest National Laboratory is operated for the DOE by Battelle under Contract DE-AC05-76RLO1830.
- Bordbar A, Lewis NE, Schellenberger J, Palsson BO, Jamshidi N: Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions. Mol Syst Biol. 2010, 6: 422-PubMed CentralView ArticlePubMedGoogle Scholar
- Balazsi G, Heath AP, Shi L, Gennaro ML: The temporal response of the Mycobacterium tuberculosis gene regulatory network during growth arrest. Mol Syst Biol. 2008, 4: 225-PubMed CentralView ArticlePubMedGoogle Scholar
- Burstein D, Zusman T, Degtyar E, Viner R, Segal G, Pupko T: Genome-scale identification of Legionella pneumophila effectors using a machine learning approach. PLoS Pathog. 2009, 5: e1000508- 10.1371/journal.ppat.1000508PubMed CentralView ArticlePubMedGoogle Scholar
- Graham SM: Salmonellosis in children in developing and developed countries and populations. Curr Opin Infect Dis. 2002, 15: 507-512. 10.1097/00001432-200210000-00009View ArticlePubMedGoogle Scholar
- Raffatellu M, Santos RL, Verhoeven DE, George MD, Wilson RP, Winter SE, Godinez I, Sankaran S, Paixao TA, Gordon MA, Kolls JK, Dandekar S, Baumler AJ: Simian immunodeficiency virus-induced mucosal interleukin-17 deficiency promotes Salmonella dissemination from the gut. Nat Med. 2008, 14: 421-428. 10.1038/nm1743PubMed CentralView ArticlePubMedGoogle Scholar
- Santos RL, Zhang S, Tsolis RM, Kingsley RA, Adams LG, Baumler AJ: Animal models of Salmonella infections: enteritis versus typhoid fever. Microbes Infect. 2001, 3: 1335-1344. 10.1016/S1286-4579(01)01495-2View ArticlePubMedGoogle Scholar
- Tsolis RM, Kingsley RA, Townsend SM, Ficht TA, Adams LG, Baumler AJ: Of mice, calves, and men. Comparison of the mouse typhoid model with other Salmonella infections. Adv Exp Med Biol. 1999, 473: 261-274.View ArticlePubMedGoogle Scholar
- Geddes K, Cruz F, Heffron F: Analysis of cells targeted by Salmonella type III secretion in vivo. PLoS Pathog. 2007, 3: e196- 10.1371/journal.ppat.0030196PubMed CentralView ArticlePubMedGoogle Scholar
- Meresse S, Steele-Mortimer O, Moreno E, Desjardins M, Finlay B, Gorvel JP: Controlling the maturation of pathogen-containing vacuoles: a matter of life and death. Nat Cell Biol. 1999, 1: E183-8. 10.1038/15620View ArticlePubMedGoogle Scholar
- Holden DW: Trafficking of the Salmonella vacuole in macrophages. Traffic. 2002, 3: 161-169. 10.1034/j.1600-0854.2002.030301.xView ArticlePubMedGoogle Scholar
- Knodler LA, Steele-Mortimer O: Taking possession: biogenesis of the Salmonella-containing vacuole. Traffic. 2003, 4: 587-599. 10.1034/j.1600-0854.2003.00118.xView ArticlePubMedGoogle Scholar
- Bang IS, Audia JP, Park YK, Foster JW: Autoinduction of the ompR response regulator by acid shock and control of the Salmonella enterica acid tolerance response. Mol Microbiol. 2002, 44: 1235-1250. 10.1046/j.1365-2958.2002.02937.xView ArticlePubMedGoogle Scholar
- Shi Y, Cromie MJ, Hsu FF, Turk J, Groisman EA: PhoP-regulated Salmonella resistance to the antimicrobial peptides magainin 2 and polymyxin B. Mol Microbiol. 2004, 53: 229-241. 10.1111/j.1365-2958.2004.04107.xView ArticlePubMedGoogle Scholar
- Buchmeier N, Bossie S, Chen CY, Fang FC, Guiney DG, Libby SJ: SlyA, a transcriptional regulator of Salmonella typhimurium, is required for resistance to oxidative stress and is expressed in the intracellular environment of macrophages. Infect Immun. 1997, 65: 3725-3730.PubMed CentralPubMedGoogle Scholar
- Hensel M: Salmonella pathogenicity island 2. Mol Microbiol. 2000, 36: 1015-1023. 10.1046/j.1365-2958.2000.01935.xView ArticlePubMedGoogle Scholar
- Adkins JN, Mottaz HM, Norbeck AD, Gustin JK, Rue J, Clauss TR, Purvine SO, Rodland KD, Heffron F, Smith RD: Analysis of the Salmonella typhimurium proteome through environmental response toward infectious conditions. Mol Cell Proteomics. 2006, 5: 1450-1461. 10.1074/mcp.M600139-MCP200View ArticlePubMedGoogle Scholar
- Ansong C, Yoon H, Norbeck AD, Gustin JK, McDermott JE, Mottaz HM, Rue J, Adkins JN, Heffron F, Smith RD: Proteomics analysis of the causative agent of typhoid fever. J Proteome Res. 2008, 7: 546-557. 10.1021/pr070434uView ArticlePubMedGoogle Scholar
- Eriksson S, Lucchini S, Thompson A, Rhen M, Hinton JC: Unravelling the biology of macrophage infection by gene expression profiling of intracellular Salmonella enterica. Mol Microbiol. 2003, 47: 103-118.View ArticlePubMedGoogle Scholar
- Faucher SP, Porwollik S, Dozois CM, McClelland M, Daigle F: Transcriptome of Salmonella enterica serovar Typhi within macrophages revealed through the selective capture of transcribed sequences. Proc Natl Acad Sci USA. 2006, 103: 1906-1911. 10.1073/pnas.0509183103PubMed CentralView ArticlePubMedGoogle Scholar
- Shi L, Adkins JN, Coleman JR, Schepmoes AA, Dohnkova A, Mottaz HM, Norbeck AD, Purvine SO, Manes NP, Smallwood HS, Wang H, Forbes J, Gros P, Uzzau S, Rodland KD, Heffron F, Smith RD, Squier TC: Proteomic analysis of Salmonella enterica serovar typhimurium isolated from RAW 264.7 macrophages: identification of a novel protein that contributes to the replication of serovar typhimurium inside macrophages. J Biol Chem. 2006, 281: 29131-29140. 10.1074/jbc.M604640200View ArticlePubMedGoogle Scholar
- Detweiler CS, Cunanan DB, Falkow S: Host microarray analysis reveals a role for the Salmonella response regulator phoP in human macrophage cell death. Proc Natl Acad Sci USA. 2001, 98: 5850-5855. 10.1073/pnas.091110098PubMed CentralView ArticlePubMedGoogle Scholar
- Coombes BK, Wickham ME, Lowden MJ, Brown NF, Finlay BB: Negative regulation of Salmonella pathogenicity island 2 is required for contextual control of virulence during typhoid. Proc Natl Acad Sci USA. 2005, 102: 17460-17465. 10.1073/pnas.0505401102PubMed CentralView ArticlePubMedGoogle Scholar
- Faith JJ, Hayete B, Thaden JT, Mogno I, Wierzbowski J, Cottarel G, Kasif S, Collins JJ, Gardner TS: Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles. PLoS Biol. 2007, 5: e8- 10.1371/journal.pbio.0050008PubMed CentralView ArticlePubMedGoogle Scholar
- Yoon H, McDermott JE, Porwollik S, McClelland M, Heffron F: Coordinated regulation of virulence during systemic infection of Salmonella enterica serovar Typhimurium. PLoS Pathog. 2009, 5: e1000306- 10.1371/journal.ppat.1000306PubMed CentralView ArticlePubMedGoogle Scholar
- Sittka A, Lucchini S, Papenfort K, Sharma CM, Rolle K, Binnewies TT, Hinton JC, Vogel J: Deep sequencing analysis of small noncoding RNA and mRNA targets of the global post-transcriptional regulator, Hfq. PLoS Genet. 2008, 4: e1000163- 10.1371/journal.pgen.1000163PubMed CentralView ArticlePubMedGoogle Scholar
- Ansong C, Yoon H, Porwollik S, Mottaz-Brewer H, Petritis BO, Jaitly N, Adkins JN, McClelland M, Heffron F, Smith RD: Global systems-level analysis of Hfq and SmpB deletion mutants in Salmonella: implications for virulence and global protein translation. PLoS One. 2009, 4: e4809- 10.1371/journal.pone.0004809PubMed CentralView ArticlePubMedGoogle Scholar
- Osborne SE, Walthers D, Tomljenovic AM, Mulder DT, Silphaduang U, Duong N, Lowden MJ, Wickham ME, Waller RF, Kenney LJ, Coombes BK: Pathogenic adaptation of intracellular bacteria by rewiring a cis-regulatory input function. Proc Natl Acad Sci USA. 2009, 106: 3982-3987. 10.1073/pnas.0811669106PubMed CentralView ArticlePubMedGoogle Scholar
- Deiwick J, Nikolaus T, Erdogan S, Hensel M: Environmental regulation of Salmonella pathogenicity island 2 gene expression. Mol Microbiol. 1999, 31: 1759-1773. 10.1046/j.1365-2958.1999.01312.xView ArticlePubMedGoogle Scholar
- Miao EA, Miller SI: A conserved amino acid sequence directing intracellular type III secretion by Salmonella typhimurium. Proc Natl Acad Sci USA. 2000, 97: 7539-7544. 10.1073/pnas.97.13.7539PubMed CentralView ArticlePubMedGoogle Scholar
- Zimmer JS, Monroe ME, Qian WJ, Smith RD: Advances in proteomics data analysis and display using an accurate mass and time tag approach. Mass Spectrom Rev. 2006, 25: 450-482. 10.1002/mas.20071PubMed CentralView ArticlePubMedGoogle Scholar
- McDermott JE, Taylor RC, Yoon H, Heffron F: Bottlenecks and hubs in inferred networks are important for virulence in Salmonella typhimurium. J Comput Biol. 2009, 16: 169-180. 10.1089/cmb.2008.04TTView ArticlePubMedGoogle Scholar
- Rappl C, Deiwick J, Hensel M: Acidic pH is required for the functional assembly of the type III secretion system encoded by Salmonella pathogenicity island 2. FEMS Microbiol Lett. 2003, 226: 363-372. 10.1016/S0378-1097(03)00638-4View ArticlePubMedGoogle Scholar
- McGhie EJ, Brawn LC, Hume PJ, Humphreys D, Koronakis V: Salmonella takes control: effector-driven manipulation of the host. Curr Opin Microbiol. 2009, 12: 117-124. 10.1016/j.mib.2008.12.001PubMed CentralView ArticlePubMedGoogle Scholar
- Niemann GS, Brown RN, Gustin JK, Stufkens A, Shaikh-Kidwai AS, Li J, McDermott JE, Brewer HM, Schepmoes A, Smith RD, Adkins JN, Heffron F: Discovery of novel secreted virulence factors from Salmonella enterica serovar Typhimurium by proteomic analysis of culture supernatants. Infect Immun. 2011, 79: 33-43. 10.1128/IAI.00771-10PubMed CentralView ArticlePubMedGoogle Scholar
- Sory MP, Boland A, Lambermont I, Cornelis GR: Identification of the YopE and YopH domains required for secretion and internalization into the cytosol of macrophages, using the cyaA gene fusion approach. Proc Natl Acad Sci USA. 1995, 92: 11998-12002. 10.1073/pnas.92.26.11998PubMed CentralView ArticlePubMedGoogle Scholar
- Ohlson MB, Fluhr K, Birmingham CL, Brumell JH, Miller SI: SseJ deacylase activity by Salmonella enterica serovar Typhimurium promotes virulence in mice. Infect Immun. 2005, 73: 6249-6259. 10.1128/IAI.73.10.6249-6259.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Belden WJ, Miller SI: Further characterization of the PhoP regulon: identification of new PhoP-activated virulence loci. Infect Immun. 1994, 62: 5095-5101.PubMed CentralPubMedGoogle Scholar
- Gunn JS, Alpuche-Aranda CM, Loomis WP, Belden WJ, Miller SI: Characterization of the Salmonella typhimurium pagC/pagD chromosomal region. J Bacteriol. 1995, 177: 5040-5047.PubMed CentralPubMedGoogle Scholar
- Gunn JS, Belden WJ, Miller SI: Identification of PhoP-PhoQ activated genes within a duplicated region of the Salmonella typhimurium chromosome. Microb Pathog. 1998, 25: 77-90. 10.1006/mpat.1998.0217View ArticlePubMedGoogle Scholar
- Miller SI, Kukral AM, Mekalanos JJ: A two-component regulatory system (phoP phoQ) controls Salmonella typhimurium virulence. Proc Natl Acad Sci USA. 1989, 86: 5054-5058. 10.1073/pnas.86.13.5054PubMed CentralView ArticlePubMedGoogle Scholar
- Geddes K, Worley M, Niemann G, Heffron F: Identification of new secreted effectors in Salmonella enterica serovar Typhimurium. Infect Immun. 2005, 73: 6260-6271. 10.1128/IAI.73.10.6260-6271.2005PubMed CentralView ArticlePubMedGoogle Scholar
- Poh J, Odendall C, Spanos A, Boyle C, Liu M, Freemont P, Holden DW: SteC is a Salmonella kinase required for SPI-2-dependent F-actin remodelling. Cell Microbiol. 2008, 10: 20-30.PubMed CentralPubMedGoogle Scholar
- Galan JE: Salmonella interactions with host cells: type III secretion at work. Annu Rev Cell Dev Biol. 2001, 17: 53-86. 10.1146/annurev.cellbio.17.1.53View ArticlePubMedGoogle Scholar
- Hansen-Wester I, Hensel M: Salmonella pathogenicity islands encoding type III secretion systems. Microbes Infect. 2001, 3: 549-559. 10.1016/S1286-4579(01)01411-3View ArticlePubMedGoogle Scholar
- Ghelardi E, Celandroni F, Salvetti S, Beecher DJ, Gominet M, Lereclus D, Wong AC, Senesi S: Requirement of flhA for swarming differentiation, flagellin export, and secretion of virulence-associated proteins in Bacillus thuringiensis. J Bacteriol. 2002, 184: 6424-6433. 10.1128/JB.184.23.6424-6433.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Young GM, Schmiel DH, Miller VL: A new pathway for the secretion of virulence factors by bacteria: the flagellar export apparatus functions as a protein-secretion system. Proc Natl Acad Sci USA. 1999, 96: 6456-6461. 10.1073/pnas.96.11.6456PubMed CentralView ArticlePubMedGoogle Scholar
- Lejona S, Aguirre A, Cabeza ML, Garcia Vescovi E, Soncini FC: Molecular characterization of the Mg2+-responsive PhoP-PhoQ regulon in Salmonella enterica. J Bacteriol. 2003, 185: 6287-6294. 10.1128/JB.185.21.6287-6294.2003PubMed CentralView ArticlePubMedGoogle Scholar
- Yamamoto K, Ogasawara H, Fujita N, Utsumi R, Ishihama A: Novel mode of transcription regulation of divergently overlapping promoters by PhoP, the regulator of two-component system sensing external magnesium availability. Mol Microbiol. 2002, 45: 423-438. 10.1046/j.1365-2958.2002.03017.xView ArticlePubMedGoogle Scholar
- Mouslim C, Hilbert F, Huang H, Groisman EA: Conflicting needs for a Salmonella hypervirulence gene in host and non-host environments. Mol Microbiol. 2002, 45: 1019-1027. 10.1046/j.1365-2958.2002.03070.xView ArticlePubMedGoogle Scholar
- Worley MJ, Ching KH, Heffron F: Salmonella SsrB activates a global regulon of horizontally acquired genes. Mol Microbiol. 2000, 36: 749-761.View ArticlePubMedGoogle Scholar
- Bomberger JM, Maceachran DP, Coutermarsh BA, Ye S, O'Toole GA, Stanton BA: Long-distance delivery of bacterial virulence factors by Pseudomonas aeruginosa outer membrane vesicles. PLoS Pathog. 2009, 5: e1000382- 10.1371/journal.ppat.1000382PubMed CentralView ArticlePubMedGoogle Scholar
- Kesty NC, Mason KM, Reedy M, Miller SE, Kuehn MJ: Enterotoxigenic Escherichia coli vesicles target toxin delivery into mammalian cells. EMBO J. 2004, 23: 4538-4549. 10.1038/sj.emboj.7600471PubMed CentralView ArticlePubMedGoogle Scholar
- Wai SN, Lindmark B, Soderblom T, Takade A, Westermark M, Oscarsson J, Jass J, Richter-Dahlfors A, Mizunoe Y, Uhlin BE: Vesicle-mediated export and assembly of pore-forming oligomers of the enterobacterial ClyA cytotoxin. Cell. 2003, 115: 25-35. 10.1016/S0092-8674(03)00754-2View ArticlePubMedGoogle Scholar
- Deatherage BL, Lara JC, Bergsbaken T, Rassoulian Barrett SL, Lara S, Cookson BT: Biogenesis of bacterial membrane vesicles. Mol Microbiol. 2009, 72: 1395-1407. 10.1111/j.1365-2958.2009.06731.xPubMed CentralView ArticlePubMedGoogle Scholar
- Kolling GL, Matthews KR: Export of virulence genes and Shiga toxin by membrane vesicles of Escherichia coli O157:H7. Appl Environ Microbiol. 1999, 65: 1843-1848.PubMed CentralPubMedGoogle Scholar
- Kitagawa R, Takaya A, Ohya M, Mizunoe Y, Takade A, Yoshida S, Isogai E, Yamamoto T: Biogenesis of Salmonella enterica serovar typhimurium membrane vesicles provoked by induction of PagC. J Bacteriol. 2010, 192: 5645-5656. 10.1128/JB.00590-10PubMed CentralView ArticlePubMedGoogle Scholar
- Yoon H, Ansong C, Adkins JN, Heffron F: Discovery of Salmonella Virulence Factors Translocated via Outer Membrane Vesicles (OMV) to Murine Macrophages. Infect Immun. 2011,Google Scholar
- Aderem A, Adkins JN, Ansong C, Galagan J, Kaiser S, Korth MJ, Law GL, McDermott JG, Proll SC, Rosenberger C, Schoolnik G, Katze MG: A systems biology approach to infectious disease research: innovating the pathogen-host research paradigm. MBio. 2011, 2:Google Scholar
- Datsenko KA, Wanner BL: One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci USA. 2000, 97: 6640-6645. 10.1073/pnas.120163297PubMed CentralView ArticlePubMedGoogle Scholar
- Porwollik S, Frye J, Florea LD, Blackmer F, McClelland M: A non-redundant microarray of genes for two related bacteria. Nucleic Acids Res. 2003, 31: 1869-1876. 10.1093/nar/gkg298PubMed CentralView ArticlePubMedGoogle Scholar
- Xia X, McClelland M, Wang Y: WebArray: an online platform for microarray data analysis. BMC Bioinformatics. 2005, 6: 306- 10.1186/1471-2105-6-306PubMed CentralView ArticlePubMedGoogle Scholar
- Jaitly N, Monroe ME, Petyuk VA, Clauss TR, Adkins JN, Smith RD: Robust algorithm for alignment of liquid chromatography-mass spectrometry analyses in an accurate mass and time tag data analysis pipeline. Anal Chem. 2006, 78: 7397-7409. 10.1021/ac052197pView ArticlePubMedGoogle Scholar
- Kiebel GR, Auberry KJ, Jaitly N, Clark DA, Monroe ME, Peterson ES, Tolic N, Anderson GA, Smith RD: PRISM: a data management system for high-throughput proteomics. Proteomics. 2006, 6: 1783-1790. 10.1002/pmic.200500500View ArticlePubMedGoogle Scholar
- Monroe ME, Shaw JL, Daly DS, Adkins JN, Smith RD: MASIC: a software program for fast quantitation and flexible visualization of chromatographic profiles from detected LC-MS(/MS) features. Comput Biol Chem. 2008, 32: 215-217. 10.1016/j.compbiolchem.2008.02.006PubMed CentralView ArticlePubMedGoogle Scholar
- Monroe ME, Tolic N, Jaitly N, Shaw JL, Adkins JN, Smith RD: VIPER: an advanced software package to support high-throughput LC-MS peptide identification. Bioinformatics. 2007, 23: 2021-2023. 10.1093/bioinformatics/btm281View ArticlePubMedGoogle Scholar
- Petritis K, Kangas LJ, Yan B, Monroe ME, Strittmatter EF, Qian WJ, Adkins JN, Moore RJ, Xu Y, Lipton MS, Camp DGn, Smith RD: Improved peptide elution time prediction for reversed-phase liquid chromatography-MS by incorporating peptide sequence information. Anal Chem. 2006, 78: 5026-5039. 10.1021/ac060143pPubMed CentralView ArticlePubMedGoogle Scholar
- Livesay EA, Tang K, Taylor BK, Buschbach MA, Hopkins DF, LaMarche BL, Zhao R, Shen Y, Orton DJ, Moore RJ, Kelly RT, Udseth HR, Smith RD: Fully automated four-column capillary LC-MS system for maximizing throughput in proteomic analyses. Anal Chem. 2008, 80: 294-302. 10.1021/ac701727rPubMed CentralView ArticlePubMedGoogle Scholar
- Kelly RT, Page JS, Tang K, Smith RD: Array of chemically etched fused-silica emitters for improving the sensitivity and quantitation of electrospray ionization mass spectrometry. Anal Chem. 2007, 79: 4192-4198. 10.1021/ac062417eView ArticlePubMedGoogle Scholar
- Page JS, Tolmachev AV, Tang K, Smith RD: Theoretical and experimental evaluation of the low m/z transmission of an electrodynamic ion funnel. J Am Soc Mass Spectrom. 2006, 17: 586-592. 10.1016/j.jasms.2005.12.013PubMed CentralView ArticlePubMedGoogle Scholar
- Hixson KK, Adkins JN, Baker SE, Moore RJ, Chromy BA, Smith RD, McCutchen-Maloney SL, Lipton MS: Biomarker candidate identification in Yersinia pestis using organism-wide semiquantitative proteomics. J Proteome Res. 2006, 5: 3008-3017. 10.1021/pr060179yView ArticlePubMedGoogle Scholar
- Polpitiya AD, Qian WJ, Jaitly N, Petyuk VA, Adkins JN, Camp DGn, Anderson GA, Smith RD: DAnTE: a statistical tool for quantitative analysis of -omics data. Bioinformatics. 2008, 24: 1556-1558. 10.1093/bioinformatics/btn217PubMed CentralView ArticlePubMedGoogle Scholar
- McQuitty LL: Similarity analysis by reciprocal pairs for discrete and continuous data. Educ Psychol Meas. 1966, 26: 825-831. 10.1177/001316446602600402.View ArticleGoogle Scholar
- Ho TD, Figueroa-Bossi N, Wang M, Uzzau S, Bossi L, Slauch JM: Identification of GtgE, a novel virulence factor encoded on the Gifsy-2 bacteriophage of Salmonella enterica serovar Typhimurium. J Bacteriol. 2002, 184: 5234-5239. 10.1128/JB.184.19.5234-5239.2002PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.