Original Research ARTICLE
In vitro Increased Respiratory Activity of Selected Oral Bacteria May Explain Competitive and Collaborative Interactions in the Oral Microbiome
- 1Center for Microbial Ecology and Technology, Ghent University, Ghent, Belgium
- 2Department of Oral Health Sciences, KU Leuven, Leuven, Belgium
- 3Colgate-Palmolive Technology Center, Piscataway, NJ, United States
Understanding the driving forces behind the shifts in the ecological balance of the oral microbiota will become essential for the future management and treatment of periodontitis. As the use of competitive approaches for modulating bacterial outgrowth is unexplored in the oral ecosystem, our study aimed to investigate both the associations among groups of functional compounds and the impact of individual substrates on selected members of the oral microbiome. We employed the Phenotype Microarray high-throughput technology to analyse the microbial cellular phenotypes of 15 oral bacteria. Multivariate statistical analysis was used to detect respiratory activity triggers and to assess similar metabolic activities. Carbon and nitrogen were relevant for the respiration of health-associated bacteria, explaining competitive interactions when grown in biofilms. Carbon, nitrogen, and peptides tended to decrease the respiratory activity of all pathobionts, but not significantly. None of the evaluated compounds significantly increased activity of pathobionts at both 24 and 48 h. Additionally, metabolite requirements of pathobionts were dissimilar, suggesting that collective modulation of their respiratory activity may be challenging. Flow cytometry indicated that the metabolic activity detected in the Biolog plates may not be a direct result of the number of bacterial cells. In addition, damage to the cell membrane may not influence overall respiratory activity. Our methodology confirmed previously reported competitive and collaborative interactions among bacterial groups, which could be used either as marker of health status or as targets for modulation of the oral environment.
The oral cavity is inhabited by a large and complex microbial community (Aas et al., 2005) Structural, metabolic, and chemical interactions within the oral microbiome are fundamental for maintaining community homeostasis and a healthy microbiome (Hajishengallis et al., 2012). In this way, oral health entails microbiota in symbiosis with the host (Marsh, 2003). The unique structural characteristics of the teeth facilitate bacterial biofilm deposition over time, but excessive accumulation at the gingival margin triggers inflammation (Armitage, 2004). Further disruptions on the microbial homeostasis may promote the expression of community virulence factors from pathobionts (e.g., Porphyromonas gingivalis, Tannerella forsythia, Treponema denticola, and Fusobacterium nucleatum) (Socransky et al., 1998), resulting in dysbioses such as caries or periodontitis (Slomka et al., 2017). Recent studies have focused on exploring the community diversity in subjects with periodontal disease (Griffen et al., 2011), and in the structure shifts from health to periodontitis (Abusleme et al., 2013). However, taxonomic information may provide only partial comprehension of the function of each member within the community. Thus, principles such as substrate utilization must be considered when evaluating the impact of selected members of the oral microbiota on the host (Ramsey and Whiteley, 2009). This increased understanding of the metabolic activities of the oral microbiome may assist the successful management and personalized treatment of oral diseases in the future.
Cell response is a useful tool for characterizing performance variations resulting from environmental influences (Campbell et al., 2013). Metabolic similarity (Anderson et al., 2002) and functional diversity (Zhang et al., 2014) of dental biofilms have been studied, while additional techniques, such as genomics, have described cellular responses involved in periodontitis (Divaris et al., 2013). Decoding the metabolic associations between bacterial species is fundamental for the development of alternative therapies for oral disease. Nevertheless, the individual cell-specific substrate preferences of disease-associated, commensal, and/or health-associated bacteria have not been fully described. Phenomics entails the use of a phenotypic microarray (PM) to study changes in the cellular respiration through the detection of H+ ions, allowing for pathogen identification with either clinical or diagnostic potential (Bochner, 2009). In addition, analyzing metabolic variations between closely related species or between treated and control strains, is possible (Greetham, 2014).
In this study, we focused on investigating the impact of groups of specific molecules on the respiratory activity of selected members of the oral microbiome. Keystone pathogens such as P. gingivalis and T. forsythia are related to dysbiosis (Hajishengallis et al., 2012), while Aggregatibacter actinomycetemcomitans is found in aggressive periodontitis (Jiao et al., 2014) and as opportunistic pathogen (Stacy et al., 2016). Other bacteria, such as Prevotella spp. have increased prevalence on disease states (Hajishengallis and Lamont, 2012). Moreover, Streptococcus mutans is associated with peri-implantitis and caries (Maddi and Scannapieco, 2013). Other species such as F. nucleatum (Jiao et al., 2014), Prevotella intermedia (Maddi and Scannapieco, 2013), and Streptococcus sobrinus (Kilian et al., 2006) may act as putative pathogens. Health-associated bacterial community profiles are not well–understood in caries (Gross et al., 2010). Significant inverse correlations between caries and S. gordonii, but not S. sanguinis were reported (Gross et al., 2010), while other studies show increases in both species in caries vs. health (Peterson et al., 2013). On the contrary, Streptococcus mitis, Streptococcus salivarius, Streptococcus oralis, Streptococcus sanguinis, Veillonella sp., Actinomyces spp., and Capnocytophaga sp. have been reported in healthy individuals or negatively correlated with disease status (Socransky et al., 1998; Aas et al., 2005; Kilian et al., 2006; Maddi and Scannapieco, 2013; Loozen et al., 2014). Furthermore, S. sanguinis and S. gordonii were found to produce H2O2 and inhibit the growth of S. mutans and A. actinomycetemcomitans (Kreth et al., 2008; Sliepen et al., 2008). Commensal populations may prevent adhesion of pathogen species (Teughels et al., 2011). For instance, commensal biofilms containing S. oralis, S. gordonii, S. mitis, and S. sanguinis inhibited the planktonic and biofilm concentrations of P. gingivalis and P. intermedia (Herrero et al., 2016). Finally, species normally associated with dental health include S. mitis, S. oralis, and S. sanguinis (Peterson et al., 2013).
The fact that there is inconsistent definition about the functional roles of different species in health and disease was a rationale to perform this study. We selected the species described above to show a proof of concept of whether respiratory activity can be selectively promoted. We hypothesized that modulating bacterial respiratory metabolism using such functional compounds may enhance the fitness of bacteria performing potentially health-associated activities. We applied a multivariate statistical analysis to detect possible triggers common to several species. Nutrient availability may ultimately impact composition of a polymicrobial community (Hajishengallis and Lamont, 2016). Thus, if the requirements for certain metabolites are shared among the selected species, this may indicate that several species have the capacity to utilize the same metabolite. This approach for analyzing functional metabolic capacities may allow for explaining competitive and collaborative interactions among bacterial groups (Slomka et al., 2017). These individual requirements may certainly vary in live oral biofilms, but they can give us an indication on the metabolic strategies that bacteria with known functions (i.e., pathobionts) employ to persist in the biofilm. In addition, increased understanding of the metabolic capacities of members of the oral microbiome may be useful to highlight individual compounds. These can be used as markers for health status or for developing targets for modulation of the oral ecosystem.
Materials and Methods
Bacterial Cultures and Collection
The following strains commonly present in the oral microbiome (Aas et al., 2005; Maddi and Scannapieco, 2013; Loozen et al., 2014) were obtained from the American Type Culture Collection (ATCC): A. actinomycetemcomitans (ATCC 43718), F. nucleatum (ATCC 10953), P. gingivalis (ATCC 33277), P. intermedia (ATCC 25611), S. mutans (ATCC 25175), S. sobrinus (ATCC 33478), T. forsythia (ATCC 43037), Actinomyces naeslundii (ATCC 51655), Capnocytophaga sputigena (ATCC 33612), Streptococcus gordonii (ATCC 49818), Actinomyces viscosus (ATCC 15987), and S. mitis (ATCC 49456). Veillonella parvula was obtained from the Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures (DSM 2007); S. sanguinis was acquired from the BCCM/LMG Bacteria Collection (LMG 14657). S. salivarius strain TOVE-R (Tanzer et al., 1985) was utilized. Bacterial strains were classified as either commensal (indigenous), potentially health-associated or disease-associated, as indicated in Table 1.
Biolog PM Assays
Biolog phenotype microarrays (PM, Biolog, Inc., Hayward, CA, USA) are a commercial high-throughput system for the analysis of microbial cellular phenotypes. The assay uses the reduction of a tetrazolium dye by NADH as a sensitive indicator of microbial respiration in response to the presence of individual nutrients or chemicals (Bochner et al., 2001). However, this methodology does not allow for estimating population size or bacterial abundance, because only a small portion of all compounds may be actively engaged in cell growth (Posch et al., 1997; Smith and Del Giorgio, 2003). PM plates 1–8, which contained 190 sole carbon sources, 95 sole nitrogen sources, 59 sole phosphate sources, 35 sole sulfur sources, 95 nutrient supplements, and 285 sole peptide-nitrogen sources were used in this analysis. All bacteria, excepting T. forsythia, were grown on blood agar plates (Blood Agar Base No 2, Oxoid, Hampshire, UK) supplemented with 5 μg/ml hemin (Sigma, Diegem, Belgium), 1 μg/ml menadione (Sigma, Diegem, Belgium), and 5% defibrinated horse blood (E&O Laboratories Limited, Bonnybridge, Scotland). T. forsythia was grown brain-heart infusion agar (Oxoid, Hampshire, UK) supplemented with 0.5% (w/v) of yeast extract (Oxoid, Hampshire, UK), 0.1% (w/v) of L-cysteine (Merck, Darmstadt, Germany), 5.0 mg/L of hemin (Sigma, Diegem, Belgium), 0.05 mg/L of menadione (Sigma, Diegem, Belgium), 0.001% (w/v) of N-acetylmuramic acid (Sigma, Diegem, Belgium), 5% Fetal Bovine Serum (Sigma, Diegem, Belgium), and 5% defibrinated horse blood (E&O Laboratories Limited, Bonnybridge, Scotland). A. actinomycetemcomitans (Aa), S. gordonii (Sgord), S. mitis (Smit), S. mutans (Smut), S. salivarius (Ssal), S. sanguinis (Ssang), and S. sobrinus (Ssob) were incubated at 37°C in a 5% CO2 environment. A. naeslundii (Anaes), A. viscosus (Avisc), C. sputigena (Csput), F. nucleatum (Fn), P. gingivalis (Pg), P. intermedia (Pi), T. forsythia (Tf), and V. parvula (Vp) were incubated at 37°C under anaerobic conditions (10% CO2, 10% H2, 80% N2). Anaerobic incubation was performed to recreate the conditions in the periodontal pockets (Marsh et al., 2015). PM assays were performed following the manufacturer's instructions using Biolog Redox Dye Mix D. Metabolic activity was estimated based on the optical density (OD), which was measured after 0, 24, and 48 h of incubation under anaerobic conditions at 37°C at 492 nm using a microplate reader (Multiskan Ascent, Thermo Fisher Scientific Oy, Vantaa, Finland). Incubations were performed in triplicate, under anaerobic conditions (10% CO2, 10% H2, 80% N2), because all bacterial strains were facultative or anaerobic.
Multivariate Statistical Analysis
All the metabolites were categorized according to the respiratory activity triggered on each bacterial species, regarding the control (no nutrient added). Mean OD was calculated for each metabolite, based on the OD readings from the 15 bacterial species. If the OD reading on a species was above the 0.975 quantile, the metabolite was considered to generate high respiratory activity (H-RA) in that strain. Metabolites whose OD value was below the 0.025 quantile were considered to trigger low respiratory activity (L-RA). If the OD was between the 0.025 and 0.975 quantiles, it was considered not significantly different from the control (NS); control OD readings were included within this category (Hernandez-Sanabria et al., 2010). The frequency of the activity triggered by each metabolite within a group of bacteria (commensal, disease-associated, or health-associated) was compared using a Chi-square analysis. In this way, 3 × 2 contingency tables of cross-classification containing the respiratory activity (H, L, or NS) were produced for each of the bacterial groups. When the count of any of the cells was below 5, Fisher Exact test (R Core Team, 2000) was used to calculate the table probability (Hernandez-Sanabria et al., 2013), and activity-associated metabolites were highlighted.
Multiple Correspondence Analysis (MCA) was applied to detect similarities in the respiratory activity of the 15 species. The coordinates best representing bacterial species on two dimensions were retained to create clusters of bacteria with similarities in the use of the PM microarray metabolites. Thus, Hierarchical Clustering of Principal Components (HCPC) was performed, followed by construction of trees using the Ward's method (double factor hierarchical clustering) (Husson et al., 2010). Further, the Calinski-Harabasz (CH) index on k-means clustering was used to evaluate the cluster validity based on the average between- and within- cluster sum of squares, using the vegan package in R (Oksanen et al., 2013). The package NbClust in R was employed to calculate the CH, Silhouette, and Dunn indices (Charrad et al., 2014). Scores were visualized with the package factoextra (Kassambara and Mundt, 2016). Multiple Factor Analysis (MFA) was employed to discriminate particular metabolites triggering respiratory activity (H-RA, L-RA, NS) on each bacterial group at both time points. The function MFA from the FactoMineR package (Lê et al., 2008) was performed in R and the square cosine values (cos2) extracted. Results were plotted in a Principal Coordinate Analysis (PCoA) (Oksanen et al., 2013) and visualized (Wickham, 2016).
Flow Cytometry Analysis
Further proof-of-principle validation was performed in selected species grown in complex medium, to determine whether activity can nevertheless be increased by specific compounds. Cell number and membrane integrity were used as indicators of activity (De Roy et al., 2012). Culturing of F. nucleatum, A. actinomycetemcomitans, and S. salivarius was in modified brain-heart infusion (Alvarez et al., 2013), using the Hungate tube method and maintaining the media in anaerobic conditions (10% CO2, 10% H2, 80% N2). Optical density of the culture used for inoculation was measured after 18 h. Growth analyses on single nutrient sources were performed in 8.0 ml of medium, supplemented with 1% v/v of either L-Asparagine, L-Aspartic Acid, or D-Sorbitol (Sigma, Diegem, Belgium), all at 25mM. Triplicate tubes were inoculated with 100 μl of a 18-h culture diluted 10 × and further incubated anaerobically at 37°C. Flow cytometry was performed at 0, 24, and 48 h. SYBR Green I (SG, Invitrogen, Merelbeke, Belgium) was used alone or in combination with Propidium Iodide (Live/dead BacLight Bacterial Viability kit, Invitrogen, Merelbeke, Belgium), to quantify the total number of bacteria and to differentiate between cells with intact and damaged cytoplasmic membranes, respectively (Berney et al., 2007). Stock staining solutions were prepared as described elsewhere (Van Nevel et al., 2013). Under anaerobic conditions, samples were placed in 96-well plates, diluted 100 × with PBS (pH 7.0) and stained with 2 μl of staining solution and 2 μl of 500 mM Na2EDTA (pH 8.0) for outer membrane permeabilization. Prior to flow cytometry analysis, the stained samples were incubated for 13 min in the dark at 37°C. Cell count analysis was performed using an Accuri C6 flow cytometer (BD Biosciences, Erembodegem, Belgium) equipped with an autoloader, at a flow rate of 66 μl min−1 in 25 μl (Van Nevel et al., 2013). Processing of the flow cytometry results was performed with the Summit v4.3 software. Differences in total and intact number of bacterial cells supplemented with different metabolites were compared using a mixed model with time as a repeated measure, in SAS version 9.3 (Sas, 2012). Statistical significance was assumed at P < 0.05.
Respiratory Activity in Response to Diverse Metabolic Triggers Was Different for Each Bacterial Group
The results of the tetrazolium dye indicator assay were used as estimates of the respiratory activity. Changes regarding the control were classified into High respiratory activity (H-RA), Low respiratory activity (L-RA), or Not Significantly different (NS) from the control using a Chi-square procedure. In this way, 3 × 2 contingency tables of cross-classification containing the respiratory activity (H, L, or NS) were produced for each of the bacterial groups (health-associated, disease-associated, or commensals). The objective of this analysis was to detect common metabolites that promote similar respiratory activity on an entire bacterial group.
Thus, Val-Ala increased the respiratory activity in five of the disease-associated species but the effect was not significantly different from the control for F. nucleatum and in P. gingivalis. Additional 17 compounds tended to decrease the respiratory activity of disease-associated at 24 h (n = 7, P < 0.10, Figure 1A and Table 2A). From these, 11 were carbon sources, one nitrogen and one nutrient source, and four peptides. The respiratory activity of all the commensal species (n = 3), was increased with sources of nitrogen and peptides. Only carbon sources triggered high respiratory activity (P < 0.0001) in all the health-associated (n = 5). At 48 h, 2-aminoethanol increased the respiratory activity of 5 disease-associated bacteria, but not in P. gingivalis and A. actinomycetemcomitans (Figure 1B). Similarly, glucuronamide stimulated respiration in all disease-associated except for F. nucleatum and P. gingivalis. A total of 14 carbon, four nitrogen sources, and three peptides tended to decrease the respiratory activity of all disease-associated, but not significantly (Table 2B). For the commensals, 20 compounds promoted their respiratory activity, mainly peptides (P < 0.0001). Three sources of carbon increased the respiratory activity of all health-associated species (P < 0.0001). Three molecules triggered activity at both 24 and 48 h in commensals, whereas α-hydroxyglutaric acid γ-lactone had similar impact in health-associated. None of the evaluated compounds significantly increased activity in disease-associated at both time points (Figure 1B). Our results indicated that the metabolite requirements of disease-associated bacteria are dissimilar, suggesting that collective modulation of their respiratory activity may be challenging. However, one molecule may consistently support the respiratory activity of all species considered health-associated.
Figure 1. Specific Biolog substrates correlated with high respiratory activity at (A) 24 h and (B) 48 h, in health-associated (n = 5), commensal (n = 3), and disease-associated bacteria (n = 7). Substrates in red significantly increased the respiratory activity in regards to the control (P < 0.05). Blue squares indicate metabolites that significantly decrease the respiratory activity and white squares are for compounds that had no significant effect on the selected bacteria. Table 2 includes the complete description of the metabolites, the significant differences in their effects on the three groups of bacteria and P-values.
Table 2. Metabolites associated with respiratory activity changes in oral bacteria at 24 and 48 h (n = 15).
Functional Similarities among Bacterial Groups Shift over Time
We first focused in detecting the metabolites that collectively increased the respiratory activity on each of the bacterial groups tested, using the frequency analysis. Then, MCA was employed to determine whether bacteria could be grouped by overall metabolic capacities and requirements, regardless of their initial classification (health-associated, disease-associated or commensal). The loadings of the first five components of the MCA were retained and clustered using the k-means hierarchical clustering (Supplementary Table 1). At 24 h, commensals and health-associated were grouped together while a second cluster included five disease-associated, one commensal, and one health-associated (Figure 2A). C. sputigena showed functional similarities with P. gingivalis at 48 h and their activity was significantly different from the bacteria grouped in the other clusters (P < 0.05, Figure 2B). S. mitis, S. sanguinis, and V. parvula were clustered with disease-associated, whereas A. naeslundii, A. viscosus, S. gordonii, and S. salivarius remained together over time. This clustering analysis suggested that there may be overall metabolic similarities among bacterial species. Following, MFA was applied to explain why bacteria considered as commensal or health-associated were grouped with disease-associated. Based on the respiratory activity triggered by all the screened metabolites, health-associated were included in dimensions 3 and 4 at 24 h (Figure 3A). Commensal bacteria were also included in Dim 3, where the highest scores from the metabolic variables corresponded to the carbon and nitrogen sources (Supplementary Table 2A). This result indicates that carbon and nitrogen sources promote respiratory activity in both health-associated and commensals. At 48 h, health-associated bacteria were spread in Dimension 4 and 5, while disease-associated were in Dimension 3 and commensals in Dim. 2. (Figure 3B). The location of health-associated bacteria in the plot, regarding disease-associated, indicated that the differences in the metabolic activity among both groups were triggered by phosphorus and sulfur (Supplementary Table 2B). These findings confirmed that metabolic requirements in health-associated species fluctuate over time but they are shared among the members of a specific bacterial community.
Figure 2. Ward's method of hierarchical unsupervised clustering based on the respiratory activity triggered by substrates contained in the Biolog Phenotype microarray at (A) 24 h and (B) 48 h. The first 2 principal components of the multiple correspondence analysis were retained for cluster analysis. The x-axis of dendrogram represents the linkage distance measured by the Euclidean distance squared. The y-axis indicates the bacterial species included on each cluster. Brackets outline the clusters statistically significant at the 0.05% level and a vertical line indicates this cut-off linkage distance. The groups discussed in the text are colored as follows: health-associated species in white, commensals in gray, and disease-associated bacteria in black circles. Validation metrics of the clustering performance are included in the lower panel. On the left, the inertia explained by each of the clusters; on the center, the Calinski-Harabasz (CH) score indicating the optimal number of clusters; and on the right, the Silhouette scores. CH scores at different number of clusters are described in Supplementary Table 1.
Figure 3. Specific groups of metabolites influenced respiratory activity of each bacterial community. Principal Coordinate Analysis assisted to explain the grouping trends in the cluster analysis, and highlighted similarities in the respiratory activity among the selected oral bacteria at 24 h (A) and 48 h (B). Aa, Aggregatibacter actinomycetemcomitans (ATCC 43718); Fn, Fusobacterium nucleatum (ATCC 10953); Pg, Porphyromonas gingivalis (ATCC 33277); Pi, Prevotella intermedia (ATCC 25611); Sm, Streptococcus mutans (ATCC 25175); Ssob, Streptococcus sobrinus (ATCC 33478); Tf, Tannerella forsythia (ATCC 43037); An, Actinomyces naeslundii (ATCC 51655); Csput, Capnocytophaga sputigena (ATCC 33612); Sgord, Streptococcus gordonii (ATCC 49818); Avisc, Actinomyces viscosus (ATCC 15987); Smitis, Streptococcus mitis (ATCC 49456); Vparv, Veillonella parvula (DSM 2007); Ssang, Streptococcus sanguinis (LMG 14657), and Ssal, Streptococcus salivarius strain TOVE-R.
Potential Competitive Interactions Influence the Fitness of the Oral Ecosystem Members
MFA was additionally used to discriminate metabolites that exclusively triggered respiratory activity in a single group of bacteria, and to uncover correlations among these metabolites. The square correlation ratios (cos2) measure the degree of association between variables and any given dimension of the factor analysis. Thus, the cos2 between the coordinates of the bacterial species and the metabolites revealed that carbon and nitrogen sources were the main metabolic triggers for health-associated bacteria at 24 h 1 (cos2 > 0.8, Supplementary Table 2A). Ala-Gly, urea, D, L-α-Amino-N-Butyric Acid, L-Phenylalanine, nitrate, L-Cysteine, α-Amino-N-Valeric Acid, L-Citrulline, L-Asparagine, D-Lysine, D-Glucosamine, L-Lysine, and ammonia had the highest contributions to dimensions 3 and 4 of the MFA, which both described health-associated bacteria. This suggested that the above were the metabolites that triggered high respiratory activity only in all the health-associated (Supplementary Table 2). After 48 h, carbon and P and S sources were the metabolites that contributed to discriminate between disease-associated and health-associated bacteria. L-Alanyl-Glycine, L-Methionine, Ala-Gly, γ-Hydroxy Butyric Acid, and D-Glucose-1-Phosphate, were described in Dimension 3 with disease-associated, while Glycyl-L-Proline, Glycyl-L-Aspartic Acid, Glycyl-L-Glutamic Acid, Gly-Asn, L-Leucine, gelatine, and L-Alanyl-Glycine, 2-Hydroxyethane Sulfonic Acid, and L-Cysteine Sulfinic Acid, were included in Dimension 4, with health-associated species (Supplementary Table 2B, P < 0.05). L-Valine, L-Isoleucine, and O-Phospho-D-Tyrosine were positively associated with both health-associated and disease-associated, indicating that both groups increase their respiratory activity with these compounds, suggesting potential competitive interactions. L-Glutamine and N-Acetyl-L-Cysteine were negatively associated with disease-associated and Adenosine-3',5'-cyclic monophosphate, and Taurine were partially associated with health-associated, because their cos2 values were high on either dimension 4 or 5 (Supplementary Table 2B). These results explain why collective modulation of respiratory activity may result challenging and thus, alternative strategies may need to be developed.
Metabolic Activity is Not Solely Attributable to Increased Cell Number
Flow cytometry was used to validate whether variations in the respiratory activity were a result of increased cell number. Bacterial cells were quantified following supplementation with compounds that either promoted respiratory activity (L-Asparagine, in health-associated) or that were not found to have effect on the respiratory activity of neither group (D-Sorbitol and L-Aspartic Acid). SYBR Green I (SG) allows for quantification of the total cell number irrespective of viability or cultivability (Zipper et al., 2004), while Propidium iodide (PI) is appropriate for investigating oxidative damage to cells. We used them in combination to differentiate between intact and damaged cells. Significant differences in total cell number in comparison with the control were identified only when L-Aspartic Acid was supplemented (P < 0.05, Table 3). PI allowed for determination of the impact that diverse metabolites may have on the membrane integrity. L-Asparagine and L-Aspartic Acid decreased the number of intact cells in F. nucleatum and S. salivarius at 48 h. However, no significant difference in cell numbers of F. nucleatum was detected with either stain at 24 h. Intact cells of S. salivarius were less when sources of amino acids were added to the medium, while the opposite effect was observed in A. actinomycetemcomitans (P < 0.05). These results indicate that the metabolic activity detected in the Biolog plates may not be a direct result of the number of bacterial cells. In addition, damage to the cell membrane may not influence overall respiratory activity.
As suggested by numerous studies, the oral ecosystem displays synergic activity (Kolenbrander et al., 2010; Tonetti et al., 2011). In the current study, we determined potential substrates that can be used by a wide range of oral bacteria. Multivariate analysis was used to support our hypothesis that specific molecules can selectively stimulate health-associated bacteria. We aimed at explaining whether the metabolic activities could clarify why bacteria considered commensal are frequently accompanying disease-associated bacteria, while some other times they are linked to oral health. Thus, we showed how different bacteria have different responses to such molecules.
Although the respiratory activity was assessed individually, our results may potentially serve as indicators of functional triggers shared by different populations. The commercially available assays of cell-specific metabolic activity in bacteria provide a direct measurement only of the responsive bacterial cells performing respiration. Thus, cell-specific respiration rates may not represent total bacteria abundance (Smith and Del Giorgio, 2003). In fact, a small portion of this abundance may potentially account for the overall activity detected (Posch et al., 1997), which was suggested by the flow cytometry results.
Few substrates increased the respiratory activity of all disease-associated bacteria, indicating that each species may follow different metabolic pathways to generate dysbiosis. Although disease-associated bacteria clustered all together at 24 h, P. gingivalis grouped separately after 48 h. V. parvula was the only health-associated bacteria in the same cluster of the disease-associated, suggesting metabolic similarities, regardless time. Previous reports revealed that A. actinomycetemcomitans grows in the presence of Veillonella sp. and/or F. nucleatum (Periasamy and Kolenbrander, 2009), while F. nucleatum enables growth of P. gingivalis by scavenging O2 (Diaz et al., 2004). Moreover, S. gordonii, Actinomyces sp., Veillonella sp. and A. actinomycetemcomitans are associated with P. gingivalis, because P. gingivalis provides pH support during supragingival colonization (Periasamy and Kolenbrander, 2009). Additionally, our analysis revealed that streptococci clustered in two groups, independently from genetic relatedness (e.g., S. gordonii and S. sanguinis are closely related, while S. mutans and S. sobrinus belong to the same subgroup of streptococci; Xu et al., 2014). Previous reports proposed that the distinction between dental health and disease may be independent of the overall abundance of Streptococcus spp. but rather linked to relative rates of acid production (Peterson et al., 2013). In this way, the fitness of periodontal pathobionts to colonize dental plaque may be shaped by both synergistic and competitive interactions with other members of the oral ecosystem (Jakubovics, 2010). Even though we separately examined the respiratory activities of 15 members of the oral plaque, the results of our clustering analyses presented associations validated in other systems (Kumar and Mason, 2015).
The relevance of carbon and nitrogen sources for the metabolism of health-associated bacteria, inferred in our analysis, (Figure 3) may aid to explain competitive interactions detected when bacteria are grown in biofilms (Jakubovics and Palmer, 2013). For instance, it has been reported that F. nucleatum (Bolstad et al., 1996), P. gingivalis (Mazumdar et al., 2009), and T. forsythia (Cross et al., 2014) prefer amino acids, while carbon is utilized by A. actinomycetemcomitans (Brown and Whiteley, 2007). Peptides are necessary for P. gingivalis, P. intermedia (Takahashi, 2005) and are the second most important energy source for F. nucleatum (Kapatral et al., 2002), while valine, leucine and isoleucine are required for P. intermedia and S. sanguinis (Takahashi, 2005). Our analysis confirmed that L-Leucine increased the respiratory activity of health-associated bacteria (Supplementary Table 2B). P. gingivalis and S. sanguinis have evolved to employ “energetically cheaper” amino acids for protein biosynthesis, as opposed to T. forsythia and P. intermedia (Cross et al., 2014). Consequently, using peptide uptake transport systems may influence the respiratory activity of the health-associated species, and possibly modulate in vivo competitive interactions, as observed with other disease-associated bacteria (Hajishengallis and Lamont, 2016).
The nitrogen released or transferred from deamination or transamination of amino acids results in carbon skeletons that can be used via glucogenic or ketogenic catabolic pathways. In this way, amino acids potentially supported the respiratory activity of health-associated bacteria. Indeed Bacilli (Streptococcus spp.), Fusobacteria, Bacteroidetes (Prevotella, Tannerella spp. and Porphyromonas spp.), Actinobacteria (Actinomyces spp.), Gammaproteobacteria (Aggregatibacter spp.), and Negativicutes (Veillonella spp.) are capable of synthesizing methionine from homoserine (Ferla and Patrick, 2014). In this way, the increase of respiratory activity when health-associated bacteria were exposed to amino acids, may explain competitive interactions observed in oral biofilms in vivo (Slomka et al., 2017). Further, we observed that L-glutamine, a precursor of leucine and glutamic acid, increased the respiratory activity of S. mutans after 48 h. The ability of this species to survive long-term low pH may be associated with the biosynthesis of branched chain amino acids from L-glutamine (Krzyściak et al., 2014). Glucuronamide and ethanolamine (2-aminoethanol) can be used for cellular membrane synthesis (Hendrickson et al., 2014) and as sources of carbon for pathogens (Chiok et al., 2013). The increased respiratory activity of F. nucleatum and P. gingivalis at 48 h when glucuronamide was supplied, suggested that their reported coaggregation interactions (Kolenbrander et al., 2010), may be based not only on localization but on the utilization of common resources and similar metabolic activities.
The boost in the respiratory activity triggered by α-hydroxyglutaric acid-γ-lactone in health-associated may be explained by its role in the pentose phosphate pathway, a major source of reducing equivalents for biosynthesis reactions. The complete genome of S. salivarius JIM8777 in the NCBI database indicated the presence of an enol-lactone hydrolase, confirming its potential ability of using lactones. Moreover, lactonizing enzymes perform “house-cleaning” functions and prevent the accumulation of toxic, non-degradable metabolites (Galperin et al., 2006). The ability of health-associated bacteria for using lactones as a potential source of energy (Mochizuki, 2001), and for membrane functions suggests that health-associated may be more efficient at sensing and communicating variations in the surrounding environment (Kumar and Mason, 2015). Later biofilm formation could be triggered by and include pathways utilizing other molecules pinpointed in our analysis. However, phenotypic variation among strains and genetic relatedness may govern how commensal relationships are ultimately transformed into pathogenic (Kilian et al., 2006; Xu et al., 2014).
Our results suggest that disease-associated bacteria inhabit distinct ecological niches defined not only by physical location and resource utilization (Socransky et al., 1998), but also by their level of metabolic activity (Maurice et al., 2013). Although the increase in respiratory activity triggered by some metabolites may be strain-dependent, we inferred metabolites that might explain the potential mechanisms that oral bacteria employ to adjust their habitat. In this way, it is possible that specific molecules can selectively stimulate health-associated bacteria, as indicated in previous work (Slomka et al., 2017). Combinations of these species in biofilms are likely to buffer the effects of key metabolic substrates identified. For this reason, future efforts to elucidate the roles of specific molecules, using whole-genome sequencing will be necessary for a complete understanding of the observed shifts. Our study provides a framework to potentially identify how the variations in the metabolic activity of individual key players may impact the later biofilm formation. These results may aid to elucidate the functional diversity and metabolic potential of the growing and metabolically active individual members of the oral biofilm. This description of potential functional targets supplies a preliminary outlook that can assist to further develop strategies to manage oral dysbiosis, and improve overall oral health.
Conceived and designed the experiments: VS, WT, and NB. Performed the experiments: VS. Methodology development: EH-S, VS, ERH, WT, and NB. Data mining, statistical analysis, results interpretation, figure, and table preparation: EH-S. Writing the code for frequency analysis in R: FK with the guidance of EH-S. Funding acquisition: NB, WT, and LZ. Contributed reagents/materials/analysis tools: NB, WT. Wrote the paper: EH-S. Reviewed the manuscript: EH-S, VS, ERH, FK, WT, LZ, and NB.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This study was supported by grants from the Research Foundation of Flanders (Fonds Wetenschappelijk Onderzoek-Vlaanderen, FWO), from the Katholieke Universiteit Leuven and from the Colgate-Palmolive Company. We acknowledge Ruben Props and Jana de Bodt for their technical assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fcimb.2017.00235/full#supplementary-material
Aas, J. A., Paster, B. J., Stokes, L. N., Olsen, I., and Dewhirst, F. E. (2005). Defining the normal bacterial flora of the oral cavity. J. Clin. Microbiol. 43, 5721–5732. doi: 10.1128/JCM.43.11.5721-5732.2005
Abusleme, L., Dupuy, A. K., Dutzan, N. Í., Silva, N., Burleson, J. A., Strausbaugh, L. D., et al. (2013). The subgingival microbiome in health and periodontitis and its relationship with community biomass and inflammation. ISME J. 7, 1016–1025. doi: 10.1038/ismej.2012.174
Alvarez, G., Gonzalez, M., Isabal, S., Blanc, V., and Leon, R. (2013). Method to quantify live and dead cells in multi-species oral biofilm by real-time PCR with propidium monoazide. AMB Express 3:1. doi: 10.1186/2191-0855-3-1
Anderson, S. A., Sissons, C. H., Coleman, M. J., and Wong, L. (2002). Application of carbon source utilization patterns to measure the metabolic similarity of complex dental plaque biofilm microcosms. Appl. Environ. Microbiol. 68, 5779–5783. doi: 10.1128/AEM.68.11.5779-5783.2002
Berney, M., Hammes, F., Bosshard, F., Weilenmann, H.-U., and Egli, T. (2007). Assessment and interpretation of bacterial viability by using the LIVE/DEAD BacLight Kit in combination with flow cytometry. Appl. Environ. Microbiol. 73, 3283–3290. doi: 10.1128/AEM.02750-06
Campbell, A. G., Campbell, J. H., Schwientek, P., Woyke, T., Sczyrba, A., Allman, S., et al. (2013). Multiple single-cell genomes provide insight into functions of uncultured Deltaproteobacteria in the human oral cavity. PLoS ONE 8:e59361. doi: 10.1371/journal.pone.0059361
Charrad, M., Ghazzali, N., Boiteau, V., and Niknafs, A. (2014). NbClust: an R Package for determining the relevant number of clusters in a data set. J. Stat. Soft. 61, 1–36. doi: 10.18637/jss.v061.i06
Chiok, K. L., Addwebi, T., Guard, J., and Shah, D. H. (2013). Dimethyl Adenosine Transferase (KsgA) deficiency in Salmonella enterica serovar enteritidis confers susceptibility to high osmolarity and virulence attenuation in chickens. Appl. Environ. Microbiol. 79, 7857–7866. doi: 10.1128/AEM.03040-13
Cross, K. J., Huq, N. L., Park, S. H., Park, J. S., Pham, P., Quah, M., et al. (2014). Bioinformatic investigation of the cost management strategies of five oral microbes. Mol. Oral Microbiol. 30, 87–96. doi: 10.1111/omi.12071
Diaz, P., Zilm, P. S., Wasinger, V., Corthals, G. L., and Rogers, A. H. (2004). Studies on NADH oxidase and alkyl hydroperoxide reductase produced by Porphyromonas gingivalis. Oral Microbiol. Immunol. 19, 137–143. doi: 10.1111/j.0902-0055.2004.00120.x
Divaris, K., Monda, K. L., North, K. E., Olshan, A. F., Reynolds, L. M., Hsueh, W. C., et al. (2013). Exploring the genetic basis of chronic periodontitis: a genome-wide association study. Hum. Mol. Genet. 22, 2312–2324. doi: 10.1093/hmg/ddt065
Griffen, A. L., Beall, C. J., Campbell, J. H., Firestone, N. D., Kumar, P. S., Yang, Z. K., et al. (2011). Distinct and complex bacterial profiles in human periodontitis and health revealed by 16S pyrosequencing. ISME J. 6, 1176–1185. doi: 10.1038/ismej.2011.191
Gross, E. L., Leys, E. J., Gasparovich, S. R., Firestone, N. D., Schwartzbaum, J. A., Janies, D. A., et al. (2010). Bacterial 16S sequence analysis of severe caries in young permanent teeth. J. Clin. Microbiol. 48, 4121–4128. doi: 10.1128/JCM.01232-10
Hajishengallis, G., and Lamont, R. J. (2016). Dancing with the stars: how choreographed bacterial interactions dictate nososymbiocity and give rise to keystone pathogens, accessory pathogens, and pathobionts. Trends Microbiol. 24, 477–489. doi: 10.1016/j.tim.2016.02.010
Hajishengallis, G., and Lamont, R. J. (2012). Beyond the red complex and into more complexity: the polymicrobial synergy and dysbiosis (PSD) model of periodontal disease etiology. Mol. Oral Microbiol. 27, 409–419. doi: 10.1111/j.2041-1014.2012.00663.x
Hendrickson, E. L., Wang, T., Beck, D. A., Dickinson, B. C., Wright, C. J., Lamont, R. J., et al. (2014). Proteomics of Fusobacterium nucleatum within a model developing oral microbial community. Microbiologyopen 3, 729–751. doi: 10.1002/mbo3.204
Hernandez-Sanabria, E., Goonewardene, L. A., Li, M., Mujibi, D. F., Stothard, P., Moore, S. S., et al. (2010). Correlation of particular bacterial PCR-denaturing gradient gel electrophoresis patterns with bovine ruminal fermentation parameters and feed efficiency traits. Appl. Environ. Microbiol. 76, 6338–6350. doi: 10.1128/AEM.01052-10
Hernandez-Sanabria, E., Goonewardene, L. A., Wang, Z., Zhou, M., and Moore, S. S. (2013). Influence of sire breed on the interplay among rumen microbial populations inhabiting the rumen liquid of the progeny in beef cattle. PLoS ONE 8:e58461. doi: 10.1371/journal.pone.0058461
Herrero, E. R., Slomka, V., Bernaerts, K., Boon, N., Hernandez-Sanabria, E., Passoni, B. B., et al. (2016). Antimicrobial effects of commensal oral species are regulated by environmental factors. J. Dent. 47, 23–33. doi: 10.1016/j.jdent.2016.02.007
Kapatral, V., Anderson, I., Ivanova, N., Reznik, G., Los, T., Lykidis, A., et al. (2002). Genome Sequence and analysis of the oral bacterium Fusobacterium nucleatum strain ATCC 25586. J. Bacteriol. 184, 2005–2018. doi: 10.1128/JB.184.7.2005-2018.2002
Kilian, M., Frandsen, E. V., Haubek, D., and Poulsen, K. (2006). The etiology of periodontal disease revisited by population genetic analysis. Periodontol. 2000 42, 158–179. doi: 10.1111/j.1600-0757.2006.00159.x
Kolenbrander, P. E., Palmer, R. J. Jr., Periasamy, S., and Jakubovics, N. S. (2010). Oral multispecies biofilm development and the key role of cell-cell distance. Nat. Rev. Microbiol. 8, 471–480. doi: 10.1038/nrmicro2381
Kreth, J., Zhang, Y., and Herzberg, M. C. (2008). Streptococcal antagonism in oral biofilms: Streptococcus sanguinis and Streptococcus gordonii interference with Streptococcus mutans. J. Bacteriol. 190, 4632–4640. doi: 10.1128/JB.00276-08
Krzyściak, W., Jurczak, A., Kościelniak, D., Bystrowska, B., and Skalniak, A. (2014). The virulence of Streptococcus mutans and the ability to form biofilms. Eur. J. Clin. Microbiol. Infect. Dis. 33, 499–515. doi: 10.1007/s10096-013-1993-7
Loozen, G., Ozcelik, O., Boon, N., De Mol, A., Schoen, C., Quirynen, M., et al. (2014). Inter-bacterial correlations in subgingival biofilms: a large-scale survey. J. Clin. Periodontol. 41, 1–10. doi: 10.1111/jcpe.12167
Maddi, A., and Scannapieco, F. (2013). Oral biofilms, oral and periodontal infections, and systemic disease. Am. J. Dent. 26, 249–254. Available online at: https://www.researchgate.net/publication/259986828_Oral_biofilms_oral_and_periodontal_infections_and_systemic_disease
Mochizuki, K. (2001). Purification and characterization of a lactonase from Burkholderia sp. R-711, that hydrolyzes (R)-5-oxo-2-tetrahydrofurancarboxylic acid. Arch. Microbiol. 175, 430–434. doi: 10.1007/s002030100283
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2013). Package ‘vegan’. Community ecology package, version 2. Available online at: http://cran.r-project.org/web/packages/vegan/index.html (Accessed April 21, 2017).
Periasamy, S., and Kolenbrander, P. E. (2009). Aggregatibacter actinomycetemcomitans builds mutualistic biofilm communities with Fusobacterium nucleatum and Veillonella species in saliva. Infect. Immun. 77, 3542–3551. doi: 10.1128/IAI.00345-09
Posch, T., Pernthaler, J., Alfreider, A., and Psenner, R. (1997). Cell-specific respiratory activity of aquatic bacteria studied with the tetrazolium reduction method, cyto-clear slides, and image analysis. Appl. Environ. Microbiol. 63, 867–873.
Ramsey, M. M., and Whiteley, M. (2009). Polymicrobial interactions stimulate resistance to host innate immunity through metabolite perception. Proc. Natl. Acad. Sci. U.S.A. 106, 1578–1583. doi: 10.1073/pnas.0809533106
Sliepen, I., Hofkens, J., Van Essche, M., Quirynen, M., and Teughels, W. (2008). Aggregatibacter actinomycetemcomitans adhesion inhibited in a flow cell. Oral Microbiol. Immunol. 23, 520–524. doi: 10.1111/j.1399-302X.2008.00456.x
Slomka, V., Hernandez-Sanabria, E., Herrero, E. R., Zaidel, L., Bernaerts, K., Boon, N., et al. (2017). Nutritional stimulation of commensal oral bacteria suppresses pathogens: the prebiotic concept. J. Clin. Periodontol. 44, 344–352. doi: 10.1111/jcpe.12700
Stacy, A., Fleming, D., Lamont, R. J., Rumbaugh, K. P., and Whiteley, M. (2016). A commensal bacterium promotes virulence of an opportunistic pathogen via cross-respiration. MBio 7:e00782-16. doi: 10.1128/mBio.00782-16
Takahashi, N. (2005). Microbial ecosystem in the oral cavity: Metabolic diversity in an ecological niche and its relationship with oral diseases. Int. Congress Ser. 1284, 103–112. doi: 10.1016/j.ics.2005.06.071
Teughels, W., Loozen, G., and Quirynen, M. (2011). Do probiotics offer opportunities to manipulate the periodontal oral microbiota? J. Clin. Periodontol. 38(Suppl. 11), 159–177. doi: 10.1111/j.1600-051X.2010.01665.x
Tonetti, M. S., and Chapple, I. L., on Behalf of Working Group 3 of the Seventh European Workshop on Periodontology (2011). Biological approaches to the development of novel periodontal therapies–consensus of the Seventh European Workshop on Periodontology. J. Clin. Periodontol. 38(Suppl. 11), 114–118. doi: 10.1111/j.1600-051X.2010.01675.x
Van Nevel, S., Koetzsch, S., Weilenmann, H. U., Boon, N., and Hammes, F. (2013). Routine bacterial analysis with automated flow cytometry. J. Microbiol. Methods 94, 73–76. doi: 10.1016/j.mimet.2013.05.007
Xu, H., Jenkinson, H. F., and Dongari-Bagtzoglou, A. (2014). Innocent until proven guilty: mechanisms and roles of Streptococcus–Candida interactions in oral health and disease. Mol. Oral Microbiol. 29, 99–116. doi: 10.1111/omi.12049
Zhang, Y., Zheng, Y., Hu, J., Du, N., and Chen, F. (2014). Functional diversity of the microbial community in healthy subjects and periodontitis patients based on sole carbon source utilization. PLoS ONE 9:e91977. doi: 10.1371/journal.pone.0091977
Zipper, H., Brunner, H., Bernhagen, J. R., and Vitzthum, F. (2004). Investigations on DNA intercalation and surface binding by SYBR Green I, its structure determination and methodological implications. Nucleic Acids Res. 32:e103. doi: 10.1093/nar/gnh101
Keywords: oral pathobionts, health-associated bacteria, commensal bacteria, phenotype microarray, respiratory activity, multivariate statistical analysis
Citation: Hernandez-Sanabria E, Slomka V, Herrero ER, Kerckhof F-M, Zaidel L, Teughels W and Boon N (2017) In vitro Increased Respiratory Activity of Selected Oral Bacteria May Explain Competitive and Collaborative Interactions in the Oral Microbiome. Front. Cell. Infect. Microbiol. 7:235. doi: 10.3389/fcimb.2017.00235
Received: 06 December 2016; Accepted: 22 May 2017;
Published: 07 June 2017.
Edited by:Lorenza Putignani, Bambino Gesù Ospedale Pediatrico (IRCCS), Italy
Reviewed by:Valerio Iebba, Sapienza Università di Roma, Italy
Nick Stephen Jakubovics, Newcastle University, United Kingdom
Parish Paymon Sedghizadeh, University of Southern California, United States
Copyright © 2017 Hernandez-Sanabria, Slomka, Herrero, Kerckhof, Zaidel, Teughels and Boon. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nico Boon, nico.boon@UGent.be
†These authors have contributed equally to this work.