- Research article
- Open Access
Comparative transcriptomic analysis of streptococcus pseudopneumoniae with viridans group streptococci
BMC Microbiology volume 12, Article number: 77 (2012)
Streptococcus pseudopneumoniae, is a novel member of the genus Streptococcus, falling close to related members like S. pneumoniae, S. mitis, and S. oralis. Its recent appearance has shed light on streptococcal infections, which has been unclear till recently. In this study, the transcriptome of S. pseudopneumoniae CCUG 49455T was analyzed using the S. pneumoniae R6 microarray platform and compared with those of S. pneumoniae KCTC 5080T, S. mitis KCTC 3556T, and S. oralis KCTC 13048T strains.
Comparative transcriptome analysis revealed the extent of genetic relatedness among the species, and implies that S. pseudopneumoniae is the most closely related to S. pneumoniae. A total of 489, 444 and 470 genes were upregulated while 347, 484 and 443 were downregulated relative to S. pneumoniae in S. pseudopneumoniae, S. oralis and S. mitis respectively. Important findings were the up-regulation of TCS (two component systems) and transposase which were found to be specific to S. pseudopneumoniae.
This study provides insight to the current understanding of the genomic content of S. pseudopneumoniae. The comparative transcriptome analysis showed hierarchical clustering of expression data of S. pseudopneumoniae with S. pneumoniae and S. mitis with S. oralis. This proves that transcriptional profiling can facilitate in elucidating the genetic distance between closely related strains.
Streptococcus pseudopneumoniae is a recently described member of the ‘S. mitis’ group of viridians streptococci, which is phenotypically and genetically close to S. pneumoniae S. mitis, and S. oralis. S. pseudopneumoniae strains characterized to date has been isolated from the lower respiratory tract [2–4]. This species is known to cause infections in patients having a history of chronic obstructive pulmonary disease or exacerbation of chronic obstructive pulmonary disease [4, 5]. However, the clinical significance of this species is currently unknown.
Streptococcus pneumoniae is the most common cause of well-defined clinical syndrome of pneumonia, bacterial meningitis, and nongonoccal urethritis in humans [6–8]. By contrast, two medically important ‘S. mitis’ group streptococci, S. mitis and S. oralis are recognized as important etiological agents for subacute endocarditis and septicaemia [9, 10]. Recently, pancreatic cancer has been associated with S. mitis, increasing the clinical relevance of this group .
The pathogenicity and the underlying genetic identity of S. pseudopneumoniae are not well characterized in relation to its phylogenetic neighbours, S. pneumoniae, S. mitis, and S. oralis. Unlike S. pneumoniae S. pseudopneumoniae is optochin resistant in the presence of 5% CO2, is bile insoluble, and lacks the pneumococcal capsule [12, 13]. The use of MLST described in this paper allowed a good differentiation between the species . In clinical studies, the phenotypic characterization of the isolates showed relatedness to the species S. pseudopneumoniae, but genotypically it was difficult to distinguish from its close neighbour S. pneumoniae. Indeed, S. pseudopneumoniae shares over 99% 16S rRNA gene homology with S. pneumoniae, S. mitis, and S. oralis showing that it has evolved from a common genetic ancestor [16–18]. In recent years, several reports have shown that S. pneumoniae share genes encoding virulence factors with S. mitis and S. oralis, providing suggestive evidence of lateral gene transfer between these species [19, 20].
Genotypic characterization of S. pseudopneumoniae in relation to its neighboring members is necessary to increase its clinical relevance. Comparative genomics or transcriptomics based on genome wide microarrays , is now the logical approach used to determine inter-species comparisons [22, 23]. Since whole-genome sequencing to elucidate the genetic content of a microorganism is considered to be expensive and time consuming, an approach used for the identification of large number of genes without the need for sequencing is the trend in present era. The entire genomes of S. pneumoniae S. mitis, and S. oralis have been fully sequenced. However, transcriptome has not been studied in these microorganisms to date, which may lead to the identification of unique virulence genes specific to the strain of interest.
Previously, we identified species-specific genes using suppressive subtractive hybridization (SSH), such as the cpsA gene for S. pneumoniae and the rgg gene for S. oralis[24–26]. In the current study, the gene expression of S. pseudopneumoniae is determined and compared with those of S. pneumoniae KCTC 5080TS. mitis KCTC 3556T and S. oralis KCTC 13048T by in silico analysis and by in vitro transcriptome microarrays experiments using open reading frame (ORF) microarrays of Streptococcus pneumoniae R6 (GenBank accession number NC_003098) platform.
Results and discussion
Statistical analysis of microarray experiments
We compared the expression profiles by hybridization to the immobilized probes on the microarray of S. pneumoniae TIGR4: NC_003028 with the total RNA of S. oralis KCTC 13048T, S. mitis KCTC 3556T, and S. pseudopneumoniae CCUG 49455T. Total RNA from the strains S. pneumoniae KCTC 5080T, S. mitis KCTC 3556T, S. oralis KCTC 13048T, and S. pseudopneumoniae CCUG 49455T was hybridized to NimbleGen S. pneumoniae TIGR4: NC_003028 Gene Expression 4x72K microarrays. Each array contains 4 sets of strains, and each strain was compared with each other strains. Interarray correlation values (Range: -1 ≤ r ≤ 1) are shown in the upper right panels and pairwise scatter plots of gene expression values (log2) are shown in the lower left panels (Figure 1). A correlation value close to 1 shows high similarity between samples. This correlation value between strains S. oralis-S. mitis was 0.609, S. oralis-S. pneumoniae was 0.365, S. oralis-S. pseudopneumoniae was 0.375, S. mitis-S. pneumoniae was 0.438, S. mitis-S. pseudopneumoniae was 0.536 and S. pneumoniae-S. pseudopneumoniae was 0.499.
Phylogenetic relatedness between streptococcal species
Based on their overall genomic profiles, there was clear delineation between each Streptococcus species. The hierarchical clustering analysis from a normalized signal grouped the isolates mainly according to their phylogenetic relationship between each Streptococcus species. The clustering of S. mitis, S. oralis and S. pneumoniae, S. pseudopneumoniae strains showed two distinct branches, placing them in two separate clades that clearly differentiated each species group (Figure 2). The map shows the expression levels of the 1,123 probes (Figure 3). A total of 444 genes were upregulated (red) and 484 genes were downregulated(green) in S. oralis KCTC 13048T, 470 genes were upregulated (red) and 443 genes were downregulated (green) in S. mitis KCTC 3556T and 489 genes were upregulated (red) and 347 genes were downregulated (green) in S. pseudopneumoniae CCUG 49455T (Figure 3). Red represents high expression; green represents low expression (Figure 4).
Identification of functional genes revealed by transcriptome analysis
Whole-genome sequence of S. pseudopneumoniae (isolate number: IS7493, GenBank accession numbers: CP002925 and CP002926) was done by Shahinas et al.. Their study shows the presence or absence of genes in the whole genome but not the functional analysis of RNA transcripts. In this study, the availability of the complete S. pneumoniae TIGR4: NC_003028 genome  allowed for the analysis of S. oralis KCTC 13048T total RNA transscripts.
About 53 genes were up regulated in S. oralis KCTC 13048T when compared with other Streptococcus species (Table 1). About 26 genes were identified as hypothetical proteins while the remaining 27 were associated with amino acid biosynthesis, transport and degenerate transposase proteins.
The 37 genes differentially regulated in S. mitis KCTC 3556T were found to function in amino acid biosynthesis, transport and were transposases, including 4’-phosphopantetheinyl transferase, ABC transporter, alcohol dehydrogenase, alkaline amylopullulanase, Smf DNA processing protein, MSM (multiple sugar metabolism) operon regulatory protein, Peptidoglycan GlcNAc deacetylase, Phosphatidate cytidylyltransferase, RecA regulator RecX, Transport protein ComB, UDP-galactose 4-epimerase, truncation, as well as other hypothetical proteins (Table 1).
The 117 upregulated genes of S. pseudopneumoniae CCUG 49455T, were found to play a role in amino acid biosynthesis and transport, such as ABC transporter ATP-binding protein, conserved hypothetical protein, D-alanine glycine permease, histidine kinase, major facilitator superfamily transporter, maltose operon transcriptional repressor, mannitol PTS EII, mannitol-1-phosphate 5-dehydrogenase, mannitol-specific enzyme IIA component, negative regulator of pho regulon for phosphate transport, peptidoglycan GlcNAc deacetylase, phosphotransferase system, positive transcriptional regulator of mutA, response regulator, riboflavin synthase, sortase and transcriptional proteins.
The degenerate transposon was significantly overexpressed in S. pseudopneumoniae compared to its expression in S. oralis and S. mitis. On the other hand, histidine kinase and response regulators associated with the two component system (TCS) were down regulated in the S. oralis and S. mitis (Table 1). Additionally pneumolysin and penicillin-binding protein were also down regulated in S. oralis and S. mitis and showed no signal in the S. pseudopneumoniae.
Upregulation of some interesting genes in the transport group was found in S. pseudopneumoniae like the ATP-binding cassette (ABC) transporters and the two component system (TCS). ABC transporters are integral membrane proteins that actively transport chemically diverse substrates across the lipid bilayers of cellular membranes. This is of clinical importance because multidrug resistance in human cancer cells is mostly the result of the over expression of ABC transporters that catalyze the extrusion of the cytotoxic compounds used in cancer therapy . Bacterial drug resistance has become an increasing problem. In bacterial cells, ABC transporters are known to contribute to multidrug and antibiotic resistance by extruding drugs or antibiotics .
The TCSs of bacteria consist of two proteins, histidine kinase and response regulators, and have received increasing attention for their potential as a novel antibacterial drug targets [31, 32]. Some TCSs regulate the expression of antibiotic resistance determinants, including drug-efflux pumps . The overexpression of response regulators of bacterial two-component signal transduction system confers drug resistance by controlling the expression of some drug transporter genes. Various TCSs ubiquitously present in bacteria regulate the transcription of different gene products. The regulation of osmolarity, nutrient uptake, redox potential, sporulation and the expression of virulence factors are under the control of TCSs. The two component system (TCS) serves as a basic stimulus–response coupling mechanism that allows organisms to sense and respond to changes in environmental conditions. The sensor kinase monitors a certain environmental condition and modulates the phosphorylation state of the response regulator that controls genes. One of the most attractive aspects of the TCS is its regulation of antimicrobial resistance factors.
In summary, based on comparative genomics/transcriptome analysis, using S. pneumoniae as the control strain, facilitated the identification of S. pseudopneumoniae transcriptome within streptococci viridans group. We postulate that transcriptional profiling with high statistical power implies the great genetic distance between each streptococci of viridans group. The correlation values by statistical analysis show the closest association between S. oralis and S.mitis. This is also clearly shown by the clustering method which placed S.oralis and S.mitis in a separate clade from S.pneumoniae and S. pseudopneumoniae revealing their genetic relatedness. Overall expression levels of 489 genes were higher in S.mitis strain when compared with the control strain. Some of the important genes identified by functional analysis at RNA level were those belonging to amino acid biosynthesis, transport and degenerate transposase proteins. One of the significant findings in this study was the upregulation of ABC transporters and TCS in S. pseudopneumoniae where the former are known to play a role multi-drug antibiotic resistance and the latter in controlling the virulence factors. Therefore, we conclude by this study that genetic relatedness and pathogenecity in S. pseudopneumoniae in comparison to viridans group was well revealed by transcriptome analysis.
Bacterial culture, RNA extraction and cDNA synthesis
S. pneumoniae KCTC 5080T was used as the reference strain for comparative microarray experiments with other viridians group of streptococci. S. pneumoniae KCTC 5080T, S. pseudopneumoniae CCUG 49455T, S. mitis KCTC 3556T, and S. oralis KCTC 13048T strains were grown on Brain Heart Infusion (BHI) agar (Difco, Detroit, MI, U.S.A.) at 37°C for 18 hours. Total RNA was isolated using a RiboPure Bacteria Kit (Ambion, UK) following manufacturer’s instructions. Extracted RNA was treated with TURBO DNase (Ambion). RNA quality was checked for purity and integrity as evaluated by OD 260/280 ratio, and analyzed on Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, USA). cDNA was synthesized according to the NimbleGen Expression protocol (Nimblegen, Madison, USA) using the SuperScript double-stranded cDNA synthesis kit (Invitrogen Life Technologies, Carlsbad, CA, U.S.A.). Briefly, 10 μg of total RNA was reverse-transcribed to cDNA using an oligo dT primer. Then second-strand cDNA was synthesized. After purification, cDNA was quantified using the ND-1000 Spectrophotometer (NanoDrop, Wilmington, USA).
Labeling and purification
cDNA was labelled using the One-Color Labelling Kit (Nimblegen) following manufacturer’s instructions. 1 μg of cDNA samples were labelled with Cy3 using Cy3-random nonamer. After purification, the labelled cDNA was quantified using the ND-1000 Spectrophotometer (NanoDrop).
Generation of microarray data
The Streptococcus pneumoniae R6 microarrays (Nimblegen) were used for the transcriptome analysis. The S. pneumoniae R6 microarray contains 2,037 genes: 4 × 72,000 probes and 5 replicates (GenBank accession numbers: NC_003098). Labelled cDNA samples of S. pseudopneumoniae S. mitis and S. oralis were hybridized onto Nimblegen Expression array (Nimblegen) for 16-20 hours at 42°C, according to manufacturer's instructions. Arrays were scanned with a NimbleGen MS 200 Microarray scanner set- at 532 nm with a resolution of 2 μm to produce images in TIFF format according to the manufacturer's instructions. Array data export processing and analysis was performed using NimbleScan (version 2.5). The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series accession number GSE37539 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE37539).
Data acquisition and statistical analysis
Raw data was extracted using NimbleScan (version 2.5, Gene Expression RMA algorithm). A single raw intensity value was determined for each gene in each array with 2535 genes by taking an average of spot replicates of all 24 probes. Gene signal value was determined by logarithmic transformation (base 2). Statistical significance of the expression data was determined using fold change. Hierarchical cluster analysis was performed using complete linkage and Euclidean distance as a measure of similarity. NimbleScan was used for quantification, image analysis of mRNA data. R scripts (‘R’ software) were used for all other analytical process.
Arbique JC, Poyart C, Trieu-Cuot P, Quesne G, Carvalho Mda G, Steigerwalt AG, Morey RE, Jackson D, Davidson RJ, Facklam RR: Accuracy of phenotypic and genotypic testing for identification of Streptococcus pneumoniae and description of Streptococcus pseudopneumoniae sp. nov. J Clin Microbiol. 2004, 42 (10): 4686-4696. 10.1128/JCM.42.10.4686-4696.2004.
Carvalho Mda G, Tondella ML, McCaustland K, Weidlich L, McGee L, Mayer LW, Steigerwalt A, Whaley M, Facklam RR, Fields B, et al.: Evaluation and improvement of real-time PCR assays targeting lytA, ply, and psaA genes for detection of pneumococcal DNA. J Clin Microbiol. 2007, 45 (8): 2460-2466. 10.1128/JCM.02498-06.
Cochetti I, Vecchi M, Mingoia M, Tili E, Catania MR, Manzin A, Varaldo PE, Montanari MP: Molecular characterization of pneumococci with efflux-mediated erythromycin resistance and identification of a novel mef gene subclass, mef(I). Antimicrob Agents Chemother. 2005, 49 (12): 4999-5006. 10.1128/AAC.49.12.4999-5006.2005.
Keith ER, Podmore RG, Anderson TP, Murdoch DR: Characteristics of Streptococcus pseudopneumoniae isolated from purulent sputum samples. J Clin Microbiol. 2006, 44 (3): 923-927. 10.1128/JCM.44.3.923-927.2006.
Harf-Monteil C, Granello C, Le Brun C, Monteil H, Riegel P: Incidence and pathogenic effect of Streptococcus pseudopneumoniae. J Clin Microbiol. 2006, 44 (6): 2240-2241. 10.1128/JCM.02643-05.
Marrie TJ, Durant H, Yates L: Community-acquired pneumonia requiring hospitalization: 5-year prospective study. Rev Infect Dis. 1989, 11 (4): 586-599. 10.1093/clinids/11.4.586.
Schmidt A, Bisle B, Kislinger T: Quantitative peptide and protein profiling by mass spectrometry. Meth Mol Biol. 2009, 492: 21-38. 10.1007/978-1-59745-493-3_2.
Fine MJ, Smith MA, Carson CA, Mutha SS, Sankey SS, Weissfeld LA, Kapoor WN: Prognosis and outcomes of patients with community-acquired pneumonia. A meta-analysis. JAMA. 1996, 275 (2): 134-141. 10.1001/jama.1996.03530260048030.
Dyson C, Barnes RA, Harrison GA: Infective endocarditis: an epidemiological review of 128 episodes. J Infect. 1999, 38 (2): 87-93. 10.1016/S0163-4453(99)90074-9.
Willcox MD, Drucker DB, Hillier VF: In-vitro adherence of oral streptococci in the presence of sucrose and its relationship to cariogenicity in the rat. Arch Oral Biol. 1988, 33 (2): 109-113. 10.1016/0003-9969(88)90053-2.
Farrell JJ, Zhang L, Zhou H, Chia D, Elashoff D, Akin D, Paster BJ, Joshipura K, Wong DT: Variations of oral microbiota are associated with pancreatic diseases including pancreatic cancer. Gut. 2012, 61 (4): 582-588. 10.1136/gutjnl-2011-300784.
Johnston C, Hinds J, Smith A, van der Linden M, Van Eldere J, Mitchell TJ: Detection of large numbers of pneumococcal virulence genes in streptococci of the mitis group. J Clin Microbiol. 2010, 48 (8): 2762-2769. 10.1128/JCM.01746-09.
Simoes AS, Sa-Leao R, Eleveld MJ, Tavares DA, Carrico JA, Bootsma HJ, Hermans PW: Highly penicillin-resistant multidrug-resistant pneumococcus-like strains colonizing children in Oeiras, Portugal: genomic characteristics and implications for surveillance. J Clin Microbiol. 2010, 48 (1): 238-246. 10.1128/JCM.01313-09.
Do T, Jolley KA, Maiden MCJ, Gilbert SC, Clark D, Wade WG DB: Population structure of Streptococcus oralis. Microbiology. 2009, 155: 2593-2602. 10.1099/mic.0.027284-0.
Suzuki N, Seki M, Nakano Y, Kiyoura Y, Maeno M, Yamashita Y: Discrimination of Streptococcus pneumoniae from viridans group streptococci by genomic subtractive hybridization. J Clin Microbiol. 2005, 43 (9): 4528-4534. 10.1128/JCM.43.9.4528-4534.2005.
Whatmore AM, Efstratiou A, Pickerill AP, Broughton K, Woodard G, Sturgeon D, George R, Dowson CG: Genetic relationships between clinical isolates of Streptococcus pneumoniae, Streptococcus oralis, and Streptococcus mitis: characterization of "Atypical" pneumococci and organisms allied to S. mitis harboring S. pneumoniae virulence factor-encoding genes. Infect Immun. 2000, 68 (3): 1374-1382. 10.1128/IAI.68.3.1374-1382.2000.
Mager DL, Ximenez-Fyvie LA, Haffajee AD, Socransky SS: Distribution of selected bacterial species on intraoral surfaces. J Clin Periodontol. 2003, 30 (7): 644-654. 10.1034/j.1600-051X.2003.00376.x.
Whiley RA, Beighton D: Current classification of the oral streptococci. Oral Microbiol Immunol. 1998, 13 (4): 195-216. 10.1111/j.1399-302X.1998.tb00698.x.
Seki M, Yamashita Y, Torigoe H, Tsuda H, Sato S, Maeno M: Loop-mediated isothermal amplification method targeting the lytA gene for detection of Streptococcus pneumoniae. J Clin Microbiol. 2005, 43 (4): 1581-1586. 10.1128/JCM.43.4.1581-1586.2005.
Verhelst R, Kaijalainen T, De Baere T, Verschraegen G, Claeys G, Van Simaey L, De Ganck C, Vaneechoutte M: Comparison of five genotypic techniques for identification of optochin-resistant pneumococcus-like isolates. J Clin Microbiol. 2003, 41 (8): 3521-3525. 10.1128/JCM.41.8.3521-3525.2003.
van Hijum SA, Baerends RJ, Zomer AL, Karsens HA, Martin-Requena V, Trelles O, Kok J, Kuipers OP: Supervised Lowess normalization of comparative genome hybridization data–application to lactococcal strain comparisons. BMC Bioinforma. 2008, 9: 93-10.1186/1471-2105-9-93.
Aguado-Urda M, Lopez-Campos GH, Fernandez-Garayzabal JF, Martin-Sanchez F, Gibello A, Dominguez L, Blanco MM: Analysis of the genome content of Lactococcus garvieae by genomic interspecies microarray hybridization. BMC Microbiol. 2010, 10: 79-10.1186/1471-2180-10-79.
Fukiya S, Mizoguchi H, Tobe T, Mori H: Extensive genomic diversity in pathogenic Escherichia coli and Shigella strains revealed by comparative genomic hybridization microarray. J Bacteriol. 2004, 186 (12): 3911-3921. 10.1128/JB.186.12.3911-3921.2004.
Park HK, Lee HJ, Jeong EG, Shin HS, Kim W: The rgg gene is a specific marker for Streptococcus oralis. J Dent Res. 2010, 89 (11): 1299-1303. 10.1177/0022034510378426.
Park HK, Lee HJ, Kim W: Real-time PCR assays for the detection and quantification of Streptococcus pneumoniae. FEMS Microbiol Lett. 2010, 310 (1): 48-53. 10.1111/j.1574-6968.2010.02044.x.
Park HK, Lee SJ, Yoon JW, Shin JW, Shin HS, Kook JK, Myung SC, Kim W: Identification of the cpsA gene as a specific marker for the discrimination of Streptococcus pneumoniae from viridans group streptococci. J Med Microbiol. 2010, 59 (10): 1146-1152. 10.1099/jmm.0.017798-0.
Shahinas D, Tamber GS, Arya G, Wong A, Lau R, Jamieson F, Ma JH, Alexander DC, Low DE, Pillai DR: Whole-genome sequence of Streptococcus pseudopneumoniae isolate IS7493. J Bacteriol. 2011, 193 (21): 6102-6103. 10.1128/JB.06075-11.
Tettelin H, Nelson KE, Paulsen IT, Eisen JA, Read TD, Peterson S, Heidelberg J, DeBoy RT, Haft DH, Dodson RJ, et al.: Complete genome sequence of a virulent isolate of Streptococcus pneumoniae. Science. 2001, 293 (5529): 498-506. 10.1126/science.1061217.
Gottesman MM, Ambudkar SV: Overview: ABC transporters and human disease. J Bioenerg Biomembr. 2001, 33 (6): 453-458. 10.1023/A:1012866803188.
Sutcliffe IC, Russell RR: Lipoproteins of gram-positive bacteria. J Bacteriol. 1995, 177 (5): 1123-1128.
Macielag MJ, Goldschmidt R: Inhibitors of bacterial two-component signalling systems. Expet Opin Investig Drugs. 2000, 9 (10): 2351-2369. 10.1517/135437188.8.131.521.
Matsushita M, Janda KD: Histidine kinases as targets for new antimicrobial agents. Bioorg Med Chem. 2002, 10 (4): 855-867. 10.1016/S0968-0896(01)00355-8.
Hirakawa H, Nishino K, Hirata T, Yamaguchi A: Comprehensive studies of drug resistance mediated by overexpression of response regulators of two-component signal transduction systems in Escherichia coli. J Bacteriol. 2003, 185 (6): 1851-1856. 10.1128/JB.185.6.1851-1856.2003.
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.
This study was supported by a grant of the Korea Healthcare Technology R&D Project, Ministry for Health & Welfare, Republic of Korea (A085138).
WK and SCM contributed to the design of experiments. HKP implemented experiments and drafted the manuscript. WK analyzed results and edited the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Park, H.K., Myung, S.C. & Kim, W. Comparative transcriptomic analysis of streptococcus pseudopneumoniae with viridans group streptococci. BMC Microbiol 12, 77 (2012). https://doi.org/10.1186/1471-2180-12-77
- Chronic Obstructive Pulmonary Disease
- Histidine Kinase
- Amino Acid Biosynthesis
- KCTC 13048T
- Comparative Transcriptome Analysis