Next Article in Journal
IL-10 Rescues CLL Survival through Repolarization of Inflammatory Nurse-like Cells
Next Article in Special Issue
Neoadjuvant Treatment Lowers the Risk of Mesopancreatic Fat Infiltration and Local Recurrence in Patients with Pancreatic Cancer
Previous Article in Journal
Systemic Triple Therapy in Metastatic Hormone-Sensitive Prostate Cancer (mHSPC): Ready for Prime Time or Still to Be Explored?
Previous Article in Special Issue
Cobaltabis(dicarbollide) ([o-COSAN]) as Multifunctional Chemotherapeutics: A Prospective Application in Boron Neutron Capture Therapy (BNCT) for Glioblastoma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Changes to the Oral Microbiome Based on 16S rRNA Sequencing among Children Treated for Cancer

by
Patrycja Proc
1,*,
Joanna Szczepańska
1,
Beata Zarzycka
2,
Małgorzata Szybka
2,
Maciej Borowiec
3,
Tomasz Płoszaj
3,
Wojciech Fendler
4,5,
Jędrzej Chrzanowski
5,
Małgorzata Zubowska
6,
Małgorzata Stolarska
7 and
Wojciech Młynarski
7,*
1
Department of Pediatric Dentistry, Medical University of Lodz, Pomorska 251, 92-213 Lodz, Poland
2
Department of Microbiology and Medical Laboratory Immunology, Medical University of Lodz, Pomorska 251, 92-213 Lodz, Poland
3
Department of Clinical Genetics Medical, University of Lodz, Pomorska 251, 92-213 Lodz, Poland
4
Department of Radiation Oncology, Dana-Farber Cancer Institute, Boston, MA 02215, USA
5
Department of Biostatistics and Translational Medicine, Medical University of Lodz, Mazowiecka 15, 92-215 Lodz, Poland
6
Departments of Pediatrics, Diabetology, Endocrinology and Nephrology, Medical University of Lodz, Sporna 36/50, 91-738 Lodz, Poland
7
Department of Pediatrics, Oncology & Hematology, Medical University of Lodz, Sporna 36/50, 91-738 Lodz, Poland
*
Authors to whom correspondence should be addressed.
Cancers 2022, 14(1), 7; https://doi.org/10.3390/cancers14010007
Submission received: 19 November 2021 / Accepted: 13 December 2021 / Published: 21 December 2021
(This article belongs to the Special Issue Feature Paper from Journal Reviewers)

Abstract

:

Simple Summary

Childhood cancer survivors suffer from many oral complications during and after primary therapy. Our study focuses on changes in the oral microbiome of cancer survivors. Using 16S rRNA sequencing, we observed global and distinct changes in oral microbiome associated with a patient’s age and therapy duration, but not antibiotic therapy or cancer type. Observed changes in the oral microbiome could differentiate patients at higher risk of long-term oral complications.

Abstract

A child’s mouth is the gateway to many species of bacteria. Changes in the oral microbiome may affect the health of the entire body. The aim of the study was to evaluate the changes in the oral microbiome of childhood cancer survivors. Saliva samples before and after anti-cancer treatment were collected from 20 patients aged 6–18 years, diagnosed de novo with cancer in 2018–2019 (7 girls and 13 boys, 7.5–19 years old at the second time point). Bacterial DNA was extracted, and the microbial community profiles were assessed by 16S rRNA sequencing. The relative abundances of Cellulosilyticum and Tannerella genera were found to significantly change throughout therapy (p = 0.043 and p = 0.036, respectively). However, no differences in the alpha-diversity were observed (p = 0.817). The unsupervised classification revealed two clusters of patients: the first with significant changes in Campylobacter and Fusobacterium abundance, and the other with change in Neisseria. These two groups of patients differed in median age (10.25 vs. 16.16 years; p = 0.004) and the length of anti-cancer therapy (19 vs. 4 months; p = 0.003), but not cancer type or antibiotic treatment.

Graphical Abstract

1. Introduction

The emergence of new molecular biology techniques, such as high-throughput sequencing techniques, has completely changed the face of clinical microbiology. Their use in metagenomics research has also increased the speed and accuracy of diagnosis [1] and revealed changes in the bacterial composition throughout the disease [2,3]. One of the latest strategies involves the use of targeted sequencing for the bacterial 16S rRNA gene. This method has already gained many positive results in examining microbiota changes in such regions as the lungs, vagina, urinary tract, or gut [4,5,6,7]. In the Human Microbiome Project (HMP), led by four sequencing centers (the Broad Institute, Baylor College of Medicine, Washington University School of Medicine, and the J. Craig Venter Institute) and the Data Analysis and Coordination Center (DACC), microbial communities of 300 healthy individuals were characterized by 16S rRNA sequencing [8]. The overall goal was to determine if there was a healthy core microbiome, so samples were collected at multiple time points from several different places on the human body: nasal passages, mouth, skin, gastrointestinal tract, and genitourinary tract.
The results of 16S rRNA sequencing indicated a close connection between the microbiota in different compartments of an organism, which were formerly treated as independent niches with their microbiomes. Several recent studies have focused on the immunological interaction between the gut and lungs and found that the composition of their microbiomes is influenced by the diseases or homeostasis of the entire organism [5]. It was also found that bacteria could penetrate from one anatomical region to another, changing the composition of the microbiome in even anatomically distinct regions. For example, pathogenic bacteria typical for the oral cavity (Fusobacterium nucleatum and Porphyromonas gingivalis) can be found in the colon of patients suffering from periodontitis. These bacteria are highly pathogenic and can seriously alter the composition of the microbiome in the whole gastrointestinal tract, which may lead to chronic inflammation and even to colorectal cancer tumorigenesis [7]. A relationship between oral microbiota and cancer was also observed in another large study, in which the presence of oral bacteria Porphyromonas gingivalis and Aggregatibacter actinomycetemcomitans were associated with a higher risk of pancreatic cancer [9].
These relationships may also have an influence in situ as the composition of the oral microbiome is responsible for the health of the oral cavity itself, as the oral cavity is the first part of the gastrointestinal tract, and in adults has a microbiome of over 700 microorganisms, including bacteria, fungi, and viruses [10]. Bacteria residing in the oral cavity of a newborn come from people in their immediate surroundings. Colonization takes place in stages of child development to achieve some stabilization in adulthood. This balance may be disturbed by childhood diseases, such as early childhood caries, celiac disease, autism, Henoch-Schönlein purpura, appendicitis in children, inflammatory bowel disease, and obstructive sleep apnea syndrome in children [10].
Changes in the oral microbiome may also be of particular importance in patients with childhood malignant neoplasms who have a number of oral complications like caries, gingivitis, stomatitis, dry mouth, and candidiasis [11]. Patients receiving chemotherapy also had a different microbiome composition than healthy controls, which is a potential risk factor for the development of oral mucositis [11,12]. It was concluded that childhood cancer therapy might affect tooth development, saliva function, craniofacial development, and temporomandibular joint function, exposing some children who have survived cancer to an increased risk of oral and dental problems [13,14].
However, the role of the oral microbiome in cancer treatment and healing remains unclear, and relatively few studies have evaluated the oral microbiome of cancer patients with the use of 16S rRNA sequencing. The present study examines whether the salivary microbiome of childhood patients is altered over the course of anti-cancer chemotherapy.

2. Materials and Methods

2.1. The Cohort

A total of 71 children newly diagnosed with various types of cancer at the Department of Pediatrics, Oncology & Hematology, between March 2018 and February 2019, were approached to participate in the study. The schema used for participant selection is presented in Figure 1. The inclusion criteria comprised the following: age below eighteen years, newly diagnosed cancer before any treatment, no primary or secondary immunodeficiency identified, and no antibiotics treatment within six months prior to the study.
A total of 30 patients did not meet the inclusion criteria or refused to take part. Therefore, the first saliva sample was taken from 41 children (26 boys and 15 girls) aged 6–18 years. During the study, six patients died, eight had a relapse, four withdrew, and three missed the second sampling. Therefore, a second sample was taken from 20 participants: 7 (35%) girls and 13 (65%) boys. Paired samples from these 20 children were included in the statistical analysis. The patients suffered from the following cancer types: brain tumor (n = 4), leukemia (n = 4), Hodgkin’s lymphoma (n = 4), non-Hodgkin’s lymphoma (n = 3), Rhabdomyosarcoma (n = 2), neuroblastoma (n = 1), ovarian tumor (n = 1), and tumor of the testis (n = 1).
The median age was 13.79 years at the first saliva collection and 15.58 years during the second collection. The median duration of cancer treatment was 10.5 months (IQR: 4.5–24 months); (Table 1). Detailed clinical characteristics is provided in Table S1.
Written consent to take part was collected from parents or guardians, as well as from the child, before participation. The study was approved by the Bioethics Committee of the Medical University of Lodz (IRB number: RNN/37/13/KE).

2.2. Saliva Collection and DNA Isolation

The first saliva sample was collected on the first day after cancer diagnosis, before any anti-cancer treatment. The second saliva sample was taken when the anti-cancer treatment was finished, i.e., three weeks after the last anti-cancer drug was given, corresponding to a median of 22 (19–24) days after oncological treatment.
The non-stimulated saliva collection process supervised by a nurse took at least 15 min. A collection of 4 mL of saliva into a 15 mL sterile plastic vial (sputum tube, Sarstedt, Germany) was performed in the morning: 2–4 h after breakfast and teeth brushing. Saliva samples were stored immediately at −20 degrees Celsius for further genomic DNA extraction. Bacterial DNA was extracted from frozen 2 mL clinical saliva samples using 200 U mutanolysin and 50 µL lysozyme with a Genomic Midi AX DNA Extraction Kit (A&A Biotechnology, catalog number: 895-20-M1, Gdynia, Poland).

2.3. Library Preparation and Sequencing

Microbial community profiles were assessed by sequencing the 16S rRNA gene. The first step was the amplification of the V3 and V4 variable fragments of the analyzed gene according to the manufacturer’s protocol (Illumina, San Diego, CA, USA). Primers with overhanging adapters compatible with Illumina indexes and sequencing adapters in paired-end sequencing techniques were used. Kapa HiFi polymerase (Roche, Cape Town, South Africa) was used to amplify a fragment with a mean length of 464 bp. Next, the specificity of the obtained products was evaluated in an agarose gel and then purified on AMPure XP magnetic beads (Beckman Coulter, Brea, CA, USA). The indexing reactions were also carried out using Kapa HiFi polymerase (Roche, Cape Town, South Africa) with the Nextera XT dual-index set (Illumina, San Diego, CA, USA). The concentration of the obtained libraries was determined using the Qubit 2.0 device (Thermo Fisher Scientific, Waltham, MA, USA) and pooled in equal concentrations. The resulting library was sequenced on the MiSeq platform (Illumina, San Diego, CA, USA) using a MiSeq Reagent Kit v3 (600 cycles).

2.4. NGS Data Processing

Raw sequencing data from the MiSeq device were uploaded to the Galaxy web platform [15] and were analyzed using the public server at usegalaxy.org. FASTQ format files were unified to Sanger FASTQ encoding with the FASTQ Groomer tool [16]. Paired-end reads were first merged using the FLASH tool [17], then the Trimmomatic algorithm was used to remove adapters and low-quality reads (below Q20 value) [18]. Operational taxonomic units (OTUs) were assigned by the Kraken algorithm [19] and filtered by the classification confidence score at a level of 0.05. Read counts for each taxonomy level were extracted into tables, and the percentage abundance of each identified bacterial taxon was calculated. Only those taxonomic units that constituted at least 1% of mappable reads in at least one of the examined patients were considered for further analysis.

2.5. Data Analysis

2.5.1. Alpha and Beta Diversity

The compositional diversity of the microbiome in the samples was analyzed at the genus level on the web Microbiome Analyst platform [20]. The Chao1 algorithm and the paired t-test were used to determine the alpha-diversity of the microbiome and the statistical significance of any differences. Beta diversity was visualized using non-metric multidimensional scaling (NMDS), and the statistical significance was tested by the PERMANOVA method.

2.5.2. Visualization of Data

Hierarchical Clustering and Heatmap were created using a Ward clustering algorithm and Euclidean distance. The results were visualized as a dendrogram using the Ward clustering algorithm and Bray–Curtis index distance measure; the dendrogram was created using the web Microbiome Analyst platform [16]. Changes in raw abundance before and after anticancer therapy were presented using stacked-bar plots using ggplot2 package in R. Differences in relative abundances between time-points in groups were obtained with k-means clustering and plotted as chard plots using the cyclize package in R.

2.5.3. Statistical Analysis

Changes in the relative abundances were tested with Wilcoxon signed-rank. Patterns of change in oral microbiota were investigated using a k-means clustering algorithm. Differences between relative abundances between time-points were compared using Wilcoxon signed-rank test. Statistical analysis was performed using R-studio, Statistica 13.1 (TIBCO Software, Palo Alto, CA, USA), and Microsoft Excel software. Two-tailed p-values lower than 0.05 were deemed statistically significant. Due to the small study group, no multiple comparisons correction was applied.

3. Results

3.1. 16S rRNA Sequencing and Microbiome Abundance Analysis

The number of all obtained reads assigned to bacteria was 5,060,942, which translated into a mean value of 126,524 reads per single analyzed sample. The detailed analysis focused mainly on the genus level because it was the lowest level that had a high percentage of assigned reads, with a mean value of 89%. In the analyzed material, 20 different types of bacteria were identified, with the most abundant being Streptococcaceae 46.4%, Prevotellaceae 12.6%, and Veillonellaceae 9.9%, which together accounted for 68.9% of all identified OTUs.
From investigated bacteria, Tannerella, Candidatus Kinetoplastibacterium, Lactobacillus, Gardnerella, and Cellulosilyticum were determined to have a detection threshold below 0.001% in more than 50% of samples. However, because the study’s main aim was to determine the change in microbiome before and after anticancer therapy, it was decided not to exclude these genera, as they may be clinically relevant, despite their low prevalence (Figure 2b). Statistically significant differences were observed for the relative abundances of Cellulosilyticum (88% increase in relative abundance; p = 0.043) and Tannerella genera (97% increase in relative abundance; p = 0.036). Both genera, however, were very rare in the evaluated material. A trend towards significance was also noted for Bifidobacterium, Prevotella, and Kribbella (Table 2).
Alpha-diversity between time points (before and after anticancer therapy) reflects a change in microbiome diversity due to the treatment. However, the observed shift in alpha-diversity was not significant (p = 0.9456 for Chao1, 0.7303 for ACE, 0.8440 for Shannon, 0.5297 for Simpson, Figure 2a), indicating conserved microbiome diversity throughout the anticancer therapy (see also Figures S1 and S2).

3.2. Differences in Microbiome in Patients before and after Cancer Treatment

Despite a non-significant change in alpha-diversity metrics in the course of anticancer therapy, relative changes in genera abundance were further investigated. The microbiological profile of each patient, before and after anticancer therapy, is demonstrated in Figure 3a. For selected genera (p < 0.15, Table 2), overall changes in relative abundances were presented using stacked bar plots (Figure 3b–f).

3.3. Patterns of Microbiome Changes during Cancer Treatment Is Associated with Clinical Outcomes

To investigate the patterns of changes within patients, we applied k-means unsupervised clustering. Two clusters were identified, with significantly different changes in microbiome. Group median values and inter-quartile ranges for genera significantly altered by treatment in at least one group were provided in Table 3. For Cluster 1 we observed a 62% increase in Campylobacter and an 80% increase in Fusobacterium abundance. As for Cluster 2, there was a 242% decrease in Neisseria abundance.
The characteristics of changes due to anticancer therapy in the clusters are illustrated via circus plots (Figure 4a,b). To focus on the specific changes, selected genera are presented as a stacked bar plot (Figure 4c–j).
These two clusters differed with the age of patients at the diagnosis (p = 0.0049) and with the length of their therapies (p = 0.0038); (Table 4 and Figure 5a,b).
Although the length of antibiotic therapy was not significantly different between clusters, an abundance of specific bacteria correlated with the duration of antibiotic treatment. For the Gardnerella, we observed a strong positive correlation between the length of antibiotic therapy and its relative abundance after anticancer therapy (R = 0.4729, p = 0.0353, Figure 6a). Moreover, Cluster 2 presented a significant correlation between the duration of antibiotic therapy and change in relative abundance pre-post anticancer therapy for Gardnerella (R = −0.7280, p = 0.0260, Figure 6b) and Prevotella (R = 0.7699, p = 0.0152, Figure 6c).

4. Discussion

Anticancer therapy among pediatric patients consists of multidrug regimens and generates severe immunosuppression; these significantly increase the risk of infection and antibiotic use and, thus, affects the microbiota profile of the patient. Our results show that treatment of pediatric malignancies impacts the salivary microbiota. The length of intensive chemotherapy increases the risk of infection or empirical antibiotic treatment. Interestingly, unsupervised clustering analysis found these clusters to be linked with certain clinical features of the patients, including age at diagnosis and length of cancer therapy. This suggests that saliva microbiota might be modified by anticancer treatment in childhood.
While some changes in components were significant and substantial heterogeneity was observed between the patients, others generally remained constant despite the treatment and immunological suppression. The major component of the microbiome was constituted by Streptococcus, which was found to constitute 46.4% of the salivary microbiome. This seems to be in line with the prevalence of the genus in the Polish population and begins with the relationships among the mode of delivery (natural or C-section), the oral microbiome pattern of the infant, and the similarity of commensal bacteria inhabiting the oral cavity of the mother and her placenta [21,22]. Streptococcus was also found to be the predominant genus in a pre-dentate mouth (62.2%), later connected with the initiation of dental caries [23]. Other genera usually found to be abundant in the infant’s mouth are Veillonella, Neisseria, Rothia, Haemophilus, Gemella, Granulicatella, Leptotrichia, and Fusobacterium. This differs from the oral microbiota of adults, in whom Haemophilus, Neisseria, Veillonella, Fusobacterium, Oribacterium, Rothia, Treponema, and Actinomyces predominated [23]. In our group of children, aged between 6 and 18 years, the second and third most abundant genera were Prevotella and Veillonella while the remaining comprised 31.1% of the microbiome.
The adult microbiome shows greater variety than those found in children, with the greatest changes taking place up to the age of three years [22]. Teeth gradually erupt, occlusions form, and the diet becomes more complex. As the process of immune postnatal development is closely related to microbial exposure in the gut, diet plays a crucial role in its proper development [24]. In general, the diet alters the structure and function of the microbiome throughout the gastrointestinal tract, which in some cases may lead to a decrease in species richness and promote the proliferation of pathogenic bacterial taxa [25]. Such complications appear when a child suffers from a condition that impairs proper nutrition, such as celiac disease or inflammatory bowel disease [26,27]. Indeed, during cancer therapy, children may often require changes in nutrition or even a need for parenteral feeding [28]. It has also been observed that oncological treatment may influence the outcome of therapy by causing deviations in gut microbiota composition [29].
In the investigated group, a significant difference between the two time-points was found only for Cellulosilyticum and Tannerella. Cellulosilyticum is a newly described genus that could efficiently hydrolyze cellulose, and it is mostly identified in the rumen of ruminantia [30]. In contrast, Tannerella is quite commonly detected in the oral cavity of humans and its presence is widely connected with periodontitis and loss of teeth in adults [31]. It was also found that Tannerella was overrepresented in the tongue coating in patients with pancreatic head cancer [32].
The abundance of Prevotella bacteria was also found to decrease over treatment. This genus is known to be associated with periodontal inflammation, and its abundance has been found to increase in patients with asthma or chronic obstructive pulmonary disease [33]. A trend toward significance was also noted in the abundance of Bifidobacterium, whose fraction increased two-fold after the treatment. These bacteria colonize the gastrointestinal tract in children and they are believed to have a range of health-promoting properties [34].
Finally, Kribbella also demonstrated considerable changes in the number during therapy. The genus is usually isolated from the soil or plants [35]. This is the first confirmation of its constant presence in the oral microbiome of humans.
To identify the existence of specific microbiome patterns, the patients were clustered with regard to the observed changes in oral microbiome. Two clusters were identified, characterized by changes in Campylobacter and Fusobacterium (Cluster 1), or Neisseria (Cluster 2).
Campylobacter is the most common cause of diarrheal illness, including acute enteritis, extraintestinal infections, and post-infectious complications [36], and its increased presence can be easily understood in the microbiome of children treated from cancer. In contrast, while Fusobacterium is considered a part of the normal microbiome of the oral cavity, gastrointestinal tract, and female genital tract, it is believed to potentially cause acute otitis media in small children [37]. Increases in the relative abundance of Fusobacterium were also confirmed in the colon of colorectal cancer patients [38] and several studies indicate a relationship between Fusobacterium, periodontal disease, and oral cancer in adults [39,40].
Regarding Cluster 2, the patients demonstrated a significant decrease in the abundance of Neisseria. Some Neisseria species are invasive pathogens inhabiting mucosal surfaces in the upper respiratory and genito-urinary tracts, but most are considered part of the commensal microbiome, which is needed to prevent infection in the oropharyngeal area [41]. However, it has been shown that high levels of Fusobacterium periodonticum and low levels of Neisseria mucosa are specific risk factors for pancreatic cancer [42].
The patients from Cluster 2 were older than children from Cluster 1 and their cancer therapy length was significantly shorter. No significant differences were observed in patients with different durations of antibiotic therapy or taking steroids. However, the one whose abundance was related to the length of antibiotic treatment was Gardnerella. When it comes to cluster analysis, there are relationships between the duration of antibiotic therapy and changes in the pre-post percentage of a given bacterium. A short course of antibiotics was correlated with an increase in Gardnerella abundance in Cluster 1, but the longer the course, the smaller the change or even increase in abundance.
This is in line with the results of others, which show that in men treated with antibiotics for urethritis, the Gardnerella abundance in oral saliva did not change significantly [43].
This was also the case for Prevotella, even though its change in abundance was not significant after clustering, the Prevotella count in Cluster 2 patients was related to the duration of antibiotic treatment and was proportional to the duration of antibiotic treatment.
There were also no changes in patients with different types of cancer or the presence of caries. The non-significance may be related to the small number of patients in the groups and this can generally be considered the main limitation of the study. As such, our results should be interpreted tentatively and further studies with a larger sample size are needed. Even so, our study is the first to alter the oral microbiome in childhood cancer survivors using 16S rRNA genotyping.

5. Conclusions

Our results clearly indicate the feasibility of measuring changes in oral microbiome in the course of chemotherapy and provide insight into the changes of oral microbiome in children in course of therapy for cancer. However, further studies involving greater cohorts are necessary, and these should include factors such as dietary habits and proven infections throughout the course of chemotherapy.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers14010007/s1, Table S1: Detailed cancer status in the study group; Figure S1: Alpha-diversity measure using Simpson at Genus level across all the samples; Figure S2: Alpha-diversity measure using Shannon at Genus level across all the samples.

Author Contributions

Conceptualization, W.M., J.S. and P.P.; methodology, W.M. and M.B.; software, W.F. and T.P.; validation, W.M. and M.B.; formal analysis, W.F. and J.C.; investigation, P.P., M.Z. and M.S. (Małgorzata Stolarska); resources, P.P., M.B., M.S. (Małgorzata Szybka) and B.Z.; data curation, T.P.; writing—original draft preparation, P.P. and W.M.; writing—review and editing, P.P.; visualization, J.C.; supervision, W.M.; project administration, P.P.; funding acquisition, J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Bioethics Committee of the Medical University of Lodz, Poland (IRB number: RNN/37/13/KE).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author and in Bioproject with Accession: PRJNA770171 https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA770171. Last access 12 December 2021.

Acknowledgments

We would like to thank Małgorzata Brauncjas who did bacteria identification.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Greninger, A.L.; Messacar, K.; Dunnebacke, T.; Naccache, S.N.; Federman, S.; Bouquet, J.; Mirsky, D.; Nomura, Y.; Yagi, S.; Glaser, C.; et al. Clinical metagenomic identification of Balamuthia mandrillarisencephalitis and assembly of the draft genome: The continuing case for reference genome sequencing. Genome Med. 2015, 7, 113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Babaev, E.A.; Balmasova, I.P.; Mkrtumyan, A.M.; Kostryukova, S.N.; Vakhitova, E.S.; Ilina, E.N.; Tsarev, V.N.; Gabibov, A.G.; Arutyunov, S.D. Metagenomic Analysis of Gingival Sulcus Microbiota and Pathogenesis of Periodontitis Associated with Type 2 Diabetes Mellitus. Bull. Exp. Biol. Med. 2017, 163, 718–721. [Google Scholar] [CrossRef] [PubMed]
  3. Peng, W.; Yi, P.; Yang, J.; Xu, P.; Wang, Y.; Zhang, Z.; Huang, S.; Wang, Z.; Zhang, C. Association of gut microbiota composition and function with a senescence-accelerated mouse model of Alzheimer’s Disease using 16S rRNA gene and metagenomic sequencing analysis. Aging 2018, 10, 4054–4065. [Google Scholar] [CrossRef] [PubMed]
  4. Virtanen, S.; Kalliala, I.; Nieminen, P.; Salonen, A. Comparative analysis of vaginal microbiota sampling using 16S rRNA gene analysis. PLoS ONE 2017, 12, e0181477. [Google Scholar] [CrossRef] [Green Version]
  5. Wang, J.; Li, F.; Tian, Z. Role of microbiota on lung homeostasis and diseases. Sci. China Life Sci. 2017, 60, 407–1415. [Google Scholar] [CrossRef]
  6. Chen, C.; Song, X.; Wei, W.; Zhong, H.; Dai, J.; Lan, Z.; Li, F.; Yu, X.; Feng, Q.; Wang, Z.; et al. The microbiota continuum along the female reproductive tract and its relation to uterine-related diseases. Nat. Commun. 2017, 17, 875. [Google Scholar] [CrossRef] [Green Version]
  7. Koliarakis, I.; Messaritakis, I.; Nikolouzakis, T.K.; Hamilos, G.; Souglakos, J.; Tsiaoussis, J. Oral Bacteria and Intestinal Dysbiosis in Colorectal Cancer. Int. J. Mol. Sci. 2019, 20, 4146. [Google Scholar] [CrossRef] [Green Version]
  8. Integrative HMP (iHMP) Research Network Consortium. The Integrative Human Microbiome Project. Nature 2019, 569, 641–648. [Google Scholar]
  9. Fan, X.; Alekseyenko, A.V.; Wu, J.; Peters, B.A.; Jacobs, E.J.; Gapstur, S.M.; Purdue, M.P.; Abnet, C.C.; Stolzenberg-Solomon, R.; Miller, G.; et al. Human oral microbiome and prospective risk for pancreatic cancer: A population-based nested case-control study. Gut 2018, 67, 120–127. [Google Scholar] [CrossRef] [Green Version]
  10. Xiao, J.; Fiscella, K.A.; Gill, S.R. Oral microbiome: Possible harbinger for children’s health. Int. J. Oral Sci. 2020, 12, 12. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, Y.; Zeng, X.; Yang, X.; Que, J.; Du, Q.; Zhang, Q.; Zou, J. Oral Health, Caries Risk Profiles, and Oral Microbiome of Pediatric Patients with Leukemia Submitted to Chemotherapy. BioMed Res. Int. 2021, 2021, 6637503. [Google Scholar] [CrossRef] [PubMed]
  12. de Farias Gabriel, A.; Silveira, F.M.; Curra, M.; Schuch, L.F.; Wagner, V.P.; Martins, M.A.T.; da Silveira Matte, U.; Siebert, M.; Botton, M.R.; Brunetto, A.T.; et al. Risk factors associated with the development of oral mucositis in pediatric oncology patients: Systematic review and meta-analysis. Oral Dis. 2021. epub ahead of print. [Google Scholar] [CrossRef]
  13. Proc, P.; Szczepańska, J.; Skiba, A.; Zubowska, M.; Fendler, W.; Młynarski, W. Dental Anomalies as Late Adverse Effect among Young Children Treated for Cancer. Cancer Res. Treat. 2016, 48, 658–667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Effinger, K.E.; Migliorati, C.A.; Hudson, M.M.; McMullen, K.P.; Kaste, S.C.; Ruble, K.; Guilcher, G.M.; Shah, A.J.; Castellino, S.M. Oral and dental late effects in survivors of childhood cancer: A Children’s Oncology Group report. Support. Care Cancer 2014, 22, 2009–2019. [Google Scholar] [CrossRef] [Green Version]
  15. Afgan, E.; Baker, D.; van den Beek, M.; Blankenberg, D.; Bouvier, D.; Cech, M.; Chilton, J.; Clements, D.; Coraor, N.; Eberhard, C.; et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic Acids Res. 2016, 44, W3–W10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Blankenberg, D.; Gordon, A.; Von Kuster, G.; Coraor, N.; Taylor, J.; Nekrutenko, A.; Galaxy, T. Manipulation of FASTQ data with Galaxy. Bioinformatics 2010, 26, 1783–1785. [Google Scholar] [CrossRef] [PubMed]
  17. Magoc, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  18. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  19. Wood, D.E.; Salzberg, S.L. Kraken: Ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014, 15, R46. [Google Scholar] [CrossRef] [Green Version]
  20. Dhariwal, A.; Chong, J.; Habib, S.; King, I.L.; Agellon, L.B.; Xia, J. Microbiome Analyst: A web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res. 2017, 45, W180–W188. [Google Scholar] [CrossRef]
  21. Bearfield, C.; Davenport, E.S.; Sivapathasundaram, V.; Allaker, R.P. Possible association between amniotic fluid micro-organism infection and microflora in the mouth. BJOG Int. J. Obstet. Gynaecol. 2002, 109, 527–533. [Google Scholar] [CrossRef]
  22. Lif Holgerson, P.; Harnevik, L.; Hernell, O.; Tanner, A.C.; Johansson, I. Mode of birth delivery affects oral microbiota in infants. J. Dent. Res. 2011, 90, 1183–1188. [Google Scholar] [CrossRef] [Green Version]
  23. Cephas, K.D.; Kim, J.; Mathai, R.A.; Barry, K.A.; Dowd, S.E.; Meline, B.S.; Swanson, K.S. Comparative analysis of salivary bacterial microbiome diversity in edentulous infants and their mothers or primary care givers using pyrosequencing. PLoS ONE 2011, 6, e23503. [Google Scholar] [CrossRef] [PubMed]
  24. Munyaka, P.M.; Khafipour, E.; Ghia, J.E. External influence of early childhood establishment of gut microbiota and subsequent health implications. Front. Pediatr. 2014, 9, 109. [Google Scholar] [CrossRef] [Green Version]
  25. Shen, T.D. Diet and Gut Microbiota in Health and Disease. Intest. Microbiome Funct. Asp. Health Dis. 2017, 88, 117–126. [Google Scholar] [CrossRef] [Green Version]
  26. Francavilla, R.; Ercolini, D.; Piccolo, M.; Vannini, L.; Siragusa, S.; De Filippis, F.; De Pasquale, I.; Di Cagno, R.; Di Toma, M.; Gozzi, G.; et al. Salivary microbiota and metabolome associated with celiac disease. Appl. Environ. Microbiol. 2014, 80, 3416–3425. [Google Scholar] [CrossRef] [Green Version]
  27. Docktor, M.J.; Paster, B.J.; Abramowicz, S.; Ingram, J.; Wang, Y.E.; Correll, M.; Jiang, H.; Cotton, S.L.; Kokaras, A.S.; Bousvaros, A. Alterations in diversity of the oral microbiome in pediatric inflammatory bowel disease. Inflamm. Bowel Dis. 2012, 18, 935–942. [Google Scholar] [CrossRef] [PubMed]
  28. Brinksma, A.; Roodbol, P.F.; Sulkers, E.; Kamps, W.A.; de Bont, E.S.; Boot, A.M.; Burgerhof, J.G.; Tamminga, R.Y.; Tissing, W.J. Changes in nutritional status in childhood cancer patients: A prospective cohort study. Clin. Nutr. 2015, 34, 66–73. [Google Scholar] [CrossRef]
  29. Aarnoutse, R.; de Vos-Geelen, J.M.P.G.M.; Penders, J.; Boerma, E.G.; Warmerdam, F.A.R.M.; Goorts, B.; Olde Damink, S.W.M.; Soons, Z.; Rensen, S.S.M.; Smidt, M.L. Study protocol on the role of intestinal microbiota in colorectal cancer treatment: A pathway to personalized medicine 2.0. Int. J. Colorectal Dis. 2017, 32, 1077–1084. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Li, S.; Shao, N.; Luo, Y.; Liu, H.; Cai, S.; Dong, X. Transcriptome and Zymogram Analyses Reveal a Cellobiose-Dose Related Reciprocal Regulatory Effect on Cellulase Synthesis in Cellulosilyticum ruminicola H1. Front. Microbiol. 2017, 12, 2497. [Google Scholar] [CrossRef] [Green Version]
  31. Chen, W.P.; Chang, S.H.; Tang, C.Y.; Liou, M.L.; Tsai, S.J.; Lin, Y.L. Composition Analysis and Feature Selection of the Oral Microbiota Associated with Periodontal Disease. BioMed Res. Int. 2018, 2018, 3130607. [Google Scholar] [CrossRef] [Green Version]
  32. Lu, H.; Ren, Z.; Li, A.; Li, J.; Xu, S.; Zhang, H.; Jiang, J.; Yang, J.; Luo, Q.; Zhou, K.; et al. Tongue coating microbiome data distinguish patients with pancreatic head cancer from healthy controls. J. Oral Microbiol. 2019, 11, 1563409. [Google Scholar] [CrossRef] [PubMed]
  33. Larsen, J.M. The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology 2017, 151, 363–374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. O’Callaghan, A.; van Sinderen, D. Bifidobacteria and Their Role as Members of the Human Gut Microbiota. Front. Microbiol. 2016, 15, 925. [Google Scholar] [CrossRef] [Green Version]
  35. Ozdemir-Kocak, F.; Saygin, H.; Saricaoglu, S.; Cetin, D.; Guven, K.; Spröer, C.; Schumann, P.; Klenk, H.P.; Sahin, N.; Isik, K. Kribbella soli sp. nov., isolated from soil. Antonie Leeuwenhoek 2017, 110, 641–649. [Google Scholar] [CrossRef] [PubMed]
  36. Fitzgerald, C. Campylobacter. Clin. Lab. Med. 2015, 35, 289–298. [Google Scholar] [CrossRef]
  37. Arane, K.; Goldman, R.D. Fusobacterium infections in children. Can. Fam. Physician 2016, 62, 813–814. [Google Scholar]
  38. Amitay, E.L.; Werner, S.; Vital, M.; Pieper, D.H.; Höfler, D.; Gierse, I.J.; Butt, J.; Balavarca, Y.; Cuk, K.; Brenner, H. Fusobacterium and colorectal cancer: Causal factor or passenger? Results from a large colorectal cancer screening study. Carcinogenesis 2017, 1, 781–788. [Google Scholar] [CrossRef] [Green Version]
  39. Gholizadeh, P.; Eslami, H.; Yousefi, M.; Asgharzadeh, M.; Aghazadeh, M.; Kafil, H.S. Role of oral microbiome on oral cancers, a review. Biomed. Pharmacother. 2016, 84, 552–558. [Google Scholar] [CrossRef]
  40. Takahashi, Y.; Park, J.; Hosomi, K.; Yamada, T.; Kobayashi, A.; Yamaguchi, Y.; Iketani, S.; Kunisawa, J.; Mizuguchi, K.; Maeda, N.; et al. Analysis of oral microbiota in Japanese oral cancer patients using 16S rRNA sequencing. J. Oral Biosci. 2019, 61, 120–128. [Google Scholar] [CrossRef]
  41. Liu, G.; Tang, C.M.; Exley, R.M. Non-pathogenic Neisseria: Members of an abundant, multi-habitat, diverse genus. Microbiology 2015, 161, 1297–1312. [Google Scholar] [CrossRef] [PubMed]
  42. Sun, H.; Zhao, X.; Zhou, Y.; Wang, J.; Ma, R.; Ren, X.; Wang, H.; Zou, L. Characterization of Oral Microbiome and Exploration of Potential Biomarkers in Patients with Pancreatic Cancer. BioMed Res. Int. 2020, 47, 12498. [Google Scholar] [CrossRef] [PubMed]
  43. Le, P.T.; Hamasuna, R.; Matsumoto, M.; Furubayashi, K.; Hatanaka, M.; Kawai, S.; Yamaguchi, T.; Uehara, K.; Murakami, N.; Yoshioka, M.; et al. The detection of microorganisms related to urethritis from the oral cavity of male patients with urethritis. J. Infect. Chemother. 2017, 23, 668–673. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Plan of the study—number of patients enrolled in each stage of the study.
Figure 1. Plan of the study—number of patients enrolled in each stage of the study.
Cancers 14 00007 g001
Figure 2. Changes in salivary microbiome before and after anticancer treatment: (a) Alpha diversity measures for matched patient samples (PAT(#ID)(before/after)) using Chao1 index (red—before, blue—after); (b) Heatmap of taxon prevalence by its relative abundance, cell color represents a relative prevalence of taxa at certain detection threshold (from 0 to 1, blue to red).
Figure 2. Changes in salivary microbiome before and after anticancer treatment: (a) Alpha diversity measures for matched patient samples (PAT(#ID)(before/after)) using Chao1 index (red—before, blue—after); (b) Heatmap of taxon prevalence by its relative abundance, cell color represents a relative prevalence of taxa at certain detection threshold (from 0 to 1, blue to red).
Cancers 14 00007 g002
Figure 3. Changes in salivary microbiome before and after anticancer treatment for specific patient: (a) paired measurements of relative abundances of genera for a specific patient, before and after anticancer therapy; (bf) stacked bar plots for genus relative abundance before and after anticancer therapy. Layer color indicated specific patient samples (across bf graphs). The demonstrated genera were significantly different between paired measures (before and after therapy), p-value < 0.15 using Wilcoxon signed-rank test.
Figure 3. Changes in salivary microbiome before and after anticancer treatment for specific patient: (a) paired measurements of relative abundances of genera for a specific patient, before and after anticancer therapy; (bf) stacked bar plots for genus relative abundance before and after anticancer therapy. Layer color indicated specific patient samples (across bf graphs). The demonstrated genera were significantly different between paired measures (before and after therapy), p-value < 0.15 using Wilcoxon signed-rank test.
Cancers 14 00007 g003
Figure 4. Clustering analysis of microbiome changes after anti-cancer therapy: (a,b) chord plots represent the relative abundance estimated for a perfect profile in Cluster 1 and Cluster 2. Values before and after therapy are right- and left side of a circle respectively. Colors for both chord plots correspond to the genus (legend “GENUS”); (cj). Stacked bar plots represent the difference in relative abundance specific for clusters (Cluster 1 = c,d,g,h; Cluster 2 = e,f,i,j). Color layers represent patients for a specific group and they are consequently used within a group.
Figure 4. Clustering analysis of microbiome changes after anti-cancer therapy: (a,b) chord plots represent the relative abundance estimated for a perfect profile in Cluster 1 and Cluster 2. Values before and after therapy are right- and left side of a circle respectively. Colors for both chord plots correspond to the genus (legend “GENUS”); (cj). Stacked bar plots represent the difference in relative abundance specific for clusters (Cluster 1 = c,d,g,h; Cluster 2 = e,f,i,j). Color layers represent patients for a specific group and they are consequently used within a group.
Cancers 14 00007 g004
Figure 5. Clinical diversity of identified clusters: (a) difference in the age at diagnosis of patients from Cluster 1 and Cluster 2; p = 0.0049. (b) Difference in the length of anti-cancer therapy between patients from Cluster 1 and Cluster 2; p = 0.0038.
Figure 5. Clinical diversity of identified clusters: (a) difference in the age at diagnosis of patients from Cluster 1 and Cluster 2; p = 0.0049. (b) Difference in the length of anti-cancer therapy between patients from Cluster 1 and Cluster 2; p = 0.0038.
Cancers 14 00007 g005
Figure 6. Spearman correlations between duration (time) of antibiotic therapy (days) and relative abundance of Gardnerella before (a), after (b), and Prevotella change (c) for Cluster 1 (continuous line, circle markers) and Cluster 2 (dashed line, square markers); 95% confidence interval for regression line was marked as dotted lines.
Figure 6. Spearman correlations between duration (time) of antibiotic therapy (days) and relative abundance of Gardnerella before (a), after (b), and Prevotella change (c) for Cluster 1 (continuous line, circle markers) and Cluster 2 (dashed line, square markers); 95% confidence interval for regression line was marked as dotted lines.
Cancers 14 00007 g006
Table 1. Clinical characteristics of the study group (n = 20).
Table 1. Clinical characteristics of the study group (n = 20).
Clinical FeatureMedian (25–75%) or % (n)
Age at diagnosis (years)13.79 (9.13–16.08)
Length of therapy (months)10.50 (4.50–24.00)
Length of antibiotic therapy (days)28.00 (8.00–96.50)
Sex (%males)65% (13/20)
Table 2. Comparison of changes in relative abundance before and after the anticancer therapy.
Table 2. Comparison of changes in relative abundance before and after the anticancer therapy.
BacteriaBefore Treatment (%)After Treatment (%)p-Value
Bifidobacterium0.37 (0.28–0.68)0.72 (0.37–1.00)0.1005
Cellulosilyticum0.0021 (0.0007–0.0029)0.0011 (0.0007–0.0012)0.0438
Kribbella0.48 (0.32–0.68)0.80 (0.50–1.21)0.1454
Prevotella12.28 (6.32–17.69)9.47 (4.93–12.94)0.1259
Tannerella0.10 (0.05–0.43)0.09 (0.01–0.27)0.0366
Table 3. Genera abundance summary after unsupervised classification (using k-means), based on patient-specific relative abundances (before and after anticancer therapy). Campylobacter, Fusobacterium, and Neisseria demonstrate a difference in abundance before and after therapy in at least one of the predefined groups. Change in abundance is compared using Wilcoxon signed-rank test, p < 0.15 is defined as significant.
Table 3. Genera abundance summary after unsupervised classification (using k-means), based on patient-specific relative abundances (before and after anticancer therapy). Campylobacter, Fusobacterium, and Neisseria demonstrate a difference in abundance before and after therapy in at least one of the predefined groups. Change in abundance is compared using Wilcoxon signed-rank test, p < 0.15 is defined as significant.
BacteriaBefore
Treatment
(Cluster 1)
After
Treatment
(Cluster 1)
p-Value Before
Treatment
(Cluster 2)
After
Treatment
(Cluster 2)
p-Value
Campylobacter0.16 (0.05–0.20)0.36 (0.20–0.63)0.01280.66 (0.45–0.76)0.38 (0.29–0.44)0.1097
Cellulosilyticum0.0014
(0.0000–0.0029)
0.0009
(0.0006–0.0012)
0.53370.0029
(0.0021–0.0030)
0.0011
(0.0009–0.0012)
0.0506
Fusobacterium0.09 (0.01–0.33)1.52 (0.15–1.81)0.02080.85 (0.26–1.61)0.71 (0.41–0.91)0.5940
Gardnerella0.01 (0.01–0.03)0.06 (0.01–0.14)0.18230.03 (0.01–0.05)0.02 (0.01–0.03)0.1386
Kribbella0.34 (0.27–0.58)0.64 (0.45–1.18)0.13070.65 (0.49–0.71)1.05 (0.69–1.23)0.1731
Megasphaera0.12 (0.01–0.42)0.08 (0.01–0.22)0.09120.49 (0.13–0.71)0.31 (0.16–0.60)0.5940
Neisseria2.28 (0.08–7.18)1.21 (0.06–5.40)0.47691.45 (0.69–2.65)0.09 (0.01–0.73)0.0284
Table 4. Clinical characteristics of the patients according to the cluster adherence.
Table 4. Clinical characteristics of the patients according to the cluster adherence.
Cluster 1Cluster 2
Clinical FeatureMedian (25–75%) or % (n)Median (25–75%) or % (n)p-Value
Age at diagnosis (years)10.25 (6.75–14.00)16.16 (13.91–17.66)0.0049
Length of anticancer therapy (months)19.00 (13.00–32.00)4.00 (4.00–7.00)0.0038
Length of antibiotic therapy (days)71.00 (9.00–108.00)26.00 (7.00–35.00)0.4467
Sex (% males)63.64% (7/11)66.67% (6/9)0.7415
Steroid therapy (Yes)27.27% (3/11)66.67% (6/9)0.1902
Diagnosis: leukemia/lymphoma45.45% (5/11)66.67% (6/9)0.8812
Diagnosis: CNS tumor27.27% (3/11)11.11% (1/9)
Diagnosis: soft tissue tumor27.27% (3/11)22.22% (2/9)
Caries (Yes)63.64% (7/11)55.56% (5/9)0.7136
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Proc, P.; Szczepańska, J.; Zarzycka, B.; Szybka, M.; Borowiec, M.; Płoszaj, T.; Fendler, W.; Chrzanowski, J.; Zubowska, M.; Stolarska, M.; et al. Evaluation of Changes to the Oral Microbiome Based on 16S rRNA Sequencing among Children Treated for Cancer. Cancers 2022, 14, 7. https://doi.org/10.3390/cancers14010007

AMA Style

Proc P, Szczepańska J, Zarzycka B, Szybka M, Borowiec M, Płoszaj T, Fendler W, Chrzanowski J, Zubowska M, Stolarska M, et al. Evaluation of Changes to the Oral Microbiome Based on 16S rRNA Sequencing among Children Treated for Cancer. Cancers. 2022; 14(1):7. https://doi.org/10.3390/cancers14010007

Chicago/Turabian Style

Proc, Patrycja, Joanna Szczepańska, Beata Zarzycka, Małgorzata Szybka, Maciej Borowiec, Tomasz Płoszaj, Wojciech Fendler, Jędrzej Chrzanowski, Małgorzata Zubowska, Małgorzata Stolarska, and et al. 2022. "Evaluation of Changes to the Oral Microbiome Based on 16S rRNA Sequencing among Children Treated for Cancer" Cancers 14, no. 1: 7. https://doi.org/10.3390/cancers14010007

APA Style

Proc, P., Szczepańska, J., Zarzycka, B., Szybka, M., Borowiec, M., Płoszaj, T., Fendler, W., Chrzanowski, J., Zubowska, M., Stolarska, M., & Młynarski, W. (2022). Evaluation of Changes to the Oral Microbiome Based on 16S rRNA Sequencing among Children Treated for Cancer. Cancers, 14(1), 7. https://doi.org/10.3390/cancers14010007

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop