Abstract
Systemic metabolic signatures of oral diseases have been rarely investigated, and prospective studies do not exist. We analyzed whether signs of current or past infectious/inflammatory oral diseases are associated with circulating metabolites. Two study populations were included: the population-based Health-2000 (
Introduction
Metabolomics refers to the large-scale study of small molecules called metabolites, which collectively form the metabolic phenotype of an individual. Diet, clinical measurements such as blood pressure, body mass index, use of medications, and the composition of gut microbiome predicted 76% of the profiled serum metabolites (Bar et al. 2020). Considering the importance of oral health in cardiovascular diseases (CVDs), its contribution to metabolic profiles in serum is a relevant target for investigation. Although the mechanisms behind the association of oral inflammations and CVD are partially unknown, systemic inflammation and endotoxemia are among the most studied (Pietiäinen et al. 2018), and they both induce multiple metabolic alterations (Khovidhunkit et al. 2004; Määttä et al. 2021; Jiménez et al. 2022).
A widely used high-throughput nuclear magnetic resonance (NMR) platform provides a broad range of clinically and analytically validated metabolic biomarker concentrations. The biomarkers represent cell metabolism and pathways, such as fatty acid balance, glycolysis, ketone bodies, inflammation, lipoprotein class/subclass profiles, and other lipids (Würtz et al. 2017). They have been used extensively to predict disease events (Würtz et al. 2015) and all-cause morbidity (Fischer et al. 2014), measure chronic inflammation (Ritchie et al. 2015), and identify metabolic risk factors, such as insulin resistance (Wang et al. 2012).
The NMR platform provides detailed lipoprotein subclass profiling. All classes—namely, very low-density lipoprotein (VLDL), intermediate-density lipoprotein (IDL), low-density lipoprotein (LDL), and high-density lipoprotein (HDL)—are divided into subclasses from XXL to XS based on particle size. The compositions comprise the particle core containing cholesterol esters and triglycerides, as well as the surface containing free cholesterol, phospholipids, and apolipoproteins. Clinical laboratories traditionally produce only total cholesterol and lipid determinations. For example, plasma triglyceride or cholesterol concentration refers to the sum of these lipids, making no differentiation between carrier lipoprotein classes, sizes, or overall composition. Detailed metabolic profiling extends the knowledge of disease mechanisms and aids to identify modifiable risk factors. For example, trials investigating predictors of future CVD events imply that inflammation is one of the most important sources of the “residual cardiovascular risk” requiring adjunctive treatments (Ridker et al. 2023).
There are only few previous studies concerning systemic metabolic profiling of oral inflammations (Barnes et al. 2014; Chen et al. 2018), and further investigations using untargeted, modern methods are warranted. In the cross-sectional part of our study, we analyzed whether signs of current or past infectious/inflammatory oral diseases are associated with circulating metabolites, including concentrations, compositions, and subclass distributions of lipoprotein particles, as well as levels of fatty acids (FAs), amino acids, ketone bodies, metabolic substrates, and various other metabolic markers determined by NMR spectroscopy-based methods in 2 cohorts. In the longitudinal part of our study, we analyzed whether signs of oral diseases at baseline are associated with metabolic changes during 11 y of follow-up.
Materials and Methods
Study Populations
The Parogene study and the Health-2000 survey are described in detail in the Appendix. The oral parameters registered included radiographically determined caries lesions, periapical infections, the number and quality of root canal fillings (RCFs), and alveolar bone loss (ABL), as well as clinically determined number of teeth with probing pocket depths (PPDs) ≥4 mm and bleeding on probing (BOP). The periodontal inflammatory burden index (PIBI) was calculated by adding the number of teeth with PPD from 4 to 5 mm to the number of teeth with PPD ≥6 mm multiplied by 2.
NMR Metabolomics
Serum samples of the participants were analyzed with a high-throughput serum NMR metabolomics platform as described previously (Würtz et al. 2017) by Nightingale Health. Concentrations of 157 metabolites were available for 6,229 and 452 participants of the Health-2000 survey and Parogene, respectively. In addition, 4,116 participants of the Health-2000 gave follow-up serum samples for metabolic profiling 11 y after the baseline examination.
Statistical Analyses
All statistical analyses were performed with R (v3.6.1 or higher, http://www.r-project.org/).
Metabolite concentrations were log(x + 1) transformed to obtain approximately normal distributions. Thereafter, the measures were scaled to a mean of zero and a standard deviation of 1. We fitted linear regression models between metabolite measures and oral parameters describing periodontal and endodontic conditions in both study populations. The results from the 2 cohorts were combined by a fixed-effect meta-analysis.
For follow-up data of Health-2000, we fitted linear regression models with metabolite concentrations at follow-up as the outcome and oral parameters at baseline as predictors. Both cross-sectional and longitudinal analyses were adjusted for number of teeth, age, sex, smoking (current/former/never), diabetes, and the level of education (only Health-2000), and longitudinal analyses were additionally adjusted for baseline metabolite concentration. Due to the correlated nature of data, we conducted a principal component analysis to define the threshold for statistical significance. The first 18 principal components explained more than 95% of the variation of the metabolomics data. Thus, the
Results
The main characteristics of the study participants are presented in Table 1. The mean age of Parogene patients was higher than that of Health-2000 participants (63.5 vs. 52.8 y), and Parogene included fewer female participants (34.5 vs. 54.6%). As expected, diabetes, hypertension, and CVD were more often present in Parogene patients. Decayed teeth (mean 0.98 vs. 0.45) and RCF (mean 2.18 vs. 1.50), ABL (76.4 vs. 29.5%), and pathologically deepened periodontal pockets (mean 7.16 vs. 4.14) were more common in the Parogene than the Health-2000 population.
Characteristics of the Participants.
BMI, body mass index; BOP, bleeding on probing; PPD, probing pocket depth; RCF, root canal filling.
Having medication for hypertension (Parogene), self-reported (Health-2000).
History of a cardiovascular disease (Health-2000), ≥50% stenosis in coronary angiography (Parogene).
BOP was recorded as the percentage of bleeding sites in Parogene and as the number of sextants with BOP in Health-2000.
Associations of 157 metabolic measures with oral health parameters were analyzed with linear regression models (Appendix Tables 1 and 2). The associations between oral parameters and metabolite concentrations from the meta-analyses are illustrated in Figure 1. Overall, we found 13 significant associations between oral parameters and metabolic measures in Parogene, while the number of significant associations was 168 in Health-2000. In the meta-analyses, the number of teeth with PPD ≥4 mm displayed the highest number of significant associations (93) followed by PPD ≥4 mm combined with BOP (88), PIBI (77), the number of teeth with RCF (47), inadequate RCF (27), and the number of decayed teeth (8). ABL was associated with concentrations of 2 metabolites, and apical rarefactions were associated with the concentration of 1 metabolite.

Strengths of associations between oral parameters and metabolite measures. From Parogene and Health-2000 cohorts, concentrations of 157 metabolites were determined by a nuclear magnetic resonance platform. Their associations with parameters from the oral examination were investigated by using linear regression analyses. The analyses were adjusted for age, gender, smoking, number of teeth, diabetes, and the level of education (Health-2000). The results from the 2 cohorts were combined in a meta-analysis, except for BOP, which could not be harmonized between the cohorts. The estimates are shown color coded and the statistically significant associations (
Associations between caries and endodontic parameters and selected metabolic measures are presented in Figure 2. The number of decayed teeth presented positive associations with LDL diameter and the concentrations of pyruvate and glycine. Negative associations were found between caries and the unsaturation degree of FAs (UNSAT), the proportion of docosahexaenoic acid (DHA%), and omega-3 FAs (omega-3%). These negative associations displayed substantial heterogeneity in the meta-analysis (

Associations between endodontic parameters and metabolite concentrations. From Parogene and Health-2000 cohorts, concentrations of 157 metabolites were determined by a nuclear magnetic resonance platform. Their associations with parameters from the oral examination were investigated by using linear regression analyses. (
Associations between periodontal parameters and selected metabolic measures are presented in Figure 3. BOP was positively associated with several lipoprotein-related parameters, including total TG concentration and VLDL and LDL triglycerides, the number of IDL particles, and the number of medium and small LDL particles. PPD ≥4 mm was associated with total cholesterol, VLDL and LDL cholesterol, TG, and the numbers of VLDL, IDL, and LDL particles; total FA concentration; MUFA%, pyruvate; citrate; isoleucine; leucine; valine; phenylalanine; tyrosine; and GlycA. In addition, PPD ≥4 mm with BOP was associated positively with apoB/apoA-I ratio and negatively with HDL diameter, HDL cholesterol, FA unsaturation degree, omega-6%, and proportion of polyunsaturated FAs (PUFA%).

Associations between periodontal parameters and metabolite concentrations. From Parogene and Health-2000 cohorts, concentrations of 157 metabolites were determined by a nuclear magnetic resonance platform. Their associations with parameters from the oral examination were investigated by using linear regression analyses. (
Participants of the Health-2000 survey (
We analyzed whether the oral health parameters at baseline were associated with changes of metabolic measures during follow-up. Number of teeth was a strong determinant of estimate levels, and thus we analyzed the results separately for those with ≥20 teeth or <20 teeth at baseline (Appendix Tables 4 and 5). Number of significant associations (≥20 teeth/<20 teeth) were for caries 30/11, for RCF 3/0, for inadequate RCF 5/0, for apical rarefactions 3/0, for PPD ≥4 mm 1/2, for BOP 8/0, for PPD ≥4 mm + BOP 1/0, for PIBI 1/2, and for ABL 0/0. Significant associations between caries, BOP, and metabolites are presented for subjects with ≥20 teeth at baseline in Table 2. In the follow-up, baseline caries was associated with changes of IDL lipids, phospholipids, free cholesterol, non-TG lipids in XL- and L-HDL, HDL3 cholesterol, apoA-I, and several FA levels. Mean changes of VLDL-, IDL-, and LDL-related parameters were positive in those with fewer decayed teeth at baseline and negative in those with more decayed teeth (Appendix Table 6). The number of XL-HDL particles and XL-HDL lipid concentration, serum total cholesterol and HDL3 cholesterol, PUFA, and omega-3 increased in individuals with fewer decayed teeth and decreased in those with more decayed teeth.
Longitudinal Associations of Oral Health Parameters with Changes of Metabolic Measures in the Follow-up in Participants with Functional Dentition (Number of Teeth ≥20) at Baseline.
BOP, bleeding on probing; DHA, docosahexaenoic acid; HDL, high-density lipoprotein; IDL, intermediate-density lipoprotein; LDL, low-density lipoprotein; MUFA, monounsaturated fatty acid; NS, nonsignificant; PUFA, polyunsaturated fatty acid; VLDL, very low-density lipoprotein.
BOP at baseline was associated with changes in VLDL particle size; composition of L-HDL particles; serum cholesterol including HDL, HDL2, and HDL3 cholesterol; apoA-I; and omega-3% in follow-up. The increase in VLDL particle size during follow-up was larger in the subgroup with BOP (Appendix Table 6).
Discussion
In this study including 2 populations, a follow-up, and altogether 6,681 participants, we show that signs of current or past infectious/inflammatory oral diseases were associated with a large number of unbeneficial circulating metabolic measures. We found that various signs of oral diseases were associated with rises in VLDL-related parameters, including particle number, size, and neutral lipids. In addition, we found marked associations between oral diseases and increased saturation degree of FAs, lower levels of omega-3 FAs, and increased levels of glycolysis-related metabolites. These metabolic profiles, which were independent of age, sex, number of teeth, diabetes, education, or smoking status, are inflammatory in nature. Furthermore, periodontal inflammation was associated with GlycA, a novel biomarker of systemic inflammation.
Systemic inflammation is an essential feature in the onset, maintenance, and exacerbation of cardiometabolic disorders. So far, metabolomics studies on the contribution of oral diseases to metabolic profiles are scarce, and prospective studies do not exist. Thus, our study is one of those in the frontline. Caries and periodontitis are putatively associated with systemic inflammation through overlapping mechanisms, although the role of caries in systemic health has been far less investigated. Caries reflects a high-carbohydrate/sugar-rich diet, which induces receptors of advanced glycation end products, promoting oxidative stress and inflammation (Sabharwal et al. 2021). Local inflammation in the oral cavity may lead to the presence of proinflammatory mediators in circulation promoting systemic inflammation (Sharma et al. 2017; Garrido et al. 2019). An additional mechanism resulting in systemic inflammation in periodontitis is the translocation of bacteria from the oral cavity directly to blood or down via the gastrointestinal tract (Kitamoto and Kamada 2022). Our results suggest that metabolic profiles in individuals with caries or periodontal inflammation resemble each other, but caries/endodontic parameters are associated with fewer metabolic markers than periodontal parameters, suggesting smaller systemic impact.
The main lipoprotein class whose compositional alterations were associated with periodontal inflammation and the number of RCF was apoB-containing, proatherogenic lipoproteins. RCFs were directly linked to VLDL diameter and the number of VLDL particles of all sizes, whereas BOP and PPD ≥4 mm were positively associated with VLDL cholesterol, apoB, and the number of smaller apoB-containing particles, such as S-VLDL, IDL, and L-LDL. The NMR platform categorizes VLDL into 6 subclasses (Wang et al. 2012); the diameter of XXL- and XL-particles exceeds 64 nm, whereas the S-sized have an average size of 37 nm. Although XX- or XL-VLDL cannot enter the intima of the vessel walls, they have been associated with cardiometabolic diseases (Lee et al. 2022). Large VLDL particle size has been associated with carbohydrate-rich diet (Feingold and Grunfeld 2022), which may explain the observed associations between RCF and VLDL parameters. The smaller VLDL, IDL, and LDL, however, can more readily cross the endothelium (Borén et al. 2000). During inflammation, activity of lipoprotein lipase decreases maintaining large VLDL size, whereas lipolysis of adipose tissue provides FAs for increased hepatic triglyceride synthesis (Khovidhunkit et al. 2004). Thus, an increased amount of triglycerides is available for formation and hepatic secretion of VLDL (Feingold and Grunfeld 2022). Our observations indicate disturbances in both VLDL production and clearance in oral diseases.
In addition to an overall decrease of HDL, inflammation leads to structural changes diminishing the antiatherogenic character of HDL (Khovidhunkit et al. 2004): HDL remodeling produces large amounts of small HDL due to replacement of apoA-I with serum amyloid A, and the particles become enriched with triglycerides. Earlier publications suggest that childhood caries and cumulative signs of oral inflammations are associated with low HDL cholesterol levels in adulthood (Pussinen et al. 2019, 2020). In addition, adult periodontitis patients have lower HDL cholesterol levels than those without periodontitis (Buhlin et al. 2003), and HDL in periodontitis patients displays decreased ability to function in reverse cholesterol transport (Pussinen et al. 2004). Also in the present study, the amount of HDL-TG increased in the presence of periodontal inflammation and in individuals with RCFs. In addition, HDL particle size was smaller and the number of large HDL particles lower in individuals with RCFs. Thereby, oral diseases were associated with HDL characteristics resembling those found in acute-phase HDL (Khovidhunkit et al. 2004).
ApoB/apoA-I ratio was positively associated with RCFs and PPD ≥4 mm. ApoB/apoA-I ratio, which constitutes from the major proteins of LDL and HDL, is regarded as an indicator of the balance of proatherogenic and antiatherogenic lipoproteins as well as inflammatory status (Wu et al. 2019). Markers of periodontal inflammation were also positively associated with a novel inflammatory biomarker, GlycA. The GlycA signal arises largely from the glycan portions of acute-phase proteins. Its levels are associated with interleukin-6, C-reactive protein, and serum amyloid A (Ritchie et al. 2015; Connelly et al. 2017), and it reflects chronic inflammation, cardiometabolic risk, and poor survival (Ritchie et al. 2015; Julkunen et al. 2023).
Unlike other periodontal parameters, ABL had mostly nonsignificant associations with metabolites. The development of ABL is usually a slow process that takes place over decades. Conventional radiographs demonstrate bone loss only after a significant amount of destruction (Zaki et al. 2015). Thus, panoramic radiographs are not optimal for detecting current disease activity, which might explain the weak association with metabolic measures. Instead, periodontal pockets are a favorable growth site for dysbiotic biofilm, which initiates and maintains the host immune-inflammatory response. It is clinically marked by BOP, which still is a gold standard marker of active periodontal inflammation (Bartold and Van Dyke 2017).
Several parameters related to FA synthesis/turnover were associated with caries, endodontic lesions, RCFs, and periodontal inflammation. BOP was associated with increased levels of saturated FAs, as well as PPD ≥4 mm and caries with a clear decrease in unsaturated FAs. The omega-3 PUFAs decreased in the presence of caries and periapical inflammations, and omega-6 PUFAs decreased in individuals with BOP and PPD ≥4 mm. Several pathways are likely to contribute to these findings: first, an oxidative response to inflammation could induce peroxidation/modification of these highly unsaturated FAs (Khovidhunkit et al. 2004). Second, the omega-3 and omega-6 FA pool is fueled into proresolving and proinflammatory lipid mediator synthesis in oral diseases (Sommakia and Baker 2017). In addition, the FA composition of serum reflects diet, which obviously has a major effect on oral health. In the recent atlas of plasma NMR metabolites, MUFA was the biomarker that was associated with the highest number of disease end points and displayed similar disease clustering as GlycA (Julkunen et al. 2023). The authors suggested that several FA-related biomarkers should be considered markers of systemic inflammation rather than markers of recent diet (Julkunen et al. 2023).
Infectious/inflammatory oral diseases were also associated with non-lipid-related metabolic measures. Pyruvate concentrations were elevated in patients with caries, endodontic lesions, and periodontal inflammation. Our findings could be indicative of increased aerobic glycolysis and disrupted carbon flux through the citric acid cycle similarly as in endotoxemia (Määttä et al. 2021). These metabolites may also reflect higher sugar content in diet (Wang et al. 2012). Amino acids, which have been consistently associated earlier with periodontitis in different biofluids (Brito et al. 2022), were in the present study upregulated in individuals with caries, inadequate RCF, and signs of periodontal inflammation.
In addition to cross-sectional associations, our observations indicate that infectious/inflammatory oral diseases may predict future metabolic profiles. In the follow-up of the Health-2000 population, baseline caries and BOP were associated with decreases of L-HDL particle concentration, their non-TG lipid content, and apoA-I, but associations with proatherogenic apoB-containing lipoproteins were also observed. It is known that serum cholesterol concentration decreases during inflammation due to its reduction in both HDL and LDL (Khovidhunkit et al. 2004). These associations were more obvious in participants with functional dentition, since the number of teeth was a strong confounder. Although no information on the number of occluding pairs was available, 20 teeth was selected as a cutoff for “adequate dentition” because functional dentition comprises ≥10 occluding pairs (Ikebe et al. 2007). Any professional intervention, such as the number of RCFs (Liljestrand et al. 2021), may be a sign of better health-consciousness and the utilization of health care services. This may be a more relevant predictor of future metabolic profiles than oral health. Unfortunately, we did not have information on potential periodontal/restorative/endodontic treatment during the follow-up, which warrants for further research.
Strengths of the study include the sample including both a population-based study with adults of all ages and a cohort of elderly patients with cardiac symptoms. The study populations were obviously different regarding the presence of CVD or their risk factors, such as hypertension, and some parameters were collected from the populations using different methods. However, virtually all significant associations observed between oral diseases and metabolic measures were similar in both study populations, strengthening the credibility of the findings. Limitations of the study include the exclusion of several potentially important determinants of the metabolic profiles: we only adjusted the analyses for diabetes, but other cardiometabolic diseases, such as CVD or hypertension, and their medications may also have an impact. In addition, reverse causality may introduce a bias in the analyses, as the results from the cross-sectional part of our study were stronger than those from the follow-up part. The contribution of unbeneficial metabolic profiles on future oral disease risk remains unknown. Diet may be reflected in virtually all metabolites investigated, which was clearly shown in the short term when comparing postprandial and fasting samples (Bermingham et al. 2023): 85% of the metabolites of the present panel changed significantly during fasting. In the present study, both caries and BOP were associated with several FA measures, such as saturation degree, suggesting an unhealthy diet. Thus, consumption of fat, sugar, vegetables, fibers, and alcohol would have been of interest but represents a topic for a separate study. Overall, the metabolic measures are affected by several lifestyle choices, diseases, medications, and even genetics (Nath et al. 2017). Despite these flaws, our study is the first to investigate the association between serum metabolomics and oral health with a large sample including 2 study populations, both caries- and periodontitis-related parameters, and with follow-up data.
We conclude that both caries- and periodontitis-related signs of oral diseases are associated with unbeneficial metabolic profiles with an inflammatory nature. Due to the high prevalence of oral diseases in most populations, their role in the total inflammatory burden is significant, and they may lead to deteriorating public health. Infectious/inflammatory oral diseases may represent a modifiable risk factor for systemic chronic inflammation and thus cardiometabolic disorders.
Author Contributions
A. Salminen, contributed to design, data acquisition, analysis, and interpretation, drafted the manuscript; A.M. Määttä, J. Leskelä, contributed to data analysis and interpretation, critically revised the manuscript; P. Mäntylä, A.L. Suominen, contributed to design, data acquisition and interpretation, critically revised the manuscript; M. Pietiäinen, contributed to conception, data interpretation, critically revised the manuscript; K. Buhlin, S. Paju, contributed to data acquisition and interpretation, critically revised the manuscript; W. Sattler, contributed to data interpretation, critically revised the manuscript; J. Sinisalo, contributed to conception, design, data acquisition and interpretation, critically revised the manuscript; P.J. Pussinen, contributed to conception, design, data acquisition and interpretation, drafted the manuscript. All authors gave final approval and agree to be accountable for all aspects of the work.
Supplemental Material
sj-docx-1-jdr-10.1177_00220345231203562 – Supplemental material for Systemic Metabolic Signatures of Oral Diseases
Supplemental material, sj-docx-1-jdr-10.1177_00220345231203562 for Systemic Metabolic Signatures of Oral Diseases by A. Salminen, A.M. Määttä, P. Mäntylä, J. Leskelä, M. Pietiäinen, K. Buhlin, A.L. Suominen, S. Paju, W. Sattler, J. Sinisalo and P.J. Pussinen in Journal of Dental Research
Footnotes
A supplemental appendix to this article is available online.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The Health-2000 survey received support from the Finnish Dental Society (Apollonia), the Finnish Dental Association, and Plandent Oyj. The research group received funding from the Sigrid Juselius foundation and the Academy of Finland (#1340750).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
