Skip to main content

PHIDIAS: a pathogen-host interaction data integration and analysis system

Abstract

The Pathogen-Host Interaction Data Integration and Analysis System (PHIDIAS) is a web-based database system that serves as a centralized source to search, compare, and analyze integrated genome sequences, conserved domains, and gene expression data related to pathogen-host interactions (PHIs) for pathogen species designated as high priority agents for public health and biological security. In addition, PHIDIAS allows submission, search and analysis of PHI genes and molecular networks curated from peer-reviewed literature. PHIDIAS is publicly available at http://www.phidias.us.

Rationale

An infectious disease is the result of an interactive relationship between a pathogen and its host. According to estimations of the World Health Organization, infectious diseases caused 14.7 million deaths in 2001, accounting for 26% of the total global mortality [1]. Integration and analysis of various data related to pathogens and pathogen-host interactions (PHIs) will yield a better understanding of, and means for, the control of infectious diseases induced by such pathogens.

Completely sequenced genomic information provides valuable information for gene and protein functions, and intra-organismic processes. Pathogen genome information also lays a foundation for the study of the interactions between host and microbial organisms. Several genome data resources, such as the National Center for Biotechnology Information (NCBI), European Bioinformatics Institute (EBI) and Swiss Institute of Bioinformatics (SIB), are available to the public. However, data obtained from these sources often are not integrated. Lack of such integration prompted us to develop the Brucella Bioinformatics Portal (BBP) [2]. This program allows integration of data from more than 20 sources including information on the Brucella genome. The same strategy can be expanded to include other pathogens, thereby enhancing our ability to conduct comparative studies. The program can be modified to include additional features not yet available in BBP. For example, protein conserved domains (distinct units of molecular evolution usually associated with particular molecular functions) could be listed. The NCBI Conserved Domain Database (CDD) mirrors several collections, including the Protein families database of alignments (Pfam) [3], Simple Modular Architecture Research Tool (SMART) [4], and Clusters of Orthologous Groups (COG) [5], and thus provides comprehensive information about conserved protein domains. Conserved domains are critical for protein functions and provide important clues about microbial pathogenesis and interactions between pathogens and hosts.

While CDD contains conserved domains derived from various eukaryotic and prokaryotic organisms [6], it is difficult to compare and analyze pathogen-specific conserved domains. The availability of a program that permits the acquisition and storage of pathogen-specific domain information in an integrated system would be extremely useful, as would the combination of such a database with BLAST search programs and other programs for the determination of sequence analyses. To facilitate comparison and better understanding of pathogens and fundamental PHI mechanisms, it is necessary to integrate genome information from publicly important pathogens with effective tools for browsing, searching, and analyzing annotated genome sequences and conserved domains. Such an integrated system would also benefit from the inclusion of large amounts of published literature data relating to pathogens and their interactions with host immune systems. To allow machine-readable data exchange of the now voluminous pathogen information, He et al. [7] developed an Extensible Markup Language (XML)-based Pathogen Information Markup Language (PIML). PIML contains comprehensive pathogen-oriented information, including pathogen taxonomy, genomic information, life cycle, epidemiology, induced diseases in host, diagnosis, treatment, and relevant laboratory analysis. A list of PIML documents addressing pathogens deemed of high priority for public health and biological defense have been created and are available on the worldwide web or through a web service [7]. However, compared to relational databases, XML databases do not efficiently support query functions and scalability. These deficiencies prompted us to design a web-based relational database system to store and query PIML data. The database system can also integrate efficiently other PHI-related data, including manually curated information related to the pathobiology and management of laboratory animals that are given high priority pathogens [8].

The molecular functions of pathogen and host genes as well as their roles in specific PHI pathways have been extensively studied. Molecules that play important roles in the virulence of pathogens and in the host immune defense are particularly important for PHI. A systematic collation from the literature of these molecules and their functions is lacking. Once PHI-related molecules are collated, the next step is to illustrate molecular interactions and pathways involving these molecules. Existing pathway databases, such as the Kyoto Encyclopedia of Genes and Genomes (KEGG) [9], BioCyc [10, 11], and Biomolecular Interaction Network Database (BIND) [12], contain pathways for various metabolic and molecular interactions of different organisms. Although richly documented, the networks of microbial and host molecular and cellular interactions that occur during pathogenic infections of hosts are underrepresented in current database systems. He and colleagues [13] developed the Molecular Interaction Network Markup Language (MINetML, previously called ProNetML) to summarize information related to microbial pathogenesis. However, MINetML cannot be exchanged with other standard data exchange formats such as the Biological Pathways Exchange format (BioPAX) [14]. This deficiency prevents active data exchange and communication with biological pathway databases. In addition, there is no effective MINetML visualization tool available.

Experimental methodologies, including microarrays and mass spectrometry, provide abundant sources of gene expression data. Publicly available gene expression data repositories, including the NCBI Gene Expression Omnibus (GEO) [15] and the EBI ArrayExpress [16] store large amounts of gene expression data, much of which is related to interactions between pathogens and hosts. Summaries of gene expression experiments and gene profiles allow querying and comparison of PHI-related gene expression patterns.

To better understand the intricate interactions between pathogens and hosts, we have now developed a web-based PHI data integration and analysis system (PHIDIAS) that permits integration and analysis of genome sequences, curated literature data for general PHI information and PHI networks, and PHI-related gene expression data. PHIDIAS currently targets 42 pathogens. These include most category A, B, and C priority pathogens identified by the National Institute of Allergy and Infectious Diseases (NIAID) and the Centers for Disease Control and Prevention (CDC) in the USA, and other pathogens deemed of high priority with regards to public health, such as the human immunodeficiency virus (HIV) and Plasmodium falciparum (Table 1).

Table 1 Forty-two pathogens included in PHIDIAS

System design

PHIDIAS is implemented using a three-tier architecture built on two Dell Poweredge 2580 servers that run the Redhat Linux operating system (Redhat Enterprise Linux ES 4). Users can submit database or analysis queries through the web. These queries are then processed using PHP/Perl/SQL (middle-tier, application server based on Apache) against a MySQL (version 5.0) relational database (back-end, database server). The result of each query is then presented to the user in the web browser. Two servers are scheduled to regularly backup each others' data.

PHIDIAS includes six components that search and analyze annotated genome sequences, curated PHI data, and PHI-related gene expression data (Figure 1a). Pathogen genomes are displayed and analyzed by PGBrowser, Pacodom, and BLAST searches. The PGBrowser has been developed to browse and analyze the gene and protein sequences of 77 genomes from 42 bacterial, viral, and parasitic pathogens (Table 1). Although PHDIAS does not include non-pathogenic species, PHIDIAS includes genomes from both pathogenic strains (for example, Escherichia coli O157:H7 strain Sakai) and non-pathogenic strains (for example, E. coli strain K12) in the same pathogen species. Pacodom is used to search and analyze conserved protein domains of the pathogen genomes. Customized BLAST programs allow users to perform similarity searches on pathogen genome sequences. Curated PHI data are separated into Phinfo, Phigen and Phinet, based on general PHI information, PHI molecules and networks, respectively. PHI gene expression experiments and gene profiles are searched through the Phix database system.

Figure 1
figure 1

PHIDIAS data flow. (a) The PHIDIAS system architecture. (b) PhiDB data flow among key elements of different PhiDB database modules. The relationships among these elements are represented by the following signs: *, zero or more; 1, one; and 2...*, two or more. For example, the labeling of a pathway with '1' and '2...*' indicates that one pathway includes two or more interactions.

PhiDB is the PHIDIAS relational database that integrates different PHIDIAS components. Figure 1b illustrates the relationship and data flow among different database modules and PHIDIAS components. PhiDB integrates PHI-related data from more than 20 public databases (Table 2) and from data curated by the PHIDIAS curation team. PhiDB contains gene information, including sequences, conserved domains from pathogen genomes as well as gene information for PHI and diagnosis of pathogen infections. The biological objects (Bio Object) in the data flow diagram are flexible, that is, they can be a gene or gene product, or any other molecular or cellular entity, including metabolites, cell membrane, mitochondria and so on. The Bio Object element also enables representation of a cluster or group of molecules such as virulent factors and protective antigens. Each interaction includes two or more Bio Objects that function as input or output objects. Each pathway contains more than one interaction. General information pertaining to each pathogenic organism and each disease is available and integrates with pathway and gene information. PHI-related gene expression experiments are also recorded. Detailed information for references, including peer-reviewed journal publications, reliable websites and databases for each of the components is also stored. Each of the PHIDIAS components focuses on different PhiDB elements. All of these components are integrated together and readily available for biomedical researchers working on different pathogens and PHI systems.

Table 2 Public databases and software programs integrated in PHIDIAS

To illustrate the features of data integration and comparative analyses using PHIDIAS, the pathogenic Brucella serves as an example and demonstrates how PHIDIAS can promote Brucella research. Brucella species are Gram-negative, facultative intracellular bacteria that cause brucellosis in humans and animals [17]. B. melitensis, B. suis, B. abortus, and B. canis are human pathogens in decreasing order of severity. Brucella species have been identified as priority agents amenable for use in biological warfare and bioterrorism and are listed as USA NIAID category B priority pathogens. The genomes of B. melitensis strain 16 M [18], B. suis strain 1330 [19], and B. abortus strain 994-1 [20] and strain 2308 [21] have been sequenced and published.

PHIDIAS components

PGBrowser: pathogen genome browser

Pathogen genomes serve as the foundation for the study of PHI in the post-genomic era. PGBrowser integrates data from more than 20 different sources, including NCBI, EBI, and The Institute for Genomic Research (TIGR) (Table 2). Currently, PGBrowser stores 77 genome sequences and 203,297 features from 42 pathogens. NCBI Entrez Programming Utilities are used to download genome information for the pathogens selected from Reference Sequences (RefSeq) and other NCBI databases. The information obtained is formatted in XML. A script has been developed to parse all the protein/gene features, including raw sequences. These are stored in the PhiDB database. Another script has also been developed to query UniProt and other EBI databases, and to download all of the protein information that relates to the 42 pathogens using the SwissProt format. The information is then parsed and stored in a database based on Locus Tag matches. The molecular weights and isoelectric points (pI) are calculated from the protein sequences using the modules (Bio::Tools::pICalculator and Bio::Tools::SeqStats) from BioPerl [22]. In order to enhance the query process, all pathogen sequences and annotation information for PGBrowser are stored in the database server instead of flat files.

The genome browser web interface of PGBrowser was developed based on the Generic Genome Browser (GBrowse) available at the Generic Software Components for Model Organism Databases (GMOD), a popular genome browser tool because of its portability, simple installation, convenient data input and easy integration with other software programs [23]. The GBrowse program has been used to display genome information about the bacterial pathogens Brucella spp. [2] and Pseudomonas aeruginosa [24]. PGBrowser modifies GBrowse and allows simultaneous query and analysis for any bacterial or viral gene across all 77 genomes of the 42 pathogens. For example, a query for sodC in PGBrowser results in 32 sodC hits from 32 genomes in 11 bacterial species, among which are four Brucella sodC genes from four Brucella genomes (Figure 2a). One can query any Brucella gene (for example, sodC) among the different Brucella genomes, analyze the gene sequences before and after a particular gene (Figure 2b), and obtain gene DNA, RNA, and protein sequences, and perform sequence analyses (for example, finding restriction enzyme digestion sites). As a feature inherited from GBrowse, PGBrowser also provides means for annotating restriction sites, finding short oligonucleotides, and downloading protein or DNA sequence files. PGBrowser can also be directly accessed from other PHIDIAS components such as Pacodom.

Figure 2
figure 2

Comparison and analyses of sodC genes in the PGBrowser. Thirty two sodC genes are found in 32 genomes from 11 bacteria species (a), including sodC from B. abortus strain 9-941 (b).

A detailed page of pathogen gene information has been developed to summarize integrative information about a specific pathogen gene, such as sodC in B. melitensis strain 16 M (Figure 3). It not only provides web links to various databases but also lists detailed protein annotation from authorized databases (for example, UniProt). Additionally, this page includes PHI specific information curated internally by the PHIDIAS curation team. A curator is also prompted to provide additional information using an online submission system. This page also provides DNA and protein sequences in FASTA format. The sequences can be directly linked to a customized BLAST search to find similar sequences from other pathogens. The references for curated PHI information are listed. A PubMed link is available for searching more related peer-reviewed articles. Figure 3 shows that Cu/Zn superoxide dismutase (SOD) encoded by the B. abortus sodC gene is required for Brucella protection from endogenous superoxide stress [25]. The B. abortus sodC mutant is attenuated in macrophages and mice [25]. Figure 3 also indicates that Brucella Cu/Zn SOD induces protective Th1 type immune responses and has been used for Brucella vaccine development [26]. For comparative purposes, one may examine sodC genes from other bacterial pathogens, such as Bacillus anthracis. Passalacqua et al. [27] recently showed that B. anthracis Cu/Zn SOD plays only a trivial role in protecting against endogenous superoxide stress. This indicates that the same gene may have different roles in microbial pathogenesis, suggesting that it is important to analyze pathogen genes individually, particularly in terms of the interactions between pathogens and hosts.

Figure 3
figure 3

Integrative pathogen gene information in PHIDIAS.

While PHIDIAS is pathogen-oriented and focuses on functional analysis of pathogen genes during PHI, host genome sequences may be requested for gene level PHI analyses. Since GBrowse-based human and mouse genome browsers are publicly available, PGBrowser contains a web interface that allows users to conveniently search the host genome sequence browsers by linking them to the websites.

Pacodom: pathogen protein conserved domains

The conserved domain data from completely sequenced pathogenic organisms provide valuable information for the identification of protein functions and for the study of PHI. Currently, the NCBI CDD database contains 12,589 position-specific score matrix (PSSM) models that are commonly used representations of motifs present in biological sequences. However, the PSSM models cover a broad range of organisms and, therefore, it is difficult to compare conserved domains from select priority pathogens. To circumvent this problem, a pathogen-specific protein conserved domains database module called Pacodom was developed. This program contains all possible conserved domains found in the 77 pathogen genomes of 42 pathogens. To build this system, a local reverse-position-specific (RPS) CDD library was constructed based on the CDD conserved domain data downloaded from NCBI [28]. The RPS BLAST program (downloaded from the NCBI toolkit distribution) [29] was run for each protein sequence against the RPS CDD library with an expectation value of 10-6. The domain alignments obtained from the RPS BLAST search are used to calculate the PSSM. A Perl script was developed to store non-redundant PSSM models [30] in the Pacodom MySQL database module. Currently, the Pacodom database contains 7,919 PSSMs found in 151,787 protein sequences. This value comprises 76.4% of a total of 198,696 proteins from all genomes available in PhiDB.

The conserved domain data from completely sequenced pathogenic organisms provide valuable information for comparative analysis of functional roles of pathogen proteins and their involvement in the interactions between host and microbial organisms. For example, conserved domain data can be used to study phagocytosis, a process where host phygocytic cells (for example, macrophages) engulf pathogen cells (for example, Brucella). A search for 'phagocytosis' in Pacodom yields 14 domains; 13 domains do not match any protein from any PhiDB pathogen genome (Figure 4a). However, one domain, 'Nramp' (pfam01566), matches 42 pathogen proteins (Figure 4b). As summarized in the Pfam description of this domain (available in Pacodom), the natural resistance-associated macrophage protein (Nramp) family consists of Nramp1 and Nramp2 in human and mouse systems. Nramp1 plays an important role in phagocytosis and the macrophage activation pathway and regulates the interphagosomal replication of bacteria. Nramp2 is a transporter of multiple divalent cations (for example, Fe2+, Mn2+ and Zn2+) and is involved in a major transferrin-independent iron uptake system in mammals. The Pfam summary does not list any related microbial Nramp proteins. However, a Pacodom search shows Nramp is very common in the bacterial pathogens listed in PHIDIAS. Those 42 proteins containing the Nramp domain come from many bacterial species, such as Brucella spp., Mycobacterium tuberculosis, and Salmonella enterica. Nramp exists in all strains from these bacteria, whether the strain is pathogenic or non-pathogenic. In contrast, Nramp does not exist in the following species: Campylobacter jejuni, Clostridium perfringens, Coxiella burnetii, Francisella tularensis, and Rickettsia prowazekii. The Nramp domain has been investigated in depth in mycobacteria [31]. Since pathogenic mycobacteria survive within phagosomes, a nutrient-restricted environment, divalent cation transporters of the Nramp family in phagosomes and mycobacteria may compete for metals that are crucial for bacterial survival [31]. However, inactivation of mycobacterial Nramp, called Mramp, does not affect virulence in mice, suggesting a sufficient redundancy in the cation acquisition systems [32]. A more recent report [33] demonstrated that the Salmonella enterica serovar typhimurium (S. typhimurium) requires both of the divalent cation transport systems, MntH (Nramp1 homolog) and SitABCD (putative ABC iron and/or manganese transporter), for full virulence in congenic Nramp1-expressing mice. These results suggest that bacterial Nramp is required for pathogenesis in S. typhimurium and probably other bacteria by synchronizing with other redundant cation transport system(s) to compete for divalent cations with host cells. The role of Brucella Nramp in pathogenesis remains unclear and deserves further analysis. This example demonstrates how Pacadom can be used to find valuable information and form testable hypotheses by comparative analysis of conserved domains.

Figure 4
figure 4

Example of Pacodom applications. (a) Pacodom search of 'phagocytosis'. (b) There are 42 Nramp protein matches from 42 pathogen genomes of 15 microbial species available in Pacodom.

It is noted that the Nramp domain (pfam01566), while found in a list of pathogens in Pacodom, is also found in many bacterial species that are not pathogens. Therefore, it may be important for investigators to cross reference PHIDIAS search results against databases that contain both pathogen and non-pathogen species. Since Pacodom includes conserved domains from both pathogenic strains and non-pathogenic strains of the same microbial species, it can be used to find domains shown in pathogenic but not in non-pathogenic strains. For example, a query of 'bacteriophage' in Pacodom results in many conserved domains being found, such as Phage_Mu_Gp45 (pfam06890) and Phage_Mu_F (pfam04233), which exist in pathogenic E. coli O157:H7 strain Sakai but not in the benign K12 strain. Such domains have previously been reported as required for pathogenesis [34].

BLAST searches

Gene or protein sequences among different pathogen genomes can be analyzed by different BLAST search approaches. PHIDIAS BLAST uses the latest web server version of BLAST obtained from NCBI [35]. It includes regular BLAST services (blastn, blastp, blastx, tblastn, tblastx), PSI/PHI BLAST, Mega BLAST, RPS BLAST, and BLAST 2 sequences. The nucleotide and protein BLAST libraries contain sequences from all the 77 genomes of the 42 pathogens (Table 1). The 7,919 PSSMs available in Pacodom are combined to form a customized RPS BLAST library specifically used for the RPS BLAST program. The sequence libraries are updated periodically to reflect newly curated annotations and the addition of new genomes.

The approaches used with BLAST greatly help comparative studies for all the genes available in PhiDB. However, some gene annotations from certain genomes are not satisfactory. Based on sequence similarity, these are readily detected with BLAST. The PHIDIAS BLAST methods can also be used to find a group of pathogen genes using a seeding DNA or protein sequence. For example, a PHIDIAS blastp search for the protein sequence of human Nramp1 (also known as SLC11A1, RefSeq#: NP_000569) yields 65 hits from 77 pathogen genomes, most of which are attributable to a single putative manganese transport protein (MntH, which belongs to the Nramp family) found in different pathogens, including four Brucella strains. A blastp search using human Nramp2 (also known as SLC11A2, RefSeq#: NP_000608) as input yields similar hits. The BLAST search results are consistent with the analysis of conserved domains as described in the section on Pacodom above.

Phinfo: curated pathogen-host interaction general information

The Phinfo database module stores pathogen and PHI information curated from the biomedical literature and other curated databases. A major source of Phinfo data are PIML documents available from Virginia Bioinformatics Institute (VBI) [7]. A Java program was developed to extract PIML documents from the ToolBus/PathPort PIML XML database via the PathInfo web service [36]. An Extensible Stylesheet Language for Transformations (XSLT) script was developed to parse the PIML documents into a text-based SQL script. This in turn was used to insert the parsed data into a pre-designed MySQL database system. Phinfo also integrates data manually curated by the PHIDIAS curation team from PubMed literature and other databases such as KEGG [9]. Phinfo links to the Hazards in Animal Research Database (HazARD). This database was developed internally at the University of Michigan [8]. Pathobiology and management of laboratory animals administered USA NIAID/CDC priority pathogens are subjects of the HazARD database and can be searched with Phinfo [8]. Currently, Phinfo includes information for 36 pathogens and corresponding PHI information supported by 2,894 references.

Phinfo provides an integrative web interface for user-friendly querying and display of curated pathogen and PHI information. Two query programs are available in Phinfo: Keyword Search and Topic Search. The Keyword Search program allows queries for specific pathogen and PHI information. Such information is displayed with the searched keywords highlighted in color. The Topic Search program searches for one or many of 47 topics listed in the hierarchical structure (Figure 5). Compared to the native PIML XML database [7], the relational Phinfo database system provides secure storage, efficient querying, and database extendibility (that is, the ability to add new data categories). In addition, Phinfo provides links to public databases (for example, NCBI taxonomy, NCBI Gene database, and PubMed). Phinfo is also integrated with other PHIDIAS components. For example, Phinfo of Brucella spp. indicates that a PCR assay based on the B. abortus gene wboA (forward primer: TTAAGCGCTGATGCCATTTCCTTCAC, reverse primer: GCCAACCAACCCAAATGCTCACAA) has been used to differentiate B. abortus vaccine strain RB51 from other Brucella strains. Either of the primer sequences can be linked directly by clicking to local nucleotide BLAST analysis. Genes found from local BLAST searches are also linked to the PHIDIAS gene table (Figure 3). The wboA genes from four Brucella genomes are always the first four hits. Other microbial genes (for example, from Vibrio and Yersinia) are also found, indicating a possible cross-reaction during PCR assays and/or functional similarities among these genes.

Figure 5
figure 5

PhiDB Topic Search. The PhiDB Topic Search web interface is shown on the left and a comparison of immunoassays for diagnosis of B. melitensis and B. anthracis is shown on the right.

Phigen: pathogen-host interaction genes

The interactions between pathogen and host genes have been extensively studied in the post-genomic era [37]. However, most databases of genes and proteins focus on sequence annotation and function in a single cell species. Phigen focuses on functional annotation of pathogen genes and their interaction with host genes during the process of pathogen-host reactions. The main source of the PHI-related gene annotation comes from literature curation and data integration. The information about genes and/or proteins required for virulence, able to induce protective immune responses in hosts, or used for diagnosis, has been annotated and stored in the Phigen system. Phigen consists of two parts, pathogen gene search and manual curation submission.

Every pathogen gene may be involved in an interaction between the pathogen and its host. The pathogen gene search interface of Phigen allows users to search for any pathogen genes from the 77 genomes of the 42 pathogens available in PhiDB (Table 1). The Phigen search has a function for simple Boolean-powered keyword searches and an advanced topic search (Figure 6). The advanced topic search allows searching for PHI-specific information and generic features, including chromosomes and chromosomal position, RefSeq identifier, GenBank accession number, locus tag and name, molecular weight, pI, and description. Searched results can also be sorted in ascending or descending order. Molecular weight and pI data obtained in each search may be used to aid the interpretation of two-dimensional mass spectrometry data for proteomics analyses.

Figure 6
figure 6

Gene search web interface in Phigen.

Phigen provides an efficient online submission system for submitting of data for curation of pathogen genes, especially their roles in PHI. The information is fully referenced from peer-reviewed publications, with direct links to PubMed paper abstracts and full texts for additional details. Submitted information is critically reviewed and verified by reviewers prior to acceptance. Currently, Phigen has manually curated and stored more than 400 genes from 42 pathogens. Instead of altering records from other public databases, the curation is currently focusing on adding PHI-related information, such as host immune responses, gene mutations and resultant pathogenic changes in the host. In addition to integrated gene information, the PHI-specific information assists researchers in surveying, comparing, and studying gene-specific PHI mechanisms.

Phinet: pathogen-host interaction network curation, data exchange, and visualization

PHI has the ability to reveal complicated networks between pathogen and host molecules. Phinet is targeted at analyzing molecular networks responsible for PHI. Phinet data are stored in PhiDB and are derived from the MINetML XML database extracted through the web service, other curated databases (for example, KEGG), and manual annotation based on literature curation. Similar to that implemented in Phinfo, a Java program was developed to extract MINetML documents from the ToolBus/PathPort MINetML XML database via the MINet web service [38]. An XSLT script was further developed to parse the MINetML documents into a text-based SQL script, which is used to insert the parsed data into a pre-designed MySQL database system. Data from the KEGG pathway database are manually curated and added to Phinet. Phinet also includes a web-based data submission system that permits internal or external curators to submit PHI-related network data. The Phinet data submission follows a similar curation policy as described for Phigen online submission above. If conflicts exist for data from different sources, those records with the strongest reference support are selected, or in some circumstances, conflicting data were included with well-documented references. Currently, Phinet includes PHI network information for 21 pathogens.

A Graphviz-based visualization software program has been developed internally to dynamically display all the biological interactions in Phinet (Figure 7). The visualization program effectively displays all pathway data for each pathogen available in Phinet. The user can select to view information about a biological object or the interaction between biological objects (Figure 7).

Figure 7
figure 7

Visualization of an E. coli pathogenesis network in Phinet. A click on each node provides detailed information about a biological object in the bottom frame. When a mouse cursor moves over a node, a brief description of the biological object will appear. An interaction between biological objects is represented by a centered gray ball and arrows between nodes. Once the centered gray ball is clicked, details about the specific interaction appear in the bottom frame. Subcellular locations of biological objects are differentiated by the node border colors. The biological object types (for example, protein or gene) are represented by a combination of the node background colors and shapes. The program also displays different interactions, such as inhibition (solid T sign), activation (solid arrow), and indirect effects (dashed line).

Data exchange among different pathway databases is critical for data sharing and integration. BioPAX is a community-supported data exchange format for biological pathway data [14]. Current BioPAX Level 2 covers metabolic pathways, molecular interactions and protein post-translational modifications. Compared to the model representation format SBML, BioPAX focuses on molecule and interaction classification schemes and database cross-referencing for pathway components. PHI networks involve complex signaling pathways and gene regulatory networks that are similar to BioPAX, although they are not supported in their entirety by the current BioPAX version. A program was developed to transform Phinet data to the closest BioPAX OWL format using the current BioPAX Level 2 format. These BioPAX documents can be used to communicate with other biological pathway databases and, additionally, provide input files for other software programs.

Phix: pathogen-host interaction gene expression

Gene expression data for pathogens and/or hosts during PHIs comprise important data for analysis of pathogen pathogenesis and host defense mechanisms. The NCBI GEO [15] and EBI ArrayExpress [39] are the two biggest repositories that store publicly available microarray and proteomics data, many of which relate to PHI. The Phix database stores all gene expression experiment records for the targeted 42 pathogens and their infected hosts from the GEO and ArrayExpress databases. Since new gene expression experiments are frequently submitted to these databases, a Linux cron job [40] was developed to check daily for any new information; if found, the new data are added to the database. The Phix module currently stores 187 GEO records and 79 ArrayExpress records. The Phix gene expression search program provides a one-step system for users to query PHI gene expression experimental data. For example, a query of 'macrophage' in Phix leads to 13 search hits representing various experimental studies involving pathogen-infected macrophages. Each hit links to detailed information in GEO or ArrayExpress. These results are particularly useful for comparing different pathogen-macrophage interaction systems. Finally, Phix also includes a gene profile search engine for query and comparison of expression profiles of specific genes from one, or all, of the pathogen genomes selected from the GEO and ArrayExpress databases. In contrast to the general GEO and ArrayExpress gene profile search engines, this program is specifically targeted to pathogen and PHI studies.

To improve further integration of different PHIDIAS components, the PHIDIAS web site contains a keyword search engine that simultaneously allows searching for information from all PHIDIAS components. All results are sorted based on the components and displayed in one page for convenient data analysis (data not shown).

Discussion

A deeper understanding of PHI is required for effectively combating infectious diseases. To efficiently analyze the ever-increasing amount of PHI data in the post-genomics era, PHIDIAS was developed. This program permits integration of PHI related data from genome sequences, the biomedical literature, curated databases, and gene expression experiments. PHIDIAS covers 42 microbial and viral pathogens of high priority for public heath and security. The gene and protein sequences from each genome are available for browsing and analysis using PGBrowser and customized BLAST searches. The conserved domains are analyzed and stored in Pacadom. PHI data extracted from existing databases, or internally manually curated, are stored in Phinfo (general PHI information), Phigen (PHI genes) and Phinet (PHI networks). PHI-related gene expression experiment records and profiles from public GEO and ArrayExpress repositories can be directly searched in Phix. The PHIDIAS components are interconnected (Figure 1). Scenarios have been used in this report to show that PHIDIAS greatly helps Brucella research by allowing users to search and analyze integrative Brucella data derived from different sources and compare these data with those from other pathogens.

Similar PHI-related biological programs exist. PHI-base is a web-accessible database devoted to the identification and presentation of information on fungal and oomycete pathogenicity genes and their host interactions [41]. PathoPlant deals with plant-pathogen interactions, signal transduction reactions, and microarray gene expression data from Arabidopsis thaliana subjected to pathogen infection and elicitor treatment [42]. In contrast to PHI-base and PathoPlant, which target the interactive relationships between pathogens and hosts, PHIDIAS includes a list of other bacterial, viral and parasitic pathogens and their interactions with hosts. Similar to PHIDIAS, PHI-base and PathoPlant contain manually curated information supported by strong experimental evidence (gene disruption experiments) and literature references. Each system allows interlinking of gene information with external data sources. However, PHIDIAS integrates more data sources for a broader scope of data integration and analysis. PHIDIAS also provides on-line submission systems for curators to submit annotated data for genes as well as genetic interactions and pathways.

Many biological systems allow systematic genome comparison. MicrobesOnline is a publicly available suite of web-based comparative genomic tools designed to facilitate multispecies comparison among prokaryotes [43]. The database PRODORIC systematically organizes information about the prokaryotic gene expression of multiple prokaryotic species, and integrates this information into regulatory networks [44]. As does PHIDIAS, these systems contain many comparative analysis and visualization tools. However, while MicrobesOnline and PRODORIC target more general prokaryotic species, PHIDIAS focuses on pathogenic bacteria as well as viral and parasitic pathogens important for biodefense and/or human health. PHIDIAS also emphasizes interactions between pathogens and hosts, which MicrobesOnline and PRODORIC currently lack. PHIDIAS also contains manually curated data for functional annotation of genes and genetic networks in pathogen genomes.

Eight Bioinformatics Resource Centers (BRCs), sponsored by the USA NIAID, provide web-based resources for organisms that are considered potential agents of biowarfare or bioterrorism or cause emerging or re-emerging diseases [45]. Each BRC is targeted to maintain and annotate genomes from a selected list of pathogens. Each BRC contains a web site to display the data and analyses for these pathogens. BRC Central [46] serves as a repository linking these eight BRCs. Many of the pathogens contained in the BRCs are also found in PHIDIAS. However, PHIDIAS also targets non-biodefense pathogens (for example, HIV) not included in the BRCs. Additionally, PHIDIAS includes not only data analysis and search functions found in the BRC resources, but also provides tighter integration of various data types. Finally, PHI and literature data curation are emphasized in PHIDIAS but not in the BRCs.

PHIDIAS is unique in that it integrates existing knowledge about a broad range of human or zoonotic priority pathogens, and focuses on efficient searching, visualization, comparison, and analysis of pathogen genes and their interactions with their hosts using genome sequences, manually curated literature data, and gene expression data from public resources. PHIDIAS utilizes online data submission systems for efficient data curation, making integrative PHI data more comprehensive. All the PHIDIAS components are scalable, and more pathogens and PHI systems may be added to the system. Due to inclusion of an ever increasing number of pathogens in PHIDIAS and in view of the dramatically increasing amount of literature information, it will be an ongoing challenge to curate all the significant genes and keep the PHI-related information in PhiDB current. Therefore, one of our future directions will be to explore ontology-based natural language processing and statistical methods for efficient literature acquisition and curation. In this regard, we have now developed a literature mining and curation system (Limix). This system has been used efficiently for literature mining and curation for four Brucella genomes [2]. Systematic curation and incorporation of Brucella-specific mutation and genetic interaction information has allowed a comprehensive investigation of Brucella pathogenesis [2]. Limix is currently being expanded to annotate literature for other pathogens and PHI systems. Finally, future plans for expanding PHIDIAS include development of a web-based database and an analysis pipeline that permit storage, processing, and modeling of PHI-related gene expression data. This approach will allow researchers to address scientific PHI questions with the ultimate goal of successfully fighting infectious diseases.

References

  1. Becker K, Hu Y, Biller-Andorno N: Infectious diseases - a global challenge. Int J Med Microbiol. 2006, 296: 179-185. 10.1016/j.ijmm.2005.12.015.

    Article  PubMed  Google Scholar 

  2. Xiang Z, Zheng W, He Y: BBP: Brucella genome annotation with literature mining and curation. BMC Bioinformatics. 2006, 7: 347-10.1186/1471-2105-7-347.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, et al: The Pfam protein families database. Nucleic Acids Res. 2004, 32: D138-141. 10.1093/nar/gkh121.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, Bork P: SMART 5: domains in the context of genomes and networks. Nucleic Acids Res. 2006, 34: D257-260. 10.1093/nar/gkj079.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, et al: The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003, 4: 41-10.1186/1471-2105-4-41.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Marchler-Bauer A, Anderson JB, Derbyshire MK, DeWeese-Scott C, Gonzales NR, Gwadz M, Hao L, He S, Hurwitz DI, Jackson JD, et al: CDD: a conserved domain database for interactive domain family analysis. Nucleic Acids Res. 2007, 35: D237-240. 10.1093/nar/gkl951.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. He Y, Vines RR, Wattam AR, Abramochkin GV, Dickerman AW, Eckart JD, Sobral BW: PIML: the Pathogen Information Markup Language. Bioinformatics. 2005, 21: 116-121. 10.1093/bioinformatics/bth462.

    Article  PubMed  CAS  Google Scholar 

  8. He Y, Rush HG, Liepman RS, Xiang Z, Colby LA: Pathobiology and management of laboratory rodents administered CDC Category A agents. Comparative Med. 2007, 57: 18-32.

    CAS  Google Scholar 

  9. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32: D277-280. 10.1093/nar/gkh063.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  10. Karp PD, Riley M, Saier M, Paulsen IT, Collado-Vides J, Paley SM, Pellegrini-Toole A, Bonavides C, Gama-Castro S: The EcoCyc Database. Nucleic Acids Res. 2002, 30: 56-58. 10.1093/nar/30.1.56.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Krieger CJ, Zhang P, Mueller LA, Wang A, Paley S, Arnaud M, Pick J, Rhee SY, Karp PD: MetaCyc: a multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res. 2004, 32: D438-442. 10.1093/nar/gkh100.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  12. Bader GD, Betel D, Hogue CW: BIND: the Biomolecular Interaction Network Database. Nucleic Acids Res. 2003, 31: 248-250. 10.1093/nar/gkg056.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. The Molecular Interaction Network Markup Language(MINetML). [http://pathport.vbi.vt.edu/xml/molecules/molecules.dtd]

  14. Stromback L, Lambrix P: Representations of molecular pathways: an evaluation of SBML, PSI MI and BioPAX. Bioinformatics. 2005, 21: 4401-4407. 10.1093/bioinformatics/bti718.

    Article  PubMed  CAS  Google Scholar 

  15. Barrett T, Suzek TO, Troup DB, Wilhite SE, Ngau WC, Ledoux P, Rudnev D, Lash AE, Fujibuchi W, Edgar R: NCBI GEO: mining millions of expression profiles - database and tools. Nucleic Acids Res. 2005, 33: D562-566. 10.1093/nar/gki022.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Parkinson H, Sarkans U, Shojatalab M, Abeygunawardena N, Contrino S, Coulson R, Farne A, Lara GG, Holloway E, Kapushesky M, et al: ArrayExpress - a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2005, 33: D553-555. 10.1093/nar/gki056.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Roop RM, Bellaire BH, Valderas MW, Cardelli JA: Adaptation of the brucellae to their intracellular niche. Mol Microbiol. 2004, 52: 621-630. 10.1111/j.1365-2958.2004.04017.x.

    Article  PubMed  CAS  Google Scholar 

  18. DelVecchio VG, Kapatral V, Redkar RJ, Patra G, Mujer C, Los T, Ivanova N, Anderson I, Bhattacharyya A, Lykidis A, et al: The genome sequence of the facultative intracellular pathogen Brucella melitensis. Proc Natl Acad Sci USA. 2002, 99: 443-448. 10.1073/pnas.221575398.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Paulsen IT, Seshadri R, Nelson KE, Eisen JA, Heidelberg JF, Read TD, Dodson RJ, Umayam L, Brinkac LM, Beanan MJ, et al: The Brucella suis genome reveals fundamental similarities between animal and plant pathogens and symbionts. Proc Natl Acad Sci USA. 2002, 99: 13148-13153. 10.1073/pnas.192319099.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Halling SM, Peterson-Burch BD, Bricker BJ, Zuerner RL, Qing Z, Li LL, Kapur V, Alt DP, Olsen SC: Completion of the genome sequence of Brucella abortus and comparison to the highly similar genomes of Brucella melitensis and Brucella suis. J Bacteriol. 2005, 187: 2715-2726. 10.1128/JB.187.8.2715-2726.2005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Chain PS, Comerci DJ, Tolmasky ME, Larimer FW, Malfatti SA, Vergez LM, Aguero F, Land ML, Ugalde RA, Garcia E: Whole-genome analyses of speciation events in pathogenic brucellae. Infect Immun. 2005, 73: 8353-8361. 10.1128/IAI.73.12.8353-8361.2005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. BioPerl. [http://www.bioperl.org]

  23. Stein LD, Mungall C, Shu S, Caudy M, Mangone M, Day A, Nickerson E, Stajich JE, Harris TW, Arva A, Lewis S: The generic genome browser: a building block for a model organism system database. Genome Res. 2002, 12: 1599-1610. 10.1101/gr.403602.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Winsor GL, Lo R, Sui SJ, Ung KS, Huang S, Cheng D, Ching WK, Hancock RE, Brinkman FS: Pseudomonas aeruginosa Genome Database and PseudoCAP: facilitating community-based, continually updated, genome annotation. Nucleic Acids Res. 2005, 33: D338-343. 10.1093/nar/gki047.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Gee JM, Valderas MW, Kovach ME, Grippe VK, Robertson GT, Ng WL, Richardson JM, Winkler ME, Roop RM: The Brucella abortus Cu, Zn superoxide dismutase is required for optimal resistance to oxidative killing by murine macrophages and wild-type virulence in experimentally infected mice. Infect Immun. 2005, 73: 2873-2880. 10.1128/IAI.73.5.2873-2880.2005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. He Y, Vemulapalli R, Schurig GG: Recombinant Ochrobactrum anthropi expressing Brucella abortus Cu, Zn superoxide dismutase protects mice against B. abortus infection only after switching of immune responses to Th1 type. Infect Immun. 2002, 70: 2535-2543. 10.1128/IAI.70.5.2535-2543.2002.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  27. Passalacqua KD, Bergman NH, Herring-Palmer A, Hanna P: The superoxide dismutases of Bacillus anthracis do not cooperatively protect against endogenous superoxide stress. J Bacteriol. 2006, 188: 3837-3848. 10.1128/JB.00239-06.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. NCBI CDD Download. [ftp://ftp.ncbi.nih.gov/pub/mmdb/cdd/cdd.tar.gz]

  29. NCBI Toolkit Download. [ftp://ftp.ncbi.nlm.nih.gov/toolbox]

  30. Marchler-Bauer A, Panchenko AR, Shoemaker BA, Thiessen PA, Geer LY, Bryant SH: CDD: a database of conserved domain alignments with links to domain three-dimensional structure. Nucleic Acids Res. 2002, 30: 281-283. 10.1093/nar/30.1.281.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  31. Agranoff D, Monahan IM, Mangan JA, Butcher PD, Krishna S: Mycobacterium tuberculosis expresses a novel pH-dependent divalent cation transporter belonging to the Nramp family. J Exp Med. 1999, 190: 717-724. 10.1084/jem.190.5.717.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Boechat N, Lagier-Roger B, Petit S, Bordat Y, Rauzier J, Hance AJ, Gicquel B, Reyrat JM: Disruption of the gene homologous to mammalian Nramp1 in Mycobacterium tuberculosis does not affect virulence in mice. Infect Immun. 2002, 70: 4124-4131. 10.1128/IAI.70.8.4124-4131.2002.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. Zaharik ML, Cullen VL, Fung AM, Libby SJ, Kujat Choy SL, Coburn B, Kehres DG, Maguire ME, Fang FC, Finlay BB: The Salmonella enterica serovar typhimurium divalent cation transport systems MntH and SitABCD are essential for virulence in an Nramp1G169 murine typhoid model. Infect Immun. 2004, 72: 5522-5525. 10.1128/IAI.72.9.5522-5525.2004.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Hayashi T, Makino K, Ohnishi M, Kurokawa K, Ishii K, Yokoyama K, Han CG, Ohtsubo E, Nakayama K, Murata T, et al: Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12. DNA Res. 2001, 8: 11-22. 10.1093/dnares/8.1.11.

    Article  PubMed  CAS  Google Scholar 

  35. NCBI BLAST Download. [http://www.ncbi.nih.gov/BLAST/download.shtml]

  36. PathInfo Web Service. [http://staff.vbi.vt.edu/pathport/services/wsdls/pathinfo.wsdl]

  37. Forst CV: Host-pathogen systems biology. Drug Discov Today. 2006, 11: 220-227. 10.1016/S1359-6446(05)03735-9.

    Article  PubMed  CAS  Google Scholar 

  38. MINet Web Service. [http://www.vbi.vt.edu/~pathport/services/wsdls/pathway.wsdl]

  39. Brazma A, Parkinson H, Sarkans U, Shojatalab M, Vilo J, Abeygunawardena N, Holloway E, Kapushesky M, Kemmeren P, Lara GG, et al: ArrayExpress - a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2003, 31: 68-71. 10.1093/nar/gkg091.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  40. Petersen R: Linux: The Complete Reference. 2000, Emeryville, CA: McGraw-Hill Osborne Media, 4

    Google Scholar 

  41. Winnenburg R, Baldwin TK, Urban M, Rawlings C, Kohler J, Hammond-Kosack KE: PHI-base: a new database for pathogen host interactions. Nucleic Acids Res. 2006, 34: D459-464. 10.1093/nar/gkj047.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  42. Bulow L, Schindler M, Hehl R: PathoPlant: a platform for microarray expression data to analyze co-regulated genes involved in plant defense responses. Nucleic Acids Res. 2007, 35: D841-845. 10.1093/nar/gkl835.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Alm EJ, Huang KH, Price MN, Koche RP, Keller K, Dubchak IL, Arkin AP: The MicrobesOnline Web site for comparative genomics. Genome Res. 2005, 15: 1015-1022. 10.1101/gr.3844805.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  44. Munch R, Hiller K, Barg H, Heldt D, Linz S, Wingender E, Jahn D: PRODORIC: prokaryotic database of gene regulation. Nucleic Acids Res. 2003, 31: 266-269. 10.1093/nar/gkg037.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  45. NIAID Bioinformatics Resource Centers for Biodefense and Emerging or Re-emerging Infectious Diseases: an Overview. [http://www.niaid.nih.gov/dmid/genomes/brc/]

  46. BRC Central. [http://www.brc-central.org]

Download references

Acknowledgements

We thank the authors of published data in various programs (for example, RefSeq, CDD, Pfam, PubMed, PathInfo, MINet, HazARD, KEGG, and so on) for making them available to the public. We also acknowledge the public availability of many open-source programs (for example, GBrowse and NCBI BLAST) that have allowed the integration and extension into PHIDIAS. The critical review and editing of this manuscript by Drs L Colby and GW Jourdian from the University of Michigan Medical School is gratefully acknowledged.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yongqun He.

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Xiang, Z., Tian, Y. & He, Y. PHIDIAS: a pathogen-host interaction data integration and analysis system. Genome Biol 8, R150 (2007). https://doi.org/10.1186/gb-2007-8-7-r150

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2007-8-7-r150

Keywords