Skip to main content

The novel inflammatory biomarker GlycA and triglyceride-rich lipoproteins are associated with the presence of subclinical myocardial dysfunction in subjects with type 1 diabetes mellitus

Abstract

Background

Subjects with Type 1 diabetes mellitus (T1DM) have an increased incidence of heart failure (HF). Several pathophysiological mechanisms have been involved in its development. The aim of this study was to analyze the potential contribution of the advanced lipoprotein profile and plasma glycosylation (GlycA) to the presence of subclinical myocardial dysfunction in subjects with T1DM.

Methods

We included subjects from a Danish cohort of T1DM subjects (Thousand & 1 study) with either diastolic and/or systolic subclinical myocardial dysfunction, and a control group without myocardial dysfunction, matched by age, sex and HbA1c. All underwent a transthoracic echocardiogram and an advanced lipoprotein profile obtained by using the NMR-based Liposcale® test. GlycA NMR signal was also analyzed. Systolic dysfunction was defined as left ventricular ejection fraction ≤ 45% and diastolic dysfunction was considered as E/e′≥12 or E/e′ 8–12 + volume of the left atrium > 34 ml/m2. To identify a metabolic profile associated with the presence of subclinical myocardial dysfunction, a multivariate supervised model of classification based on least squares regression (PLS-DA regression) was performed.

Results

One-hundred forty-six subjects had diastolic dysfunction and 18 systolic dysfunction. Compared to the control group, patients with myocardial dysfunction had longer duration of diabetes (p = 0.005), and higher BMI (p = 0.013), serum NTproBNP concentration (p = 0.001), systolic blood pressure (p < 0.001), albuminuria (p < 0.001), and incidence of advanced retinopathy (p < 0.001). The supervised classification model identified a specific pattern associated with myocardial dysfunction, with a capacity to discriminate patients with myocardial dysfunction from controls. PLS-DA showed that triglyceride-rich lipoproteins (TGRLs), such as VLDL (total VLDL particles, large VLDL subclass and VLDL-TG content) and IDL (IDL cholesterol content), as well as the plasma concentration of GlycA, were associated with the presence of subclinical myocardial dysfunction.

Conclusion

Proatherogenic TGRLs and the proinflammatory biomarker Glyc A are strongly associated to myocardial dysfunction in T1DM. These findings suggest a pivotal role of TGRLs and systemic inflammation in the development of subclinical myocardial dysfunction in T1DM.

Background

Diabetes mellitus (DM) and heart failure (HF) are two multifaceted entities that involve high morbidity and mortality when both conditions coexist [1, 2]. The risk of HF is increased both in subjects with type 2 DM (T2DM) [3] and type 1 DM (T1DM) [4]. Indeed, DM is highly prevalent amongst patients with HF [5, 6], especially those with HF and preserved ejection fraction (HFpEF) [7, 8].

In population-based studies, the risk of HF in patients with diabetes (particularly T2DM) is significantly increased following adjustment for well-established HF risk factors [9]. The resulting specific form of cardiomyopathy is known as “diabetic cardiomyopathy” [10]. Although the concept of diabetic cardiomyopathy is often considered in individuals affected by T2DM, a metabolically-induced cardiomyopathy, independent of hypertension, nephropathy or ischemic heart disease is also evident in individuals with T1DM [11, 12].

Several pathophysiological mechanisms directly affecting the structure and function of the myocardium have been proposed to contribute to the development of diabetic cardiomyopathy. Among those described are hyperglycemia, hyperinsulinemia, inflammation and increased levels of circulating fatty acids (FAs) and triglycerides (TGs) [13,14,15,16,17]. Systemic inflammation plays a key role in HF etiopathogenesis [18, 19]. In that sense, GlycA has been described as a “composite biomarker of systemic inflammation” since its signal on nuclear magnetic resonance (NMR) spectra represents both the levels and degree of glycosylation of various acute phase proteins. GlycA NMR signal has been reported to be associated with increased risk of CV events, peripheral arterial disease, and mortality even after adjusting for other inflammatory markers [20, 21].

Lipotoxicity and cardiac lipid accumulation are other factors that have been related to the etiopathogenesis of diabetic cardiomyopathy [22]. In this line, myocardial metabolism studies have shown a reduced myocardial glucose uptake and an increased uptake of FAs in subjects with T1DM [23]. In T1DM, insulin deficiency promotes the mobilization of FAs from fat pads which results in an increased availability of excess FAs in different tissues, including the myocardium. When the capacity for storage and oxidation of mobilized FAs is exceeded, they can be transformed in other reactive species that further potentiates myocardial lipotoxicity. This cause of non-ischemic and non-hypertensive cardiomyopathy is often referred to as diabetic or “lipotoxic” cardiomyopathy.

Diabetic dyslipidaemia may also contribute to the diabetic myocardial dysfunction. Particularly, the excess flux of mobilized FAs to the liver promotes overproduction of TG-rich lipoproteins (TGRLs) and their remnants. Elevations in circulating TGRLs are frequently associated with increased concentrations of remnant cholesterol and with reduced high-density lipoproteins (HDL) cholesterol, and all contribute to the development of ischemic heart disease [24]. However, their contribution, if any, on non-ischemic cardiomyopathy remains poorly explored.

To the best of our knowledge, no previous studies in T1DM patients have studied the relationship between metabolic advanced profile with subclinical HF, defined as presence of impaired cardiac diastolic and/or systolic function without previous clinical manifestation of HF. Thus, the present study aims to analyze the contribution of inflammation (GlycA) and an advanced lipoprotein profile to the presence of subclinical myocardial dysfunction in a well-defined cohort of T1DM patients.

Methods

Study population

Our study population was selected from the Thousand & 1 study cohort study. This study was carried out at the Steno Diabetes Center in Copenhagen (SDCC) with cardiology examinations conducted at the Department of Cardiology, Copenhagen University Hospital Herlev-Gentofte. It was conducted between April 1st, 2010, and April 1st, 2012, and was based on a large cohort of 1.093 patients with T1DM without known heart disease. Patients were included if they were 18 years of age or older, attending the outpatient clinic at the SDCC, diagnosed with T1DM, without known heart disease (defined as any known HF, coronary artery disease [including previous myocardial infarction, stable angina, previous percutaneous coronary intervention, or coronary artery bypass surgery], atrial fibrillation or atrial flutter, left bundle branch block, congenital heart disease, pacemaker or implantable cardioverter defibrillator implantation), and if willing to participate in the study. The study population has been described in detail elsewhere [25]. Briefly, from the total study population included in the study, 15.5% (n = 169) of the participants had grossly abnormal systolic or diastolic function.

From the original cohort of Thousand & 1 study, we finally included a subgroup of 304 patients with T1DM, comprising 154 patients with myocardial dysfunction (the myocardial dysfunction group) and 150 controls. For the myocardial dysfunction group, we selected all T1DM patients who had either diastolic and/or systolic myocardial dysfunction and enough serum sample for conducting the NMR spectroscopy analysis. Additionally, we included 150 patients from the same cohort who had no echocardiographic alterations (the control group) matched by age, glycated hemoglobin (HbA1c) and gender.

Ethical consideration

The original study was performed in accordance with the second Helsinki declaration and approved by the regional ethics committee (H-3-2009-139) and the Danish Data Protection Agency (00934-Geh-2010-003). All subjects gave written informed consent.

Study visit

Prior to the echocardiographic examination, all patients received study information, signed the consent form and filled out a questionnaire with information about lifestyle factors, including smoking, exercise, alcohol consumption and cardiorespiratory symptoms. The use of cardiovascular treatments, such as lipid lowering medication (statins) and antihypertensive medication (beta blockers, calcium antagonists, diuretics, angiotensin-converting enzyme inhibitors or angiotensin II receptor antagonists), was also recorded. Blood pressure was measured in the supine position.

Echocardiogram

Echocardiography was performed with a General Electric, Vivid 7 Dimension imaging system device (GE Vingmed Ultrasound AS, Horten, Norway) with a 3.5-MHz transducer in accordance with the recommendations from the European Association of Echocardiography/American Society of Echocardiography [26]. Echocardiographic examinations were read and analysed using the General Electric EchoPAC software (BT11), recording three consecutive heart cycles. Left ventricle ejection fraction (LVEF) was determined by Simpson’s biplane method. Left atrial volume (LAV) was determined by the recommended biplane area-length method and indexed for body surface area. LV mass was determined by the linear method and indexed for body surface area.

Subclinical myocardial dysfunction was defined when patients had systolic and/or diastolic myocardial dysfunction on the echocardiogram in the absence of HF symptoms. Systolic dysfunction was defined as LVEF ≤ 45% determined by Simpson´s biplane method and diastolic dysfunction was considered if there was evidence of long-standing LV filling pressure defined as E/e′ >12 (where E is diastolic mitral early inflow velocity and e′ is pulsed-wave early diastolic tissue doppler velocity) or E/e′ 8–12 and LAV > 34 ml/m2.

Biochemistry

Information about biochemistry such as HbA1c, p-creatinine and albuminuric status was collected from electronic patient files at the SDCC from the ambulatory visit closest to study inclusion, which was maximally ± 4 months from inclusion. This information was collected after analyzing the echocardiography.

The urinary albumin excretion rate (UAER) was measured in 24-h sterile urine collections by enzyme immunoassay. Patients were categorized as normoalbuminuric if the UAER, in two out of three consecutive measurements, was < 30 mg/24 h, microalbuminuric if the UAER was between 30 and 300 mg/24 h, and macroalbuminuric if the UAER > 300 mg/24 h. HbA1c was measured by high-performance liquid chromatography (normal range: 21–46 mmol/mol [4.1–6.4%]; Variant; Bio-Rad Laboratories, Munich, Germany), and serum creatinine concentration was measured by an enzymatic method (Hitachi 912; Roche Diagnostics, Mannheim, Germany). The estimated glomerular filtration rate (eGFR) was calculated by the MDRD method.

Lipoprotein and glycoprotein analysis by NMR spectroscopy (advanced profile)

Serum samples were shipped on dry ice from the SDCC to the Biosfer Teslab facilities (Reus, Spain) for Liposcale® lipoprotein and glycoprotein analysis. Samples were kept at −80°C until the NMR analysis. 200 µl of serum was diluted with 50 µl deuterated water and 300 µl of 50 mM phosphate buffer solution at pH 7.4. 1H-NMR spectra were recorded at 306 K on a Bruker Avance III 600 spectrometer operating at a proton frequency of 600.20 MHz (14.1 T).

Lipoprotein analysis

Lipoprotein profiling was obtained by using the Liposcale® test (IVD-CE), a previously reported method based on a two-dimensional 1H-NMR diffusion-ordered spectroscopy (DOSY) approach for lipoprotein profile characterization including lipid content (cholesterol and triglyceride concentration), size and particle number of the main lipoprotein classes [27]. The methyl signal was deconvoluted by using 9 lorentzian functions to determine the lipid concentration of the large, medium and small subclasses of the main lipoprotein classes: very low-density lipoproteins (VLDL), low-density lipoproteins (LDL) and HDL, and their size associated diffusion coefficients. Then, the lipid concentration was combined with their associated particle volume in order to quantify the number of particles required to transport the measured lipid concentration of each lipoprotein subclass. Weighted average VLDL, LDL and HDL particle sizes were calculated from various subclass concentrations by summing the known diameter of each subclass multiplied by its relative percentage of subclass particle number. The variation coefficients for the particle numbers were between 2% and 4%, and for the particle sizes they were lower than 0.3%.

Glycoprotein analysis

The region of the 1H-NMR spectrum where the glycoproteins resonate (2.15–1.90 ppm) was deconvoluted using several analytical functions associated with specific plasmatic sugar-protein bonds according to a previously published procedure [28]. For each function, the total area (proportional to concentration), height, position and bandwidth were determined. The area of the specific glycoprotein signal defined as the GlycA NMR signal arose from the acetyl groups of N-acetylglucosamine and N-acetylgalactosamine bonded to plasmatic proteins [29]. Consistent with that, a larger GlycA area reflects a higher level of plasmatic glycosylation. The height-to-width (H/W) ratio of the GlycA signal, associated with its molecular aggregation state, was also reported. The H/W parameter, increased during inflammatory processes, reflects the sugar-protein bond flexibility, indicating glycosylation in accessible regions of the proteins. Height is defined as the difference from baseline to the maximum of the corresponding NMR peaks and the width value corresponds to the peak width at half height.

Statistical analysis

We used widely described chemometric methods [30] in order to identify a specific lipoprotein/glycoprotein profile associated with myocardial dysfunction. Briefly, multivariate statistical analyses were computed in MATLAB, Ver. 7.10.0 using PLS-Toolbox, Ver. 5.2.2 (Eigenvector Research Inc., Manson, WA, United States) after application of a genetic algorithm (GA) for variable selection to optimize the predictive ability of the model. Partial least squares discriminant analysis (PLS-DA) models were used as a supervised classification method between the study groups. This is a well established and widely used method in chemometrics- and metabolomics-based analyses when a large number of variables are used to classify two-stage conditions. PLS-DA relates the X matrix (input data including NMR derived parameters, clinical and anthropometric variables), and the Y matrix (presence or absence of myocardial disfunction) to find the maximum discrimination between classes (study groups). The PLS-DA method reduces the dimensionality of the initial dataset (X matrix), creating a new multidimensional dataset for each individual maximizing the total variance of data using just a few components (the latent scores), obtained from the specific contribution of each variable (loadings) to the new multidimensional axis [31].

On the other hand, to avoid overfitting and correct the multiple testing effects, we auto scaled and cross-validated by the permutation Venetian Blinds cross validation method by using 10 splits of the input data. The area under the curve (AUC) was used to evaluate the capacity of the NMR and clinical variables to distinguish between the two groups (with and without subclinical myocardial dysfunction).

The following subset of 19 clinical and biochemical variables was selected after a GA approach: diabetes duration, body mass index (BMI), hemoglobin, intermediate-density lipoproteins (IDL) cholesterol content (IDL-C), LDL cholesterol content (LDL-C), VLDL TG content (VLDL-TG), total VLDL particles (VLDL Particles), Large VLDL, total LDL particles (LDL Particles), total HDL particles (HDL Particles); Large HDL, Medium HDL, LDL size (LDL-Z), HDL size (HDL-Z), GlycA area, H/W GlycA ratio (H/W GlycA), HDL-TG/HDL-C ratio (HDL ratio), Small VLDL/total VLDL particles ratio (% Small VLDL) and Small LDL/total LDL particles ratio (% Small LDL).

Finally, we assessed the discrimination capacity of NMR biomarkers to predict the presence of myocardial dysfunction when these variables were added to the model including only traditional clinical variables (age, sex, eGFR, NTproBNP, BMI, diabetes duration and systolic blood pressure) by using logistic regression.

Results

Clinical characteristics of our study population

Our study population comprised a total of 304 T1DM subjects (53.6% women, 50.7% subjects with myocardial dysfunction) with a median [min;max] age of 62.1 [22.5;87.0] years at inclusion, median diabetes duration of 33 years and a median HbA1c value of 8.0%. Clinical characteristics by group are shown in Table 1. Briefly, there were no differences regarding gender, age, HbA1c value, smoking habit, total cholesterol, LDL cholesterol and TGs values between the groups. From a total of 154 subjects with myocardial dysfunction, 146 (94.8%) and 18 (11.7%) had diastolic and systolic dysfunction, respectively. Subjects with myocardial dysfunction showed, in comparison with controls, a longer diabetes duration (35.1 ± 14.9 years vs. 30.1 ± 15.5 years; p = 0.005), a higher BMI (26.1 ± 3.9 kg/m2 vs. 25.0 ± 3.7 kg/m2; p = 0.013), a higher systolic blood pressure (143 mmHg vs. 136 mmHg; p < 0.001) and a lower eGFR (75.4 ± 26.2 mL/min/1.73m2 vs. 83.7 ± 21.0 mL/min/1.73m2; p = 0.003). Differences in echocardiographic parameters and serum cardiac biomarkers of HF were also observed between the two groups. In summary, subjects with myocardial dysfunction had lower mean (SD) ejection fraction (55.8 ± 7.5% vs. 58.6 ± 5.2%, p < 0.001), higher LAV (34.1 ± 7.6 ml/m2 vs. 28.9 ± 5.5 ml/m2, p < 0.001) and higher median [25th;75th] concentrations of N-terminal fragment of pro-B-type natriuretic peptide (NTproBNP) (352 [163;591] pg/mL vs. 249 [116;449] pg/mL, p = 0.001), and mid-regional pro atrial natriuretic peptide (MRproANP) (100 [69.2;144] pmol/L vs. 82.3 [56.4;112] pmol/L, p < 0.001). Also, the percentage of subjects with E/e′ >8 on echocardiography was higher in subjects with myocardial dysfunction compared with those without myocardial dysfunction (94.8% vs. 38.7%, p < 0.001). Additionally, the use of statin and antihypertensive therapies was more frequent in subjects with myocardial dysfunction than in those without (p = 0.039 and p < 0.001 respectively). Finally, subjects with myocardial dysfunction were more likely to have advanced stages of retinopathy and albuminuria (p < 0.001 for both comparisons) than those without.

Table 1 Descriptive analysis of clinical variables by group

Lipoprotein and glycoprotein NMR advanced profile in subjects with myocardial dysfunction and controls

Results from the lipoprotein and glycoprotein profile are shown in Table 2. Most remarkably, compared with controls, subjects with myocardial dysfunction presented a more proatherogenic and pro-inflammatory profile, associated with increased IDL cholesterol and TG content (p = 0.004 and 0.003 respectively) and a greater GlycA area (p < 0.001). The LDL lipid composition, defined as LDL-TG / LDL-C, was also increased in the myocardial dysfunction group, reflecting increased remnant levels and neutral lipid heteroexchange (TGs and cholesteryl esters) from TG-rich particles to LDL particles, and vice versa, by the cholesteryl ester transfer protein.

Table 2 Lipoprotein and glycoprotein NMR advanced profile in cases and controls

Lipoprotein and glycoprotein signature associated with the presence of myocardial dysfunction

We used a multivariant classification approach based on PLS-DA in order to identify a specific lipoprotein and glycoprotein profile associated with myocardial dysfunction in subjects with T1DM. Figure 1 shows the good performance of the classification method by using a ROC curve analysis, and the contribution of each variable to the model, summarized in the loadings plot of the principal two latent multidimensional variables (LV1 and LV2).

Fig. 1
figure 1

 A: Estimated and cross-validated ROC curve of the PLS-DA classification model. B: Contribution of each variable to the multidimensional latent variables LV1 and LV2. PLS-DA relates the X matrix (input data including NMR derived parameters, clinical and anthropometric variables) and the Y matrix (presence or absence of myocardial disfunction) to find the maximum discrimination between both groups. AUROC: area under the ROC curve, PLS-DA: partial least squares discriminant analysis, LV: latent variable, NMR: nuclear magnetic resonance

This supervised classification model succeeded in identifying a specific pattern associated with myocardial dysfunction, with a capacity to discriminate diabetic patients with myocardial dysfunction from the rest of the individuals without CVD, in a modest but significant way in relation to fortuitously obtaining the correct classification (area under the ROC curve 0.63, p < 1.1801e-012).

In that sense, we used the PLS-DA classification approach to investigate which variables discriminated best between subjects with myocardial dysfunction and controls (non-myocardial dysfunction subjects). LV1 showed that the variables with the most significant contribution to explain the presence of myocardial dysfunction in T1DM subjects beyond well-established risk factors included: NMR-determined GlycA area and H/W GlycA ratio as well as TGRLs, such as VLDL (total VLDL particles, large VLDL subclass and VLDL-TG content) and IDL (IDL cholesterol content). In the second latent variable (i.e. after the first discriminating analysis based on inflammatory parameters and VLDL and IDL related parameters was accomplished), other clinical variables including diabetes duration and LDL related characteristics (LDL-Z) became discriminative, whereas VLDL related variables were no longer discriminant. Furthermore, we compared the predicted probability of the presence of myocardial dysfunction when NMR variables were added together with traditional risk factors by using logistic regressions. Supplementary Fig. 1 shows the predicted probability of myocardial dysfunction in a model (model 1) including only classical risk variables (age, sex, eGFR, NTproBNP, BMI, diabetes duration and systolic blood pressure > 140mmHg). We built a second model (model 2) including both classical variables and the NMR-assessed biomarkers. The inclusion of the NMR variables in this second model significantly increased the area under the ROC curve from 0.62 [0.56–0.68] to 0.67 [0.61–0.73], with a net reclassification analysis improvement considering NMR-assessed parameters of 21%.

Discussion

To our knowledge, this is the first study to date evaluating the association between a metabolic advanced profile with subclinical myocardial dysfunction assessed with echocardiography in subjects with T1DM. The analysis has been performed in a well characterized large cohort of T1DM patients without prior CVD. Our data showed that the inflammation biomarker GlycA as well as VLDL-related variables (total VLDL particles, large VLDL subclass and VLDL-TG content) and IDL-related variables (IDL cholesterol content) were associated with the presence of myocardial dysfunction in T1DM subjects.

Although the underlying mechanisms explaining the myocardial dysfunction remain still obscure, accumulating evidence supports the role of chronic inflammation in its manifestation and progress to HF [3233]. In support of this, our data suggest that the circulating levels of GlycA are related to myocardial dysfunction in T1DM subjects. Our finding is consistent with the significance of circulating GlycA, as it integrates in a single metrics the plasma concentrations [34] and glycosylation states [29, 35] of several abundant acute-phase inflammatory proteins. Notably, GlycA exhibits lower intraindividual variability and greater analytical precision than other inflammatory markers [34], thereby suggesting a role as a potential novel biomarker linked to systemic inflammation and CVD assessment [36,37,38]. In line with this, the serum levels of GlycA are directly associated with those of C-reactive protein (CRP) [39], which is a well-established diagnostic biomarker of inflammation used in clinical practical clinical. Noteworthy, the advantage offered by GlycA over CRP is that it may integrate multiple inflammatory pathways by capturing the global signal of several proteins and, therefore, better captures the degree of systemic inflammation [40]. In addition, the measurement of GlycA is associated with higher reliability and lower intra-individual variability than CRP because it is detected similarly in both serum and plasma samples, in fasting and non-fasting states, after short or long-term storage [34], and also after repeated intraindividual determinations over time [41].

A main finding of our study was that GlycA values allowed identification of T1DM subjects with and without subclinical myocardial dysfunction. Our data are in line with other study that support the concept that GlycA is associated with an increased risk of HF, particularly HFpEF, which is independent of traditional CVD risk factors and other inflammatory markers [42]. To our knowledge, data about GlycA concentrations in subjects with T1DM are lacking and for the first time our results regarding GlycA further suggest a role for inflammatory mechanisms in the pathogenesis of myocardial dysfunction in these patients.

Our data also showed that some characteristics of VLDL-related particles, together with GlycA, could help to better discriminate between T1DM subjects with and without subclinical myocardial dysfunction. Regarding T1DM, the advanced characteristics of VLDL has been positively and independently associated with arterial stiffness [43]. Likewise, similar associations have also been described in T1DM women with previous pre-eclampsia with carotid atherosclerosis [44].

Additionally, the role of TGRLs in HF has been previously established. In two prospective studies of 113,554 individuals from the general population in Denmark, a higher risk of HF for stepwise higher non-fasting TGs (i.e., TGRLs, which also include the chylomicron remnants) was found [45].

Although evidence supporting a relationship between cholesterol remnants and non-ischemic causes of HF is limited, the induction of VLDL receptor (VLDLR) abundance and TGRL uptake by cardiomyocytes in an experimental postprandial setting [46] suggest a role for VLDLR in atrial cardiomyopathy and ventricular dysfunction [47, 48]. Indeed, an association of atrial fibrillation with abnormal VLDL-related lipid metabolism has been suggested by some authors [48] and an increased risk of developing atrial fibrillation among patients with T1DM compared with the general population has been reported [49].

According to previous studies, the circulating baseline levels of GlycA exhibit significant positive associations with incident CVD event rates and associated mortality [37]. Remarkably, such correlations remain significant even after adjusting for several other established CVD risk factors. In our study, the discriminating capability of our PLS-DA model was not improved when including either the 24 h-urinary albumin excretion rate or NTproBNP as input variables, maintaining GlycA as the best discriminating variable for the presence of myocardial disfunction in subjects with T1DM (data not shown).

The multivariant approach used in the present study by using the PLS-DA methodology helped in identifying concomitant lipoprotein characteristics (TGRLs) that, together with the proinflammatory GlycA, could mainly account for myocardial dysfunction in T1DM subjects on top of the commonly used traditional and well-established risk factors, including BMI, gender, HbA1c or diabetes duration. Inflammation showed the higher discriminant ability distinguishing myocardial dysfunction among T1DM patients, with the GlycA area being the highest contributor to the first multidimensional latent variable (LV1). Of note, diabetes duration helped to discriminate myocardial dysfunction from non-myocardial dysfunction T1DM patients only once the inflammatory signature was considered, as shown in the second latent variable (LV2) of the multivariable classification analysis. Furthermore, this second multidimensional latent variable showed a loss of association between isolated increased TG-rich parameters and myocardial dysfunction. Therefore, our data revealed an interaction between GlycA and TGRLs in T1DM subjects with myocardial dysfunction. This association was lost in those patients without myocardial dysfunction.

Our study had some limitations that should be taken into account. We matched individuals with and without subclinical myocardial dysfunction by age, sex and glycemic control. This may have resulted in a biased comparison because, according to original Thousand & 1 study results, patients with myocardial dysfunction were likely to be older and have longer diabetes duration than subjects with normal function. Therefore, by selecting subjects who were similar for these criteria, at least one group may be unrepresentative of the population from which it came, and the abnormalities observed in the NMR profile may be different from a truly representative population. Furthermore, considering the condition of hypertension as a risk factor for myocardial dysfunction, is important to emphasize that the use of antihypertensive therapies was significantly more frequent in the group of subjects with myocardial dysfunction. Nevertheless, due to the fact that these drugs can be used both as an antihypertensive treatment as well as a nephroprotective intervention, we consider this as a limitation to assess the condition of hypertension associated with myocardial dysfunction. Additionally, the small number of subjects who presented systolic dysfunction within subjects with myocardial dysfunction (n = 18, 11.7%) implies that the results obtained in the present study should be primordially referred to diastolic dysfunction. Finally, since this study is cross-sectional we cannot infer causality between the metabolic profile and presence of myocardial dysfunction. Possible strengths should also be mentioned, since the current population is selected from a large study of ambulatory T1DM patients without previous known heart disease. This cohort had the advantage of a relatively low HbA1c and may therefore be at a lower risk of HF than typical T1DM patients and that may potentially have strengthened our findings.

Conclusion

The present study uncovered associations between subclinical myocardial dysfunction and advanced NMR metabolic characteristics in patients with T1DM that were hidden in conventional analyses. The GlycA area and the H/W GlycA ratio, as well as TGRLs (VLDL and IDL) related variables, were revealed as strong contributors of subclinical myocardial dysfunction. According to the aforementioned results, we propose a pivotal role of TGRLs characteristics and systemic inflammation reflected by the GlycA biomarker in subclinical myocardial dysfunction in T1DM patients.

Estimation of the presence of subclinical myocardial dysfunction among T1DM patients and its relationship with lipoprotein and glycoprotein characteristics revealed in this study still deserves more attention in the future. Therefore, further studies will be required to clarify the potential clinical applications of these findings as well as to investigate their biological basis.

Data availability

Primary material is held by the authors.

Abbreviations

DM:

Diabetes mellitus.

HF:

Heart failure.

T2DM:

Type 2 diabetes mellitus.

T1DM:

Type 1 diabetes mellitus.

HFpEF:

Heart failure with preserved ejection fraction.

CVD:

Cardiovascular disease.

FAs:

Fatty acids.

TGs:

Triglycerides.

NMR:

Nuclear magnetic resonance.

TGRLs:

Triglycerides-rich lipoproteins.

HDL:

High density lipoproteins.

SDCC:

Steno diabetes center in Copenhagen.

HbA1c:

Glycated hemoglobin.

LVEF:

Left ventricle ejection fraction.

LAV:

Left atrial volume.

UAER:

Urinary albumin excretion rate.

eGFR:

Estimated glomerular filtration rate.

DOSY:

Diffusion-ordered spectroscopy.

VLDL:

Very low-density lipoproteins.

LDL:

Low-density lipoproteins.

GA:

Genetic algorithm.

PLS-DA:

Partial least square discriminant analysis.

AUC:

Area under the curve.

BMI:

Body mass index.

IDL:

Intermediate-density lipoprotein.

NTproBNP:

N-terminal fragment of pro-B-type natriuretic peptide.

MRproANP:

Mid-regional pro atrial natriuretic peptide.

CRP:

C-reactive protein.

VLDLR:

VLDL receptor.

hsCRP:

High-sensitivity CRP.

References

  1. Formiga F, Camafort M, Carrasco Sánchez FJ. Heart failure and diabetes: The confrontation of two major epidemics of the 21st century. Rev Clin Esp. 2020 Mar;220(2):135–8.

    Article  CAS  PubMed  Google Scholar 

  2. Lindman BR, Dávila-Román VG, Mann DL, McNulty S, Semigran MJ, Lewis GD, et al. Cardiovascular phenotype in HFpEF patients with or without diabetes: a RELAX trial ancillary study. J Am Coll Cardiol. 2014;64(6):541–9.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Shah AD, Langenberg C, Rapsomaniki E, Denaxas S, Pujades-Rodriguez M, Gale CP, et al. Type 2 diabetes and incidence of a wide range of cardiovascular diseases: a cohort study in 1?9 million people. Lancet. 2015;26(385 Suppl 1):86.

    Article  Google Scholar 

  4. Avogaro A, Azzolina D, Fadini GP, Baldi I. Incidence of heart failure in patients with type 1 diabetes: a systematic review of observational studies. J Endocrinol Invest. 2021;44(4):745–53.

    Article  CAS  PubMed  Google Scholar 

  5. Nichols GA, Hillier TA, Erbey JR, Brown JB. Congestive heart failure in type 2 diabetes: prevalence, incidence, and risk factors. Diabetes Care. 2001;24:1614–9.

    Article  CAS  PubMed  Google Scholar 

  6. From AM, Leibson CL, Bursi F, Redfield MM, Weston SA, Jacobsen SJ, et al. Diabetes in heart failure: prevalence and impact on outcome in the population. Am J Med. 2006;119:591–9.

    Article  PubMed  Google Scholar 

  7. McHugh K, DeVore AD, Wu J, Matsouaka RA, Fonarow GC, Heidenreich PA, et al. Heart failure with preserved ejection fraction and diabetes: JACC state-of-the-art review. J Am Coll Cardiol. 2019;73(5):602–11.

    Article  PubMed  Google Scholar 

  8. Fonarow GC, Stough WG, Abraham WT, Albert NM, Gheorghiade M, Greenberg BH, et al. Characteristics, treatments, and outcomes of patients with preserved systolic function hospitalized for heart failure: a report from the OPTIMIZE-HF Registry. J Am Coll Cardiol. 2007;50(8):768–77.

    Article  PubMed  Google Scholar 

  9. Bugger H, Abel ED. Molecular mechanisms of diabetic cardiomyopathy. Diabetologia. 2014;57(4):660–71.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Rydén L, Grant PJ, Anker SD, Berne C, Cosentino F, Danchin N, et al. ESC Guidelines on diabetes, pre-diabetes, and cardiovascular diseases developed in collaboration with the EASD: the Task Force on diabetes, pre-diabetes, and cardiovascular diseases of the European Society of Cardiology (ESC) and developed in collaboration with the European Association for the Study of Diabetes (EASD). Eur Heart J. 2013;34(39):3035–87.

    Article  PubMed  Google Scholar 

  11. Gøtzsche O, Darwish A, Gøtzsche L, Hansen LP, Sørensen KE. Incipient cardiomyopathy in young insulin-dependent diabetic patients: a seven-year prospective Doppler echocardiographic study. Diabet Med. 1996;13(9):834–40.

    Article  PubMed  Google Scholar 

  12. Rosengren A, Vestberg D, Svensson AM, Kosiborod M, Clements M, Rawshani A, et al. Long-term excess risk of heart failure in people with type 1 diabetes: a prospective case-control study. Lancet Diabetes Endocrinol. 2015;3(11):876–85.

    Article  PubMed  Google Scholar 

  13. Evangelista I, Nuti R, Picchioni T, Dotta F, PalazzuoliInt A. Molecular Dysfunction and Phenotypic Derangement in Diabetic Cardiomyopathy. J Mol Sci. 2019;20(13):3264.

    Article  CAS  Google Scholar 

  14. Borghetti G, von Lewinski D, Eaton DM, Sourij H, Houser SR, Wallner M. Diabetic Cardiomyopathy: Current and Future Therapies. Beyond Glycemic Control. Front Physiol. 2018;9:1514.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Tan Y, Zhang Z, Zheng C, Wintergerst KA, Keller BB, Cai L. Mechanisms of diabetic cardiomyopathy and potential therapeutic strategies: preclinical and clinical evidence. Nat Rev Cardiol. 2020;17:585–607.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Spillmann F, Linthout SV, Tschöpe C. Cardiac effects of HDL and its components on diabetic cardiomyopathy. Endocr Metab Immune Disord Drug Targets. 2012;12(2):132–47.

    Article  CAS  PubMed  Google Scholar 

  17. Xu X, Luo Z, He Y, Shan J, Guo J, Li J. Application of untargeted lipidomics based on UHPLC-high resolution tandem MS analysis to profile the lipid metabolic disturbances in the heart of diabetic cardiomyopathy mice. J Pharm Biomed Anal. 2020;190:113525.

    Article  CAS  PubMed  Google Scholar 

  18. Shirazi LF, Bissett J, Romeo F, Mehta JL. Role of inflammation in heart failure. Curr Atheroscler Rep. 2017;19(6):27.

    Article  PubMed  Google Scholar 

  19. Murphy SP, Kakkar R, McCarthy CP, Januzzi JL Jr. Inflammation in heart failure: JACC state-of-the-art review. J Am Coll Cardiol. 2020;75(11):1324–40.

    Article  PubMed  Google Scholar 

  20. Duprez DA, Otvos J, Sanchez OA, Mackey RH, Tracy R, Jacobs DR Jr. Comparison of the predictive value of GlycA and other biomarkers of inflammation for total death, incident cardiovascular events, noncardiovascular and noncancer inflammatory-related events, and total cancer events. Clin Chem. 2016;62(7):1020–31.

    Article  CAS  PubMed  Google Scholar 

  21. Akinkuolie AO, Buring JE, Ridker PM, Mora S. A novel protein glycan biomarker and future cardiovascular disease events. J Am Heart Assoc. 2014;3(5):e001221.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Ritchie RH, Zerenturk EJ, Prakoso D, Calkin AC. Lipid metabolism and its implications for type 1 diabetes-associated cardiomyopathy. J Mol Endocrinol. 2017;58(4):R225–40.

    Article  CAS  PubMed  Google Scholar 

  23. Herrero P, Peterson LR, McGill JB, Matthew S, Lesniak D, Dence C, et al. Increased myocardial fatty acid metabolism in patients with type 1 diabetes mellitus. JACC. 2006;47(3):598–604.

    Article  CAS  PubMed  Google Scholar 

  24. Varbo A, Nordestgaard BG. Remnant cholesterol and ischemic heart disease. Curr Opin Lipidol. 2014;25(4):266–73.

    Article  CAS  PubMed  Google Scholar 

  25. Jensen MT, Sogaard P, Andersen HU, Bech J, Hansen TF, Galatius S, et al. Prevalence of systolic and diastolic dysfunction in patients with type 1 diabetes without known heart disease: the Thousand & 1 Study. Diabetologia. 2014;57(4):672–80.

    Article  CAS  PubMed  Google Scholar 

  26. Lang RM, Bierig M, Devereux RB, Flachskampf FA, Foster E, Pellikka PA, et al. Recommendations for chamber quantification: a report from the American Society of Echocardiography’s Guidelines and Standards Committee and the Chamber Quantification Writing Group, developed in conjunction with the European Association of Echocardiography, a branch of the European Society of Cardiology. J Am Soc Echocardiogr. 2005;18(12):1440–63.

    Article  PubMed  Google Scholar 

  27. Mallol R, Amigó N, Rodríguez MA, Heras M, Vinaixa M, Plana N, et al. Liposcale: a novel advanced lipoprotein test based on 2D diffusion-ordered 1H NMR spectroscopy. J Lipid Res. 2015;56(3):737–46.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Fuertes-Martín R, Taverner D, Vallvé JC, Paredes S, Masana L, Correig X, et al. Characterization of 1H NMR Plasma Glycoproteins as a New Strategy To Identify Inflammatory Patterns in Rheumatoid Arthritis. J Proteome Res. 2018;17(11):3730–9.

    Article  PubMed  Google Scholar 

  29. Fuertes-Martín R, Moncayo S, Insenser M, Martínez-García MA, Luque-Ramírez M, Amigó Grau N, et al. Glycoprotein A and B Height-to-Width Ratios as Obesity-Independent Novel Biomarkers of Low-Grade Chronic Inflammation in Women with Polycystic Ovary Syndrome (PCOS). J Proteome Res. 2019;18(11):4038–4045.

  30. Bro R, Smilde AK. Principal component analysis. Anal Methods. 2014;6:2812–31.

    Article  CAS  Google Scholar 

  31. Lee LC, Liong CY, Jemain AA. Partial least squares-discriminant analysis (PLS-DA) for classification of high-dimensional (HD) data: a review of contemporary practice strategies and knowledge gaps. Analyst. 2018;143(15):3526–39.

    Article  CAS  PubMed  Google Scholar 

  32. Riehle C, Bauersachs J. Of mice and men: models and mechanisms of diabetic cardiomyopathy. Basic Res Cardiol. 2018;114(1):2.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Riehle C, Bauersachs J. Key inflammatory mechanisms underlying heart failure. Herz. 2019:96–106.

  34. Otvos JD, Shalaurova I, Wolak-Dinsmore J, Connelly MA, Mackey RH, Stein JH, et al. GlycA: a composite nuclear magnetic resonance biomarker of systemic inflammation. Clin Chem. 2015;61(5):714–23.

    Article  CAS  PubMed  Google Scholar 

  35. Fuertes-Martín R, Correig X, Vallvé JC, Amigó N. Title: human serum/plasma glycoprotein analysis by 1H-NMR, an emerging method of inflammatory assessment. J Clin Med. 2020;9(2):354.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Ballout RA, Remaley AT. GlycA: a new biomarker for systemic inflammation and cardiovascular disease (CVD) risk assessment. J Lab Precis Med. 2020;5:17.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Rodríguez-Carrio J, Alperi-López M, López P, Pérez-Álvarez AI, Gil-Serret M, Amigó N, et al. GlycA levels during the earliest stages of rheumatoid arthritis: potential use as a biomarker of subclinical cardiovascular disease. J Clin Med. 2020;9(8):2472.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Carmona-Maurici J, Amigó N, Cuello E, Bermúdez M, Baena-Fustegueras JA, Peinado-Onsurbe J. Bariatric surgery decreases oxidative stress and protein glycosylation in patients with morbid obesity. Eur J Clin Invest. 2020;50(11):e13320.

    Article  CAS  PubMed  Google Scholar 

  39. Mehta NN, Dey AK, Maddineni R, Kraus WE, Huffman KM. GlycA measured by NMR spectroscopy is associated with disease activity and cardiovascular disease risk in chronic inflammatory diseases. Am J Prev Cardiol. 2020;4:100120.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Connelly MA, Otvos JD, Shalaurova I, Playford MP, Mehta NN. GlycA, a novel biomarker of systemic inflammation and cardiovascular disease risk. J Transl Med. 2017;15:219.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Moncayo S, Insenser M, Martínez-García MA, Fuertes-Martín R, Amigó-Grau N, Álvarez-Blasco F, et al. Acute-phase glycoprotein profile responses to different oral macronutrient challenges: Influence of sex, functional hyperandrogenism and obesity. Clin Nutr. 2021;40(3):1241–6.

    Article  CAS  PubMed  Google Scholar 

  42. Jang S, Ogunmoroti O, Ndumele CE, Zhao D, Rao VN, Fashanu OE, et al. Association of the novel inflammatory marker GlycA and incident heart failure and its subtypes of preserved and reduced ejection fraction: the multi-ethnic study of atherosclerosis. Circ Heart Fail. 2020;13(8):e007067.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Llauradó G, Amigó N, Cano A, Ballesta S, Albert L, Mazarico I, et al. Specific nuclear magnetic resonance lipoprotein subclass profiles and central arterial stiffness in Type 1 diabetes mellitus: a case control study. J Clin Med. 2019;8(11):1875.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Amor AJ, Vinagre I, Valverde M, Urquizu X, Meler E, López E, et al. Nuclear magnetic resonance lipoproteins are associated with carotid atherosclerosis in type 1 diabetes and pre-eclampsia. Diabetes Metab Res Rev. 2021;37(1):e3362.

    Article  CAS  PubMed  Google Scholar 

  45. Varbo A, Nordestgaard BG. Nonfasting triglycerides, low-density lipoprotein cholesterol, and heart failure risk: two cohort studies of 113 554 individuals. Arterioscler Thromb Vasc Biol. 2018;38(2):464–72.

    Article  CAS  PubMed  Google Scholar 

  46. Nakajima K, Tokita Y, Tanaka A, Takahashi S. The VLDL receptor plays a key role in the metabolism of postprandial remnant lipoproteins. Clin Chim Acta. 2019;495:382–93.

    Article  CAS  PubMed  Google Scholar 

  47. Perman JC, Boström P, Lindbom M, Lidberg U, StÅhlman M, Hägg D, et al. The VLDL receptor promotes lipotoxicity and increases mortality in mice following an acute myocardial infarction. J Clin Invest. 2011;121(7):2625–40.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Lee HC, Lin YH. The Pathogenic Role of Very Low Density Lipoprotein on Atrial Remodeling in the Metabolic Syndrome. Int J Mol Sci. 2020;21(3):891.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Dahlqvist S, Rosengren A, Gudbjörnsdottir S, Pivodic A, Wedel H, Kosiborod M, et al. Risk of atrial fibrillation in people with type 1 diabetes compared with matched controls from the general population: a prospective case-control study. Lancet Diabetes Endocrinol. 2017;5(10):799–807.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We want to particularly acknowledge the participants of the Thousand & 1 study and all the collaboration provided from SDCC.

Funding

The support for the original Thousand & 1 Study has been provided by the European Foundation for the Study of Diabetes/Pfizer European Programme 2010 for Research into Cardiovascular Risk Reduction in Patients with Diabetes; and The Danish Heart Foundation (number 12-04-R90-A3840-22725). This work was supported by Grants from Fundació La MARATÓ de TV3 (201602- 30-31 to NAL and JJ), Instituto de Salud Carlos III (PI17- 00232 to JJ, PI17-01362 to NAL) and by CIBER -Consorcio Centro de Investigación Biomédica en Red- (CB15/00071), Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación, Spain, with co-funding from the European Regional Development Fund (ERDF). JJ is supported by funds provided by ISCIII (Grant CPII18/00004, Miguel Servet II program). This research was partially funded by the Agencia de Gestió d’Ajuts Universitaris I de Recerca AGAUR, grant number 2020-PANDE-00083.

Author information

Authors and Affiliations

Authors

Contributions

CPJ, JJ, PR, DM and NAL contributed to the study design, data analysis and interpretation, and the writing of the manuscript. HUA, PR and MTJe contributed to the acquisition, analysis and interpretation of data. EC, MTJu, NAM and TSA contributed to data interpretation and discussion. PR, DM and NAL contributed to the study coordination, conduct of the study and writing of the manuscript. All authors critically reviewed the manuscript and agreed with the final version of the manuscript.

Corresponding author

Correspondence to Núria Alonso.

Ethics declarations

Ethics approval and consent to participate

The original study was approved by the regional ethics committee (H-3-2009-139), and the Danish Data Protection Agency (00934-Geh-2010-003). The investigation has been conducted according to the principles expressed in the Declaration of Helsinki, and all the included participants provided written informed consents.

Consent for publication

Not applicable.

Competing interests

NAM is stockholders of and serve on the board of director of Biosfer Teslab, a diagnostic laboratory company that performed the lipoprotein and glycoprotein analyses described herein. The remaining authors declared not having anything to disclose regarding conflict of interest with respect to this manuscript.

Additional information

Publisher’s Note

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

Dídac Mauricio, Magnus T Jensen and Núria Alonso share senior authorship

Electronic supplementary material

Below is the link to the electronic supplementary material.

12933_2022_1652_MOESM1_ESM.docx

Supplementary Material 1: Supplementary Fig. 1: Comparison of the distribution of the predicted probabilities showing the classification performance of the presence of MCD between two different models, being the x-axis the predicted probabilities for both classes and the y-axis the count of observations. A: Model 1 includes classical risk variables (age, sex, eGFR, NTproBNP, BMI, diabetes duration and systolic blood pressure >140mmHg). B: Model 2 includes classical risk variables and the NMR-assessed biomarkers. The inclusion of NMR parameters significantly increased the AUROC from 0.62 [0.56–0.68] to 0.67 [0.61–0.73], with a NRI considering NMR-assessed parameters of 21%. MCD: myocardial dysfunction, BMI: body mass index, NMR: nuclear magnetic resonance, AUROC: area under the ROC curve, NRI: net reclassification improvement.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Puig-Jové, C., Julve, J., Castelblanco, E. et al. The novel inflammatory biomarker GlycA and triglyceride-rich lipoproteins are associated with the presence of subclinical myocardial dysfunction in subjects with type 1 diabetes mellitus. Cardiovasc Diabetol 21, 257 (2022). https://doi.org/10.1186/s12933-022-01652-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12933-022-01652-z

Keywords