Next Article in Journal
Involvement of Il-33 in the Pathogenesis and Prognosis of Major Respiratory Viral Infections: Future Perspectives for Personalized Therapy
Next Article in Special Issue
Gelsolin as a Potential Biomarker for Endoscopic Activity and Mucosal Healing in Ulcerative Colitis
Previous Article in Journal
Natural Compounds Targeting Cancer-Associated Fibroblasts against Digestive System Tumor Progression: Therapeutic Insights
Previous Article in Special Issue
Frequency and Clinical Significance of Elevated IgG4 in Rheumatoid Arthritis: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Periodontal Disease Augments Cardiovascular Disease Risk Biomarkers in Rheumatoid Arthritis

1
Section of Oral Health and Periodontology, Division of Oral Diseases, Department of Dental Medicine, Karolinska Institutet, 14104 Huddinge, Sweden
2
Department of Microbiology, Faculty of Life Sciences and Informatics, Balochistan University of Information Technology, Engineering and Management Sciences (BUITEMS), Quetta 87300, Pakistan
3
Department of Applied Oral Sciences, The Forsyth Institute, Cambridge, MA 02142, USA
4
Rheumatology Clinic, Habib Medical Centre, Karachi 75950, Pakistan
5
Department of Periodontology, Altamash Institute of Dental Medicine, Karachi 75500, Pakistan
6
Clinical Chemistry, Department of Medical Sciences, Uppsala University, 75236 Uppsala, Sweden
7
Department of Oral Medicine, Infection and Immunity, Harvard Faculty of Medicine, Boston, MA 02115, USA
*
Author to whom correspondence should be addressed.
Biomedicines 2022, 10(3), 714; https://doi.org/10.3390/biomedicines10030714
Submission received: 15 February 2022 / Revised: 15 March 2022 / Accepted: 16 March 2022 / Published: 19 March 2022

Abstract

:
Objectives: Periodontal disease (PD) and rheumatoid arthritis (RA) are known chronic conditions with sustained inflammation leading to osteolysis. Cardiovascular diseases (CVD) are frequent comorbidities that may arise from sustained inflammation associated with both PD and RA. In order to determine CVD risk, alterations at the molecular level need to be identified. The objective of this study, therefore, was to assess the relationship of CVD associated biomarkers in RA patients and how it is influenced by PD. Methods: The study consisted of patient (26 RA with PD, 21 RA without PD, 51 patients with PD only) and systemically and periodontally healthy control (n = 20) groups. Periodontal parameters bleeding on probing, probing pocket depth, and marginal bone loss were determined to characterize the patient groups. Proteomic analysis of 92 CVD-related protein biomarkers was performed using a multiplex proximity extension assay. Biomarkers were clustered using the search tool for retrieval of interacting genes (STRING) to determine protein–protein interaction (PPI) networks. Results: RA patients with PD had higher detection levels for 47% of the measured markers (ANGPT1, BOC, CCL17, CCL3, CD4, CD84, CTRC, FGF-21, FGF-23, GLO1, HAOX1, HB-EGF, hOSCAR, HSP 27, IL16, IL-17D, IL18, IL-27, IL6, LEP, LPL, MERTK, MMP12, MMP7, NEMO, PAPPA, PAR-1, PARP-1, PD-L2, PGF, PIgR, PRELP, RAGE, SCF, SLAMF7, SRC, THBS2, THPO, TNFRSF13B, TRAIL-R2, VEGFD, VSIG2, and XCL1) as compared to RA without PD. Furthermore, a strong biological network was identified amongst these proteins (clustering coefficient = 0.52, PPI enrichment p-value < 0.0001). Coefficients for protein clusters involved in CVD (0.59), metabolic (0.53), and skeletal (0.51) diseases were strongest in the PD group. Conclusion: Periodontal disease augments CVD-related biomarkers in RA through shared pathological clusters, concurrently enhancing metabolic and skeletal disease protein interactions, independent of autoimmune status.

1. Introduction

Chronic inflammation stems from persistent acute inflammation due to the failure to resolve the acute phase, often associated with the inability to remove the inducing agent or stimulus. Several diseases that acquire such chronicity due to a dysregulated immune response include atherosclerosis, type 2 diabetes, rheumatoid arthritis (RA), and periodontal disease (PD) [1].
Cardiovascular diseases (CVD) are the leading cause of global mortality with over 75% of cases in low- and middle-income countries [2]. CVD comprises a group of disorders that involve the heart muscle and blood vessels. The most common pathogenic pathway that leads to CVD is atherosclerosis [3]. Risk factors such as smoking, diabetes, hypertension, and obesity are transduced into atherosclerotic events via complex interactions between endothelial adhesion molecules and inflammatory cells including macrophages and T lymphocytes. The inflammatory response also has an autoimmune component as low-density lipoprotein (LDL) cholesterol, one of the retained lipids in atherosclerotic plaques, is antigenic resulting in production of high affinity antibodies [4].
PD is an independent risk factor for the development of CVD [5]. Systematic reviews have shown a consistent association between CVD and PD which may be partially attributed to shared risk factors and the dissemination of periodontal pathogens into the bloodstream or an increase in systemic inflammation [6].
RA is an autoimmune disease characterized by synovial inflammation and destruction due to immune mediated inflammation. This sustained inflammation in RA promotes cardiovascular pathology to such an extent that it remains the leading cause of mortality in RA patients [7]. The overall increased CVD risk in RA has been attributed less to traditional CVD risk factors and more to underlying autoimmunity and inflammatory burden [8].
The influence of PD and RA combined may increase the menacing effects of inflammation and raise an individual’s risk of developing CVD even further. This can be evaluated through exploration of emergent biomarkers involved in CVD initiation and pathology and, therefore, the aim of this study is to assess CVD related biomarkers in RA patients and how it is influenced by PD.

2. Materials and Methods

The study was performed at the Department of Periodontology, the Altamash Institute of Dental Medicine, between October 2012 and August 2017 in Karachi, Pakistan. Upon obtaining informed consent, a detailed questionnaire was used to acquire information pertaining to medical and dental history. The minimum sample size was calculated using the Epitools Epidemiological Calculators [9] with the assumptions of a power of 80% and a confidence interval (CI) of 95%. The parameters reported in the literature pertaining to South Asian population were used [10]: (1) frequency of 60% for PD among RA patients and (2) frequency of 28% for PD without RA.

2.1. Rheumatoid Arthritis Group

A total of 47 RA patients were recruited via consecutive sampling from the Department of Rheumatology at the Habib Medical Centre in Karachi, Pakistan. These patients were established RA cases diagnosed by a rheumatologist (AG) using current ACR/EULAR classification criteria [11]. All patients were receiving disease-modifying anti-rheumatic drugs (DMARDs), corticosteroids, non-steroidal anti-inflammatory drugs (NSAIDs), or a biologic DMARD (Rituximab) at the time of the examination. Based on their periodontal status, the patients were divided into two groups: RA with PD (n = 26) and RA without PD (n = 21).

2.2. Periodontal Disease Group

A group of 51 participants diagnosed with PD, but exhibiting no signs of RA, gout, or osteoarthritis were also included. Individuals with a history of treatment for PD during the last six months and/or treatment with antibiotics in the last three months were excluded.
We used twenty controls for comparison and better characterization of the groups. All controls had clinically healthy periodontium and no systemic disease. Blood samples were drawn from all participants and prepared sera were stored at −80 °C until the time of analyses.

2.3. Periodontal Examination

Periodontal examination was carried out for all teeth except for the third molars by a single examiner (JP). The defining criteria for PD was probing pocket depth (PPD) of ≥5 mm (to the nearest millimeter) in at least three different sites using a periodontal probe (Hu-Friedy manufacturing, Chicago, IL, USA). Pockets measuring ≥5 mm were added to calculate the sum of deep pockets representing PD affected sites. Our group has designed continuous periodontal indices to gauge the severity of PD as a continuous variable rather than dichotomous. These parameters were used to assess the severity of PD:
  • Bleeding on probing (BOP)
  • ΣPPD Total
  • ΣPPD Disease
  • Adjusted PPD Total
  • Adjusted PPD Diseased sites
  • Σ Marginal bone loss (MBL)
  • Adjusted ΣMBL
Details of the parameters and their measurements are described in our previous publications [12,13].

2.4. Anthropometric Measures

Body weight was measured to the nearest kg. Using non-stretchable measuring tape, height and waist were measured to the nearest cm. Waist was measured in the horizontal plane at the midpoint between the lowest rib and the iliac crest. Body mass index (BMI) was calculated from weight and height measurements (kg/m2). Anthropometric measurements were recorded for all four groups.

2.5. Glycated Hemoglobin (HbA1c)

Glycated hemoglobin levels were determined for all four groups after collecting four milliliters of whole blood into spray-coated EDTA tubes (lavender top, Becton, Dickinson, Franklin Lakes, NJ, USA). The samples were analyzed on the same day using the ion-exchange high-performance liquid chromatography system Bio-Rad D-10 Hemoglobin Testing System (Bio-Rad Laboratories, Hercules, CA, USA). The HbA1c values are standardized according to the National Glycohemoglobin Standardization Program (NGSP) system [14].

2.6. Proteomic Profiling

Serum samples were analyzed at the SciLifeLab Affinity Proteomics Uppsala, Uppsala University, (Uppsala, Sweden) using proximity extension assay (PEA) technology (Olink Proteomics, Uppsala, Sweden). Levels of 92 proteins from the Olink Target96 CVD II panel were measured (Supplementary File Table S1). The PEA technology utilize pairs of antibodies equipped with DNA reporter molecules [15]. When binding in close proximity to their correct targets, the antibody pairs give rise to new DNA amplicons each ID-barcoding their respective antigens. The amplicons are subsequently quantified using the Fluidigm BioMark™ HD real-time PCR platform (South San Francisco, CA, USA). For data analysis and quality control Olink NPX Manager Software was used and the inter-plate variability was adjusted by intensity normalization. The final protein values are expressed as Normalized Protein eXpression (NPX) values which are on a log2 scale and one unit higher NPX represents a doubling of the measured protein concentration. Data quality was controlled and normalized using an internal and an interpolate control. Assay validation data for all proteins from the panel are available (www.olink.com).

2.7. Protein–Protein Interaction (PPI) Network Analysis

The Protein–Protein Interactions (PPI) Network analysis was performed using the search tool for retrieval of interacting genes (STRING) (https://string-db.org, accessed on 21 February 2021). The STRING database interaction evidence is thematically grouped into ‘channels’ (such as text mining, co-expression, and lab experiments) and limited to Homo sapiens. An interaction score > 0.4 was applied to construct the PPI networks. STRING performed identifier mapping (test the proteins of each known pathway for any nonrandom skew within the user-provided input values, and report statistically significant pathways) and displayed a network with all of the mapped proteins and their interconnections [16]. In the networks, the nodes correspond to the proteins and the edges represent the interactions. STRING was employed to seek potential interactions among proteins. The clustering coefficient, where 0 represents the absence of connections and 1 a fully connected network, was calculated quantifying the abundance of connected nodes in a PPI network. PPI enrichment p-value is used to indicate that the nodes are not random and that the observed number of edges is significant.

2.8. Statistical Analyses

All analyses were performed test using GraphPad Prism version 9.0. for Windows (GraphPad Software, San Diego, CA, USA). Patient characteristics were analyzed using one-way ANOVA and Kruskal–Wallis tests depending on the normality of the data to identify group wise differences. Inter-group differences in biomarker distributions were analyzed using the Mann–Whitney U test. The relationship between each of the 92 protein biomarkers and periodontal parameters was assessed using Spearman correlation analyses. To control the false discovery rate (FDR), the Benjamini–Hochberg procedure was applied to adjust p-values from multiple testing [17]. The significance level was defined at p ≤ 0.05.
For exploration of biomarker patterns within disease groups, principal component analysis (PCA) was performed. PCA is a powerful exploratory model statistically used for data exploration and simplification. The technique is based on generating principal components (latent variable) from the original dataset. The relationship of the principal components to the samples is referred to as ‘scores’, and that to the variables is called ‘loadings’. A threshold of 0.5 was deemed significant for variable loadings.

3. Results

3.1. Characteristics of Study Groups

The characteristics for the disease groups (RA with PD, RA without PD, and PD only) and controls are shown in Table 1. There were no differences in age amongst the four groups. The number of females were higher in disease groups as compared to controls. The clinical status comprised self-reported hypertension and diabetes confirmed by medication and prescription. The frequency of both conditions were similar amongst the disease groups. Periodontal parameters, waist circumference and HbA1c differed significantly amongst the groups with the highest medians in PD patients. The median for HbA1c value in the PD group classifies them as pre-diabetes overall (Supplementary File Table S3).

3.2. Group-Wise Biomarker Distribution

The distribution of 92 CVD biomarkers was assessed among the four groups. Two samples from the PD group were excluded due to unacceptable technical variations. Biomarkers with significantly increased levels in RA with PD groups as compared to RA without PD are shown in Figure 1. Higher NPX values were noted for 47% (43/92) of the markers which were: ANGPT1, BOC, CCL17, CCL3, CD4, CD84, CTRC, FGF-21, FGF-23, GLO1, HAOX1, HB-EGF, hOSCAR, HSP 27, IL16, IL-17D, IL18, IL-27, IL6, LEP, LPL, MERTK, MMP12, MMP7, NEMO, PAPPA, PAR-1, PARP-1, PD-L2, PGF, PIgR, PRELP, RAGE, SCF, SLAMF7, SRC, THBS2, THPO, TNFRSF13B, TRAIL-R2, VEGFD, VSIG2 and XCL1. For 32 of these biomarkers (BOC, CCL17, CCL3, CD84, CTRC, FGF-21, GLO1, HAOX1, HB-EGF, hOSCAR, HSP 27, IL-16, IL-17D, IL-18, IL-27, IL-6, LEP, LPL, MERTK, MMP12, NEMO, PAR-1, PARP-1, PD-L2, PRELP, RAGE, SCF, SLAMF7, THBS2, TNFRSF13B, TRAIL-R2, and XCL1) PD and RA with PD groups exhibited no differences (Supplementary File Table S2).

3.3. Correlation of CVD Biomarkers with Periodontal Parameters

The correlation between periodontal parameters and CVD-related biomarkers are shown in Table 2. The highest frequency of significant correlations was seen in the PD group for all parameters except for adjusted MBL. Anti-inflammatory marker IL-4RA was inversely related with three out of five indices for inflammation and pocketing. The Proto-oncogene tyrosine-protein kinase Src (SRC) was inversely correlated with four out of five indices for inflammation and pocketing. All correlations were direct amongst the RA with PD group, except for ACE-2. The most frequent and moderately strong correlations were noted with adjusted MBL. Least frequent correlations were noted in the RA without PD group. Dickkopf-related protein 1 (Dkk-1) and thrombospondin 2 (THBS2) were directly associated with Adjusted PPD Total. There was no overlap between the associated biomarkers amongst the three groups.

3.4. PCA

PCA was performed using standardized data and PC selection via parallel analysis. In the initial PCA output, selected component PC1 with loading structure >0.5 are shown for all disease groups (Table 3). The individual values show the correlation between the specific biomarker and the PC 1 for which the loading is calculated for. For RA with PD group, 64 biomarkers exceeded the 0.5 threshold of loading significance. Similarly, RA without PD had 65 biomarkers exceeding the threshold whereas PD group showed 55 biomarkers exceeding the threshold.
Visual representation of PC loadings plot (Figure 2) shows how the biomarkers are clustered closely together. In the disease groups, the majority of the biomarkers not only correlated strongly with each other but also with PC1 as most of the values were close to 1. The clustering pattern was more similar between the PD and RA without PD groups as some biomarkers showed stronger correlation with PC2. Using controls for reference, the loading plot showed weaker correlations between the biomarkers themselves and PC 1 and 2. The PC score plots reveal the variation in the dimensionality of the four groups.

3.5. Protein–Protein Interaction Network

The Protein–Protein interaction (PPI) network analysis of 43 proteins discriminating RA with PD from RA without PD is shown in Figure 3. The potential interactions between ANGPT1, BOC, CCL17, CCL3, CD4, CD84, FGF-21, FGF-23, HB-EGF, hOSCAR, HSP 27, IL16, IL-17D, IL18, IL-27, IL6, LEP, LPL, MMP12, MMP7, NEMO, PAPPA, PAR-1, PARP-1, PD-L2, PGF, RAGE, SCF, SRC, THBS2, THPO, TNFRSF13B, TRAIL-R2, VEGFD, and XCL1 yielded a clustering coefficient of 0.52, with a PPI enrichment p value < 0.0001. Markers from the CVD panel that also play a significant role in metabolic and skeletal disease areas were identified from the PC1 results for each disease group based on bioinformatic databases, including UniProt, Human Protein Atlas, Gene Ontology (GO), and DisGeNET. PPI network analysis was performed for three disease areas per group. These results are shown in Figure 4. The clustering coefficient was strongest for the PD group in all three disease areas when compared to the RA groups. The metabolic disease proteins were identical in clustering strength in both RA groups, uninfluenced by PD status.

4. Discussion

In this report, we identified 43 markers with a strong interactive network in patients suffering from PD, with and without RA. The risk of CVD exists in both PD and RA through shared pathological clusters. Several markers also increase associated metabolic and skeletal disease risk, independent of autoimmune status. In order to prevent CVD related morbidity and mortality in chronic inflammatory conditions, it is crucial to identify and study CVD risk biomarkers in the early stages of inflammatory disease. Studying a vast array of biomarkers that are significant in CVD development is an advantage offered by protein profiling using proteomic techniques. The biological mechanisms can be better understood with identification of early stage biomarkers which predispose RA and PD independently or combined to risk of CVD.
In our study, we examined an array of 92 biomarkers related to cardiovascular dysfunction and inflammation in RA patients with or without PD and PD patients alone. The disease groups showed a higher number of women of a relatively young age (<50 years). The gender dominance of women was expected since they are affected more by RA and seek dental care more frequently as compared to men [18]. In young women, being affected with RA is a risk for CVD [19]. RA female patients are 2.6-times more likely to develop CVD as compared to the general population. Our findings in relation to the age of the present cohort are, therefore, relevant.
An overall dysregulated level of HbA1c and increased waist circumference, a measure of central obesity, in PD patients has been confirmed previously as well [20]. Periodontal parameters were less severe in RA patients with PD and attributable to the use of disease modifying anti-rheumatic drugs (DMARDs) by the former group.
For direct comparisons, CVD biomarker distribution was assessed in all groups. Based on biological processes, the frequency of PD relevant biomarkers represented immune response (47%), cell adhesion (40%), intracellular mitogen-activated protein kinase (MAPK) signaling cascade (35%), inflammation (30%), catabolic process (23%), and proteolysis (19%). MAPKs are implicated as key regulators of inflammatory cytokines like IL-6 and TNF, thus transducing inflammation [21]. One of the contributing factors to CVD is endothelial dysfunction which is brought about by over expression of adhesion molecules due to inflammatory mediators [22].
The association of periodontal parameters with CVD biomarkers was also examined per disease group. In RA with PD, the associated biomarkers for periodontal pocketing spanned from enzymes (ACE-2) and membrane proteins (LOX-1) to plasma proteins (PTX3). The inverse relationship between ACE-2 levels and PPD Total scores reflect a pro-atherogenic state as ACE-2 levels have been detected in RA patients with a negative correlation with intima media thickness of carotid arteries [23]. Diseased probing sites correlated moderately with PTX3, also a pro-atherogenic inflammatory marker expressed by vascular endothelium known to modify angiogenesis and atherosclerotic lesion development [24]. Oxidized low density lipoproteins (ox-LDL) have been directly implicated in the pathogenesis of RA through signaling via the lectin-like ox-LDL receptor 1 (LOX-1) in the joint synovium [25]. LOX-1 activates downstream pathways that enhance atherosclerosis via endothelial dysfunction. LOX-1 is also expressed in platelets, where it enhances platelet activation and adhesion to endothelial cells [26]. Both LOX-1 and PTX3 associations with PPD Disease were moderately strong suggesting that PD contributes to a pro-atherosclerotic milieu in RA.
In PD patients only, three biomarkers (IL-4RA, SRC, MMP-12) conveyed a consistent pattern associated with deep pocketing. They reflect an unbalanced state with low anti-inflammatory IL-4RA levels confirming previous findings [27]. These findings corroborate a defect in the regulatory involvement of SRC and MMP-12 with phagocytosis and host defense mechanisms in PD patients. Low-MMP12 levels in periodontal tissues may be a risk factor underlying excessive pro-inflammatory IFN-γ macrophage activation in disease [28].
FGF-23, a bone-derived hormone, can also drive an increased production of pro-inflammatory cytokines [29]. Dkk-1 is known to play a pathophysiological role in bone erosion and joint remodeling in RA patients [30]. It negatively regulates the function of the Wnt pathway which is involved in the differentiation of osteoblasts. Thrombospondin 2 (THBS2) a matricellular protein, has been demonstrated as an endogenous regulator of angiogenesis and inflammation in the RA synovium [31].
High levels of LEP (leptin) associated with increased MBL in RA with PD patients further enforce previous findings of increased LEP levels found in dysfunctional immune phenotype including insulin resistance, inflammation, and disturbances in hemostatic factors [32]. TNSFR13B and its association with MBL reflects an increased B-cell proliferative and surviving capacity via its receptor BAFF (B-cell-activating factor). BAFF are up-regulated in RA synovial joints as well as early stages of PD [33,34]. IL-27′s enhancement of TNF-α mediated upregulation of adhesion molecules and pro-inflammatory IL-6 in blood monocytes of patients with acute myocardial infarction (MI) makes it high CVD risk associated [35]. The TGM2 levels in RA without PD correlate inversely with the total sum of MBL which aligns with previous findings that TGM2 correlates with RANKL production in human periodontal ligament cells as part of the inflammatory response in PD [36].
Protein biomarkers with high loadings (>0.5) on PC1 exceeded 50% of the total biomarkers analyzed in all disease groups. The biomarkers contributing to the greatest variance were similar in all three groups. Based on their disease–gene associations, these biomarkers are involved in vascular inflammation (HO-1, LPL, PAPPA, ADAMTS13, ADM, PGF, and GDF-2), hypertension, and arterial disease (PAPPA, ADAMTS13, ADM, and PGF). The underlying gene ontology represents upregulation of chemotaxis (XCL1 and CCL3), T helper 1 cytokines (SLAMF1, IL-18, XCL1, and IL-1ra), T helper 2 cytokines (XCL1), negative regulation of vasoconstriction (ADM and LEP), hematopoietic stem cell proliferation (THPO and ATXN) and increased bone loss (TNFSF11A and TF) [37].
The additional analyses of PPI networking for PC1 markers in CVD, metabolic and skeletal disease areas was performed as osseous and metabolic disturbances, especially insulin resistance, are highly frequent co-morbidities in both PD and RA [38,39,40,41]. The clustering coefficients displayed by PD group PC 1 biomarkers reflect a greater involvement of disease related proteins that make them a group with the highest risk for developing CVD, insulin resistance and skeletal diseases. The dampening of inflammatory circuits due to the use of NSAIDs and DMARDs in RA groups are to have some impact on the level of engagement amongst these proteins. Future studies are required to identify and validate markers of diagnostic and therapeutic relevance that may enhance the ‘treat-to-target’ strategy for RA and, hopefully, PD.
The limitations to our study pertain to the limited size of samples and the exploration of proteins which have been associated with cardiovascular diseases. Due to the exploratory nature of our study and the low prevalence of RA (~1%), we used a non-probability sampling method in which groups were not sex-matched. Despite these limitations, our findings have identified a direction for the exploration of other pathways in order to understand molecular alterations responsible for increased risk of CVD development in RA and PD.

5. Conclusions

We identified 43 markers with a strong interactive network in patients suffering from PD, with and without RA. In addition, several of these markers also increase associated metabolic and skeletal disease risk, independent of autoimmune status.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines10030714/s1, Table S1, Olink Cardiovascular II panel; Table S2; Table S3, Clinical Parameters.

Author Contributions

Conceptualization, J.P., M.Å., T.E.V.D., A.L. and P.-E.E.; methodology, J.P., A.G., M.A., M.Å., A.L. and P.-E.E.; software, J.P., M.Å. and A.L.; validation, J.P., A.G., M.A., M.Å., T.E.V.D., A.L. and P.-E.E.; formal analysis, J.P., M.Å. and A.L.; investigation, J.P., A.G., M.A., M.Å. and A.L.; resources, A.G., M.A. and P.-E.E.; data curation, J.P.; writing—original draft preparation, J.P.; writing—review and editing, J.P., A.G., M.A., M.Å., T.E.V.D., A.L. and P.-E.E.; visualization, J.P., T.E.V.D., A.L. and P.-E.E.; supervision, J.P., M.A., T.E.V.D., A.L. and P.-E.E.; project administration, J.P., A.G., M.A., T.E.V.D., A.L. and P.-E.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study protocol was approved by the regional research ethical committees in Karachi, Pakistan (26 September 2012, 30 September 2016) and Stockholm, Sweden (2016/296-31/1).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All relevant data are provided as Supporting Information in Tables S1–S3.

Acknowledgments

The authors would like to acknowledge support of SciLifeLab Affinity Proteomics Uppsala, (Uppsala University, Sweden) for aiding in protein analyses. All authors gave their final approval and agree to be accountable for all aspects of the work.

Conflicts of Interest

The authors declare no conflict of interest. Van Dyke is supported by USPHS, NIH, NIDCR DE025020.

References

  1. Hajishengallis, G.; Moutsopoulos, N.M.; Hajishengallis, E.; Chavakis, T. Immune and regulatory functions of neutrophils in inflammatory bone loss. Semin. Immunol. 2016, 28, 146–158. [Google Scholar] [CrossRef] [Green Version]
  2. Li, Z.; Lin, L.; Wu, H.; Yan, L.; Wang, H.; Yang, H.; Li, H. Global, Regional, and National Death, and Disability-Adjusted Life-Years (DALYs) for Cardiovascular Disease in 2017 and Trends and Risk Analysis from 1990 to 2017 Using the Global Burden of Disease Study and Implications for Prevention. Public Health Front. 2021, 9, 559751. [Google Scholar] [CrossRef]
  3. Gallino, A.; Aboyans, V.; Diehm, C.; Cosentino, F.; Stricker, H.; Falk, E.; Schouten, O.; Lekakis, J.; Amann-Vesti, B.; Siclari, F.; et al. Non-coronary atherosclerosis. Eur. Heart J. 2014, 35, 1112–1119. [Google Scholar] [CrossRef] [Green Version]
  4. Kobiyama, K.; Ley, K. Atherosclerosis. Circ. Res. 2018, 123, 1118–1120. [Google Scholar] [CrossRef]
  5. Beukers, N.G.; van der Heijden, G.J.; van Wijk, A.J.; Loos, B.G. Periodontitis is an independent risk indicator for atherosclerotic cardiovascular diseases among 60174 participants in a large dental school in the Netherlands. J. Epidemiol. Community Health 2017, 71, 37–42. [Google Scholar] [CrossRef] [Green Version]
  6. Makkar, H.; Reynolds, M.A.; Wadhawan, A.; Dagdag, A.; Merchant, A.T.; Postolache, T.T. Periodontal, metabolic, and cardiovascular disease: Exploring the role of inflammation and mental health. Pteridines 2018, 29, 124–163. [Google Scholar] [CrossRef]
  7. Aviña-Zubieta, J.A.; Choi, H.K.; Sadatsafavi, M.; Etminan, M.; Esdaile, J.M.; Lacaille, D. Risk of cardiovascular mortality in patients with rheumatoid arthritis: A meta-analysis of observational studies. Arthritis Care Res. 2008, 59, 1690–1697. [Google Scholar] [CrossRef]
  8. Santos, M.J.; Cordeiro, A.C.; Gil, V.M. Biomarkers of Cardiometabolic Risk, Inflammation and Disease; Springer: Cham, Switzerland, 2015. [Google Scholar]
  9. Sergeant, ESG, 2018. Epitools Epidemiological Calculators. Ausvet. Available online: http://epitools.ausvet.com (accessed on 20 May 2021).
  10. Potikuri, D.; Dannana, K.C.; Kanchinadam, S.; Agrawal, S.; Kancharla, A.; Rajasekhar, L.; Pothuraju, S.; Gumdal, N. Periodontal disease is significantly higher in non-smoking treatment-naive rheumatoid arthritis patients: Results from a case-control study. Ann. Rheum. Dis. 2012, 71, 1541–1544. [Google Scholar] [CrossRef] [PubMed]
  11. Aletaha, D.; Neogi, T.; Silman, A.J.; Funovits, J.; Felson, D.T.; Bingham, C.O., III; Birnbaum, N.S.; Burmester, G.R.; Bykerk, V.P.; Cohen, M.D.; et al. 2010 rheumatoid arthritis classification criteria: An American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum 2010, 62, 2569–2581. [Google Scholar] [CrossRef] [PubMed]
  12. Panezai, J.; Ali, A.; Ghaffar, A.; Benchimol, D.; Altamash, M.; Klinge, B.; Engström, P.E.; Larsson, A. Upregulation of circulating inflammatory biomarkers under the influence of periodontal disease in rheumatoid arthritis patients. Cytokine 2020, 131, 155117. [Google Scholar] [CrossRef] [PubMed]
  13. Panezai, J.; Ghaffar, A.; Altamash, M.; Sundqvist, K.G.; Engström, P.E.; Larsson, A. Correlation of serum cytokines, chemokines, growth factors and enzymes with periodontal disease parameters. PLoS ONE 2017, 12, e0188945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Little, R.R.; Rohlfing, C.L.; Sacks, D.B. National Glycohemoglobin Standardization Program (NGSP) Steering Committee Status of hemoglobin A1c measurement and goals for improvement: From chaos to order for improving diabetes care. Clin. Chem. 2011, 57, 205–214. [Google Scholar] [CrossRef] [PubMed]
  15. Assarsson, E.; Lundberg, M.; Holmquist, G.; Bjorkesten, J.; Thorsen, S.B.; Ekman, D.; Eriksson, A.; Rennel Dickens, E.; Ohlsson, S.; Edfeldt, G.; et al. Homogenous 96-plex PEA immunoassay exhibiting high sensitivity, specificity, and excellent scalability. PLoS ONE 2014, 9, e95192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021, customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef]
  17. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R Stat. Soc. Ser. B Stat. Methodol. 1995, 1, 289–300. [Google Scholar] [CrossRef]
  18. Lipsky, M.S.; Su, S.; Crespo, C.J.; Hung, M. Men and Oral Health: A Review of Sex and Gender Differences. Am. J. Men’s Health 2021, 15, 15579883211016361. [Google Scholar] [CrossRef]
  19. Fransen, J.; Kazemi-Bajestani, S.M.; Bredie, S.J.; Popa, C.D. Rheumatoid Arthritis Disadvantages Younger Patients for Cardiovascular Diseases: A Meta-Analysis. PLoS ONE 2016, 11, e0157360. [Google Scholar] [CrossRef]
  20. Lindström, J.; Tuomilehto, J. The diabetes risk score: A practical tool to predict type 2 diabetes risk. Diabetes Care 2003, 26, 725–731. [Google Scholar] [CrossRef] [Green Version]
  21. Thalhamer, T.; McGrath, M.A.; Harnett, M.M. MAPKs and their relevance to arthritis and inflammation. Rheumatology 2008, 47, 409–414. [Google Scholar] [CrossRef] [Green Version]
  22. Theofilis, P.; Sagris, M.; Oikonomou, E.; Antonopoulos, A.S.; Siasos, G.; Tsioufis, C.; Tousoulis, D. Inflammatory Mechanisms Contributing to Endothelial Dysfunction. Biomedicines 2021, 9, 781. [Google Scholar] [CrossRef]
  23. Braz, N.F.; Pinto, M.R.; Vieira, É.L.; Souza, A.J.; Teixeira, A.L.; Simoes-e-Silva, A.C.; Kakehasi, A.M. Renin–angiotensin system molecules are associated with subclinical atherosclerosis and disease activity in rheumatoid arthritis. Mod. Rheumatol. 2021, 31, 119–126. [Google Scholar] [CrossRef] [PubMed]
  24. CieŚlik, P.; Hrycek, A. Long pentraxin 3 (PTX3) in the light of its structure, mechanism of action and clinical implications. Autoimmunity 2012, 45, 119–128. [Google Scholar] [CrossRef] [PubMed]
  25. Ishikawa, M.; Ito, H.; Akiyoshi, M.; Kume, N.; Yoshitomi, H.; Mitsuoka, H.; Tanida, S.; Murata, K.; Shibuya, H.; Kasahara, T.; et al. Lectin-like oxidized low-density lipoprotein receptor 1 signal is a potent biomarker and therapeutic target for human rheumatoid arthritis. Arthritis Rheum. 2012, 64, 1024–1034. [Google Scholar] [CrossRef] [PubMed]
  26. Barreto, J.; Karathanasis, S.K.; Remaley, A.; Sposito, A.C. Role of LOX-1 (lectin-like oxidized low-density lipoprotein receptor 1) as a cardiovascular risk predictor: Mechanistic insight and potential clinical use. Arter. Thromb. Vasc. Biol. 2021, 41, 153–166. [Google Scholar] [CrossRef]
  27. Pradeep, A.R.; Roopa, Y.; Swati, P.P. Interleukin-4, a T-helper 2 cell cytokine, is associated with the remission of periodontal disease. J. Periodontal Res. 2008, 43, 712–716. [Google Scholar] [CrossRef]
  28. Dufour, A.; Bellac, C.L.; Eckhard, U.; Solis, N.; Klein, T.; Kappelhoff, R.; Fortelny, N.; Jobin, P.; Rozmus, J.; Mark, J.; et al. C-terminal truncation of IFN-γ inhibits proinflammatory macrophage responses and is deficient in autoimmune disease. Nature Commun. 2018, 9, 2416. [Google Scholar] [CrossRef]
  29. Czaya, B.; Faul, C. The role of fibroblast growth factor 23 in inflammation and anemia. Int. J. Mol. Sci. 2019, 20, 4195. [Google Scholar] [CrossRef] [Green Version]
  30. Rossini, M.; Viapiana, O.; Adami, S.; Fracassi, E.; Idolazzi, L.; Dartizio, C.; Povino, M.R.; Orsolini, G.; Gatti, D. In patients with rheumatoid arthritis Dickkopf-1 serum levels are correlated with parathyroid hormone, bone erosions and bone mineral density. Clin. Exp. Rheumatol. 2015, 33, 77–83. [Google Scholar]
  31. Park, Y.W.; Kang, Y.M.; Butterfield, J.; Detmar, M.; Goronzy, J.J.; Weyand, C.M. Thrombospondin 2 functions as an endogenous regulator of angiogenesis and inflammation in rheumatoid arthritis. Am. J. Pathol. 2004, 165, 2087–2098. [Google Scholar] [CrossRef] [Green Version]
  32. Wannamethee, S.G.; Tchernova, J.; Whincup, P.; Lowe, G.D.; Kelly, A.; Rumley, A.; Wallace, A.M.; Sattar, N. Plasma leptin: Associations with metabolic, inflammatory and haemostatic risk factors for cardiovascular disease. Atherosclerosis 2007, 191, 418–426. [Google Scholar] [CrossRef]
  33. Abe, T.; AlSarhan, M.; Benakanakere, M.R.; Maekawa, T.; Kinane, D.F.; Cancro, M.P.; Korostoff, J.M.; Hajishengallis, G. The B cell-stimulatory cytokines BLyS and APRIL are elevated in human periodontitis and are required for B cell-dependent bone loss in experimental murine periodontitis. J. Immunol. 2015, 195, 1427–1435. [Google Scholar] [CrossRef] [PubMed]
  34. Wei, F.; Chang, Y.; Wei, W. The role of BAFF in the progression of rheumatoid arthritis. Cytokine. 2015, 76, 537–544. [Google Scholar] [CrossRef] [PubMed]
  35. Al Shahi, H.; Shimada, K.; Miyauchi, K.; Yoshihara, T.; Sai, E.; Shiozawa, T.; Naito, R.; Aikawa, T.; Ouchi, S.; Kadoguchi, T.; et al. Elevated circulating levels of inflammatory markers in patients with acute coronary syndrome. Int. J. Vasc. Med. 2015, 2015, 805375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Matarese, G.; Currò, M.; Isola, G.; Caccamo, D.; Vecchio, M.; Giunta, M.L.; Ramaglia, L.; Cordasco, G.; Williams, R.C.; Ientile, R. Transglutaminase 2 up-regulation is associated with RANKL/OPG pathway in cultured HPDL cells and THP-1-differentiated macrophages. Amino Acids 2015, 47, 2447–2455. [Google Scholar] [CrossRef] [PubMed]
  37. Carbon, S.; Ireland, A.; Mungall, C.J.; Shu, S.; Marshall, B.; Lewis, S. AmiGO Hub, Web Presence Working Group. AmiGO: Online access to ontology and annotation data. Bioinformatics 2009, 25, 288–289. [Google Scholar] [CrossRef]
  38. Song, I.S.; Han, K.; Park, Y.M.; Ji, S.; Jun, S.H.; Ryu, J.J.; Park, J.B. Severe periodontitis is associated with insulin resistance in non-abdominal obese adults. J. Clin. Endocrinol. Metab. 2016, 101, 4251–4259. [Google Scholar] [CrossRef] [PubMed]
  39. Llorente, I.; García-Castañeda, N.; Valero, C.; González-Álvaro, I.; Castañeda, S. Osteoporosis in rheumatoid arthritis: Dangerous liaisons. Front. Med. 2020, 7, 802. [Google Scholar] [CrossRef]
  40. Avouac, J.; Elhai, M.; Forien, M.; Sellam, J.; Eymard, F.; Molto, A.; Banal, F.; Damiano, J.; Dieudé, P.; Larger, E.; et al. Influence of inflammatory and non-inflammatory rheumatic disorders on the clinical and biological profile of type-2 diabetes. Rheumatology 2021, 60, 3598–3606. [Google Scholar] [CrossRef]
  41. Hong, S.J.; Yang, B.E.; Yoo, D.M.; Kim, S.J.; Choi, H.G.; Byun, S.H. Analysis of the relationship between periodontitis and osteoporosis/fractures: A cross-sectional study. BMC Oral Health 2021, 21, 125. [Google Scholar] [CrossRef]
Figure 1. Group wise analyses for CVD-related biomarkers. Graphs 1–43 showing higher detection levels in RA with PD as compared to RA without PD. (1) ANGPT1, (2) BOC, (3) CCL17, (4) CCL3, (5) CD4, (6) CD84, (7) CTRC, (8) FGF-21, (9) FGF-23, (10) GLO1, (11) HAOX1, (12) HB-EGF, (13) hOSCAR (14) HSP 27, (15) IL-16, (16) IL-17D, (17) IL-18 (18) IL-27, (19) IL-6, (20) LEP, (21) LPL, (22) MERTK, (23) MMP7, (24), MMP12, (25) NEMO, (26) PAPPA, (27) PAR-1, (28) PARP-1, (29) PD-L2, (30) PGF, (31) PIgR, (32) PRELP, (33) RAGE, (34) SCF, (35) SLAMF7, (36) SRC, (37) THBS2, (38) THPO, (39) TNFRSF13B, (40) TRAIL-R2 (41) VEGFD, (42) VSIG2, and (43) XCL1. Data are presented as median with interquartile range. Group differences were calculated using Mann–Whitney U test. * p value ≤ 0.05, ** p value < 0.01, *** p value < 0.001, **** p value < 0.0001.
Figure 1. Group wise analyses for CVD-related biomarkers. Graphs 1–43 showing higher detection levels in RA with PD as compared to RA without PD. (1) ANGPT1, (2) BOC, (3) CCL17, (4) CCL3, (5) CD4, (6) CD84, (7) CTRC, (8) FGF-21, (9) FGF-23, (10) GLO1, (11) HAOX1, (12) HB-EGF, (13) hOSCAR (14) HSP 27, (15) IL-16, (16) IL-17D, (17) IL-18 (18) IL-27, (19) IL-6, (20) LEP, (21) LPL, (22) MERTK, (23) MMP7, (24), MMP12, (25) NEMO, (26) PAPPA, (27) PAR-1, (28) PARP-1, (29) PD-L2, (30) PGF, (31) PIgR, (32) PRELP, (33) RAGE, (34) SCF, (35) SLAMF7, (36) SRC, (37) THBS2, (38) THPO, (39) TNFRSF13B, (40) TRAIL-R2 (41) VEGFD, (42) VSIG2, and (43) XCL1. Data are presented as median with interquartile range. Group differences were calculated using Mann–Whitney U test. * p value ≤ 0.05, ** p value < 0.01, *** p value < 0.001, **** p value < 0.0001.
Biomedicines 10 00714 g001
Figure 2. PCA analysis. The principal component analysis showing loadings (left side) and scores (right side) RA with PD (panel A), PD (panel B), RA without PD (panel C), and controls (panel D).
Figure 2. PCA analysis. The principal component analysis showing loadings (left side) and scores (right side) RA with PD (panel A), PD (panel B), RA without PD (panel C), and controls (panel D).
Biomedicines 10 00714 g002aBiomedicines 10 00714 g002b
Figure 3. Protein–Protein interactions (PPI) showing networking of 43 CVD related biomarkers identified to be increased in RA with PD patients. The cluster shows frequent and strong interactions (represented by the same color of the nodes).
Figure 3. Protein–Protein interactions (PPI) showing networking of 43 CVD related biomarkers identified to be increased in RA with PD patients. The cluster shows frequent and strong interactions (represented by the same color of the nodes).
Biomedicines 10 00714 g003
Figure 4. PC1 biomarkers and their protein network analysis according to disease area in RA with PD, PD, and RA without PD groups. The network nodes represent proteins with red colored nodes denoting first shell interactors and green color showing second shell of interactors. All cluster coefficients (CC) have a PPI enrichment value of <0.0001.
Figure 4. PC1 biomarkers and their protein network analysis according to disease area in RA with PD, PD, and RA without PD groups. The network nodes represent proteins with red colored nodes denoting first shell interactors and green color showing second shell of interactors. All cluster coefficients (CC) have a PPI enrichment value of <0.0001.
Biomedicines 10 00714 g004
Table 1. Characteristics for RA and PD (disease) groups and controls.
Table 1. Characteristics for RA and PD (disease) groups and controls.
CharacteristicsDisease GroupsControlp-Value
RA with PD
(n = 26)
RA without PD
(n = 21)
PD
(n = 51)
(n = 20)
Age (years) a48.5 (8.8)43.1 (13.3)47 (9.5)43 (6.3)0.11
Female sex, n (%) b21 (81)20 (95)33 (65)8 (40)<0.001
Clinical Status, n b
  – Hypertension978-0.14
– Diabetes2610-
BOP c23 (57)43 (60)77 (56)15 (32)<0.0001
PPD Total c301(81)276 (86)384 (113)191 (24)<0.0001
PPD Disease c107.5 (104)0 (2.5)229 (136)0 (5)<0.0001
Adjusted PPD Total c11.6 (2.9)10.8 (2.8)15.5 (4.2)6.8 (1)<0.0001
Adjusted PPD Disease c8 (4.3)0 (0)10.4 (4)0 (0)<0.0001
∑MBL c27.4 (10.8) d13.5 (12.9) e34.2 (15.4)8.8 (17.5) f<0.0001
Adjusted ∑MBL c4.57 (1)3.02 (0.9)5.24 (2)2.88 (0.8)<0.0001
Body mass index (kg/m2) c24.2 (5)24.1 (6.2)25.2 (4)23.9 (4.6)0.35
Waist circumference (cm) c102 (30)97 (23)109 (19)86 (17)<0.0001
HbA1c % c5.0 (1)5.0 (2)5.7 (1.2)4.5 (0.8)<0.0001
BOP = bleeding on probing, PPD = probing pocket depth, MBL = marginal bone loss, HbA1c = glycated hemoglobin. a Differences in means were tested using one-way ANOVA test (testing overall difference among the three groups). b Differences in frequency were tested using χ² (chi-squared) test (testing overall difference among the three groups). c Differences in medians were tested using Kruskal–Wallis test (testing overall difference among the three groups). d Missing data (n = 5) was excluded in the analyses. e Missing data (n = 1) was excluded in the analyses. f Missing data (n = 1) was excluded in the analyses.
Table 2. Correlations of CVD risk biomarkers with periodontal pocketing and marginal bone loss.
Table 2. Correlations of CVD risk biomarkers with periodontal pocketing and marginal bone loss.
Periodontal Pocketing and InflammationMarginal Bone Loss
BOPPPD TotalPPD DiseaseAdj. PPD TotalAdj. PPD Disease∑MBLAdj. MBL
AnalyterAnalyterAnalyterAnalyterAnalyterAnalyterAnalyter
RA with PD ACE-2−0.42LOX-10.41 PTX30.44ANGPT1
PGF
0.47
0.45
LEP
TNFRSF13B
IL-27
0.48
0.48
0.46
PDCXCL1
SRC
XCL1
−0.31
−0.32
−0.36
IL-4RA−0.29IL-4RA
MERTK
SRC
−0.37
−0.28
−0.29
IL-4RA
MMP-12
SRC
−0.33
−0.31
−0.38
MMP-12
SRC
−0.34
−0.42
ADAM-TS13−0.29
RA without PDFGF-23−0.52 Dkk-1
THBS2
0.47
0.49
CD40-L
TGM2
−0.47
−0.45
Spearman rank correlation was used to identify correlations. All coefficients show biomarkers with adjusted p-values ≤ 0.05 after using the Benjamini–Hochberg procedure for multiple testing.
Table 3. PC loadings for disease groups.
Table 3. PC loadings for disease groups.
RA with PDPDRA without PD
VariablePC1PC2VariablePC1PC2VariablePC1PC2
C22D40.98−0.08PDGF subunit B0.87−0.12PDGF subunit B0.920.02
SCF0.960.05SOD20.87−0.05CD840.90−0.14
IL-17D0.950.00MMP70.860.12SCF0.90−0.31
PAR-10.93−0.16CD40.85−0.28BOC0.90−0.32
BOC0.93−0.16hOSCAR0.840.08CXCL10.900.11
PIgR0.93−0.02CCL170.84−0.08PD-L20.890.21
VEGFD0.930.20IL160.83−0.24MERTK0.870.22
IL160.930.07HB-EGF0.830.05VEGFD0.870.06
MMP70.920.15CCL30.830.11VSIG20.870.20
SPON20.91−0.14PIgR0.80−0.30MMP70.860.09
PDGF subunit B0.910.19VEGFD0.800.01THBS20.86−0.11
THPO0.91−0.18RAGE0.790.07BNP0.85−0.06
CD840.900.21SCF0.78−0.41PIgR0.84−0.33
hOSCAR0.900.25HSP 270.78−0.21PARP-10.82−0.36
LPL0.90−0.18IL-17D0.770.04HO-10.820.16
FABP20.89−0.23THBS20.77−0.02hOSCAR0.820.27
FGF-210.89−0.30CD840.76−0.14CD40.82−0.26
THBS20.88−0.06ADAM-TS130.750.17CCL170.820.28
CCL170.880.28PD-L20.750.05IL-17D0.82−0.20
CXCL10.880.30FGF-210.74−0.06DECR10.81−0.17
PARP-10.88−0.10HO-10.73−0.22FABP20.81−0.41
CTRC0.870.26Dkk-10.72−0.22GDF-20.810.15
ANGPT10.86−0.04BOC0.72−0.36RAGE0.800.17
PRELP0.85−0.44PAPPA0.710.03SORT10.800.20
MERTK0.85−0.06VSIG20.700.44PGF0.800.28
TM0.85−0.28SORT10.700.32ADAM-TS130.800.23
CCL30.840.18MERTK0.690.07CCL30.80−0.23
BMP-60.84−0.36GDF-20.680.33FGF-210.80−0.31
SOD20.830.02TNFRSF13B0.670.11HB-EGF0.780.48
SORT10.830.36CXCL10.660.07THPO0.78−0.39
STK40.81−0.07GLO10.660.03SLAMF70.770.39
VSIG20.81−0.04NEMO0.660.02LEP0.770.24
SRC0.810.10CTRC0.66−0.02Dkk-10.77−0.22
RAGE0.810.24SPON20.65−0.11TF0.76−0.01
PGF0.81−0.31CTSL10.650.32AGRP0.76−0.14
HO-10.81−0.06IL-1ra0.630.16IL-1ra0.760.18
FGF-230.80−0.28AGRP0.61−0.19NEMO0.75−0.06
IL180.79−0.18FGF-230.60−0.29TNFRSF13B0.740.19
Dkk-10.790.01IL1RL20.600.08PAR-10.74−0.57
AGRP0.78−0.27THPO0.60−0.37PSGL-10.73−0.04
XCL10.76−0.50DCN0.59−0.05IL160.73−0.39
HB-EGF0.760.53IL-270.570.18DCN0.730.25
HSP 270.760.20TNFRSF11A0.560.61HSP 270.72−0.01
TNFRSF13B0.75−0.15PTX30.560.08PAPPA0.72−0.22
PD-L20.740.50BNP0.56−0.05SOD20.71−0.38
NEMO0.740.39IL60.550.37TNFRSF11A0.710.56
BNP0.74−0.08TGM20.540.13TM0.700.32
LEP0.730.08PGF0.540.63LPL0.70−0.62
PAPPA0.730.42LPL0.54−0.67FGF-230.69−0.40
IL-1ra0.72−0.17PRSS270.540.08GLO10.660.16
TNFRSF11A0.70−0.24FABP20.53−0.53SPON20.650.05
DCN0.700.50PAR-10.52−0.66SRC0.64−0.33
IL1RL20.690.04TIE20.520.36CTRC0.640.10
PSGL-10.68−0.05ANGPT10.51−0.44IL1RL20.64−0.06
TF0.67−0.05IDUA0.500.04IL-270.630.47
ADM0.66−0.30IL180.500.33PRSS270.620.38
ADAM-TS130.660.55SRC0.49−0.39BMP-60.62−0.49
TNFRSF10A0.59−0.53LEP0.480.07IL180.610.30
MARCO0.57−0.24TM0.470.64TIE20.600.52
CTSL10.570.04CEACAM80.470.09SERPINA120.570.25
GDF-20.570.29STK40.47−0.49ANGPT10.55−0.48
PRSS270.560.35PARP-10.47−0.40STK40.54−0.53
TRAIL-R20.51−0.61HAOX10.47−0.04ITGB1BP20.54−0.20
IL-270.50−0.30DECR10.430.33REN0.530.31
GLO10.500.28TRAIL-R20.420.72PRELP0.51−0.77
IgG Fc receptor II-b0.44−0.31TF0.420.56IL60.490.01
SERPINA120.390.17TNFRSF10A0.400.63GH0.470.15
HAOX10.36−0.32ITGB1BP20.38−0.10GIF0.46−0.03
IL60.33−0.18PRELP0.37−0.63CTSL10.460.25
DECR10.310.26LOX-10.360.08TGM20.440.39
GH0.300.20BMP-60.36−0.58TNFRSF10A0.420.54
FS0.29−0.21XCL10.34−0.56XCL10.38−0.79
SLAMF70.29−0.35GH0.34−0.28TRAIL-R20.380.56
KIM10.25−0.46Gal-90.330.24PTX30.36−0.33
Gal-90.22−0.33MARCO0.30−0.14LOX-10.36−0.18
ITGB1BP20.200.36ADM0.29−0.38KIM10.340.68
PTX30.180.12AMBP0.260.64CD40-L0.320.07
TGM20.150.27CA5A0.250.34IgG Fc receptor II-b0.32−0.02
GIF0.12−0.18KIM10.240.74GT0.310.42
MMP120.11−0.54PSGL-10.190.23CA5A0.290.26
CEACAM80.090.25SERPINA120.180.15MARCO0.29−0.01
ACE20.08−0.49SLAMF70.170.10HAOX10.27−0.29
AMBP0.060.26GIF0.170.36IDUA0.250.36
TIE2−0.030.58REN0.150.35AMBP0.230.63
CA5A−0.04−0.44ACE20.130.54PRSS80.140.81
LOX-1−0.040.37IL-4RA0.120.31MMP120.110.34
IDUA−0.10−0.03CD40-L0.080.02Gal-90.100.35
REN−0.10−0.21IgG Fc receptor II-b0.030.24CEACAM80.060.09
GT−0.11−0.66PRSS80.020.81ADM0.03−0.25
CD40-L−0.20−0.04MMP12−0.040.32FS0.010.23
PRSS8−0.44−0.17FS−0.060.34ACE20.000.16
IL-4RA−0.46−0.76GT−0.130.16IL-4RA−0.010.26
Proportion of variance 47.0%8.9% 32.9%11.6% 43.2%11.3%
Cumulative proportion of variance 47.0%55.9% 32.9%44.5% 43.2%54.5%
PC = principal component. All loadings > 0.5 are in bold. The variance represented by two principal components in proportion and cumulatively are shown as percentages.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Panezai, J.; Ghaffar, A.; Altamash, M.; Åberg, M.; Van Dyke, T.E.; Larsson, A.; Engström, P.-E. Periodontal Disease Augments Cardiovascular Disease Risk Biomarkers in Rheumatoid Arthritis. Biomedicines 2022, 10, 714. https://doi.org/10.3390/biomedicines10030714

AMA Style

Panezai J, Ghaffar A, Altamash M, Åberg M, Van Dyke TE, Larsson A, Engström P-E. Periodontal Disease Augments Cardiovascular Disease Risk Biomarkers in Rheumatoid Arthritis. Biomedicines. 2022; 10(3):714. https://doi.org/10.3390/biomedicines10030714

Chicago/Turabian Style

Panezai, Jeneen, Ambereen Ghaffar, Mohammad Altamash, Mikael Åberg, Thomas E. Van Dyke, Anders Larsson, and Per-Erik Engström. 2022. "Periodontal Disease Augments Cardiovascular Disease Risk Biomarkers in Rheumatoid Arthritis" Biomedicines 10, no. 3: 714. https://doi.org/10.3390/biomedicines10030714

APA Style

Panezai, J., Ghaffar, A., Altamash, M., Åberg, M., Van Dyke, T. E., Larsson, A., & Engström, P. -E. (2022). Periodontal Disease Augments Cardiovascular Disease Risk Biomarkers in Rheumatoid Arthritis. Biomedicines, 10(3), 714. https://doi.org/10.3390/biomedicines10030714

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop