- Open Access
Chemotherapy-induced oral mucositis is associated with detrimental bacterial dysbiosis
Microbiomevolume 7, Article number: 66 (2019)
Gastrointestinal mucosal injury (mucositis), commonly affecting the oral cavity, is a clinically significant yet incompletely understood complication of cancer chemotherapy. Although antineoplastic cytotoxicity constitutes the primary injury trigger, the interaction of oral microbial commensals with mucosal tissues could modify the response. It is not clear, however, whether chemotherapy and its associated treatments affect oral microbial communities disrupting the homeostatic balance between resident microorganisms and the adjacent mucosa and if such alterations are associated with mucositis. To gain knowledge on the pathophysiology of oral mucositis, 49 subjects receiving 5-fluorouracil (5-FU) or doxorubicin-based chemotherapy were evaluated longitudinally during one cycle, assessing clinical outcomes, bacterial and fungal oral microbiome changes, and epithelial transcriptome responses. As a control for microbiome stability, 30 non-cancer subjects were longitudinally assessed. Through complementary in vitro assays, we also evaluated the antibacterial potential of 5-FU on oral microorganisms and the interaction of commensals with oral epithelial tissues.
Oral mucositis severity was associated with 5-FU, increased salivary flow, and higher oral granulocyte counts. The oral bacteriome was disrupted during chemotherapy and while antibiotic and acid inhibitor intake contributed to these changes, bacteriome disruptions were also correlated with antineoplastics and independently and strongly associated with oral mucositis severity. Mucositis-associated bacteriome shifts included depletion of common health-associated commensals from the genera Streptococcus, Actinomyces, Gemella, Granulicatella, and Veillonella and enrichment of Gram-negative bacteria such as Fusobacterium nucleatum and Prevotella oris. Shifts could not be explained by a direct antibacterial effect of 5-FU, but rather resembled the inflammation-associated dysbiotic shifts seen in other oral conditions. Epithelial transcriptional responses during chemotherapy included upregulation of genes involved in innate immunity and apoptosis. Using a multilayer epithelial construct, we show mucositis-associated dysbiotic shifts may contribute to aggravate mucosal damage since the mucositis-depleted Streptococcus salivarius was tolerated as a commensal, while the mucositis-enriched F. nucleatum displayed pro-inflammatory and pro-apoptotic capacity.
Altogether, our work reveals that chemotherapy-induced oral mucositis is associated with bacterial dysbiosis and demonstrates the potential for dysbiotic shifts to aggravate antineoplastic-induced epithelial injury. These findings suggest that control of oral bacterial dysbiosis could represent a novel preventive approach to ameliorate oral mucositis.
One of the most frequent complications of chemotherapy is oral mucositis, reported to affect about 75% of patients receiving high-dose conditioning chemotherapy prior to hematopoietic cell transplantation and 20 to 60% of individuals treated for solid tumors [1, 2]. Lesions present as erythema and ulceration of the non-keratinized mucosa. Although the condition is self-limiting, it can impact the delivery of optimal cancer treatment, since it is associated with clinically significant pain, compromised nutrition, prolonged hospitalization, bloodstream infections and antineoplastic dose reductions [3, 4]. Different approaches to prevent or treat oral mucositis have been largely ineffective as there is insufficient knowledge on its pathophysiology .
Oral mucositis primarily results from the cytotoxic effects of chemotherapeutics on the rapidly dividing oral epithelium, with other areas along the gastrointestinal tract also commonly affected. Antineoplastics exert their cytotoxic effects by a variety of mechanisms leading to impaired DNA replication and repair, cell-cycle arrest, DNA damage, and cell death [6, 7]. However, the specific downstream cellular events mediating oral epithelial damage remain poorly characterized. It has been proposed that oral microbiome communities, which live in constant cross-talk with the adjacent mucosal tissues, could contribute to the development of mucositis [8, 9]. Evidence from animal models of intestinal mucositis suggests it is plausible that epithelial-microbiome cross-talk along the gastrointestinal track modulates susceptibility to mucositis, with resident microbial commensals shown to be necessary for irinotecan-dependent intestinal mucosal injury . Also, mice deficient in the innate immunity activator Toll-like receptor 2 show increased susceptibility to methotrexate-induced intestinal mucositis, an indication that microbial signaling could modify the severity of the response .
It is not completely clear how chemotherapy alters the oral environment. Various studies in humans using microbiological techniques with limited power of detection indicate the oral microbiota changes during cancer treatment [12, 13]. However, there is a paucity of longitudinal, well-controlled studies that use highly sensitive high throughput sequencing to characterize the oral microbiome during chemotherapy. Disruption of microbiome communities, allowing growth of pathobionts, could negatively impact the ability of mucosal tissues to remain intact during an antineoplastic challenge. The plausibility of this hypothesis was demonstrated in an animal model of intestinal mucositis, in which the antineoplastic cisplatin induced significant changes in the intestinal microbiome, but restoration of the gut microbiota through fecal-pellet gavage promoted healing of cisplatin-induced mucositis . It is possible that similar events take place in the oral cavity, where chemotherapy could induce microbiome disruptions that affect the susceptibility of oral tissues to mucositis.
Alterations in the oral microbiome during chemotherapy could result from impairment of the different constituents involved in maintaining these communities in a commensal state. It has been suggested that chemotherapy compromises the flow of saliva, one of the most important defense mechanisms in the oral cavity [15, 16]. Chemotherapy-associated myelosuppression could also compromise the availability of oral neutrophils, resulting in a more permissive environment for detrimental oral microorganisms. Other chemotherapy-associated treatments such as antibiotics could also affect the oral microbiome. In addition, antineoplastics have been shown to have antimicrobial activity and could potentially contribute to microbiome changes during chemotherapy [17, 18].
Here, we present results of a comprehensive longitudinal evaluation of subjects undergoing 5-fluorouracil (5-FU) or doxorubicin-based chemotherapy, characterizing changes in the bacterial and fungal oral microbiome, salivary flow rate (SFR), and oral granulocyte availability during one treatment cycle and in relation to the development of oral mucositis. As a control for microbiome stability, non-cancer subjects were also longitudinally sampled. We found that oral mucositis correlated with the dose of 5-FU, docetaxel, and cisplatin, increased salivary flow, and greater oral granulocyte presence. Chemotherapy disrupted the oral bacteriome with these shifts showing a strong correlation with oral mucositis severity. Shifts were not specific to mucositis lesions but rather affected the oral cavity as a whole. Microbiome shifts associated with mucositis were independent from the effect of antibiotics and acid inhibitors. Although shifts correlated with antineoplastic doses, and 5-FU showed potential to modify directly the microbiome, mucositis-associated shifts were unlikely the result of a selective antibacterial action of 5-FU, as commensals depleted during mucositis were resistant to physiologically relevant drug concentrations and Fusobacterium nucleatum, a species enriched during mucositis, was sensitive to 5-FU. Bacteria with a longitudinal change in abundance that negatively correlated with mucositis severity included species of Streptococcus, Actinomyces, Gemella, Granulicatella, and Veillonella, all common in oral health, while more severe mucositis correlated with enrichment of Gram-negatives previously associated with other oral inflammatory conditions . These results are consistent with mucositis-associated dysbiosis as a result of alterations in nutritional resources due to inflammation, echoing other oral conditions in which inflammatory products promote pathobiont growth, which in turn alters host responses . To better understand mucosal responses during chemotherapy, we longitudinally characterized oral epithelial transcriptome changes in a subset of subjects, finding genes involved in the innate immune response and apoptosis, among others, to be upregulated. The potential for a dysbiotic microbiome to modify mucosal responses was evaluated in a 3D oral mucosa multilayered construct. In this model, we show that Streptococcus salivarius, a potential symbiont depleted during mucositis, is tolerated by the oral mucosa. In contrast, Fusobacterium nucleatum subsp. vincentii, associated with more severe mucositis, has pro-inflammatory and pro-apoptotic effects on oral epithelial cells. To our knowledge, this study provides the first comprehensive characterization of oral mucositis pathophysiology in humans demonstrating bacterial dysbiosis is associated with lesion severity, with such bacteriome changes showing clear potential to contribute to mucositis by promoting further epithelial damage.
Oral mucositis is associated with 5-FU exposure
To better understand the pathophysiology of oral mucositis we followed 49 subjects during one cycle of 5-FU- or doxorubicin-based chemotherapy (Additional file 1: Figure S1 and Additional file 2: Tables S1 and S2) and recorded the incidence and severity of oral mucositis. Fig. 1a shows a typical oral mucositis clinical presentation with erythema and ulceration affecting mostly the non-keratinized mucosa. As seen in Fig. 1b, mucositis occurred in 65 to 76% of subjects as measured via the WHO and OMAS scales, respectively, with mucositis severity peaking at V3 (Fig. 1c, d). Mucositis incidence was higher in subjects on the 5-FU-based regimen in comparison to those taking doxorubicin (Fig. 1e, WHO p = 0.016 and OMAS p = 0.019), and a moderate but highly significant positive correlation was seen between the dose of 5-FU and mucositis severity (Fig. 1f). The doses of docetaxel and cisplatin, which were commonly given in combination with 5-FU (Additional file 2: Table S2), also positively correlated with more severe mucositis; while other drugs given with 5-FU such as carboplatin, oxaliplatin and the chemoprotectant leucovorin showed no correlation or a negative correlation with mucositis severity (Fig. 1f).
Oral mucositis is associated with increased salivary flow and greater oral neutrophil presence
A longitudinal analysis was next conducted to better understand if disruptions in two important oral homeostatic mechanisms, saliva and neutrophil surveillance, occurred during chemotherapy and were associated with mucositis. Changes in these variables were evaluated by comparing values at each visit to baseline and by modelling the change over four visits with linear or quadratic polynomial contrasts. In contrast to previous reports , we found that SFR increased during chemotherapy (Fig. 2b, p = 0.06 at V3; p = 0.028 at V4; and p = 0.002 for linear change). Increased SFR correlated with more severe mucositis, as shown by significant correlations between the linear increase in SFR during the cycle and the change in mucositis symptoms modelled either quadratically (Fig. 2c, p = 0.001) or linearly (Additional file 2: Table S3, p = 0.003). These findings suggested mucositis was not associated with decreased saliva, but instead, mucosal inflammation and ulceration occurred concomitant with and/or were followed by increased SFR.
Since neutrophil surveillance could affect oral homeostasis by controlling the growth of oral commensals, which could in turn affect mucosal health, we evaluated the effect of chemotherapy on blood and oral neutrophil counts. As expected due to antineoplastic cytotoxicity, blood neutrophil counts in the cancer group decreased at V3 in comparison to baseline levels (p = 0.00009); however, at V4 neutrophil counts rebounded (Fig. 2e). Decreased blood neutrophils in the middle of the cycle (positive quadratic change) correlated with a linear increase in mucositis severity (Fig. 2f, p = 0.001 and Additional file 2: Table S3, p = 0.001). Oral neutrophils enter the oral cavity through the thin junctional epithelium at the gingiva/tooth interface and via transmucosal migration . Although chemotherapy decreased blood neutrophils, oral granulocytes did not follow a similar pattern, only decreasing slightly at V2 (Fig. 2h, p = 0.015). A linear change in oral granulocytes positively correlated with a linear change in mucositis severity (Fig. 2i, p = 0.007). Taken together, these findings indicate that despite decreased blood neutrophils as a consequence of chemotherapy, these cells were still being recruited into the oral tissues, with mucositis lesions promoting increased efflux of granulocytes into the oral space.
Additional file 2: Table S3 shows other demographic and clinical characteristics that correlated with the longitudinal clinical progression of oral mucositis. The negative quadratic change in mucositis severity (low, high, low) was associated with 5-FU (p = 0.004) but also with number of prosthetic teeth (p = 0.003) and wearing a removable oral prosthesis (p = 0.002). A linear change in mucositis severity, indicating subjects with more severe mucositis late in the cycle (delayed healing), positively correlated with smoking (p = 0.007), steroid intake (p = 0.002), and cisplatin dose (p = 0.002), while it negatively correlated with doxorubicin dose (p = 0.00001).
Disruption of salivary bacterial communities during chemotherapy correlates with mucositis severity
The oral microbiome was longitudinally evaluated by sequencing of salivary and mucosal 16S rRNA gene amplicons and salivary internal transcribed spacer (ITS)-1 DNA amplicons. Chemotherapy disrupted the oral microbiome with the most severe changes seen in salivary bacterial communities. Salivary bacterial alpha-diversity decreased during the cycle (p = 0.01 at V3 and p = 0.009 for linear change) and at the visit with the most severe mucositis (Fig. 3a, p = 0.001). Mucosal bacterial community diversity was only minimally affected, with a small increase seen at V2 (p = 0.042), while salivary fungal community diversity did not change during chemotherapy (Fig. 3a). Changes over time in salivary bacterial diversity modelled with orthogonal polynomial contrasts showed a correlation of the linear change with the dose of 5-FU and intake of a multi-dose antibiotic, and of the quadratic change with docetaxel and cisplatin doses (Additional file 2: Table S3). To better understand the contributions of different predictors to changes in salivary bacterial diversity and evaluate their effect at each time point during the cycle, we also conducted a by visit analysis, calculating the change in diversity from each visit to baseline (Fig. 3b). This analysis showed that at V2, there was a trend for a correlation between decreased salivary bacterial diversity with higher cisplatin dose and intake of acid inhibitors. At V3, decreased diversity correlated with higher OMAS scores at the same visit (p = 0.003), a higher docetaxel dose (p = 0.008), and use of acid inhibitors (p = 0.004). The use of a multi-dose antibiotic and the dose of cisplatin showed a trend for a negative correlation with diversity. At V4, the use of a multi-dose antibiotic (p = 0.0005) was associated with a decrease in diversity. A linear regression analysis was run considering the change in diversity at V3 as dependent variable and OMAS score, the use of multi-dose antibiotic and use of acid inhibitors as predictors. This analysis showed that OMAS (β = − 0.300, p = 0.046) and acid inhibitor use (β = − 0.301, p = 0.045), but not antibiotic use, remained significantly associated with the change in diversity at V3. In summary, different factors appeared to have influenced salivary bacterial alpha-diversity with changes at V3 associated with mucositis severity and acid inhibitor use and at V4 with the intake of a multi-dose antibiotic.
Echoing the changes seen in alpha-diversity, salivary bacterial community structure (beta-diversity) was also disturbed during chemotherapy (Fig. 3c). Mucosal bacterial communities followed a similar trend as their salivary counterparts but their changes were not greater than those of non-cancer controls. The structure of salivary fungal communities was not disturbed by chemotherapy (Fig. 3c).
To inquire which patient factors were associated with disruptions in community structure, the ThetaYC distances between baseline to each visit, indicating the magnitude of change, were evaluated for a correlation with clinical characteristics (Fig. 3d). This analysis showed that the degree of salivary bacterial community disruption at V2 correlated with cisplatin dose (p = 0.004). At V3, salivary bacterial ThetaYC distances positively correlated with the degree of mucositis severity (p = 0.000005), with the doses of 5-FU (p = 0.00007) and docetaxel (p = 0.0001) and with the use of acid inhibitors (p = 0.002). Linear regression using salivary bacterial ThetaYC distances at V3 as the dependent variable and OMAS scores and acid inhibitor use as predictors showed both variables remained independently associated with community structure change (OMAS β = 0.527, p = 0.0002 and acid inhibitor β = 0.287, p = 0.033). At V4, salivary bacterial community disruption correlated with mucositis severity (p = 0.002) and with 5-FU (p = 0.005) and docetaxel (p = 0.007) doses. Surprisingly, antibiotic use only showed a trend for a correlation with salivary bacterial community disruption, and this occurred only at V4 (p values below 0.05 but not significant after adjustment for multiple testing). A linear regression analysis evaluating the effect of OMAS, multi-dose antibiotics, and acid inhibitors on ThetaYC distances at V4 showed that when the effect of each of the three predictors was adjusted for the effect of the others, OMAS (β = 0.347, p = 0.014) and antibiotics (β = 0.281, p = 0.048) remained significant below a 0.05 p level. For mucosal communities, the only significant correlation seen was between ThetaYC distances at V4 and antibiotic use. A correlation trend was seen between mucosal ThetaYC distances at V3 with cisplatin dose and antibiotics. Overall, these results show that mucositis severity and acid inhibitor use were the main predictors associated with salivary bacterial community structure disruption during the cycle. Antineoplastic doses were also associated with ThetaYC distances, but it should be noted that their independent association could not be evaluated as they are highly correlated with OMAS scores. Antibiotics seemed to only affect the microbiome at V4 having a greater effect on mucosal communities than on their salivary counterparts.
We next asked if disruption of salivary bacterial community structure occurred in a consistent direction among subjects as visualized by PCoA plots (Additional file 1: Figure S2). At V2, there was no significant separation of samples from baseline communities. At V3 and V4, there was a significant separation of the corresponding data clouds from baseline samples, although not all microbiomes changed in a uniform direction (Additional file 1: Figure S2A–B, p = 0.006 at V3 and p = 0.003 at V4). Although microbiome changes were heterogeneous among subjects, microbiome variability along axis 2 in both the V3 and V4 plots (Additional file 1: Figure S2A–B) correlated with OMAS score, doses of 5-FU, and docetaxel, with receiving a multi-dose antibiotic prior to V3, and with changes in alpha-diversity.
Changes in the proportions of individual taxa during chemotherapy were also evaluated. Additional file 1: Figure S2C shows the number of taxa that significantly changed at each visit, with the greatest number of altered taxa seen in salivary bacterial communities at V3 and V4. Additional file 1: Figure S3 and Additional file 3: Table S4 show the identity of these taxa. At V3, all changes in saliva were due to significant decreases in relative abundance. Several Gram-negative species from the genera Prevotella, Selenomonas, Leptotrichia, Tannerella, Campylobacter, Capnocytophaga, and Parvimonas showed a trend to increase in relative abundance at V3, but no species reached significance, highlighting the heterogeneity among individuals in the species that replaced those that decreased. At V4, we observed both significant increases and decreases in abundance.
Taken together, these results show that the salivary bacterial microbiome was profoundly altered during chemotherapy. While intake of multi-dose antibiotics and acid inhibitor use during the cycle were associated with some of the observed microbiome modifications, disruption of salivary bacterial communities was also independently associated with oral mucositis severity. The total doses of 5-FU, docetaxel, and cisplatin, associated with more severe mucositis, also correlated with greater salivary bacteriome disruption.
Oral mucositis is associated with oral bacteriome dysbiosis
To evaluate the longitudinal covariation of microbial taxa and clinical variables during chemotherapy, we conducted a multilevel multivariate sparse partial least square analysis focusing on bacterial species enriched or depleted in accordance with mucositis severity. Figure 4a shows the correlations between the longitudinal change in salivary bacterial proportions and clinical variables, including OMAS scores. Mucositis severity positively correlated with enrichment of 3 salivary Gram-negative species (taxa 1 to 3 in Fig. 4a), namely Fusobacterium nucleatum subsp. vincentii, an uncultured Clostridiales, and Treponema maltophilum, all taxa associated with other oral inflammatory conditions . Mucositis severity also correlated with a decrease in proportions of 24 salivary bacterial species (taxa 4 to 27 in Fig. 4a), among them abundant commensals from genera typically associated with oral health such as Streptococcus, Actinomyces, Gemella, Granulicatella, and Veillonella [19, 22]. Consistent with previous analyses (Fig. 2), mucositis severity positively correlated with SFR and oral neutrophils (Fig. 4a). It should also be noted that the intake of a multi-dose antibiotic did not correlate with changes in the same salivary bacteria associated with OMAS. This result further confirmed antibiotic treatment was not the cause of mucositis-associated dysbiosis.
A similar analysis of mucosal bacterial communities found that changes in only a few taxa correlated with mucositis severity (Fig. 4b). Changes, however, were consistent with those seen in saliva, with enrichment of the Gram-negative, inflammation-associated Prevotella oris  and lower abundance of Streptococcus, Veillonella, and Actinomyces species as OMAS scores increased (Fig. 4b).
Taken together, these results revealed that mucositis is associated with depletion of prevalent health-associated commensals (likely symbionts) and enrichment of Gram-negative species typically enriched during oral inflammation (likely pathobionts). These changes are reminiscent of dysbiotic shifts seen in other oral inflammatory conditions such as gingivitis and periodontitis , although the changes are not directly related to these conditions as mucositis did not correlate with the presence of severe periodontitis.
The antimicrobial activity of antineoplastics is unlikely to explain chemotherapy-associated dysbiosis
We found that the change in relative abundance of many species associated with mucositis (in Fig. 4a, b) also correlated with the use and dose of specific antineoplastics (Additional file 1: Figure S4A). For instance, depletion of Streptococcus salivarius during mucositis correlated with higher doses of 5-FU and docetaxel (p = 0.004 in both cases). Additional file 1: Figures S4B and S4C show taxa depleted during mucositis in subjects taking 5-FU, but that did not change in subjects on doxorubicin who also presented with mucositis, albeit of less severity. These results raised the possibility that the antimicrobial activity of specific antineoplastics was contributing to the selective depletion or enrichment of oral bacteria during chemotherapy-induced mucositis. Indeed, when we exposed salivary communities to physiological concentrations of 5-FU or 5-FU and docetaxel, we observed a decrease in viability of the total culturable community after a 2-h incubation suggesting these drugs have the potential to modify the microbiome (Additional file 1: Figure S4D). We hypothesized that if 5-FU was responsible for the dysbiotic changes seen during chemotherapy, in particular those shifts associated with mucositis, then depleted taxa should be susceptible to the drug while enriched species would be resistant. To test this, bacterial strains of the depleted S. salivarius, Streptococcus parasanguinis, Veillonella atypica, and Veillonella rogosae and the enriched Fusobacterium nucleatum were exposed to 5-FU. After a 2-h incubation, however, only Veillonella atypica showed decreased viability at the same physiologic concentration (7.7 μM) tested with whole saliva (Additional file 1: Figure S4E). In fact, many of the commensals depleted during mucositis were resistant to killing by 5-FU even at drug concentrations likely to be above the physiologic range (e.g., 770 μM). The ability of 5-FU to inhibit the growth of S. salivarius and F. nucleatum was also tested with these evaluations showing that even high concentrations of 5-FU did not completely inhibit S. salivarius, while the growth of F. nucleatum was inhibited by all concentrations tested (Additional file 1: Figure S4F–G). Since taxa negatively correlated with 5-FU and mucositis were not consistently sensitive to 5-FU exposure and a species enriched during mucositis was not resistant to the drug, we conclude that an antimicrobial effect of 5-FU is unlikely to explain the microbiome dysbiosis associated with mucositis. Chemotherapeutics, however, could still temporarily modify the oral microbiome, especially while salivary drug concentrations are in the high range, but these changes do not seem to be associated with mucositis-specific shifts.
Mucosal lesions did not differ in their microbiome from healthy mucosa
The analysis of mucosal communities used combined libraries from sites with different clinical appearance (healthy, erythematous, ulcerated) obtained from the same subject. An analysis using separate libraries (see list of libraries in Additional file 2: Table S5) was also conducted to explore whether mucosal lesions harbored a microbiome distinct from their subject-matched healthy site counterparts. We did not observe differences in either alpha-diversity (Additional file 2: Table S6) or in the abundance of individual taxa when comparing visit-matched, within-subject, mucosal communities with different clinical appearance (healthy vs erythematous or healthy vs ulcerated). This result indicates that mucosal surface changes were not selecting for a specific microbiome colonization at the lesion sites but rather mucositis-associated dysbiosis was affecting the mouth as a whole.
The oral epithelium responds to chemotherapy by upregulating genes involved in innate immune responses and apoptosis
Since little information exists on genes and pathways mediating oral mucosal injury, we conducted a gene expression analysis in a subset of 14 subjects comparing oral epithelial transcript levels between baseline and V3. The clinical characteristics of the subjects analyzed are shown in Additional file 2: Table S7. SAM analysis showed 158 upregulated and two downregulated genes at V3 (> 1.5-fold) (Additional file 4: Table S8). Although there was substantial variability in the response of subjects, this was not related to the type of chemotherapeutic regimens subjects were taking according to clustering analysis. Figure 5a shows a summary of upregulated gene ontology (GO) term categories (for a complete list see Additional file 5: Table S9), which included genes involved in epidermis development, apoptosis and cell death, proteolysis, immune responses, and response to stress. Figure 5b shows the main genes involved in innate immune responses that were upregulated. These included TNF, which is a key inflammatory mediator that leads to activation of nuclear factor κB (NF-κB) but could also activate apoptosis via the extrinsic pathway or induce necroptosis . Other genes involved in epithelial innate immune responses to microbial commensals such as IL17C, CCL20, CXCL2, DEFB4A, and DEFB103A were also upregulated suggesting a microbiome influence on mucosal responses during chemotherapy. Another category of genes upregulated that could have a role in mediating mucosal damage was genes related to apoptosis (Fig. 5c). Of particular interest was the upregulation of PMAIP1, the gene encoding the proapoptotic Bcl2 homology 3 (BH3)-only protein NOXA, which activates the intrinsic apoptosis pathway , and has been shown to mediate cisplatin-induced cell death .
Potential for bacterial dysbiotic changes to affect epithelial responses
We next explored the potential for dysbiotic microbiome changes associated with mucositis severity to affect epithelial responses. As Fig. 5d shows, exposure of a 3D multilayer oral epithelial construct to the potential symbiont S. salivarius, depleted during mucositis, provoked minimal induction of innate immune response genes. In contrast, challenge of epithelial tissues with a similar load of the potential pathobiont F. nucleatum, associated with more severe mucositis, produced a dramatic upregulation of TNF, IL17C, CCL20, and CXCL2. We also observed that F. nucleatum induced an upregulation of PMAIP1 and, consistent with this effect, tissues displayed death of cells in the apical layer (Fig. 5e). These results suggest that dysbiosis triggered by chemotherapy could affect epithelial responses and may play a role in the course of oral mucositis lesions.
Despite the adverse clinical impact of oral mucositis, there is insufficient knowledge on its pathophysiology. Enhanced understanding of events that mediate oral mucosal damage and whether microbial signaling modifies the response to antineoplastics is important for the future development of novel therapies. To our knowledge, our work represents the most comprehensive characterization to date of changes in the human oral environment during chemotherapy in relation to oral mucositis. Our study cohort included patients undergoing treatment based on either 5-FU or doxorubicin. Although both groups developed oral mucositis, its incidence was higher in the 5-FU group with a highly significant correlation seen between the total dose of 5-FU and the severity of lesions. 5-FU is an anti-metabolite drug widely used either alone or as a foundational therapeutic in combination treatment regimens for a range of cancers, including colorectal, breast, and cancers of the aerodigestive tract . A recently developed animal model showed that 5-FU induces atrophy of the oral mucosa with reduced proliferation, increased cell death, upregulation of pro-inflammatory mediators, and compromised barrier integrity . The mechanisms by which 5-FU alone or in combination with adjunct chemotherapeutics exerts its toxic effects on the human oral mucosa, however, remain incompletely understood, with the role of the microbiome in oral mucositis almost entirely unexplored.
Our work demonstrates that the oral microbiome is disrupted by chemotherapy. Although we expected chemotherapy to affect the oral mycobiome due to myelosuppression, we did not observe major changes in the composition of fungal communities. Instead, it was the oral bacteriome that displayed major modifications. These changes were of greater magnitude in salivary communities, but were also detected in mucosal samples. Shifts in the oral bacteriome were strongly correlated with mucositis severity. Bacteriome shifts, however, affected the mouth as a whole and were not specific to the surfaces of lesions as we did not detect any differences in the mucosal bacterial communities present at healthy, erythematous, or ulcerated sites. It should be noted, however, that since swabs were used to sample mucosal lesions, there could have been subsurface microorganisms in lesions not captured by our sampling approach.
An exploration of possible causes of bacteriome shifts found that the use of a multiple dose antibiotic correlated with decreased salivary bacterial diversity but showed only a weak correlation with changes in salivary bacterial community structure. Antibiotic intake mostly affected the oral microbiome at the last visit (V4) and its effect was independent from mucositis-associated dysbiotic shifts. Furthermore, a single-dose prophylactic antibiotic exposure did not correlate with oral microbiome changes. This lack of a major effect of antibiotics on the oral microbiome is not surprising, as previous studies indicate the oral microbiome is resistant to change after systemic antibiotic administration, in contrast to gut communities, which are easily disrupted by antibiotics . Microbiome changes during chemotherapy and mucositis could be neither explained by impaired salivary flow, which in fact was increased during chemotherapy and positively associated with mucositis severity. This result contrasted with previous reports in which patients undergoing repeated chemotherapy had decreased SFR . The discrepancy with our results is probably due to the fact that our cohort was comprised of patients who were either chemotherapy-naïve or in their first cycles of treatment. It is possible that repeated exposure to antineoplastics could cause induction of temporary glandular damage therefore affecting SFR. Similarly, dysbiotic microbiome shifts were not related to decreased oral granulocyte availability, which in fact was positively correlated with more severe mucositis. Increased oral neutrophils during mucositis may be the result of upregulation in chemoattractants such as CXCL2 by the epithelium, as seen in our transcriptome data, in combination with barrier defects through which neutrophils escape tissues. Furthermore, evaluation of a direct antibacterial effect of chemotherapeutic drugs did not support the possibility that chemotherapeutics, in particular 5-FU, were responsible for microbiome shifts associated with mucositis, although 5-FU showed antimicrobial activity against whole salivary communities and against a few of the species tested. Salivary concentrations of 5-FU have been measured at approximately 77 μM 12 min after the start of 5-FU administration, but concentrations rapidly decrease to about 7 μM after 3 h and become undetectable minutes after the end of a 22-h infusion .Based on our results, it is thus possible that for a short period of time during 5-FU administration, the viability of selected oral bacteria is affected by the drug, but salivary 5-FU concentrations rapidly decline allowing oral commensals to recover. In addition, it should be noted that most study subjects received 5-FU within a maximum 5-day window, while most of the microbiome changes associated with mucositis were seen to occur at visits 3 and 4, which were days apart from the last infusion. Based on the salivary drug kinetics described above, it is therefore unlikely that mucositis-associated microbiome shifts were due to a direct antibacterial effect of 5-FU. In support of this, we did not observe that microorganisms whose depletion correlated with high 5-FU doses and more severe mucositis were consistently sensitive to 5-FU, while the growth of F. nucleatum, a microorganism enriched during mucositis, was inhibited by the drug.
The correlation of microbiome shifts and mucositis is in agreement with a previous study that evaluated changes in the oral mucosal bacterial communities in children undergoing chemotherapy, finding greater microbiome modifications in children that developed oral mucositis . This prior study, however, did not define specific species associated with lesion development. We found that shifts in bacterial communities during mucositis consist of depletion of abundant species from genera associated with health and enrichment of Gram-negative species linked to other oral inflammatory conditions . Mucositis-associated microbiome shifts may have resulted from the release of inflammatory products into the oral environment after chemotherapy-induced damage and inflammation. Bacterial species from genera found enriched in mucositis are known to thrive using inflammatory products as nutrients [30, 31], and therefore, the observed mucositis-associated microbiome shifts could have resulted from alterations in oral nutritional resources that promote the growth of these Gram-negative species.
To better understand epithelial cellular events altered by chemotherapy and potential contributors to oral mucosal injury, we performed a global gene expression analysis of changes in the oral epithelium from baseline to V3. Although the specific sites sampled were not affected by mucositis, most subjects presented with lesions in other areas of the mouth at the time of sampling. One upregulated category was that of genes involved in the innate immune response, including here TNF, IL17C, CXCL2, and CCL20, which are cytokines typically induced in epithelial tissues by stimulation with microorganisms, but not part of the response of cell lines after an antineoplastic challenge [32,33,34,35,36]. This suggested that some of the responses seen in vivo were the result of microbial-host cross talk. In vitro challenge of a 3D oral epithelial construct with F. nucleatum, a microorganism associated with more severe mucositis, confirmed the potential for a dysbiotic microbiome to contribute to upregulation of these innate immune response genes during chemotherapy. In contrast, S. salivarius, a microorganism depleted during mucositis, only promoted a slight upregulation or had no effect on innate immunity mediators. The role microbially induced innate immune responses could play in amplifying oral epithelial damage during chemotherapy warrants further investigation.
An expected category of upregulated genes in the oral epithelium during chemotherapy was that of apoptosis-related genes. Among them, the pro-apoptotic gene PMAIP1 (NOXA) is known to mediate antineoplastic-induced cytotoxicity . NOXA is a BH3s death sentinel that acts by directly inducing stepwise, bimodal activation of BAX–BAK, which mediate permeabilization of mitochondria, followed by the release of cytochrome c and caspase activation . NOXA also binds to the anti-apoptotic protein MCL-1 and may induce apoptosis by preventing its activity [24, 37, 38]. PMAIP1 is upregulated in gastric epithelial cells upon Helicobacter pylori infection, which suggests a link between microbial sensing and PMAIP1 transcription . Unexpectedly, we also observed that in vitro challenge of oral epithelial cells with F. nucleatum upregulated expression of PMAIP1 and consistent with this, F. nucleatum-challenged tissues displayed cell death markers. This suggested that the upregulation of signals leading to cellular death could also occur in response to oral dysbiosis, which may serve to aggravate epithelial injury.
Another factor identified as associated with mucositis was wearing a removable prosthesis. Contact of a denture with oral tissues may constitute a mechanical insult that aggravates antineoplastic-driven epithelial injury. Dentures are also known to promote the overgrowth of oral commensals, which in turn could modify mucosal damage . Additional factors associated with a linear increase in mucositis severity were smoking and steroid use during the cycle. Both smoking and steroids are likely to modify innate immune responses of mucosal tissues that are important for controlling dysbiosis [41, 42], and may play a role in mucositis pathophysiology.
In summary, our work shows that the severity of oral mucositis lesions correlates with the degree of oral microbiome dysbiosis, with enrichment of potential Gram-negative pathobionts and depletion of health-associated symbionts in severe mucositis. These microbiome alterations resemble the dysbiotic shifts associated with other oral diseases in which changes in the environment and nutrient resources, due to inflammation, induce the outgrowth of low abundant species, which in turn dysregulate host defenses leading to damage . Indeed, F. nucleatum, a potential pathobiont enriched as mucositis became more severe, showed an ability to induce pro-inflammatory and apoptotic responses in a 3D oral epithelium model, while S. salivarius, a potential symbiont depleted in severe mucositis, was tolerated. While antineoplastics represent the primary initiators of mucositis, our observations also suggest that antineoplastic-triggered inflammation may induce microbiome disruptions and that such dysbiotic microbiome could play a role in the clinical course of lesions aggravating epithelial injury. Further studies are warranted to elucidate whether interfering with dysbiotic events could reduce oral mucositis incidence and/or severity.
This observational prospective study with two cohorts was approved by the Institutional Review Board at UConn Health (IRB number IE-11-037 J-2) and conformed with STROBE guidelines. Written informed consent was received from all participants. Additional file 1: Figure S1 shows the study design and inclusion and exclusion criteria. Subjects planned to or already receiving chemotherapeutic treatment for a solid tumor either on a 5-fluorouracil (5-FU)- or a doxorubicin-based regimen were recruited at the Neag Comprehensive Cancer Center of UConn Health and at the Helen and Harry Gray Cancer Center at Hartford Hospital. Every effort was made to recruit chemotherapy-naïve subjects. From the 49 enrollments in the cancer cohort, 32 were chemotherapy-naïve while 17 had already completed at least one previous cycle. Thirty non-cancer controls were enrolled from the general local population through the UConn Health Dental Clinical Research Center. Cancer subjects were seen at four visits which included a baseline visit completed prior to subjects commencing chemotherapy, and three more visits (V2–V4) within 14 ± 2 days after commencing treatment. Control subjects were seen at two visits, including a baseline and a visit within 14 ± 2 days (V4).
Collection of demographic and medical data and oral evaluation
Medical information was obtained via questionnaires and from medical charts. All subjects received at baseline an oral evaluation including assessment of the presence of periodontitis via the Community Periodontal Index of Treatment Needs (CPITN) , presence and type of prosthetic restorations, number of teeth, and presence of visible cavitated carious lesions. In addition, at all visits subjects were evaluated for the presence and severity of mucositis per the World Health Organization (WHO) scale, which assesses oral soreness, erythema, ulceration, and ability to eat using a 0–4 scale, giving a single score for the entire oral cavity . A score of 0 indicates no signs or symptoms; 1 indicates erythema and soreness; 2 indicates ulcers, able to eat solids; 3 indicates ulcers, requires liquid diet (due to mucositis); and 4 indicates ulcers, alimentation not possible (due to mucositis). A second mucositis assessment method, the Oral Mucositis Assessment Scale (OMAS) was also used . The OMAS provides an objective assessment of oral mucositis based on an evaluation of the appearance and the extent of redness and ulceration in various areas of the mouth. Nine oral mucosal sites were evaluated for presence and severity of mucositis ulcerative lesions assigning to each site a score from 0 to 3 depending on lesion size, and for severity of erythema, assigning each site a score from 0 to 2. These scores were then aggregated per site and therefore each site could have a score ranging from 0 to 5. Scores from all nine sites were then aggregated for a total mouth score ranging from 0 to 45. Aggregated total mouth scores were used for analysis. At all visits, a subjective evaluation of oral dryness was also conducted , and an assessment of the amount of plaque on teeth, as determined via the Plaque Index of Silness and Löe .
Salivary flow rate determination and collection of samples for microbiome evaluation
At each visit, unstimulated saliva and mucosal swabs were collected for microbiome evaluation. Participants were instructed to avoid eating or drinking anything other than water for 1 h prior to each study visit. Unstimulated saliva was collected for 5 min by having subjects lean forward over a sterile funnel attached to a sterile vial placed on ice. Saliva was weighted to determine salivary flow rate, and then aliquoted and centrifuged at 6000×g for 10 min. Supernatants were removed and pellets stored at − 80 °C. Mucosal swab samples were collected in subjects presenting mucosal health by swabbing the entire right and left buccal mucosa for 10 s each with a single CatchAll™ swab. In subjects presenting with mucosal lesions, separate swabs were collected from healthy mucosal areas, erythematous tissue, and areas with ulceration. Swabs were immediately swirled in a tube containing 500 μL of TE buffer (20 mM Tris pH 7.0, 2 mM EDTA) pressing against the tube walls to release the swab contents. Samples were then stored at − 80 °C.
Collection of samples for oral epithelium transcriptome analysis
Oral mucosal cytology samples were collected using a Cytobrush™, which was pressed and rotated against a healthy mucosal area for 60 s. If the procedure lead to mucosal bleeding, the sample was discarded to avoid contamination by cells beyond the epithelium. After collection, the brush was swirled in 1 mL of RNAProtect (Qiagen), pressing against the tube walls to release the brush contents. Cells were pelleted at 8000×g for 20 min, snap-frozen, and stored in liquid nitrogen.
Evaluation of blood and oral neutrophil counts
Total counts and percentages of neutrophils were obtained from complete blood counts performed by the hospital laboratory using an automated Beckman Coulter analyzer. Oral cells were collected at all study visits by asking subjects to rinse with 10 mL of a bicarbonate solution for 30 s. Oral rinse samples were centrifuged at 1258×g for 15 min and pellets were resuspended in Hanks’ Balanced Salt solution supplemented with 2 μg/ml acridine orange. Granulocytes, most of which are neutrophils, characteristically stained as multilobulated cells with this nuclear stain, were enumerated in a hemacytometer under a fluorescence microscope .
DNA extraction for microbiome evaluation
DNA was extracted separately for bacterial and fungal microbiome evaluation. For bacteria, we followed a previously described procedure using lysozyme and proteinase K treatment and the DNeasy Blood and Tissue kit (Qiagen) . For fungi, the DNA extraction protocol involved bead beating with a matrix containing Lysing Matrix B (MP Biomedicals) and very high density 0.5 mm yttrium-stabilized zirconium oxide (95% ZrO2 + 5% Y2O3) grinding media (YSZ) (Glen Mills Inc, Clifton, NJ), followed by extraction using the FastDNA SPIN KIT (MP Biomedicals), as previously described .
Amplification and sequencing of 16S rRNA gene and ITS-1 DNA
Amplicon libraries of the 16S rRNA gene V1–V2 hypervariable regions were generated in triplicate using fusion primers which included universal primers 8F AGAGTTTGATCMTGGCTCAG or 361R CYIACTGCTGCCTCCCGTAG . PCR conditions have been described previously . For mycobiome characterization, fusion primers containing fungal-specific ITS1F forward primer (CTTGGTCATTTAGAGGAAGTAA) or ITS2 reverse primer (GCTGCGTTCTTCATCGATGC) were used to amplify ITS-1 sequences in triplicate, as previously described . Combined libraries were further purified and sequenced using 454 Titanium chemistry and the 454-GS-FLX sequencing platform (454 Life Sciences, Branford, CT). Since mucosal samples did not consistently yield ITS-1 PCR products, probably due to low fungal load, only saliva samples were evaluated via ITS-1 sequencing.
Processing of sequences and taxonomic classification
16S rRNA gene reads were processed in mothur . Primers and barcodes were trimmed followed by removal of sequences shorter than 200 bp, with homopolymers greater than eight nucleotides or ambiguous base calls. Sequences were then filtered using a 50 bp sliding window approach and an average quality score threshold of 35 . Chimeric sequences were removed with UChime , in mothur. Sequences were then classified to species level by using the classify.seqs command and the Human Oral Microbiome Database (HOMD) V14.5 as reference. Parameters used were: method = knn, search = blast, gapopen = − 5, gapextend = − 5, match = 4, mismatch = − 5, numwanted = 1, following recommendations by Al-Hebshi et al. . This taxonomy assignment algorithm was validated by classifying HOMD reference sequences trimmed to include only the V1–V2 region against the HOMD full-length reference sequence database. With a few exceptions, all V1–V2 short sequences were correctly classified in the validation test. However, the following species could not be discriminated from each other: Lactobacillus casei and Lactobacillus rhamnosus; Veillonella parvula and Veillonella dispar; Streptococcus mitis, Streptococcus pneumoniae, and Streptococcus sp. HOT423; and Neisseria flavescens and Neisseria subflava. Counts for species that could not be correctly identified were aggregated.
ITS-1 reads were also processed in mothur . Sequences were trimmed and quality filtered and chimeras removed as described for 16S rRNA sequences. Sequences were then classified to genus level by using the classify.seqs command and a modified version of the Findley et al. database  as a reference, employing the recommended parameters (https://www.mothur.org/wiki/Findley_ITS_Database). Prior to this, a curation of the Findley et al. database was conducted to join synonym taxa under one preferred name . To accomplish this, reference sequences from synonym taxa were compared via BLAST against NCBI’s nucleotide type strain database and against the Fungal Metagenomics Project database. After confirming their genus identity, Gueomyces sequences were included under Trichosporon; Lewia under Alternaria; Valsa under Cytospora; Coprinellus and Coprinopsis under Coprinus; Erythrobasidium under Rhodotorula; Cochliobolus under Curvularia; Filobasidium, Cystofilobasidium and Dioszegia under Cryptococcus; and Emericella under Aspergillus. In addition, sequences of Cladosporium, Toxicocladosporium, Aureobasidium, Kabatiella, Scleroconidioma, and Candida were added as they were underrepresented in the database. Malassezia and Candida ITS-1 sequences were further classified to species level, using curated and aligned reference libraries and the software package pplacer , as previously described . A total of 691 saliva- and mucosa-derived libraries were sequenced. All sequences are available at the Short Reads Archive (Accession number PRJNA399163).
Alpha and beta-diversity estimates for 16S rRNA gene and ITS-1 amplicon data
16S rRNA gene libraries were subsampled at 3074 reads and ITS-1 reads at 1338 reads. Rarefaction curves were generated via the R package vegan and can be found in Additional file 1: Figure S5. Most curves were starting to become asymptotic at the thresholds used for subsampling. Community diversity was evaluated via the non-parametric Shannon Index and community structure via the Yue-Clayton Theta Index (ThetaYC) as calculated in mothur. These metrics were constructed based on species-level taxonomic units for 16S rRNA gene reads and genus-level for ITS-1 (with the exception of Malassezia and Candida that were speciated). Mucosal communities were analyzed in two different manners, including a separate analysis of samples according to the clinical status of the sampled site (healthy, erythematous, or ulcerated) and an analysis in which data from different sites within a subject were combined and analyzed as one mucosal unit.
Evaluation of the antimicrobial activity of chemotherapeutics
To evaluate the effect of antineoplastics on the whole cultivable salivary microbiota, unstimulated saliva from eight subjects was collected in sterile tubes. Aliquots were incubated with 5-FU (7.7 μM), 5-FU (7.7 μM) + docetaxel (1.2 μM), or a vehicle control, at 37 °C for 2 h. Samples were serially diluted and plated on solid media containing tryptic soy agar (20 g/L), brain heart infusion agar (26 g/L), yeast extract (10 g/L), sheep blood (5% vol/vol), hemin (5 μg/mL), menadione (0.3 μg/mL), and N-acetyl muramic acid (10 ug/mL) for enumeration of salivary bacteria . Sets of duplicate plates were incubated either anaerobically or aerobically at 37 °C for 5 days.
To evaluate the effect of 5-FU on pure bacterial cultures, microorganisms were grown in appropriate medium and atmosphere under late logarithmic phase and this culture was diluted in fresh medium to an OD(600 nm) = 0.4. Normalized culture aliquots were incubated with different concentrations of 5-FU or with a vehicle control for 2 h, at 37 °C, under an appropriate atmosphere, after which cultures were serially diluted and plated for colony forming unit determination. The effect of 5-FU on the growth of S. salivarius ATCC 9222 and F. nucleatum ATCC 49256 was determined by adding 5-FU or vehicle control to starter cultures diluted to an OD = 0.1. Growth was followed until early stationary phase. Inhibition was determined by comparing the area under the curve of test conditions compared to control vehicle. All experiments were repeated three times. Concentrations tested of 5-FU and docetaxel in the above experiments were based on their bioavailability in plasma and oral fluids [59, 60].
Changes in global oral epithelial gene expression during chemotherapy
RNA was isolated from cytologic smears using the miRNeasy Mini Kit (Qiagen). RNA quantity and quality were measured using the Experion High Sensitivity RNA Kit. As previously observed by others , RNA quantity and quality from oral epithelium varies with some samples showing very low yields and high degradation. Therefore, we employed the whole genome cDNA-mediated Annealing Selection Extension and Ligation assay (wgDASL) in combination with the HumanHT-12 v4 Expression Bead Chip (Illumina), to evaluate gene expression. This assay performs well with partially degraded RNA as it does not depend on an intact poly-A tail for cDNA synthesis . Only subjects with matched baseline and V3 samples that yielded at least 500 ng of RNA were included in the analysis. The wgDASL assay was performed with 250 ng of RNA, including two technical replicates. Initial data analysis to eliminate samples that did not perform well in the assay was carried out using GenomeStudio (Illumina). Data were first normalized using the quantile method, followed by filtering out probes without a detection p value < 0.05 in at least one sample. Outlier samples for which technical replicates did not have a correlation of at least r = 0.95 were removed from the dataset. Differential gene expression analysis in 14 subject-paired samples from baseline to V3 was then performed using Significance Analysis of Microarrays (SAM) . Gene ontology enrichment analysis was performed using DAVID , and redundant GO terms summarized and visualized using REVIGO .
Effect of commensals on a 3D oral mucosa analog tissue
To evaluate the response of the oral mucosa to microbial commensals, 3D multi-layered epithelial cell constructs derived from primary keratinocytes (EpiOral™, MatTek Corporation) were challenged with 108 bacterial cells per cm2. Microorganisms were placed apically and incubated with tissues for 24 h. Membranes containing tissues were harvested and used for histologic sectioning or RNA extraction. Cell death was evaluated in tissue sections via the DeadEnd™ Fluorometric TUNEL System (Promega). RNA was extracted via Qiagen’s RNeasy mini kit, reverse transcribed using the High Capacity RNA-to-cDNA Kit (Applied Biosystems), and evaluated via real-time PCR using TaqMan Assays (Applied Biosystems). The 18S rRNA gene was used as endogenous control.
Baseline data were compared between groups by either chi-square or Fisher’s exact test, the latter applied when more than 20% of table cell counts had expected frequencies < 5. Continuous data were evaluated via independent sample t tests or Mann-Whitney U tests, according to data distribution. The distribution of continuous data was tested for normality using measures of Skewness and Kurtosis and the Shapiro-Wilk test for normality in SPSS. Differences in mucositis incidence across regimens were evaluated via chi-square. Longitudinal data were evaluated in subject-matched samples. Of interest was how changes in one variable were associated with changes in another variable or with one time-point measures. Specifically, we were interested in the longitudinal analysis of changes in mucositis scores and the relationship of these changes with other variables. We therefore modelled the change in each longitudinal variable using linear or quadratic orthogonal polynomial contrasts. For this, each data point was multiplied by a standard coefficient (for a linear or quadratic 4 level contrast) and data aggregated by subject. Correlations between linear or quadratic longitudinal change and other variables were determined via Spearman Rank tests. Differences in clinical variables and microbiome diversity from baseline to each visit were calculated using paired-Wilcoxon Rank tests. The change in alpha-diversity (measured with the non-parametric Shannon Index) was also evaluated by subtracting values at each visit from the respective baseline values. Correlations between the difference in non-parametric Shannon Index values and other variables were evaluated via Spearman Rank tests. To quantify changes in community structure during chemotherapy, we calculated the ThetaYC distance of each sample to its corresponding baseline and compared distances to control subjects via Mann-Whitney Rank tests. ThetaYC distances were also evaluated for their correlation with demographic and clinical variables via Spearman Rank tests. Linear regression analysis was used to evaluate the contribution of different predictors to alpha- and beta-diversity microbiome change. Changes in salivary community structure during chemotherapy and correlation of these changes with clinical variables and microbiome diversity were also visualized in Principal Coordinates Analysis (PCoA) plots based on the ThetaYC distances. Significant separation between data clouds was tested via AMOVA and the correlations of principal components with metadata were evaluated via Spearman Rank tests. Differences in relative abundances of individual microbial taxa from baseline to each visit were evaluated via Wilcoxon signed rank tests. The significance threshold for all statistical tests mentioned above was adjusted using the Benjamini Hochberg false discovery rate method. Evaluation of the longitudinal covariation in microbial relative abundances and clinical parameters measured at the four visits of the study was performed using the R package MixOmics . For this analysis, we used the multilevel function and compared microbial and clinical variables using sparse Partial Least Squares (sPLS). This multilevel multivariate approach takes into account the repeated measures within subjects and by using sPLS reveals covariation patterns in all variables in an unsupervised manner. Differences in colony forming units (CFUs) recovered after treatment of saliva or pure bacterial cultures with antineoplastics were determined via paired Wilcoxon signed rank tests. Differences in growth of bacteria in the presence and absence of 5-FU were determined by comparing the area under the growth curve for different conditions via t tests.
Community Periodontal Index of Treatment Needs
Human Oral Microbiome Database
Oral Mucositis Assessment Scale
Salivary flow rate
Sparse Partial Least Squares analysis
- ThetaYC :
Yue-Clayton Theta Index
World Health Organization
Vera-Llonch M, Oster G, Ford CM, Lu J, Sonis S. Oral mucositis and outcomes of allogeneic hematopoietic stem-cell transplantation in patients with hematologic malignancies. Support Care Cancer. 2007;15(5):491–6.
Keefe DM, Schubert MM, Elting LS, Sonis ST, Epstein JB, Raber-Durlacher JE, et al. Updated clinical practice guidelines for the prevention and treatment of mucositis. Cancer. 2007;109(5):820–31.
Murphy BA, Beaumont JL, Isitt J, Garden AS, Gwede CK, Trotti AM, et al. Mucositis-related morbidity and resource utilization in head and neck cancer patients receiving radiation therapy with or without chemotherapy. J Pain Symptom Manage. 2009;38(4):522–32.
Rosenthal DI. Consequences of mucositis-induced treatment breaks and dose reductions on head and neck cancer treatment outcomes. J Support Oncol. 2007;5(9 Suppl 4):23–31.
Peterson DE, Srivastava R, Lalla RV. Oral mucosal injury in oncology patients: perspectives on maturation of a field. Oral Dis. 2015;21(2):133–41.
Dasari S, Tchounwou PB. Cisplatin in cancer therapy: molecular mechanisms of action. Eur J Pharmacol. 2014;740:364–78.
Wilson PM, Danenberg PV, Johnston PG, Lenz HJ, Ladner RD. Standing the test of time: targeting thymidylate biosynthesis in cancer therapy. Nat Rev Clin Oncol. 2014;11(5):282–98.
Sonis ST. The Chicken or the Egg? Changes in Oral Microbiota as Cause or Consequence of Mucositis During Radiation Therapy. EBioMedicine. 2017;18:7–8.
Vasconcelos RM, Sanfilippo N, Paster BJ, Kerr AR, Li Y, Ramalho L, et al. Host-microbiome cross-talk in oral mucositis. J Dent Res. 2016;95(7):725–33.
Pedroso SH, Vieira AT, Bastos RW, Oliveira JS, Cartelle CT, Arantes RM, et al. Evaluation of mucositis induced by irinotecan after microbial colonization in germ-free mice. Microbiol. 2015;161(10):1950–60.
Frank M, Hennenberg EM, Eyking A, Runzi M, Gerken G, Scott P, et al. TLR signaling modulates side effects of anticancer therapy in the small intestine. J Immunol. 2015;194(4):1983–95.
Napenas JJ, Brennan MT, Bahrani-Mougeot FK, Fox PC, Lockhart PB. Relationship between mucositis and changes in oral microflora during cancer chemotherapy. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2007;103(1):48–59.
de Mendonca RM, de Araujo M, Levy CE, Morari J, Silva RA, Yunes JA, et al. Prospective evaluation of HSV, Candida spp., and oral bacteria on the severity of oral mucositis in pediatric acute lymphoblastic leukemia. Support Care Cancer. 2012;20(5):1101–7.
Perales-Puchalt A, Perez-Sanz J, Payne KK, Svoronos N, Allegrezza MJ, Chaurio RA, et al. Frontline Science: Microbiota reconstitution restores intestinal integrity after cisplatin therapy. J Leukoc Biol. 2018;103(5):799–805.
Jensen SB, Mouridsen HT, Reibel J, Brünner N, Nauntofte B. Adjuvant chemotherapy in breast cancer patients induces temporary salivary gland hypofunction. Oral Oncol. 2008;44(2):162–73.
Marsh PD, Do T, Beighton D, Devine DA. Influence of saliva on the oral microbiota. Periodontol 2000. 2016;70(1):80–92.
Maier L, Pruteanu M, Kuhn M, Zeller G, Telzerow A, Anderson EE, et al. Extensive impact of non-antibiotic drugs on human gut bacteria. Nature. 2018;555(7698):623–8.
Vanlancker E, Vanhoecke B, Smet R, Props R, Van de Wiele T. 5-Fluorouracil sensitivity varies among oral micro-organisms. J Med Microbiol. 2016;65(8):775–83.
Diaz PI, Hoare A, Hong BY. Subgingival microbiome shifts and community dynamics in periodontal diseases. J Calif Dent Assoc. 2016;44(7):421–35.
Lamont RJ, Hajishengallis G. Polymicrobial synergy and dysbiosis in inflammatory disease. Trends Mol Med. 2015;21(3):172–83.
Rijkschroeff P, Loos BG, Nicu EA. Impaired polymorphonuclear neutrophils in the oral cavity of edentulous individuals. Eur J Oral Sci. 2017;125(5):371–8.
Diaz PI, Chalmers NI, Rickard AH, Kong C, Milburn CL, Palmer RJ Jr, et al. Molecular characterization of subject-specific oral microflora during initial colonization of enamel. Appl Environ Microbiol. 2006;72(4):2837–48.
Wang X, Lin Y. Tumor necrosis factor and cancer, buddies or foes? Acta Pharmacol Sin. 2008;29(11):1275–88.
Chen HC, Kanai M, Inoue-Yamauchi A, Tu HC, Huang Y, Ren D, et al. An interconnected hierarchical model of cell death regulation by the BCL-2 family. Nat Cell Biol. 2015;17(10):1270–81.
Sheridan C, Brumatti G, Elgendy M, Brunet M, Martin SJ. An ERK-dependent pathway to Noxa expression regulates apoptosis by platinum-based chemotherapeutic drugs. Oncogene. 2010;29(49):6428–41.
Bertolini M, Sobue T, Thompson A, Dongari-Bagtzoglou A. Chemotherapy induces oral mucositis in mice without additional noxious stimuli. Transl Oncol. 2017;10(4):612–20.
Zaura E, Brandt BW, Teixeira de Mattos MJ, Buijs MJ, Caspers MP, Rashid MU, et al. Same exposure but two radically different responses to antibiotics: resilience of the salivary microbiome versus long-term microbial shifts in feces. mBio. 2015;6(6):e01693–15.
Joulia JM, Pinguet F, Ychou M, Duffour J, Astre C, Bressolle F. Plasma and salivary pharmacokinetics of 5-fluorouracil (5-FU) in patients with metastatic colorectal cancer receiving 5-FU bolus plus continuous infusion with high-dose folinic acid. Eur J Cancer. 1999;35(2):296–301.
Ye Y, Carlsson G, Agholme MB, Wilson JA, Roos A, Henriques-Normark B, et al. Oral bacterial community dynamics in paediatric patients with malignancies in relation to chemotherapy-related oral mucositis: a prospective study. Clin Microbiol Infect. 2013;19(12):E559–67.
ter Steeg PF, Van der Hoeven JS, de Jong MH, van Munster PJ, Jansen MJ. Enrichment of subgingival microflora on human serum leading to accumulation of Bacteroides species, Peptostreptococci and Fusobacteria. Antonie van Leeuwenhoek. 1987;53(4):261–72.
ter Steeg PF, van der Hoeven JS, de Jong MH, van Munster PJJ, Jansen MJH. Modelling the gingival pocket by enrichment of subgingival microflora in human serum in chemostats. Microb Ecol Health Dis. 1988;1:73–84.
Ramirez-Carrozzi V, Sambandam A, Luis E, Lin Z, Jeet S, Lesch J, et al. IL-17C regulates the innate immune function of epithelial cells in an autocrine manner. Nat Immunol. 2011;12(12):1159–66.
Dickinson BC, Moffatt CE, Hagerty D, Whitmore SE, Brown TA, Graves DT, et al. Interaction of oral bacteria with gingival epithelial cell multilayers. Mol Oral Microbiol. 2011;26(3):210–20.
Conti HR, Whibley N, Coleman BM, Garg AV, Jaycox JR, Gaffen SL. Signaling through IL-17C/IL-17RE is dispensable for immunity to systemic, oral and cutaneous candidiasis. PloS One. 2015;10(4):e0122807.
Maxwell PJ, Longley DB, Latif T, Boyer J, Allen W, Lynch M, et al. Identification of 5-fluorouracil-inducible target genes using cDNA microarray profiling. Cancer Res. 2003;63(15):4602–6.
Cho YE, Singh TS, Lee HC, Moon PG, Lee JE, Lee MH, et al. In-depth identification of pathways related to cisplatin-induced hepatotoxicity through an integrative method based on an informatics-assisted label-free protein quantitation and microarray gene expression approach. Mol Cel Proteomics. 2012;11(1):M111 010884.
Czabotar PE, Lee EF, van Delft MF, Day CL, Smith BJ, Huang DC, et al. Structural insights into the degradation of Mcl-1 induced by BH3 domains. PNAS. 2007;104(15):6217–22.
Nakajima W, Hicks MA, Tanaka N, Krystal GW, Harada H. Noxa determines localization and stability of MCL-1 and consequently ABT-737 sensitivity in small cell lung cancer. Cell Death Dis. 2014;5:e1052.
Rath S, Das L, Kokate SB, Pratheek BM, Chattopadhyay S, Goswami C, et al. Regulation of Noxa-mediated apoptosis in Helicobacter pylori-infected gastric epithelial cells. FASEB J. 2015;29(3):796–806.
Narhi TO, Kurki N, Ainamo A. Saliva, salivary micro-organisms, and oral health in the home-dwelling old elderly--a five-year longitudinal study. J Dent Res. 1999;78(10):1640–6.
Qian YJ, Wang X, Gao YF, Duan N, Huang XF, Sun FF, et al. Cigarette smoke modulates NOD1 signal pathway and human beta defensins expression in human oral mucosa. Cell Physiol Biochem. 2015;36(2):457–73.
Conti HR, Shen F, Nayyar N, Stocum E, Sun JN, Lindemann MJ, et al. Th17 cells and IL-17 receptor signaling are essential for mucosal host defense against oral candidiasis. J Exp Med. 2009;206(2):299–311.
Ainamo J, Barmes D, Beagrie G, Cutress T, Martin J, Sardo-Infirri J. Development of the World Health Organization (WHO) community periodontal index of treatment needs (CPITN). Int Dent J. 1982;32(3):281–91.
Peterson DE, Boers-Doets CB, Bensadoun RJ, Herrstedt J, Committee EG. Management of oral and gastrointestinal mucosal injury: ESMO Clinical Practice Guidelines for diagnosis, treatment, and follow-up. Ann Oncol. 2015;26(Suppl 5):v139–51.
Sonis ST, Eilers JP, Epstein JB, LeVeque FG, Liggett WH Jr, Mulagha MT, et al. Validation of a new scoring system for the assessment of clinical trial research of oral mucositis induced by radiation or chemotherapy. Mucositis Study Group. Cancer. 1999;85(10):2103–13.
Fall-Dickson JM, Mitchell SA, Marden S, Ramsay ES, Guadagnini JP, Wu T, et al. Oral symptom intensity, health-related quality of life, and correlative salivary cytokines in adult survivors of hematopoietic stem cell transplantation with oral chronic graft-versus-host disease. Biol Blood Marrow Transplant. 2010;2010:4.
Silness J, Loe H. Periodontal disease in pregnancy. II. Correlation between oral hygiene and periodontal condition. Acta Odontol Scand. 1964;22:121–35.
Cheretakis C, Dror Y, Glogauer M, Leung R, Sun CX. A noninvasive oral rinse assay to monitor engraftment, neutrophil tissue delivery and susceptibility to infection following HSCT in pediatric patients. Bone Marrow Transplant. 2005;36(3):227–32.
Diaz PI, Dupuy AK, Abusleme L, Reese B, Obergfell C, Choquette L, et al. Using high throughput sequencing to explore the biodiversity in oral bacterial communities. Mol Oral Microbiol. 2012;27(3):182–201.
Dupuy AK, David MS, Li L, Heider TN, Peterson JD, Montano EA, et al. Redefining the human oral mycobiome with improved practices in amplicon-based taxonomy: discovery of Malassezia as a prominent commensal. PloS One. 2014;9(3):e90899.
Sundquist A, Bigdeli S, Jalili R, Druzin ML, Waller S, Pullen KM, et al. Bacterial flora-typing with targeted, chip-based pyrosequencing. BMC Microbiol. 2007;7:108.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Schloss PD, Gevers D, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PloS One. 2011;6(12):e27310.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Al-Hebshi NN, Nasher AT, Idris AM, Chen T. Robust species taxonomy assignment algorithm for 16S rRNA NGS reads: application to oral carcinoma samples. J Oral Microbiol. 2015;7:28934.
Findley K, Oh J, Yang J, Conlan S, Deming C, Meyer JA, et al. Topographic diversity of fungal and bacterial communities in human skin. Nature. 2013;498(7454):367–70.
Diaz PI, Hong BY, Dupuy AK, Strausbaugh LD. Mining the oral mycobiome: methods, components, and meaning. Virulence. 2017;8(3):313–23.
Matsen FA, Kodner RB, Armbrust EV. pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree. BMC Bioinformatics. 2010;11:538.
Feres M, Haffajee AD, Allard K, Som S, Goodson JM, Socransky SS. Antibiotic resistance of subgingival species during and after antibiotic therapy. J Clin Periodontol. 2002;29(8):724–35.
Rosing H, Lustig V, van Warmerdam LJ, Huizing MT, ten Bokkel Huinink WW, Schellens JH, et al. Pharmacokinetics and metabolism of docetaxel administered as a 1-h intravenous infusion. Cancer Chemother Pharmacol. 2000;45(3):213–8.
Spira A, Beane J, Schembri F, Liu G, Ding C, Gilman S, et al. Noninvasive method for obtaining RNA from buccal mucosa epithelial cells for gene expression profiling. Biotechniques. 2004;36(3):484–7.
Waddell N, Cocciardi S, Johnson J, Healey S, Marsh A, Riley J, et al. Gene expression profiling of formalin-fixed, paraffin-embedded familial breast tumours using the whole genome-DASL assay. J Pathol. 2010;221(4):452–61.
Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. PNAS. 2001;98(9):5116–21.
Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature Prot. 2009;4(1):44–57.
Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PloS One. 2011;6(7):e21800.
Rohart F, Gautier B, Singh A, Le Cao KA. mixOmics: An R package for ‘omics feature selection and multiple data integration. PLoS Comput Biol. 2017;13(11):e1005752.
The authors would like to acknowledge Dr. Julia Oh (Jackson Laboratory for Genomic Medicine) for kindly providing us with the scripts and database to classify Malassezia species. The authors would also like to acknowledge Dr. Niki M. Moutsopoulos (NIDCR, NIH) for reading the manuscript and providing helpful comments.
This work was supported by grant R01DE021578 from The National Institute of Dental and Craniofacial Research (NIDCR), National Institutes of Health (NIH). We also acknowledge the support of pilot grant M01RR006192 from UConn Health General Clinical Research Center.
Availability of data and materials
All 16S rRNA gene and ITS-1 sequences are available at the Short Reads Archive (Accession number PRJNA399163).
Ethics approval and consent to participate
This study was approved by the Institutional Review Board at UConn Health (IRB number IE-11-037 J-2). Written informed consent was received from all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Study design. Figure S2. Effect of chemotherapy on the oral microbiome community structure. Figure S3. Taxa that changed in relative abundance at each study visit compared to baseline levels. Figure S4. Effect of antineoplastics on the oral microbiome. Figure S5. Rarefaction curves for all bacterial and fungal amplicon libraries. (PDF 9205 kb)
Table S1. Demographic and baseline clinical characteristics of study participants. Table S2. Descriptive statistics specific to the cancer cohort. Table S3. Correlations between change in oral mucositis severity (OMAS) or change in salivary bacterial diversity with demographic and clinical characteristics. Table S5. Sequencing summary. Table S6. Within subject differences in alpha-diversity (non-parametric Shannon index) between mucosal sites with different clinical appearance. Table S7. Main clinical characteristics of subjects included in the mucosal gene expression analyses. (PDF 221 kb)
Table S4. Identity and proportions of taxa with significantly different relative abundance at any study visit compared to baseline. (XLSX 15 kb)
Table S8. SAM analysis of differential gene expression between baseline and V3 in 14 subjects, showing 158 upregulated and 2 downregulated genes at V3 (> 1.5-fold). (XLSX 24 kb)
Table S9. Summary of upregulated gene ontology (GO) term categories. (XLSX 13 kb)