Shuvashis
Das Gupta‡
a,
Markus
Killenberger‡
b,
Tarja
Tanner
bd,
Lassi
Rieppo
a,
Simo
Saarakkala
ac,
Jarkko
Heikkilä
b,
Vuokko
Anttonen
bd and
Mikko A. J.
Finnilä
*a
aResearch Unit of Medical Imaging, Physics and Technology, University of Oulu, 90220 Oulu, Finland. E-mail: shuvashis.dasgupta@oulu.fi; lassi.rieppo@oulu.fi; simo.saarakkala@oulu.fi; mikko.finnila@oulu.fi
bResearch Unit of Oral Health Sciences, Department of Cariology, Endodontology and Pediatric Dentistry, University of Oulu, 90220 Oulu, Finland. E-mail: markus.killenberger@student.oulu.fi; tarja.Tanner@oulu.fi; jarkko.heikkila@student.oulu.fi; vuokko.anttonen@oulu.fi
cDepartment of Diagnostic Radiology, Oulu University Hospital, 90220 Oulu, Finland
dMedical Research Center, Oulu University Hospital and University of Oulu, 90220 Oulu, Finland
First published on 18th November 2020
Dental caries is the most common oral disease that causes demineralization of the enamel and later of the dentin. Depth-wise assessment of the demineralization process could be used to help in treatment planning. In this study, we aimed to provide baseline information for the development of a Raman probe by characterizing the mineral composition of the dental tissues from large composition maps (6 × 3 mm2 with 15 μm step size) using Raman microspectroscopy. Ten human wisdom teeth with different stages of dental caries lesions were examined. All of the teeth were cut in half at representative locations of the caries lesions and then imaged with a Raman imaging microscope. The pre-processed spectral maps were combined into a single data matrix, and the spectra of the enamel, dentin, and caries were identified by K-means cluster analysis. Our results showed that unsupervised identification of dental caries is possible with the K-means clustering. The compositional analysis revealed that the carious lesions are less mineralized than the healthy enamel, and when the lesions extend into the dentin, they are even less mineralized. Furthermore, there were more carbonate imperfections in the mineral crystal lattice of the caries tissues than in healthy tissues. Interestingly, we observed gradients in the sound enamel showing higher mineralization and greater mineral crystal perfection towards the tooth surface. To conclude, our results provide a baseline for the methodological development aimed at clinical diagnostics for the early detection of active caries lesions.
The traditional treatment for cavitated caries lesions is to remove the infected enamel and dentin and place a filling. However, the prognosis is substantially better if caries lesions are detected at an early stage and their progression can be halted non-invasively.5 Caries is traditionally detected with a visual and tactile inspection aided by radiography (X-ray).6 Unfortunately, this system has a relatively low sensitivity yet high specificity.7 Conventional radiography tends to underestimate caries-induced tissue damage and demineralization. Additionally, it usually only detects lesions that have already extended beyond the enamel–dentin junction.8 Therefore, modern methods to detect early lesions are needed.
The healthy enamel is approximately 95% mineral, 1% organic matter, and 4–5% water by weight.9 The mineral crystals in the enamel are arranged in a tightly packed rod-like pattern that forms the enamel prisms. Caries forms when the crystals begin to dissolve in the acidic environment of the dental biofilm. When the minerals are dissolved, the intercrystalline space increases, and the surface of the enamel becomes softer and more porous.3 The lesion will grow substantially quicker when it reaches the dentin as it is less mineralized. The pulp and dentin react to the stimuli from a caries lesion by sclerosis of the dentinal tubules and the formation of reactive dentin.
Raman spectroscopy uses monochromatic light (mostly in the near-infrared, visible, or UV range) that interacts with chemical bonds within the specimen.10 Most of the light is scattered elastically as Rayleigh scattering having the same wavelength as the original light. Some of the light (low probability: around 1 in 108 photons), however, is scattered inelastically as Raman scattering with a different wavelength due to the change in the polarizability of molecules. The Raman scattering spectrum provides quantitative information of the molecular composition of the specimen. Raman spectroscopy is particularly well suited for studying biomaterials because of its high specificity for biomolecules and low interference with water.10 Moreover, it is a non-destructive technique and does not require molecular labeling.
Akkus et al. used Raman spectroscopy to study healthy enamel from extracted human incisors.11 They reported that the mineral content was lower in the cervical region than in the rest of the crown, and that the overall mineral content had substantial variance between individuals. Carious lesions have also been studied by Raman spectroscopy;12–16 the clinical and subclinical enamel lesions were found to have a lower mineral content in these studies.14,15 Wide-field Raman spectroscopy has also been used to determine regions of hypo- and hyper-mineralized lesions.17 Moreover, Raman polarization anisotropy has been used to detect early carious lesions.18 This method was based on the fact that there is a reduction in the polarization anisotropy of certain peaks in the Raman spectra from the caries lesion.
In our previous proof-of-concept study with a time-resolved complementary-metal–oxide-semiconductor (CMOS) single-photon-avalanche-diode (SPAD) based Raman spectrometer, we have demonstrated that the fluorescence-suppressed Raman spectrometers have the potential to be used in clinical settings.19 Although time-resolved Raman spectroscopy holds great potential, detailed mapping has remained challenging. In continuation of these efforts, we aimed to use Raman microspectroscopy in this in vitro study to map the mineral composition of the dental tissues and caries lesions in a detailed manner. We hypothesized that enamel mineralization is not homogeneous and that a simple multivariate clustering technique can be used to detect even small carious lesions from the Raman maps.
Sample number | ICDAS and activity grading | ICCMS radiographic grading |
---|---|---|
1 | 2+ | 3 |
2 | 3+ | 3 |
3 | 5+ | 4 |
4 | 1+ | 3/2 |
5 | 5+ | 3 |
6 | 2+ | 3 |
7 | 3+ | 2/3 |
8 | 6+ | 3 |
9 | 3+ | 3 |
10 | 2− | 3 |
A 785 nm laser (30 mW) and a 50 μm confocal pinhole aperture were used to excite the Raman scattering. The Raman spectra were collected for 0.025 s and averaged for 5 scans. A wide-range grating (spectral resolution of 5 cm−1 in the range of 50–3250 cm−1) was used for the Raman measurements. A rectangular region-of-interest containing the enamel, dentin, and caries lesions was selected for Raman mapping based on visual inspection of the full mosaic image (dark-field image using an optical microscope) for each sample. The measurement regions were mapped with a step size of 15 μm.
The pre-processed Raman maps of all samples were combined into a single data matrix for the cluster analysis to reduce the sample-specific bias and increase the interpretability of the clusters. The number of clusters was increased iteratively, and after each iteration, a pseudo-colored cluster image was generated and “subjectively” matched with the microscopic images until the tissue-specific clusters matched the respective tissue locations from the microscopic images. After the final iteration, all clusters were annotated using the microscopic images as references. Subsequently, binary masks were generated for the enamel, dentin, and caries tissues by combining the clusters assigned to the same tissue type. Speckles were removed from the binary masks by performing area-opening operations (less than 3 pixels). Finally, the mean cluster spectra of each tissue type were determined.
We used the V1PO43− band (932–980 cm−1) to represent the mineral content.17,21,22 This band is strongly polarization dependent, but as we used a depolarized laser for the excitation of the Raman scattering, we expected that this effect would not influence our measurements. The amide III band (1215–1300 cm−1) (less susceptible to the tissue orientation effects23) was used to represent the organic matrix.
The type-B carbonate symmetric stretch band (1055–1090 cm−1) was normalized by the mineral band to measure the type-B carbonate substitution in the apatite lattice.24,25 Moreover, the full width at half maximum (FWHM) of the V1PO43− band was calculated, as it inversely correlates with the degree of mineral crystallinity (a measure of the mineral crystal size and/or perfection).22
Finally, 2-D chemical maps were constructed by calculating the ratios of the areas under the Raman bands (calculated for the local baseline within the wavenumber range) for each pixel.
The mean spectra of the clusters for the enamel, dentin, and caries are shown in Fig. 3B. The mean spectra of all twenty-four clusters are shown in ESI Fig. 1.† As expected, compared to the mean spectra of the dentin, the mean spectra of the enamel have higher values in the major mineral peaks and lower values in the organic peaks (ESI Fig. 2†). Moreover, the mean cluster spectra of both enamel and dentin have higher values in the V1PO43− band compared to the mean spectra of the caries. Additionally, the raw Raman spectra (average of 40 spectra) of the enamel, dentin, and caries (approximate location) collected from a sample with caries lesions extended into the dentin are shown in Fig. 3A.
The enamel lesions (Fig. 4C) were observed to be less mineralized and they have more substituted carbonate in the crystal lattice compared to the sound enamel (Fig. 4D). When the lesion extends into the dentin (Fig. 4A), the lesion tissues were observed to have even less mineral than the sound dentin, but many mineral apatites of caries were more crystalline and had less substituted carbonate in the lattice.
The chemical maps show that the sound (healthy) enamel tissues have a greater degree of mineralization and mineral crystallinity and less carbonate in the crystal lattice compared to the sound dentin. The mineral composition is not homogeneous within the sound enamel tissue, but more crystallinity (greater mineral crystal size) and greater mineral crystal perfection are observed towards the tooth surface, whereas the mineral apatites are less crystalline and have more carbonate imperfection in the crystal lattice near the enamel–dentin junction. These gradients of the mineral crystal size and imperfection in the crystal lattice from the tooth surface to the junction of the sound enamel tissues are shown in ESI Fig. 3.† As there are many undefined mineral-to-matrix ratio values present in the sound enamel tissue, it is challenging to profile the gradient of the degree of mineralization. However, in Fig. 4D, a higher degree of mineralization can be observed towards the tooth surface.
In Fig. 5, the violin plots show the distributions of the Raman parameters in different tissue types for all samples, and the descriptive statistics (median and median absolute deviation) are presented in Table 2. The Kruskal–Wallis test shows that there is a statistically significant difference between the degree of mineralization of the different tissue types (chi-squared value = 376554, p < 0.001). The enamel has a higher (p < 0.001) degree of mineralization than the dentin and the caries lesions have a lower (p < 0.001) degree of mineralization than the enamel and dentin.
Raman parameters | Tissue | Median | Median absolute deviation |
---|---|---|---|
Degree of mineralization | Enamel | 130 | 127 |
Dentin | 7.32 | 1.81 | |
Caries | 4.02 | 4.02 | |
Mineral crystallinity | Enamel | 0.0751 | 0.005 |
Dentin | 0.0581 | 0.00531 | |
Caries | 0.0697 | 0.016 | |
Type-B carbonate substitution | Enamel | 0.0564 | 0.0113 |
Dentin | 0.0844 | 0.0142 | |
Caries | 0.0965 | 0.0531 |
There is a statistically significant difference also between the mineral crystallinity of the different tissue types (chi-squared value = 277515, p < 0.001). The enamel is more crystalline than the dentin (p < 0.001). The caries lesions have higher crystallinity (p < 0.001) than the dentin but lower crystallinity than the enamel (p < 0.001).
Finally, the differences in the type-B carbonate substitution in the mineral crystal lattice between the different tissue types are statistically significant (chi-squared value = 264836, p < 0.001). The enamel has a less (p < 0.001) substituted lattice compared to the dentin, indicating the presence of mineral crystals with greater stoichiometric perfection in the enamel than in the dentin. The caries lesions have a higher (p < 0.001) carbonate substitution in the lattice than the enamel and dentin, which indicates the presence of crystals with less stoichiometric perfection in the caries lesion.
Previous studies have also used Raman spectroscopy to determine the mineral composition of dental hard tissues.14,15,17,26 Raman composition maps (using a step size of 50 μm) have been used previously to identify enamel demineralization caused by artificial caries using the intensity and FWHM of the V1PO43− band, depolarization ratio, and polarization anisotropy,26 where, based on the chemical analysis, it was concluded that demineralized enamel tissues have a disordered structure at very early stages.
Here, we applied K-means cluster analysis to identify dental caries from the Raman maps. K-Means cluster analysis is an unsupervised technique, but it depends on the user-defined number of clusters. We iteratively increased the number of clusters and visually compared the cluster images with the microscopic images after each iteration. Since we did not have the tissue-specific “true” labels for the microscopic images, objective validation of clustering could not be performed. To increase the interpretability of cluster analysis, we performed the clustering on the combined data matrix. This combined data matrix contains the Raman spectra of the enamel, dentin, caries with different severities, and resin, and other background noises due to aqueous measurement. This might be the reason why as many as twenty-four clusters were needed to segment all the tissue types. Overall, we observed that the cluster analysis was sensitive to the degree of mineralization. For instance, in the sample with severe caries lesions (sample 3, Fig. 2B), we found more than six clusters in the caries region. However, it should be possible to identify the different dental tissues with a smaller number of clusters if the clustering was performed on the Raman map of each sample individually. Furthermore, the limited sample set (10 samples) prevented us from associating specific clusters to caries progression.
The locations of the caries clusters were well correlated with the locations of the visually identified lesions (Fig. 2 and 4). The carious tissue clusters were located in the areas where caries are most typically found: either on the surface of the enamel or inside the enamel or dentin. In samples with severe clinical lesions, caries was also found grossly invading the dentin. Moreover, thin caries clusters were identified on the enamel surface in most of the samples that were assessed clinically as having active caries lesions. We suggest that those caries clusters represent the demineralization happening on the tooth surface. On tooth surfaces, there is a continuous demineralization–remineralization process going on. If the balance turns towards demineralization then sub-clinical and later clinically detectable lesions occur known as white spot lesions. Here, these surface clusters indicate that the oral circumstances of these samples having active caries lesions are favoring the demineralization process.
We found that the enamel carious lesions are less mineralized than the healthy enamel. This local demineralization might be due to the lactic acid produced by oral bacteria. When the demineralization process exceeds the remineralization process by saliva or clinical fluoride repair, the dental caries lesion originates at certain anatomical predilection sites on the teeth.5 We found that when the lesions extend into the dentin, they are even less mineralized than the dentin. Moreover, the values of all the analyzed Raman parameters are widely spread for the caries lesions. One of the reasons for this might be the combining of the values from the caries with different progression stages (surface caries, enamel caries, and dentin caries) in a single group.
In dental tissues, the type-B carbonate substitution happens when CO32− ions (impurity) occupy the PO43− sites of the hydroxyapatite crystal lattice. Carbonated hydroxyapatite is less stable, has lower hardness, and is more acid-soluble than non-carbonated hydroxyapatite.3,27 Here, caries lesions were found to be rich in substituted carbonate, and, therefore, they have stoichiometrically imperfect crystals. Stoichiometrically imperfect crystals are common in the demineralization process. Precipitations with high amounts of carbonate content are associated with caries lesions.28,29 Based on our mineral crystallinity results, these crystals of the caries tissues were thicker than those of the dentin tissues but thinner than the crystals of the enamel tissues (Table 2).
Previously, Fourier transform Raman microspectroscopic mapping was used to examine the distribution of phosphate (960 cm−1), carbonate (1070 cm−1), and C–H stretch bands (2700–2880 cm−1) in a cross-section of a human tooth.30 An increase of carbonate ions from the outside of the enamel towards the enamel–dentin junction was reported, while the phosphate ions showed the opposite gradient. Calcium hydroxyapatite and magnesium phosphate are the main sources of calcium and magnesium, respectively, in tooth enamel. In a recent study,31 an increase of calcium and magnesium contents towards the enamel–dentin junction in the central upper incisor teeth was reported using atomic absorption spectrometry. However, a decrease of the calcium content from the enamel surface towards the dentin has been also suggested.32,33 As magnesium ions act as a mineral growth inhibitor in calcium and phosphorus solutions,31 it is reasonable to assume that the magnesium content would increase with the decrease of the degree of mineralization towards the enamel–dentin junction.
This study used a relatively small laser step size (15 μm) for the measurements, which gave a detailed spatial picture of the chemical properties of hydroxyapatite and allowed us to visualize the mineral composition gradient in the sound enamel. We observed gradients of decreasing degree of mineralization, mineral crystallinity, and lattice perfection (Fig. 4 and ESI Fig. 3†) from the occlusal surface towards the enamel–dentin junction. It could be that the dentin integrates into the enamel more coronally than towards the enamel–dentin junction or that the enamel contains more impurities as it approaches the dentin. Since the previous literature on enamel mineralization is contradictory, our results provide the much-needed validity to these observations.
In our proof-of-concept study, we imaged human teeth using a time-resolved CMOS SPAD-based Raman spectrometer and analyzed the Raman chemical image maps (step size: 250 μm).19 A fluorescence suppressed time-resolved Raman spectrometer can be used to rapidly acquire comparable spectra with continuous wave Raman spectroscopy. Furthermore, a time-resolved CMOS SPAD-based Raman spectrometer can be used to collect depth-resolved chemical information.34 Our current study shows a natural decline in enamel mineralization which is less severe as in caries-associated demineralization. In continuation of the proof-of-concept study, the findings of this study provide baseline characteristics of caries lesions and depth-wise mineralization to enable the methodological development using a time-resolved Raman spectrometer aimed at clinical diagnostics for the early detection of active caries lesions.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0an01938k |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2021 |