- Research article
- Open Access
Encyclopedia of bacterial gene circuits whose presence or absence correlate with pathogenicity – a large-scale system analysis of decoded bacterial genomes
BMC Genomics volume 16, Article number: 773 (2015)
Bacterial infections comprise a global health challenge as the incidences of antibiotic resistance increase. Pathogenic potential of bacteria has been shown to be context dependent, varying in response to environment and even within the strains of the same genus.
We used the KEGG repository and extensive literature searches to identify among the 2527 bacterial genomes in the literature those implicated as pathogenic to the host, including those which show pathogenicity in a context dependent manner. Using data on the gene contents of these genomes, we identified sets of genes highly abundant in pathogenic but relatively absent in commensal strains and vice versa. In addition, we carried out genome comparison within a genus for the seventeen largest genera in our genome collection. We projected the resultant lists of ortholog genes onto KEGG bacterial pathways to identify clusters and circuits, which can be linked to either pathogenicity or synergy. Gene circuits relatively abundant in nonpathogenic bacteria often mediated biosynthesis of antibiotics. Other synergy-linked circuits reduced drug-induced toxicity. Pathogen-abundant gene circuits included modules in one-carbon folate, two-component system, type-3 secretion system, and peptidoglycan biosynthesis. Antibiotics-resistant bacterial strains possessed genes modulating phagocytosis, vesicle trafficking, cytoskeletal reorganization, and regulation of the inflammatory response. Our study also identified bacterial genera containing a circuit, elements of which were previously linked to Alzheimers disease.
Present study produces for the first time, a signature, in the form of a robust list of gene circuitry whose presence or absence could potentially define the pathogenicity of a microbiome. Extensive literature search substantiated a bulk majority of the commensal and pathogenic circuitry in our predicted list. Scanning microbiome libraries for these circuitry motifs will provide further insights into the complex and context dependent pathogenicity of bacteria.
Microbiology experiments identified a large number of bacterial virulence mechanisms conserved through evolution . Pathogenic bacterial strains could be defined as those with capacity to harm the host and cause disease . A bacterial strain may appear as asymptotic or pathogenic depending on the state of the immune system of the host , composition of the microbiome , presence and absence of elicitors , and other environmental factors . Synthesis of secreted virulence factors are transcriptionally regulated by environmental stimuli . Production of enzymes that degrade host cytoskeleton and cause damage is linked to bacterial strain density and quorum sensing . It is clear that pathogenicity is context dependent.
The large research literature on bacterial pathogenicity was recently curated into an open access web platform, the Virulence Factor Database or VFDB [2, 6, 9]. The database presents genes and gene groups (virulence factors) associated with infectious disease. Since the metagenome of a human microbiome contains over three million genes, the extent of our knowledge on virulence genes may not be close to saturation. Nevertheless, not only the VFDB database is valuable for identifying genes linked to pathogenicity in specified bacteria (even strains), but it also serves as a benchmark for studies predicting pathogenicity linked gene clusters.
Virulence factors may cause disease in multiple hosts . Yet, specificity of virulence factors to certain bacterial and strains is poorly understood. A recent study involving 50 bacterial genomes found some virulence factors to be exclusive to pathogenic bacterial strains in this small sample . Virulence genes abundant in pathogenic as well as nonpathogenic strains included those facilitating coding of translocation proteins, apparatus proteins, and chaperons. The findings of the study suggest the presence of pathogenic gene circuits in which some genes belong exclusively to pathogenic strains whereas others can also be abundant in bacterial strains synergistic or commensal to the host.
The drivers of pathogenicity may not only be virulence factors but also the absence of antivirulence factors in bacterial genomes [11–13]. Some bacterial species contain strains known as commensal to a host and other strains that exhibit pathogenicity. Transformation to pathogenic state maybe due to exchange of DNA fragments between bacteria . Literature points out to the acquisitions of gene clusters and pathogenicity islands via horizontal gene transfer [15–17]. Deletion of genes or loss of gene function through mutation appears to be part of the adaptation to newly acquired pathogenesis .
Recent identification of the elicitors for activating antibiotics-synthesizing bacterial gene circuits  provides a new dimension in our understanding of the environmental forces affecting pathogenicity. Although, the research literature contains many examples of bacterial genes linked to synthesis of antibiotics [19–21], a large-scale study exposing diversity of such genes and gene circuits is yet to be carried out. The same is true for pathogenic gene circuits. That is the reason we set out in this study a systematic approach to identify and annotate prokaryotic ortholog gene circuits whose presence or absence are linked to phenotypes of antagonism to the host. Comparative genomics has been utilized before in the identification of drug and vaccine targets in Staphylococcus aureus  and mycobacterial peptidoglycan remodeling enzymes linked to pathogenicity . Ours, however, takes it one step further, employing comparisons among 2527 distinct genomes. Restriction of genome comparison to ortholog groupings reduced dimension of this meta-scale analysis. It also expanded the reach of findings within the context of evolution. Aspects of pathogenicity, which are universal across eukaryotic hosts, will likely emerge in this approach .
Our methodology has multiple steps. First, we annotate existing genome sequences of bacterial strains as pathogenic or otherwise, based on literature curation. Next, we identify the presence and absence of orthologs in the genomes of these bacterial strains. Thirdly, we determine the relative abundance or absence of these orthologs in the pathogenic and non pathogenic strains within and across genera. Projection onto cellular pathways result in annotation of gene circuits linked to pathogenicity. Extensive comparison with experimental literature provides biological context to our findings. Our study in effect creates an encyclopedia for pathogenicity, built on big data on genome sequences and literature on phenotypes of bacterial strains.
A. Orthologs linked either to pathogenicity or synergy with the host
Annotation of pathogenic bacterial strains
Our literature search detailed in the Methods section identified 949 bacterial strains as pathogenic, meaning they had been reported in the literature as pathogenic to an animal host at least once. The label pathogenic, in the way we use it, doesn’t mean that the bacterial strain will cause disease to an animal host under all circumstances.
Additional file 1 presents the list of bacterial strains deemed as pathogenic, with evidence provided in the file in the form of references or database citations. This supplemental file also contains labeling of pathogenic strains found as antibiotic resistant in the Antibiotic Resistance Genes Database (ARDB) . Table 1 shows a sample of the bacterial genera possessing both pathogenic and nonpathogenic strains. Strains of the same bacterial genus often separated into pathogenic and nonpathogenic clusters.
There was no KEGG- or literature-recorded evidence of pathogenicity for the remaining 1578 decoded bacterial strain genomes in the KEGG database. Thus, they were deemed for this study as nonpathogenic. Additional file 2 presents the list of the 1578 bacterial strains deemed nonpathogenic in our study.
Orthology contents of pathogenic and nonpathogenic strains
Next, we created two sets of matrices for ortholog genes, one set for the genomes of pathogenic and the other for nonpathogenic strains. The columns of each matrix identified the bacterial strain whereas rows identified whether an ortholog was present (1) or absent (0) in that strain. It turned out that 7194 different orthologs were present in at least one of the genomes of the 2527 bacterial strains under study. These large matrices are used for the abundance computations presented in this study and hence included as Additional file 3.
For a given ortholog, fractions of pathogenic and nonpathogenic bacterial strains presenting the ortholog in their genomes are represented by symbols Ap and Anp, respectively. The scatter diagram shown in Fig. 1 presents the Ap and Anp values for the 7194 orthologs present in bacterial strains. It appears that most of the orthologs have comparable presence in all bacterial strains whereas a small portion (green and red dots for pathogenic and nonpathogenic strains, respectively) is biased towards one of the two phenotypes.
The histogram shown in Fig. 2 is another view of the data presented in the scatter diagram in Fig. 1. Here, we plotted the frequency of occurrence against the pathogen abundance score log PA for all the orthologs under consideration. The parameter PA = Ap / (Anp + 0.0001) is a measure of relative abundance of the ortholog in pathogenic strains. In cases where Anp equaled zero, the equation still enables division due to the presence of 0.0001 in the denominator. The two tail ends of the distribution indicate those orthologs abundant in pathogenic but rarely found in nonpathogenic (red) and vice versa (green). The cutoff values we used (PA > 4, and PA < ¼), although somehow arbitrary, were placed at the inner edges of the tails of the histograms.
Shown in Fig. 3a are the overall characteristics of the PA distribution among the orthologs. In brief, there were 229 pathogenic only and an additional 379 orthologs for which PA > 4, representing about 8 percent of the orthologs found in our bacterial strain library. Taken together, we deemed this group as pathogen abundant or pathogen-linked. Total number of genes in the 608 pathogen-abundant orthologs was 18,982, indicating their presence in a diverse set of bacterial species. Pathogen exclusive orthologs comprised only 1,518 of this set of genes, suggesting most genes previously linked to pathogenicity is not exclusive to disease-causing bacterial strains.
The number of orthologs exclusive to nonpathogenic strains was much larger at 879, and an additional 485 had PA < ¼. The rest, a total of 5222 orthologs, were commonly present among pathogenic and nonpathogenic strains. It is expected that these numbers will change as the number of available bacterial genomes in the literature increase from thousands to tens of thousands.
Next, we identified those orthologs in our list, a total of 1308, which were also present in the Virulence Factor Database, VFDB, by matching either the gene names or gene descriptions. As indicated in Fig. 3b, VFDB orthologs are significantly biased towards pathogen abundant orthologs. Additional file 4 lists the pathogen-abundant and non-pathogen-abundant orthologs in accordance with the PA ranking, along with VFDB labeling if present in that database. Overall, our study indicates the absence of one-to-one match between known virulence factors and pathogen-abundant orthologs.
Orthologs enriched in pathogenic, antibiotic resistant, and nonpathogenic bacterial strains
Statistical enrichment of KEGG pathways was conducted based on the hypergeometric test for ortholog sets PA > 4, and PA < ¼, respectively. Results are shown in Fig. 4. Orthologs abundant in pathogenic strains crowd pathogen-linked cellular pathways: Staphyloccus aureus, Leigonellosis, Pertussis, Salmonella, Shigellosis, and Escherichia coli infections, as well as epithelial signaling in H. pylori infection. Pathogen abundant orthologs are also found in pathways involving bacterial secretion systems, Nod like receptor signaling, bacterial invasion of epithelial cells, and plant-pathogen interactions.
Orthologs found exclusively in nonpathogenic strains occupy nodes in metabolic pathways (Fig. 4). These pathways include biosynthesis of peptidoglycans, microlides, carotenoids, ansamycins, and nonribiosomal peptides. The GO cell compartment investigations not shown in the figure indicate that pathogen abundant orthologs are enriched in crosstalk positions of contact with the host whereas nonpathogen exclusive orthologs code proteins involving in events in the cell interior.
Next, we looked at molecular function enrichments of pathogenic associated orthologs and compared the results with corresponding enrichments obtained using the VFDB database. Figure 5 shows that both our annotation and VFDB contain roughly equal amounts of orthologs in secretion, toxins, peptidase, and pilin categories. However, pathogen abundant list of the present study has significantly more abundance in enzyme categories such as oxidoreductases, transferases, hydrolases, lyases, and ligases. Some of the orthologs in our list are also enriched in regulatory function, particularly involving transcription and translation.
Pathway enrichment was also conducted within the population of pathogenic strains for the subset identified as antibiotic resistant using ARDB Database . The hypergeometric test revealed the pathways shown in Table 2 as particularly enriched in antibiotic resistant strains. These included sphingolipid metabolism, producing bioactive metabolites that regulate cell function , PI3K-Akt signaling pathway, an intracellular pathway important in apoptosis , and Aminoacyl-tRNA biosynthesis . Some of the modules in the enriched pathways also appear in eukaryotic processes for drug resistance against chemotherapy. One must caution, however, that the results could potentially change with the updating of ARDB, even if the p values in these enrichments are vanishingly small.
B. Gene circuits linked either to pathogenicity or synergy
This section presents results on genetic circuits statistically enriched in pathogenic and nonpathogenic bacteria. We used two different types of comparison to achieve our results: a) analyzing the entire set of genomes partitioned into pathogenic and nonpathogenic phenotypes; and b) conducting the same operation within genera for the 17 genera identified in Table 1 with a star. In the first approach, we mapped the list of ortholog genes linked to pathogenicity (PA > 4) and non-pathogenicity onto KEGG reference pathways and identified, based on KEGG repository, those multiply connected clusters of genes (gene circuits) containing at least three pathogen-linked or nonpathogen-linked orthologs. Results are shown in Tables 3 and 4, respectively. Additional file 5 presents corresponding results for within-genera comparisons, both for pathogen- and nonpathogen-linked circuits along with the genera containing such circuits.
Gene clusters more common in pathogenic strains
The Table 3 presents a set of gene circuits with the ortholog genes linked to pathogenicity and also indicates the pathway to which the gene circuits belong. Along with Table 3, comes Fig. 6, in which the wiring diagrams for the gene circuitry are shown in the form of cutouts from the KEGG Reference pathways. The circuits in the figure have the same ordering number used in Table 3. Note also that the actual circuits contain orthologs not only pathogen-linked (shown in pink and orange) but also others, some preferentially found in pathogenic strains and others more ubiquitous. The p value through hypergeometric test for a bacterial strain containing at least one pathogen-ortholog in a circuit shown in Table 3 was less than 0.01.
Gene circuits for bacterial secretion and invasion pathways:
The type III bacterial secretion system pathway mediates toxin and protein delivery to host cells. Table 3 shows the existence of multiple clusters of pathogen-linked orthologs in this pathway. Consistent with our findings, the type III pathway is listed in the literature as modulating pathogenic interactions with host organisms including animals and plants [29–32].
Also shown in Fig. 6 are examples of pathogen-linked circuits in the secretion system. One such circuit contains pathogen-linked orthologs yscF, yscO, yscP, yscX, yscC, and yscW. Subsets of pathogen-linked orthologs of the cluster are present in 242 pathogenic and 106 nonpathogenic strains, resulting in vanishing p values in hypergeometric test. Moreover, the bias towards pathogenicity increases dramatically with increasing number of pathogen-linked orthologs in this cluster in the genome of the bacterial strain.
Another secretion-linked pathway is that of type IV gene circuit, for which some of its genes exist in both pathogenic and nonpathogenic strains. The circuit functions in translocation of DNA and protein substrates to target cells via direct cell-to-cell contact. In our study, the complete circuit is preferentially present in pathogenic strains. Consistent with these observations, recent investigations uncovered a role for pathogenicity for this circuit [33–35].
Pathogen-linked gene clusters in the two-component System:
The two-component regulatory system is a stimulus–response coupling pathway, which enables bacteria to sense and respond to changes in its environment [36–40]. Membrane-bound histidine kinases are major building blocks of the pathway.
These signal transduction systems modulate crosstalk between species within the microbiome. The Table 3contains multiple gene clusters (circuits) in the two-component system containing orthologs linked to pathogenicity: cluster (devS, nreA, nreB, narT) involved in hypoxia, oxygen, and nitrogen assimilation; cluster (uhpC, uhpA, uhpT) modulating hexose phosphate uptake; and the cluster (pagC, pagO, pagD, pagK, pgtE) involved in Mg2+ starvation, and others. See also Fig. 6 for the wiring diagrams of these clusters. Elements of the metabolite assimilation cluster have been linked in the literature to pathogenicity of mycobacterium tuberculosis [41, 42]. The second cluster in the list in the two-component system, mediating hexose phosphate uptake, plays an important role in the sodium-dependent D-glucose transport protein of Helicobacter pylori . This gene circuit is involved in Mg2+ starvation and was shown to play a role in the pathogenicity of Salmonella enterica . Mg2+ starvation is also involved in quorum sensing of Pseudomonas fluorescens  and in biosynthesis of complex lipids needed for virulence of mycobacterium tuberculosis .
Metabolic circuits linked to pathogenicity:
A metabolic gene circuit whose genes are commonly found in pathogenic strains is the CMP-Pse metabolism circuit cluster belonging to the amino nucleotide sugar mechanism. Pathogen-linked ortholog genes in this circuit consist of pseC, pseH, pseF, and UAP1 (Table 3). This circuit is linked to the synthesis of glycoconjugates, which are typically expressed on the surfaces of pathogenic bacteria. The protein products of the circuit have already been identified as virulence factors in the VFDB database and in the literature [47–49].
Nodal elements of the Peptidoglycan biosynthesis circuit cluster shown in Table 3 are also preferentially present in pathogenic strains. Pathogen-linked orthologs in this gene circuit consist of the genes sgtA, sgtB, femA, pbpA, femB, pbp3, femX, and fmhB. Peptidoglycans are polymers consisting of sugars and amino acids forming a mesh scaffold external to the plasma membrane. Recent studies in the literature point to the role of peptidoglycans in the pathogen phenotype of different bacteria [50–52].
Sorbose to Sorbose 1-phosphate circuit of the Phosphotransferase (PTS) system also shown in Table 3 contains pathogen-linked genes PTS-Sor-EIIC, sorA, PTS-Sor-EIID, sorM, PTS-Sor-EIIA, sorF, PTS-Sor-EIIB, sorB. PTS circuit codes a group translocation process present in many bacteria, transporting sugars from the environment into the bacterial cell. The circuit has been linked in the literature to Streptococcus invasion .
Our statistical computations based on hypergeometric test indicate that the likelihood of pathogenic identification of a strain increases dramatically with increasing numbers of the circuit genes expressed in the strain’s genome. Reflecting this finding, the Table 3 contains 15 gene circuits for which bacterial strains containing at least 75 percent of the circuit elements are always pathogenic. Hence a signature for pathogenicity may be derived from the study of clusters of pathogen-linked orthologs in bacterial strains.
Figure 6 presents other examples of ortholog groupings listed in Table 3 and acting in tandem in pathogenic processes. One such circuitry shown in the figure is involved in the biosynthesis of siderophore group of nonribosomal peptides. These are high affinity iron binding compounds  and were found to play an important role in virulent bacterial infections . Also shown in the Figure is pathogen associated ortholog circuit clusters crowding the bacterial secretion system not discussed above in detail. As noted in the literature, the secretion system facilitates transport, injection, and release of effector compounds including enzymes, and toxins in bacteria [56, 57].
Additional pathogen-linked circuitry identified through comparisons of genomes belonging to the same genera:
Additional circuits linked to pathogenicity could be identified using within-genera genome comparisons. We have conducted comparisons of ortholog contents of strains belonging to the same genera for the seventeen genera with most number of strains in our dataset, shown in Table 1. Again, the clusters of pathogen-linked orthogs forming on KEGG reference pathways were identified. However, in this case, we relaxed the pathogen-linkage evaluation from PA > 4 to PA > 2 since genomes belonging to the same genera are more or less similar. In addition, we are looking here for circuitry common across genera.
Results of these computations are presented in Additional file 5, identified in rows 1 to 21 for pathogen-linked circuits. The table shows not only the circuitry but also the genera associated with the specified circuitry. The circuit clusters most common across genera in this Table lie in the pathways for glycine, serine and threonine metabolism and sulfur metabolism (Fig. 7). These pathways have been implicated in playing important roles in pathogenicity [58–61]. Also in this category, is the gene circuit in Additional file 5 row 1 linked to Alzheimer's disease via Amyloid B and Mitochondial Disfunction.
Another pathogen-linked gene circuitry that comes out in genera-specific comparisons is the Amyloid B and Mitochondrial Dysfunction circuitry in the KEGG Alzheimer’s pathway (Additional file 5). The pathogen-linked orthologs in this circuitry (UQCRFS1, RIP1, petA, MME, IDE, ide, CALM, NDUFV2) are found in 16 of the 17 genera under consideration. This observation suggests the diversity of a bacterial infection that could be linked as a possible modulator of the Alzheimer’s disease [62–64].
Gene circuits found in nonpathogenic strains
Circuits linked to nonpathogenicity are shown in Tables 4 and Additional file 5, respectively, for across genera and within genera comparisons. Our detailed results shown in Table 4 and Fig. 8 are summarized below.
Antibiotics and metabolite producing circuits:
Presence of clusters of nonpathogen-linked orthologs in metabolic circuits such as steroid biosynthesis, arginine and proline metabolism, and the Insulin signaling pathway indicate the importance of these pathways in establishing synergy with the host in all the major genera considered (Table 4). Some of the orthologs in the nonpathogen-linked bacterial gene circuits have orthologs in the human. Other synergy circuits in Table 4 are involved in radiation survival . KEGG reference metabolic pathways contain large numbers of nodes creating thousands of clusters for bacterial species, and hence the relative lack of literature for some of the clusters shown in Table 4.
The polyketide circuit shown in Fig. 8 and presented in Table 4 facilitates the synthesis of common antibiotics . Polyketides are complex organic compounds, which are highly active biologically. Many pharmaceuticals are derived from or inspired by polyketides. In addition to the polyketide circuits, the circuit shown in Fig. 8 composed of scyllo-inosamine orthologs (stsE, strB1, E126.96.36.199, strK, K12570, aphD, strA) is involved in the biosynthesis of streptomycin and similar anti-mycobacteria antibiotics. Recent studies show bacterial virulence factors in type 3-secretion pathway as targeted by virulence inhibitors such as those illustrated in Fig. 8 . Also shown in the figure is the one-carbon pool by folate pathway, activating one-carbon units for biosynthesis . It plays a major role in amino acid metabolism . It has been shown to affect proof reading of DNA replication, DNA methylation, and chromatin structure [72–74]. The list for commensal circuitry presented in Table 4 is not complete, but representative of the diversity of commensal circuits found in bacterial strains.
Genera-specific genome comparisons reveal additional circuitry clusters found almost exclusively in nonpathogenic bacterial strains. Shown in Additional file 5 in rows 22 to 24 are clusters for benzoate degradation, and the cluster for dopamine circuitry in Isoquinoline alkaloid biosynthesis. Benzoate degradation is an important factor in reducing drug-induced toxicity . It is not clear how dopamine inducing bacterial gene circuits drive synergy with the host, yet modulations in dopamine circuitry in bacteria was previously linked to Alzheimer’s disease stage progressions via Borrelia infection .
Pathogenicity is context dependent. A survey of the literature indicates many dimensions of complexity in defining and investigating pathogenic behavior of bacterial strains and microbiomes in relation to host [75–78]. Consider a microbiome composed of thousands of bacterial strains. The combined genome is large, containing millions of genes. It is clear in this scenario the need for obtaining a signature for pathogenic or commensal phenotypes for dimension reduction. This study begins with such a task, namely deducing a list of genes which are highly abundant in pathogenic and relatively absent in commensal strains and vice versa. Two distinct datasets were used in our analysis: (a) list of bacterial strains observed to be pathogenic at least in a context dependent manner and (b) bacterial genomes annotated with ortholog genes. The literature search yielded 949 decoded bacterial genomes deemed as pathogenic and another 1578 as nonpathogenic or commensal due to lack of evidence in the literature for pathogenicity. Additional file 1, presenting the list of pathogenic strains along with literature-curated evidence for each, will prove to be a useful resource for microbiology.
The ortholog contents of genomes in KEGG Orthology (KO) Database were recently used in developing an algorithm (PICRUSt) for analysis of genes in microbiomes in health and disease . In the present case, KO database yielded two large matrices with columns identifying pathogenic or nonpathogenic strains and rows representing 7194 orthologs present in at least one bacterial strain. We created similar matrices for within genera comparison of pathogenic and nonpathogenic strains. From this point on, we could compare the two matrices and identify, for each ortholog, the relative abundance in pathogenic and nonpathogenic strains. The resulting ortholog based files compiled in Additional file 2 is easy to use, as for each ortholog, we provide its NCBI identified symbol, full ortholog name, and KEGG ID. In essence, the focus on orthology enabled genome comparison at a meta-scale, and enabled us to integrate discrete pieces of data in the literature into a system wide portrait.
Statistical enrichment processes we employed allowed us to investigate antibiotic resistance within the context of pathogenicity. Specifically, we considered pathway enrichment of orthologs present in 268 antibiotic resistant bacterial strains with respect to the 681 pathogenic strains with no documented antibiotic resistance. Some of the orthologs associated in the literature with antibiotic resistance  were indeed preferentially abundant in antibiotic resistant strains. Others were not but they co-localized with orthologs highly enriched in antibiotic resistant strains in cellular pathways. The pathways enriched by orthologs abundant in antibiotic resistant strains included P13K-Akt signaling pathway of eukaryotic hosts. The effector bacterial orthologs in this pathway divert host cell signaling pathways to the benefit of the pathogen and target kinase signaling cascades present in P13K-Akt , resulting in pathogenic infections . The other most enriched pathway in antibiotics resistant bacerial strains was peptidoglycan biosynthesis, previously linked to biofilm production . The pathway mediating the metabolism of Sphingolipids, a class of lipids, was occupied by orthologs found in antibiotic resistant bacteria. Sphingolipids play significant roles in membrane biology and provide many bioactive metabolites that regulate cell function . It has already been linked in the literature to antibiotic resistance [85, 86].
Pathway enrichment protocols provided us with an overall portrait of ortholog sets linked either to commensal or pathogenic behavior to the host. Additional annotation was needed however in order to get a signature for synergy or pathogenicity. Visually, when we projected pathogen- or synergy-linked ortholog lists onto KEGG pathways, we could see orthologs from our lists forming multiply connected clusters along the bacterial pathways. Thus, our next task was identifying circuits along the pathways that contained such clusters and then to conduct extensive literature searches to annotate our discoveries.
The circuits identified included toxin-producing units for pathogenic strains and the gene circuits synthesizing antibiotics for commensal strains. We observed that the list of circuitry linked to synergy contained among others, biosynthesis modules of type II polyketide products, such as Erythromycin, and Doxycycline. Other examples of nonpathogen-linked ortholog clusters in gene circuits could be found in VEGF signaling pathway and in metabolic pathways.
The ortholog modules highly abundant in pathogenic bacteria included gene circuits found in bacterial secretion mechanism type III. Also in this category were gene circuits involved in peptidoglycan biosynthesis. In general, pathogen abundant orthologs were statistically enriched in KEGG pathways for pathogen interactions and in signaling pathways. Comparison of pathogenic strains with nonpathogenic strains of the same genera identified one carbon pool by folate pathway as highly abundant in pathogenic strains. This pathway mediates the activation of one-carbon units in the biosynthesis and metabolism of amino acids. It affects DNA methylation as well as DNA replication. Recent research implicates it in folic acid-mediated degeneration of the brain–blood barrier of the host . In addition, our system approach identified a bacterial gene circuit whose genes are preferentially found in pathogenic strains and linked in the literature to Alzheimer’s disease. We find the gene clusters in this circuit to be abundant in a diverse set of genera, possibly providing new avenues of research for linkage between bacterial infections and Alzheimer’s disease.
The work presented here is a first draft of an encyclopedia for bacteria – host interactions. The study uncovers identity of ortholog clusters as possible signatures of pathogenicity or synergy in a mixture of bacteria. We show that gene clusters formed on bacterial pathways to be much stronger determinants of pathogenicity than a list of virulence and/or anti-virulence factors. The signature we derive in the form of gene circuits may not be complete as our results relied on the decoded genomes of bacterial strains currently available in the literature. Moreover, our reliance on the KEGG reference pathways in identifying clusters of orthologs preferentially found in pathogenic and nonpathogenic bacteria introduces additional limitations due to incompleteness of the KEGG pathway models. Nevertheless, the gene circuitry signature we discovered for synergy and pathogenicity is comprehensive enough to derive from it biomarker clusters identifying pathogenic phenotypes in bacterial strains isolated or in a mixture.
This study presents a system approach for identifying gene clusters either preferentially present or absent in pathogenic bacterial strains. The study utilized 2527 fully sequenced bacterial strain genomes available in the public domain. Literature search identified 949 of these genomes to belong to strains with demonstrated pathogenic potential. Ortholog abundance comparisons between pathogenic and nonpathogenic strains within and across genera revealed signatures for pathogenic and commensal bacteria. Known virulence factors were highly enriched in the genomes of pathogenic strains.
Projection of ortholog gene signatures onto cellular pathways revealed gene circuits linked either to synergy or pathogenicity. The pathogenicity related gene circuits included those in bacterial two-component system, biosynthesis of siderophores, and one-carbon pool by folate pathways. Circuits belonging to shingolipid metabolism, P13-Akt pathway, and tRNA synthesis were particularly enriched with orthologs preferentially expressed in antibiotic-resistant bacteria. Genes preferentially expressed by nonpathogenic bacteria also formed circuits, among which were those linked to the synthesis of antibiotics. The study comprises an important step towards addressing crosstalk between host, virus, prokaryotes, and environment using a system approach [88–90].
Annotation of bacterial strains exhibiting pathogenicity
The list of pathogenic bacterial strains was obtained via literature search using three complementary and partially redundant approaches. First, we text mined the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [91, 92] for pathogenicity label for all decoded bacterial strains in KEGG. Second, we identified as pathogenic those bacterial strains with pathogenic citation in at least one of the following web tools: Virulence Factor Database (VFDB) [2, 6, 9], High-quality Automated and Manual Annotation of Proteins (HAMAP) , and the Interactive Atlas for Exploring Bacterial Genomes (BacMap) . Third, we text mined PUBMED article abstracts with a) the name of the bacterial strain and b) one of the pathogen-related code words (pathogen, virulence, pathogenic, virulent). All hits were then verified for pathogenicity by reading the articles. We labeled bacterial strains as nonpathogenic if the process described above did not yield evidence for pathogenesis. We have used the Antibiotic Resistance Database (ABRD)  to identify antibiotic resistant bacterial strains within the set of pathogenic bacterial strains.
Orthology content of decoded bacterial strains
KEGG programming interface was used to obtain gene orthology information for the bacterial organisms found in KEGG. The individual organism orthology content was translated into a logical content vector that describes the orthology information for that particular organism in binary form, “1” meaning that the ortholog is present and “0” meaning it is absent, with respect to all collective orthologs found in the KEGG repository. The individual feature vectors were than accumulated into a logic matrix, with rows representing the presence/absence of orthologs for each decoded bacterial strain. The resulting array contained 2527 sequenced bacterial strain genomes potentially expressing 7194 orthologs. The number of genes that were accounted for in this array totaled to over 3.2 million. The logical array was then split into pathogenic vs. non-pathogenic logical arrays with equal number of columns representing the orthologs.
Abundance scores for orthologs in bacterial strain genomes
An abundance score for each ortholog was generated for both the pathogenic and nonpathogenic bacterial strain arrays based on the sum of all elements in the ortholog column vector divided by vector length (number of bacterial strains). These scores were denoted as Ap and Anp for pathogenic and nonpathogenic strains, respectively. They represent the fraction of genomes expressing the ortholog in pathogenic and nonpathogenic logical matrices. Next we defined a pathogenicity abundance score as PA = Ap / (Anp + 0.0001). If Anp turned out to be “0” for that ortholog, meaning that ortholog was absent in nonpathogenic bacteria, the equation yields PA as equal to Ap / 10−4. The list of orthologs could then be ordered with respect to the PA value, creating a histogram.
The comparisons of contents of genome orthology of bacterial strains belonging to the same genera are identical to the one for all the 2527 content comparison described above. For this operation, we used only strains belonging to the genera shown in Table 2, one genus at a time. Also, in this case we reduced the cutoffs for pathogen association and defined orthologs as pathogen-abundant for PA > 2, and nonpathogen abundant for PA < 1/2. This was necessitated by the relative similarity of genomes within a genus. However, to compensate, we studied only those clusters in KEGG Reference pathways abundant in pathogen or nonpathogen strains in multiple genera.
Statistical enrichment of cellular pathways with pathogenic and nonpathogen-abundant orthologs
We identified orthologs as pathogen exclusive if they were expressed in some of the pathogenic strains but not at all in nonpathogenic strains. Nonpathogen exclusive orthologs were similarly defined. We defined orthologs as pathogen-abundant for PA > 4, and nonpathogen abundant for PA < 1/4. These were also called pathogenic orthologs and vice versa.
For KEGG cellular pathway enrichments, we created a score matrix, similar to that described in the making of the orthology database. Next, the enrichment was carried out via hypergeometric test [63, 95] using total number of orthologs in bacterial strains as the population size, pathogenic orthologs as the number of success states in the population; number of orthologs in the pathways as the number of draws; and the number of pathogenic orthologs in the pathways as the number of successes in the draw.
Subpathways containing pathogenic ortholog circuits were identified via manual KEGG orthology mapping and screening for clustering. Cellular pathways were drawn using the KEGG web tool for pathways. Although the cutoff values for PA for generating pathogenic and nonpathogenic ortholog lists appear arbitrary, they capture the tails of the PA distribution for orthologs in the thousands of decoded genomes under consideration. Preliminary studies involving the perturbation of the cutoff did not yield variation in the enriched pathways.
Ortholog clusters relevant to pathogenicity in KEGG reference cellular pathways
For this purpose, we mapped the list of pathogen-abundant and nonpathogen abundant orthologs to all available KEGG Reference pathways, and manually curated the clusters (sub circuitry) they form. The criterion for cluster was the minimum of three orthologs from the list to be multiply-connected with each other and have connections to other neighboring orthologs. Results shown in Table 3 to Additional file 5 indicate only the gene symbols of those abundant in pathogenic or vice versa. As presented in these tables, most clusters contained many more orthologs then three. However, the Reference pathways typically contain genes not found in all strains and hence in almost all cases we studied 75 % presence of pathogene-abundant orthologs in a strain as indicating that the strain was deemed pathogenic in the literature. The p values listed in these tables were obtained using hypergeometric test based on the following scheme: a) there are 949 strains deemed pathogenic in a pool of 2527 strains, and b) among M strains, m cluster orthologs were found only on N strains. In this set up, M, N, and m are positive integers.
Antibiotics resistance genes database
Kyoto encyclopedia of genes and genomes
Vascular endothelial growth factor
Virulence factor database
Rahme LG et al. Plants and animals share functionally common bacterial virulence factors. Proc Natl Acad Sci U S A. 2000;97(16):8815–21.
Chen L, Xiong Z, Sun L, Yang J, Jin Q. VFDB 2012 update: toward the genetic diversity and molecular evolution of bacterial virulence factors. Nucleic Acids Res. 2012;40(Database issue):D641–5.
Long TM, Nisa S, Donnenberg MS, Hassel BA. Enteropathogenic Escherichia coli inhibits type I interferon- and RNase-L-mediated host defense to disrupt intestinal epithelial cell barrier function. Infect Immun. 2014;82(7):2802–14.
Jorth P, Turner KH, Gumus P, Nizam N, Buduneli N, Whiteley M. Metatranscriptomics of the human oral microbiome during health and disease. MBio. 2014;5(2):e01012–4.
Seyedsayamdost MR. High-throughput platform for the discovery of elicitors of silent bacterial gene clusters. Proc Natl Acad Sci U S A. 2014;111(20):7266–71.
Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y et al. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2005;33(Database issue):D325–8.
Passador L, Cook JM, Gambello MJ, Rust L, Iglewski BH. Expression of Pseudomonas aeruginosa virulence genes requires cell-to-cell communication. Science. 1993;260(5111):1127–30.
de Kievit TR, Kakai Y, Register JK, Pesci EC, Iglewski BH. Role of the Pseudomonas aeruginosa las and rhl quorum-sensing systems in rhlI regulation. FEMS Microbiol Lett. 2002;212(1):101–6.
Yang J, Chen L, Yu J, Jin Q. VFDB 2008 release: an enhanced web-based resource for comparative pathogenomics. Nucleic Acids Res. 2008;36(Database issue):D539–42.
Niu C, Yu D, Wang Y, Ren H, Jin Y, Zhou W et al. Common and pathogen-specific virulence factors are different in function and structure. Virulence. 2013;4(6):473–82.
Marvig RL, Johansen HK, Molin S, Jelsbak L. Genome analysis of a transmissible lineage of pseudomonas aeruginosa reveals pathoadaptive mutations and distinct evolutionary paths of hypermutators. PLoS Genet. 2013;9(9):e1003741.
Koskiniemi S, Gibbons HS, Sandegren L, Anwar N, Ouellette G, Broomall S et al. Pathoadaptive mutations in Salmonella enterica isolated after serial passage in mice. PLoS One. 2013;8(7):e70147.
Prosseda G, Di Martino ML, Campilongo R, Fioravanti R, Micheli G, Casalino M et al. Shedding of genes that interfere with the pathogenic lifestyle: the Shigella model. Res Microbiol. 2012;163(6–7):399–406.
Juhász J, Kertész-farkas A, Szabó D, Pongor S. Emergence of collective territorial defense in bacterial communities: horizontal gene transfer can stabilize microbiomes. PLoS One. 2014;9(4):e95511.
Maurelli AT. Black holes, antivirulence genes, and gene inactivation in the evolution of bacterial pathogens. FEMS Microbiol Lett. 2007;267(1):1–8.
Swan BK, Chaffin MD, Martinez-garcia M, Morrison HG, Field EK, Poulton NJ et al. Genomic and metabolic diversity of marine group I thaumarchaeota in the mesopelagic of two subtropical gyres. PLoS One. 2014;9(4):e95380.
Yang X, Brisbin J, Yu H, Wang Q, Yin F, Zhang Y et al. Selected lactic Acid-producing bacterial isolates with the capacity to reduce salmonella translocation and virulence gene expression in chickens. PLoS One. 2014;9(4):e93022.
Miskinyte M, Sousa A, Ramiro RS, de Sousa JA, Kotlinowski J, Caramalho I et al. The genetic basis of Escherichia coli pathoadaptation to macrophages. PLoS Pathog. 2013;9(12):e1003802.
Podoll JD, Liu Y, Chang L, Walls S, Wang W, Wang X. Bio-inspired synthesis yields a tricyclic indoline that selectively resensitizes methicillin-resistant Staphylococcus aureus (MRSA) to beta-lactam antibiotics. Proc Natl Acad Sci U S A. 2013;110(39):15573–8.
Uda N, Matoba Y, Kumagai T, Oda K, Noda M, Sugiyama M. Establishment of an in vitro D-cycloserine-synthesizing system by using O-ureido-L-serine synthase and D-cycloserine synthetase found in the biosynthetic pathway. Antimicrob Agents Chemother. 2013;57(6):2603–12.
Paredes D, Ortiz C, Torres R. Synthesis, characterization, and evaluation of antibacterial effect of Ag nanoparticles against Escherichia coli O157:H7 and methicillin-resistant Staphylococcus aureus (MRSA). Int J Nanomedicine. 2014;9:1717–29.
Ghosh S, Prava J, Samal HB, Suar M, Mahapatra RK. Comparative genomics study for the identification of drug and vaccine targets in Staphylococcus aureus: MurA ligase enzyme as a proposed candidate. J Microbiol Methods. 2014;101:1–8.
Machowski EE, Senzani S, Ealand C, Kana BD. Comparative genomics for mycobacterial peptidoglycan remodelling enzymes reveals extensive genetic multiplicity. BMC Microbiol. 2014;14:75.
Wichmann F, Udikovic-kolic N, Andrew S, Handelsman J. Diverse antibiotic resistance genes in dairy cow manure. MBio. 2014;5(2):e01017.
Liu B, Pop M. ARDB–Antibiotic Resistance Genes Database. Nucleic Acids Res. 2009;37(Database issue):D443–7.
Gouaze-andersson V, Cabot MC. Sphingolipid metabolism and drug resistance in hematological malignancies. Anticancer Agents Med Chem. 2011;11(9):891–903.
West KA, Castillo SS, Dennis PA. Activation of the PI3K/Akt pathway and chemotherapeutic resistance. Drug Resist Updat. 2002;5(6):234–48.
Hurdle JG, O’neill AJ, Chopra I. Prospects for aminoacyl-tRNA synthetase inhibitors as new antimicrobial agents. Antimicrob Agents Chemother. 2005;49(12):4821–33.
Macho AP, Zipfel C. Targeting of plant pattern recognition receptor-triggered immunity by bacterial type-III secretion system effectors. Curr Opin Microbiol. 2014;23C:14–22.
Dai W, Li Z. Conserved type III secretion system exerts important roles in Chlamydia trachomatis. Int J Clin Exp Pathol. 2014;7(9):5404–14.
Patel S, Mccormick BA. Mucosal Inflammatory Response to Salmonella typhimurium Infection. Front Immunol. 2014;5:311.
Galán JE, Lara-tejero M, Marlovits TC, Wagner S. Bacterial type III secretion systems: specialized nanomachines for protein delivery into target cells. Annu Rev Microbiol. 2014;68:415–38.
Christie PJ, Whitaker N, González-rivera C. Mechanism and structure of the bacterial type IV secretion systems. Biochim Biophys Acta. 2014;1843(8):1578–91.
Wang Y, Wei X, Bao H, Liu SL. Prediction of bacterial type IV secreted effectors by C-terminal features. BMC Genomics. 2014;15:50.
Cabezón E, Ripoll-rozada J, Peña A, De La Cruz F, Arechaga I. Towards an integrated model of bacterial conjugation. FEMS Microbiol Rev. 2014;39(1):81–95.
Stephenson K, Hoch JA. Two-component and phosphorelay signal-transduction systems as therapeutic targets. Curr Opin Pharmacol. 2002;2(5):507–12.
Bakhtiar SM, Leblanc JG, Salvucci E, Ali A, Martin R, Langella P et al. Implications of the human microbiome in inflammatory bowel diseases. FEMS Microbiol Lett. 2013;342(1):10–7.
Michail S, Durbin M, Turner D, Griffiths AM, Mack DR, Hyams J et al. Alterations in the gut microbiome of children with severe ulcerative colitis. Inflamm Bowel Dis. 2012;18(10):1799–808.
Honda K, Littman DR. The microbiome in infectious disease and inflammation. Annu Rev Immunol. 2012;30:759–95.
Weiser JN, Lindberg AA, Manning EJ, Hansen EJ, Moxon ER. Identification of a chromosomal locus for expression of lipopolysaccharide epitopes in Haemophilus influenzae. Infect Immun. 1989;57(10):3045–52.
Neilands JB. Siderophores: structure and function of microbial iron transport compounds. J Biol Chem. 1995;270(45):26723–6.
Chen T, He L, Deng W, Xie J. The Mycobacterium DosR regulon structure and diversity revealed by comparative genomic analysis. J Cell Biochem. 2013;114(1):1–6.
Saini DK, Malhotra V, Dey D, Pant N, Das TK, Tyagi JS. DevR-DevS is a bona fide two-component system of Mycobacterium tuberculosis that is hypoxia-responsive in the absence of the DNA-binding domain of DevR. Microbiology (Reading, Engl). 2004;150(Pt 4):865–75.
Psakis G, Saidijam M, Shibayama K, Polaczek J, Bettaney KE, Baldwin JM et al. The sodium-dependent D-glucose transport protein of Helicobacter pylori. Mol Microbiol. 2009;71(2):391–403.
Choi E, Lee KY, Shin D. The MgtR regulatory peptide negatively controls expression of the MgtA Mg2+ transporter in Salmonella enterica serovar Typhimurium. Biochem Biophys Res Commun. 2012;417(1):318–23.
Yan Q, Gao W, Wu XG, Zhang LQ. Regulation of the PcoI/PcoR quorum-sensing system in Pseudomonas fluorescens 2P24 by the PhoP/PhoQ two-component system. Microbiology (Reading, Engl). 2009;155(Pt 1):124–33.
Walters SB, Dubnau E, Kolesnikova I, Laval F, Daffe M, Smith I. The Mycobacterium tuberculosis PhoPR two-component system regulates genes essential for virulence and complex lipid biosynthesis. Mol Microbiol. 2006;60(2):312–30.
Rokosz-chudziak N, Rastawicki W. Selected mechanisms of pathogenicity of Campylobacter jejuni. Med Dosw Mikrobiol. 2014;66(1):47–58.
Agarwal S, Vasudhev S, Deoliveira RB, Ram S. Inhibition of the classical pathway of complement by meningococcal capsular polysaccharides. J Immunol. 2014;193(4):1855–63.
Tan Y, Kagan JC. A cross-disciplinary perspective on the innate immune responses to bacterial lipopolysaccharide. Mol Cell. 2014;54(2):212–23.
Boneca IG. The role of peptidoglycan in pathogenesis. Curr Opin Microbiol. 2005;8(1):46–53.
Hassan SS, Tiwari S, Guimarães LC, Jamal SB, Folador E, Sharma NB et al. Proteome scale comparative modeling for conserved drug and vaccine targets identification in Corynebacterium pseudotuberculosis. BMC Genomics. 2014;15 Suppl 7:S3.
Landau GM, Parida L, Weimann O. Gene proximity analysis across whole genomes via PQ trees. J Comput Biol. 2005;12(10):1289-306.
Molinari G, Chhatwal GS. Streptococcal invasion. Curr Opin Microbiol. 1999;2(1):56–61.
Miethke M, Marahiel MA. Siderophore-based iron acquisition and pathogen control. Microbiol Mol Biol Rev. 2007;71(3):413–51.
Anderson LL. Discovery of the ‘porosome’; the universal secretory machinery in cells. J Cell Mol Med. 2006;10(1):126–31.
Galán JE, Wolf-watz H. Protein delivery into eukaryotic cells by type III secretion machines. Nature. 2006;444(7119):567–73.
Li J, Lv C, Sun W, Li Z, Han X, Li Y et al. Cytosporone B, an inhibitor of the type III secretion system of Salmonella enterica serovar Typhimurium. Antimicrob Agents Chemother. 2013;57(5):2191–8.
Nirmala J, Drader T, Lawrence PK, Yin C, Hulbert S, Steber CM et al. Concerted action of two avirulent spore effectors activates Reaction to Puccinia graminis 1 (Rpg1)-mediated cereal stem rust resistance. Proc Natl Acad Sci U S A. 2011;108(35):14676–81.
Lemos EG, Alves LM, Campanharo JC. Genomics-based design of defined growth media for the plant pathogen Xylella fastidiosa. FEMS Microbiol Lett. 2003;219(1):39–45.
Williams MJ, Shanley CA, Zilavy A, Peixoto B, Manca C, Kaplan G et al. bis-Molybdopterin guanine dinucleotide is required for persistence of Mycobacterium tuberculosis in guinea pigs. Infect Immun. 2014;83(2):544–50.
Tang Q, Li X, Zou T, Zhang H, Wang Y, Gao R et al. Mycobacterium smegmatis BioQ defines a new regulatory network for biotin metabolism. Mol Microbiol. 2014; doi: 10.1111/mmi.12817
Kountouras J, Gavalas E, Polyzos SA, Deretzi G, Kouklakis G, Grigoriadis S et al. Association between Helicobacter pylori burden and Alzheimer’s disease. Eur J Neurol. 2014;21(12):e100.
Bibi F, Yasir M, Sohrab SS, Azhar EI, Al-Qahtani MH, Abuzenadah AM et al. Link between chronic bacterial inflammation and Alzheimer disease. CNS Neurol Disord: Drug Targets. 2014;13(7):1140–7.
Maheshwari P, Eslick GD. Bacterial infection and Alzheimer’s disease: a meta-analysis. J Alzheimers Dis. 2015;43(3):957–66.
Nuyts S, Van Mellaert L, Theys J, Landuyt W, Lambin P, Anné J. The use of radiation-induced bacterial promoters in anaerobic conditions: a means to control gene expression in clostridium-mediated therapy for cancer. Radiat Res. 2001;155(5):716–23.
Isoda K, Nozawa T, Tezuka M, Ishida I. Toxicity of 50-nm polystyrene particles co-administered to mice with acetaminophen, 5-aminosalicylic acid or tetracycline. Pharmazie. 2014;69(9):676–9.
Macdonald AB. Alzheimer’s disease Braak Stage progressions: reexamined and redefined as Borrelia infection transmission through neural circuits. Med Hypotheses. 2007;68(5):1059–64.
Aiello D, Williams JD, Majgier-Baranowska H, Patel I, Peet NP, Huang J et al. Discovery and characterization of inhibitors of Pseudomonas aeruginosa type III secretion. Antimicrob Agents Chemother. 2010;54(5):1988–99.
Fox JT, Stover PJ. Folate-mediated one-carbon metabolism. Vitam Horm. 2008;79:1–44.
Stover PJ. One-carbon metabolism-genome interactions in folate-associated pathologies. J Nutr. 2009;139(12):2402–5.
Kamen B. Folate and antifolate pharmacology. Semin Oncol. 1997;24(5 Suppl 18):S18–30. S18-39.
Crider KS, Yang TP, Berry RJ, Bailey LB. Folate and DNA methylation: a review of molecular mechanisms and the evidence for folate’s role. Adv Nutr. 2012;3(1):21–38.
Zhou HR, Zhang FF, Ma ZY, Huang HW, Jiang L, Cai T et al. Folate polyglutamylation is involved in chromatin silencing by maintaining global DNA methylation and histone H3K9 dimethylation in Arabidopsis. Plant Cell. 2013;25(7):2545–59.
Lijek RS, Weiser JN. Co-infection subverts mucosal immunity in the upper respiratory tract. Current Opin Immunol. 2012;24:417–23.
Branchu P, Matrat S, Vareille M, et al. NsrR, GadE, and GadX interplay in repressing expression of the Escherichia coli O157:H7 LEE pathogenicity island in response to nitric oxide. PLoS Pathog. 2014;10(1):e1003874.
Marullo S, Coureuil M. Arrestins in host-pathogen interactions. Handb Exp Pharmacol. 2014;219:361–74.
Tack AJ, Horns F, Laine AL. The impact of spatial scale and habitat configuration on patterns of trait variation and local adaptation in a wild plant parasite. Evolution. 2014;68(1):176–89.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.
Nesme J, Cécillon S, Delmont TO, Monier JM, Vogel TM, Simonet P. Large-scale metagonomic-based study of antibiotic resistance in the environment. Current Biol. 2014;24:1–5.
Krachler AM, Woolery AR, Orth K. Manipulation of kinase signaling by bacterial pathogens. J Cell Biol. 2011;195(7):1083–92.
Oviedo-boyso J, Cortés-vieyra R, Huante-mendoza A, Yu HB, Valdez-Alarcón JJ, Bravo-Patiño A et al. The phosphoinositide-3-kinase-Akt signaling pathway is important for Staphylococcus aureus internalization by endothelial cells. Infect Immun. 2011;79(11):4569–77.
Li L, Zhu J, Yang K, Xu Z, Liu Z, Zhou R. Changes in gene expression of Actinobacillus pleuropneumoniae in response to anaerobic stress reveal induction of central metabolism and biofilm formation. J Microbiol. 2014;52(6):473–81.
Gault CR, Obeid LM, Hannun YA. An overview of sphingolipid metabolism: from synthesis to breakdown. Adv Exp Med Biol. 2010;688:1–23.
Morad SA, Cabot MC. Tamoxifen regulation of sphingolipid metabolism--Therapeutic implications. Biochim Biophys Acta. 2015;1851(9):1134–45.
Hill JA, Ammar R, Torti D, Nislow C, Cowen LE. Genetic and genomic architecture of the evolution of resistance to antifungal drug combinations. PLoS Genet. 2013;9(4):e1003390.
Kalani A, Kamat PK, Givvimani S, Brown K, Metreveli N, Tyagi SC et al. Nutri-epigenetics ameliorates blood–brain barrier damage and neurodegeneration in hyperhomocysteinemia: role of folic acid. J Mol Neurosci. 2014;52(2):202–15.
Freeman S, Bhatt A, Pedamallu C, King S, Duke F, Jung J et al. A121: in search of infectious triggers of periodic Fever, aphthous stomatitis, pharyngitis and adenitis syndrome. Arthritis Rheumatol. 2014;66(11):S158.
Xu G, Strong MJ, Lacey MR, Baribault C, Flemington EK, Taylor CM. RNA CoMPASS: a dual approach for pathogen and host transcriptome analysis of RNA-seq datasets. PLoS One. 2014;9(2):e89445.
Minot S, Bryson A, Chehoud C, Wu GD, Lewis JD, Bushman FD. Rapid evolution of the human gut virome. Proc Natl Acad Sci U S A. 2013;110(30):12450–5.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28:27–30.
Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42:D199–205.
Pedruzzi I, Rivoire C, Auchincloss AH, Coudert E, Keller G, de Castro E et al. HAMAP in 2013, new developments in the protein family classification and annotation system. Nucleic Acids Res. 2013;41(Database issue):D584–9.
Cruz J, Liu Y, Liang Y, Zhou Y, Wilson M, Dennis JJ, et al. BacMap: an up-to-date electronic atlas of annotated bacterial genomes. Nucleic Acids Res. 2012;40(Database issue):D599–604.
Rivals I, Personnaz L, Taing L, Potier MC. Enrichment or depletion of a GO category within a class of genes: which test? Bioinformatics. 2007;23(4):401–7.
MS was supported by the US Department of Education Doctoral Training Fellowship Program (GAANN), entitled Bioinformatics of Infectious Diseases. AT is the PI and the Director of the GAANN project.
This study utilized open access data, exclusively. No human or animal subjects were involved. Supporting files provide all the information required to reproduce our results.
The authors declare that they have no competing interests.
AT and MS conceptualized and designed the study. MS wrote the algorithms for the study and implemented the research with continuous input and discussions with AT and SO. AT and MS co-wrote the manuscript. All authors read and approved the final manuscript.
Availability of Data and Materials
List of 949 decoded bacterial genomes present in the KEGG web platform deemed as pathogenic along with literature citations for pathogenicity and antibiotic resistance phenotype. (XLSX 51 kb)
List of 1578 decoded bacterial genomes present in the KEGG web platform with no present citation implicating them as pathogenic. (XLSX 54 kb)
Bacterial strain orthology logic array indicating the presence of given orthologs (7194 of them) in the 2527 bacterial strain genomes considered in this study. (RAR 1397 kb)
List of 7194 orthologs ordered according to pathogen association score PA. The orthologs shaded in pink represent ones present in VFDB. (XLSX 570 kb)
Pathogen- and nonpathogen gene circuits identified via within genera comparisons. The genera exhibiting a pathogen-linked circuit are painted in pink whereas nonpathogen-linked circuits are painted in green.
About this article
Cite this article
Shestov, M., Ontañón, S. & Tozeren, A. Encyclopedia of bacterial gene circuits whose presence or absence correlate with pathogenicity – a large-scale system analysis of decoded bacterial genomes. BMC Genomics 16, 773 (2015). https://doi.org/10.1186/s12864-015-1957-7
- Gene circuits
- Virulence factors