- Research article
- Open Access
- Open Peer Review
Comparing dental plaque microbiome diversity of extrinsic black stain in the primary dentition using Illumina MiSeq sequencing technique
BMC Oral Health volume 19, Article number: 269 (2019)
Extrinsic black stain (EBS) is characterized by discrete dark dots or lines on the tooth surface. The relationship between EBS and oral microbiota in children remains elusive. The aim of this study was to compare dental plaque microbiome in EBS children with that in EBS-free children in the primary dentition.
The Illumina MiSeq sequencing technique was utilized in the cross-sectional pilot study to investigate the diversity and composition of the supragingival plaque microbiota from 10 EBS-positive and 10 EBS-free children. The results were analysed with nonparametric Mann-Whitney U test, Pearson Chi-Square test, Fisher’s Exact test and one-way ANOVA tests.
We identified 13 different phyla, 22 classes, 33 orders, 54 families, 105 genera, and 227 species from a total of 52,646 high-quality sequences. Between two groups, no statistical differences were observed in the estimators of community richness and diversity at 97% similarity, as well as in the Unweighted Unifrac principal co-ordinates analysis (PCoA). At the species level, higher level of relative abundance of Actinomyces naeslundii and lower level of relative abundance of a species belonging to Candidate_division_TM7 was observed in dental plaque of EBS-positive subjects, compared to dental plaque of EBS-free subjects (P < 0.05). This indicated that some species might be involved in the EBS process.
Changes in dental plaque microbiota is possibly relevant to the process of EBS in the primary dentition.
Extrinsic black tooth stain (EBS, Fig. 1a), which is defined as dark pigmented extrinsic substance in lines or dots parallel to the third cervical line of the tooth crown in the primary and permanent teeth , is often associated with clinical and aesthetic problems .
The prevalence of EBS varies between 2.4 and 18% among different populations . Based on the extension of the tooth surface area affected, Gasparetto et al. introduced the following diagnostic scoring system: score 1 corresponded to the presence of pigmented dots or thin lines with incomplete coalescence parallel to gingival margin; score 2 corresponded to continuous pigmented lines that could easily be observed and were limited to half of the cervical third of the tooth surface; score 3 corresponded to the presence of pigmented stains that extended beyond half of the cervical third of the tooth surface .
The aetiology of EBS is multi-factorial with chromogens derived from dietary and pigmented elements. Possible causes of EBS include the intake of particular foods and beverages, including tea, coffee, and red wine . In addition, the consumption of fruits, dairy products, and dark colour seasonings, such as soy sauce were reported to promote the development of EBS . Other causes of EBS may include chlorhexidine treatment and the presence of metals such as iron .
Previous studies of the microflora of EBS have indicated that black-pigmented bacteria are linked to EBS [8, 9]. It has been reported that Gram-positive rods, especially Actinomyces israelii and Actinomyces naeslundii, are the predominant cultivable microorganisms of EBS . Saba et al. found that Actinomyces spp, had four-fold higher probability of EBS compared to the subjects without these bacteria . Interestingly, Porphyromonas gingivalis and Prevotella melaninogenica were formerly regarded as the main etiologic cause of these pigmentations, however, they were absent in both EBS and control groups. While studying the dental biofilms of patients with or without EBS by multiplex PCR, similar bacterial compositions of EBS-positive and EBS-free groups were identified. However, another research showed similar prevalence of Actinomyces spp in the EBS and non-EBS group (19.2% vs. 11.5%, respectively) , which did not match the findings reported by Saba et al. The results differ between studies, suggesting that the composition of the EBS-associated microbiota has not been definitively identified.
The oral cavity harbours one of the most complex microbiomes in the body . Traditional culture-dependent studies that focused on mono-species or certain types of microbiome underestimated the mixed species oral microbial communities , therefore, efforts to investigate microbial diversity increasingly relied on cultivation-independent, molecular approaches . DAN-based methods, such as next-generation FLX+ and, Illumina pyrosequencing methods, combined with curated gene databases, such as the Human Oral Microbiome Database , are comprehensive methods to analyse the microbiota. MiSeq, Illumina’s next generation sequencing techonology, has been wildly used nowadays as it is intuitively multi-function, including amplification, sequencing and data analysis. By using MiSeq, it is capable of identifying the bacteria composition from complex environments in a culture-free way .
The aim of the present study was to evaluate the plaque microbiota in children with high prevalence of EBS and who had no EBS by using the Illumina MiSeq, thereby identifying a potential role for the microbiome in EBS among children. While we expand our efforts to identify the dental plaque microbiome diversity of EBS in children with primary dentition will greatly aid in our understanding of the potential role of bacteria in the process of EBS as limited research has been conducted regarding microbiology by next-generation sequencing methods.
Materials and methods
A total of 100 children in primary dentition were recruited for this cross-sectional pilot study by the Department of Paediatric Dentistry, West China Hospital of Stomatology, Sichuan University (Chengdu, China). The oral health status of all participants was determined by a dentist who performed a full-mouth clinical examination that included inspection of the teeth, oral mucosa, and periodontal tissues.
Based on the clinical examination according to criteria defined by Gasparetto et al., 10 children who presented EBS (sore 2–3) pigmented stains that were limited to or extended beyond half of the cervical third of the tooth surface and were negative to other extrinsic discoloration (e.g. brown, green or orange stain), were included in this study (Group BP).
For each child with EBS, a gender-, age-, and dentition status-matched EBS-free counterpart (10 children in total, Group BFP) was recruited as a control. All children didn’t have antibiotics or any type of polishing within 3 months prior to the start of the study. Children with other oral disease, such caries, periodontitis or intrinsic BS, or systematic diseases were excluded from the study. Inclusion and exclusion criteria are listed in Table 1.
Plaque samples were collected from supragingival tooth surfaces using sterilized curettes. For Group BP, samples were collected from the labial or buccal surfaces of EBS teeth and pooled. For Group BFP, specimens were obtained from the corresponding positions and surfaces of the teeth. Pooled plaque samples were removed from the curette by agitation in 700 ml of Tries-EDTA (TE) buffer (10 mM Tris-Cl [pH 7.5] and 1 mM EDTA) and stored at − 80 °C for further studies.
DNA extraction and PCR amplification
The total bacterial genomic DNA was extracted from pooled plaque samples using the QIAamp DNA Micro Kit (Qiagen, Hilden, Germany) according to manufacturer’s protocols. The V3-V4 hypervariable region of the bacteria 16S ribosomal RNA gene were amplified by PCR using primers 338F 5′-barcode-ACTCCTACGGGAGGCAGCAG-3′ and 806R 5′-GGACTACHVGGGTWTCTAAT − 3′, in which the barcode represented an eight-base sequence that was unique to each sample . PCR reactions were performed in triplicate in a total of 20 μl mixture, containing 4 μl of 5 × FastPfu Buffer, 2 μl of 2.5 mM dNTPs, 0.8 μl of each primer (5 μM), 0.4 μl of FastPfu Polymerase, and 10 ng of template DNA. The PCR conditions were as follows: 94 °C for 3 min for DNA denaturation, followed by 27 cycles (94 °C for 30 s, 50 °C for 30 s, and 72 °C for 45 s), and a final extension for 10 min at 72 °C. The amplicons of three replicates were pooled to obtain sufficient PCR products.
Illumina MiSeq sequencing
After PCR, amplicons were separated using 2% agarose gels and extracted and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA). Then, amplicons were quantified using QuantiFluor™ -ST (Promega, Madison, WI, USA) according to the manufacturer’s instructions. Purified amplicons were pooled in equimolar and paired-end sequenced (2 × 250) on an Illumina MiSeq platform (Majorbio, Shanghai, China) according to the standard protocols. The raw data were deposited into the NCBI Sequence Read Archive (SRA) database for evaluation.
Processing of sequencing data
Raw FASTQ files were de-multiplexed, quality-filtered using QIIME (version 1.9.1) using the following criteria: (i) 300 bp reads were truncated at any site receiving an average quality score < 20 over a 50 bp sliding window. Truncated reads shorter than 50 bp were eliminated; (ii) exact barcode matching reads (2 nucleotide mismatch in primer matching) containing ambiguous characters were removed; (iii) based to their overlap sequence, only sequences with an overlap of > 10 bp were assembled. Reads that could not be assembled were eliminated.
UPARSE software (version 7.1 http://drive5.com/uparse/) was used to cluster quality-checked sequences into Operational Taxonomic Units (OTUs) at a similarity threshold of 97% for production of OTUs and removed using UCHIME .
The taxonomy of each 16S rRNA gene sequence was analysed by Ribosomal Database Project (RDP) Classifier (http://rdp.cme.msu.edu/) against the SILVA (SSU123) 16S rRNA database using a confidence threshold of 70% . The relative abundance of bacterial genera was calculated using QIIME software , and rarefaction curve analysis, community richness (ACE, Chao1), and community diversity indices (Simpson, Shannon diversity index) were determined by the MOTHUR program at a 97% similarity level. To examine dissimilarities in community composition, principal co-ordinates analysis (PCoA) in QIIME was performed, which was used to compare groups of samples based on weighted and weighted UniFrac distance metrics . Using Vegan packages in R, heat map profiles were generated.
The independent sample nonparametric Mann-Whitney U test, Pearson Chi-Square test, Fisher’s Exact test and one-way ANOVA test were performed using SPSS version 23.0 (SPSS Inc., Chicago, IL, USA).
Sequencing results and diversity indices
In this study, a total of 552,646 high-quality sequences were produced, with an average of 27,632 sequences per sample. Summary information is shown in Table 2, and detailed characteristics of each sample are presented in Additional file 1: Table S1. The numbers of OTUs (3% dissimilarity) and detailed characteristics for each sample are shown in Additional file 1: Table S1. The average OTUs of Group BP and BFP were 153.7 ± 32.2 and 156.4 ± 24.6, respectively (P > 0.05). The rarefaction curves of Group BP and BFP reached a saturation plateau, indicating that the data volume of sequenced reads was reasonable (Additional file 2: Figure S1). Comparisons of the estimators of community richness (Chao1 and Ace) and diversity (Shannon index and Simpson indices) at 97% similarity between two groups are shown in Table 3. One-way ANOVA analysis did not show significant differences in α-diversity between Group BP and BFP (P > 0.05).
Taxonomy-based comparisons of oral microbiota between group BP and BFP
To investigate whether EBS was correlated with potential changes in the abundance of specific bacterial taxa, the relative abundance of taxa was compared.
Sequences from these samples were classified into 13 different phyla, 22 classes, 33 orders, 54 families, 105 genera, and 227 species. The overall microbiota structure for each group at the phylum level is presented in Fig. 2a. Out of a total of 13 phonotypes, six different phyla (Formicates, Fusobacteria, Proteobacteria, Bacteroidetes, Actinobacteria, and Candidate_division_TM7) were common, and comprised a considerable proportion of 99.5%. Firmicutes were strongly enriched, accounting for 28.9% of the total sequences in Group BP, followed by Fusobacteria (24.4%), Proteobacteria (18.3%), Bacteroidetes (16.2%), Actinobacteria (8.3%), and Candidate_division_TM7 (3.6%). Moreover, microbiota in Group BFP primarily included Firmicutes (29.7%), Fusobacteria (25.3%), Bacteroidetes (16.1%), Candidate_division_TM7 (11.6%), Proteobacteria (11.4%), and Actinobacteria (5.3%).
Statistical differences between the two groups were found only for Candidate_division_TM7 (P < 0.05). The relative abundance of the dominant genera identified in the two groups is presented in Fig. 2b (only those with a relative abundance of > 1% are listed). In Group BP, Actinomyces, Candidate_division_TM7_norank, Capnocytophaga, Centipeda, Corynebacterium, Porphyromonas, and Fusobacterium exhibited a relatively higher abundance compared to Group BFP. Gemella, Haemophilus, Lautropia, Leptotrichia, Neisseria, and Neisseriaceae_unclassified were relatively more abundant in Group BFP. Statistically significant differences between Group BP and BFP were only found in genus Candidate_division_TM7_norank (P < 0.05).
Figure 2c showed the relative abundance of the major bacterial community at the species levels (only those above 1% in relative abundance are listed). Our data indicated that the bacterial community composition showed a few differences between Group BP and BFP. For example, species Actinomyces naeslundii were significantly more abundant in Group BP compared to Group BFP, whereas species TM7_phylum_sp._canine_oral_taxon_308 was significantly less abundant in Group BP compared to Group BFP (P < 0.05).
The top 100 abundant genera among the 20 samples were displayed based on the heat map of the genus-level plaque composition (Fig. 4), and indicated that between Group BP and BFP most genera were shared. In both groups, the dominant genera were Neisseria, Streptococcus, Fusobacterium, Capnocytophaga, and Veilonella. Moreover, minor genera, such as Gemella, Lactobacillales, and Stomatobaculum were found in both groups. However, some genera, including genoceanbacillus, myroidesu, brochotrix, were predominantly found in Group BFP.
Comparison of beta diversity between groups
In this study, PCoA based on weighted UniFrac matrix was used to identify community structure differences (Fig. 3), and showed that the oral microbial structure of Group BP differed from that of Group BFP. On the primary axes, the PCoA ordination did not reveal strong grouping of disease and health, however, a segregation trend for EBS-positive patients and EBS-free subjects was observed.
The human mouth harbours a highly diverse bacterial community, however, only about 50% of oral microorganisms can be cultivated . In previous studies, the microbiology of EBS has been investigated. However, these studies were impeded by the limitation of conventional culture-dependent approaches , or only focused on a few selected bacteria [10, 22], without characterizing the overall community of oral microbiota in EBS patients. Accordingly, we used high-throughput DNA sequencing technique, Illumina Miseq to investigate the holistic pattern of dental biofilm microbiome. Moreover, bacteria with low abundance or not-yet-cultivated members could be investigated in children with or without EBS .
No statistical differences were found from the Chao1 and Ace estimator between Group BP and BFP. In addition, a high similarity of Shannon and Simpson index was observed between the two groups, demonstrating that similar dental biofilm microbial community richness and diversity in these two groups. These findings were in line with the recent findings reported by Yue Li, et al. . The difference was not statistically significant, which was likely due to the small sample size of our pilot study, however a segregation trend was observed in the Unweighted Unifrac PCoA analysis. Additional research, with larger study groups, is required to verify the role of EBS in shaping the overall pattern of the dental biofilm microbiome community.
In our study, we haven’t observed significant differences in alpha and beta diversity, however, we did identify potential microbiota associated with EBS in children. Our findings indicated that Candidate_division_TM7 was less abundant in the supragingival plaque of Group BP compared to Group BFP from the phylum down to the species level. To our knowledge, this is the first study that reported significant differences between the prevalence of Candidate_division_TM7 in children with or without EBS.
Candidate_division_TM7, initially discovered from the German peat bog by environmental sequencing data, represents one of uncultivable bacterial divisions . Using next-generation sequencing strategies, highly divergent members of this division were identified in a range of cosmopolitan habitats, including terrestrial, aquatic, and human body sites . Recently, the ability to modulate the oral biofilm formation of TM7 was reported. The morphology of TM7 in dual-species biofilms with different types of bacteria varied, as did the morphological characteristics of biofilm. It was speculated that the interactions between TM7 and other microorganisms affected the oral microbial community and the process of infectious oral diseases . In present study, Candidate_division_TM7 was reduced in Group BP compared to Group BFP, however the mechanism remains elusive. We hypothesized a possible correlation between the decrease in Candidate_division_TM7 and the development of EBS through complicated bacterial interactions. Further studies are required to clarify the underlying mechanism.
Furthermore, our study is consistent with previous reports [22, 24] that Actinomyces naeslundii was more prevalent in supragingival plaque of EBS patients. A. naeslundii is one of the resident microbiota in the human mouth and plays a considerable role in dental biofilm formation . Sarkonen et al.  demonstrated that other Actinomyces species, including A. odontolyticus, could form brown to black pigmentation. In addition, Reid et al. considered that the EBS was formed as a result of the reaction between ferric ions in saliva and the hydrogen sulfide produced by microorganisms (for example Actinomycetes) .
Prevotella spp. was identified as the predominant bacteria of black pigmentations in earlier studies , while, our study showed no significant differences regarding of the prevalence of Prevotella spp. between Group BP and BFP, which is consistent with most recent studies [9, 10, 22].
Based on the results of the current study, our data demonstrated that children in primary dentition with or without EBS share a similar dental biofilm microbial community. Candidate_division_TM7 and Actinomyces naeslundii may possibly be involved in the presence of EBS. Although the sample size in our study was relatively small and further experiments need to be done, our findings represented a comprehensive picture of the microbial structure of children with EBS.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Extrinsic black stain
Plaque samples from EBS-positive subjects
Plaque samples from EBS-free subjects
Operational Taxonomic Units
Principal co-ordinates analysis
Watts A, Addy M. Tooth discolouration and staining: a review of the literature. Br Dent J. 2001;190(6):309–16.
Ronay V, Attin T. Black stain - a review. Oral Health Prev Dent. 2011;9(1):37–45.
Żyła T, Kawala B, Antoszewska-Smith J, Kawala M. Black stain and dental caries: a review of the literature. Biomed Res Int. 2015;2015:469392. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4354720/.
Gasparetto A, Conrado CA, Maciel SM, Miyamoto EY, Chicarelli M, Zanata RL. Prevalence of black tooth stains and dental caries in Brazilian schoolchildren. Braz Dent J. 2003;14(3):157–61.
Pearson D. The chemical analysis of foods: Longman Group Ltd; 1976. http://trove.nla.gov.au/work/3564386?q&versionId=45250339.
Chen X, Zhan JY, Lu HX, Ye W, Zhang W, Yang WJ, et al. Factors associated with black tooth stain in Chinese preschool children. Clin Oral Investig. 2014;18(9):2059–66.
Addy M, Moran J, Griffiths A, Wills-Wood N. Extrinsic tooth discoloration by metals and chlorhexidine. I. Surface protein denaturation or dietary precipitation? Br Dent J. 1985;159(9):281–5.
Slots J. The microflora of black stain on human primary teeth. Scand J Dent Res. 1974;82(7):484–90.
Saba C, Solidani M, Berlutti F, Vestri A, Ottolenghi L, Polimeni A. Black stains in the mixed dentition: a PCR microbiological study of the etiopathogenic bacteria. J Clin Pediatr Dent. 2006;30(3):219–24.
Costa MT, Dorta ML, Ribeiro-Dias F, Pimenta FC. Biofilms of black tooth stains: PCR analysis reveals presence of Streptococcus mutans. Braz Dent J. 2012;23(5):555–8.
Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, et al. The human oral microbiome. J Bacteriol. 2010;192(19):5002–17.
Amann RI, Ludwig W, Schleifer KH. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol Rev. 1995;59(1):143–69.
Schloss PD, Handelsman J. Status of the microbial census. Microbiol Mol Biol Rev. 2004;68(4):686–91.
Chen T, Yu WH, Izard J, Baranova OV, Lakshmanan A, Dewhirst FE. The human oral microbiome database: a web accessible resource for investigating oral microbe taxonomic and genomic information. Database (Oxford). 2010;2010. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2911848/.
Ravi RK, Walton K, Khosroheidari M. MiSeq: a next generation sequencing platform for genomic analysis. Methods Mol Biol. 1706;2018:223–32.
Dennis KL, Wang Y, Blatner NR, Wang S, Saadalla A, Trudeau E, et al. Adenomatous polyps are driven by microbe-instigated focal inflammation and are controlled by IL-10 producing T-cells. Cancer Res. 2013;73(19):5905–13.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Amato KR, Yeoman CJ, Kent A, Righini N, Carbonero F, Estrada A, et al. Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes. ISME J. 2013;7(7):1344–53.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.
Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011;5(2):169–72.
Paster BJ, Olsen I, Aas JA, Dewhirst FE. The breadth of bacterial diversity in the human periodontal pocket and other oral sites. Periodontol 2000. 2006;42:80–7.
Heinrich-Weltzien R, Bartsch B, Eick S. Dental caries and microbiota in children with black stain and non-discoloured dental plaque. Caries Res. 2014;48(2):118–25.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq illumina sequencing platform. Appl Environ Microbiol. 2013;79(17):5112–20.
Li Y, Zhang Q, Zhang F, Liu R, Liu H, Chen F. Analysis of the microbiota of black stain in the primary dentition. PLoS One. 2015;10(9):e0137030.
Rheims H, Spröer C, Rainey FA, Stackebrandt E. Molecular biological evidence for the occurrence of uncultured members of the actinomycete line of descent in different environments and geographical locations. Microbiology. 1996;142(10):2863–70.
Hugenholtz P, Tyson GW, Webb RI, Wagner AM, Blackall LL. Investigation of candidate division TM7, a recently recognized major lineage of the domain bacteria with no known pure-culture representatives. Appl Environ Microbiol. 2001;67(1):411–9.
Soro V, Dutton LC, Sprague SV, Nobbs AH, Ireland AJ, Sandy JR, et al. Axenic culture of a candidate division TM7 bacterium from the human oral cavity and biofilm interactions with other oral bacteria. Appl Environ Microbiol. 2014;80(20):6480–9.
Henssge U, Do T, Radford DR, Gilbert SC, Clark D, Beighton D. Emended description of Actinomyces naeslundii and descriptions of Actinomyces oris sp. nov. and Actinomyces johnsonii sp. nov., previously identified as Actinomyces naeslundii genospecies 1, 2 and WVA 963. Int J Syst Evol Microbiol. 2009;59(3):509–16.
Sarkonen N, Könönen E, Summanen P, Könönen M, Jousimies-Somer H. Phenotypic identification of Actinomyces and related species isolated from human sources. J Clin Microbiol. 2001;39(11):3955–61.
Reid JS, Beeley JA, MacDonald DG. Investigations into black extrinsic tooth stain. J Dent Res. 1977;56(8):895–9.
Theilade J. Development of bacterial plaque in the oral cavity. J Clin Periodontol. 1977;4(5):1–12.
The study design, the collection, analysis, and interpretation of data and writing the manuscript were funded by the National Natural Science Foundation of China (81400502 to QZ, 81470035 to JZ and 81600864 to YW).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of State Key Laboratory of Oral Diseases, Sichuan University (Chengdu, China) (WCHSIRB-D-2013-128-R1). Written informed consent was obtained from the parents prior to the start of the study. This study was in accordance with the Ethical Guidelines of the Declaration of Helsinki (First revision, 1975). The authors confirm that all ongoing and related trials for these experiments are registered. All experiments were performed in accordance with relevant guidelines, and informed written consent was obtained from all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, L., Zhang, Q., Wang, Y. et al. Comparing dental plaque microbiome diversity of extrinsic black stain in the primary dentition using Illumina MiSeq sequencing technique. BMC Oral Health 19, 269 (2019) doi:10.1186/s12903-019-0960-9
- Extrinsic black stain
- Dental biofilm
- Primary dentition
- Illumina MiSeq