Apical periodontitis: preliminary assessment of microbiota by 16S rRNA high throughput amplicon target sequencing
BMC Oral Health volume 18, Article number: 55 (2018)
Apical periodontitis includes periapical granulomas and radicular cysts, which are histologically distinguished by the absence and the presence of an epithelial lining, respectively. The main cause of apical periodontitis is the bacterial colonization of the root canal space. This research aimed at assessing whether and how periapical granulomas and radicular cysts differ in terms of microbiota using high throughput amplicon target sequencing (HTS) techniques.
This study included 5 cases of Periapical Granulomas (PGs) and 5 cases of Radicular Cysts (RCs) selected on the base of histology out of 37 patients from January 2015 to February 2016. Complete medical history, panoramic radiograms (OPTs) and histologic records of each patient were assessed. Only lesions greater than 1 cm in diameter and developed in proximity to teeth with bad prognosis were included. The microbiota present in periapical granulomas and radicular cysts thus retrieved was finely characterized by pyrosequencing of the 16S rRNA genes.
The core of OTUs shared between periapical granulomas and radicular cysts was dominated by the presence of facultative anaerobes taxa such as: Lactococcus lactis, Propionibacterium acnes, Staphylococcus warneri, Acinetobacter johnsonii and Gemellales. L. lactis, the main OTUs of the entire datasets, was associated with periapical granuloma samples. Consistently with literature, the anaerobic taxa detected were most abundant in radicular cyst samples. Indeed, a higher abundance of presumptive predicted metabolic pathways related to Lipopolysaccharide biosynthesis was found in radicular cyst samples.
The present pilot study confirmed the different microbial characterization of the two main apical periodontitis types and shade light on the possible role of L. lactis in periapical granulomas.
Apical periodontitis (AP) is associated with endodontically involved teeth [1, 2]. In most cases, it is impossible to distinguish between periapical granulomas (PGs) and radicular cysts (RCs), without recurring to biopsy . The occurrence of periapical granulomas ranges between 9.3 and 87.1% . Radicular cysts are believed to form by proliferation of the epithelial cell rests of Malassez in inflamed periradicular tissues . Whether they be pocket cysts (cavity open to the root canal) or true cysts (completely enclosed by lining epithelium) , their reported incidence among periapical lesions varies from 6% to 55% . As radiography for determining APs has been questioned for scientific investigations , differential diagnosis is possible only with histopathological examination . Whatever the diagnosis, root canal debridement is the first choice treatment of APs [10, 11].
Bacterial colonization of the root canal space has been demonstrated as the main etiologic factor of APs [12, 13]. In two paradigmatic studies by Ricucci & Siqueira [13, 14], bacterial biofilm varied and no unique pattern for endodontic infections was identified. Bacteria could modify the severity and prognoses of APs and yet, surprisingly, little information is available in the scientific literature comparing the microbiota within PGs and RCs. The application of high throughput amplicon target sequencing (HTS) to study the microbial ecology has been witnessed over the past couple of years aimed at estimating the microbial diversity in different ecosystems using 16S rRNA gene as the target. The HTS provides an unprecedented greater sampling depth and allows the detection not only of the dominant community members but also of low-abundance taxa .
Flurry of research has been carried out in past decades to assess the microbiota of the human oral cavity, as well as the endodontic microbiome [14, 16, 17]. Hence, the purpose of this study was to finely characterize the microbiota present in histologically determined PGs and RCs by pyrosequencing their 16S rRNA genes. An in depth analysis of the PGs and RCs microbiota is required for a better understanding of the bacterial taxa involved with the inflammation process.
Study design and patients
The study was planned and performed in accordance with the Declaration of Helsinki and was approved by the Ethics committee of the Dental School, University of Turin. From January 2015 to February 2016, 121 patients with apical periodontitis were referred to the Triage of the Dental School of the University of Turin. Complete medical history and panoramic radiograms (OPTs) of each patient were assessed seeking large Apical periodontits (APs) that clearly had developed in proximity to teeth with bad prognosis. In the Oral Surgery Department, 37 patients were selected by applying the following exclusion criteria: systemic or local disease or condition (hematologic diseases, uncontrolled diabetes, serious coagulopathies, history of intravenous therapy with bisphosphonates, and/or diseases of the immune system) possibly precluding oral surgical intervention; immunosuppression; HIV+, HCV+, HBV+, TBC+, corticosteroid treatment, pregnancy, radiotherapy to the head or neck region within 12 months before surgery. A further restriction was achieved excluding teeth a) periodontally compromised, b) with endodontic communication to the oral cavity and c) vertical tooth fractures. Finally, after giving their informed consent, 10 patients underwent surgery. By careful manipulation, the periapical lesions were harvested sterilely and fixed in 4% formalin, while the hopeless teeth were extracted. Based on histology and a dimensional cut-off (lesion greater than 1 cm in diameter), among the apical periodontal lesions analyzed, 5 Radicular Cysts (RCs) and 5 Periapical Granulomas (PGs) were retrieved, dependently on the clear presence or absence of cavity lining epithelium, respectively.
The histological specimens retrieved were fixed in 4% formalin for 24 h and subsequently embedded with paraffin wax and cut into 3 μm thick sections, using a motorized microtome. Polylysine coated slides were used to enhance the adhesion of the tissue section during staining procedures. The histological structure of the lesions was assessed by traditional hematoxylin and eosin staining for optical microscopy.
DNA analysis by pyrosequencing
Two slices of formalin-fixed, paraffin-embedded (FFPE) tissue samples (about 10 mg of tissue) were used for total genomic DNA extraction. Samples were pre-treated at 55 °C for a minimum of 1 h with the dissolving buffer and Proteinase K (20 mg/ml) according to the manufacturer’s instructions (BiOstic® FFPE Tissue DNA Isolation Kit Mobio, Carlsbad, CA, USA). DNA was used to study the microbial diversity by pyrosequencing of the amplified V1–V3 region of the 16S rRNA gene, recurring to the primers Gray28F (5’-TTTGATCNTGGCTCAG) and Gray519r (5’-GTNTTACNGCGGCKGCTG) that amplify a fragment of 520 bp, following PCR conditions previously reported . PCR products were purified twice with Agencourt AMPure purification kit (Beckman Coulter, Milan, Italy), and quantified using the PlateReader AF2200 (Eppendorf, Hamburg, Germany) with PicoGreen assay and an equimolar pool was obtained prior to further processing. Due to poor DNA quality, PG_8 sample was excluded. The amplicon pool was used for pyrosequencing on a GS Junior platform (454 Life Sciences, Roche, Monza, Italy) according to the manufacturer’s instructions by using Titanium chemistry.
Bioinformatics analysis and metagenomic prediction
QIIME 1.9.0 software was employed to analyze 16S rRNA data : OTUs (operational taxonomic units) were picked at 99% of similarity by means of UCLUST clustering methods . Representative sequences from each cluster were used to assign taxonomy through matching against the Greengenes 16S rRNA gene database version 2013 by the RDP classifier. R environment (www.r-project.org) was adopted to elaborate statistics as well as plotting. To calculate the microbiota alpha diversity the authors chose the “diversity” function of the vegan package of R. Weighted UniFrac distance matrices as obtained through QIIME were imported in R to generate PCoA (Principal Coordinates Analysis) plots. Weighted UniFrac distance matrices were also used to perform ADONIS and ANOSIM statistical tests owing to compare_categories.py script of QIIME. OTU tables, which were filtered at 0.2% abundance in at least two samples, were used to compare each OTU based on the passed sample groupings (PGs and RCs) through the group_significance.py script of QIIME. OTUs co-occurrence co-exclusion was carried out by the psych package of R (www.r-project.org) and it was further visualized through the corrplot package of R . In order to predict the inferred metagenome, the authors used PICRUSt  so as to predict abundances of gene families based on 16S rRNA sequences data . Briefly, the pick OTUs step was re-performed at 97% similarity against the Greengenes database and the resulting KEGG orthologs table was then collapsed at level 3 of the KEGG annotations in order to display the inferred metabolic pathways. The resulting table was imported in the GAGE Bioconductor package  to identify biological pathways overrepresented or underrepresented between PGs and RCs samples. To characterize the accuracy of PICRUSt, the Nearest Sequenced Taxon Indexes (NSTI) were calculated . All the sequencing data were deposited at the Sequence Read Archive of the National Center for Biotechnology Information (accession numberSRP096711).
APs were subdivided into PGs and RCs on the base of the histological analysis performed on the biopsy samples. The relevant data concerning patients’ age and gender are reported in Table 1, along with dimensions and region of the lesion. Also any previous canal root therapy was recorded. It is to be noted that all patients were older than 18 years and were in good health conditions as per inclusion criteria (Table 1).
A total of 163.832 raw reads were obtained after the sequencing. 65.070 reads passed the filters applied through QIIME, with an average value of 7.230 reads/sample and a sequence length of 486 bp. The rarefaction analysis and the estimated sample coverage (Table 2) indicated that there was a satisfactory coverage for all the samples (ESC > 96%). The richness of the samples varied from a minimum of 22 to a maximum of 202 OTUs. (Table 2) Alpha-diversity indices (Table 2) showed no difference on the level of complexity (P > 0.05) of RCs samples compared to PGs.
In Fig. 1, the box plot shows the OTUs with a relative abundance of 0.2% in at least two samples (Fig. 1). The data showed a varied microbiota composition characterized by the presence of 40% of facultative microbes , 37% of anaerobes and 22% of aerobe microbes (Table 3). In details, the samples were characterized by the predominance of Lactococcus lactis (55% of the relative abundance), Propionibacterium acnes (18%), Corynebacterium matruchotii (5.5%), Staphylococcus warneri (5%), Gemellales (2%), Actinomyces johnsonii (2%), and Lactobacillus zeae (2.5%). Through principal coordinate analysis (PCoA) with a weighted UniFrac distance matrix, it was possible to show that samples from RCs grouped together and that they were well separated from PGs on the basis of their microbiota (Fig. 2); ADONIS and ANOSIM statistical tests confirmed this difference (P < 0.001). The differential abundance analysis showed a higher abundance (Bonferroni corrected P value of < 0.001) of several minor OTUs in RCs compared to PGs samples. In particular, it was possible to observe the most abundant presence of several facultative anaerobes or anaerobe OTUs such P. acnes, Gemellales, Capnocytophaga ochracea, Paracoccus, Fusobacterium nucleatum, Prevotella intermedia and Rothia dentocariosa. However L. lactis was found significantly more abundant (P < 0.001) in PGs samples than in RCs samples.
The OTU co-occurrence was investigated by considering the species-level taxonomic assignment and significant correlations at false-discovery rate [FDR], < 0.05. (Fig. 3). P. intermedia showed the highest number of positive correlations including those with Streptococcus mitis and F. nucleatum and a co-exclusion with Sphyngomonas sp. Moreover the most significant OTUs in cyst samples such as Gemellales, C. ochracea, showed the highest number of positive correlation. The core OTUs L. lactis co-exclude the presence of Acinetobacter lwoffii, while P. acnes co-exclude the presence of S. mitis. (Fig. 3).
Regarding the predicted metagenomes, the weighted nearest-sequenced-taxon index (NSTI) for the samples, expressed as the mean SD, was 0.042 ± 0.004. Thus, a NSTI score of 0.042 indicates a satisfactory accuracy for all of the samples (96%). The pathway enrichment analysis (performed by GAGE) of the predicted metagenomes showed an enrichment of metabolic pathways such as Biosynthesis of amino acids (ko01230), Pyruvate metabolism (ko00620), Propanoate metabolism (ko00640) in PGs samples compared to RCs samples (data not shown). In contrast, from RCs samples only pathways involved in cellular processes, biosynthesis of secondary metabolites, and genes involved in Lipopolysaccharide biosynthesis (ko00540) were found.
The oral cavity is exposed to the external environment and is, therefore, one of the most important ways of microbial entry into the human body . By invading the adjacent tissues, bacteria may induce an immune response resulting in inflammatory manifestations such as apical periodontitis [24, 25]. The presence of bacteria in PGs and RCs was previously confirmed . Recently, by culture dependent methods , RCs were clearly demonstrated to possess a great variety of anaerobic and facultative anaerobic microbial taxa. Our results showed that the core of OTUs shared between PGs and RCs was dominated by the presence of facultative anaerobes taxa such as: L. lactis, P. acnes, S. warneri, A. johnsonii and Gemellales. In particular, P. acnes, reported as the most commonly detected bacterium, has been studied due to its capacity to induce the differentiation of T lymphocytes into CD25 regulatory bright cells with a potentially inhibitory effect on the immune response . Actinomyces species have been implicated frequently as a cause of endodontic failure because of their ability to persist in periapical tissues [27,28,29,30].
These species are all normal commensals of the human oral cavity and were isolated in radicular cyst . Beta diversity calculation as well as ADONIS and ANOSIM statistical tests display a degree of separation of the samples due to the relative abundance of the minor OTUs. Of course, within the oral cavity, many initial interactions between food microbes and human microbiota occur. Our results revealed the presence of taxa clearly derived from food like L. lactis, a non pathogenic taxon, usually not associated with the oral microbiota. Although the genus Lactococcus had already been isolated by culture dependent methods from PGs , here we show, for the first time, that L. lactis was the main OTUs of the entire datasets and it was associated with PGs samples.
It may be noteworthy that the anaerobic taxa detected were most abundant in RCs samples. F. nucleatum, P. intermedia and R. dentocariosa were observed to be statistically more abundant in those samples and, as previously reported, the presence of these taxa suggested the onset of secondary infection . In addition, Iatrou et al.  determined in their study that the isolated bacteria were mostly anaerobic. Aerobe and facultative anaerobic bacteria growth was seen in 10.8% of the cases. The OTU co-occurrence analysis displays the strong correlation between the presence of F. nucleatum and P. intermedia. Fusobacteria are present in apical abscesses because they constitute an important part of the apical biofilm [27, 32]. In particular F. nucleatum was frequently isolated and cultured from teeth with apical periodontitis and its virulence is greatly enhanced in presence of P. intermedia .
Predicted metagenomes confirmed differences between the two types of samples and indicated that the RCs samples displayed a higher abundance of presumptive predicted metabolic pathways related to Lipopolysaccharide biosynthesis. This metabolic pathway is indicative of the presence of Gram negative bacteria and it is widely accepted as a subclinical pro inflammation marker . The putative role of bacterial endotoxins in supporting epithelial proliferation typical of RCs has already been reported . According to Meghji et al. , when epithelial cell proliferation assays were performed, the Lipopolysaccharides derived form three different bacteria displayed mitogenic effects, which were even enhanced if cyst fibroblast culture media were used.
The present research, albeit preliminary, may contribute to the progress of the existing knowledge concerning the microbiota within the two main apical periodontal lesions. The use of sophisticated and sensitive techniques allowed unprecedented results and could help elucidating the possible etiopathologic role of a complex microbial environment in either promoting or refraining the epithelial proliferation.
High throughput amplicon target sequencing
Operational taxonomic unit
Principal coordinate analysis
Lin LM, Huang GT-J, Rosenberg PA. Proliferation of epithelial cell rests, formation of apical cysts, and regression of apical cysts after periapical wound healing. J Endod. 2007;33:908–16.
Poswar FDO, Farias LC, Fraga CADC, Bambirra W, Brito-Júnior M, Sousa-Neto MD, et al. Bioinformatics, interaction network analysis, and neural networks to characterize gene expression of radicular cyst and periapical granuloma. J Endod. 2015;41:877–83.
Xia W, Zhu Y, Wang X. Six cases report of differential diagnosis of periapical diseases. Int J Oral Sci. 2011;3:153–9.
Schulz M, von Arx T, Altermatt HJ, Bosshardt D. Histology of periapical lesions obtained during apical surgery. J Endod. 2009;35:634–42.
Shear M. Cysts of the jaws: Recent advances. J Oral Pathol. 1985;14:43–59.
Simon JHS. Incidence of periapical cysts in relation to the root canal. J Endod. 1980;6:845–8.
Ramachandran Nair PN, Pajarola G, Schroeder HE. Types and incidence of human periapical lesions obtained with extracted teeth. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 1996;81:93–102.
López FU, Kopper PMP, Cucco C, Della Bona A, de Figueiredo JAP, Vier-Pelisser FV. Accuracy of cone-beam computed tomography and periapical radiography in apical periodontitis diagnosis. J Endod. 2014;40:2057–60.
Alcantara BAR, de Carli ML, Beijo LA, Pereira AAC, Hanemann JAC. Correlation between inflammatory infiltrate and epithelial lining in 214 cases of periapical cysts. Braz Oral Res. 2013;27:490–5.
Tavares PBL, Bonte E, Boukpessi T, Siqueira JF, Lasfargues J-J. Prevalence of apical periodontitis in root canal–treated teeth from an urban french population: Influence of the quality of root canal fillings and coronal restorations. J Endod. 2009;35:810–3.
Song M, Kim H-C, Lee W, Kim E. Analysis of the cause of failure in nonsurgical endodontic treatment by microscopic inspection during endodontic microsurgery. J Endod. 2011;37:1516–9.
Kakehashi S, Stanley HR, Fitzgerald RJ. The effects of surgical exposures of dental pulps in germ-free and conventional laboratory rats. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 1965;20:340–9.
Ricucci D, Pascon EA, Pitt Ford TR, Langeland K. Epithelium and bacteria in periapical lesions. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2006;101:239–49.
Siqueira JF, Rôças IN, Alves FRFF, Silva MG. Bacteria in the apical root canal of teeth with primary apical periodontitis. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2009;107:721–6.
De Filippis F, Parente E, Ercolini D. Metagenomics insights into food fermentations. Microb Biotechnol. 2017;10:91–102.
Li L, Hsiao WWL, Nandakumar R, Barbuto SM, Mongodin EF, Paster BJ, et al. Analyzing endodontic infections by deep coverage pyrosequencing. J Dent Res. 2010;89:980–4.
De Filippis F, Vannini L, La Storia A, Laghi L, Piombino P, Stellato G, et al. The same microbiota and a potentially discriminant metabolome in the saliva of omnivore, ovo-lacto-vegetarian and vegan individuals. PLoS One. 2014;9(11):e112373.
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:335–6.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.
Friendly M. Corrgrams. Am Stat. 2002;56:316–24.
Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes J, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.
Ferrocino I, Greppi A, La Storia A, Rantsiou K, Ercolini D, Cocolin L. Impact of nisin-activated packaging on microbiota of beef burgers during storage. Appl Environ Microbiol. 2016;82:549–59.
Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ. GAGE: Generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009;10:161.
Kopitar AN, Ihan Hren NIA. Commensal oral bacteria antigens prime human dendritic cells to induce Th1, Th2 or T reg differentiation. Oral Microbiol Immunol. 2006;21:1–5.
Tseng SK, Tsai YL, Li UM, Jeng JH. Radicular cyst with actinomycotic infection in an upper anterior tooth. J Formos Med Assoc. 2009;108:808–13.
Wayman BE, Murata SM, Almeida RJ, Fowler CB. A bacteriological and histological evaluation of 58 periapical lesions. J Endod. 1992;18:152–5.
Tek M, Metin M, Sener I, Bereket C, Tokac M, Kazancioglu HO, et al. The predominant bacteria isolated from radicular cysts. Head Face Med. 2013;9:25.
Gomes NR, Diniz MG, Pereira TD, Estrela C, de Macedo Farias L, de BAB A, et al. Actinomyces israelii in radicular cysts: A molecular study. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2017;123:586–90.
Shibasaki M, Iwai T, Chikumaru H, Mitsudo K, Inayama Y, Tohnai I. Actinomyces-associated calcifications in a Dentigerous cyst of the mandible. J Craniofac Surg. 2013;24:e311–4.
Kaplan I, Anavi K, Anavi Y, Calderon S, Schwartz-Arad D, Teicher S, et al. The clinical spectrum of Actinomyces-associated lesions of the oral mucosa and jawbones: Correlations with histomorphometric analysis. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2009;108:738–46.
Iatrou IA, Legakis N, Ioannidou E, Patrikiou A. Anaerobic bacteria in jaw cysts. Br J Oral Maxillofac Surg. 1988;26:62–9.
Nair PNR, Sundqvist G, Sjögren U. Experimental evidence supports the abscess theory of development of radicular cysts. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2008;106:294–303.
Zehnder M, Belibasakis GN. On the dynamics of root canal infections -what we understand and what we don’t. Virulence. 2015;6(3):216–22.
Ordiz MI, May TD, Mihindukulasuriya K, Martin J, Crowley J, Tarr PI, et al. The effect of dietary resistant starch type 2 on the microbiota and markers of gut inflammation in rural Malawi children. Microbiome. 2015;3:37.
Meghji S, Qureshi W, Henderson B, Harris M. The role of endotoxin and cytokines in the pathogenesis of odontogenic cysts. Arch Oral Biol. 1996;41:523–31.
The Authors wish to heartily thank collectively all the members of Department of Oral Surgery.
The research was funded with Ricerca Locale 2015 (PI F.M.) of the University of Turin. No involvement of the funding body ever occurred in any phase of the study.
Availability of data and materials
The datasets generated and analysed during the current study are available in the NCBI repository [https://www.ncbi.nlm.nih.gov/SRA/SRP096711].
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the Dental School, University of Turin. The cystic formations removed during surgery were subjected to histology as part of the standard care. Written informed consent to participate was obtained from all study 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
Mussano, F., Ferrocino, I., Gavrilova, N. et al. Apical periodontitis: preliminary assessment of microbiota by 16S rRNA high throughput amplicon target sequencing. BMC Oral Health 18, 55 (2018). https://doi.org/10.1186/s12903-018-0520-8