Social isolation, loneliness, and incident type 2 diabetes mellitus: results from two large prospective cohorts in Europe and East Asia and Mendelian randomization.

Introduction

Type 2 diabetes mellitus (T2DM) is a widespread public health concern with high global prevalence. According to the International Diabetes Federation Diabetes Atlas Tenth Edition,1 the worldwide diabetes population quadrupled between 1980 and 2021, reaching 537 million,2 of which nearly 95% is T2DM.3 This number is projected to rise to 578 million (9.3% global prevalence) by 2030 and 700 million (10.2%) by 2045.2,4 Social isolation and loneliness are pervasive global concerns, intertwined yet distinct in their concepts. Social isolation is defined as an objective deficiency in the quantity of relationships and the frequency of interactions with family, friends, and the broader community,5,6 whereas loneliness is defined as the subjective adverse emotion of solitude or the distressing sensation of being socially isolated, regardless of objectively being socially isolated or not.7 The estimated prevalence of social isolation is 25% in community-dwelling older adults,8 and the estimated prevalence of loneliness, as reported in a recent meta-analysis, ranged from 1.8% to 12.0% in young or middle-aged adults and 4.2%–24.2% in older adults in Europe.9 Studies have shown that loneliness and social isolation are associated with activation of the hypothalamic–pituitary–adrenal axis, leading to enhanced sympathetic nervous system activity, increased release of inflammatory cytokines, and increased oxidative stress,10,11 which directly contribute to the development of T2DM. In addition, social isolation and loneliness are associated with poor weight control,12 high tobacco uses,13,14 low physical activity,15 and poor sleep quality,16,17 which indirectly increase the risk of T2DM.18,19

A significant association between loneliness and a heightened risk of T2DM has consistently been seen in recent longitudinal studies, including in cohorts in England (n = 4,112),20 Denmark (n = 24,687),21 and Norway (n = 24,024).22 Conversely, findings on the relation between social isolation and risk of T2DM have been less consistent, since the significant association between social isolation and an increased risk of T2DM was only observed in a Danish cohort study,21 but not reported in other European cohort research.20,23,24 In general, although recent research has demonstrated a positive connection between loneliness and the risk of T2DM, the relationship between social isolation and T2DM risk remains a subject of debate. The existing studies have provided promising insights into the associations of social isolation and loneliness with T2DM risk. However, the available data have several limitations, including a focus on European populations, which limits generalizability; reliance on a traditional cohort study design, which limits causal inference; and a lack of consideration for genetic predisposition to glucose metabolism or T2DM, which may be important confounders or effect modifiers. Furthermore, weaknesses such as insufficient follow-up (absence of specific event time) and less accurate T2DM diagnostic methods, may also introduce potential bias. Also, evidence from Mendelian randomization (MR) analyses, which could help to establish the causal associations between social isolation, loneliness, and the risk of T2DM, is lacking. The role played by genetic susceptibility to T2DM in the relations between social isolation, loneliness, and T2DM risk has yet to be thoroughly elucidated. Hence, despite the range of studies investigating the correlation between social isolation, loneliness, and T2DM, the existing evidence remains insufficient, demanding further evaluation.

To address these gaps, we analyzed the data on social isolation, loneliness, and incident T2DM in two large-scale prospective cohorts from Europe and Asia and conducted two-sample MR analysis to address causality of the associations. We integrated genetic susceptibilities to T2DM, insulin resistance, and beta cell function in the analysis, with the aim of bolstering the study's robustness through comprehensive consideration of genetic factors our overarching aim is to reveal connections between social isolation, loneliness, and T2DM risk, thereby enriching the landscape of T2DM prevention through a nuanced consideration of social dimensions.

Methods

Two prospective cohort studies are based on data from the UK Biobank (UKB) and the China Health and Retirement Longitudinal Study (CHARLS) databases, respectively. The UKB is a comprehensive biomedical database of population health and genetic research resources. Between 2006 and 2010, more than 500,000 participants aged 37–73 years were recruited from 22 assessment centers across the UK. Participants completed touchscreen questionnaires and physical measurements and provided blood samples at recruitment.25 Multiple follow-ups were conducted, and the database is frequently updated. The full UKB study protocol is available at https://www.ukbiobank.ac.uk/media/gnkeyh2q/study-rationale.pdf. The study was approved by the North West Multicenter Research Ethics Committee (REC reference for UK Biobank 11/NW/0382), and all participants provided written informed consent.25 The CHARLS is a continuous nationwide survey that seeks to gather comprehensive data on middle-aged and older adults to support research and inform policy decisions related to aging in China.26 In CHARLS, individuals aged 45 years and older from 150 county-level units and 450 village-level units completed baseline interviews in 2011. Follow-up interviews were conducted using face-to-face computer-assisted personal interviews (CAPI) in 2013, 2015, and 2018. Ethical approval for CHARLS was granted by the Institutional Review Board at Peking University (IRB00001052-11015), and all participants provided written informed consent.26 Additional details can be found at http://charls.pku.edu.cn/.

In the current study, 502,369 UK adults and 25,533 Chinese were initially screened. In the UKB cohort, we excluded those with missing data on loneliness or social isolation (n = 35,707), a prior diagnosis of diabetes at baseline (n = 32,562), and missing data on the polygenic risk score (PRS) for T2DM (n = 10,597), resulting in 423,503 participants in our analyses (Supplementary Figure S1). Supplementary Table S1 and Supplementary Table S2 present a comparison of baseline characteristics between enrolled participants and excluded individuals, along with those with missing exposure data, in the UKB cohort. In general, participants excluded from the study, including those with missing exposure data, tended to be older, belong to the ethnic group of white, have a higher percentage of males, and exhibit lower socioeconomic statuses (including lower household income, lower percentage of college degree attainment, and unemployment). Additionally, this group also tended to have unhealthy lifestyle patterns (such as smoking, low levels of physical activity, unhealthy diet, and poor sleep patterns), and experience poorer health conditions (with a higher prevalence of prediabetes, cardiovascular disease [CVD], hypertension, high cholesterol, chronic lung diseases [CLD], cancer, and obesity) compared with the enrolled subjects. Within the CHARLS cohort, a total of 13,800 participants were included in the final analysis after excluding individuals with missing data on social isolation or loneliness (n = 10,034), baseline diagnosis of diabetes mellitus (n = 1,105), with censoring of those lost to follow-up (n = 594) (Supplementary Figure S1). Further information on the participant selection process can be found in Appendix 1.

The exposures of interest are social isolation and loneliness. The evaluation of social isolation primarily centers on objective social connections, encompassing factors such as frequency of social interactions or engagement in social activities,5 whereas the assessment of loneliness concentrates on the subjective perception of isolation.7 To maintain consistency with previous studies conducted within the UKB and CHARLS cohorts, distinct assessment methods were employed between the two cohorts. In the current study, social isolation was evaluated based on the number of people in the household, frequency of social activities, contact with others, and marriage status (CHARLS only); whereas loneliness was evaluated by subjective feelings of loneliness and the willingness to confide in others (UKB only).

In the UKB cohort, social isolation was evaluated through responses to three questions on questionnaires: 1) “How often do you visit friends or family or have them visit you?” (1 point for the response of “once a month”, “once every few months”, “never or almost never”, or “no friends or family outside the household” and 0 point for the response of “once a week”, “2–4 times a week”, and “Almost daily”); 2) “Which of the following leisure or social activities do you engage in once a week or more often? You may select more than one of them: sports club or gym, pub or social club, religious group, adult education class, or other group activities” (1 point for the response of “none of above” and 0 point for the response of one of activities mentioned above); and 3) “Including yourself, how many people live in your household? Include those who usually live in the house such as students living away from home during term time, and partners in the armed forces or in professions such as pilots” (1 point for the response of “0” and 0 point for other numbers higher than 0). The social isolation score was calculated by summing up the points of three questions (ranging from 0 to 3). Participants were allocated into three groups according to the social isolation score: least isolated group (0 point), moderate isolated group (1 point), and most isolated group (2 or 3 points). Loneliness was assessed by two questions from UKB questionnaires: “Do you often feel lonely?” (1 point for the response of “yes” and 0 point for “no”) and “How often are you able to confide in someone close to you?” (1 point for the response of “never or almost never” and 0 point for the response of “once every few months”, “once a month”, “once a week”, “2–4 times a week”, or “almost daily”). The loneliness score was determined by adding up the scores from two questions, resulting in a range of 0–2 points. Individuals were allocated into two groups according to the loneliness score: the no-loneliness group (0 or 1 point) and the loneliness group (2 points). The field ID of the questions above is presented in Supplementary Table S5.

In the CHARLS cohort, social isolation and loneliness were also assessed through self-reported questionnaires. For the social isolation evaluation, participants received 1 point for each of the following: being unmarried (including separated, divorced, widowed, or never married), living alone, having less than weekly contact with their children (via phone, in person, or email), and not participating in any social activities over the past month (including interacting with friends, playing chess or cards, attending sports, social, or other clubs). Scores on the social isolation index ranged from 0 to 4, with higher scores indicating a higher level of social isolation. Participants with ≥2 points were allocated to the social isolation group.27 For the loneliness assessment, a single question from the Center for Epidemiological Studies Depression Scale (CESD), “In the past week, how often did you feel lonely?”, was used. Participants were classified as lonely if they said they felt lonely occasionally (1–2 days per week), frequently (3–4 days per week), or most of the time (5–7 days per week). Those respondents who said they felt lonely infrequently or never (less than one day) were categorized as not feeling lonely.28,29

Although the assessment methods were not the same in the two cohorts, the assessment methods for loneliness and social isolation in both cohorts were consistent with published studies focusing on the topic of loneliness and social isolation within the UKB and CHARLS cohorts. In the UKB cohort, researchers utilized methods above to document the association of loneliness and social isolation with heart failure,30 myocardial infarction,31 stroke,31 and in-hospital infection.32 The assessment methods utilized in CHARLS have also been applied in prior research that explored the associations of loneliness and social isolation with cognitive function27 and functional disability.29

T2DM was the primary outcome in both cohorts. In the UKB cohort, the diagnosis of T2DM was based on hospital admission data and cause of death registry records. The diagnosis of T2DM was based on hospital admission data from the Hospital Episode Statistics for England (up to 31 October 2022), Scottish Morbidity Record data for Scotland (up to 31 July 2021), and Patient Episode Database for Wales (up to 28 February 2018), and the death cause registry records from the National Health Service (NHS) Information Centre (England and Wales, up to 30 November 2022) and the NHS Central Register, National Records of Scotland (Scotland, up to 30 November 2022). T2DM was confirmed when there was a medical diagnosis in the hospital admission data or when T2DM was listed as a cause of death in the death register. We collected the diagnostic data by following the International Classification of Diseases-Tenth Revision (ICD-10); the code for T2DM was ICD E11.33 In the CHARLS cohort, DM was confirmed by the questions “Have you been diagnosed with diabetes by a doctor?” and “Are you now taking any of the following treatments to treat or control your diabetes?”. Participants who responded affirmatively to either of the two questions were considered to have a diagnosis of DM.

Covariates for this study were collected at baseline. In the UKB cohort, covariates included age, sex, race/ethnicity, assessment center, household income, educational level, employment status, smoking, alcohol drinking, physical activity, sun exposure time, healthy diet score, healthy sleep pattern, body mass index (BMI), self-reported medical conditions (history of prediabetes, CVD, hypertension, dyslipidemia, CLD, and cancer), and glucose or lipid metabolism measures from laboratory tests on blood samples. Physical activity was measured by the metabolic equivalent task (MET).30 A healthy sleep pattern was defined as a healthy sleep score ≥4.34 Prediabetes was defined as a fasting blood glucose (FBG) of 5.6–6.9 mmol/L or HbA1c of 5.7–6.4%.35 Family history of T2DM was defined by the presence of T2DM diagnosis in at least one first-degree relative. Grades of genetic risk (low, moderate, and high genetic risk) were defined according to the tertiles of the PRS for T2DM, insulin resistance, and beta cell function. The PRS for T2DM from the UKB database has been previously described.36 The calculation of PRS for insulin resistance and beta cell function is presented in Appendix 2. Information on single nucleotide polymorphisms (SNPs) utilized in the calculation of PRS for insulin resistance and beta cell function is presented in Supplementary Tables S3 and S4, respectively.

In the CHARLS cohort, covariates included demographic characteristics (age, sex, ethnicity, and residential location), socioeconomic information (household income, educational level, and employment), lifestyle (smoking, alcohol intake, physical activity), BMI, medical conditions (hypertension, dyslipidemia, heart disease, and stroke), and treatments including anti-hypertension and anti-hyperlipidemia therapies.

Specific information on covariates in the two cohorts is presented in Appendix 3 and the field ID of the questions above in the UKB is presented in Supplementary Table S5.

Our study utilized a two-sample Mendelian randomization (MR) approach using summary-level data from genome-wide association studies (GWAS) conducted on individuals of European ancestry. The exposures of MR analysis included all items of social isolation and loneliness, with instrumental variables derived from the UKB GWAS. For the outcome of T2DM, instrumental variables were identified from two independent GWAS: UKB and European Bioinformatics Institute (EBI). Detailed information on the MR design is presented in Appendix 4.

Baseline variates were presented as means ± standard error or median (interquartile range) for continuous variables and frequency (percentages) for categorical variables. Continuous variables were assessed for statistical differences using two-sample T-tests, ANOVA tests, or Mann–Whitney U tests. Categorical variables were evaluated for differences using the χ2 test. For missing data, we imputed the median for variables with a missingness rate <5%, and treated missing data as a separate category labeled “unknown” for variables with a missingness rate ≥5%. Detailed information on the covariates with a missingness rate ≥5% in the UKB and CHARLS cohorts is presented in Supplementary Table S6.

The association of exposure with incident T2DM was assessed in multivariable Cox proportional hazards and presented as hazard ratios (HR) with 95% confidence intervals (CI). The association of social isolation and loneliness with T2DM risk was evaluated in the main-effect and interaction models. Three models were generated for the analysis in the UKB cohort: Model 1 adjusted for age, sex, and race/ethnicity; Model 2 further adjusted for BMI (continuous), grip strength, assessment center, household income, college/university degree, employed, currently smoking, current alcohol drinking, physical activity, sun exposure time, healthy diet score, and healthy sleep score. Model 3 further adjusted for family history of T2DM, prediabetes, history of CVD, history of hypertension, dyslipidemia, chronic lung disease, cancer, and PRS for T2DM. In the CHARLS cohort, we also built three models for the Cox regression: Model 1 adjusted for age, sex, and ethnicity; Model 2 further adjusted for location, income, educational level, currently employed, current smoking, drinking frequency, physical activity, and BMI; Model 3 further adjusted for heart problems, hypertension, hyperlipidemia, stroke, and medication for hypertension and hyperlipidemia. The assessment of the proportional hazard's assumption for the Cox models in both cohorts is provided in Supplementary Table S7, and the P values for all covariates within the models, as well as the P-values for the global models, exceeded 0.05. The population-attributable fraction (PAF) was calculated to reflect the proportion of the events that could be avoided by eliminating the exposure. In addition, the associations of loneliness and social isolation with T2DM were further explored in participants with different levels of genetic susceptibility to T2DM (low, moderate, and high genetic risk).

Beyond the analyses above, we conducted several supplementary analyses to enhance the comprehensiveness of our study. In supplementary analyses, the cumulative risk of incident T2DM among different groups is presented in Kaplan–Meier (KM) curves with a log-rank test. As a supplement to the baseline information, baseline characteristics of the CHARLS cohort and the comparisons between participants who developed T2DM during the follow-up and those who did not in the UKB cohort are presented. Moreover, we delved into the associations of specific items of social isolation and loneliness with T2DM risk and examined the trends in T2DM risk associated with increasing scores of loneliness and social isolation. Furthermore, we explored potential interactions between social isolation and loneliness in relation to T2DM risk. In addition to the whole population, the separate and joint associations of social isolation and loneliness with the risk of T2DM, the interaction between social isolation and loneliness in the risk of T2DM, and the validation for their associations among participants with different levels of genetic susceptibility to T2DM are further investigated in participants with different glucose metabolism statuses (normal glucose tolerance and prediabetes). Furthermore, we performed subgroups for the separate and joint associations of social isolation and loneliness with T2DM risk. Subgroup analyses were performed by stratifying by age (<65 years old or ≥65 years old), sex (male or female), race/ethnicity (white or others), BMI(<25 kg/m2, 25–29.9 kg/m2, and ≥30 kg/m2), high household income (yes or no), educational level (college/university degree or others), employment status (currently employed or not), smoking status (currently smoking or not), drinking status (currently drinking or not), physical activity (MET Score less than the cohort median or above the median), healthy diet (healthy diet score ≥4 or not), healthy sleep pattern (yes or no), less sun exposure (sun exposure time <3 h or not), CVD (yes or no), hypertension (yes or no), dyslipidemia (yes or no), chronic lung disease (yes or no), cancer (yes or no), PRS for insulin resistance (low, intermediate, and high), and PRS for beta cell function (low, intermediate, and high). Sensitivity analyses were further performed to validate the results, including the independent and joint associations of social isolation and loneliness with T2DM risk and the interaction between social isolation and loneliness concerning T2DM risk. In sensitivity analyses, we validated results in new models with further adjustment for physical measurements of blood pressure (systolic and diastolic blood pressure) and laboratory tests on glucose metabolism (FBG and HbA1C), and lipid metabolism (total cholesterol (TC), total triglycerides (TG), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and lipoprotein a (Lp (a))) to evaluate the independence of the findings from the established biological pathways for T2DM. Further, we excluded participants with a family history of T2DM to eliminate hereditary effects and those with a personal history of cancer due to shortened survival commonly associated with cancer. At last, we excluded those diagnosed with T2DM within the first year of follow-up to mitigate the potential impact of reverse causality. Also, we further utilized the Fine–Gray sub-distribution hazard model to account for the potential impact of mortality as a competing event and performed mediation analyses with bootstrapping (2000 replications) to confirm significant mediating covariates between social isolation, loneliness, and T2DM incidence.

In MR analysis, we used odds ratios (OR) and 95% CI to estimate the T2DM risk caused by increased levels of loneliness and social isolation. For the primary MR analysis, we utilized the random-effects inverse-variance weighted (RE-IVW) method to estimate the causal effect.37 Several other MR analyses, including IVW, fixed effects IVW, MR-Egger regression, weighted median (WM), and penalized WM, were performed to assess the robustness of the results. We used the MR-Egger intercept test to assess for horizontal pleiotropy38 and Cochran's Q statistic for population heterogeneity.39 To confirm the causal effect of any single SNP, we conducted a leave-one-out (LOO) analysis by discarding each exposure-associated SNP and repeatedly performing IVW analysis. We also performed a bidirectional two-sample Mendelian randomization analysis from the reverse direction. Additional details on the MR analyses are presented in Appendix 5.

We performed analyses with the R version (R Foundation for Statistical Computing, Vienna, Austria). P < 0.05 was used to indicate statistical significance in 2-sided statistical testing. Because of the absence of adjustment for multiple comparisons, secondary and subgroup analyses should be interpreted cautiously.

The funders played no part in the design of the study, collection, analysis, interpretation of data, or the writing of the report.

Results

In the UKB cohort, the average age of the population was 56.3 years, and 44.7% of the participants were male. In the CHARLS cohort, the population had an average age of 58.3 years, with 46.7% of the participants being male. The baseline characteristics of participants from the UKB and CHARLS databases are presented in Table 1 and Supplementary Table S5, respectively. In the UKB cohort, individuals with higher levels of social isolation or loneliness tended to be male and had unhealthier lifestyle behaviors, including smoking, less physical activity, lower healthy diet scores, and unhealthier sleep patterns. They were also more likely to have obesity, prediabetes, and a history of cardiovascular disease, hypertension, dyslipidemia, chronic liver disease, or cancer (Table 1). In the CHARLS cohort, participants with social isolation and loneliness were more likely to be female, live in urban regions, have lower socioeconomic status and less physical activity, have obesity, and have comorbidities such as hypertension, heart disease, and stroke (Supplementary Table S8).Table 1Baseline characteristics grouped by social isolation and loneliness in the UKB cohort.All (n = 423,503)Social isolationLonelinessLeast isolated (n = 193,693)Moderately isolated (n = 170,641)Most isolated (n = 59,169)PvalueNo loneliness (n = 404,016)Loneliness (n = 19,487)PvalueAge, years56.3 ± 8.156.6 ± 8.156.2 ± 8.155.8 ± 7.9<0.00156.3 ± 8.155.7 ± 8.0<0.001Male, n (%)189,117 (44.7)84,575 (43.7)75,405 (44.2)29,137 (49.2)<0.001179,461 (44.4)9,656 (49.6)<0.001Ethnicity, white, n (%)406,435 (96.0)187,319 (96.7)163,380 (95.7)55,736 (94.2)<0.00114,968 (3.7)963 (4.9)<0.001Assessment center<0.0010.006 English, n (%)374,869 (88.5)170,758 (88.2)151,296 (88.7)52,815 (89.3)357,757 (88.6)17,112 (87.8) Scotland, n (%)30,878 (7.3)14,321 (7.4)12,353 (7.2)4,204 (7.1)29,383 (7.3)1,495 (7.7) Wales, n (%)17,756 (4.2)8,614 (4.4)6,992 (4.1)2,150 (3.6)16,876 (4.2)880 (4.5)Household income<0.001<0.001 Low, n (%)78,417 (18.5)25,035 (12.9)36,237 (21.2)17,145 (29.0)72,315 (17.9)6,102 (31.3) Medium, n (%)191,139 (45.1)90,955 (47)75,803 (44.4)24,381 (41.2)183,097 (45.3)8,042 (41.3) High, n (%)99,623 (23.5)51,258 (26.5)37,357 (21.9)11,008 (18.6)96,973 (24)2,650 (13.6) Unknown, n (%)54,324 (12.8)26,445 (13.7)21,244 (12.4)6,635 (11.2)51,631 (12.8)2,693 (13.8)College/university degree, n (%)142,166 (33.6)65,182 (33.7)56,671 (33.2)20,313 (34.3)<0.001137,434 (34)4,732 (24.3)<0.001Employed, n (%)250,564 (59.2)110,163 (56.9)103,166 (60.5)37,235 (62.9)<0.001239,835 (59.4)10,729 (55.1)<0.001Current smoker, n (%)43,721 (10.3)15,365 (7.9)18,901 (11.1)9,455 (16.0)<0.00140,218 (10.0)3,503 (18.0)<0.001Current drinker, n (%)387,809 (91.6)177,485 (91.6)156,150 (91.5)54,174 (91.6)0.400369,968 (91.6)17,841 (91.6)0.925Physical activity (MET)10.7 ± 4.811.2 ± 4.610.5 ± 4.99.7 ± 5.0<0.00110.7 ± 4.810.0 ± 5.2<0.001Sun exposure time, hours/day2.8 ± 1.92.8 ± 1.92.7 ± 1.92.6 ± 2.0<0.0012.8 ± 1.92.8 ± 2.10.725Healthy diet score3.0 (2.0, 3.0)3.0 (2.0, 3.0)2.0 (2.0, 3.0)2.0 (2.0, 3.0)<0.0013.0 (2.0, 3.0)2.0 (2.0, 3.0)<0.001Healthy sleep patterna, n (%)246,372 (58.2)11,7784 (60.8)97,823 (57.3)30,765 (52.0)<0.001238,629 (59.1)7,743 (39.7)<0.001Family history of T2DM, n (%)84,692 (20.0)38,160 (19.7)34,488 (20.2)12,044 (20.4)<0.00180,457 (19.9)4,235 (21.7)<0.001Prediabetes, n (%)85,864 (20.3)37,526 (19.4)35,222 (20.6)13,116 (22.2)<0.00181,428 (20.2)4,436 (22.8)<0.001History of CVD, n (%)20,547 (4.9)9,023 (4.7)8,288 (4.9)3,236 (5.5)<0.00119,022 (4.7)1,525 (7.8)<0.001History of hypertension, n (%)108,671 (25.7)48,360 (25)44,189 (25.9)16,122 (27.2)<0.001102,784 (25.4)5,887 (30.2)<0.001History of dyslipidemia, n (%)182,085 (43.0)83,200 (43.0)73,316 (43.0)25,569 (43.2)0.511173,546 (43.0)8,539 (43.8)0.017History of CLDa, n (%)6,451 (1.5)2,353 (1.2)2,778 (1.6)1,320 (2.2)<0.0015,832 (1.4)619 (3.2)<0.001History of cancer, n (%)31,968 (7.6)14,557 (7.5)13,006 (7.6)4,405 (7.5)0.28330,552 (7.6)1,416 (7.3)0.154Obesitya, n (%)93,748 (22.1)40,506 (20.9)38,849 (22.8)14,393 (24.3)<0.00188,106 (21.8)5,642 (29.0)<0.001Grip strength, kg30.8 ± 11.030.9 ± 10.930.7 ± 11.030.9 ± 11.1<0.00130.8 ± 11.030.3 ± 11.4<0.001SBP, mmHg139.4 ± 19.7139.6 ± 19.7139.2 ± 19.6138.9 ± 19.6<0.001139.4 ± 19.7138.2 ± 19.2<0.001DBP, mmHg82.2 ± 10.782.1 ± 10.682.2 ± 10.782.5 ± 10.9<0.00182.2 ± 10.782.2 ± 11.00.789FBG, mmol/L4.9 (4.6, 5.2)4.9 (4.6, 5.2)4.9 (4.6, 5.2)4.9 (4.6, 5.2)<0.0014.9 (4.6, 5.2)4.9 (4.6, 5.2)0.026HbA1C, %5.3 (5.1, 5.6)5.3 (5.1, 5.6)5.4 (5.1, 5.6)5.4 (5.1, 5.6)<0.0015.3 (5.1, 5.6)5.4 (5.2, 5.6)<0.001TC, mmol/L5.7 (4.9, 6.4)5.6 (4.9, 6.4)5.7 (4.9, 6.4)5.6 (4.9, 6.4)0.1345.6 (4.9, 6.4)5.7 (4.9, 6.4)0.427TG, mmol/L1.5 (1.0, 2.1)1.4 (1.0, 2.1)1.5 (1.0, 2.1)1.5 (1.1, 2.2)<0.0011.4 (1.0, 2.1)1.6 (1.1, 2.3)<0.001LDL, mmol/L1.4 (1.2, 1.7)1.4 (1.2, 1.7)1.4 (1.2, 1.7)1.4 (1.2, 1.6)0.4681.4 (1.2, 1.7)1.3 (1.1, 1.6)0.454HDL, mmol/L3.6 (3.0, 4.2)3.6 (3.0, 4.1)3.6 (3.0, 4.2)3.6 (3.0, 4.2)<0.0013.6 (3.0, 4.1)3.6 (3.0, 4.2)<0.001Lp (a), nmol/L21.0 (9.6, 61.6)20.9 (9.6, 61.5)21.1 (9.6, 61.9)21.0 (9.5, 61.3)0.14821.0 (9.6, 61.7)21.3 (9.6, 60.6)0.705Abbreviations: CVD, cardiovascular disease; CLD, chronic lung disease; DM, diabetes mellitus; DBP, diastolic blood pressure; FBG, fasting blood glucose; HbA1c, glycosylated hemoglobin A1c; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; Lipoprotein (a), Lp (a); MET, metabolic equivalent of task; SBP, systolic blood pressure; T2DM, type 2 diabetes mellitus; TC, total cholesterol; TG, Triglyceride; UKB, UK biobank.MET: Range 0–21; positively correlated with weekly physical activity.Sun exposure time: Range 0–24 h/day.Healthy diet score: Range 0–5 points; positively linked with the level of adherence to a healthy diet.Healthy sleep score: Range 0–5 points; positively connected with the degree of adherence to a healthy sleep pattern.aHealthy sleep pattern: the healthy sleep score ≥4; Chronic lung disease: chronic bronchitis, emphysema, chronic obstructive pulmonary disease; Obesity: BMI ≥30 kg/m2.

Baseline characteristics grouped by social isolation and loneliness in the UKB cohort.

Abbreviations: CVD, cardiovascular disease; CLD, chronic lung disease; DM, diabetes mellitus; DBP, diastolic blood pressure; FBG, fasting blood glucose; HbA1c, glycosylated hemoglobin A1c; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; Lipoprotein (a), Lp (a); MET, metabolic equivalent of task; SBP, systolic blood pressure; T2DM, type 2 diabetes mellitus; TC, total cholesterol; TG, Triglyceride; UKB, UK biobank.

MET: Range 0–21; positively correlated with weekly physical activity.

Healthy diet score: Range 0–5 points; positively linked with the level of adherence to a healthy diet.

Healthy sleep score: Range 0–5 points; positively connected with the degree of adherence to a healthy sleep pattern.

Healthy sleep pattern: the healthy sleep score ≥4; Chronic lung disease: chronic bronchitis, emphysema, chronic obstructive pulmonary disease; Obesity: BMI ≥30 kg/m2.

In the UKB cohort, 15,072 study participants were diagnosed with T2DM during a mean follow-up period of 13.5 years. Compared with those without a T2DM diagnosis, those with T2DM exhibited higher levels of social isolation and loneliness, as well as more unhealthy lifestyle factors such as smoking, low levels of physical activity, and unhealthy sleep patterns (Supplementary Table S9). In comparison to participants in the least isolated group, those in the moderately and most isolated groups had a significantly higher risk of T2DM (moderately isolated: HR 1.06, 95% CI 1.02–1.11; most isolated: HR 1.17, 95% CI 1.11–1.23) (Table 2). Similarly, individuals in the loneliness group had a significantly higher risk of T2DM compared with those in the no-loneliness group (HR 1.21, 95% CI 1.13–1.30). With increasing levels of both social isolation and loneliness, participants had a progressively and significantly higher risk of T2DM (P trend < 0.001) (Table 2). Consistent outcomes were obtained from Kaplan–Meier plots in supplementary materials (social isolation: Supplementary Figure S2; Loneliness: Supplementary Figure S3; Joint effects: Supplementary Figure S4). Further, we identified a positive dose–gradient association of social isolation score and loneliness score with T2DM risk (P trend < 0.001), and found that participants with less social contact (HR 1.14, 95% CI 1.07–1.21), living alone (HR 1.08, 95% CI 1.04–1.13), feeling lonely (HR 1.22, 95% CI 1.17–1.27), and being unwilling to confide in others (HR 1.06, 95% CI 1.01–1.11) were at significantly enhanced risk of T2DM (Supplementary Table S10). The loneliness and social isolation scores were positively associated with the risk of T2DM (P trend < 0.001) (Supplementary Table S10). A significant positive gradient of association between social isolation and T2DM risk was found in participants who did not experience loneliness (P trend < 0.001) but not in participants with loneliness (P trend = 0.391), whereas the association of loneliness with increased T2DM risk was significant regardless of social isolation level (P trend < 0.001) (Supplementary Table S11). The significant association between social isolation, loneliness, and increased risk of T2DM was observed across strata of genetic susceptibility for T2DM (social isolation: P value = 0.003 [low], <0.001 [moderate and high]; loneliness: P value = 0.003 [low], <0.001 [moderate], and 0.003 [high]), as depicted in Fig. 1, and these association was not modified by genetic susceptibility (P interaction = 0.420 for social isolation and 0.070 for loneliness).Table 2Separate and joint association of social isolation and loneliness with long-term risk of T2DM in the UKB cohort.NCases/Person-yearsModel 1 HR (95% CI)Model 2 HR (95% CI)Model 3 HR (95% CI)PAF (%) (95% CI [%])Separate effectsSocial isolation12.61 (10.93–14.28) Least isolated193,6936,051/2,632,0081.00 [Reference]1.00 [Reference]1.00 [Reference] Moderate isolated170,6416,311/2,297,8121.21 (1.17–1.25)1.10 (1.05–1.14)1.06 (1.02–1.11) Most isolated59,1692,710/784,5681.50 (1.43–1.56)1.22 (1.16–1.29)1.17 (1.11–1.23) Ptrend<0.001<0.001<0.001Loneliness8.47 (7.64–9.29) No loneliness404,01613,853/5,456,4681.00 [Reference]1.00 [Reference]1.00 [Reference] Loneliness19,4871,219/257,9201.97 (1.86–2.09)1.29 (1.20–1.38)1.21 (1.13–1.30)Joint effects14.39 (12.68–16.09)No loneliness Least isolated188,7255,776/2,565,1501.00 [Reference]1.00 [Reference]1.00 [Reference] Moderate isolated162,2975,805/2,187,0491.19 (1.15–1.24)1.09 (1.05–1.14)1.06 (1.02–1.10) Most isolated52,9942,272/704,2691.43 (1.36–1.50)1.21 (1.15–1.28)1.16 (1.10–1.23)Loneliness Least isolated4,968275/66,8581.90 (1.68–2.14)1.37 (1.19–1.57)1.29 (1.12–1.49) Moderate isolated8,344506/110,7622.07 (1.89–2.26)1.36 (1.23–1.52)1.23 (1.11–1.37) Most isolated6,175438/80,2992.38 (2.16–2.63)1.43 (1.28–1.61)1.34 (1.19–1.50)Ptrend<0.001<0.001<0.001Abbreviations: CVD, cardiovascular disease; CLD, chronic lung disease; DM, diabetes mellitus; PAF, population attributed fraction; PRS, polygenic risk score; T2DM, type 2 diabetes mellitus; UKB, UK biobank.Model 1: adjusted for age, sex, and ethnicity.Model 2: adjusted for age, sex, ethnicity, BMI, grip strength, assessment center, household income, college/university degree, employed, currently smoking, currently drinking, physical activity, sun exposure time, healthy diet score, and healthy sleep score.Model 3: age, sex, ethnicity, BMI, grip strength, assessment center, household income, college/university degree, employed, currently smoking, currently drinking, physical activity, sun exposure time, healthy diet score, healthy sleep score, family history of T2DM, prediabetes, history of CVD, history of hypertension, history of dyslipidemia, history of CLD,a history of cancer, and PRS for T2DM.aChronic lung disease: chronic bronchitis, emphysema, and chronic obstructive pulmonary disease.Fig. 1Separate (A) and joint (B) association of social isolation and loneliness with the long-term risk of T2DM among subjects with different levels of PRS for T2DM in the UKB. P interaction of genic risk subgroups for social isolation: 0.420. P interaction of genic risk subgroups for loneliness: 0.070. The analysis was performed in Model 3 (adjusted with age, sex, ethnicity, BMI, grip strength, assessment center, household income, college/university degree, employed, healthy sleep score, current smoking, current drinking, physical activity, sun exposure, healthy diet score, family history of T2DM, prediabetes, history of CVD, history of hypertension, history of dyslipidemia, history of CLD, and history of cancer). Abbreviations: BMI, body mass index; CVD, cardiovascular disease; CLD, chronic lung disease; PRS, polygenic risk score; T2DM, type 2 diabetes mellitus; UKB, UK biobank.

Separate and joint association of social isolation and loneliness with long-term risk of T2DM in the UKB cohort.

Abbreviations: CVD, cardiovascular disease; CLD, chronic lung disease; DM, diabetes mellitus; PAF, population attributed fraction; PRS, polygenic risk score; T2DM, type 2 diabetes mellitus; UKB, UK biobank.

Model 1: adjusted for age, sex, and ethnicity.

Model 2: adjusted for age, sex, ethnicity, BMI, grip strength, assessment center, household income, college/university degree, employed, currently smoking, currently drinking, physical activity, sun exposure time, healthy diet score, and healthy sleep score.

Model 3: age, sex, ethnicity, BMI, grip strength, assessment center, household income, college/university degree, employed, currently smoking, currently drinking, physical activity, sun exposure time, healthy diet score, healthy sleep score, family history of T2DM, prediabetes, history of CVD, history of hypertension, history of dyslipidemia, history of CLD,a history of cancer, and PRS for T2DM.

Chronic lung disease: chronic bronchitis, emphysema, and chronic obstructive pulmonary disease.

Separate (A) and joint (B) association of social isolation and loneliness with the long-term risk of T2DM among subjects with different levels of PRS for T2DM in the UKB. P interaction of genic risk subgroups for social isolation: 0.420. P interaction of genic risk subgroups for loneliness: 0.070. The analysis was performed in Model 3 (adjusted with age, sex, ethnicity, BMI, grip strength, assessment center, household income, college/university degree, employed, healthy sleep score, current smoking, current drinking, physical activity, sun exposure, healthy diet score, family history of T2DM, prediabetes, history of CVD, history of hypertension, history of dyslipidemia, history of CLD, and history of cancer). Abbreviations: BMI, body mass index; CVD, cardiovascular disease; CLD, chronic lung disease; PRS, polygenic risk score; T2DM, type 2 diabetes mellitus; UKB, UK biobank.

In the CHARLS cohort, 1,247 T2DM cases were documented during the mean follow-up of 5.8 years. Consistent with the UKB cohort, participants with social isolation (HR 1.22, 95% CI 1.06–1.40) and loneliness (HR 1.21, 95% CI 1.07–1.36) had a significantly higher risk of T2DM when compared with those without social isolation/loneliness (Table 3). Similar to the UKB cohort, a progressively increased risk of T2DM was observed with increasing levels of social isolation and loneliness (P trend < 0.001), and no interaction between social isolation and loneliness was observed (P interaction = 0.758) (Table 3).Table 3Separate and joint association of social isolation and loneliness with long-term risk of T2DM in the CHARLS cohort.NCases/Person-yearsModel 1 HR (95% CI)Model 2 HR (95% CI)Model 3 HR (95% CI)PAF (%) (95% CI [%])Separate effectsSocial isolation3.90 (1.03–6.76) No isolated10,816947/63,1451.00 [Reference]1.00 [Reference]1.00 [Reference] Isolated2,984302/16,9571.17 (1.02–1.33)1.19 (1.04–1.36)1.22 (1.06–1.40)Loneliness6.65 (3.17–10.12) No loneliness9,733826/56,8031.00 [Reference]1.00 [Reference]1.00 [Reference] Loneliness4,067423/23,2991.22 (1.09–1.38)1.24 (1.10–1.40)1.21 (1.07–1.36)Pinteraction0.6100.6270.758Joint effects8.41 (3.99–12.83)No loneliness No isolated8,060674/47,2021.00 [Reference]1.00 [Reference]1.00 [Reference] Isolated1,673152/9,6011.09 (0.91–1.31)1.12 (0.93–1.34)1.17 (0.97–1.40)Loneliness No isolated2,756273/15,9431.18 (1.02–1.36)1.20 (1.04–1.38)1.17 (1.02–1.35) Isolated1,311150/7,3561.39 (1.16–1.67)1.44 (1.19–1.73)1.42 (1.18–1.71)Ptrend<0.001<0.001<0.001Abbreviations: BMI, body mass index; CHARLS, China Health and Retirement Longitudinal Study; DM, diabetes mellitus; PAF, population attributed fraction; T2DM, type 2 diabetes mellitus.Model 1: adjusted for age, sex, ethnicity.Model 2: adjusted for age, sex, ethnicity, location, income, educational level, currently employed, currently smoking, drinking frequency, physical activity, and BMI.Model 3: adjusted for age, sex, ethnicity, location, income, educational level, currently employed, currently smoking, drinking frequency, physical activity, BMI, heart problems, hypertension, hyperlipidemia, stroke, and medication for hypertension and hyperlipidemia.

Separate and joint association of social isolation and loneliness with long-term risk of T2DM in the CHARLS cohort.

Abbreviations: BMI, body mass index; CHARLS, China Health and Retirement Longitudinal Study; DM, diabetes mellitus; PAF, population attributed fraction; T2DM, type 2 diabetes mellitus.

Model 1: adjusted for age, sex, ethnicity.

Model 2: adjusted for age, sex, ethnicity, location, income, educational level, currently employed, currently smoking, drinking frequency, physical activity, and BMI.

Model 3: adjusted for age, sex, ethnicity, location, income, educational level, currently employed, currently smoking, drinking frequency, physical activity, BMI, heart problems, hypertension, hyperlipidemia, stroke, and medication for hypertension and hyperlipidemia.

Regardless of participants' glucose status, the association between social isolation, loneliness, and increased risk of T2DM remained statistically significant (social isolation: P value < 0.001 [normal glucose tolerance], = 0.001 [prediabetes]; loneliness: P value < 0.001 [normal glucose tolerance], = 0.005 [prediabetes]) (Supplementary Figure S5A, Supplementary Table S12). Glucose status had significant interactions with the effects of social isolation and loneliness in T2DM risk (P interaction = 0.009 for social isolation and <0.001 for loneliness) (Supplementary Table S12). No interaction was found between social isolation and loneliness for T2DM risk among participants with normal glucose tolerance or prediabetes (Supplementary Tables S13 and S14). Joint analysis for social isolation and loneliness showed a positive gradient association with T2DM risk in both the normal glucose tolerance and prediabetes groups (P trend < 0.001) (Supplementary Figure S5B, Supplementary Table S15). Among participants with normal glucose tolerance, the association of social isolation with increased T2DM risk was significant regardless of genetic risk (P value = 0.032 [low], 0.014 [moderate], and <0.001 [high]), whereas, among participants with prediabetes, the association was significant only in those with low and intermediate genetic risk (P value = 0.016 [low], 0.006 [moderate], and 0.180 [high]) (Supplementary Table S16). Loneliness was significantly associated with enhanced T2DM risk in each group of genetic risk among participants with normal glucose tolerance (P value = 0.002 [low], 0.001 [moderate], and 0.002 [high]), but only in the low genetic risk group among participants with prediabetes (P value = 0.011 [low], 0.162 [moderate], and 0.115 [high]) (Supplementary Table S17).

Subgroup analyses for the separate and joint associations of social isolation and loneliness with T2DM risk are presented in Supplementary Tables S18–S20. No interaction was found between most of the covariates and social isolation or loneliness in the risk of T2DM, except for ethnicity (P interaction = 0.001 for social isolation and 0.014 for loneliness). The results were similar in the models that further adjusted for blood pressure, glucose metabolic, and lipid metabolic biomarkers (Supplementary Tables S21–S23). In the sensitivity analyses excluding participants with a family history of T2DM (Supplementary Tables S24–S26), cancers (Supplemental Tables S27–S29), and those diagnosed with T2DM within the first year (Supplementary Tables S30–S32), the association between loneliness and social isolation and T2DM remained significant. Moreover, in the analysis of the Fine–Gray models (Supplemental Tables S33 and S34), the results were consistent with the main findings. In mediation analyses, we confirmed that BMI, current smoking, low physical activity, and an unhealthy sleep pattern were significant mediators between social isolation, loneliness, and incident T2DM (all P value < 0.001) (Supplemental Table S35).

In the two-sample MR analyses based on UKB GWAS for T2DM, feeling lonely (OR 1.04, 95% CI 1.02–1.06, P value < 0.001), and engaging in less leisure/social activities (OR 1.03, 95% CI 1.02–1.05, P value < 0.001) were associated with an increased risk of T2DM, while more contact with friends or family (OR 0.99, 95% CI 0.98–0.99, P value = 0.046) was associated with a reduced risk of T2DM (Fig. 2 A). Furthermore, the analysis identified that participating in more sports clubs or gym (OR 0.95, 95% CI 0.93–0.97, P value < 0.001), religious activities (OR 0.98, 95% CI 0.96–0.99, P value = 0.022), and/or adult education classes (OR 0.95, 95% CI 0.90–0.99, P value = 0.042) were each associated with a reduced risk of T2DM (Fig. 2A). After validating the causal associations above in the EBI GWAS for T2DM, the causal impacts of feeling lonely (OR 2.23, 95% CI 1.15–4.32, P value = 0.017), taking less leisure/social activities (OR 3.02, 95% CI 1.61–5.64, P value = 0.001), more sports club or gym (OR 0.36, 95% CI 0.17–0.77, P value = 0.008) and more religious activities (OR 0.39, 95% CI 0.20–0.76, P value = 0.006) on the risk of T2DM were still robust (Fig. 2B). Scatter plots and LOO plots are presented in Supplementary Figures S6–S25. Estimates from MR analyses with RE-IVW, MR-Egger intercept test for horizontal pleiotropy, and heterogeneity tests with Cochran's Q statistic are presented in Supplementary Table S36 (UKB GWAS) and Supplementary Table S37 (EBI GWAS). In the reverse direction, we found evidence of a causal effect of T2DM on the tendency to participate in less leisure/social activities (pleiotropy-corrected MR-PRESSO OR 1.18, 95% CI 1.024–1.32, P value = 0.009) and a lower tendency to participate in sports club or gym activities (pleiotropy-corrected MR-PRESSO OR 0.85, 95% CI 0.74–0.96, P value = 0.014), suggesting potential bidirectional causality between the traits (Supplementary Table S38).Fig. 2MR analyses for the causal effects of loneliness and social isolation on the risk of T2DM by using IVW method based on GWAS of UKB (A) and EBI (B). Abbreviations: GWAS, genome-wide association studies; EBI, European Bioinformatics Institute; MR, Mendelian randomization; IVW, inverse-variance weighted; T2DM, type 2 diabetes mellitus; UKB, UK biobank.

MR analyses for the causal effects of loneliness and social isolation on the risk of T2DM by using IVW method based on GWAS of UKB (A) and EBI (B). Abbreviations: GWAS, genome-wide association studies; EBI, European Bioinformatics Institute; MR, Mendelian randomization; IVW, inverse-variance weighted; T2DM, type 2 diabetes mellitus; UKB, UK biobank.

Discussion

The key findings from this prospective cohort study can be summarized as follows: Loneliness and social isolation were found to be independently and jointly associated with a higher risk of developing long-term T2DM. Two-sample MR analysis provided evidence to support a causal link between social isolation, loneliness, and the incidence of T2DM. Associations between social isolation, loneliness, and increased T2DM risk remained significant regardless of genetic susceptibility to T2DM. The link between loneliness, social isolation, and an elevated risk of T2DM was seen across participants with varying baseline glucose metabolism status. There was no interaction between loneliness and social isolation with respect to T2DM risk.

The association of loneliness and social isolation with the risk of T2DM has been previously explored. In 2020, Hackett et al. conducted a cohort study involving 4,112 participants who were initially free from DM, using data from the English Longitudinal Study of Ageing database. Their study concluded that loneliness, but not social isolation, was associated with an increased risk of developing T2DM.20 Christiansen et al. explored the association between loneliness, social isolation, and T2DM risk in a Danish cohort study involving 24,687 participants and found a significant increase in T2DM risk among individuals experiencing loneliness over a 5-year follow-up period.21 In contrast to Hackett et al.'s study, Christiansen et al. also reported a significant association between social isolation and an elevated risk of T2DM. Recently, Henriksen et al. conducted a cohort study involving 24,024 participants from the Trøndelag Health Study and found that individuals experiencing loneliness had twice the risk of developing T2DM compared with those who did not report feelings of loneliness.22 However, the association between social isolation and T2DM risk was not analyzed. Overall, the research on the association between loneliness and T2DM risk appears to be more consistent than the research on social isolation and T2DM risk. Although existing studies provided novel insights into the association of loneliness and social isolation with risk of T2DM, some limitations remain. For example, previous studies mainly focused on European populations, which may constrain the broader applicability of their findings. Moreover, previous studies mainly used cohort settings, limiting causal inference, and did not consider genetic predisposition, a critical determinant of T2DM incidence.40 Additional concerns, such as inadequate follow-up periods (lacking specific event date) and less accurate diagnostic criteria for T2DM (relying solely on immediate blood glucose biomarkers or self-report questionnaires), may further weaken the solidity of previous conclusions. In the current study, we aimed to address these limitations. We incorporated both European and East Asian populations to improve generalizability of the results. We used two-sample MR analyses to allow for investigation of the causal relationship between loneliness, social isolation, and T2DM risk. Additionally, we incorporated PRS for T2DM, insulin resistance, and beta cell function into our analyses, thus factoring in genetic determinants of glucose metabolism. Furthermore, use of the UKB cohort, with its large sample size (n > 400,000), prolonged follow-up (mean duration of 13.5 years, with specific event dates), and precise T2DM diagnosis (relying on inpatient hospital records and death registers), substantially reinforced the validity of our findings. Beyond these robust methodological enhancements, our study results also exhibit several compelling features. First, our research provides additional evidence to support a link between social isolation and increased T2DM risk. Second, the MR analyses provide novel evidence for a causal relationship between loneliness, social isolation, and the heightened risk of T2DM. Third, the research found significant positive associations between loneliness, social isolation, and T2DM risk regardless of participants’ underlying genetic susceptibility to T2DM. Fourth, the study is the first to demonstrate significant associations between loneliness, social isolation, and risk of T2DM across participants with diverse baseline glucose metabolism statuses. Finally, the study showed that social isolation and loneliness have joint positive effects on T2DM risk, while also confirming the absence of an interaction effect between them in relation to T2DM risk. In summary, our study provides a more nuanced understanding of the association between social isolation, loneliness, and T2DM risk.

Distinct facets of social isolation (reduced social contact and solitary living) and loneliness (emotional loneliness and hesitance to confide) were examined separately. Intriguingly, within the UKB cohort, each of these components—living alone, reduced social contact, emotional loneliness, and reluctance to confide–exhibited a significant connection with heightened T2DM risk. These findings generally align with prior research. In 2009, a study of 8,804 participants in Germany underscored that living alone was an independent predictor of T2DM, particularly among males.41 Similarly, a Swedish cohort study involving 4,963 participants highlighted that individuals engaged in social contact and activities were less prone to T2DM development.42 The significance of feeling lonely in relation to increased T2DM risk was demonstrated in the cohort study by Hackett et al., as described earlier.20 Furthermore, a study involving 139,924 postmenopausal women revealed that social strain, characterized by unease in social interactions and reduced confiding, correlated with elevated T2DM risk.43 Beyond these indications, the current study employed the factors of social isolation and loneliness discussed above to compute the scores for both loneliness and social isolation. The outcomes revealed that as the scores for loneliness and social isolation increased, a corresponding upward trend in T2DM risk was evident. This outcome underscores a direct correlation between the degree of loneliness and social isolation and the susceptibility to T2DM, once more highlighting the potential utility of addressing these factors to mitigate T2DM risk.

To the best of our knowledge, this study is the first to use a two-sample MR analysis to evaluate the causal relationship between loneliness, social isolation, and the risk of T2DM, further reinforcing the strength of these associations. We also explored differential associations of various social activities with T2DM risk in MR analyses. Specifically, engagement in sports clubs or gyms and religious activities was associated with significant risk reductions in T2DM, whereas other social activities did not significantly correlate with the long-term risk of this outcome. This distinction may be attributed to the multifaceted nature of sports clubs or gyms, which involve not only social interactions but also physical activities known to be effective in T2DM prevention. Religious activities have been previously hypothesized to provide preventive and protective roles against the development of T2DM. In a qualitative study in British Bangladeshis, religious leaders highlighted the significant role of religious education in T2DM prevention, emphasizing the T2DM-preventive support from religious teaching regarding dietary and physical activity principles.44 In 2022, a cross-sectional study of participants in Ghana revealed that a higher frequency of engagement in religious activities was linked to reduced HbA1C levels in individuals with T2DM, providing additional evidence for the beneficial impact of religious activities on blood glucose control.45 As suggested by earlier studies carried out in Ghana, involvement in religious activities may have the potential to diminish stress and counteract the influence of stress-triggering factors on health.46,47 Furthermore, consistent engagement in religious activities may foster a sense of hope, which may lead to better-coping strategies to prevent and manage T2DM.48 In addition, it is worth noting that the additional bidirectional MR results suggest that increasing leisure/social activities, such as sports club or gym activities, are inversely related to T2DM risk, and that reducing incident T2DM may lead to additional health benefits by increasing sports/gym activities, such as weight loss.

Loneliness and social isolation were found to have no significant interaction in relation to T2DM risk. Notably, the association between loneliness and increased T2DM risk remained significant regardless of social isolation, whereas the positive link between social isolation and T2DM risk was contingent on the presence of non-loneliness. This outcome highlights the pivotal role of the subjective feeling of social disconnection compared with the objective state of social isolation. Consistent conclusions on the interaction between loneliness and social isolation have been drawn in relation to the risk of heart failure30 and mental disorders.49 This may be attributed to loneliness exerting a stronger adverse psychological impact,50 potentially leading to stress or depression that directly impacts health.51 Conversely, social isolation is an objective condition and may not evoke negative emotions in all individuals. Consequently, prioritizing psychological intervention, beyond mere alleviation of objective social isolation, is critical.

The precise mechanisms underlying the behavioral and emotional alterations caused by social isolation, and loneliness, and their subsequent impact on health outcomes are not yet fully understood.52 A previous study determined that T2DM patients with loneliness had significantly higher levels of monocyte chemoattractant protein-1, indicating a potential association between loneliness and increased inflammation.11 Research has also found that both loneliness and social isolation can trigger the activation of the hypothalamic–pituitary–adrenal axis, enhance sympathetic activity, impair parasympathetic function, and stimulate a pro-inflammatory immune response and oxidative stress.10 These are potential pathways, beyond unfavorable changes in lifestyle behaviors, through which social isolation and loneliness increase the risk of T2DM and metabolic disease.10,53

According to our study, loneliness and social isolation pose significant public health challenges that mental health professionals and clinicians should consider. Unfortunately, there is a lack of evidence-based interventions available to address and improve loneliness and social isolation,54 and preventive measures for T2DM do not currently include strategies for reducing loneliness and social isolation.55 Tackling loneliness requires a collective effort from society as a whole to provide protective interventions and support for individuals who feel lonely and socially isolated. Effective societal strategies include the establishment of community engagement initiatives aimed at fostering social interactions in public domains, such as local events, workshops, and clubs.56 Furthermore, the implementation of elderly care services, including home visits, transportation aid, and senior centers, can provide older individuals with avenues to connect with peers and engage in activities.57 Additionally, the promotion of information and communication technology can be crucial for individuals experiencing social isolation or loneliness.6 This approach has the potential to transcend social and geographical barriers by enabling convenient and affordable communication through various means regardless of time and location, especially among older people with limited mobility.6 As we navigate the post-pandemic recovery, it is crucial to establish a comprehensive response to address these challenges, supported by appropriate national policies.58

Certain limitations in this study need to be declared. First, there were limitations in the assessment of loneliness and social isolation. Specifically, there were differences between the UKB and CHARLS studies in the assessment of these exposures. CHARLS incorporated an evaluation of marital status for assessing social isolation, a component not present in the UKB. Conversely, the UKB evaluated loneliness based on the willingness to confide in others, a criterion not included in CHARLS. Although these assessment methods have been previously used, the differences might have introduced bias. Hence, we recommend that future studies employ consistent assessment methods for loneliness and social isolation across diverse populations. In addition, the exposure assessments did not rely on validated scales. Although these methods have been employed in previous research, the lack of validation introduces the possibility of bias. Second, there were limitations in the assessment of T2DM. The diagnostic methods for assessing T2DM in the UKB and CHARLS cohorts were not consistent, and each method had certain shortcomings. In the UKB cohort, the diagnosis of T2DM was based on in-patient data and the death register. Individuals diagnosed with T2DM in outpatient settings and who underwent outpatient treatment during follow-up might not have been captured, potentially leading to an underestimation of T2DM incidence. Future studies should utilize diagnostic methods that capture data from both outpatient and inpatient systems. In the CHARLS cohort, the diagnosis of DM was made via self-report on questionnaires, without subsequent corroboration from medical records. Thus, the diagnostic accuracy was less than ideal, and the precise type of DM could not be confirmed. Moreover, the DM diagnosis questionnaire in the CHARLS database was administered biennially and did not ask about a specific date of diagnosis during the follow-up. This design allowed for time-to-event data to be presented only in the form of even-numbered years within the follow-up, potentially introducing bias into the analysis. Additionally, the variation in outcome definitions between the two cohorts could introduce potential bias. Third, there was a large discrepancy in sample sizes between the two cohorts. The UKB cohort included more than 400,000 participants, whereas the CHARLS cohort included 13,800 participants. Although associations between loneliness, social isolation, and heightened T2DM risk were seen in both cohorts, the robustness of evidence within the Asian population remains comparatively limited. Future research should include larger-scale Asian cohorts. Fourth, the baseline disparities between the included and excluded participants in the UKB cohort, and not using survey weights in the CHARLS cohort analysis might limit the representation of the population and introduce potential bias. Fifth, the study did not capture changes in social isolation and loneliness during the follow-up period in either cohort, which could lead to underestimates of the observed associations. Sixth, although we controlled for several potential covariates in the analysis, residual confounding may still be present. Finally, the absence of randomized interventions limits the strength of evidence in drawing causal conclusions. Future research should include randomized controlled trials of interventions to mitigate loneliness and social isolation among individuals experiencing or at high risk of these conditions. This approach will allow for more definitive answers regarding the causality of the relationships between loneliness, social isolation, and T2DM.

In conclusion, our findings demonstrate a significant association between social isolation, loneliness, and an increased risk of T2DM. Implementing strategies aimed at mitigating these factors warrants further investigation.

Contributors

Prof. Kefei Dou attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted. Yanjun Song and Chen Zhu have directly accessed and verified the underlying data. All the authors were responsible for the decision to submit the manuscript. Concept and design: Kefei Dou and Yanjun Song. Acquisition of data: Yanjun Song (UKB and CHARLS), Chen Zhu (CHARLS), Chenxi Song (UKB), and Zhen'ge Chang (CHARLS). Formal analysis and interpretation of data: Yanjun Song, Chen Zhu, Boqun Shi, Chenxi Song, Kongyong Cui, Zhen'ge Chang, Guofeng Gao, Lei Jia, Rui Fu, Qiuting Dong, Lei Feng, Chenggang Zhu, Dong Yin. Statistical analysis: Yanjun Song and Chen Zhu. Drafting of the manuscript: Yanjun Song, Chen Zhu, and Boqun Shi. Critical revision of the manuscript for important intellectual content: JoAnn E. Manson and Kefei Dou. Obtained funding: Kefei Dou and Chen Zhu. Administrative, technical, or material support: Kefei Dou. Supervision: JoAnn E. Manson and Kefei Dou.

Data sharing statement

Researchers interested in accessing the data used in this study can apply for access to the UK Biobank by visiting their website (https://www.ukbiobank.ac.uk/enable-your-research/apply-for-access) and submitting an application that includes a research protocol summary and requested data fields. Upon approval by the UK Biobank management team and payment of applicable fees, researchers will be granted access to the database. As for the CHARLS resource, it is freely available to the public, and researchers can directly apply for data access on the website (http://charls.pku.edu.cn/). Please note that no additional data beyond what was used in this study are available.