Skip to main content
  • More from ADA
    • Diabetes Care
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care in Diabetes
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care
  • Subscribe
  • Log in
  • My Cart
  • Follow ada on Twitter
  • RSS
  • Visit ada on Facebook
Diabetes

Advanced Search

Main menu

  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
  • More from ADA
    • Diabetes Care
    • Clinical Diabetes
    • Diabetes Spectrum
    • ADA Standards of Medical Care in Diabetes
    • ADA Scientific Sessions Abstracts
    • BMJ Open Diabetes Research & Care

User menu

  • Subscribe
  • Log in
  • My Cart

Search

  • Advanced search
Diabetes
  • Home
  • Current
    • Current Issue
    • Online Ahead of Print
    • ADA Scientific Sessions Abstracts
  • Browse
    • By Topic
    • Issue Archive
    • Saved Searches
    • ADA Scientific Sessions Abstracts
    • Diabetes COVID-19 Article Collection
    • Diabetes Symposium 2020
  • Info
    • About the Journal
    • About the Editors
    • ADA Journal Policies
    • Instructions for Authors
    • Guidance for Reviewers
  • Reprints/Reuse
  • Advertising
  • Subscriptions
    • Individual Subscriptions
    • Institutional Subscriptions and Site Licenses
    • Access Institutional Usage Reports
    • Purchase Single Issues
  • Alerts
    • E­mail Alerts
    • RSS Feeds
  • Podcasts
    • Diabetes Core Update
    • Special Podcast Series: Therapeutic Inertia
    • Special Podcast Series: Influenza Podcasts
    • Special Podcast Series: SGLT2 Inhibitors
    • Special Podcast Series: COVID-19
  • Submit
    • Submit a Manuscript
    • Submit Cover Art
    • ADA Journal Policies
    • Instructions for Authors
    • ADA Peer Review
Genetics/Genomes/Proteomics/Metabolomics

Circulating Metabolites and Lipids Are Associated to Diabetic Retinopathy in Individuals With Type 1 Diabetes

  1. Viktor Rotbain Curovic1⇑,
  2. Tommi Suvitaival1,
  3. Ismo Mattila1,
  4. Linda Ahonen1,
  5. Kajetan Trošt1,
  6. Simone Theilade1,
  7. Tine W. Hansen1,
  8. Cristina Legido-Quigley1,2 and
  9. Peter Rossing1,3
  1. 1Steno Diabetes Center Copenhagen, Gentofte, Denmark
  2. 2Institute of Pharmaceutical Science, Faculty of Life Sciences & Medicine, King’s College London, London, U.K.
  3. 3University of Copenhagen, Copenhagen, Denmark
  1. Corresponding author: Viktor Rotbain Curovic, viktor.rotbain.curovic{at}regionh.dk
  1. V.R.C., T.S., C.L.-Q., and P.R. contributed equally.

Diabetes 2020 Oct; 69(10): 2217-2226. https://doi.org/10.2337/db20-0104
PreviousNext
  • Article
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF
Loading

Abstract

Omics-based methods may provide new markers associated to diabetic retinopathy (DR). We investigated a wide omics panel of metabolites and lipids related to DR in type 1 diabetes. Metabolomic analyses were performed using two-dimensional gas chromatography with time-of-flight mass spectrometry and lipidomic analyses using an ultra-high-performance liquid chromatography quadruple time-of-flight mass spectrometry method in 648 individuals with type 1 diabetes. Subjects were subdivided into no DR, mild nonproliferative DR (NPDR), moderate NPDR, proliferative DR, and proliferative DR with fibrosis. End points were any progression of DR, onset of DR, and progression from mild to severe DR tracked from standard ambulatory care and investigated using Cox models. The cohort consisted of 648 participants aged a mean of 54.4 ± 12.8 years, 55.5% were men, and follow-up was 5.1–5.5 years. Cross-sectionally, 2,4-dihydroxybutyric acid (DHBA), 3,4-DHBA, ribonic acid, ribitol, and the triglycerides 50:1 and 50:2 significantly correlated (P < 0.042) to DR stage. Longitudinally, higher 3,4-DHBA was a risk marker for progression of DR (n = 133) after adjustment (P = 0.033). We demonstrated multiple metabolites being positively correlated to a higher grade of DR in type 1 diabetes and several triglycerides being negatively correlated. Furthermore, higher 3,4-DHBA was an independent risk marker for progression of DR; however, confirmation is required.

Introduction

One of the most frequent and debilitating complications of diabetes is diabetic retinopathy (DR). Despite radical improvements in diagnosis and treatment throughout the last decades (1), DR is still the primary cause of blindness in individuals with diabetes aged 20–74 years, and the prevalence was 34.6% in a large meta-analysis including 22,896 individuals with diabetes (2). The advances in diagnosis and treatment have particularly been made for later stages of the disease, and robust and specific risk markers for onset and early progression of DR are still lacking.

Novel omics methods have been developed, allowing for simultaneous evaluation of large panels of metabolites using mass spectrometry–based approaches that allow for comprehensive study of metabolic pathways compared with more traditional single-biomarker approaches. Omics facilitates advanced and detailed analysis faster than standard methods. More specifically, metabolomics and lipidomics is the analysis and categorization of circulating metabolites and lipids using this method and can provide a more comprehensive view of the biological effects of various metabolites (3,4).

Applying this method to DR is an intriguing concept. Omics allows for a unique ability to understand the biological pathways for DR as well as facilitate the discovery of novel biomarkers associated to DR development and progression. Few studies have evaluated the association between circulating metabolites and the presence of DR. In one cross-sectional study, several metabolites were associated with DR in type 2 diabetes, among them, hydroxyl fatty acids, such as 3,4 dihydroxybutyric acid (3,4-DHBA), and sugar derivatives such as lactose, maltose, and ribose (5). Likewise, other studies have identified plasma metabolites associated to arginine-, pyrimidine- and fatty acid-related pathways, as well as several amino acids highly related to insulin resistance, in relation to DR (6–8).

In the current study, we investigated the predictive qualities of a wide panel of metabolites and lipids in plasma in relation to the presence, onset, and progression of DR in individuals with type 1 diabetes (T1D).

Research Design and Methods

Study Population

Between 2009 and 2011, 648 individuals with T1D and a large range of albuminuria were recruited from the outpatient clinic at Steno Diabetes Center Copenhagen. The details of the cohort have previously been described (9). Participants were subdivided by stages of albuminuria (normo-, micro-, and macroalbuminuria). End-stage kidney disease, defined as receiving dialysis, renal transplantation, or an estimated glomerular filtration rate (eGFR) <15 mL/min/1.73 m2 at baseline was an exclusion criterion. In the current study, metabolomics and lipidomics data along with information on retinopathy status was available for 601 (serum metabolomics) and 648 (plasma lipidomics) participants, respectively.

The study was conducted in compliance with the Declaration of Helsinki and was approved by the ethics committee for the Capital Region of Denmark (Hillerød, Denmark). All participants have given informed written consent.

Baseline Clinical Analyses

Serum creatinine, plasma LDL cholesterol, triglycerides, and HbA1c were measured using standardized methods from venous samples. The urinary albumin excretion rate (UAER) was analyzed by enzyme immunoassay based on three consecutive 24-h urine collections. eGFR was calculated based on serum creatinine using the Chronic Kidney Disease Epidemiology Collaboration equation. An automated validated device was used to measure sitting brachial blood pressure after a 10-min rest.

Baseline DR stage was classified using in-house algorithms on a 0–4 scale based on regular retinopathy screenings at Steno Diabetes Center Copenhagen performed by specifically trained and certified nursing staff under the supervision of ophthalmologists. Five mydriatic nonstereoscopic fundus photos are taken; one macula-centered and four peripheral 45° fundus images. These are combined into a mosaic, grading the macula and periphery separately, according to a modified version of the international classification of DR disease severity scale (10). The presence of microaneurysms, intraretinal hemorrhages, hard or soft exudates, and proliferations are recorded and quantified. Likewise, intraretinal microvascular abnormalities and venous beading are recorded. A weighted quantitation of the distinct retinal pathologies was used to perform staging. Overall stage was defined as the highest stage diagnosed in either eye. Stage 0 is defined as no DR in any eye, stage 1 as mild nonproliferative DR (NPDR), stage 2 as moderate NPDR, stage 3 as proliferative DR (PDR), and stage 4 as PDR with fibrosis. Blind subjects are not screened for retinopathy at Steno Diabetes Center Copenhagen and were therefore excluded. Blindness was defined as visual acuity of less than 1/60, lack of ability to count fingers in front of a white screen at 1 m distance, and lack of ability to see hand motion in front of a white screen at 1 m distance.

Sample Quantification and Identification

The metabolomics analysis is detailed in Tofte et al (11), and the lipidomics analysis is detailed in Tofte et al. (12). For completeness, the analyses are outlined here as follows: Serum samples, stored at −80°C, were analyzed by two different analytical methods. Metabolomics samples were analyzed using a two-dimensional gas chromatography with time-of-flight mass spectrometry. Peak-picking from the raw data was performed with ChromaTOF, and the resulting features were aligned with Guineu (13).

Lipidomics samples were prepared using a modified Folch extraction procedure (14) and analyzed by a previously presented ultra-high-performance liquid chromatography quadrupole time-of-flight mass spectrometry method (15). The raw data were preprocessed with MZmine 2 (16). A complete list of identified metabolites is available in Tofte et al. (11). Finally, the metabolomics and lipidomics data were postprocessed in R software, as described previously (11,12). Lipid species are defined as the number of carbon atoms (indicating total fatty acid chain length) and number double-bonds for the specific species. They are presented as “species (number of carbon atoms:number of double-bonds).”

Within the coverage of the two mass spectrometry platforms, the inclusion of metabolites and lipids in subsequent data analysis was solely based on the certainty of identification and the level of technical precision, thereby not restricting to any particular pathway or prior hypothesis.

Follow-up

Data regarding retinopathy were obtained using local electronic records from Steno Diabetes Center Copenhagen up to 31 December 2016 and were available for 563 subjects. The end points were defined as 1) progression from any stage to any other stage of DR (any progression); 2) onset of DR; and 3) progression from stage 1–2 to stage 3–4 (progression from mild to severe DR). In the case of participants experiencing multiple end points, only the first occurrence was included.

Statistical Analysis

Continuous variables are presented as mean (SD), if normally distributed, and as median (interquartile range [IQR]) if skewed. Before all analyses, skewed variables were log2 transformed, including all metabolites, to achieve normal distribution. Categorical variables are presented as total number (%). Baseline clinical characteristics were compared across baseline DR status using ANOVA and the χ2 test for continuous and categorical variables, respectively.

Metabolites and lipid species were analyzed using a narrowing-down approach in relation to DR stages and outcomes as follows: Cross-sectional relationships between single metabolites or lipid species and baseline DR stages were assessed using multivariate linear regression models adjusted for relevant clinical variables. Thereafter, the single measures were cross-sectionally associated to categories of DR stage and tested using ANCOVA. The Benjamini-Hochberg (BH) method (PBH) (17,18) was used to correct for multiple testing for presented P values throughout the analysis. Metabolites with PBH < 0.05 and lipids with PBH < 0.1 in the adjusted cross-sectional model were included in survival analysis with the Cox proportional hazards model for end points. All hazard ratios (HRs) are reported per doubling of the metabolite or lipid.

Clinical variables in the adjusted models were age, sex, HbA1c, systolic blood pressure, smoking, BMI, statin treatment, triglycerides, LDL cholesterol, and prescribed antihypertensive medication.

Partial correlation network analysis was done with the R package “huge” (19) using the graphical lasso algorithm (20) coupled with the extended Bayesian information criterion (21). All metabolites were included in inferring the network. Subsequently, the subnetwork of metabolites, which were immediately connected to the four retinopathy-associated metabolites, was visualized with the R package qgraph (22). Edges or the network were colored by the respective partial correlation in the graphical lasso model, and nodes were colored with respective Spearman correlation to the top candidate biomarkers. In both color annotations, red and blue refer to positive and negative correlation, respectively. Size of nodes refers to the respective degree (i.e., the number of associations to other nodes). Statistical analysis and data visualization were performed using R (version 3.4.2).

Data and Resource Availability

The data sets generated during and/or analyzed during the current study are available from the corresponding authors on reasonable request.

Results

Baseline Characteristics

Baseline characteristics for the participants divided according to baseline DR stages are summarized in Table 1. Most had no DR (n = 141 [22%]) or moderate NPDR (n = 186 [29%]). Compared with no DR, participants with PDR with fibrosis had higher mean (SD) BMI (24.7 [3.5] vs. 26.8 [11.2] kg/m2), UAER (11.5 [IQR 7.6–22.4] vs. 49.0 [IQR 14.9–231.8] mg/24 h), and systolic blood pressure (127 [(46.1] vs. 135 [20.3] mmHg), and lower eGFR (94.8 [24.7] vs. 64.5 [30.6] mL/min/1.73 m2). Across groups, age, systolic blood pressure, and UAER were higher with higher DR stage, and the eGFR was lower. Also, the frequency of treatment with antihypertensive drugs and statins was higher with higher DR stage.

View this table:
  • View inline
  • View popup
Table 1

Baseline characteristics across DR stage

Metabolomic Cross-sectional Analyses

A total of 75 metabolite species were identified and passed quality control (Supplementary Table 1). These metabolites were included in the multivariate linear regression models. Four metabolites were positively correlated to baseline DR stage after adjustment for clinical variables and correction for multiple testing: 2,4-DHBA (PBH < 0.001), ribonic acid (PBH = 0.017), ribitol (PBH = 0.032), and 3,4-DHBA (PBH = 0.036). Thereafter, all 75 metabolites were included in ANCOVA for baseline DR stages. The same four metabolite levels were significantly increased by higher DR stage after adjustment (ribonic acid, PBH < 0.001; 2,4-DHBA, PBH < 0.001; ribitol, PBH = 0.013; and 3,4-DHBA, PBH = 0.041) (Table 2). Figure 1 illustrates the distribution and relative levels of these four metabolites. It is apparent that the levels of all four metabolites increase with higher DR stage. Both 2,4-DHBA and 3,4-DHBA (Fig. 1B and D) follow the same pattern of lower levels of the metabolite for no DR and NPDR and a nonlinear increase toward higher DR stages and PDR with fibrosis, in particular. In contrast, the levels of ribonic acid and ribitol (Fig. 1A and C) increase linearly with increasing DR stage.

View this table:
  • View inline
  • View popup
Table 2

Cross-sectional association between metabolites and lipids and baseline DR grade

Figure 1
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1

Violin plots show metabolite levels across baseline DR stages for four metabolites—ribonic acid (A), 2,4-DHBA (B), ribitol (C), and 3,4-DHBA (D)—where a difference between the stages was detected using ANCOVA. Observations are shown as dots, and their distribution in each of the stages of retinopathy as a violin geom. Pairwise differences between the stages are indicated with P values at the top of the figure. The results are from models with adjustment for age, sex, HbA1c, systolic blood pressure, smoking, BMI, statin treatment, triglycerides, LDL cholesterol, and prescribed antihypertensive medication.

The partial correlation network of the metabolome connected to the four metabolites is shown in Fig. 2. This network included kidney function–related metabolites (e.g., creatinine and myo-inositol), glucose metabolism–related compounds (e.g., citric acid and glycine), fatty acids (e.g., fumaric acid and malic acid), and amino acids (e.g., alanine and serine). Particularly, creatinine and myo-inositol, but also glyceryl-glycoside, 4-hydroxybenzeneacetic acid, and fumaric acid, were highly associated with the hub of the four retinopathy-associated metabolites. Except for the amino acids, most of the compounds in the network were positively correlated with the four highlighted metabolites, as indicated by the red color in Fig. 2.

Figure 2
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2

Partial correlation network shows associations in the measured metabolome. Metabolites with cross-sectional association to retinopathy (ribitol, ribonic acid, 2,4-DHBA and 3,4-DHBA) are indicated as diamond-shaped nodes. The node color of the other metabolites (circular nodes) indicates the Spearman correlation to ribonic acid. Partial correlations (i.e., independent associations) are indicated by lines between the metabolites, where thickness and color of the line, respectively, indicate the strength and the sign of the association (red: positive, blue: inverse). Additionally, the size of the node indicates the degree of the node (i.e., the number of associations with the metabolite).

Lipidomic Cross-sectional Analyses

After identification and quality control, 104 lipid species from the following five major lipid classes were included in the analyses: diacyl-phosphatidylcholines (PCs), alkyl-acyl-phosphatidylcholines, lyso-phosphatidylcholines (LPCs), triacylglycerols (TGs), and sphingomyelins. Lipid species are defined as number of carbon atoms (indicating total fatty acid chain length) and number double-bonds for the specific species. They are presented as “species (number of carbon atoms:number of double-bonds).” The investigated lipids are listed in Supplementary Table 2.

TG(50:1) and TG(50:2) were inversely associated with DR grade at baseline after adjustment for clinical covariates, when testing with linear regression analysis (PBH < 0.05). Furthermore, LPC(16:1), PC(32:1), TG(14:0/16:0/18:1), TG(50:3), and PC(32:2) were inversely associated with baseline DR grade at a higher false discovery rate of 10%.

When investigating with ANOVA, a difference in the lipid level between the DR grades was detected in LPC(16:0) (PBH < 0.1). This association was lost after adjustment for the clinical covariates (ANCOVA). In addition to LPC(16:0), medium-sized unsaturated TGs and small LPCs had stronger indicative associations with DR grade than other lipids, as shown in the lipidome-wide heatmap of the ANCOVA F statistics (Fig. 3). In particular, LPCs(16:1), as well as TG(49:3), TG(50:1), and TG(50:2), emerged with an indicative association with the DR grade.

Figure 3
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3

Heatmaps of the ANCOVA model F statistics across the entire lipidomic panel in relation to DR stage. The results are from models with adjustment for age, sex, HbA1c, systolic blood pressure, smoking, BMI, statin treatment, triglycerides, LDL cholesterol, and prescribed antihypertensive medication. Lipid species are grouped according to the lipid classes in each panel. Each cell represents one lipid species. On the y-axis is number of double-bonds for the specific species (indicating level of saturation), and on the x-axis is the number of carbon atoms (indicating total fatty acid chain length). PC-O/P, alkyl-acyl-phosphatidylcholines; SM, sphingomyelins.

Longitudinal Analyses

Metabolites and lipids identified in cross-sectional analyses were thereafter analyzed with compound-specific Cox proportional hazards models for association to any progression, onset of DR, and progression from mild to severe DR. Median follow-up ranged between 5.1 and 5.5 years depending on the end point. The number of events was 133, 47, and 29 for any progression, onset of DR, and progression from mild to severe DR, respectively. For the any progression end point, higher 3,4-DHBA exhibited significance after adjustment for clinical covariates and multiple testing (HR 1.55, 95% CI 1.12–2.15, P = 0.033). The other metabolites were not associated with any of the end points neither before nor after adjustment. Although not statistically significant, 2,4-DHBA showed a high HR for progression from mild to severe DR (HR 1.92, 95% CI 0.94–3.93, P = 0.290). Unadjusted and adjusted HRs for the metabolites are presented as a forest plot in Fig. 4.

Figure 4
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4

Forest plot of multiple testing–corrected HRs of the selected metabolites for three retinopathy events: any retinopathy, onset of retinopathy, and progression from mild to severe retinopathy (top, middle, and bottom panels, respectively). Adjusted model (right) is adjusted for age, sex, HbA1c, systolic blood pressure, smoking, BMI, statin treatment, triglycerides, LDL cholesterol, and prescribed antihypertensive medication. The crude model (left) is without these adjustments.

Unlike the metabolites, no lipids were independent risk factors for any of the end points (P > 0.05).

Discussion

The current study illustrates an exciting new avenue in characterizing individuals with T1D with DR. In the present cohort we have investigated individuals with long diabetes duration and a broad range of albuminuria, leading to a high proportion of subjects with more severe DR than would be expected in a general clinical population with T1D. We identified four metabolites associated with presence of DR as well as higher 3,4-DHBA as an independent risk marker for progression of DR. Our results were independent of a panel of metabolic risk factors traditionally used for risk stratification of DR in T1D. Therefore, we now argue for the need for further investigation of omics-based risk stratification of DR. Using omics in relation to DR is a relatively new venture, and clinical studies assessing its viability are sparse, especially in subpopulations such as individuals with T1D, and using longitudinal data.

The metabolites identified in this study mainly stem from two etiopathogenic factors, namely, hyperglycemia (ribitol and ribonic acid) and dyslipidemia (2,4- and 3,4-DHBA) (23). Ribitol and ribonic acid are derivatives from ribose, which is highly active in the pentose phosphate pathway and in the production of nucleotides and nucleic acids. Furthermore, sugar alcohols, such as sorbitol, are active in the polyol pathway, which has been identified as a crucial insulin-independent pathway relevant in the onset of DR, and have been suggested as a possible therapeutic target in the treatment of DR (24). The fructose created in this pathway becomes further phosphorylated, resulting in the formation of advanced glycation end products (AGEs), which in turn bind to receptors for AGEs (RAGE)—a known facilitator of DR (25).

DHBAs, on the other hand, which are closely related to the ketone body hydroxybutyric acid, have not been directly associated with any major pathways implicated in the onset of DR or T1D. Other diseases, such as succinic semialdehyde dehydrogenase deficiency, an autosomal recessive genetic disease, leads to 4-DHBA aggregation, and, in turn, is associated to severe neurological complications and symptoms (26). In general, ketone bodies are associated with dyslipidemia and high-fat diets (27), which are highly relevant risk factors in the development of diabetic complications. A theory that has drawn recent attention is that DHBAs could be implicated in the butyrate metabolism by the gut microbiota (28), although at present this still calls for more investigation. Similarly, Sumarriva et al. (6) demonstrated that higher plasma carnitine, a metabolite highly present in food that contains meat, was associated to the presence of PDR compared with NPDR. Carnitine is further metabolized by the gut microbiota into trimethylamine-N-oxide, which has been associated to cardiovascular and metabolic diseases (29).

Interestingly, visualized in Fig. 2, the subnetwork of metabolites associated with DR shows that all metabolites that correlated significantly with DR in this study all seem strongly associated to myo-inositol. Despite myo-inositol itself not being associated to DR, in this study, the association to the other metabolites could propose another pathway of DR etiology. myo-Inositol is a sugar alcohol, the gastrointestinal absorption and intracellular transport of which has been shown to be impaired in individuals with diabetes (during hyperglycemia) (33). Furthermore, while previously thought mainly to be expressed in renal tissue, recent studies show evidence of myo-inositol activity in extrarenal tissue, such as retinal and lens epithelium as well (34,35), and has additionally been associated to existing DR in individuals with type 2 diabetes (36).

Our metabolomic results stand well in relation to a study by Chen et al. (5), describing similar metabolites associated to DR in a cross-sectional study including individuals with type 2 diabetes with DR (n = 40) or without DR (n = 40) in the Singapore Indian Eye Study. They demonstrated that 3,4-DHBA as well as ribose were significantly higher in the DR group compared with the group without DR. Furthermore, similar results were found for 2-deoxyribonic acid as well as for several other sugars and sugar derivates (5). As such, we can partly validate these results in our larger population of individuals with T1D, and additionally, we have shown that 3,4-DHBA was a risk marker for progression of DR during follow-up. The inherent issue with omics discovery studies being explorative is that replication of results is necessary across populations and cohorts, but is often difficult, especially due to platform heterogeneity. The study by Chen et al. (5) strengthens the findings in our study suggesting that 3,4-DHBA and ribose derivatives are valid risk markers of DR. However, Lin et al. (8) showed that branched-chain, aromatic, and glucogenic amino acids, such as leucine, valine, tyrosine, and alanine, were positively associated with diabetic microangiopathy in type 2 diabetes. These amino acids were also investigated in our study; however, we could not confirm the results, possibly due to the heterogeneity between type 1 and type 2 diabetes.

In the case of lipids, some, such as LDL cholesterol and TGs, have throughout many years been comprehensively studied and have proved to be robust risk markers for vascular disease in T1D. In addition, studies targeting dyslipidemia with fibrates have found beneficial effects on retinopathy (37,38). Therefore, the concept of applying lipidomic strategies into finding novel markers to strengthen the identification and prediction of vascular risk is not implausible. In our panel of lipids, we were only able to identify three TGs negatively correlated with DR stage in our linear models using a 5% α-level, adjusted for, among others, baseline LDL cholesterol, triglycerides, and BMI. However, these associations could not be replicated in ANCOVA models, arguably in part due to nonlinear trends for different lipids across DR stages.

Comparing the present results in relation to other studies is difficult because very few have investigated DR. One study found stearic acid, trans oleic acid, linoleic acid, arachidonic acid, and free cholesterol in circulation to be significant differentiators between preclinical DR, NPDR, and PDR (39). Likewise, Schwartzman et al. (40) identified several free-fatty acid autacoids in vitreous humor as markers of PDR in T1D.

PDR is demonstrated to be associated to an impaired blood-brain barrier (BBB), and early damage to the BBB is hypothesized to be a predictor of progression to more advanced stages of DR (41). However, the mechanisms for the association surrounding BBB impairment and DR are not understood, and the association to circulating biomarkers has not been described. Hogan et al. (42) demonstrated that various polyunsaturated fatty acids and sphingolipids were associated with impaired BBB in traumatic brain injury in rats, but our results do not support the association between these and DR.

Moving outside of purely metabolomic and lipidomic studies, a substantial amount of research has been performed on the proteome and its effect on DR risk and risk progression, however, primarily in small samples with largely nonreplicated results (43).

This study is not without limitations. The DR staging during follow-up did not take laser surgery or vascular endothelial growth factor injections into account, because it was based only on changes in DR stage from baseline. Furthermore, no data on concomitant medication throughout the follow-up were available, and as such, there are no data on how treatment with statins, antihypertensive medication, or insulin has changed during follow-up. In addition, no information on lifestyle parameters, which could have influence on lipid composition, was available at baseline. Finally, the lack of a validation population is another limiting factor. Nonetheless, the sizable strengths of this study are, firstly, the large, well-defined cohort of individuals with T1D, including 7 years of longitudinal data, and secondly, a comprehensive metabolomic and lipidomic analysis regarding presence of and changes in DR in T1D.

In summary, we identified four metabolites and three lipids with an association to the DR stage: ribonic acid, ribitol, and two DHBAs were associated with DR stage, and three triglycerides were negatively correlated with the DR stage. Furthermore, we have identified 3,4-DHBA as an independent risk marker for progression in DR stage. Our results may serve as a basis for further studies regarding sugar metabolism, hydroxy acids, and lipids in relation to diabetic complications such as retinopathy, because more investigative studies are needed before these markers can be clinically applied.

Article Information

Acknowledgments. The authors wish to acknowledge the excellent technical assistance of Steno Diabetes Center Copenhagen staff Tina R. Juhl, Anne G. Lundgaard, Berit R. Jensen, Jessie Hermann, and Ulla M. Smidt.

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

Author Contributions. V.R.C. wrote the manuscript and is responsible for the integrity of the work as a whole. V.R.C., T.S., I.M., L.A., S.T., T.W.H., C.L.-Q., and P.R. conceived and designed the research. V.R.C., T.S., I.M., L.A., S.T., T.W.H., C.L.-Q., and P.R. analyzed and interpreted the data. V.R.C., T.S., I.M., L.A., S.T., T.W.H., C.L.-Q., and P.R. critically revised the manuscript for key intellectual content. T.S. performed the statistical analysis. P.R. obtained funding and supervised the study. All authors approved the final version of the manuscript. V.R.C. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Prior Presentation. Parts of this study were presented as two abstracts at the 56th Eastern Association for the Study of Diabetes Annual Meeting, which was held online, 21–25 September 2020.

Footnotes

  • This article contains supplementary material online at https://doi.org/10.2337/figshare.12733202.

  • L.A. is currently affiliated with Biosyntia ApS, Copenhagen, Denmark.

  • Received January 27, 2020.
  • Accepted July 29, 2020.
  • © 2020 by the American Diabetes Association
https://www.diabetesjournals.org/content/license

Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. More information is available at https://www.diabetesjournals.org/content/license.

References

  1. ↵
    1. Klein R,
    2. Lee KE,
    3. Gangnon RE,
    4. Klein BE
    . The 25-year incidence of visual impairment in type 1 diabetes mellitus the Wisconsin epidemiologic study of diabetic retinopathy. Ophthalmology 2010;117:63–70
    OpenUrlCrossRefPubMedWeb of Science
  2. ↵
    1. Yau JWY,
    2. Rogers SL,
    3. Kawasaki R, et al.; Meta-Analysis for Eye Disease (META-EYE) Study Group
    . Global prevalence and major risk factors of diabetic retinopathy. Diabetes Care 2012;35:556–564
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Wenk MR
    . The emerging field of lipidomics. Nat Rev Drug Discov 2005;4:594–610
    OpenUrlCrossRefPubMedWeb of Science
  4. ↵
    1. Hyötyläinen T,
    2. Orešič M
    . Analytical lipidomics in metabolic and clinical research. Trends Endocrinol Metab 2015;26:671–673
    OpenUrl
  5. ↵
    1. Chen L,
    2. Cheng CY,
    3. Choi H, et al
    . Plasma metabonomic profiling of diabetic retinopathy. Diabetes 2016;65:1099–1108
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Sumarriva K,
    2. Uppal K,
    3. Ma C, et al
    . Arginine and carnitine metabolites are altered in diabetic retinopathy. Invest Ophthalmol Vis Sci 2019;60:3119–3126
    OpenUrl
    1. Zhu X-R,
    2. Yang F-Y,
    3. Lu J, et al
    . Plasma metabolomic profiling of proliferative diabetic retinopathy. Nutr Metab (Lond) 2019;16:37
    OpenUrl
  7. ↵
    1. Lin HT,
    2. Cheng ML,
    3. Lo CJ, et al
    . 1H nuclear magnetic resonance (NMR)-based cerebrospinal fluid and plasma metabolomic analysis in type 2 diabetic patients and risk prediction for diabetic microangiopathy. J Clin Med 2019;8:874
    OpenUrl
  8. ↵
    1. Theilade S,
    2. Lajer M,
    3. Hansen TW,
    4. Rossing P
    . Pulse wave reflection is associated with diabetes duration, albuminuria and cardiovascular disease in type 1 diabetes. Acta Diabetol 2014;51:973–980
    OpenUrlPubMed
  9. ↵
    1. Wilkinson CP,
    2. Ferris FL III.,
    3. Klein RE, et al.; Global Diabetic Retinopathy Project Group
    . Proposed international clinical diabetic retinopathy and diabetic macular edema disease severity scales. Ophthalmology 2003;110:1677–1682
    OpenUrlCrossRefPubMedWeb of Science
  10. ↵
    1. Tofte N,
    2. Suvitaival T,
    3. Trost K, et al
    . Metabolomic assessment reveals alteration in polyols and branched chain amino acids associated with present and future renal impairment in a discovery cohort of 637 persons with type 1 diabetes. Front Endocrinol (Lausanne) 2019;10:818
    OpenUrl
  11. ↵
    1. Tofte N,
    2. Suvitaival T,
    3. Ahonen L, et al
    . Lipidomic analysis reveals sphingomyelin and phosphatidylcholine species associated with renal impairment and all-cause mortality in type 1 diabetes. Sci Rep 2019;9:16398
    OpenUrlCrossRef
  12. ↵
    1. Castillo S,
    2. Mattila I,
    3. Miettinen J,
    4. Orešič M,
    5. Hyötyläinen T
    . Data analysis tool for comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry. Anal Chem 2011;83:3058–3067
    OpenUrlCrossRefPubMed
  13. ↵
    1. Folch J,
    2. Lees M,
    3. Sloane Stanley GH
    . A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem 1957;226:497–509
    OpenUrlFREE Full Text
  14. ↵
    1. O’Gorman A,
    2. Suvitaival T,
    3. Ahonen L, et al
    . Identification of a plasma signature of psychotic disorder in children and adolescents from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. Transl Psychiatry 2017;7:e1240
    OpenUrl
  15. ↵
    1. Pluskal T,
    2. Castillo S,
    3. Villar-Briones A,
    4. Oresic M
    . MZmine 2: modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics 2010;11:395
    OpenUrlCrossRefPubMed
  16. ↵
    1. Darshi M,
    2. Van Espen B,
    3. Sharma K
    . Metabolomics in diabetic kidney disease: unraveling the biochemistry of a silent killer. Am J Nephrol 2016;44:92–103
    OpenUrlPubMed
  17. ↵
    1. Storey JD,
    2. Tibshirani R
    . Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 2003;100:9440–9445
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Zhao T,
    2. Liu H,
    3. Roeder K,
    4. Lafferty J,
    5. Wasserman L
    . The huge package for high-dimensional undirected graph estimation in R. J Mach Learn Res 2012;13:1059–1062
    OpenUrlPubMed
  19. ↵
    1. Friedman J,
    2. Hastie T,
    3. Tibshirani R
    . Sparse inverse covariance estimation with the graphical lasso. Biostatistics 2008;9:432–441
    OpenUrlCrossRefPubMedWeb of Science
  20. ↵
    1. Chen J,
    2. Chen Z
    . Extended Bayesian information criteria for model selection with large model spaces. Biometrika 2008;95:759–771
    OpenUrlCrossRefWeb of Science
  21. ↵
    1. Epskamp S,
    2. Cramer AOJ,
    3. Waldorp LJ,
    4. Schmittmann VD,
    5. Borsboom D
    . qgraph: network visualizations of relationships in psychometric data. J Stat Softw 2012;48:18
    OpenUrl
  22. ↵
    1. Ahsan H
    . Diabetic retinopathy--biomolecules and multiple pathophysiology. Diabetes Metab Syndr 2015;9:51–54
    OpenUrl
  23. ↵
    1. Kinoshita JH
    . A thirty year journey in the polyol pathway. Exp Eye Res 1990;50:567–573
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    1. Mahajan N,
    2. Arora P,
    3. Sandhir R
    . Perturbed biochemical pathways and associated oxidative stress lead to vascular dysfunctions in diabetic retinopathy. Oxid Med Cell Longev 2019;2019:8458472
    OpenUrl
  25. ↵
    1. Pearl PL,
    2. Novotny EJ,
    3. Acosta MT,
    4. Jakobs C,
    5. Gibson KM
    . Succinic semialdehyde dehydrogenase deficiency in children and adults. Ann Neurol 2003;54(Suppl. 6):S73–S80
    OpenUrlPubMed
  26. ↵
    1. Sikder K,
    2. Shukla SK,
    3. Patel N,
    4. Singh H,
    5. Rafiq K
    . High fat diet upregulates fatty acid oxidation and ketogenesis via intervention of PPAR-γ. Cell Physiol Biochem 2018;48:1317–1331
    OpenUrlCrossRef
  27. ↵
    1. Jain A,
    2. Li XH,
    3. Chen WN
    . An untargeted fecal and urine metabolomics analysis of the interplay between the gut microbiome, diet and human metabolism in Indian and Chinese adults. Sci Rep 2019;9:9191
    OpenUrl
  28. ↵
    1. Papandreou C,
    2. Moré M,
    3. Bellamine A
    . Trimethylamine N-oxide in relation to cardiometabolic health-cause or effect? Nutrients 2020;12:1330
    OpenUrl
    1. Boden G
    . Role of fatty acids in the pathogenesis of insulin resistance and NIDDM. Diabetes 1997;46:3–10
    OpenUrlAbstract/FREE Full Text
    1. McGarry JD,
    2. Dobbins RL
    . Fatty acids, lipotoxicity and insulin secretion. Diabetologia 1999;42:128–138
    OpenUrlCrossRefPubMedWeb of Science
    1. Yore MM,
    2. Syed I,
    3. Moraes-Vieira PM, et al
    . Discovery of a class of endogenous mammalian lipids with anti-diabetic and anti-inflammatory effects. Cell 2014;159:318–332
    OpenUrlCrossRefPubMed
  29. ↵
    1. Clements RS Jr..,
    2. Reynertson R
    . Myoinositol metabolism in diabetes mellitus: effect of insulin treatment. Diabetes 1977;26:215–221
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Arner RJ,
    2. Prabhu KS,
    3. Krishnan V,
    4. Johnson MC,
    5. Reddy CC
    . Expression of myo-inositol oxygenase in tissues susceptible to diabetic complications. Biochem Biophys Res Commun 2006;339:816–820
    OpenUrlCrossRefPubMedWeb of Science
  31. ↵
    1. Li WY,
    2. Zhou Q,
    3. Qin M,
    4. Tao L,
    5. Lou M,
    6. Hu TS
    . Reduced absolute rate of myo-inositol biosynthesis of cultured bovine retinal capillary pericytes in high glucose. Exp Eye Res 1991;52:569–573
    OpenUrlCrossRefPubMed
  32. ↵
    1. Tong J,
    2. Geng H,
    3. Zhang Z, et al
    . Brain metabolite alterations demonstrated by proton magnetic resonance spectroscopy in diabetic patients with retinopathy. Magn Reson Imaging 2014;32:1037–1042
    OpenUrlCrossRefPubMed
  33. ↵
    1. Keech AC,
    2. Mitchell P,
    3. Summanen PA, et al.; FIELD study investigators
    . Effect of fenofibrate on the need for laser treatment for diabetic retinopathy (FIELD study): a randomised controlled trial. Lancet 2007;370:1687–1697
    OpenUrlCrossRefPubMedWeb of Science
  34. ↵
    1. Chew EY,
    2. Ambrosius WT,
    3. Davis MD, et al.; ACCORD Study Group; ACCORD Eye Study Group
    . Effects of medical therapies on retinopathy progression in type 2 diabetes. N Engl J Med 2010;363:233–244
    OpenUrlCrossRefPubMedWeb of Science
  35. ↵
    1. Li X,
    2. Luo X,
    3. Lu X,
    4. Duan J,
    5. Xu G
    . Metabolomics study of diabetic retinopathy using gas chromatography-mass spectrometry: a comparison of stages and subtypes diagnosed by Western and Chinese medicine. Mol Biosyst 2011;7:2228–2237
    OpenUrlCrossRefPubMed
  36. ↵
    1. Schwartzman ML,
    2. Iserovich P,
    3. Gotlinger K, et al
    . Profile of lipid and protein autacoids in diabetic vitreous correlates with the progression of diabetic retinopathy. Diabetes 2010;59:1780–1788
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Serlin Y,
    2. Levy J,
    3. Shalev H
    . Vascular pathology and blood-brain barrier disruption in cognitive and psychiatric complications of type 2 diabetes mellitus. Cardiovasc Psychiatry Neurol 2011;2011:609202
    OpenUrlPubMed
  38. ↵
    1. Hogan SR,
    2. Phan JH,
    3. Alvarado-Velez M, et al
    . Discovery of lipidome alterations following traumatic brain injury via high-resolution metabolomics. J Proteome Res 2018;17:2131–2143
    OpenUrlCrossRef
  39. ↵
    1. Jenkins AJ,
    2. Joglekar MV,
    3. Hardikar AA,
    4. Keech AC,
    5. O’Neal DN,
    6. Januszewski AS
    . Biomarkers in diabetic retinopathy. Rev Diabet Stud 2015;12:159–195
    OpenUrl
PreviousNext
Back to top
Diabetes: 69 (10)

In this Issue

October 2020, 69(10)
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by Author
  • Masthead (PDF)
Sign up to receive current issue alerts
View Selected Citations (0)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word about Diabetes.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Circulating Metabolites and Lipids Are Associated to Diabetic Retinopathy in Individuals With Type 1 Diabetes
(Your Name) has forwarded a page to you from Diabetes
(Your Name) thought you would like to see this page from the Diabetes web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Circulating Metabolites and Lipids Are Associated to Diabetic Retinopathy in Individuals With Type 1 Diabetes
Viktor Rotbain Curovic, Tommi Suvitaival, Ismo Mattila, Linda Ahonen, Kajetan Trošt, Simone Theilade, Tine W. Hansen, Cristina Legido-Quigley, Peter Rossing
Diabetes Oct 2020, 69 (10) 2217-2226; DOI: 10.2337/db20-0104

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Add to Selected Citations
Share

Circulating Metabolites and Lipids Are Associated to Diabetic Retinopathy in Individuals With Type 1 Diabetes
Viktor Rotbain Curovic, Tommi Suvitaival, Ismo Mattila, Linda Ahonen, Kajetan Trošt, Simone Theilade, Tine W. Hansen, Cristina Legido-Quigley, Peter Rossing
Diabetes Oct 2020, 69 (10) 2217-2226; DOI: 10.2337/db20-0104
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Research Design and Methods
    • Results
    • Discussion
    • Article Information
    • Footnotes
    • References
  • Figures & Tables
  • Suppl Material
  • Info & Metrics
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Low-Frequency Genetic Variant in the Hepatic Glucokinase Gene Is Associated With Type 2 Diabetes and Insulin Resistance in Chinese Population
  • Baseline Assessment of Circulating MicroRNAs Near Diagnosis of Type 1 Diabetes Predicts Future Stimulated Insulin Secretion
  • Polycystic Ovary Syndrome and Risk of Type 2 Diabetes, Coronary Heart Disease, and Stroke
Show more Genetics/Genomes/Proteomics/Metabolomics

Similar Articles

Subjects

  • Complications-Retinopathy

Navigate

  • Current Issue
  • Online Ahead of Print
  • Scientific Sessions Abstracts
  • Collections
  • Archives
  • Submit
  • Subscribe
  • Email Alerts
  • RSS Feeds

More Information

  • About the Journal
  • Instructions for Authors
  • Journal Policies
  • Reprints and Permissions
  • Advertising
  • Privacy Policy: ADA Journals
  • Copyright Notice/Public Access Policy
  • Contact Us

Other ADA Resources

  • Diabetes Care
  • Clinical Diabetes
  • Diabetes Spectrum
  • Scientific Sessions Abstracts
  • Standards of Medical Care in Diabetes
  • BMJ Open - Diabetes Research & Care
  • Professional Books
  • Diabetes Forecast

 

  • DiabetesJournals.org
  • Diabetes Core Update
  • ADA's DiabetesPro
  • ADA Member Directory
  • Diabetes.org

© 2021 by the American Diabetes Association. Diabetes Print ISSN: 0012-1797, Online ISSN: 1939-327X.