Original research

Exploring the gender difference in type 2 diabetes incidence in a Swiss cohort using latent class analysis: an intersectional approach

Abstract

Introduction Type 2 diabetes is multifactorial and influenced by the intersection of gender-related variables and other determinants of health. The aim of this study was to highlight the intersectional social position of the participants and disentangle its role from administrative sex in predicting the development of type 2 diabetes.

Methods Using CoLaus|PsyCoLaus study, a Swiss single-centre prospective cohort initiated in 2003 and including 6733 participants (age 35–75 years; 54% women) at baseline, we conducted latent class analyses using gender-related variables (eg, risk-taking behaviours, gender roles represented by employment status, etc) and socioeconomic determinants at baseline (2003–2006) to construct intersectional classes and we tested their association with the development of type 2 diabetes at follow-up (2018–2021).

Results Of the 6733 participants enrolled at baseline, 3409 were included in our analyses (50.6%). Over a median follow-up time of 14.5 years, 255 (7.5%) participants developed type 2 diabetes, of which 158 men (62.0%). We identified seven latent classes highlighting different intersectional social position groups (ie, young, fit, educated men (N=413), non-White physically inactive men and women (N=170), highly qualified men, former or current smokers (N=557), working women living alone (N=914), low qualified working men with overweight (N=445), women with obesity, low education and low qualified job or housewives (N=329), low educated retired participants (N=581)). Using the class labelled as ‘young, fit, educated men’ as reference, the risk of incident type 2 diabetes was higher in all other classes (adjusted OR values between 4.22 and 13.47). Classes mostly feminine had a more unfavourable intersectional social position than that of the predominantly masculine classes. The corresponding OR increased in sex-adjusted regressions analyses.

Conclusions We observe cumulative intersectional effects across behavioural and socioeconomic profiles with different risks of developing type 2 diabetes emphasising the deleterious effect of a feminine gender profile. These patterns are only partly captured by traditional sex-stratified analyses.

What is already known on this topic

  • Sex differences in epidemiology, treatment and outcomes of type 2 diabetes exist.

  • Environmental exposure and socioeconomic position differences lead to disparities in health and health-related behaviours and influence the clinical presentation, development and predisposition of type 2 diabetes.

  • Socioeconomic position and health-related behaviours are not equally distributed among women and men.

What this study adds

  • The intersectional social position influences the risk of developing type 2 diabetes beyond administrative sex and unfavourable intersectional social position among predominantly feminine classes emphasises a deleterious effect of the feminine gender profile.

How this study might affect research, practice or policy

  • The portion of risk attributable to the gender profile incorporating the notion of intersectionality, beyond the administrative sex, is to be included into prevention strategies and risk score development.

Introduction

In 2019, 9.3% of the adult population aged 20–79 years were living with type 2 diabetes worldwide1 with a rapid prevalence increase in regions facing an epidemiological transition. Three major risk factor categories for type 2 diabetes are recognised: biological risk factors such as unfavourable genetic and epigenetic modifications (eg, related to maternal or paternal obesity) and hormonal status (eg, high testosterone levels in women or low sex hormone-binding globulin levels),2 3 a low socioeconomic status2 4 5 and cardiometabolic risk factors (eg, obesity,6 smoking,7 physical inactivity,8 depression9). In vivo, these categories are not mutually exclusive, and individuals find themselves at their intersections.10

Previous literature reported sex differences in type 2 diabetes epidemiology: worldwide and in high-income countries, men display an approximately 1-3-1.5-fold higher prevalence of type 2 diabetes than women.11–13In high-income Western countries in 2019, an age-standardised prevalence is of 7.3% in men and 5.3% in women,14 although women are predominant among youth-onset diabetic patients.3 Sex-disaggregated data on type 2 diabetes in Switzerland are scarce, but one study reported an age-adjusted prevalence of 7.8% in men and 5.7% in women.11 Concerning socioeconomic status, women tend to have a less favourable position than men (eg, lower educational level and job qualification, more often living alone and/or raising children alone)15 16 and their professional stress increased in the last decade.17 Regarding exposure to cardiometabolic risk factors, women tend to display a more favourable cardiometabolic profile and healthier health-related behaviours patterns (eg, smoking, drinking alcohol, alimentation) than men.15 18 Moreover, the probability of developing type 2 diabetes after cardiometabolic2 19–21 and socioeconomic4 risk factor exposition is higher for women and they have an excess risk of CVD compared with men exposed to the same risk factors,2 3 supporting the hypothesis of women developing type 2 diabetes at worse metabolic states than men.22 Nevertheless, available research on type 2 diabetes epidemiology mainly set hypotheses on sex differences a posteriori and study biological, cardiometabolic and socioeconomic risk factors separately failing to provide thorough explanations of the combined effects of these different categories.

Restricting research to the man/woman variable may be limiting as it entangles potential biological and social factors on one hand and prevents the integration of the other social dimensions and systemic power relations that modulate the intersectional social positions of women and men on the other hand. More precisely, intersectionality posits that individual identities and social locations such as gender, race, and class intersect and represents unique experiences that are overlooked by focusing on one identity over another.23 We assume this conceptualisation of gender as one aspect of the social positions shaping the life experience of individuals. Gender medicine research has highlighted how gender influences risk exposure, health-related behaviours and access to healthcare.24 It also defined three different levels of the gender dimension (ie, individual, interactional and institutional): as an example, risk-taking behaviours are proxy for the conformity to (masculine) gender norms on the individual level and job-related physical intensity for gender roles on the interactional level (figure 1). In recent years, this new focus on medical research challenged the sex dichotomy in how epidemiological science and knowledge are conceived and different research methods on how to integrate gender in clinical research are being developed.24–28 The authors advocate for disentangling sex and gender, illustrating how neglecting gender in its predefined sense reinforces health disparities, and arguing for robust methods to improve the reproducibility of these emerging approaches. Nevertheless, operationalisation of gender as an intersectional sociological concept remains a challenge.29 To the best of our knowledge, only a limited number of studies have delved into the multidimensional impact of gender on acute coronary syndrome30 and, more recently, metabolic syndrome.31

Figure 1
Figure 1

Gendered-related variables representing the three levels of the gender dimension.

The originality of this study is its contribution to explore the added value of a latent class analyses (LCA) approach to describe and understand the role of the intersectional social position (including multidimensional gender, sociodemographic and health-related behaviour variables) in contributing to the differences observed between women and men related to their risk of developing type 2 diabetes.

Methods

Study design

This project is a data analysis of the CoLaus|PsyCoLaus study, a single-centre prospective cohort on determinants of cardiovascular and mental disease, initiated in 2003 in Lausanne, Switzerland. Its detailed protocol has been described previously.15 32 Between 2003 and 2006, 6733 subjects (age range 35–75 years, 54% women) were randomly recruited from the population of Lausanne, located in the French-speaking part of Switzerland. Periodic resurveys of the whole cohort were conducted over an 18-year follow-up.

Patient and public involvement

No patient or public involvement in the study design.

Inclusion and exclusion criteria

The CoLaus|PsyCoLaus study initially included participants aged 35–75 who provided written informed consent and had French language ability. For this secondary analysis, exclusion criteria comprised baseline diabetes (type 1 or 2), missing information on diabetes at baseline, and missing information on diabetes at the third follow-up.

Definitions

Outcome definition

Incident type 2 diabetes was defined as having a fasting plasma glucose ≥7 mmol and/or reporting an antidiabetic drug treatment (ie, oral and/or parenteral) at third follow-up without fulfilling these criteria at baseline. Glycosylated haemoglobin (HbA1c) measurement was not available at baseline and, therefore not used as an outcome criterion.

Gender-related factors

Defining which variables available at baseline would be considered as gendered-related factors was based on the gender toolbox we developed33 to represent the three levels of gender (individual, interactional and institutionalised). These variables were also expected to differ between men and women and to have an influence on health outcomes (ie, incidence of type 2 diabetes).28

At the individual level, conformity to gender norms (use of antidepressant drug treatment) and risk-taking behaviours (including alcohol consumption, smoking, physical inactivity) was selected. At the interactional level, gender roles (represented by employment status, current job type and job-related physical intensity) and gender relations (represented by current domestic situation) were selected. Receiving social help and educational level were selected to represent the institutionalised gender (ie, the institutional level) (figure 1).

Confounding factors

We considered cardiometabolic risk factors (age categories, cardiovascular disease at baseline, abdominal obesity, high blood pressure (BP), dyslipidaemia, familiar history of type 2 diabetes, polycystic ovary syndrome (PCOS), history of gestational diabetes and menopause status) and sex as confounding factors for regression analyses.

Data collection

Following an overnight fast, participants visited Lausanne University Hospital for a physical examination, a 50 mL blood sample, and an interview with a trained nurse.

The physical examination in light clothes and barefoot included measures of weight (in kilograms to the nearest 100 g using a Seca scale (Hamburg, Germany); height (to the nearest 5 mm using a Seca (Hamburg, Germany)) height gauge); waist circumference (ie, the average of two measurements executed midway between the lowest rib and the iliac crest), and BP (measured three times using an Omron HEM-907 (Matsusaka, Japan) automated oscillometric sphygmomanometer after at least a 10 min rest in a seated position). For BP, we used the average of the last two measurements.

Overweight was defined as body mass index (BMI) ≥25 and <30 kg/m2 and obesity as BMI ≥30 kg/m2. Abdominal obesity was defined as waist circumference >102 cm for men and >88 cm for women.34

Glucose, triglycerides and high-density lipoprotein (HDL)-cholesterol were measured with a Modular P apparatus (Roche Diagnostics, Switzerland) at the clinical laboratory of the Lausanne University Hospital within 2 hours of blood collection. Low-density lipoprotein (LDL)-cholesterol levels were assessed using the Friedewald formula. Low HDL-cholesterol was defined as HDL-cholesterol <1 mmol/L for men and <1.3 mmol/L for women; high-LDL was defined as LDL-cholesterol ≥3.4 mmol/L; hypertriglyceridemia was defined as triglycerides ≥1.7 mmol/L. High BP was defined as systolic BP ≥140 and/or diastolic BP ≥90 mm Hg and/or presence of an antihypertensive drug treatment. Dyslipidaemia was defined as low-HDL and/or high-LDL and/or hypertriglyceridemia and/or presence of a hypolipidemic drug treatment. Prescribed and over-the-counter medication was collected by questionnaire.

Demographic, cardiometabolic history and lifestyle data were gathered through a questionnaire, including information on adoption status and place of birth (Switzerland or elsewhere). Other sociodemographic information retrieved were ‘White’ as self-reported race; current domestic situation (alone, monoparental family, couple living without children, couple living with children); education (‘high’ for university; ‘middle’ for secondary and high school; ‘low’ for compulsory education, apprenticeship or none); receiving social help (yes, no or does not know); current professional status (‘working at least 50%’, ‘not working or working at ≤50%’ or ‘staying at home’); current job type (‘high qualification’ for entrepreneur, liberal profession and senior management; ‘middle qualification’ for independent worker, middle management and qualified worker; ‘low qualification’ for employed worker, farmer, unqualified worker, manoeuvre; ‘not working’) and job-related physical activity (sitting, standing, carrying light load, carrying heavy load). Baseline cardiovascular disease (ie, stroke and acute coronary syndrome history) was adjudicated. Positive family history of diabetes was noted if either parent had diabetes. Relevant gynaecologic pathologies (eg, PCOS, gestational diabetes, menopause status) were documented. Lifestyle data encompassed weekly alcohol consumption (considered at risk if ≥28 units for men, ≥14 units for women), smoking status (never, former, current) and physical activity (≥20 min two times per week).

Statistical methods

For all descriptive analyses, we reported categorical variables as frequency and percentage and continuous variables as mean (SD) for normally distributed data and median (IQR) for skewed data. We used independent samples t-tests (for continuous variables) and χ2 tests (for categorical variables) to compare the distribution of baseline characteristics. We used non-parametric equivalent tests in non-normal distributions. We performed LCA, a finite mixture model identifying homogeneous groups in a diverse population using selected indicators.35 36 Variables within each latent class are independent, resulting in consistent profiles across different categorical subgroups (eg, in men, women and all age categories in this study). We used a theory-driven approach and retained as class-defining indicators the gender-related variables mentioned above and BMI as it is strongly related to socioeconomic status and lifestyle37 and cannot be reduced to a pure biological variable (online supplemental methods M1). We did not introduce the outcome (ie, incident type 2 diabetes cases) in the LCA design but indicators represented assumed determinants of its development. We fitted a series of latent class models starting from k=1 onward (where k represents the number of classes). We ensured that the smallest class size was >1.5% of the study sample as in previous research38 and assigned each individual to the class for which he or she had the highest posterior probability.36 We selected the optimal model (ie, the optimal value of k) based on model fit indices and clinical interpretability. The selected indices were the Akaike information criterion (AIC), the Bayesian information criterion (BIC), and entropy where lower values for Akaike and BIC as well as higher values for entropy indicate a better fitter model. We stopped fitting the model (ie, adding a new class) when AIC and BIC increased at the addition of a new class. The research team evaluated the interpretability and clinical coherence of the classes. For each variable category, the ratio of the class prevalence to the overall prevalence was colour coded in a heat map graphic representation and the most important differences gave their label to the classes (see online supp. table S4). To assess the relationship between class membership and the incidence of type 2 diabetes, we conducted logistic regressions analyses in a three-step process: first without adjustment variables (model 1), then adjusted for cardiometabolic risk factors (model 2), and finally adjusted for variables included in model 2 and sex (model 3).

We also conducted univariate and multivariate (non-adjusted and adjusted) logistic regression analyses of incident type 2 diabetes and sex, abdominal obesity and several socioeconomic variables to explore the magnitude of these associations compared with the relationship between class membership and incident type 2 diabetes (online supplemental table S3). As the highest rate of missing value is 2% (for physical activity) and no information gain can be expected from imputation with missing data rates below 5%,39 list wise deletion was applied. We set statistical significance at p value <0.05 and conducted statistical analyses with STATA and R softwares.40 41

Results

Baseline characteristics

From the 6733 participants who participated to the CoLaus|PsyCoLaus study at baseline, 3409 were included (50.62%). We excluded 18 (0.26%) participants with missing data for type 2 diabetes or type 1 diabetes at baseline, 119 (1.77%) participants with type 2 diabetes at baseline and 3187 (47.33%) participants with missing data for type 2 diabetes at third follow-up (due to loss to follow-up) (figure 2). Compared with included participants, participants excluded from analyses were older, more frequently men and social help recipients, and had less commonly any professional activity or a high or middle education level (online supplemental table S1).

Figure 2
Figure 2

Study flowchart.

In the final sample, 1893 participants were women (55.53%), and the mean age was 50.30 years (SD 9.75). Concerning gender-related factors, women were more frequently living alone with or without children, had more frequently a middle or low education level, received more social assistance and worked more often part time or not at all than men. Women had mostly middle qualified job, in standing position. They drank less alcohol, smoked less often, were less physically inactive and took more antidepressant drug treatment as their male counterparts. Regarding cardiometabolic risk factors, women had less frequently cardiovascular disease, high BP, dyslipidaemia, high-LDL or hypertriglyceridemia, and they had a lower BMI than men did. However, they had more frequently low-HDL levels and abdominal obesity (table 1).

Table 1
|
Baseline characteristics of participants overall and according to sex.

Type 2 diabetes cumulative incidence and relative risks

Overall, 255 (7.48%) participants developed type 2 diabetes over a median follow-up time of 14.53 years (IQR 14.40–14.77 years). The sex-specific incidence was 10.42% in men and 5.12% in women, p<0.001, and the relative risk was reduced by two-third for women compared with men (OR 0.30, 95% CI 0.19 to 0.49, p<0.001). Living in couple also had a protective effect (OR 0.67, 95% CI 0.49 to 0.91, p=0.012), whereas abdominal obesity (OR 2.48, 95% CI 1.70 to 3.63, p<0.001) or lower educational level (OR 1.86, 95% CI 1.17 to 2.96, p=0.009) increased the probability of developing type 2 diabetes (online supplemental table S2).

Latent class modelling

A seven-class model was identified as optimal according to statistical indices (AIC 58989.9, BIC 59879.4, entropy 0.83) (online supplemental table S3). Descriptive characteristics of the latent classes are available in table 2. For 1805 (52.95%) participants, the probability of belonging to the class they were assigned to was >0.85, while for 332 (9.74%) participants, this probability was <0.55, indicating more ambiguous membership. The median posterior probability ranged from 0.71 (IQR 0.58–0.82) in class 3 to 0.99 (IQR 0.92–1.00) in class 5 (data not shown). These classes were considered as clinically relevant by the authors. The ratio (r) of the prevalence of each variable category within each class to the overall prevalence is available in online supplemental table S4 (ie, higher r meaning a higher prevalence in the class than in the overall sample). According to the smallest and highest r (which identifies variables whose distribution is the most different from the whole sample), classes were labelled as follows to represent their most specific characteristics: class 1 as non-White physically inactive men and women, class 2 as highly qualified men, former or current smokers, class 3 as young, fit, educated men, class 4 as working women living alone, class 5 as low qualified working men with overweight, class 6 as women with obesity, low education and low qualified job or housewives and class 7 as low educated retired participants. We observed a different constellation of socioeconomic and behavioural factors in every latent class generated and according to the predominant sex represented in the classes.

Table 2
|
Prevalence of each categorical variable within each of the seven classes generated by latent class analysis

Regressions analyses of incident type 2 diabetes and latent classes

Class 3 (young, fit, educated men) was defined as the reference group for all regression analyses due to its favourable socioeconomic and behavioural profile compared with the other classes in relation to the risk of developing diabetes. In model 1 (ie, without adjustment), ORs were very high for the other classes (eg, OR 20.32, 95% CI 5.99 to 68.87, for class 1 non-White physically inactive men and women; OR 20.89, 95% CI 6.49 to 67.19, for class 5 low qualified working men with overweight) (figure 3). In model 2 (ie, after adjustment for cardiometabolic risk factors without integrating administrative sex), the magnitude of the ORs decreased by a twofold to threefold factor but remained statistically significant, except for class 4 working women living alone whose odds were no longer significantly different from reference group (figure 3). Model 3 (ie, adjusted as in model 2 plus for administrative sex) allows the interpretation of the ORs as associations of gender profiles with the probability of developing type 2 diabetes, independently of sex, age and other cardiometabolic risk factors. In this model, each class had significant higher OR than class 3 young, fit, educated men. More precisely and compared with model 2, the ORs for class 2 highly qualified men, former or current smokers and class 5 low qualified working men with overweight (both classes including almost only men) were attenuated while ORs increased for all other classes (containing more women than men) (figure 3).

Figure 3
Figure 3

Non-adjusted and adjusted OR with 95% CIs for incident type 2 diabetes by latent class. Model 1: non-adjusted ORs; model 2: ORs adjusted for cardiometabolic risk factors; model 3: ORs adjusted for cardiometabolic risk factors and administrative sex.

Discussion

This study uniquely assessed the role of gender as an intersectional sociological concept on the incidence of type 2 diabetes through an exploratory methodology using LCA. Overall, the cumulative incidence of type 2 diabetes in our sample was 7.48% with a 70% increased likelihood of developing diabetes for men compared with women. However, we observed a gendered distribution of the intersectional social position and a deleterious effect of the feminine gender profile.

The higher incidence of type 2 diabetes among men in our study was consistent with known European data on prevalence and burden of disease across the same period1 11 12 with a high relative risk difference. As observed in other studies,17 42 women without type 2 diabetes had a healthier cardiometabolic profile at baseline but a more unfavourable socioeconomic profile than men without type 2 diabetes. Unsurprisingly in view of these cohesive results, the seven latent classes identified correspond to social groups encountered in our regional clinical practice, with a distribution of socioeconomic risk factors reflecting their gendered distribution. As an example, the class with more former or current smokers (representing the conformity to gender norms on the individual level of the gender dimension) and high work qualification (representing gender roles on the interactional level of the gender dimension) included more men. Men were also more numerous in the class with a higher proportion of highly educated (representing institutionalised gender), physically active people of normal weight (representing conformity to gender norms on the individual level of the gender dimension). On the contrary, women outnumbered men in the class with a majority living alone while working mainly part time and in the class with a predominance of individuals with low education, low qualified jobs or housewives and with obesity (representing gender roles and gender relations on the interactional level of gender).

However, due to the independence of the variables within a class, the same gender profile was found in any member of one class, regardless of their administrative sex or age category. This allowed us to compare gender profiles across classes: each class had significant higher OR than the class of young, fit, educated men with the ORs for classes containing almost only men (ie, class 2 Highly qualified men, former or current smokers and class 5 low qualified working men with overweight) decreasing when adjusting for sex in addition to other cardiometabolic risk factors, including age. The opposite was true for women-dominated classes, reflecting the unfavourable biological profile of men. However, classes mostly feminine (class 4: OR=4.99, 95% CI 1.44 to 17.26, class 6: OR 8.91, 95% CI 2.46 to 32.31) showed globally higher ORs than male dominated classes (class 3: OR 1 (reference), class 2: OR 4.22 95% CI 1.28 to 13.94, class 5: OR 7.15 2.18 to 23.42). Compared with socially advantaged groups, people with disadvantaged intersectional positions have a higher overall risk of chronic diseases. Several hypotheses support this association: chronic stress,43 the concept of embodiment44 (how gender oppression might ‘get under the skin’ to affect the health of women and gender minorities 24), and multidirectional links between several factors including depression and obesity.42 The latter is also more prevalent in lower socioeconomic environments and a strong independent risk factor for type 2 diabetes with an analogy to countries with low sociodemographic index.1 These elements were reflected in our results: classes with a higher risk of developing type 2 diabetes than the reference represented socially disadvantaged groups. These dynamics explained the unfavourable gender profile of the predominantly feminine classes: for instance, ‘working women living alone’ represented a low socioeconomic position group in Switzerland as women are more exposed (and increasingly so) to ‘non-traditional’ risk factors (eg, stress at work, problems arranging work with family duties, major depression, etc).17

Consequently, we conclude that LCA can be an effective method for integrating gender into epidemiological data with an intersectional perspective.45

Strengths and limitations

Strengths of this study are precise data on socioeconomic status and health-related behaviours collected longitudinally. It allowed an intersectional approach that went beyond traditional sex-stratified analyses or gender scores. Furthermore, LCA revealed higher risk groups (ie, high ORs) than regression analyses incorporating administrative sex, socioeconomic and cardiometabolic risk factors and health behaviours separately. This study has several limitations. First, its single-centre design with important loss to follow-up and limited sample size prevents causal and generalisable conclusions, as the local context is inextricably linked to type 2 diabetes epidemiology. Survival analyses could not be carried out with a single follow-up, although they could have refined the analyses. However, important variables such as education, living situation and work qualification show little change in older adults and fundamental changes in results are unlikely. Second, the primary study was developed more than 20 years ago, and decisions—such as the exclusion of non-French speakers—are debatable. This study has a potential representativeness bias with more privileged participants included in the original study as shown by the characteristics of the patients excluded from this secondary analysis. Third, the absence of HbA1c measurement at baseline and systematic oral glucose tolerance test (oGTT) reflecting postprandial insulin resistance may underestimate the incidence of type 2 diabetes, particularly in women regarding oGTT.2 Finally, the lack of standardised gender measurement tool limits comparison with existing studies.

Perspectives

Better theoretical framework and operationalisation guidelines are required to improve gender-sensitive analyses in bio-medical research. The exploration of this operationalisation is essential to the integration of gender as a variable, leading to improved quality and equity of care.24 A discussion on the epistemological framework in which this research is embedded is also necessary since our beliefs about gender affects what kinds of knowledge scientists produce about sex in the first place.46 Our study also highlights uncertainty in the optimal segmentation methodology for populations with type 2 diabetes.47 Examining gender as a segmentation method can help to recognise the interconnectedness of demographics, socioeconomic factors, and health behaviours, especially in lifestyle-related chronic diseases.

Further exploration is needed for applying this approach in prevention and clinical practice, especially in underprivileged populations. The unexpected OR magnitudes underscore the necessity of integrating an intersectional approach in diverse populations/databases for comparison.

Conclusions

LCA allow the operationalisation of an intersectional approach of gender as an epidemiological risk factor for type 2 diabetes incidence beyond traditional sex-stratified analyses. Cumulative intersectional effects across behavioural and socioeconomic profiles emphasise on the deleterious effect of a feminine gender profile. Considering multifactorial aspects of gender in the evaluation of epidemiological risk factors seem to be a promising approach to better understand complex diseases such as type 2 diabetes. Prevention strategies should also account for gender to better approach unprivileged groups of the population.