Skip to main content
Download PDF

Low serum IGF1 is associated with hypertension and predicts early cardiovascular events in women with rheumatoid arthritis

BMC Medicine volume 17, Article number: 141 (2019) Cite this article

Abstract

Objectives

Since low insulin-like growth factor (IGF) 1 is often linked to inflammation, we analyze whether serum levels of IGF1 are associated with cardiovascular disease (CVD) in rheumatoid arthritis (RA) in a longitudinal observational study.

Methods

A CVD risk was estimated (eCVR) in 184 female RA patients (mean age 52 years) and in 132 female patients after ischemic stroke (mean age 56 years) with no rheumatic disease, using the Framingham algorithm. The median level of IGF1 divided the cohorts in IGF1high and IGF1low groups. A 5-year prospective follow-up for new CVD events was completed in all RA patients. The Mantel-Cox analysis and event-free survival curves were prepared. Unsupervised clustering of proteins within the IGF1 signaling pathway was employed to identify their association with eCVR.

Results

Low IGF1 resulted in a higher eCVR in RA patients (7.2% and 3.3%, p = 0.0063) and in stroke (9.3% and 7.1%, p = 0.033). RA had higher rate for new CVD events at prospective follow-up (OR 4.96, p = 0.028). Hypertension was the major risk factor associated with low IGF1 in RA and stroke. In hypertension, IGF1 was no longer responsible for intracellular activation and lost its correlation to IRS1/2 adaptor proteins. The clustering analysis confirmed that combination of low IGF1 and IRS1/2 with high IL6, insulin, and glucose predisposed to high eCVR and emphasized the functional role of serum IGF1.

Conclusions

Low serum IGF1 precedes and predicts development of early CVD events in female RA patients. Hypertension and aberrant IGF1 receptor signaling are highlighted as the important contributors to IGF1-related CVD events.

Peer Review reports

Background

Hypertension is a significant health threat and an independent predictor of CV events including coronary heart disease, stroke, heart failure, and dementia in the general population [1] and in RA [2]. The prevalence of hypertension in rheumatoid arthritis (RA) is high, but varies between 4 and 72% producing inconsistent results in the repeated meta-analysis [3,4,5]. Studies in RA population often indicated high rates of hypertension in patients with first stroke [6,7,8,9] and the highest stroke risk in young patients < 50 years [10]. Chronic inflammation, the hallmark of RA and an important contributor to the high CVD rate in these patients, is viewed as a major pathogenic process behind significant vascular dysfunction. Inflammation reduces artery elasticity and increases systemic vascular resistance, which consequently leads to an increased arterial blood pressure (BP) [11,12,13]. Elevated C-reactive protein levels are associated with hypertension in the general population by reducing the production of endothelial nitric oxide and by triggering platelet activation and leukocyte adhesion. Pro-inflammatory cytokines TNF-a and IL6 maintain hypertension by inducing endothelial proliferation and increasing vascular permeability and blood volume [14, 15].

Solid biological evidence connects insulin-like growth factor 1 (IGF1) with the regulation of endothelial cell function. Vascular endothelial and smooth muscle cells express IGF1 receptor (IGF1R), which mediates angio-protective effects of IGF1. Locally produced IGF1 supports proliferation and migration of endothelial progenitors essential for blood vessel reparation and controls oxidative stress triggered by inflammation [16,17,18]. Circulating IGF1 is known to induce vasodilation, which contributes to the regulation of arterial BP and vascular tone, while deletion of IGF1 in mice enhanced mechanisms of vasoconstriction leading to hypertension. Clinically, both a direct and inverse relation between IGF1 levels and BP has been reported. Early studies with focus on patients with congenital or acquired abnormalities in GH/IGF1 production, indicated an association between high IGF1 production and hypertension [19, 20]. Recent studies of general population with low prevalence of obvious endocrinological conditions reported an inverse association between IGF1 and BP [21,22,23,24]. Low serum IGF1 is tightly linked to vascular changes in aging [25] and obesity [26] and to the increased cardiovascular morbidity and mortality [27,28,29,30].

Disturbances in IGF1/IGF1R signaling are notable for RA pathology. It contributes to joint inflammation and affects homeostasis in chondrocytes, leukocytes and synovial fibroblasts [31,32,33]. It has been shown that inflammation is directly related to high expression of IGF1R in leukocytes [34], which might be responsible for higher pain perception in RA [33, 35]. In a study combining two independent RA cohorts, we demonstrated that smoking and inflammation had additive suppressive effects on serum IGF1 levels [36]. These reports imply IGF1/IGF1R signaling to be a natural protective response to inflammation and pain. Little is known whether IGF1 contributes to the prevalence of CVD in RA. Therefore, the objective of this study was to determine the importance of low IGF1 and alterations in IGF1/IGF1R signaling for CVD in RA and to assess whether similarities exist in patients without RA afflicted by ischemic stroke.

Methods

The RA cohort

This study involved 184 female RA patients consecutively enrolled at the Rheumatology Clinic of Sahlgrenska University Hospital, Gothenburg, and the Northern Älvsborg County Hospital, Uddevalla [33, 36]. All the patients fulfilled the classification criteria for RA (ACR1987) and had median disease duration of 7 years. The dominating majority (89%) were RF and/or ACPA positive. Active anti-rheumatic treatment had 175 patients. Among those, monotherapy with MTX had 82 patients (47%), other DMARDs 3 patients (1.7%); combination therapy with DMARDS 45 patients (26%), TNF-a-inhibitors 30 patients (17%), and other biologics 15 patients (8.6%). The study is registered at the Clinical Trials.gov with ID NCT03449589.

The ischemic stroke (IS) cohort

We use clinical and serological data collected at 3 months after index IS of 132 consecutive female patients recruited at the Sahlgrenska University Hospital in the frame of the Sahlgrenska Academy Study on Ischemic Stroke (SALHSIS) [30, 37].

The studies are approved by the Swedish Ethical Review Authority (RA cohort, Dnr.659-11 and SALHSIS cohort, Dnr.Ö469-99) and are conducted in accordance with the ethical principles of the Helsinki Declaration and Good Clinical Practice. All patients provided written informed consent before enrolment into the study.

Cardiovascular risk

A 10-year risk of developing CVD was estimated (eCVR) using a digital version of the Framingham algorithm [38] and included sex, age, systolic BP, treatment for hypertension, current smoking, diabetes mellitus, HDL, and TC.

Prospective follow-up

Five years after enrollment the RA patients were contacted for a structured telephone interview. The questions were asked for any CVD events, cases of type 2 diabetes (T2D) and current medication including anti-hypertensive drugs, anticoagulants, anti-diabetic drugs, and statins. The reported CVD events and changes in medications were then controlled against medical records and the Swedish National Patient Registry.

Hypertension was recorded if an incidental measure of systolic BP was > 140 mmHg, or diastolic BP was > 90 mmHg or pharmacological treatment for hypertension was used. The BP burden was calculated as a sum of systolic and diastolic BPs.

Disease activity score (DAS28) of RA was calculated based on assessment of 28 tender and swollen joints and ESR.

Blood sampling and storage

The samples were collected between 7 and 10 o’clock in the morning after overnight fasting. For serum preparation, the blood was obtained from the cubital vein into vacuum containers (BD Vacutainer) and for RNA preparation into PAXgene protection tubes (Becton Dickinson, Franklin Lakes, NJ, USA). Serum samples were stored at − 70 °C and PAXtubes in − 20 °C until use.

Serological measurements

In RA samples, serum IGF1, total cholesterol (TC), triglycerides, high-density lipoprotein (HDL), and low-density lipoprotein (LDL) were measured by photometry on Cobas 8000 (Roche Diagnostics, Switzerland). In the IS samples, serum IGF1 was measured with a radio-immune assay (Mediagnost, Reutlingen, Germany) [39]. Plasma glucose levels were measured using FreeStyle Lite (Abbott Diabetes Care Ltd., Oxon, UK). Sandwich ELISAs were used to measure insulin (DY8056, R&D Systems, Minneapolis, MN, USA) and IL6 and IL1b (M9316 and M1934, respectively; Sanquin, Amsterdam, the Netherlands).

Gene expression analysis

Total mRNA was prepared using PAXgene Blood RNA kit (Qiagen). Complementary DNA was synthesized using High Capacity cDNA Reverse Transcription kit (Applied Biosystems, Foster City, CA). Amplification of the gene product was attained on a ViiA™7 Real-Time PCR (Applied Biosystems) using SYBR Green qPCR Mastermix (SA Biosciences, Qiagen) and primer pairs as reported [33, 40]. Gene expression levels were calculated by the ddCt method and presented as relative quantity to the average expression in the IGF1hi group.

Statistical analysis

The SPSS v.25 (IBMSPSS, Chicago, IL), GraphPad prism v.7, www.open-epi.com, R v.3.3.0 (R Core Team, 2018) and R studio v.1.1.447 (http://www.rstudio.com/) were used for the analysis. Data are presented as mean ± SD, median [IQR], or in absolute numbers. Missing data for BP (5%) and ESR (8%) were imputed using the linear regression (SPSS). The study cohorts were dichotomized into IGF1hi and IGF1low groups by the median level. Continuous data were analyzed using the Mann-Whitney U test, the Kruskal-Wallis test followed by Dunn’s post hoc test, and Spearman’s correlation test. Relative risk prediction was done using the area under the receiver operative characteristic (ROC) curve. The Kaplan-Meier curves and the Mantel-Cox analysis were used to compare the groups. For the clustering analysis, the data were log normalized and ranked by row. Heatmaps and hierarchical clustering were performed in R, using the stats and gplots packages, Spearman correlation-based distances, and Ward2 linkage. All tests were two-tailed and conducted with 95% confidence.

Results

Low IGF1 levels are associated with higher CVD risk in RA patients

Consistent with a clinically relevant IGF1-deficiency, the IGF1low group was lower in height compared to IGF1hi group (Table 1). The IGF1low group had significantly higher eCVR compared to the IGF1hi group (7.0% vs. 3.2%, p = 0.0063), and the absolute IGF1 levels correlated negatively to the eCVR (Fig. 1a, b, d). Since both eCVR and IGF1 levels were age dependent, the analysis was performed separately within patients of different age. Expectedly, the eCVR had a gradual increase with age, both in the IGF1low and IGF1hi groups (Fig. 1b). The increase in eCVR was predominantly attributed to IGF1low patients < 50 years. Notably, these patients were associated with significantly higher BP burden, BMI, TC, and LDL, compared to the age-matched IGF1hi group (Fig. 1d). Notably, the IGF1low RA patients < 50 years reached a low threshold of IGF1 levels with no further decline (Fig. 1b). This low IGF1 was associated with significantly higher BP burden, BMI, TC, and LDL, compared to the age-matched IGF1hi group (Fig. 1e).

Table 1 Clinical characteristics of female RA and stroke cohorts
Full size table
Fig. 1
figure 1

Low serum levels of IGF1 are associated with higher estimated cardiovascular risk (eCVR) in RA patients. a eCVR was calculated in 184 female RA patients and in 132 female incidental ischemic stroke using the Framingham lipid algorithm. The median level of IGF1 formed the IGF1high and IGF1low groups. Box plots show eCVR separately for IGF1low and IGF1hi groups stratified by age. Embedded box plots show eCVR in the total cohorts. b Box plots show levels of IGF1 separately for IGF1low and IGF1hi groups stratified by age. Box plots present median, interquartile range. P values are calculated with the Mann-Whitney U test. c The forest plot shows the difference in traditional CVR factors between IGF1low and IGF1hi groups is shown as odds ratio (OR) with 95% confidence interval (CI). Filled circles indicate RA cohort (n = 184), open circles Ischemic stroke cohort (n = 132). The p values are obtained by chi-square statistics. BMI, body mass index; TC, total cholesterol; T2D, type 2 diabetes. d The Spearman correlation between eCVR and serum levels of IGF1 in the RA cohort is shown as a dot plot. The dotted line indicates the median IGF1 level separating IGF1low and IGF1hi groups. Red dots indicate the patients, who developed CV events during the prospective follow-up. e The comparison of CVR factors between IGF1low (n = 18) and IGF1hi (n = 50) patients within the youngest age group < 50 years. Box plots present median, interquartile range. P-values are calculated with the Mann-Whitney U test. BPB, blood pressure burden (a sum of systolic and diastolic BP); BMI, body mass index; DD, disease duration; TC, total cholesterol; LDL, low-density lipoproteins; HDL, high-density lipoproteins

Full size image

The comparison of the traditional CVR factors between IGF1low and IGF1hi groups showed that hypertension was the most prominent distinctive feature (Fig. 1c). IGF1low group had higher rate of pharmacological treatment for hypertension. The incidental measures of systolic and diastolic BP in IGF1low groups were also higher compared to IGF1hi (Table 1). Hypercholesterolemia and higher levels of BMI, LDL, and HDL were also characteristic for the IGF1low group, while there was no association with diabetes or smoking (Fig. 1c, Table 1).

Low IGF1 levels and CVD risk in IS patients

Age and serum IGF1 levels in the IS cohort were comparable with RA patients (Table 1, Fig. 1b). The IGF1low IS patients had higher diastolic BP and eCVR (9.3% vs. 7.0%, p = 0.033) (Fig. 1a). However, there was no significant difference between IGF1low and IGF1hi IS groups with respect to the individual CVR factors (Table 1, Fig. 1c). Compared to the RA cohort, the rate of hypertension, smoking, and diabetes was higher in the IS cohort, whereas levels of TC, HDL, and LDL were lower. Pharmacological treatment of hypertension and hyperlipidemia, which is prescribed actively during the early post-stroke phase, will probably explain this difference.

CVD risk in relation to RA disease and inflammation

Clinical characteristics of RA patients are shown in Fig. 2a. Analysis of the RA-related CVR including RA-specific antibodies, active RA with DAS28 ≥ 3.2, and disease duration > 10 years, identified no differences between the IGF1low and IGF1hi groups (Fig. 2b). Despite the absence of obvious differences in clinical disease activity and similar numbers of swollen and tender joints, IGF1low group had the higher frequency and absolute measures of ESR and IL6 (Fig. 2b, d). This led us to ask if IGF1low and IGF1hi patients needed a similar amount of treatment to control the RA disease. To investigate this, we compared the frequency of MTX monotherapy and combination therapy in the groups (Fig. 2c). The IGF1low group was significantly more often treated with MTX monotherapy (OR 2.26, p = 0.007), and the dose of MTX was similar to IGF1hi group (18.7 vs. 15.0 mg/week, ns), which implied insufficient inflammation control. The combination of MTX with other DMARDs or TNF-a inhibitors was not prevalent in the IGF1hi group (Fig. 2c).

Fig. 2
figure 2

Cardiovascular risk factors attributed to rheumatoid arthritis. a The table shows patient characteristics of the RA cohort. Forest plots show the difference between IGF1low (n = 96) and IGF1hi (n = 88) groups in RA-specific parameters (b) and in anti-rheumatic treatment (c). Lines indicate odds ratio (OR) with 95% confidence interval (CI). The p-values are obtained by chi-square statistics. d The comparison of clinical and serological parameters between IGF1low (filled boxes) and IGF1hi (open boxes) groups. Box plots show median, interquartile range. P-values are calculated with the Mann-Whitney U test. RA, rheumatoid arthritis; ESR, erythrocyte sedimentation rate; DD, disease duration, DAS28, disease activity score with assessment of 28 joints; RF, rheumatoid factor; ACPA, antibodies to cyclic citrullinated peptides; MTX, methotrexate; DMARDs, disease modifying anti-rheumatic drugs; TNFi, tumor necrosis factor inhibitors

Full size image

Low IGF1 associates with new CVD events

At the 5-year follow-up, new CVD events occurred in 12 of 184 RA women (6.5%). The baseline characteristics of patients with new CVD events are shown in Fig. 3a. We learned that new CVD events often belong to established complications of hypertension and include 4 cases of ischemic stroke, 3 atrial fibrillations, and 2 aorta aneurysms. Patients with new CVD events had longer disease (14 vs. 7 years, p = 0.006) than the average of the RA cohort. They were slightly older (57.5 vs. 53 years, not significant) and had the eCVR of 9.1%, similar to the IS cohort.

Fig. 3
figure 3

Predictive performance of serum IGF1 levels for development of cardiovascular events in rheumatoid arthritis. a The table shows individual baseline characteristics of 12 patients, who developed CV events during the 5-year follow-up. Asterisk indicates the case of death. b The receiver operator characteristic (ROC) curves depict associations between CV events developed during the 5-year follow-up and age (red line) and absolute serum levels of IGF1 (blue line) at baseline in 184 female RA patients. AUC, area under the ROC curve. c The Kaplan-Meier curves show the development of new CV events in IGF1hi (n = 88, red line) or IGF1low (n = 96, blue line) groups. d The column diagrams show the frequency of treatment at baseline (BL, open columns) and at the 5-year follow-up (5Y, filled columns) in IGF1low and IGF1hi groups stratified by age, < 50 years (n = 68), 50–60 years (n = 53), > 60 years (n = 63). The p values are obtained by chi-square statistics. eCVR, CV risk estimated by the Framingham algorithm; DD, disease duration; IS, ischemic stroke; TIA, transitory ischemic attack; DVT, deep venous thrombosis; LE, lung embolism; MI, myocardial infarction; MTX, methotrexate, RTX, rituximab; ETC, etanercept; SLZ, sulfasalazine; AZA, azathioprine; TOCI, tocilizumab; OR, odds ratio; CI, confidence interval; T2D, type 2 diabetes

Full size image

To assess the clinical importance of serum IGF1 for development of new CVD events, we performed the ROC analysis of new CVD events using the absolute levels of IGF1 and age as independent variables. The ROC curve identified a strong negative relation between IGF1 and CVD events (Fig. 3b) and confirmed higher eCVR in IGF1low group (Fig. 1a). The ROC curve with age showed no association to CVD events and confirmed independent of age input of IGF1 to the developed CVD events. Indeed, the absolute CVD risk in IGF1hi group was 2.3% and in IGF1low group, it was 10.4% giving the relative risk of CVD event 3.47 times higher for the patients in the IGF1low group. To assess the probability of future CVD events, we constructed the event-free survival curves for IGF1low and IGF1hi groups (Fig. 3c). The Mantel-Cox analysis of the curves revealed an almost 5 times higher probability of CVD event in IGF1low group (OR 4.96, p = 0.028).

During the 5-year period, we observed a remarkable increase in treatment for hypertension to occur in IGF1low but not in the IGF1hi group (Fig. 3d). This increase was significant among patients < 50 years and ≥ 60 years, which could have partly limited development of new CVD events in young RA patients.

Hypertension is associated with an imbalance in the signaling through IGF1 receptor

To assess how low serum IGF1 affects IGF1R signaling in leukocytes, we compared the transcription of individual proteins within this pathway between IGF1low and IGF1hi groups in RA. We observed that IGF1low group had a significantly higher mRNA of IGF1R and AKT1, and low mRNA of the adaptor proteins IRS1 and IRS2 (Fig. 4a). To study how these changes are related to hypertension, the clinical hallmark of low IGF1, we compared internal correlations between the individual components in normotensive and hypertensive patients (Fig. 4b). We observed that the negative correlations between IGF1 and IRS1, IRS2 and plasma glucose, which were present in the normotensive patients, could no longer be found in the hypertensive patients. Instead, novel correlations between IGF1R, IRS1, IRS2, and insulin appeared with hypertension demonstrating that IGF1 is no longer responsible for activation of intracellular processes through IGF1R. Finally, we performed an unsupervised clustering of RA patients by individual components of the IGF1 signaling pathway aiming to identify associations with the increased eCVR. Unsupervised clustering analysis of 60 patients with the complete set of the individual components formed four independent clusters (Fig. 4c). At visual inspection, cluster 4 had highest eCVR (mean 8.51) combined with lower IGF1, IRS1, and IRS2, while serum IL6, insulin, and glucose was high. This cluster differed significantly from the other clusters by eCVR (p = 0.0024), DAS28 (p = 0.047) and age (p = 0.0006), but had similarities in eCVR with patients developed new CVD events (Fig. 3a). Cluster 2, opposite to cluster 4, had the lowest eCVR, high levels of IGF1 and IRS1/2 and low IL6, insulin, and glucose. Clusters 1 and 3 had intermediate eCVR and low levels of IGF1, IRS1/2, IL6, insulin, and glucose (Fig. 4c).

Fig. 4
figure 4

The IGF1 signaling pathway is related to hypertension and cardiovascular risk. a Box plots show mRNA levels in leukocytes of patients within IGF1hi or IGF1low groups. b Correlation matrix shows Spearman’s R values in color, as indicated the color key code, in patients with normal blood pressure (BP) or hypertension. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. c Heatmap shows individual components related to IGF1 signaling of 60 RA patients. Unsupervised clustering analysis was done using the Spearman’s correlation-based distances and Ward2 linkage and formed 4 clusters. The table represents the mean and standard deviation of the parameters by cluster. The difference between clusters 1–3 and 4 was calculated with the t-test. RQ, relative quantity; MTX, methotrexate; TNFi, tumor necrosis factor inhibitors; DMARD, disease modifying anti-rheumatic drugs; DAS28, disease activity score with an assessment of 28 joints; DD, disease duration

Full size image

Discussion

In this cross-sectional study, we explored longitudinal associations between serum IGF1 and early CVD in RA. We showed harmful consequences of normal-low levels of IGF1 for the increase in the eCVR in female RA patients < 50 years and premature development of new CVD events. The profile of new CVD events consisted of the established complications of hypertension including ischemic stroke, atrial fibrillation, and aorta aneurysm, which indicated that hypertension was an essential clinical sign of IGF1-related CVD of the studied patients. Indeed, the patients with low IGF1 demonstrated higher BP measures at baseline and a significant increase in treatment for hypertension at the prospective follow-up. These findings in RA patients are consolidated by the reference cohort of IS females of similar age, where a significant association between low IGF1, hypertension, and eCVR is also observed.

The cumulative knowledge gained in a recent meta-analysis supports a direct causality between low IGF1 and hypertension [22]. Although the studied cohorts had no patients with extreme levels of IGF1, we observed a pronounced negative effect of normal-low IGF1 levels on the development of CVD and on its arterial localization. Our findings offer an appropriate molecular context to the high rate of stroke previously reported in young RA patients [10]. We found that the IGF1-related increase in eCVR was most noticeable in young women < 50 years. This implies that the vascular effects related to IGF1 may precede severe atherosclerotic change in the vessel wall and predispose to its early development. Indeed, previous studies in RA connected early CV events with endothelial dysfunction [11, 12]. Direct angioprotective properties of IGF1 for vascular endothelial cell function have been extensively studied and is experimentally verified [16]. In clinical studies, substitution with IGF1 was often followed by resolution of hypertension and improvement of the intima-media thickness of the arterial vessel wall [41,42,43]. Importantly, IGF1 has an immediate effect on repairing tissue damage after acute ischemic cerebral and myocardial injury [39, 44, 45]. In RA, the factors enhancing IGF1 levels attracted less attention. In this study, we observed that the use of TNF-a inhibitors was associated with higher IGF1 levels, while patients with low IGF1 received mostly MTX monotherapy. Treatment with TNF-a inhibitors is suggested to improve hypertension through an endothelium-dependent mechanism [46, 47], which could indicate a contribution of increasing IGF1 levels. We would expect that anti-rheumatic treatment assisting production and bioavailability of IGF1 to be an important tool to recover cardiovascular health in RA patients. Physical inactivity, a condition accumulated in long-lasting RA and in ischemic stroke [8], is frequently associated with lower serum IGF1 [48] that could be increased with physical exercise [49]. Changing lifestyle to improve serum IGF1 levels could be an attractive alternative for RA patients [50].

In the present study, we investigated how the difference in serum IGF1 changed the IGF1R signaling and its relation to eCVR. From analyzing leukocytes, we confirmed that serum IGF1 controlled the relative expression of proteins within the IGF1R pathway. Importantly, the relation between these proteins changes with hypertension, demonstrating its intimate connection with low IGF1. The unsupervised clustering of the IGF1-related proteins proposed the high CVR signature, which combined the low expression of IGF1, IRS1, and IRS2 with high serum IL6, insulin and plasma glucose and emphasized the functional role of serum IGF1 in the development of early CVD in RA. These findings are in line with our recent report in experimental arthritis, which puts forward the role of excessive inhibitory phosphorylation in IRS1 to create a condition of insulin/IGF1 resistance [35]. The situation seems more complex in patients, where low IGF1 is concurrent with high IGF1R and low adaptors IRS1/2. Being an ancestor of pro-insulin, IGF1 modulates carbohydrate metabolism, which stimulates glucose transport and inhibits insulin sensitivity [18]. In the general population and experimentally, IGF1 deficiency leads to glucose intolerance and T2D, the conditions rare among the studied RA patients. This was despite the fact that metabolic deviations as obesity, hyperlipidemia and hyperinsulinemia were expectedly frequent in IGF1low group.

There are notable strengths and limitations in our study. Among strengths, this is to our knowledge the largest study of the relation between IGF1 levels and cardiovascular morbidity in RA, and the only one to explore the longitudinal association with the development of CVD events. An additional strength is the similarity of the RA and IS female cohorts, which comprise women of the same age. This study is done on the female cohort of RA patients and allows no extrapolation of the results on male RA patients. Additionally, the study is restricted to the risk of CV events and does not cover cerebrovascular health in female RA patients, since the vascular wall has not been evaluated. Hence, it is outlined among the limitations of our study. We believe that the chosen methodology does not affect the results obtained in the between-group comparison.

Conclusion

Taken together, the results of this study attract attention to serum levels of IGF1 as a valuable parameter for estimating CVD risk in female RA patients. Low serum IGF1 precedes and predicts the development of early CVD events due to a tight connection with hypertension, which is important for future management of RA patients.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ACPA:

Anti-citrullinated protein antibodies

AZA:

Azathioprine

BMI:

Body mass index

BP:

Blood pressure

BPB:

Blood pressure burden

CI:

Confidence interval

CVD:

Cardiovascular disease

CVR:

Cardiovascular risk

DAS28:

Disease activity score based on 28 joints and ESR

DD:

Disease duration

DMARD:

Disease modifying anti-rheumatic drugs

DVT:

Deep venous thrombosis

eCVR:

Estimated cardiovascular risk

ESR:

Erythrocyte sedimentation rate

ETC:

Etanercept

GH:

Growth hormone

HDL:

High-density lipoprotein

IGF1:

Insulin-like growth factor 1

IGF1R:

IGF1 receptor

IL:

Interleukin

IQR:

Inter-quartile range

IRS:

Insulin receptor substrate

IS:

Ischemic stroke

LDL:

Low-density lipoprotein

LE:

Lung embolism

MI:

Myocardial infarction

MTX:

Methotrexate

OR:

Odds ratio

RA:

Rheumatoid arthritis

RF:

Rheumatoid factor

ROC:

Receiver operator characteristics

RTX:

Rituximab

SD:

Standard deviation

SLZ:

Sulfasalazine

T2D:

Type 2 diabetes mellitus

TC:

Total cholesterol

TIA:

Transitory ischemic attack

TNF:

Tumor necrosis factor

TOCI:

Tocilizumab

References

  1. Vasan RS. High blood pressure in young adulthood and risk of premature cardiovascular disease: calibrating treatment benefits to potential harm. JAMA. 2018;320(17):1760–3. https://doi.org/10.1001/jama.2018.16068 [published Online First: 2018/11/07].

    Article  PubMed  Google Scholar 

  2. Crowson CS, Rollefstad S, Ikdahl E, et al. Impact of risk factors associated with cardiovascular outcomes in patients with rheumatoid arthritis. Ann Rheum Dis. 2018;77(1):48–54. https://doi.org/10.1136/annrheumdis-2017-211735 published Online First: 2017/09/08.

    Article  CAS  PubMed  Google Scholar 

  3. Panoulas VF, Metsios GS, Pace AV, et al. Hypertension in rheumatoid arthritis. Rheumatology (Oxford). 2008;47(9):1286–98. https://doi.org/10.1093/rheumatology/ken159 published Online First: 2008/05/10.

    Article  CAS  Google Scholar 

  4. Boyer JF, Gourraud PA, Cantagrel A, et al. Traditional cardiovascular risk factors in rheumatoid arthritis: a meta-analysis. Joint Bone Spine. 2011;78(2):179–83. https://doi.org/10.1016/j.jbspin.2010.07.016 published Online First: 2010/09/21.

    Article  PubMed  Google Scholar 

  5. Baghdadi LR, Woodman RJ, Shanahan EM, et al. The impact of traditional cardiovascular risk factors on cardiovascular outcomes in patients with rheumatoid arthritis: a systematic review and meta-analysis. PLoS One. 2015;10(2):e0117952. https://doi.org/10.1371/journal.pone.0117952 published Online First: 2015/02/18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Tiosano S, Yavne Y, Gendelman O, et al. Stroke among rheumatoid arthritis patients: does age matter? A real-life study. Neuroepidemiology. 2017;49(3–4):99–105. https://doi.org/10.1159/000481992 published Online First: 2017/11/15.

    Article  PubMed  Google Scholar 

  7. Meissner Y, Richter A, Manger B, et al. Serious adverse events and the risk of stroke in patients with rheumatoid arthritis: results from the German RABBIT cohort. Ann Rheum Dis. 2017;76(9):1583–90. https://doi.org/10.1136/annrheumdis-2017-211209 [published Online First: 2017/05/10].

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Holmqvist M, Gransmark E, Mantel A, et al. Occurrence and relative risk of stroke in incident and prevalent contemporary rheumatoid arthritis. Ann Rheum Dis. 2013;72(4):541–6. https://doi.org/10.1136/annrheumdis-2012-201387 [published Online First: 2012/05/23].

    Article  PubMed  Google Scholar 

  9. Nadareishvili Z, Michaud K, Hallenbeck JM, et al. Cardiovascular, rheumatologic, and pharmacologic predictors of stroke in patients with rheumatoid arthritis: a nested, case-control study. Arthritis Rheum. 2008;59(8):1090–6. https://doi.org/10.1002/art.23935 [published Online First: 2008/08/01].

    Article  PubMed  PubMed Central  Google Scholar 

  10. Wiseman SJ, Ralston SH, Wardlaw JM. Cerebrovascular disease in rheumatic diseases: a systematic review and meta-analysis. Stroke. 2016;47(4):943–50. https://doi.org/10.1161/STROKEAHA.115.012052 [published Online First: 2016/02/27].

    Article  PubMed  Google Scholar 

  11. Solomon DH, Karlson EW, Rimm EB, et al. Cardiovascular morbidity and mortality in women diagnosed with rheumatoid arthritis. Circulation. 2003;107(9):1303–7 published Online First: 2003/03/12.

    Article  Google Scholar 

  12. DeMarco VG, Aroor AR, Sowers JR. The pathophysiology of hypertension in patients with obesity. Nat Rev Endocrinol. 2014;10(6):364–76. https://doi.org/10.1038/nrendo.2014.44 published Online First: 2014/04/16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Skeoch S, Bruce IN. Atherosclerosis in rheumatoid arthritis: is it all about inflammation? Nat Rev Rheumatol. 2015;11(7):390–400. https://doi.org/10.1038/nrrheum.2015.40 published Online First: 2015/04/01.

    Article  CAS  PubMed  Google Scholar 

  14. Lauper K, Gabay C. Cardiovascular risk in patients with rheumatoid arthritis. Semin Immunopathol. 2017;39(4):447–59. https://doi.org/10.1007/s00281-017-0632-2 published Online First: 2017/04/30.

    Article  CAS  PubMed  Google Scholar 

  15. Low A, Mak E, Rowe JB, et al. Inflammation and cerebral small vessel disease: a systematic review. Ageing Res Rev. 2019:100916. https://doi.org/10.1016/j.arr.2019.100916 published Online First: 2019/06/11.

    Article  CAS  Google Scholar 

  16. Higashi Y, Sukhanov S, Anwar A, et al. Aging, atherosclerosis, and IGF-1. J Gerontol A Biol Sci Med Sci. 2012;67(6):626–39. https://doi.org/10.1093/gerona/gls102 [published Online First: 2012/04/12].

    Article  CAS  PubMed  Google Scholar 

  17. Ungvari Z, Csiszar A. The emerging role of IGF-1 deficiency in cardiovascular aging: recent advances. J Gerontol A Biol Sci Med Sci. 2012;67(6):599–610. https://doi.org/10.1093/gerona/gls072 [published Online First: 2012/03/28].

    Article  CAS  PubMed  Google Scholar 

  18. Clemmons DR. Metabolic actions of insulin-like growth factor-I in normal physiology and diabetes. Endocrinol Metab Clin N Am 2012;41(2):425–443, vii-viii. doi: https://doi.org/10.1016/j.ecl.2012.04.017 [published Online First: 2012/06/12].

  19. Bondanelli M, Ambrosio MR, degli Uberti EC. Pathogenesis and prevalence of hypertension in acromegaly. Pituitary. 2001;4(4):239–49 published Online First: 2002/12/28.

    Article  CAS  Google Scholar 

  20. Vitale G, Pivonello R, Auriemma RS, et al. Hypertension in acromegaly and in the normal population: prevalence and determinants. Clin Endocrinol. 2005;63(4):470–6. https://doi.org/10.1111/j.1365-2265.2005.02370.x published Online First: 2005/09/27.

    Article  Google Scholar 

  21. Conti E, Carrozza C, Capoluongo E, et al. Insulin-like growth factor-1 as a vascular protective factor. Circulation. 2004;110(15):2260–5. https://doi.org/10.1161/01.CIR.0000144309.87183.FB published Online First: 2004/10/13.

    Article  PubMed  Google Scholar 

  22. Schutte AE, Volpe M, Tocci G, et al. Revisiting the relationship between blood pressure and insulin-like growth factor-1. Hypertension. 2014;63(5):1070–7. https://doi.org/10.1161/HYPERTENSIONAHA.113.03057 published Online First: 2014/02/26.

    Article  CAS  PubMed  Google Scholar 

  23. Hunt KJ, Lukanova A, Rinaldi S, et al. A potential inverse association between insulin-like growth factor I and hypertension in a cross-sectional study. Ann Epidemiol. 2006;16(7):563–71. https://doi.org/10.1016/j.annepidem.2005.10.007 published Online First: 2006/01/25.

    Article  PubMed  Google Scholar 

  24. Colao A, Di Somma C, Cascella T, et al. Relationships between serum IGF1 levels, blood pressure, and glucose tolerance: an observational, exploratory study in 404 subjects. Eur J Endocrinol. 2008;159(4):389–97. https://doi.org/10.1530/EJE-08-0201 published Online First: 2008/07/08.

    Article  CAS  PubMed  Google Scholar 

  25. Denti L, Annoni V, Cattadori E, et al. Insulin-like growth factor 1 as a predictor of ischemic stroke outcome in the elderly. Am J Med. 2004;117(5):312–7. https://doi.org/10.1016/j.amjmed.2004.02.049 published Online First: 2004/09/01.

    Article  CAS  PubMed  Google Scholar 

  26. Parekh N, Roberts CB, Vadiveloo M, et al. Lifestyle, anthropometric, and obesity-related physiologic determinants of insulin-like growth factor-1 in the third National Health and Nutrition Examination Survey (1988-1994). Ann Epidemiol. 2010;20(3):182–93. https://doi.org/10.1016/j.annepidem.2009.11.008 published Online First: 2010/02/18.

    Article  PubMed  Google Scholar 

  27. Kaplan RC, McGinn AP, Pollak MN, et al. High insulinlike growth factor binding protein 1 level predicts incident congestive heart failure in the elderly. Am Heart J. 2008;155(6):1006–12. https://doi.org/10.1016/j.ahj.2007.12.031 published Online First: 2008/06/03.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Laughlin GA, Barrett-Connor E, Criqui MH, et al. The prospective association of serum insulin-like growth factor I (IGF-I) and IGF-binding protein-1 levels with all cause and cardiovascular disease mortality in older adults: the Rancho Bernardo Study. J Clin Endocrinol Metab. 2004;89(1):114–20. https://doi.org/10.1210/jc.2003-030967 published Online First: 2004/01/13.

    Article  CAS  PubMed  Google Scholar 

  29. Lawlor DA, Ebrahim S, Smith GD, et al. The association of insulin-like-growth factor 1 (IGF-1) with incident coronary heart disease in women: findings from the prospective British Women’s Heart and Health Study. Atherosclerosis. 2008;201(1):198–204. https://doi.org/10.1016/j.atherosclerosis.2007.12.061 published Online First: 2008/02/26.

    Article  CAS  PubMed  Google Scholar 

  30. Aberg ND, Aberg D, Jood K, et al. Altered levels of circulating insulin-like growth factor I (IGF-I) following ischemic stroke are associated with outcome - a prospective observational study. BMC Neurol. 2018;18(1):106. https://doi.org/10.1186/s12883-018-1107-3 [published Online First: 2018/08/08].

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Koshy PJ, Henderson N, Logan C, et al. Interleukin 17 induces cartilage collagen breakdown: novel synergistic effects in combination with proinflammatory cytokines. Ann Rheum Dis. 2002;61(8):704–13 published Online First: 2002/07/16.

    Article  CAS  Google Scholar 

  32. Bostrom EA, Svensson M, Andersson S, et al. Resistin and insulin/insulin-like growth factor signaling in rheumatoid arthritis. Arthritis Rheum. 2011;63(10):2894–904. https://doi.org/10.1002/art.30527 published Online First: 2011/07/09.

    Article  CAS  PubMed  Google Scholar 

  33. Erlandsson MC, Toyra Silfversward S, Nadali M, et al. IGF-1R signalling contributes to IL-6 production and T cell dependent inflammation in rheumatoid arthritis. Biochim Biophys Acta Mol basis Dis. 2017;1863(9):2158–70. https://doi.org/10.1016/j.bbadis.2017.06.002 [published Online First: 2017/06/07].

    Article  CAS  PubMed  Google Scholar 

  34. Laurberg TB, Ellingsen T, Thorsen J, et al. Insulin-like growth factor I receptor density on CD4+T-lymphocytes from active early steroid- and DMARD-naive rheumatoid arthritis patients is up-regulated and not influenced by 1 year of clinically effective treatment. Rheumatol Int. 2012;32(2):501–4. https://doi.org/10.1007/s00296-010-1683-z published Online First: 2011/01/20.

    Article  CAS  PubMed  Google Scholar 

  35. Andersson KME, Wasen C, Juzokaite L, et al. Inflammation in the hippocampus affects IGF1 receptor signaling and contributes to neurological sequelae in rheumatoid arthritis. Proc Natl Acad Sci U S A. 2018;115(51):E12063–E72. https://doi.org/10.1073/pnas.1810553115 [published Online First: 2018/12/05].

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Erlandsson MC, Doria Medina R, Toyra Silfversward S, et al. Smoking functions as a negative regulator of IGF1 and impairs adipokine network in patients with rheumatoid arthritis. Mediat Inflamm. 2016;2016:3082820. https://doi.org/10.1155/2016/3082820 published Online First: 2016/04/05.

    Article  CAS  Google Scholar 

  37. Jood K, Ladenvall C, Rosengren A, et al. Family history in ischemic stroke before 70 years of age: the Sahlgrenska Academy study on ischemic stroke. Stroke. 2005;36(7):1383–7. https://doi.org/10.1161/01.STR.0000169944.46025.09 [published Online First: 2005/06/04].

    Article  PubMed  Google Scholar 

  38. D'Agostino RB Sr, Vasan RS, Pencina MJ, et al. General cardiovascular risk profile for use in primary care: the Framingham heart study. Circulation. 2008;117(6):743–53. https://doi.org/10.1161/CIRCULATIONAHA.107.699579 [published Online First: 2008/01/24].

    Article  PubMed  Google Scholar 

  39. Aberg D, Jood K, Blomstrand C, et al. Serum IGF-I levels correlate to improvement of functional outcome after ischemic stroke. J Clin Endocrinol Metab. 2011;96(7):E1055–64. https://doi.org/10.1210/jc.2010-2802 [published Online First: 2011/04/22].

    Article  PubMed  Google Scholar 

  40. Nadali M, Pullerits R, Andersson KME, et al. High expression of STAT3 in subcutaneous adipose tissue associates with cardiovascular risk in women with rheumatoid arthritis. Int J Mol Sci. 2017;18(11). https://doi.org/10.3390/ijms18112410 published Online First: 2017/11/16.

    Article  Google Scholar 

  41. Maison P, Griffin S, Nicoue-Beglah M, et al. Impact of growth hormone (GH) treatment on cardiovascular risk factors in GH-deficient adults: a metaanalysis of blinded, randomized, placebo-controlled trials. J Clin Endocrinol Metab. 2004;89(5):2192–9. https://doi.org/10.1210/jc.2003-030840 published Online First: 2004/05/06.

    Article  CAS  PubMed  Google Scholar 

  42. Colao A. The GH-IGF-I axis and the cardiovascular system: clinical implications. Clin Endocrinol. 2008;69(3):347–58. https://doi.org/10.1111/j.1365-2265.2008.03292 x [published Online First: 2008/05/09].

    Article  CAS  Google Scholar 

  43. Gazzaruso C, Gola M, Karamouzis I, et al. Cardiovascular risk in adult patients with growth hormone (GH) deficiency and following substitution with GH--an update. J Clin Endocrinol Metab. 2014;99(1):18–29. https://doi.org/10.1210/jc.2013-2394 published Online First: 2013/11/13.

    Article  CAS  PubMed  Google Scholar 

  44. Caplice NM, DeVoe MC, Choi J, et al. Randomized placebo controlled trial evaluating the safety and efficacy of single low-dose intracoronary insulin-like growth factor following percutaneous coronary intervention in acute myocardial infarction (RESUS-AMI). Am Heart J. 2018;200:110–7. https://doi.org/10.1016/j.ahj.2018.03.018 published Online First: 2018/06/15.

    Article  CAS  PubMed  Google Scholar 

  45. Heinen A, Nederlof R, Panjwani P, et al. IGF1 treatment improves cardiac remodeling after infarction by targeting myeloid cells. Mol Ther. 2019;27(1):46–58. https://doi.org/10.1016/j.ymthe.2018.10.020 published Online First: 2018/12/12.

    Article  CAS  PubMed  Google Scholar 

  46. Sandoo A, Panoulas VF, Toms TE, et al. Anti-TNFalpha therapy may lead to blood pressure reductions through improved endothelium-dependent microvascular function in patients with rheumatoid arthritis. J Hum Hypertens. 2011;25(11):699–702. https://doi.org/10.1038/jhh.2011.36 published Online First: 2011/04/29.

    Article  CAS  PubMed  Google Scholar 

  47. Desai RJ, Solomon DH, Schneeweiss S, et al. Tumor necrosis factor-alpha inhibitor use and the risk of incident hypertension in patients with rheumatoid arthritis. Epidemiology. 2016;27(3):414–22. https://doi.org/10.1097/EDE.0000000000000446 published Online First: 2016/01/26.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Spartano NL, Stevenson MD, Xanthakis V, et al. Associations of objective physical activity with insulin sensitivity and circulating adipokine profile: the Framingham heart study. Clin Obes. 2017;7(2):59–69. https://doi.org/10.1111/cob.12177 published Online First: 2017/01/24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Bjersing JL, Erlandsson M, Bokarewa MI, et al. Exercise and obesity in fibromyalgia: beneficial roles of IGF-1 and resistin? Arthritis Res Ther. 2013;15(1):R34. https://doi.org/10.1186/ar4187 published Online First: 2013/03/01.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Siegel P, Tencza M, Apodaca B, et al. Effectiveness of occupational therapy interventions for adults with rheumatoid arthritis: a systematic review. Am J Occup Ther 2017;71(1):7101180050p1-710118050p11. doi: https://doi.org/10.5014/ajot.2017.023176 [published Online First: 2016/12/28].

Download references

Acknowledgments

The authors appreciate the assistance of Drs. Mats Dehlin, Lovisa Leifsdottir, and Jan Bjersing, Rheumatology Clinics, Sahlgrenska University Hospital, and Dan Norberg, Rheumatology Unit, Uddevalla Hospital, for clinical evaluation of the patients.

Funding

This work has been funded by grants from the Swedish Research Council (MB, 521-2014-02637, 521-2017-03025), the Medical Society of Göteborg (MB), the Swedish Association against Rheumatism (MB, R-566961, R-477321; RP), the King Gustaf V:s 80-year Foundation (MB), Professor Nanna Svartz Foundation (RP), the Torsten Söderberg’s Foundation (MB), Rune and Ulla Amlövs Trust, the Lundberg’s Foundation (MB), the University of Gothenburg, the Regional agreement on medical training and clinical research between the Western Götaland county council and the University of Gothenburg (LUA/ALF, MB).

Author information

Authors and Affiliations

  1. Department of Rheumatology and Inflammation Research, Institute of Medicine, the Sahlgrenska Academy at University of Gothenburg, Guldhedsgatan 10A, SE-41345, Gothenburg, Sweden

    Malin C. Erlandsson, Lovisa Lyngfelt, Caroline Wasén, Rachelle A. Espino, Sofia Töyrä Silfverswärd, Mitra Nadali, Karin M.E. Andersson, Rille Pullerits & Maria I. Bokarewa

  2. Rheumatology Clinic, the Sahlgrenska University Hospital, Gothenburg, Region of West Götaland, Sweden

    Malin C. Erlandsson, Mitra Nadali, Rille Pullerits & Maria I. Bokarewa

  3. Department of Internal Medicine, Institute of Medicine, the Sahlgrenska Academy at University of Gothenburg, Gothenburg, Sweden

    N. David Åberg

  4. Keele University, Keele, UK

    Rachelle A. Espino

  5. Department of Clinical Neuroscience, Institute of Neuroscience and Physiology, the Sahlgrenska Academy at University of Gothenburg, Gothenburg, Sweden

    Katharina Jood

  6. Department of Neurology, the Sahlgrenska University Hospital, Gothenburg, Sweden

    Katharina Jood

  7. Department of Clinical Immunology and Transfusion Medicine, the Sahlgrenska University Hospital, Gothenburg, Sweden

    Rille Pullerits

Authors
  1. Malin C. Erlandsson
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lovisa Lyngfelt
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. N. David Åberg
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Caroline Wasén
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Rachelle A. Espino
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Sofia Töyrä Silfverswärd
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Mitra Nadali
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Katharina Jood
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Karin M.E. Andersson
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Rille Pullerits
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Maria I. Bokarewa
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

MB conceived the study. MB, MN, RP, STS, LL, KJ, and DÅ collected patient material. ME, KA, CW, and RP did the laboratory analyses. ME, CW, DÅ, and MB did the statistical analysis. ME, CW, RP, DÅ, KJ, and MB prepared the manuscript. All authors have read and approved the final version of the manuscript.

Corresponding authors

Correspondence to Malin C. Erlandsson or Maria I. Bokarewa.

Ethics declarations

Ethics approval and consent to participate

The study protocol was reviewed and approved by the Ethical Review Boards of Gothenburg. All methods used in this study were carried out in accordance with relevant Swedish guidelines and regulations and following the Good Clinical Practice. The informed consent was obtained from all subjects.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Erlandsson, M.C., Lyngfelt, L., Åberg, N.D. et al. Low serum IGF1 is associated with hypertension and predicts early cardiovascular events in women with rheumatoid arthritis. BMC Med 17, 141 (2019). https://doi.org/10.1186/s12916-019-1374-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-019-1374-x

Keywords

BMC Medicine

ISSN: 1741-7015

Contact us