- Research article
- Open Access
A functional genomics catalogue of activated transcription factors during pathogenesis of pneumococcal disease
- Layla K Mahdi†1,
- Tahereh Deihimi†2,
- Fatemeh Zamansani†2,
- Mario Fruzangohar3,
- David L Adelson4,
- James C Paton1,
- Abiodun D Ogunniyi1Email author and
- Esmaeil Ebrahimie5
© Mahdi et al.; licensee BioMed Central Ltd. 2014
- Received: 25 July 2014
- Accepted: 3 September 2014
- Published: 8 September 2014
Streptococcus pneumoniae (the pneumococcus) is the world’s foremost microbial pathogen, killing more people each year than HIV, TB or malaria. The capacity to penetrate deeper host tissues contributes substantially to the ability of this organism to cause disease. Here we investigated, for the first time, functional genomics modulation of 3 pneumococcal strains (serotype 2 [D39], serotype 4 [WCH43] and serotype 6A [WCH16]) during transition from the nasopharynx to lungs to blood and to brain of mice at both promoter and domain activation levels.
We found 7 highly activated transcription factors (TFs) [argR, codY, hup, rpoD, rr02, scrR and smrC] capable of binding to a large number of up-regulated genes, potentially constituting the regulatory backbone of pneumococcal pathogenesis. Strain D39 showed a distinct profile in employing a large number of TFs during blood infection. Interestingly, the same highly activated TFs used by D39 in blood are also used by WCH16 and WCH43 during brain infection. This indicates that different pneumococcal strains might activate a similar set of TFs and regulatory elements depending on the final site of infection. Hierarchical clustering analysis showed that all the highly activated TFs, except rpoD, clustered together with a high level of similarity in all 3 strains, which might suggest redundancy in the regulatory roles of these TFs during infection. Discriminant function analysis of the TFs in various niches highlights differential regulatory backgrounds of the 3 strains, and pathogenesis data confirms codY as the most significant predictor discriminating between these strains in various niches, particularly in the blood. Moreover, the predicted TF and domain activation profiles of the 3 strains correspond with their distinct pathogenicity characteristics.
Our findings suggest that the pneumococcus changes the short binding sites in the promoter regions of genes in a niche-specific manner to enhance its ability to disseminate from one host niche to another. This study provides a framework for an improved understanding of the dynamics of pneumococcal pathogenesis, and opens a new avenue into similar investigations in other pathogenic bacteria.
- Discriminant Function Analysis
- Pneumococcal Disease
- Pneumococcal Strain
- Host Niche
- Pneumococcal Virulence
Streptococcus pneumoniae (the pneumococcus) continues to cause high morbidity and mortality worldwide, in spite of the availability of vaccines and antimicrobial therapies [1, 2]. Resistance of virulent S. pneumoniae to multiple antibiotics, particularly against beta-lactams is due to alterations in the structure of six penicillin-binding proteins, while macrolide resistance is mediated through the erm(B), mefA or mefE genes . This could partly be explained by the fact that the organism deploys an efficient virulence machinery during infection. Consequently, research efforts have been geared towards understanding the molecular mechanisms underlying the pathogenesis of pneumococcal disease [4–6]. However, due to technical difficulties associated with harvesting RNA for analysis of in vivo transcription patterns of pneumococci during pathogenesis, most studies have either studied the earlier stages of infection (nasopharynx and lungs) or used in vitro surrogates [7–10]. Recently, these challenges have been largely overcome by studies involving transcriptomic comparisons of differentially-regulated genes during penetration of deeper host tissues [11–13].
The study of activated domains of expressed proteins in addition to promoter activation profiles can result in a better understanding of functional genomics. Compared to TFs, more information is available on pneumococcal pathogenesis at the level of expressed protein domains because of its use in vaccine design. The role of ATP-dependent transport and DNA binding domains in pneumococcal competence , histidine kinase in sensing environmental cues and virulence , and the role of PsaA in virulence [18, 28], are well documented. Genome wide analysis of histidine kinases and response regulators has led to the discovery of new response regulator loci , while mutation in the DNA gyrase domain has resulted in resistance to quinolones .
In the present study, we used functional genomics tools to examine the underlying molecular mechanisms that underpin the different pathogenicity characteristics of 3 pneumococcal strains at both promoter (transcriptomic) and domain (proteomic) activation levels. To our knowledge, this combined approach of characterizing in vivo microarray data has not been reported previously and can be used as a model for unraveling functional genomics attributes of other pathogenic bacteria.
For the purposes of this investigation, we have used 3 well-characterised pneumococcal strains: D39 (serotype 2), WCH43 (serotype 4) and WCH16 (serotype 6A). Mouse intranasal challenge experiments have shown that D39 and WCH43 are more virulent than WCH16. However, D39 causes severe pneumonia and high-grade bacteremia, while WCH16 and WCH43 have a propensity to translocate to the brain of infected mice [21, 30] (see Additional file 1: Figure S1). Nevertheless, WCH43 infection of mice demonstrates the “classical” disease progression from the nasopharynx to the lungs and dissemination to blood and then to the brain while WCH16 seems to progress directly to the brain with minimal lung and blood involvement [21, 31]. These striking differences in the pathogenicity and virulence characteristics of these strains make them ideal for comprehensive functional genomics analyses.
Pneumococcal strains activate a variety of TFs during pathogenesis
Our initial bioinformatic analysis of microarray data for D39, WCH16 and WCH43 at 72 h post-infection of mice suggests that WCH16 up-regulates many genes in the lungs and brain, WCH43 up-regulates many genes in the brain, whereas D39 has significant genome activation in the blood. The list of up-regulated genes of D39, WCH16 and WCH43 during pathogenesis is presented in Additional files 2, 3, 4: Tables S1-S3.
Transcription factors (TFs) with the highest number of binding sites on the promoter regions of up-regulated genes during pneumococcal pathogenesis
Lungs vs nasopharynx
Blood vs lungs
Brain vs blood
Activated TFs with the highest number of TFBs cluster together
Activated TFs exhibit different behavior between pneumococcal strains
We used discriminant function analysis to compare the relative impact (weight) of TFs between the 3 strains in various niches. Comparison of D39 and WCH16 showed that both strains have relatively similar discriminant models in the lungs. In the blood, a large number of TFs have opposite signs between D39 and WCH16. This highlights differential regulatory backgrounds of both strains in the blood (see Additional file 8: Table S7). However, the discriminant modelling algorithm indicate that D39 and WCH43 have similar regulatory backgrounds in the lungs but are different in the blood. The discriminant models of WCH16 and WCH43 in the lungs and brain show that the activated TFs have similar coefficient values and signs in both strains (Additional file 8: Table S7). Interestingly, these strains choose different models during progression from lungs to blood. In this instance, of the 7 highly activated TFs, codY was the most significant predictor discriminating WCH16 from WCH43 in the blood, where its coefficient was -38 in WCH16, compared to +13 in WCH43.
Given the predicted discriminant function for codY between these 3 strains, we hypothesized that a mutant of codY will show distinct pathogenic profiles between these strains in vivo. However, we noted that a previous manuscript reported that codY inactivation in D39 is difficult due to its essentiality  and suggested further work in other clinical strains to establish if this is a general feature of pneumococci. Therefore, we attempted to delete codY in D39, WCH16 and WCH43 by targeted mutagenesis using the overlap extension strategy, as described in Methods. Surprisingly, we obtained putative mutants in all 3 strains, which were verified to be correct by PCR and sequencing using primers flanking the codY ORF. We examined if the putative codY mutants carried the permissive amino acid (Q166 → Stop) mutation in FatC (SP_1870) and/or amino acids (D480 → Y; D487 → Y) and amino acids 1–81 deletion mutations in AmiC (SP_1890) reported previously . This was carried out by PCR and sequencing of the amiC and fatC genes of 2 independent mutants from each of the 3 strains. However, we could not detect any of these reported changes. Next, we tested the fitness of each of these putative codY mutants in an in vitro competition with the isogenic wild type over 4.5 hrs. Surprisingly, we found that the mutant of all 3 strains was completely out-competed by the wild type. Output ratios were 1:1300, 1:655 and 1:1900, for D39ΔcodY, WCH16ΔcodY and WCH43ΔcodY mutants, respectively at 3 h. We concluded that while the codY mutation might have been tolerated in these strains, the mutation is unstable, in agreement with the finding of Caymaris and colleagues regarding the essentiality of CodY . We also found that during in vitro growth, the antibiotic resistance selection marker in the codY mutant strains is lost in the absence of selection. This was verified by replica-plating cells grown on antibiotic plates onto plates with or without antibiotic; the relative rate of loss of selection was more rapid in the WCH16ΔcodY (1:6.5), compared to 1:5, and 1:1.4, for the WCH43ΔcodY and D39ΔcodY putative mutants, respectively.
Up-regulated genes are predicted to be co-regulated by similar sets of TFs in various niches
Our analysis also included the identification of genes that are commonly up-regulated by each of the 7 highly activated TFs in the lungs, blood or brain for the 3 strains. We restricted our analysis to include the genes that harbor a minimum of 4 TFBs in their promoter regions in at least 2 of the 3 strains. This showed that the same set of genes (except SP_0698) are commonly up-regulated in the lungs and brain for WCH16 and WCH43 (see Additional file 10: Table S9). However, these genes were not detected during lung infection by D39. By contrast, only one gene (SP_1329), a putative N-acetylneuraminate lyase, was up-regulated in WCH43 and D39 in the blood. These results are complementary to data presented in Additional file 8: Table S7.
Genes harboring higher numbers of regulatory elements constitute a novel index of quality-based gene selection
We developed a new quality-based selection index based on mining of the number of regulatory elements on the promoter regions of up-regulated genes controlled by the 7 highly activated TFs. Specifically, we hypothesized that genes containing higher numbers of TFBs are likely to play key roles in pneumococcal pathogenesis because of their ability to host more TFs. Of the up-regulated genes during infection of lungs and brain by WCH16 and WCH43, SP_0698 (ABC-2 type transport system permease protein) was predicted to contain the highest number of TFBs (13) in its promoter region (see Additional file 11: Table S10). In the blood, WCH43 up-regulates SP_2182 (a bacteriocin_IIc-type protein) with 18 TFBs in its promoter region. However, in D39, SP_0124 and SP_0125 (bacteriocin_IIc-type proteins), SP_0874 and SP_1036 (hypothetical proteins), and SP_1608 (a DadA-family oxidoreductase) (with 23, 26, 25, and 13, TFBs, respectively) were up-regulated.
GO analysis of up-regulated genes containing higher number of TFBs in lungs vs nasopharynx suggests that these genes are involved in key metabolic processes such as NADPH activity, fatty acid biosynthetic process, sugar phosphotransferase system, membrane development, NAD binding, and proteolysis. These processes are important during colonization of the nasopharynx and translocation to lungs by S. pneumoniae. In the blood vs lungs comparisons, genes containing high number of TFBs enrich GO processes of ion binding (magnesium ion binding/iron ion binding), transporter activity, metabolic process, as well as ion transmembrane transport. These processes are vital for survival of S. pneumoniae in the blood, highlighting the important roles of such up-regulated genes with high TFBs in their promoter regions in the blood. GO analysis of up-regulated genes with high TFBs in the brain vs blood comparison were found to be involved in important processes such as NADPH activity, sugar:hydrogen symporter activity, glutamine biosynthetic process, fatty acid biosynthetic process, sugar phosphotransferase activity, kinase activity, carbohydrate transmembrane transport, ion transmembrane transport, and NAD binding. These processes facilitate the translocation, survival and adaptation of WCH16 and WCH43 in the brain.
Domain activation profiles of different pneumococcal strains during pathogenesis
Number of activated protein domains by virulence machinery of different S. pneumoniae strains during movement from the nasopharynx → lungs → blood → brain
Lungs vs nasopharynx
Blood vs lungs
Brain vs blood
Comparison of the activated domain profiles of WCH16 and WCH43 in the brain revealed that the number of activated domains for both strains are essentially similar (105 and 110, respectively; Table 2). The only slight difference is that in WCH43 more activated domains including energy-producing (such as ATP-grasp), two-component regulatory system, tRNA, ribosomal and PTS systems domains were found.
In this work, we have used a functional genomics-based approach to provide a comprehensive prediction of TF activation profiles and domain expression patterns of 3 different S. pneumoniae strains during pathogenesis. It has been suggested that the TF and promoter activation map of an organism can be obtained by comparison with the fully annotated map of an ortholog [22–24] given the fact that functional specificity of proteins is conserved among orthologs . Therefore, we exploited the available E. coli transcriptional regulation data to gain information on the orthologous TF and promoter activation map of 3 S. pneumoniae strains (WCH16, WCH43 and D39) during pathogenesis. However, it has been argued that such methods have shortcomings when applied to distantly related organisms . We attempted to address this potential drawback by using an approach that predicts the function of genes based on the similarity/difference in the pattern of TFBs. The identified potential TFBs and their organization modules provide another way to understand gene expression and regulation during pathogenesis. Taking into account these approaches, we present a strategy for selection of novel antibacterial targets and vaccine candidates.
In a recent investigation , we introduced the concept of bacterial niche-specific virulence gene expression during pathogenesis by qualitative comparisons of transcriptomic data of one niche versus the previous niche. In the present study, we carried out an in-depth analysis of the TF-regulatory network that govern this transition in order to shed light on the virulence factors that S. pneumoniae employs to breach host tissue barriers. Our analysis showed that 7 TFs (argR, codY, hup, rpoD, rr02, scrR and smrC) were highly activated in all 3 strains in vivo and were capable of binding to a large number of up-regulated genes. We suggest that these TFs might constitute the regulatory backbone that underpins pneumococcal pathogenesis. We also found that these TFs were used by D39 in blood and by WCH16 and WCH43 during brain infection.
A close examination of activated TFs of the 3 strains suggests that S. pneumoniae adapts to each host niche during infection by increasing the number of TFs deployed (and their TFBs). As a corollary, the number of activated TFs is reduced if the bacteria are unable to adapt to a particular niche (e.g. WCH16 in the blood). Not surprisingly, rpoD was the most active TF for all 3 strains in all niches, consistent with its role in promoting the attachment of RNA polymerase to specific initiation sites during transcription. The observed similarity between the predicted TF activation profiles and their clustering patterns between distinct niches suggests that different S. pneumoniae strains might activate the same set of TFs and regulatory elements at different sites of infection. This is further supported by discriminant function analysis of the TFs in various niches, which highlights differential regulatory backgrounds of the 3 strains in the blood. We also provide experimental evidence that validates the predicted discriminant function for CodY between the 3 strains by showing distinct pathogenic profiles of their otherwise isogenic codY induced derivatives in various anatomical niches of mice. The results are consistent with our in vivo microarray and real-time RT-PCR data showing differential expression of CodY-regulated virulence genes such as aliA, ilvH and piuA in various niches of mice , and are also in agreement with results from other workers . Furthermore, the predicted domain activation catalogues of these strains during pathogenesis suggests that their functional domain profiles differed mainly in the blood, with only one domain (DNA methyltransferase) activated for WCH16, 14 for WCH43 and 470 domains for D39. These results are consistent with gene ontology classifications for WCH16 and WCH43 [21, 39], and in agreement with our published observations showing that these 3 strains display distinct pathogenicity characteristics [21, 31, 39].
Our analyses also showed that some genes could potentially be regulated by multiple TFs. For example, SP_1647 (a metallo-endopeptidase) was predicted to be regulated by argR, hup, rpoD and rr02 in the lungs of D39-infected mice. Notably, this gene has recently been characterized to have a role in pneumococcal attachment and internalization in host epithelial (A549) and endothelial (HUVEC) cell lines . We also developed a new quality-based strategy for selection of genes important for pathogenesis, based on mining of the number of TFBs in their promoter regions. GO analysis of such genes allowed us to propose that number of TFBs in the promoter regions of up-regulated genes is an important and novel index for selection of genes that play key roles at distinct stages of pneumococcal pathogenesis. In this context, we found, amongst other genes, a group of class IIc bacteriocins that harbor the highest number of TFBs in WCH43 and D39 harvested from the blood. Bacteriocins have been shown to be involved in fratricide in S. pneumoniae, a phenomenon that favors natural competence and pathogenesis in this organism [41, 42]. Additionally, a putative N-acetylneuraminate lyase (SP_1329) was up-regulated by both WCH43 and D39 in blood, suggesting a need for this enzyme for metabolic adaptation of these strains in this niche.
It has been suggested that regulation of pneumococcal virulence proteins is very complex and multifactorial, and likely involves overlapping regulatory mechanisms . We found that all the highly activated TFs argR, codY, hup, rr02, scrR and smrC), except rpoD, clustered together in all 3 strains. This is a significant finding, implying redundancy in the regulatory roles of these TFs during infection, as the functions of these TFs in pathogenesis and regulation of virulence gene expression have hitherto been studied in isolation [7, 33, 34, 44–47]. Indeed, this notion is exemplified by our recent work, which showed that a defined mutation of smrC in a WCH43 background did not completely abrogate virulence . In the same context, we showed in Figure 4 that the same groups of up-regulated genes are predicted to be co-regulated by similar sets of TFs. These findings might explain a fundamental problem associated with prevention and control of pneumococcal disease, as it appears that the task of identifying the principal targets of intervention would be challenging. Therefore, a thorough analysis of the important regulatory pathways employed by S. pneumoniae in different in vivo niches would likely improve our knowledge of pneumococcal pathogenesis and identify novel targets for intervention.
Our findings suggest a possible evolutionary strategy for the pneumococcus in which the evolution occurs in non-coding promoter regions of genes during infection. Specifically, the pneumococcus changes the short binding sites in the promoter regions of genes instead of alteration of coding sequences in a niche-specific manner to generate more virulent strains. There are two important caveats of this study. Firstly, we utilized the fully annotated E. coli promoter listings to pull out the orthologous TFs in S. pneumoniae, which might not represent the full array of active promoters in the pneumococcus. Secondly, the microarray slides used in this study cannot detect the expression of unique WCH16-specific genes that are absent in TIGR4 and R6 genomes represented on the slides. Nevertheless, this study is the first such investigation; it provides a framework towards a better understanding of the dynamics of pneumococcal pathogenesis, and opens new avenues into similar investigations in other pathogenic bacteria.
Bacterial strains and growth conditions
The pneumococcal strains used in this study were D39 (serotype 2; Sequence Type 595), clinical blood isolates WCH43 (serotype 4; Sequence Type 205) and WCH16 (serotype 6A; Sequence Type 4966). Previous mouse intranasal challenge experiments in our laboratory with these strains indicated that D39 causes fulminant bacteraemia; WCH43 demonstrates the “classical” disease progression from the nasopharynx to the lungs, followed by dissemination to the blood and then to the brain, while WCH16 demonstrates minimal lung and blood involvement before translocation to the brain  Serotype-specific capsule production was confirmed by Quellung reaction, as described previously . Opaque-phase variants of the strains, selected on Todd-Hewitt broth supplemented with 1% yeast extract (THY)-catalase plates , were used in all animal experiments. Before infection, the bacteria were grown statically at 37°C in serum broth (SB) to A600 of 0.16 (equivalent to approx. 5 × 107 CFU/ml).
Microarray analysis of RNA extracted from different mouse tissues during pathogenesis
For this investigation, we analyzed 72 h-microarray data of in vivo-derived RNA samples obtained from our previous studies with WCH16 and WCH43 [12, 13] as well as that generated for D39 [12, 13, 31, 50]. Our previous studies have also established that after intranasal infection of CD-1 mice, the number of pneumococci recovered from various tissues peaked at 72 hr after which the mice succumbed to infection quite rapidly . Therefore, we focused on this time point for further differential gene expression analyses. This was carried out by a modified two-color microarray analysis where the relative expression of each gene in one niche was calculated in comparison to expression in the previous niche when bacteria moves from the nasopharynx to lungs to blood and then to brain. Differential gene expression values with p < 0.05 were considered as statistically significant (one-sample t-test). In the case of multiple comparisons, T- Hotelling using Minitab 16 package (http://www.minitab.com) was used to adjust p-values and decrease type I error. The goal of this analysis was to unravel pneumococcal gene expression patterns in the nasopharynx, lungs, blood and brain of mice.
Promoter extraction of up-regulated genes
Comparative analysis of TF activation in different niches
Based on promoter analysis of up-regulated genes, a map of TF activation was estimated in different niches for the 3 strains. To obtain a clear understanding of similarities and/or differences in TF activation, TFs were clustered according to the number of TFBs on promoter regions of up-regulated genes during infection by hierarchical clustering analysis (HCA). This was achieved by count data (number of possible binding sites of each TF to the promoter regions of up-regulated genes in each niche comparison) for all strains. When a particular TF does not have a potential binding site on promoter regions of up-regulated genes, a value of “0” was placed in the data. For each strain, similarities/differences in different TF activation patterns (in terms of potential binding to promoter regions of up-regulated genes) were calculated by average linkage method based on Euclidean distance measure . Dendograms were designed based on calculated similarity matrixes by Minitab 16 package (http://www.minitab.com) to visualise differential binding patterns. In addition, clustering of overall TF activation map in all strains and niches was carried out as described above. The hierarchical clustering analyses were validated using c1Valid R Package , by internal and stability cluster validation measures. Internal measures represent the compactness, connectedness, and separation of the cluster partitions by Connectivity, Silhouette Width and Dunn Index, where lower Connectivity but higher Dunn and Silhouette represent better scores. Stability measures assess the results from the full clustering data based on removing each column, one at a time. The Stability measures evaluated were: Average proportion of non-overlap (APN), Average distance (AD), Average distance between means (ADM) and Figure of merit (FOM); lower APN, AD, ADM, and FOM represent better scores. We then specified best threshold line to generate 2 clusters.
Discriminant function analysis
We used linear discriminant statistical analysis to determine the functions of TFs and used this information to discriminate between the impact of a particular TF in each niche between two groups (e.g. WCH16 vs WCH43, WCH16 vs D39, WCH43 vs D39). This predictive algorithm provides an accurate comparative index of TF function, whereby TFs with relative higher coefficient values (either positive or negative) have significant influence/impact in a particular niche between the two groups being compared.
GO analysis of up-regulated genes with high number of TFBs in their promoter regions
Up-regulated genes containing higher number of TFBs during pneumococcal pathogenesis were subjected to GO classification for their functional association categories with respect to molecular function and biological process, using our recently developed bacterial GO web application [39, 57].
TF–based network construction of pneumococcal pathogenesis
A database was made based on the results of promoter analysis linking the genes to TFs. TF-based regulatory networks were designed in different niches for each strain based on the developed database. The database relationships were visualised using PathwayStudio 9 (Elsevier, USA). The database is available upon request.
Extraction of domain profile of up-regulated genes
Protein sequences of up-regulated genes of WCH16, WCH43, and D39 were extracted from PIR (Protein Information Resource) databases (http://pir.georgetown.edu/cgi-bin/batch.pl). Then, domains of each protein sequence were extracted from pfam database (http://pfam.sanger.ac.uk/) and CLC Genomics Workbench 6 (http://www.clcbio.com/). For domain prediction, E-value = 0.0000001was used to extract the significant domains.
Construction of mutants and assessment of bacterial growth
S. pneumoniae derivatives with marked mutations in codY were constructed in the 3 strains. Mutants were constructed by overlap extension PCR [58, 59] and validated by PCR and sequencing to be in-frame deletion mutation replacements. All PCR procedures were performed with the Phusion High Fidelity Kit (FINNZYMES). The primer pairs used for construction and validation of the mutants are listed in Additional file 13: Table S12. A recombinant plasmid over-expressing codY was constructed by driving its expression from ami promoter using pAL2 plasmid  in which the luxABCDE cassette had been deleted (pAL2ΔluxABCDE). The recombinant construct was transformed into the 3 strains, generating D39[codY]ind, WCH16[codY]ind, and D39[codY]ind, respectively. Induction of codY in these strains was confirmed by real-time RT-PCR analysis.
In vitro competition experiments were performed as described previously . In this assay, mutant (or codYind) and wild-type bacteria were grown to A600 in SB and then mixed at an input ratio of 1:1 in SB. At 1.5 and 3 h post incubation, an aliquot of each sample was serially diluted in SB and plated on blood agar and blood agar with a selective antibiotic to determine the ratio of mutant (or codYind) to wild-type bacteria. Each competition experiment was repeated at least twice. Competitive indices were calculated as the ratio (±SEM) of mutant (or codYind) to wild-type bacteria recovered at each time point adjusted by the input ratio.
Mixed infection experiments
Competition experiments were carried out essentially as described previously . Briefly, 8 mice were anesthetized by intraperitoneal injection of pentobarbital sodium (Nembutal; Rhone-Merieux) at a dose of 66 mg per g of body weight and separately challenged i.n. with 50 μl suspension containing approx. 2.5 × 106 CFU of either wild-type or the isogenic codYind strain. At 36 h post-challenge, mice from each separate infection experiment were sacrificed, bacteria were enumerated from the nasopharynx, lungs blood and brain, by plating on blood agar and blood agar with a selective antibiotic to determine the ratio of the codYind strain to wild-type bacteria. Each competition experiment was repeated at least twice. Competitive indices were calculated as the ratio (±SEM) of codYind strain to wild-type bacteria recovered in each niche adjusted by the input ratio.
Outbred 5- to 6-week-old female CD1 (Swiss) mice were used in all experiments. The Animal Ethics Committee of The University of Adelaide approved all animal experiments (Project Number: S-2013-053). The study was conducted in compliance with the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes (7th Edition 2004) and the South Australian Animal Welfare Act 1985.
The data reported in this paper are archived at the following databases: BμG@Sbase (http://bugs.sgul.ac.uk/E-BUGS-130 and http://bugs.sgul.ac.uk/E-BUGS-133, and also ArrayExpress (accession number E-BUGS-130, and E-BUGS-133).
James C Paton, Abiodun D Ogunniyi and Esmaeil Ebrahimie Joint Senior authors.
This work was supported by the Meningitis Research Foundation (UK) Research Grant 0802.0 to ADO, JCP and LKM, the National Health and Medical Research Council of Australia (NHMRC) Project Grant 627142 to JCP and ADO, and NHMRC Program Grant 565526 to JCP. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We thank Kim LeMessurier for help with D39 microarray, and Manije Mohammadi-dehcheshmeh for help with preparation of figures. We also thank Adam Potter for construction of the pAL2ΔluxABCDE plasmid. We acknowledge BμG@S (the Bacterial Microarray Group at St George’s, University of London) for supply of the microarray slides and advice and The Wellcome Trust for funding the multi-collaborative microbial pathogen microarray facility under its Functional Genomics Resources Initiative.
- Paton JC, Boslego JW: Protein vaccines. Pneumococcal vaccines: The impact of conjugate vaccines. Edited by: Siber GR, Klugman KP, Makela PH. 2008, Washington DC: ASM Press, 421-436.View ArticleGoogle Scholar
- O'Brien KL, Wolfson LJ, Watt JP, Henkle E, Deloria-Knoll M, McCall N, Lee E, Mulholland K, Levine OS, Cherian T: Burden of disease caused by Streptococcus pneumoniae in children younger than 5 years: global estimates. Lancet. 2009, 374 (9693): 893-902. 10.1016/S0140-6736(09)61204-6.PubMedView ArticleGoogle Scholar
- Reinert RR: The antimicrobial resistance profile of Streptococcus pneumoniae. Clin Microbiol Infect. 2009, 15 (Suppl 3): 7-11.PubMedView ArticleGoogle Scholar
- Bogaert D, De Groot R, Hermans PW: Streptococcus pneumoniae colonisation: the key to pneumococcal disease. Lancet Infect Dis. 2004, 4 (3): 144-154. 10.1016/S1473-3099(04)00938-7.PubMedView ArticleGoogle Scholar
- Kadioglu A, Weiser JN, Paton JC, Andrew PW: The role of Streptococcus pneumoniae virulence factors in host respiratory colonization and disease. Nat Rev Microbiol. 2008, 6 (4): 288-301. 10.1038/nrmicro1871.PubMedView ArticleGoogle Scholar
- van der Poll T, Opal SM: Pathogenesis, treatment, and prevention of pneumococcal pneumonia. Lancet. 2009, 374 (9700): 1543-1556. 10.1016/S0140-6736(09)61114-4.PubMedView ArticleGoogle Scholar
- Hava DL, Camilli A: Large-scale identification of serotype 4 Streptococcus pneumoniae virulence factors. Mol Microbiol. 2002, 45 (5): 1389-1406.PubMed CentralPubMedGoogle Scholar
- Marra A, Asundi J, Bartilson M, Lawson S, Fang F, Christine J, Wiesner C, Brigham D, Schneider WP, Hromockyj AE: Differential fluorescence induction analysis of Streptococcus pneumoniae identifies genes involved in pathogenesis. Infect Immun. 2002, 70 (3): 1422-1433. 10.1128/IAI.70.3.1422-1433.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Orihuela CJ, Radin JN, Sublett JE, Gao G, Kaushal D, Tuomanen EI: Microarray analysis of pneumococcal gene expression during invasive disease. Infect Immun. 2004, 72 (10): 5582-5596. 10.1128/IAI.72.10.5582-5596.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- van Opijnen T, Camilli A: A fine scale phenotype-genotype virulence map of a bacterial pathogen. Genome Res. 2012, 22 (12): 2541-2551. 10.1101/gr.137430.112.PubMed CentralPubMedView ArticleGoogle Scholar
- Molzen TE, Burghout P, Bootsma HJ, Brandt CT, van der Gaast-de Jongh CE, Eleveld MJ, Verbeek MM, Frimodt-Moller N, Ostergaard C, Hermans PW: Genome-wide identification of Streptococcus pneumoniae genes essential for bacterial replication during experimental meningitis. Infect Immun. 2011, 79 (1): 288-297. 10.1128/IAI.00631-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Mahdi LK, Wang H, Van der Hoek MB, Paton JC, Ogunniyi AD: Identification of a novel pneumococcal vaccine antigen preferentially expressed during meningitis in mice. J Clin Invest. 2012, 122 (6): 2208-2220. 10.1172/JCI45850.PubMed CentralPubMedView ArticleGoogle Scholar
- Ogunniyi AD, Mahdi LK, Trappetti C, Verhoeven N, Mermans D, Van der Hoek MB, Plumptre CD, Paton JC: Identification of genes that contribute to the pathogenesis of invasive pneumococcal disease by in vivo transcriptomic analysis. Infect Immun. 2012, 80 (9): 3268-3278. 10.1128/IAI.00295-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Deihimi T, Niazi A, Ebrahimi M, Kajbaf K, Fanaee S, Bakhtiarizadeh MR, Ebrahimie E: Finding the undiscovered roles of genes: an approach using mutual ranking of coexpressed genes and promoter architecture-case study: dual roles of thaumatin like proteins in biotic and abiotic stresses. SpringerPlus. 2012, 1: 30-10.1186/2193-1801-1-30.PubMed CentralPubMedView ArticleGoogle Scholar
- Schoolnik GK: Microarray Analysis of Bacterial Pathogenicity. Advances in Microbial Physiology. Edited by: Poole RK. 2002, Elsevier USA: Elsevier, 46: 1-45.Google Scholar
- Balleza E, Lopez-Bojorquez LN, Martinez-Antonio A, Resendis-Antonio O, Lozada-Chavez I, Balderas-Martinez YI, Encarnacion S, Collado-Vides J: Regulation by transcription factors in bacteria: beyond description. FEMS Microbiol Rev. 2009, 33 (1): 133-151. 10.1111/j.1574-6976.2008.00145.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Cases I, de Lorenzo V, Ouzounis CA: Transcription regulation and environmental adaptation in bacteria. Trends Microbiol. 2003, 11 (6): 248-253. 10.1016/S0966-842X(03)00103-3.PubMedView ArticleGoogle Scholar
- McAllister LJ, Tseng HJ, Ogunniyi AD, Jennings MP, McEwan AG, Paton JC: Molecular analysis of the psa permease complex of Streptococcus pneumoniae. Mol Microbiol. 2004, 53 (3): 889-901. 10.1111/j.1365-2958.2004.04164.x.PubMedView ArticleGoogle Scholar
- Ogunniyi AD, Mahdi LK, Jennings MP, McEwan AG, McDevitt CA, Van der Hoek MB, Bagley CJ, Hoffmann P, Gould KA, Paton JC: Central role of manganese in regulation of stress responses, physiology, and metabolism in Streptococcus pneumoniae. J Bacteriol. 2010, 192 (17): 4489-4497. 10.1128/JB.00064-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Marra A, Lawson S, Asundi JS, Brigham D, Hromockyj AE: In vivo characterization of the psa genes from Streptococcus pneumoniae in multiple models of infection. Microbiology. 2002, 148 (Pt 5): 1483-1491.PubMedView ArticleGoogle Scholar
- Mahdi LK, Ebrahimie E, Adelson DL, Paton JC, Ogunniyi AD: A transcription factor contributes to pathogenesis and virulence in Streptococcus pneumoniae. PLoS One. 2013, 8 (8): e70862-10.1371/journal.pone.0070862.PubMed CentralPubMedView ArticleGoogle Scholar
- Covert MW, Knight EM, Reed JL, Herrgard MJ, Palsson BO: Integrating high-throughput and computational data elucidates bacterial networks. Nature. 2004, 429 (6987): 92-96. 10.1038/nature02456.PubMedView ArticleGoogle Scholar
- van Schaik W, van der Voort M, Molenaar D, Moezelaar R, de Vos WM, Abee T: Identification of the σBregulon of Bacillus cereusand conservation of σB-regulated genes in low-GC-content gram-positive bacteria. J Bacteriol. 2007, 189 (12): 4384-4390. 10.1128/JB.00313-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Lintner RE, Mishra PK, Srivastava P, Martinez-Vaz BM, Khodursky AB, Blumenthal RM: Limited functional conservation of a global regulator among related bacterial genera: Lrp in Escherichia, Proteus and Vibrio. BMC Microbiol. 2008, 8: 60-10.1186/1471-2180-8-60.PubMed CentralPubMedView ArticleGoogle Scholar
- Mirny LA, Gelfand MS: Using orthologous and paralogous proteins to identify specificity-determining residues in bacterial transcription factors. J Mol Biol. 2002, 321 (1): 7-20. 10.1016/S0022-2836(02)00587-9.PubMedView ArticleGoogle Scholar
- Hui FM, Morrison DA: Genetic transformation in Streptococcus pneumoniae: nucleotide sequence analysis shows comA, a gene required for competence induction, to be a member of the bacterial ATP-dependent transport protein family. J Bacteriol. 1991, 173 (1): 372-381.PubMed CentralPubMedGoogle Scholar
- Lange R, Wagner C, de Saizieu A, Flint N, Molnos J, Stieger M, Caspers P, Kamber M, Keck W, Amrein KE: Domain organization and molecular characterization of 13 two-component systems identified by genome sequencing of Streptococcus pneumoniae. Gene. 1999, 237 (1): 223-234. 10.1016/S0378-1119(99)00266-8.PubMedView ArticleGoogle Scholar
- Berry AM, Paton JC: Sequence heterogeneity of PsaA, a 37-kilodalton putative adhesin essential for virulence of Streptococcus pneumoniae. Infect Immun. 1996, 64 (12): 5255-5262.PubMed CentralPubMedGoogle Scholar
- Pan XS, Ambler J, Mehtar S, Fisher LM: Involvement of topoisomerase IV and DNA gyrase as ciprofloxacin targets in Streptococcus pneumoniae. Antimicrob Agents Chemother. 1996, 40 (10): 2321-2326.PubMed CentralPubMedGoogle Scholar
- Orihuela CJ, Gao G, McGee M, Yu J, Francis KP, Tuomanen E: Organ-specific models of Streptococcus pneumoniae disease. Scand J Infect Dis. 2003, 35 (9): 647-652. 10.1080/00365540310015854.PubMedView ArticleGoogle Scholar
- Mahdi LK, Ogunniyi AD, LeMessurier KS, Paton JC: Pneumococcal virulence gene expression and host cytokine profiles during pathogenesis of invasive disease. Infect Immun. 2008, 76 (2): 646-657. 10.1128/IAI.01161-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Brock G, Pihur V, Datta S, Datta S: clValid, an R package for cluster validation. J Stat Soft. 2008, 25 (4): 1-22.View ArticleGoogle Scholar
- Caymaris S, Bootsma HJ, Martin B, Hermans PW, Prudhomme M, Claverys JP: The global nutritional regulator CodY is an essential protein in the human pathogen Streptococcus pneumoniae. Mol Microbiol. 2010, 78 (2): 344-360. 10.1111/j.1365-2958.2010.07339.x.PubMedView ArticleGoogle Scholar
- Hendriksen WT, Bootsma HJ, Estevao S, Hoogenboezem T, de Jong A, de Groot R, Kuipers OP, Hermans PW: CodY of Streptococcus pneumoniae: link between nutritional gene regulation and colonization. J Bacteriol. 2008, 190 (2): 590-601. 10.1128/JB.00917-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Martin RG, Rosner JL: Genomics of the marA/soxS/rob regulon of Escherichia coli: identification of directly activated promoters by application of molecular genetics and informatics to microarray data. Mol Microbiol. 2002, 44 (6): 1611-1624. 10.1046/j.1365-2958.2002.02985.x.PubMedView ArticleGoogle Scholar
- Chaturongakul S, Raengpradub S, Palmer ME, Bergholz TM, Orsi RH, Hu Y, Ollinger J, Wiedmann M, Boor KJ: Transcriptomic and phenotypic analyses identify coregulated, overlapping regulons among PrfA, CtsR, HrcA, and the alternative sigma factors σB, σC, σH, and σLin Listeria monocytogenes. Appl Environ Microbiol. 2011, 77 (1): 187-200. 10.1128/AEM.00952-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang X, Fredrick KL, Helmann JD: Promoter recognition by Bacillus subtilis sigmaW: autoregulation and partial overlap with the sigmaX regulon. J Bacteriol. 1998, 180 (15): 3765-3770.PubMed CentralPubMedGoogle Scholar
- Price MN, Dehal PS, Arkin AP: Orthologous transcription factors in bacteria have different functions and regulate different genes. Plos Comput Biol. 2007, 3 (9): 1739-1750.PubMedView ArticleGoogle Scholar
- Fruzangohar M, Ebrahimie E, Ogunniyi AD, Mahdi LK, Paton JC, Adelson DL: Comparative GO: a web application for comparative gene ontology and gene ontology-based gene selection in bacteria. PLoS One. 2013, 8 (3): e58759-10.1371/journal.pone.0058759.PubMed CentralPubMedView ArticleGoogle Scholar
- Agarwal V, Kuchipudi A, Fulde M, Riesbeck K, Bergmann S, Blom AM: Streptococcus pneumoniae Endopeptidase O (PepO) is a multifunctional plasminogen- and fibronectin-binding protein, facilitating evasion of innate immunity and invasion of host cells. J Biol Chem. 2013, 288 (10): 6849-6863. 10.1074/jbc.M112.405530.PubMed CentralPubMedView ArticleGoogle Scholar
- Guiral S, Mitchell TJ, Martin B, Claverys JP: Competence-programmed predation of noncompetent cells in the human pathogen Streptococcus pneumoniae: genetic requirements. Proc Natl Acad Sci U S A. 2005, 102 (24): 8710-8715. 10.1073/pnas.0500879102.PubMed CentralPubMedView ArticleGoogle Scholar
- Weng L, Piotrowski A, Morrison DA: Exit from competence for genetic transformation in Streptococcus pneumoniae is regulated at multiple levels. PLoS One. 2013, 8 (5): e64197-10.1371/journal.pone.0064197.PubMed CentralPubMedView ArticleGoogle Scholar
- Ogunniyi AD, Paton JC: Regulation of pneumococcal surface proteins and capsule. Regulation of Bacterial Virulence. Edited by: Vasil ML, Darwin AJ. 2012, Washington DC: ASM Press, 490-208.Google Scholar
- Kadioglu A, Echenique J, Manco S, Trombe MC, Andrew PW: The MicAB two-component signaling system is involved in virulence of Streptococcus pneumoniae. Infect Immun. 2003, 71 (11): 6676-6679. 10.1128/IAI.71.11.6676-6679.2003.PubMed CentralPubMedView ArticleGoogle Scholar
- Ng WL, Tsui HC, Winkler ME: Regulation of the pspA virulence factor and essential pcsB murein biosynthetic genes by the phosphorylated VicR (YycF) response regulator in Streptococcus pneumoniae. J Bacteriol. 2005, 187 (21): 7444-7459. 10.1128/JB.187.21.7444-7459.2005.PubMed CentralPubMedView ArticleGoogle Scholar
- Iyer R, Camilli A: Sucrose metabolism contributes to in vivo fitness of Streptococcus pneumoniae. Mol Microbiol. 2007, 66 (1): 1-13. 10.1111/j.1365-2958.2007.05878.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Kloosterman TG, Kuipers OP: Regulation of arginine acquisition and virulence gene expression in the human pathogen Streptococcus pneumoniae by transcription regulators ArgR1 and AhrC. J Biol Chem. 2011, 286 (52): 44594-44605. 10.1074/jbc.M111.295832.PubMed CentralPubMedView ArticleGoogle Scholar
- Berry AM, Paton JC: Additive attenuation of virulence of Streptococcus pneumoniae by mutation of the genes encoding pneumolysin and other putative pneumococcal virulence proteins. Infect Immun. 2000, 68 (1): 133-140. 10.1128/IAI.68.1.133-140.2000.PubMed CentralPubMedView ArticleGoogle Scholar
- Weiser JN, Austrian R, Sreenivasan PK, Masure HR: Phase variation in pneumococcal opacity: relationship between colonial morphology and nasopharyngeal colonization. Infect Immun. 1994, 62 (6): 2582-2589.PubMed CentralPubMedGoogle Scholar
- LeMessurier KS, Ogunniyi AD, Paton JC: Differential expression of key pneumococcal virulence genes in vivo. Microbiology. 2006, 152 (Pt 2): 305-311.PubMedView ArticleGoogle Scholar
- Dehal PS, Joachimiak MP, Price MN, Bates JT, Baumohl JK, Chivian D, Friedland GD, Huang KH, Keller K, Novichkov PS, Dubchak IL, Alm EJ, Arkin AP: MicrobesOnline: an integrated portal for comparative and functional genomics. Nucleic Acids Res. 2010, 38 (Suppl 1): D396-D400.PubMed CentralPubMedView ArticleGoogle Scholar
- 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 (7): 1015-1022. 10.1101/gr.3844805.PubMed CentralPubMedView ArticleGoogle Scholar
- Solovyev V, Salamov A: Automatic annotation of microbial genomes and metagenomic sequences. Metagenomics and its applications in agriculture, biomedicine and environmental studies. Edited by: Li RW. 2011, New York: Nova Science Publishers, 61-78.Google Scholar
- Lee W, Chen SL: Genome-Tools: A Flexible Package for Genome Sequence Analysis. Biotechniques. 2002, 33 (6): 1334-1341.PubMedGoogle Scholar
- Novichkov PS, Kazakov AE, Ravcheev DA, Leyn SA, Kovaleva GY, Sutormin RA, Kazanov MD, Riehl W, Arkin AP, Dubchak I, Rodionov DA: RegPrecise 3.0 - A resource for genome-scale exploration of transcriptional regulation in bacteria. BMC Genomics. 2013, 14: 745-10.1186/1471-2164-14-745.PubMed CentralPubMedView ArticleGoogle Scholar
- Manly BFJ: Multivariate statistical methods: a primer. 2005, Florida: CRC Press, 214-3Google Scholar
- Fruzangohar M, Ebrahimie E, Adelson DL: Application of global transcriptome data in gene ontology classification and construction of a gene ontology interaction network. bioRxiv. 2014, doi: http://dx.doi.org/10.1101/004911Google Scholar
- Horton RM, Ho SN, Pullen JK, Hunt HD, Cai Z, Pease LR: Gene splicing by overlap extension. Methods Enzymol. 1993, 217: 270-279.PubMedView ArticleGoogle Scholar
- Ogunniyi AD, Grabowicz M, Mahdi LK, Cook J, Gordon DL, Sadlon TA, Paton JC: Pneumococcal histidine triad proteins are regulated by the Zn2+-dependent repressor AdcR and inhibit complement deposition through the recruitment of complement factor H. FASEB J. 2009, 23 (3): 731-738. 10.1096/fj.08-119537.PubMedView ArticleGoogle Scholar
- Beard SJ, Salisbury V, Lewis RJ, Sharpe JA, MacGowan AP: Expression of lux genes in a clinical isolate of Streptococcus pneumoniae: using bioluminescence to monitor gemifloxacin activity. Antimicrob Agents Chemother. 2002, 46 (2): 538-542. 10.1128/AAC.46.2.538-542.2002.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.