Skip to main content
  • Research article
  • Open access
  • Published:

Comparison of the salivary and dentinal microbiome of children with severe-early childhood caries to the salivary microbiome of caries-free children

Abstract

Background

The main objectives of this study were to describe and compare the microbiota of 1) deep dentinal lesions of deciduous teeth of children affected with severe early childhood caries (S-ECC) and 2) the unstimulated saliva of these children and 3) the unstimulated saliva of caries-free children, and to compare microbiota compositional differences and diversity of taxa in these sampled sites.

Methods

Children with S-ECC and without S-ECC were recruited. The saliva of all children with and without S-ECC was sampled along with the deep dentinal microbiota from children affected by S-ECC. The salivary microbiota of children affected by S-ECC (n = 68) was compared to that of caries-free children (n = 70), by Illumina MiSeq sequencing of 16S rRNA amplicons. Finally, the caries microbiota of deep dentinal lesions of those children with S-ECC was investigated.

Results

Using two beta diversity metrics (Bray Curtis dissimilarity and UniFrac distance), the caries microbiota was found to be distinct from that of either of the saliva groups (caries-free & caries-active) when bacterial abundance was taken into account. However, when the comparison was made by measuring only presence and absence of bacterial taxa, all three microbiota types separated. While the alpha diversity of the caries microbiota was lowest, the diversity difference between the caries samples and saliva samples was statistically significant (p < 0.001). The major phyla of the caries active dentinal microbiota were Firmicutes (median abundance value 33.5%) and Bacteroidetes (23.2%), with Neisseria (10.3%) being the most abundant genus, followed by Prevotella (10%). The caries-active salivary microbiota was dominated by Proteobacteria (median abundance value 38.2%) and Bacteroidetes (27.8%) with the most abundant genus being Neisseria (16.3%), followed by Porphyromonas (9.5%). Caries microbiota samples were characterized by high relative abundance of Streptococcus mutans, Prevotella spp., Bifidobacterium and Scardovia spp.

Conclusions

Distinct differences between the caries microbiota and saliva microbiota were identified, with separation of both salivary groups (caries-active and caries-free) whereby rare taxa were highlighted. While the caries microbiota was less diverse than the salivary microbiota, the presence of these rare taxa could be the difference between health and disease in these children.

Peer Review reports

Background

Dental caries is the “single most common chronic disease of childhood” [1] and affects 60–90% of all school children [2]. Severe early childhood caries (S-ECC) is an aggressive form of dental caries and is classified by the presence of a decayed, missing (due to caries), or filled tooth (dmft) index score of ≥4 (age 3), ≥5 (age 4), or ≥ 6 (age 5) [3]. S-ECC is destructive [4], and when it progresses it can cause acute pain and sepsis, and potential tooth loss [5]. Poor dental health in early childhood can interfere with the child’s quality of life, nutrition and school participation [4, 6]. Because of the young age of the children, S-ECC is difficult to treat successfully in the dental chair [7,8,9] and these children frequently require treatment under general anaesthesia [10] which increases treatment costs [11]. S-ECC is a risk factor for caries of permanent teeth [12,13,14] and affected children are at a higher risk of developing recurrent caries [8, 15, 16].

Dental caries has been investigated for many years using selective culture-based methods, and the role of Streptococcus and the presence of Lactobacillus have both long been recognized as playing substantial roles in dental caries [17,18,19,20,21,22,23]. Other non-microbiological risk factors that can play a significant role in S-ECC are host factors, environmental factors and diet [24, 25]. Dental caries is a product of loss of tooth structure and driven by acid produced by certain oral bacteria that ferment carbohydrate substrates, so diet plays a major role in the abundances of relevant microbes [26]. The impact of diet on the gut microbiome has been described [27, 28], but its impact on the oral cavity is not so well understood. Recent advances in molecular methods have allowed scientists to study the microbiology of oral disease with greater power, with technologies [29] such as 16S rRNA gene amplicon sequencing (including MiSeq and HiSeq platforms) showing that the microbiology of dental caries is much richer than previously believed.

The oral cavity has been shown to harbour more than 700 bacterial taxa, with one third of these described as non-culturable in vitro [30, 31], and two thirds belonging to cultivable species [32]. Although there is compositional variation between sample sites taken from the oral cavity, a ‘core’ microbiome in health has been identified [33]. Studies have also demonstrated that oral disease is not due to an isolated organism such as Streptococcus mutans causing caries, but is more polymicrobial in nature [23, 29,30,31,32,33,34]. Studies have identified Bifidobacterium, Veillonella, Granulicatetta, Scardovia, Fusobacterium, Prevotella and Actinomyces as potential contributors to ECC evidenced by their altered abundance in the caries microbiota [8, 12, 29, 30, 35,36,37,38,39].

To understand the microbiology of dental caries, it is helpful to use the combined findings from molecular and culture-based studies [37, 40,41,42,43,44] because molecular methods, when compared to culture-dependent methods have been shown to underestimate the proportions of certain phyla such as Actinobacteria [41, 45]. When examined in more detail, primer design was shown not to influence this bias [45]. 16S rRNA is highly conserved and the variable regions are high in G + C content. These high G + C regions in the DNA of Actinobacteria can result in an interruption of Taq polymerase processivity during the PCR process [41]. Primers used for this study targeted the V4 - V5 region, which has shown high classification accuracy when compared to other regional primer sites [46]. Studies have shown differences in the oral microbiota in children with caries and those that are caries free [21, 35, 37, 40]. Streptococcus mutans is found at higher levels at early stages of caries [47], while Lactobacillus sp. are associated with disease progression of caries and Scardovia sp. have been isolated from dentinal caries and have previously been associated with having a role in the advancement of deep caries in S-ECC [40, 41, 48, 49].

Limited data are available by molecular methods on the microbiota of advanced deep dentinal caries and pulpal infections of deciduous teeth, knowledge of which could aid in the development of anti-bacterial medicaments in pulp therapy of these affected teeth. Gram-negative species have been identified in deep dentinal caries of ECC affected teeth [43, 44], and are present in deep pulpal infections of primary teeth [50,51,52,53,54,55]. The microbiota of exposed vital pulp chambers of carious deciduous teeth was found to be dominated by phyla Firmicutes and Actinobacteria using anaerobic culture and by comparison with taxa in the Human Oral Microbiome Database [56]. Rôças et al. [57] using Illumina Miseq sequencing identified these same phyla to dominate the microbiota of deep dentinal caries of permanent teeth with irreversible pulpitis.

In the present study, Illumina Miseq sequencing was employed to compare the microbiota of the deep dentinal lesions of S-ECC affected deciduous teeth, and saliva of these caries-active children, with the saliva of caries-free children. S-ECC is an aggressive form of caries, and we analysed the microbiota of the deep dentinal caries of deciduous teeth to determine if the salivary microbiota was a reservoir or source of taxa linked with this form of caries. The main objectives were to describe and compare the microbiota of 1) deep dentinal lesions of deciduous teeth of children affected with S-ECC and 2) the unstimulated saliva of these children and finally, 3) the unstimulated saliva of caries-free children and compared compositional differences and diversity of taxa in these sampled sites.

Materials and methods

Study design, ethics and recruitment

The study design was to recruit two cohorts of children under the age of 60 months. These two groups were categorized into those with S-ECC (caries-active), and a caries-free cohort, all medically healthy. Sample size calculations were estimated based on previous similar studies [36, 47, 58, 59]. In total, we recruited 68 caries-active and 70 caries-free children. The deep dentinal lesion microbiota were labelled as Caries-active cavity (CAC) and the salivary microbiota of these caries-active children labelled as Caries-active saliva (CAS), while the saliva of the caries-free children were labelled as Caries-free saliva (CFS). Of the caries-active, all were S-ECC affected deciduous teeth, and the CAC and CAS are paired samples, each from the same subject. Ethical approval was obtained from the Teaching Hospitals Clinical Research Ethics Committee (Cork, Ireland) for the recruitment and sampling of these Cohorts of children.

Recruiting of the children affected by S-ECC was performed at a Hospital Dental Treatment Centre. All children were referred to the Hospital Dental Treatment Centre (Cork, Ireland) where clinical examinations were performed and they were then scheduled for extraction of their carious teeth, under general anaesthetic. All these children referred to the Hospital Dental Treatment Centre for general anaesthetic had radiographs taken as part of the examination prior to referral. The caries-free cohort was recruited from various crèches, where a paediatric dentist travelled to each of the crèches and examined the children’s teeth.

Inclusion criteria applied to both caries-free and caries-affected groups were that they were medically healthy, had no antibiotic intake in the 3 months prior to sampling, and were under the age of 60 months.

Diet and lifestyle data collection

Before undergoing the dental examination, informed consent was obtained from the parent/guardians of the children. Habitual dietary data was also collected using a validated Food Frequency Questionnaire (FFQ) [27] which was provided to each parent/guardian, to record food intake (Additional file 1). A detailed questionnaire was given to each parent/guardian. Data collected included antibiotic treatment history, general medical history and dental history, fluoridation status of home (well/public water), feeding practices in infancy, birth mode, and oral health related quality of life.

Oral examination

Oral examinations for both caries-free and caries-active groups were performed by a trained Paediatric dentist after a full dental examination was completed by a dental surgeon in the clinic. Children in crèches and in the Hospital Dental Treatment Centre were examined in a quiet area with their parents present. Children’s teeth were wiped with a cotton wool roll and sterile gauze square to remove plaque and debris prior to examination which was carried out under natural light, using a standard size 4 mirror and ball ended CPI ‘C’ probe [60]. The mouth was illuminated with a Promed Penlight, which consists of bright concentrated halogen light when natural light was insufficient. For the caries-active group, caries was recorded at the level of cavitation into dentine (cavitation level), using the WHO criteria [60, 61], with the addition of visible non cavitated dentine caries as referenced by Whelton et al. [62]. The International Caries and Detection Assessment System (ICDAS) code for the caries affected teeth were within codes 5 and 6 [63]. The dmft score was recorded along with the dmfs score, and sample collection performed. For the caries-free group, caries was recorded at the level of cavitation into dentine (cavitation level), using the WHO criteria [61]. The dmft/dmfs was measured, and sample collection was performed. Caries-free children did not show clinical evidence of early pre-cavitation of caries or white spot lesions and had no history of treatment on any tooth surfaces, as defined [64].

Sampling

The same-trained paediatric dentist took all samples after the teeth were examined. All children were instructed not to brush their teeth the evening and the morning before sampling. A CatchAll™ collection swab, with hard pack for storage after collection was used (Cambio UK) [65] (See image: Additional file 2). For the caries-active S-ECC group, both a carious lesion sample and a saliva sample were taken. After a full dental examination, and pre-general anaesthetic, the CatchAll™ collection swab was used to collect pooled unstimulated saliva in the floor of the mouth for 1–2 min. This unstimulated saliva sample is recognized as a representation of the whole oral ecosystem [33, 66, 67]. The swab was placed back in the collection tube, and stored at − 80 °C. To sample the carious lesions, while the child was under general anaesthetic, the carious deciduous tooth was extracted and under isolation, the tooth was irrigated with saline. Under care, by a paediatric dentist, the superficial carious dentin was excavated with a sterile spoon excavator and the next layer of deep dentinal caries was excavated using a new separate sterile spoon excavator and the sample was pooled in a sterile 1.5-ml micro-centrifuge tube with 1 ml of TE buffer (50 Mm Tris-HCL, 1 Mm EDTA). The samples were placed in a sterile 1.5-ml micro-centrifuge tube and transported to the laboratory, where they were frozen until further analysis and stored at − 80 °C. For the caries-free group, after a full dental examination, the CatchAll™ collection swab was used to collect pooled unstimulated saliva in the floor of the mouth for 1–2 min. The swab was placed back in the collection tube, and stored at − 80 °C.

DNA extraction

Extraction of DNA from all samples was carried out with the MO BIO PowerLyzer® 24 homogenizer following some initial optimisation for extraction from an oral catch-all swab rather than a soil sample as previously described [68]. The saliva sample was contained in a catch-all swab in the end of a collection tube. The tube was cut 1 cm above this swab, and this was inserted into the PowerBead tubes, to which 60 μl of solution C1 had been added. Tubes were incubated at 65 °C for 10 min and then shaken horizontally at maximum speed for 2 min, using the MO BIO vortex adapter. The remainder of the protocol was followed as per manufacturer’s instructions. For the caries sample, the tubes were incubated at 65 °C for 10 min and then shaken horizontally at maximum speed for 4 min, using the MO BIO vortex adapter. The remainder of the protocol was followed as per manufacturer’s instructions. DNA was visualised on a 0.8% agarose gel and quantified using the Nanodrop 1000 (Thermo Scientific, Ireland). DNA was then stored at − 80 °C.

16 s rRNA gene amplification primers

Primers used for PCR amplification were the V4 - V5 region primers 520F (AYTGGGYDTAAAGNG) and 926R (CCGTCAATTYYTTTRAGTTT). Initial primers for Illumina sequencing contain the sequencing primer binding sites, forward or reverse 16S rRNA gene specific primer, and a 10 nt in-line multiplexing identifier (MID). Dual separate MIDs were attached to both ends of the PCR product.

The V4 - V5 amplicons for Illumina sequencing were generated using a two-step amplification procedure. The first step reaction mix contained 50 μl BIO-X-ACT™ Short Mix (BIOLINE), 10 μl of 2 nM forward and reverse primers, 50 ng genomic DNA, and ddH20 to give a final volume of 100 μl. Cycling conditions were: an initial 95 °C, 5-min denaturation step; 30 cycles of 95 °C for 15 s, 42 °C for 15 s, and 72 °C for 30s; and a final 10-min extension at 72 °C. The products were purified using SPRIselect beads (Beckman Coulter, Indianapolis IN) as per manufacturer’s instructions, using a 0.9:1 volume ratio of beads to product. The purified PCR products were eluted in 40 μl of ddH20. DNA quantity was assessed via Quant-iT™ PicoGreen® dsDNA Assay Kit (Invitrogen™). The samples were pooled in equimolar amounts and then sequenced by Eurofins Genomics (Eurofins Genetic Services Ltd., I54 Business Park, Valiant way Wolverhampton WV9 5GB, UK) using Illumina MiSeq 2 × 300 bp paired end technology. Nextflex Rapid library preparation was carried out by the company to attach bridge adaptors necessary for clustering. Sequencing of 16S DNA was carried out on the V4/V5 region using a Miseq (301 bp paired-end reads). Sequence data were stored on a Linux server and backed up on external hard-drives.

Bioinformatic analysis

Sequence processing, OTU clustering and taxonomy assignment

The software, flash (v1.2.8), was used to join paired-end reads. Paired-end reads with more than 25% incorrect bases in their region of overlap were excluded from subsequent steps. Qiime (v1.9.1) was used to extract barcodes (extract_barcodes.py) and for demultiplexing (split_libraries_fastq.py).

The USEARCH (v8.0.1623) pipeline was used for the following steps: de-replication of reads (identical reads are represented by a single sequence), exclusion of reads shorter than 350 bp and longer than 370 bp, exclusion of unique reads, chimera filtering, OTU clustering at 97% identity and calculation of representative OTU sequences. Using USEARCH, all reads (including unique reads) were then mapped back to the representative OTU sequences to give the final OTU read count for each sample. The software fastQC (v0.11.3) was used after each filtering step to assess read quality. The median read count for the samples after sequence processing was 44,400. The sample number after sequence processing was 206.

Part of the mothur (v1.36.1) [69] pipeline was used to run the RDP classifier using a filtered version of the RDP database in order to assign taxonomy down to genus level. The software SPINGO (v1.3) [70], was used to assign taxonomy at species level. For both mothur/RDP and SPINGO, confidence cut-offs of 80% were used.

Alpha and beta diversity analysis

Alpha and beta diversity metrics were calculated in Qiime (v1.9.1) [71]. To calculate diversity metrics, several additional steps were carried out (also in Qiime). The OTU table was rarefied (single_rarefaction.py) at 10,540 reads (the lowest read count in the dataset). Representative OTU sequences were aligned using pyNAST (align_seqs.py) and filtered to remove columns that do not contribute to phylogenetic signal (filter_alignment.py). A phylogenetic tree was generated using FastTree (make_phylogeny.py). This tree is necessary for phylogenetic alpha and beta diversity metrics. The rarefied OTU table was used in the calculation of all diversity metrics.

The following alpha diversity metrics were calculated: chao1, Shannon (Shannon’s index), Simpson (Simpson’s index), Observed species (OTU count) and Phylogenetic (PD whole tree). The following beta diversity metrics were calculated: weighted and un-weighted unifrac distances, and Bray-Curtis dissimilarity.

Statistics and data visualisation

All statistics and data visualisation were carried out in R (v3.2.3) [72]. Paired Mann-Whitney tests were used to compare microbiota of saliva samples (CAS) (n = 68) with that of caries samples (CAC) (n = 68) taken from the caries group. Both CAC and CAS are paired since both samples from the same subject. Un-paired Mann-Whitney tests were used to compare the saliva and caries samples from the caries group with a control group of caries-free saliva (CFS) (n = 70) of caries-free individuals. Benjamini and Hochberg correction [73] was used to adjust p-values for multiple testing. Level of significance was set at p < 0.05.

Results

A distinct microbiota in caries lesions but not saliva in children with S-ECC

Given that S-ECC is such an acute disease, it seemed possible that it resulted from a global microbiota change in the oral cavity. To investigate the relatedness of microbiome composition in CAC (caries-active caries), CAS (caries-active saliva) and CFS (caries-free saliva), we generated PCoA (principle co-ordinates) plots showing relatedness by two established metrics, Bray Curtis dissimilarity, and UniFrac distances. The Bray Curtis plot (Fig. 1a) shows separation of the caries lesion samples (CAC) and the two saliva sample types (CFS & CAS) based on PCoA axes 1 and 2. The microbiota of the two saliva groups (CAS & CFS) considerably overlap, indicating a similarity in the general composition of microbial taxa. For weighted UniFrac, (Fig. 1b) the caries microbiota group is again separated from the two saliva microbiota groups, with only minimal overlap of some samples. Combined with the Bray Curtis analysis, this shows convincingly that there is no major separation between the microbiota of CAS and CFS groups, even though this metric is very sensitive to the differences in the presence/absence and abundance of OTUs/samples.

Fig. 1
figure 1

PCoA (principle co-ordinates) plots showing relatedness by two established metrics, Bray Curtis dissimilarity, and UniFrac distances, while unweighted UniFrac illustrates separation between the three groups (CAC, CAS & CFS). a. Plot of principle co-ordinates using Bray-Curtis dissimilarity. Points are coloured according to group and ellipses describe the distribution of points for each group. Percentage variation explained: PCA 1 (22.3%) and PCA 2 (7.7%). b. Plot of principle co-ordinates using weighted unifrac distance. Points are coloured according to group and ellipses describe the distribution of points for each group. Percentage variation explained: PCA 1 (45%) and PCA 2 (11.7%). c. Plot of principle co-ordinates using un-weighted unifrac distance. Points are coloured according to group and ellipses describe the distribution of points for each group. Percentage variation explained: PCA 1 (16.6%) and PCA 2 (4.9%)

Plotting the second UniFrac metric, unweighted UniFrac distances (Fig. 1c), illustrates separation between the three groups (CAC, CAS & CFS). This index measures the presence and absence of taxa only and does not adjust the distance metric according to taxon abundance, so unlike weighted analysis, it reflects the contribution of rare taxa (that get overwhelmed in a weighted analysis). When presence/absence of taxa is the dominant parameter used to calculate distance, separation between all three groups occurs. The CAS microbiota was closer to the CAC microbiota than the CFS microbiota, suggesting rare taxa are shared between the former two samples.

Microbiota diversity in caries lesions is lower than that of saliva from children with or without S-ECC

To study the diversity of the microbiota from the caries lesions (CAC) and saliva samples of both caries-free children (CFS) and caries-active children (CAS), a series of alpha diversity metrics was determined: the Chao index, phylogenetic diversity (PD whole tree), observed species (OTU count), the Simpson index and the Shannon index (Fig. 2). All metrics illustrate that the diversity of caries microbiota was the lowest of these sample types. The difference in the diversity values of the caries samples and saliva samples was significant (p < 0.001). The Chao diversity measurement, as illustrated in Fig. 2a, estimates the number of species from observed data, and the diversity of the low abundance taxa, and according to this metric, the CFS versus CAS comparison was significantly different with a p-value of < 0.05. All other alpha diversity comparisons (Fig. 2b, c, d) indicated that the diversity difference between the caries samples and saliva samples was significant (p < 0.001).

Fig. 2
figure 2

Pairwise alpha diversity comparisons of saliva and caries microbiota. a. Boxplot of chao1 diversity in the three groups. Outliers are represented by black points. Significant differences between groups are shown by arrows and the following notation: p < 0.05 (*), p < 0.01 (**) and p < 0.001 (***). b. Boxplot of observed species in the three groups. Outliers are represented by black points. Significant differences between groups are shown by arrows and the following notation: p < 0.05 (*), p < 0.01 (**) and p < 0.001 (***). c. Boxplot of Phylogenetic diversity in the three groups. Outliers are represented by black points. Significant differences between groups are shown by arrows and the following notation: p < 0.05 (*), p < 0.01 (**) and p < 0.001 (***). d. Boxplot of Shannon diversity in the three groups. Outliers are represented by black points. Significant differences between groups are shown by arrows and the following notation: p < 0.05 (*), p < 0.01 (**) and p < 0.001 (***)

Habitual diet is not significantly different in children with or without S-ECC

Diet can have a profound impact on oral health and caries risk. We found minor technical inconsistencies within the FFQ recording of data between patients and this could have a negative impact on the significance of diet and its frequency on the oral microbiota of these children. With the data that was recorded (Additional file 3), no significant differences were found when the habitual intake of each food group was compared from the FFQ data derived from the caries-active and caries-free subjects. No food group was consumed at significantly different frequency when we tested either for unequal presence/absence of food groups in the diet using the Fisher’s test, or different frequencies of food groups in the diet using the Mann-Whitney test. Any trends were not supported by significant p-values (< 0.05).

Differentially abundant taxa in CAC compared to saliva microbiota of both CAS & CFS

Differential bacterial taxon abundance in compared microbiota datasets may be graphically demonstrated by hierarchical clustering, whereby samples are grouped based on similarity of the taxa in their microbiota. These relatedness levels between samples, and their constituent microbial taxa, are represented by vertical and horizontal dendrograms incorporating a “heatmap” colour scale to convey abundance levels. Figure 3 illustrates the abundance of bacterial taxa at family level. The three groups are split into two main branches visible in the horizontal dendrogram above the colour bar in Fig. 3. CAC clusters on one branch (red bar) and the two saliva groups (CAS in green & CFS in blue) cluster on the other. There is clustering of the CAS and the CFS within this branch, showing considerable similarity between the two groups at family level. When compared against the patient metadata, there was no obvious variable that convincingly separates CAS from CFS (data not shown).

Fig. 3
figure 3

Hierarchical clustering of microbiota data at bacterial family level. Abundances are colour-coded according to the colour key on the top left with grey representing a value of zero. Euclidean distance and complete linkage were used to cluster the rows and columns of the heatmap. The colour bar on top of the heatmap corresponds to sample type: CAC red, CAS green and CFS blue. All taxa present at less than 1% in all three groups are excluded from the heatmap

The microbiota of the CAC samples was characterized by high relative abundance of Prevotellaceae, Veillonellaceae, Bifidobacteriacae and Streptococcaceae, and by low relative abundance of Corynebacteriaceae, Carnobacteriaceae, Aerococcaceae, and Micrococcaceae. Both saliva sample types (CAS & CFS) showed higher abundances than caries samples of Leptotrichiaceae, Porphyromonadaceae and Flavobacteriaceae and of Neisseriaceae and Pasteurellaceae, while illustrating a very low abundance compared to CAC of Spirochaetaceae, Bifidobacteriaceae and Lactobacillaceae.

A more nuanced picture emerges when differentially abundant taxa were analysed at genus level (Fig. 4). The samples again separate laterally into caries versus saliva (with both caries-active and caries-free clustering together). The microbiota content appears to be split vertically in two groups of differentially abundant genera as revealed by the dendogram on the Y-axis based on bacterial abundance. The top branch is divided, with CAS and CAC showing higher abundance of Streptococcus and Prevotella, and with Neisseria at a higher abundance in both. In the saliva samples (CAS & CFS), Leptotrichia, Porphyromonas and Haemophilus are in higher abundance, with Leptotrichia at higher abundance in CFS than CAS.

Fig. 4
figure 4

Hierarchical clustering of microbiota data at bacterial genus level. Abundances are colour-coded according to the colour key on the top left with grey representing a value of zero. Euclidean distance and complete linkage were used to cluster the rows and columns of the heatmap. The colour bar on top of the heatmap is coloured according to sample type: CAC red, CAS green and CFS blue. All taxa present at less than 1% in all three groups are excluded from the heatmap

The lower branch shows clear low-abundance of the genera Lactobacillus, Treponema, Scardovia, and Parascardovia in the CFS & CAS. (Fig. 4), while low abundance taxa in CAC are Gemella and Granulicatella.

When samples in the analysis with at least one species with a median value of ≥0.5% were included, fewer outliers were identified, with only one CAS sample within the CAC branch, and two CAC within the second branch of the saliva samples. There was clearer separation of the samples, with CAC on the first branch, with high abundance of Streptococcus mutans, compared to the saliva samples. The two saliva samples split into two branch points, with CFS illustrating clustering within the first branch, with higher abundance of Tannerella forsythia, Capnocytophaga gingivalis, and Leptotrichia buccalis. Species with a clear low abundance in the majority of saliva samples included Scardovia wiggsiae, Parascardovia denticolens, Prevotella denticola and Prevotella oris, where these were present at higher abundance in CAC.

Broad and fine detail compositional differences distinguish caries microbiota from paired and healthy control saliva samples

At phylum level (Fig. 5a), the CAC microbiota was dominated by Firmicutes (median abundance value 33.45%), while the CAS microbiota was dominated by Proteobacteria (median abundance value 38.18%; p < 0.0001). At phylum level, both the CAS and CFS microbiota composition is quite similar (Fig. 5a), with CAS and CFS microbiota dominated by Proteobacteria at 38.2 and 36.1% respectively. The main phylum difference between saliva of CAS versus CFS children was Fusobacteria. Its abundance has a median value of 13.4% in the saliva microbiota of CFS children, but 7.37% in CAS (p < 0.0001).

Fig. 5
figure 5

Broad and fine detail compositional differences at Genus, phylum and species level. a. Microbiota composition at phylum level. Percentages for each taxon represent the median abundance values for the sample types. b. Barplot of percentage abundance at genus level. Percentages for each taxon represent the median values for the groups. c. Barplot of percentage abundance at species level. Percentages for each taxon represent the median values for the groups

Fig. 6
figure 6

Hierarchical clustering of microbiota data at bacterial species level. Abundances are colour-coded according to the colour key on the top left with grey representing a value of zero. Euclidean distance and complete linkage were used to cluster the rows and columns of the heatmap. The colour bar on top of the heatmap is coloured according to sample type: CAC red, CAS green and CFS blue. All taxa present with at least one species with a median value ≥0.5% in all three groups are included

Excluding the “other” and “unassigned” categories, the three microbiota types (CAC, CFS & CAS) are dominated by taxa Neisseria, Capnocytophaga, Porphyromonas, Streptococcus, Prevotella, Leptotrichia, and Haemophilus. Streptococcus, Neisseria, Prevotella, Capnocytophaga, dominate the CAC sample microbiota and at lower levels, taxa Scardovia, Parascardovia, Selenomonas and Lactobacillus. CAC sample microbiota include numerous species of higher relative abundance: Streptococcus mutans, Alloprevotella denticola, Prevotella histicola, Scardovia wiggsiae, Parascardovia denticolens, Prevotella tannerae and Bifidobacterium dentium.

The CFS microbiota composition differs significantly from the saliva of CAS group by the presence of the following microbial genera: Leptotrichia, Bifidobacterium, Corynebacterium, Alloprevotella, Cardiobacterium and Veillonella (p < 0.0001). The abundance of all six genera was significantly higher in CFS (p < 0.0001).

The CAS microbiota was dominated by Neisseria, Porphyromonas, Streptococcus and Haemophilus, and species included Streptococcus mutans, Prevotella histicola, Prevotella melaninogenica, Porphyromonas catoniae and Prevotella salivae. CFS samples were dominated by Leptotrichia, Capnocytophaga, Neisseria, Haemophilus, Streptococcus and Porphyromonas at genus level, while at species level, the CFS group include species Haemophilus haemolyticus, Haemophilus parainfluenzae, Rothia mucilaginosa, Porphyromonas catoniae and Streptococcus sanguinis.

Microbial taxa showing statistically significant differential abundance between CFS and CAS children included Streptococcus mutans, Haemophilus parainfluenzae, Prevotella histicola, Leptotrichia buccalis, Veillonella dispar, Alloprevotella tannerae and Prevotella salivae. Interestingly, Streptococcus mutans, Prevotella histicola and Veillonella dispar were present at higher abundance in CAS than the saliva of CFS.

Discussion

In this study, next generation sequencing of 16S amplicons was used to explore the microbiota of deep dentinal carious lesions and saliva of Irish children affected with S-ECC and the salivary microbiota of those that are caries-free.

Using two well established metrics (Bray Curtis dissimilarity and UniFrac distance) the caries dentine microbiota was found to be distinct from that of either CFS or CAS, illustrating, when abundance is taken into account, that CAC has considerably different proportions of certain high-abundance taxa. However, when measured using unweighted UniFrac, which measures presence and absence of taxa only, all three groups were found to be distinct, showing that each group is different in terms of rare or low-abundance taxa. Of interest was the closer microbiota relatedness of CAS to CAC, with the presence of certain CAS outliers overlapping with CAC. This suggests that some samples share similar low-abundance taxa between CAC and CAS or that CAS harbours some taxa which may have stimulated the increased caries rate compared to CFS, or there is some shedding of taxa from CAC into the CAS. When a leave-one-out strategy was implemented with dominant taxa such as Neisseria and Streptococcus, these outliers changed position, sometimes clustering within their own group and a few samples that clustered within their group in the full dataset became outliers with the reduced dataset. When these outliers were studied in detail, these samples were found not to be influenced by any metadata, suggesting that factors not accounted for in this study are responsible for unusual taxon composition in a subset of samples.

Focusing on the dominant taxa and their abundance illustrated graphically by hierarchical clustering on the heatmap, samples were identified and grouped based on microbiota similarity, but also of interest was the partial agreement between the outliers identifiable in the PCoA plots and in the heatmaps. The identity of these outliers could not be readily correlated with patient metadata and when the taxon abundance level was investigated in detail, ar genus and family levels ther is splitting of CAC with both the saliva groups (CFS & CAS), with the CFS samples clustering mostly together, with some intermixing of both saliva microbiota types. There is clear clustering of CAC in branch one to left, with CFS at next branch in blue (Figs. 3, 4), with CAS lastly split on a third branch. This is more apparent at species level by hierarchical clustering, with less intermixing of both salivary microbiota types (Fig. 6).

This incomplete separation of samples between saliva groups suggests that the salivary microbiota is not specific enough to be used as an identifier for caries risk in children. The oral cavity is an entry point for colonisation of microbial species and saliva is a reservoir for a multitude of bacteria, with its microbial and nutritional composition being shaped by food intake, reflux, environment and other influences [38, 67, 74]. In this study, no association was found between factors like habitual diet, brushing habits or fluid intake with microbiota composition, although it is possible these metadata are not sufficiently granular in the FFQ data and patient questionnaire. Furthermore, previous studies have shown variability between sites in the oral cavity itself, with niches among the tongue, soft and hard palates, supra- and sub-gingival surfaces of teeth and saliva each demonstrating microbiota variability [32, 66]. The flow rate, buffering capacity, and molecules within saliva which can aid attachment of bacterial cells, all play a role in both the compositional balance of the oral microbiome [75]. While some studies have also found an association between the microbiota and disease in plaque samples, but not within the saliva samples [37, 76] our findings support these findings, with saliva and caries representing two distinct habitats.

As caries lesions progress and become more severe, the diversity of caries microbiota decreases [34, 35, 47]. In this study, alpha diversity of caries microbiota was lowest, and differences in the diversity of the caries samples and saliva samples were significant (p < 0.001). However the Chao diversity index for CFS was significantly higher than CAS with a p-value of < 0.05, again supporting the previous data, that when low-abundance or rare taxa are given an equal weighting to higher-abundance taxa, differences in diversity between the two saliva groups become apparent. This suggests that even at low abundance, certain taxa such as the acid-producing lactobacilli can play a strong role in caries progression. This highlights an important possibility, that taxon abundance of cavity-causing microbes may not be strongly correlated with progression of caries; low-abundance taxa at abundance levels that typically do not feature in microbiome summary data might be the main indicator of future tooth decay because, for instance, a small number of acid-producing or biofilm-producing species may have a disproportionate impact on oral health. In addition, the potential affect that the removal or exclusion of these rare taxa may have on caries prevention and general oral health makes a solid case for their identification.

As caries progresses to a more advanced state, the bacteria that dominate this cavity are less diverse, because aciduric organisms have been selected and enriched, and we found that the main genera that dominated the CAC lesion were Neisseria, Streptococcus and Prevotella, while the species that dominate the caries lesion (CAC) include Streptococcus mutans, Prevotella sp, Scardovia sp. and Bifidobacterium dentium. Neisseria, Streptococcus, Prevotella and Porphyromonas have all been strongly associated with caries in past studies [29, 35, 38, 40, 41, 44, 58]. Neisseria and Streptococcus produce acid which lowers the pH of the mouth and leads to increased demineralisation of enamel [47] while Prevotella has a known role in caries progression and endodontic infections [7, 50]. The high level of Streptococcus mutans in the carious lesion is consistent with previous studies, and its presence is a strong indicator for caries [58]. Streptococcus mutans aids in caries initiation by adhering to the enamel, forming a cariogenic biofilm via glucan synthesis on the surface, aiding the binding of other species. While being aciduric and acidogenic, it can often be part of a more complex community of microorganisms working together [40, 44, 47, 77], and is a risk factor for caries progression [34].

Scardovia is documented as having a role as a cariogenic bacterium involved in the later stages of S-ECC [40]. Scardovia wiggsiae is significantly associated with S-ECC, based on a culture study of plaque from children [36], and in adults with caries [78]. Prevotella species have been shown to play an important role in endodontic infections [50], and Prevotella tannerae, Prevotella histicola (isolated from human oral mucosa [37]) and Alloprevotella denticola [41, 74] have all been shown to be associated with dental caries. Lactobacillus, which is notably associated with caries progression [22, 39, 44, 48, 79] was found at very low levels compared to other genera (0.675% (CAC) and 0.031% (CAS)). These low levels were also reported in previous studies [34, 38, 47] and this interesting finding supports the idea that when certain acid producers are at low levels, their acidogenic properties may nevertheless be strong enough to allow other acid producers to take their place, such as Neisseria, Selenomonas and Streptococcus mitis [47]. Neisseria spp. have the ability to metabolize glucose to produce lactic acid and this genus was found at high levels in CAC (10.29%) and in CAS (16.28%) and may have an active cariogenic role. Lactobacillus is found at low levels in endodontic infections with deep caries [80, 81]. It has also been suggested by Rôças et al. that altered Lactobacillus abundance may be due to the change from cariogenic microbiota to a microbiota that stimulates progression into pulpal tissue causing infection. Shifts in the microbiota composition at the outermost pulpal layer can be affected by saliva and diet to the outermost pulpal layer, while the inner deeper layer has a different environment in comparison [57]. The replacement of Lactobacillus with other taxa could be linked to the degree of pain, duration of pain, length of caries destruction, connection with pulp, diet and environmental factors, and this fine detail could reveal reasons for the low levels of Lactobacillus in these teeth.

Within CFS there were higher levels of Capnocytophaga (10.9%) and Leptotrichia (8.1% CFS), and lower levels of Porphyromonas and Neisseria. There is evidence to suggest that Capnocytophaga and Leptotrichia are health-associated species [36, 44, 58, 77] and Capnocytophaga has been found at higher levels in caries-free subjects [37, 47, 58]. Lower levels of Porphyromonas were found (9.5% CAS & 6.6% CFS) and Neisseria (16.3% CAS & 12.6% CFS) in CFS compared to CAS. Some Neisseria sp. have been shown to play a role in acid production (e.g. N. gonorrhoeae and N. meningitidis) [47], while other Neisseria species such as Neisseria flavescens have been shown to elicit higher signal of probes when targeted in caries-free children [67] together with the Porphyromonas gram-negative species, Porphyromonas catoniae. The association of these species with a caries-free oral status has been further supported by Nyvad et al. [29]. A higher abundance of Porphyromonas catoniae was found within the CAS (9.5%). Studies have reported children with high levels of caries harbouring higher levels of Porphyromonas in their saliva [78], by culture study [40], and it has been detected in root canals of necrotic deciduous teeth [79] and our findings demonstrate a positive association of Porphyromonas with S-ECC.

Conclusion

In conclusion, we identified distinct differences between the caries microbiota and saliva microbiota, with separation of both salivary groups (caries-active and caries-free), showing a clear separation when low abundance taxa were considered. While the microbiota diversity in the caries dentinal microbiota was lower than both salivary microbiota types, there were clear compositional differences between all groups from phylum to species. Firmicutes dominated the CAC, while Proteobacteria dominated the CAS and CFS salivary microbiota, and genera that dominated the CAC lesion were Neisseria, Streptococcus and Prevotella. The taxa present in the dentinal lesions could be potential instigators that drive migration of infection to the pulp, while the saliva microbiota in health and disease may be associated with caries-active or caries-free status in these children.

Abbreviations

BLAST:

Basic local alignment search tool

CAC:

Caries-active caries

CAS:

Caries-active saliva

CFS:

Caries-free saliva

dmfs:

Decayed, missing (due to caries), or filled tooth surfaces in primary teeth

dmft:

Decayed, missing (due to caries), or filled tooth in primary teeth

ECC:

Early childhood caries

FFQ:

Food frequency questionnaire

FLASH:

Fast Length Adjustment of SHort reads to improve genome assemblies

MID:

Multiplexing identifier

OTU:

Operational taxonomic unit

PCoA plot:

Principle co-ordinates plot

QIIME:

Quantitative Insights Into Microbial Ecology

rRNA:

ribosomal Ribonucleic acid

S- ECC:

Severe Early childhood caries

References

  1. National Call to Action to Promote Oral Health. U.S. Department of Health and Human Services, Public Health Service NI of HNI of D and CR. Rockville: U.S. Department of Health and Human Services; 2003.

    Google Scholar 

  2. World Health Organisation. Oral disease burdens and common risk factors: WHO; 2012. http://www.who.int/oral_health/disease_burden/global/en/. Accessed 6 Mar 2017

  3. American Academy of Pediatric Dentistry(AAPD). Policy on early childhood caries (ECC): classifications, consequences, and preventive strategies. Pediatr Dent. 2008;30(7 Suppl):40–3.

    Google Scholar 

  4. Gussy MG, Waters EG, Walsh O, Kilpatrick NM. Early childhood caries: current evidence for aetiology and prevention. J Paediatr Child Health. 2006;42:37–43.

    Article  PubMed  Google Scholar 

  5. Finucane D. Rationale for restoration of carious primary teeth: a review. Eur Arch Paediatr Dent. 2012;13:281–92.

    Article  PubMed  Google Scholar 

  6. Martins-Júnior PA, Vieira-Andrade RG, Corrêa-Faria P, Oliveira-Ferreira F, Marques LS, Ramos-Jorge ML. Impact of early childhood caries on the oral health-related quality of life of preschool children and their parents. Caries Res. 2013;47:211–8.

    Article  PubMed  Google Scholar 

  7. Tanner ACR, Kent RL, Holgerson PL, Hughes CV, Loo CY, Kanasi E, et al. Microbiota of severe early childhood caries before and after therapy. J Dent Res. 2011;90:1298–305.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Almeida AG, Roseman MM, Sheff M, Huntington N, Hughes CV. Future caries susceptibility in children with early childhood caries following treatment under general anesthesia. Pediatr Dent. 2000;22:302–6.

    PubMed  Google Scholar 

  9. Graves CE, Berkowitz RJ, Proskin HM, Chase I, Weinstein P, Billings R. Clinical outcomes for early childhood caries: influence of aggressive dental surgery. J Dent Child (Chic). 2004;71:114–7.

    Google Scholar 

  10. Albadri SS, Lee S, Lee GT, Llewelyn R, Blinkhorn AS, Mackie IC. The use of general anaesthesia for the extraction of children’s teeth. Results from two UK dental hospitals. Eur Arch Paediatr Dent. 2006;7:110–5.

    PubMed  Google Scholar 

  11. McAuliffe U, Kinirons M, Woods NHM. A retrospective investigation of the oral health records of a cohort of prechool children who received extractions under general anaesthesis including cost analysis of treatment. J Irish Dent Assoc. 2017;63:38–44.

    Google Scholar 

  12. Li Y, Wang W. Predicting caries in permanent teeth from caries in primary teeth: an eight-year cohort study. J Dent Res. 2002;81:561–6.

    Article  PubMed  Google Scholar 

  13. Alm A, Wendt LK, Koch G, Birkhed D. Oral hygiene and parent-related factors during early childhood in relation to approximal caries at 15 years of age. Caries Res. 2008;42:28–36.

    Article  PubMed  Google Scholar 

  14. Dülgergil CT, Colak H. Do the more caries in early primary dentition indicate the more caries in permanent dentition? Results of a 5-years follow-up study in rural-district. J Int Soc Prev Community Dent. 2012;2:48–52.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Berkowitz RJ, Amante A, Kopycka-Kedzierawski DT, Billings RJ, Feng C. Dental caries recurrence following clinical treatment for severe early childhood caries. Pediatr Dent. 2011;33:510–4.

    PubMed  Google Scholar 

  16. Ng MW, Ramos-Gomez F, Lieberman M, Lee JY, Scoville R, Hannon C, et al. Disease Management of Early Childhood Caries: ECC collaborative project. Int J Dent. 2014;2014:327801.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Van Houte J. Role of micro-organisms in caries etiology. J Dent Res. 1994;73:672–81.

    Article  PubMed  Google Scholar 

  18. Borgström MK, Sullivan A, Granath L, Nilsson G. On the pH-lowering potential of lactobacilli and mutans streptococci from dental plaque related to the prevalence of caries. Community Dent Oral Epidemiol. 1997;25:165–9.

    Article  PubMed  Google Scholar 

  19. Hughes CV, Dahlan M, Papadopolou E, Loo CY, Pradhan NS, Lu SC, et al. Aciduric microbiota and mutans streptococci in severe and recurrent severe early childhood caries. Pediatr Dent. 2012;34:e16–23.

    PubMed  Google Scholar 

  20. Pearce C, Bowden GH, Evans M, Fitzsimmons SP, Johnson J, Sheridan MJ, et al. Identification of pioneer viridans streptococci in the oral cavity of human neonates. J Med Microbiol. 1995;42:67–72.

    Article  PubMed  Google Scholar 

  21. Beighton D. The complex oral microflora of high-risk individuals and groups and its role in the caries process. Community Dent Oral Epidemiol. 2005;33:248–55.

    Article  PubMed  Google Scholar 

  22. Yang R, Argimon S, Li Y, Gu H, Zhou X, Caufield PW. Determining the genetic diversity of lactobacilli from the oral cavity. J Microbiol Methods. 2010;82:163–9.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Rosebury T. Acid production and tolerance of lactobacilli from dental caries and other sources as measured by the glass electrode. J Bacteriol. 1932;24:321–34.

    PubMed  PubMed Central  Google Scholar 

  24. Warren JJ, Weber-Gasparoni K, Marshall TA, Drake DR, Dehkordi-Vakil F, Dawson DV, et al. A longitudinal study of dental caries risk among very young low SES children. Community Dent Oral Epidemiol. 2009;37(2):116–22.

    Article  PubMed  Google Scholar 

  25. Mangla RG, Kapur R, Dhindsa A, Madan M. Prevalence and associated risk factors of severe early childhood caries in 12- to 36-month-old children of Sirmaur District, Himachal Pradesh. India Int J Clin Pediatr Dent. 2017;10:183–7. https://doi.org/10.5005/jp-journals-10005-1431.

    Article  Google Scholar 

  26. Wade WG. The oral microbiome in health and disease. Pharmacol Res. 2013;69:137–43. https://doi.org/10.1016/j.phrs.2012.11.006.

    Article  PubMed  Google Scholar 

  27. Claesson MJ, Jeffery IB, Conde S, Power SE, O’Connor EM, Cusack S, et al. Gut microbiota composition correlates with diet and health in the elderly. Nature. 2012;488:178–84.

    Article  PubMed  Google Scholar 

  28. Conlon MA, Bird AR. The impact of diet and lifestyle on gut microbiota and human health. Nutrients. 2015;7:17–44. https://doi.org/10.3390/nu7010017.

    Article  Google Scholar 

  29. Nyvad B, Crielaard W, Mira A, Takahashi N, Beighton D. Dental caries from a molecular microbiological perspective. Caries Res. 2013;47:89–102.

    Article  PubMed  Google Scholar 

  30. Paster BJ, Olsen I, Aas JA, Dewhirst FE. The breadth of bacterial diversity in the human periodontal pocket and other oral sites. Periodontol. 2006;42:80–7.

    Article  Google Scholar 

  31. Paster BJ, Dewhirst FE. Molecular microbial diagnosis. Periodontol. 2009;51:38–44.

    Article  Google Scholar 

  32. Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner ACR, Yu WH, et al. The human oral microbiome. J Bacteriol. 2010;192:5002–17.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Zaura E, Keijser BJF, Huse SM, Crielaard W. Defining the healthy “core microbiome” of oral microbial communities. BMC Microbiol. 2009;9:259.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Gross EEL, Beall CJC, Kutsch SR, Firestone NDN, Leys EJE, Griffen AAL, et al. Beyond Streptococcus mutans: dental caries onset linked to multiple species by 16S rRNA community analysis. PLoS One. 2012;7:1–11.

    Google Scholar 

  35. Li Y, Ge Y, Saxena D, Caufield PW. Genetic profiling of the oral microbiota associated with severe early-childhood caries. J Clin Microbiol. 2007;45:81–7.

    Article  PubMed  Google Scholar 

  36. Xu H, Hao W, Zhou Q, Wang W, Xia Z, Liu C, et al. Plaque bacterial microbiome diversity in children younger than 30 months with or without caries prior to eruption of second primary molars. PLoS One. 2014;9:e89269.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Ling Z, Kong J, Jia P, Wei C, Wang Y, Pan Z, et al. Analysis of oral microbiota in children with dental caries by PCR-DGGE and barcoded pyrosequencing. Microb Ecol. 2010;60:677–90.

    Article  PubMed  Google Scholar 

  38. Nasidze I, Li J, Quinque D, Tang K, Stoneking M. Global diversity in the human salivary microbiome. Genome Res. 2009;19:636–43.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Obata J, Takeshita T, Shibata Y, Yamanaka W, Unemori M, Akamine A, et al. Identification of the microbiota in carious dentin lesions using 16S rRNA gene sequencing. PLoS One. 2014;9:e103712.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Tanner ACR, Mathney JMJ, Kent RL, Chalmers NI, Hughes CV, Loo CY, et al. Cultivable anaerobic microbiota of severe early childhood caries. J Clin Microbiol. 2011;49:1464–74.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Munson MA, Banerjee A, Watson TF, Wade WG. Molecular analysis of the microflora associated with dental caries. J Clin Microbiol. 2004;42:3023–9.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Tanner ACR. Anaerobic culture to detect periodontal and caries pathogens. J Oral Biosci. 2015;57:18–26.

    Article  PubMed  Google Scholar 

  43. Corby PM, Lyons-Weiler J, Bretz WA, Hart TC, Aas JA, Boumenna T, et al. Microbial risk indicators of early childhood caries. J Clin Microbiol. 2005;43:5753–9.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Aas JA, Griffen AL, Dardis SR, Lee AM, Olsen I, Dewhirst FE, et al. Bacteria of dental caries in primary and permanent teeth in children and young adults. J Clin Microbiol. 2008;46:1407–17.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Schulze-Schweifing K, Banerjee A, Wade WG. Comparison of bacterial culture and 16S rRNA community profiling by clonal analysis and pyrosequencing for the characterization of the dentine caries-associated microbiome. Front Cell Infect Microbiol. 2014;4:164.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Claesson MJ, Wang Q, O’Sullivan O, Greene-Diniz R, Cole JR, Ross RP, et al. Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions. Nucleic Acids Res. 2010;38.

  47. Gross EL, Leys EJ, Gasparovich SR, Firestone ND, Schwartzbaum JA, Janies DA, et al. Bacterial 16S sequence analysis of severe caries in young permanent teeth. J Clin Microbiol. 2010;48:4121–8.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Becker MR, Paster BJ, Leys EJ, Melvin L, Kenyon SG, Galvin JL, et al. Molecular analysis of bacterial species associated with childhood caries. J Clin Microbiol. 2002;40:1001–9.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Mantzourani M, Gilbert SC, Sulong HNH, Sheehy EC, Tank S, Fenlon M, et al. The isolation of bifidobacteria from occlusal carious lesions in children and adults. Caries Res. 2009;43:308–13.

    Article  PubMed  Google Scholar 

  50. Tavares WLF, Neves de Brito LC, Teles RP, MLA M, Ribeiro Sobrinho AP, Haffajee AD, et al. Microbiota of deciduous endodontic infections analysed by MDA and Checkerboard DNA-DNA hybridization. Int Endod J. 2011;44:225–35.

    Article  PubMed  Google Scholar 

  51. Ruviére DB, Leonardo MR, da Silva LAB, Ito IY, Nelson-Filho P. Assessment of the microbiota in root canals of human primary teeth by checkerboard DNA-DNA hybridization. J Dent Child (Chic). 2007;74:118–23.

    Google Scholar 

  52. Triches TC, de Figueiredo LC, Feres M, de Freitas SFT, Zimmermann GS, Cordeiro MMR. Microbial profile of root canals of primary teeth with pulp necrosis and periradicular lesion. J Dent Child (Chic). 2014;81:14–9.

    Google Scholar 

  53. Gomes GB, Sarkis-Onofre R, Bonow MLM, Etges A, Jacinto RC. An investigation of the presence of specific anaerobic species in necrotic primary teeth. Braz Oral Res. 2013;27:149–55.

    Article  PubMed  Google Scholar 

  54. Topcuoglu N, Bozdoğan E, Kulekci G. Presence of oral bacterial species in primary endodontic infections of primary teeth. J Clin Pediatr Dent. 2013;38:155–60.

    Article  PubMed  Google Scholar 

  55. Hoshino E. Predominant obligate anaerobes in human carious dentin. J Dent Res. 1985;64:1195–8.

    Article  PubMed  Google Scholar 

  56. Chalmers NI, Oh K, Hughes CV, Pradhan N, Kanasi E, Ehrlich Y, et al. Pulp and plaque microbiotas of children with severe early childhood caries. J Oral Microbiol. 2015;7:1–8.

    Article  Google Scholar 

  57. Rôças IN, Alves FRF, Rachid CTCC, Lima KC, Assunção IV, Gomes PN, et al. Microbiome of deep dentinal caries lesions in teeth with symptomatic irreversible pulpitis. PLoS One. 2016;11:e0154653.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Jiang W, Zhang J, Chen H. Pyrosequencing analysis of oral microbiota in children with severe early childhood dental caries. Curr Microbiol. 2013;67:537–42.

    Article  PubMed  Google Scholar 

  59. Ma C, Chen F, Zhang Y, Sun X, Tong P, Si Y, et al. Comparison of oral microbial profiles between children with severe early childhood caries and caries-free children using the human oral microbe identification microarray. PLoS One. 2015;10:e0122075. https://doi.org/10.1371/journal.pone.0122075.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Whelton H, O’Mullane D, Harding M, Guiney H, Cronin M, Flannery E, et al. North South Survey of Children’s Oral Health in Ireland 2002: DublinBrunswick Press Ltd; 2006.

  61. World Health Organization. Oral Health Surveys – Basic Methods. 4th edn. Geneva: World Health Organization; 1997. https://www.who.int/oral_health/publications/9789241548649/en/.

  62. Whelton H, Crowley E, O’Mullane D, Donaldson M, Cronin M, Kelleher V. Dental caries and enamel fluorosis among the fluoridated population in the Republic of Ireland and non fluoridated population in Northern Ireland in 2002. Community Dent Health. 2006;23:37–43.

    PubMed  Google Scholar 

  63. Gugnani N, et al. International Caries Detection and Assess- ment System ( ICDAS ): A New Concept. Int J Clin Pediatr Dent. 2011;4:93–100.

    Article  PubMed  Google Scholar 

  64. Dye BA, Tan S, Smith V, Lewis BG, Barker LK, Thornton-Evans G, et al. Trends in oral health status: United States, 1988–1994 and 1999–2004. Vital Health Stat 11. 2007;248:1–92.

    Google Scholar 

  65. Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI. The human microbiome project. Nature. 2007;449:804–10.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Aas JA, Paster BJ, Stokes LN, Olsen I, Dewhirst FE. Defining the normal bacterial flora of the oral cavity. J Clin Microbiol. 2005;43:5721–32.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Crielaard W, Zaura E, Schuller AA, Huse SM, Montijn RC, BJF K. Exploring the oral microbiota of children at various developmental stages of their dentition in the relation to their oral health. BMC Med Genomics. 2011;4:22.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Dominguez-Bello MG, Costello EK, Contreras M, Magris M, Hidalgo G, Fierer N, et al. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A. 2010;107:11971–5.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Schloss PD. A high-throughput DNA sequence aligner for microbial ecology studies. PLoS One. 2009;4:e8230.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Allard G, Ryan FJ, Jeffery IB, Claesson MJ. SPINGO: a rapid species-classifier for microbial amplicon sequences. BMC Bioinformatics. 2015;16:324. https://doi.org/10.1186/s12859-015-0747-1.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high- throughput community sequencing data intensity normalization improves color calling in SOLiD sequencing. Nat Publ Gr. 2010;7:335–6.

    Google Scholar 

  72. Statistical RCT (2016). RA language and environment for statistical computing. RF for, computing, Vienna A. No Title. 2016. https://www.r-project.org/. Accessed 27 Feb 2017.

  73. Benjamini Y, Hochberg J. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc B. 1995;57(1):289–300.

    Google Scholar 

  74. Lazarevic V, Whiteson K, Hernandez D, François P, Schrenzel J. Study of inter- and intra-individual variations in the salivary microbiota. BMC Genomics. 2010;11:523.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Marsh PD, Do T, Beighton D, Devine DA. Influence of saliva on the oral microbiota. Periodontol. 2016;70:80–92.

    Article  Google Scholar 

  76. Huang S, Yang F, Zeng X, Chen J, Li R, Wen T, et al. Preliminary characterization of the oral microbiota of Chinese adults with and without gingivitis. BMC Oral Health. 2011;11:33.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Lif Holgerson P, Öhman C, Rönnlund A, Johansson I. Maturation of Oral microbiota in children with or without dental caries. PLoS One. 2015;10:1–20.

    Article  Google Scholar 

  78. Gomar-Vercher S, Cabrera-Rubio R, Mira A, Montiel-Company JM, Almerich-Silla JM. Relationship of children’s salivary microbiota with their caries status: a pyrosequencing study. Clin Oral Investig. 2014;18:2087–94.

    Article  PubMed  Google Scholar 

  79. Fabris AS, Nakano V, Avila-Campos MJ. Bacteriological analysis of necrotic pulp and fistulae in primary teeth. J Appl Oral Sci. 2014;22:118–24.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Santos AL, Siqueira JF, Rôças IN, Jesus EC, Rosado AS, Tiedje JM, et al. Comparing the bacterial diversity of acute and chronic dental root canal infections. PLoS One. 2011;6:e28088.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Hong B-Y, Lee T-K, Lim S-M, Chang SW, Park J, Han SH, et al. Microbial analysis in primary and persistent endodontic infections by using pyrosequencing. J Endod. 2013;39:1136–40.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thanks all the children, parents and staff at the crèches and dental clinic at Cork University Hospital for their help and participation in this study.

Funding

This study was supported by a grant from the Health Research Board (HRB) of Ireland (HRA_POR/2012/123) to PWOT, to support the ORALMET study. The HRB had no role in the collection of data, analysis of data and writing the manuscript. The HRB supported the funding of the study only.

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

All authors have read and approved the manuscript. EH designed the study, performed the experiments, carried out the analyses, and drafted the initial manuscript. HMBH designed the bioinformatic analysis pipeline, supervised and performed the analyses, and reviewed and revised the manuscript. MPJB performed some experiments. CAR, MK and CS critically reviewed the manuscript, and approved the final manuscript as submitted. HW and PWOT conceptualized and designed the study, reviewed and revised the manuscript, and approved the final manuscript as submitted.

Corresponding author

Correspondence to Paul W. O’Toole.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained from the Cork Teaching Hospitals Clinical Research Ethics Committee for the recruitment and sampling of these Cohorts of children. Ethical approval reference: ECM 3 (cc) 01/07/14. Informed and written consent was obtained from the parent or guardian of each child whom participated in this study.

Consent for publication

Not applicable

Competing interests

The authors have no financial relationships relevant to this article to disclose. The authors have no conflicts of interest to disclose.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Food Frequency Questionnaire (FFQ) template used. (XLS 48 kb)

Additional file 2:

Image of Catch-all™ swab used for the collection of oral saliva samples. (DOCX 18 kb)

Additional file 3:

Data collected from FFQ. (XLSX 775 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hurley, E., Barrett, M.P.J., Kinirons, M. et al. Comparison of the salivary and dentinal microbiome of children with severe-early childhood caries to the salivary microbiome of caries-free children. BMC Oral Health 19, 13 (2019). https://doi.org/10.1186/s12903-018-0693-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12903-018-0693-1

Keywords