INTRODUCTION
Parkinson’s disease (PD) is the second most frequently diagnosed neurodegenerative condition [
1,
2]. It predominantly affects people over the age of 60, with approximately 1% of this age group affected globally. The likelihood of developing PD increases with age. This disease progresses over time and is characterized by the deterioration of dopaminergic neurons in the nigrostriatal pathway, which leads to distinctive motor function–related symptoms such as rigidity, tremors, and bradykinesia. There is increasing evidence that the development of PD can begin up to two decades before motor symptoms become evident. Imaging studies and pathological examinations indicate that degeneration in the nigrostriatal region can be identified 5 to 10 years before motor symptoms become clinically apparent [
3,
4]. In the phase leading to noticeable neurodegeneration, PD is believed to undergo several stages that collectively form a molecular prodrome [
4,
5]. Despite the demonstrated efficacy of numerous compounds in laboratory or animal models, to date, none have been shown to be effective in altering the progression of PD in clinical trials.
Insulin resistance is a common metabolic disorder that is commonly linked with type 2 diabetes mellitus (DM) [
6]. This condition occurs when the body’s cells become less responsive to insulin, a hormone essential for regulating blood sugar levels [
7]. The implications of insulin resistance extend beyond DM and can affect a wide range of metabolic issues; it is closely associated with various conditions or poor prognoses, such as hypertension, dyslipidemia, liver diseases, cardiovascular diseases, certain types of cancer, obesity, and inflammatory and infectious diseases [
7-
11]. In recent years, the association between insulin and PD has attracted increasing research attention. Insulin, which is widely known for its role in managing blood glucose levels, also appears to have a protective effect on the brain. Insulin receptors are present in brain regions such as the basal ganglia and substantia nigra [
12]. Growing evidence indicates that insulin sensitivity and resistance are crucial for maintaining neuronal health and growth, supporting dopamine-related neural transmission, and preserving synaptic connections in the brain [
13].
The triglyceride-glucose (TyG) index, which is calculated by fasting triglyceride (TG) and fasting blood glucose (FBG) levels, serves as a simple and practical surrogate marker for insulin resistance [
14]. This index is easy to use and cost-effective, particularly in settings where more direct and complex measurements of insulin resistance are not readily available. It is a valuable tool for assessing metabolic health in various clinical settings and for identifying individuals at risk of developing complications associated with insulin resistance [
9,
15]. However, the association between insulin resistance and the risk of incident PD has rarely been investigated, and studies using repeatedly measured parameters in the general population are limited. We hypothesize that a higher TyG index score is associated with the development of PD. Our study aimed to investigate the association of the TyG index with the incidence risk of PD in a longitudinal setting in the general population.
MATERIALS & METHODS
- Data source
This study sourced data from the National Health Insurance Service (NHIS) Health Screening Cohort (HEALS) database. As a government program, the NHIS provides health insurance to nearly 97% of the Korean population. It also provides a nationwide free health screening program every 2 years for all South Korean adults aged 40 years and over. The Medical Aid program, an affiliate of the NHIS, attends to 3% of the population not covered by the NHIS. Our study used the NHIS-HEALS cohort database for 2002–2019 [
16]. The NHIS-HEALS database includes measurements of blood pressure, body mass index (BMI), and blood biochemistry; the results of a self-administered questionnaire on medical history; and lifestyle records for smoking, alcohol consumption, and physical activity. Health claim data covering all hospital visits, diagnoses, surgeries, medical procedures, and prescriptions of participants from 2002 to 2019 are also included. Diagnoses at each hospital visit were recorded based on the International Classification of Disease, Tenth Revision (ICD-10). Demographic information such as sex, age, and household income was provided, and data regarding participants’ health claims, insurance coverage maintenance, and deaths were available up to December 31, 2019.
- Study population
From the NHIS-HEALS database, we included participants aged 40 years and older who participated in the national health screening program during the baseline years of 2009 and 2010. Among these 362,285 potential participants, those for whom data on demographic details, lifestyle, and laboratory findings were missing (
n=9,047) were excluded from the study. The washout period extended from 2002 to the index date, during which patients with a history of PD were excluded (
n=2,388). Participants with a follow-up duration of less than 1 year (
n=56) (to avoid reverse causality or association) and those with fewer than 3 repeated measurements (
n=40,773) were excluded. After applying these inclusion and exclusion criteria, the final cohort for analysis comprised 310,021 participants (
Figure 1).
- Data collection and definitions
Details of the participants’ age, sex, BMI, waist circumference, household income, and lifestyle habits (smoking status, alcohol consumption, and regular physical activity) were collected through self-report questionnaires. BMI was calculated as weight (kg) divided by the square of height (m
2). Household income was categorized using quantiles of individuals’ health insurance premiums, with those in the ninth decile and above considered high income. Smoking status was categorized into never, former, and current smokers. The frequency of alcohol consumption was defined as the number of times alcohol was consumed per week (0, 1–2 times, 3–4 times, and ≥5 times). The frequency of regular physical activity was based on the number of days the participants engaged in exercise per week (0, 1–4 days, and ≥5 days). Biochemical measurements, including liver enzyme, lipid, and FBG levels, were collected from the health screening laboratory results. Hypertension, DM, dyslipidemia, renal disease, and liver disease were considered comorbidities, and the Charlson comorbidity index (CCI) was used to determine the burden of covariates. Information on the use of statins and antidiabetic medications (alpha-glucosidase inhibitors, sulfonylureas or meglitinides, dipeptidyl peptidase-4 [DPP-4] inhibitors, glucagonlike peptide-1 [GLP-1] agonists, sodium-glucose cotransporter-2 [SGLT2] inhibitors, biguanides, insulin, and thiazolidinedione) was collected. Detailed definitions for these can be found in the
Supplementary Materials (in the online-only Data Supplement) [
17-
21].
- Calculation of the TyG index
The TyG index was calculated as ln ([TG level]×[FBG level]/2) [
9,
22]. In this study, the TyG index was considered a time-dependent covariate throughout the follow-up period. For further analysis, the average of repeated measures of the TyG index, calculated using values from at least three repeated measurements to reduce bias in the average value, was also utilized as a variable.
- Outcome
Individuals were identified as having PD based on at least two or more related claims, with the initial date of diagnosis being noted. A diagnosis of PD was determined via the ICD-10 code of G20 and a reimbursement code of V124 for rare intractable diseases (RIDs) used by neurologists, neurosurgeons, or specialists in rehabilitation medicine, with a minimum of one annual claim for both hospital admissions and outpatient visits, along with a prescription history for any PD medication, such as amantadine, anticholinergics, catechol-O-methyltransferase inhibitors, dopamine agonists, or carbidopa/levodopa, selegiline, and rasagiline. To exclude cases of secondary parkinsonism, individuals with diagnoses of both PD (G20) and parkinsonism (G21–26) were not included as an incidence of PD [
23]. Participants who have an RID must have their diagnoses confirmed by a physician using the standard diagnostic criteria provided by the NHIS. Following a physician’s evaluation, the health care facility also examines the diagnosis prior to submitting it to the NHIS. This structured procedure guarantees the reliability of the data related to RIDs. The date of diagnosis for PD was considered the first prescription of anti-PD medication referred to from relevant ICD-10 codes on the claim record. Follow-up was carried out until December 31, 2019, death, or the first occurrence of PD.
- Statistical analysis
Comparisons between groups based on quartiles of the TyG index were made via one-way analysis of variance for continu-ous variables and chi-square tests (or Fisher’s exact test) for categorical variables. Survival curves for the time-to-event outcomes were plotted via Kaplan‒Meier curves, and a log-rank test was used to compare the survival curves across TyG index groups. To explore the linear relationship between the TyG index per standard deviation (x-axis) and the incidence of PD (y-axis), restricted cubic splines were applied. The optimal change point in the spline curve analysis was estimated via a regression model with piecewise linear relationships.
To evaluate the incidence risk of PD in relation to repeated measurements of the TyG index during the follow-up period, a time-dependent Cox proportional hazards model was applied. The participants were divided into 3 groups based on tertiles (T1, T2, and T3) of the average TyG index during the follow-up period. To determine the risk of PD according to quartile groups, a conventional Cox proportional hazards model was used. The proportionality of the hazard assumption was evaluated via the Grambsch and Therneau test of Schoenfeld residuals, which yielded satisfactory results.
The results of time-dependent Cox regression and conventional Cox regression analyses are presented as hazard ratios (HRs) and 95% confidence intervals (CIs) for an unadjusted model, Model 1, and Model 2, depending on the adjustment of covariates. Model 1 was adjusted for age and sex, whereas Model 2 was adjusted for Model 1 plus BMI, household income, smoking status, alcohol consumption, regular physical activity, hypertension, DM, renal disease, liver disease, and CCI. Antidiabetic medication was further adjusted only for the DM cohort. Blood biomarkers, such as aspartate transaminase and alanine transaminase levels and liver disease, were not adjusted for in Model 2 because of multiple collinearity factors. For covariates, in cases where participants underwent multiple health check-ups from 2009 to 2019, data from their latest examination were used in the statistical analysis. Subgroup analyses were performed according to the presence of DM. Sensitivity analyses regarding the association of the TyG index with PD were performed according to demographic data, lifestyle habits, and covariates, suggesting p values for interactions. To confirm the association between anti-diabetic medication and the risk of PD, further analysis was conducted in the cohort that received anti-diabetic medication at least once during the follow-up period from the time of inclusion. In the case of combination antidiabetic drugs, each drug was considered to have been taken. In addition, a time-dependent Cox proportional hazards model was applied according to the duration of each antidiabetic medication. All the statistical analyses were conducted via SAS version 9.4 (SAS Inc.) and R version 4.2.1 (R Foundation for Statistical Computing), with statistical significance defined as a two-sided p value <0.05.
RESULTS
- Baseline characteristics of the participants
The number of measurements repeated during the follow-up period is described in
Supplementary Table 1 (in the onlineonly Data Supplement), and the characteristics of the variables for each year are described in
Supplementary Table 2 (in the online-only Data Supplement).
Table 1 presents the baseline characteristics of the entire cohort divided into 3 groups on the basis of the tertiles of the average TyG index score (T1: <9.129; T2: 9.129–9.539; and T3: >9.539). The members of the T2 group were older than those in the other groups were. The members of the T3 group were more likely to be male and more likely to be obese. The income level of the T3 group was lower than that of the other groups. Members of the T3 group were also more likely to be smokers and consumers of alcohol and exercise less frequently (≥ 5 days/week). With respect to laboratory data, the levels of liver enzymes, total cholesterol, triglycerides, and FBG were highest in the T3 group, and the proportions of those with comorbidities, including hypertension, DM, dyslipidemia, renal disease, liver disease, and a CCI score of 2 or more, were significantly greater in the T3 group (
Table 1).
- Relationship of the TyG index with incidence risk for PD
During a median of 9.64 years (interquartile range 8.72–10.53), 4,587 individuals (1.5%) experienced PD. Survival curves depicting the incidence of PD across tertiles of the average TyG index score are presented in
Figure 2. The incidence of PD depended on the TyG index tertiles in the entire cohort (logrank test:
p<0.001) and the non-DM cohort (
p=0.002). In contrast, the incidence of PD did not depend on TyG index tertiles in the DM cohort (
p=0.122).
For the multivariable time-dependent Cox proportional hazards model with repeated measures of average TyG index scores, a per-unit increase in TyG index score significantly increased the risk of PD in the entire cohort (HR: 1.062; 95% CI 1.007–1.119). In the sensitivity analysis, repeated measures of average TyG index scores were associated with the incidence risk of PD in the non-DM cohort (HR: 1.093; 95% CI 1.025–1.165). Repeated measures of average TyG index scores were not associated with the incidence risk of PD in the DM cohort (HR: 0.990; 95% CI 0.902–1.087) in fully adjusted multivariable models (
Table 2 and
Supplementary Table 3 in the online-only Data Supplement). In a subgroup analysis, no statistically significant interaction was found according to demographic data or covariates except for sex (
Figure 3). Additionally, the multivariable Cox proportional model for average TyG index tertiles during follow-up is detailed in
Table 3 and
Supplementary Table 4 (in the onlineonly Data Supplement). The highest tertiles of the TyG index were positively associated with the incidence of PD in the entire cohort and non-DM cohort but not in the DM cohort.
Restricted cubic spline analysis (
Figure 4) revealed a clear, nonlinear and increasing trend (U- or J-shaped) in the risk of PD as measured by the TyG index per standard deviation in the entire cohort and in the DM and non-DM cohorts.
In the sensitivity analysis, participants whose TyG index was less than three were analyzed (
Supplementary Tables 5 and
6 in the online-only Data Supplement). A multivariate time-dependent Cox proportional hazards model of average TyG index scores and TyG index tertiles revealed a consistent association of the risk of PD with the TyG index in all participants and non-DM participants with more than one measurement. However, there was no correlation between the TyG index and the risk of PD in patients whose PD was measured less than 2 times.
Detailed information on DM medication is described in
Supplementary Tables 7–
10 (in the online-only Data Supplement). Among participants taking antidiabetic medications, the risk of PD was greater in the insulin group (HR: 1.161; 95% CI 1.005–1.340) than in the sulfonylurea group (
Supplementary Table 11 in the online-only Data Supplement).
DISCUSSION
The key findings of our study were that the TyG index was associated with the incidence risk of PD in a general population based on a time-dependent analysis of the TyG index and a conventional Cox regression analysis with averages of the re-peatedly measured values of the index. This association was evident for the entire cohort and the non-DM cohort, and the relationship between the TyG index and the incidence risk of PD exhibited a J shape regardless of the accompanying DM history.
The TyG index is linked to several health conditions with respect to disease presence and progression and related adverse events. For example, increased TyG index scores are correlated with a heightened incidence of coronary artery disease, cerebrovascular disease, and peripheral arterial disease [
24]. Elevated TyG index scores can predict the progression of coronary artery atherosclerosis and calcification [
15]. Moreover, a previous study established a significant association between the TyG index and all-cause and cardiovascular mortality, particularly in young and middle-aged patients [
25]. Our study presents additional information regarding the association between the TyG index as an indicator of insulin resistance and the incidence risk of PD in a general population, given the large sample size and longitudinal setting.
Our study revealed that the TyG index was associated with an increased risk of PD in the entire population and the non-DM population but not in the DM population. Although DM is a representative disease accompanied by insulin resistance, the relationship between the presence of DM and the incidence risk of PD remains controversial. Prospective studies indicate that the link between DM and PD may be less strong, with type 2 DM patients showing an approximately 40% greater likelihood of developing PD [
26,
27]. Additionally, case-control studies conducted in Scandinavian and Asian populations suggest that type 2 DM is associated with an increased risk of PD [
28,
29]. While most research supports this connection, several studies have reported either no link [
30,
31] or an inverse relationship between type 2 DM and PD [
32]. Our study revealed no association between TyG index scores and the incidence risk of PD in the DM population. These inconsistencies may be the result of differences in research methods and residual confounding factors, such as how PD diagnoses are obtained, the use of other medications, and the presence of additional medical conditions that are common among people with DM. For the non-DM population, even in those without DM, PD patients with dementia were significantly more likely to exhibit insulin resistance than PD patients without dementia were [
33]. In contrast, previous research based on the Nurses’ Health Study and Health Professionals Follow-up Study revealed that plasma levels of insulin resistance–related metabolites did not contribute to the risk of PD [
34]. Our study, which was conducted in a longitudinal setting of a general population, suggests that insulin resistance is positively associated with the incidence risk of PD in a non-DM population.
Our study revealed evidence of a nonlinear association (a J-shaped trend) between the TyG index and the incidence risk of PD. This relationship was not uniform across all the ranges. In a previous study, the association between the TyG index and the incidence of atrial fibrillation in a general population without known cardiovascular disease showed U- or J-shaped trends [
35]. In a nationwide cohort study, the TyG index was found to have a U- or J-shaped relationship with all-cause and cardiovascular mortality in patients with DM [
36]. The results of these previous studies can also be applied to our findings regarding the incidence risk of PD. In other words, because the TyG index is composed of both TG and FBG levels, it is difficult to rule out the possibility that if the TyG index score is very low, it may be associated with a relatively poor health condition. Specifically, low TyG index scores may signify optimal metabolic health, characterized by robust insulin sensitivity and reduced lipid levels. However, excessively low TyG index scores could signal an underlying health issue, such as malnutrition or a genetic predispo-sition, which could paradoxically increase cardiovascular risk. In support of these hypotheses, lower TG levels are associated with increased motor performance in PD patients [
37]. According to a meta-analysis, high TG levels are protective factors for the pathogenesis of PD [
38]. The J-shaped trends in our study therefore support the findings of previous studies.
The relationship between antidiabetic medications and PD has been of increasing interest in recent years. A populationbased cohort study with diabetic patients revealed that the use of DPP-4 inhibitors and/or GLP-1 mimetics was associated with a lower risk of PD [
39]. In another large population-based study, thiazolidinediones, meglitinides, GLP-1 analogs, DPP-4 inhibitors, and SGLT2 inhibitors were associated with a lower risk of PD than metformin was [
40]. The possible mechanisms are anti-inflammatory effects and neuroprotective effects of antidiabetic drugs [
12,
39-
41]. In our study, in the group of patients with diabetes, the use of diabetes medications—alpha-glucosidase inhibitors, sulfonylurea or meglitinide, DPP-4 inhibitors, or GLP-1 agonists, or SGLT2 inhibitors, biguanides, and insulin—did not affect the risk of PD compared with the risk in the group not taking diabetes medications. On the other hand, insulin was associated with an increased risk of PD compared with sulfonylurea. However, evaluating the association of diabetes medication with PD may have limited the design of this study.
Although our study was not mechanistic, several plausible hypotheses can be made regarding the association between the TyG index and the incidence risk of PD. Age is a primary risk factor for PD, and aging typically involves a reduction in the sensitivity of insulin receptors outside of the brain. Research has shown that the mRNA levels of insulin receptors in the brain, specifically in the hypothalamus, cortex, and hippocampus, also decrease with age. This reduction contributes to a condition known as chronic secondary hyperinsulinemia [
42]. However, a natural decline in insulin signaling with age may be more pronounced in individuals with PD. Research has revealed a significant reduction in insulin receptor mRNA in the substantia nigra pars compacta of PD patients, along with greater insulin resistance, than in individuals of the same age without PD [
43]. Alpha-synuclein may contribute to impaired insulin signaling in PD by improperly activating the PI3K/AKT/mTORC1 pathway and triggering the activation of c-Jun N-terminal kinase. This mechanism is compounded by a weakening of normal responses to insulin and insulin-like growth factor 1 (insulin resistance), leading to disrupted AKT homeostasis and diminished protective effects from FoxO activation and glycogen synthase kinase-3B inactivation [
12]. As a result, insulin resistance may be associated with intensification of alpha-synuclein pathology and the loss of dopaminergic neurons. Other studies indicate that alpha-synuclein–induced insulin resistance can promote further aggregation of alpha-synuclein, creating a vicious cycle of worsening pathology that is also observed in Alzheimer’s disease models with amyloid-beta and tau proteins [
12].
We acknowledge several limitations of our study. First, our findings may not be generalizable to different ethnic groups, as the study exclusively involved a Korean population. Second, despite multiple assessments of the TyG index to enhance reliability, the retrospective nature of the study limits the establishment of a causal relationship. Third, participants without at least three TyG index measurements during the follow-up period were excluded, which could bias the study results. Further analysis including patients with one or two TyG index measures also revealed an association between the TyG index and PD risk, with a slightly greater HR than when only participants with three or more measures were analyzed. We cannot explain why this difference occurred in this study, but it may be due to differences in the population. Fourth, the reliance on health screening data from a general population means that key PD-related imaging biomarkers, such as the results of a beta-CIT–PET study, were not included. Fifth, hemoglobin A1C and homeostatic model assessment of insulin resistance data were not available for our study cohort. Sixth, the accuracy of blood-based measurements of triglycerides and glucose in this cohort may be inaccurate because of a lack of standardization among laboratories. Additionally, long-term dietary or nutritional status can also impact these tests. Finally, because this cohort was constructed between 2002 and 2019 and the median follow-up was 9.64 years, recent antidiabetic medications, including DPP-4 inhibitors, GLP-1 agonists or SGLT2 inhibitors, were not commonly used in the baseline dataset.
- Conclusion
Our study demonstrated that increased TyG index scores were nonlinearly associated with the incidence risk of PD in the general population. Future research with direct measurement of insulin resistance combined with nutritional data in a large population longitudinal study would be helpful to further investigate the impact of insulin resistance on the occurrence of PD.