Skip to main content

Familial co-aggregation and shared genetics of cardiometabolic disorders and traits: data from the multi-generational Lifelines Cohort Study

Abstract

Background

It is unclear to what extent genetics explain the familial clustering and the co-occurrence of distinct cardiometabolic disorders in the general population. We therefore aimed to quantify the familial (co-)aggregation of various cardiometabolic disorders and to estimate the heritability of cardiometabolic traits and their genetic correlations using the large, multi-generational Lifelines Cohort Study.

Methods

We used baseline data of 162,416 participants from Lifelines. Cardiometabolic disorders including type 2 diabetes (T2D), cardiovascular diseases, hypertension, obesity, hypercholesterolemia, and metabolic syndrome (MetS), were defined in adult participants. Fifteen additional cardiometabolic traits indexing obesity, blood pressure, inflammation, glucose regulation, and lipid levels were measured in all included participants. Recurrence risk ratios (λR) for first-degree relatives (FDR) indexed familial (co-)aggregation of cardiometabolic disorders using modified conditional Cox proportional hazards models and were compared to those of spouses. Heritability (h2), shared environment, and genetic correlation (rg) were estimated using restricted maximum likelihood variance decomposition methods, adjusted for age, age2, and sex.

Results

Individuals with a first-degree relative with a cardiometabolic disorder had a higher risk of the same disorder, ranging from λFDR of 1.23 (95% CI 1.20–1.25) for hypertension to λFDR of 2.48 (95% CI 2.15–2.86) for T2D. Most of these were higher than in spouses (λSpouses < λFDR), except for obesity which was slightly higher in spouses. We found moderate heritability for cardiometabolic traits (from h2CRP: 0.26 to h2HDL: 0.50). Cardiometabolic disorders showed positive familial co-aggregation, particularly between T2D, MetS, and obesity (from λFDR obesity-MetS: 1.28 (95% CI 1.24–1.32) to λFDR MetS-T2D: 1.61 (95% CI 1.52–1.70)), consistent with the genetic correlations between continuous intermediate traits (ranging from rg HDL-Triglycerides: − 0.53 to rg LDL-Apolipoprotein B: 0.94).

Conclusions

There is positive familial (co-)aggregation of cardiometabolic disorder, moderate heritability of intermediate traits, and moderate genetic correlations between traits. These results indicate that shared genetics and common genetic architecture contribute to cardiometabolic disease.

Background

Cardiometabolic disorders, such as type 2 diabetes (T2D), cardiovascular diseases (CVD), obesity, hypercholesterolemia, and hypertension are interrelated conditions, responsible for high mortality and disability rates worldwide. The global burden of these disorders is continuously increasing: the global population with CVD doubled from 271 million in 1990 to 523 million in 2019 [1], and the case-number of diabetes (20–79 years) is predicted to increase from 537 million in 2021 to 783 million in 2045 [2]. These cardiometabolic disorders often co-occur within individuals, suggesting a co-pathogenesis of metabolic abnormalities among various cardiometabolic disorders, which is commonly described as metabolic syndrome (MetS), a highly prevalent, multifaceted cluster of metabolic abnormalities [3,4,5,6,7,8,9].

A role of genetics is likely: several studies have observed familial aggregation of CVD [10, 11] and T2D [12], obesity [13, 14], and MetS [15,16,17], while cardiometabolic traits, such as blood pressure, fasting blood glucose, and total cholesterol, were shown to be heritable [18,19,20]. Some evidence exists that different cardiometabolic disorders co-aggregate within families, i.e., a family history of a specific cardiometabolic disorder associates with elevated risk of another cardiometabolic disorder [21,22,23,24]. However, familial aggregation, and especially co-aggregation of a full spectrum of cardiometabolic disorders, has not yet been investigated comprehensively within a single study. Furthermore, the accuracy and generalizability of most family studies are limited by various factors, such as modest sample size ranging from 2302 to 17,954 individuals, specific founder populations, different family relationship included (e.g., only siblings or parent-offspring), or the use of self-reported family history not validated by objective measures.

For above reasons, it remains uncertain to what extent cardiometabolic disorders (co-)aggregate in families in the general population, and to what extent the correlation between cardiometabolic disorders and traits can be explained by genetics. Bridging this knowledge gap may help risk stratification and early detection of cardiometabolic disorders. Furthermore, knowledge of shared genetics between disorders and traits may help advance our understanding of pathophysiology. Therefore, we aimed to: (1) quantify the familial (co-)aggregation of various cardiometabolic disorders; (2) estimate the heritability of a wide array of underlying cardiometabolic traits; and (3) estimate genetic correlations between cardiometabolic traits, by using extensive data from Lifelines, a large multi-generational family study representative of the general Dutch population.

Methods

Study population

Lifelines is a multi-disciplinary prospective population-based cohort study examining in a unique three-generation design the health and health-related behaviours of 167,729 persons living in the North of the Netherlands. It employs a broad range of investigative procedures in assessing the biomedical, socio-demographic, behavioural, physical and psychological factors which contribute to the health and disease of the general population, with a special focus on multi-morbidity and complex genetics [25]. The recruitment of participants and their families, and how we define household ID are detailed in the Supplementary Methods (Additional file 1). Briefly, kinship was registered by questionnaires and verified where possible in participants with genetic data (n ~ 80,000). Participants who lived in the same house, which was determined based on registered postal codes, shared the same household ID. In total, there were 30,914 families (of size ≥ 2) of up to four generations and 40,496 singletons (i.e., participants without any relative participating in Lifelines). Physical measurements and collection of biological samples were performed in participants aged 8 years and older [25].

A total of 162,416 participants aged 8 to 93 (152,723 adults and 9693 children) were included in the current cross-sectional study. Data on six cardiometabolic disorders and fifteen cardiometabolic traits were extracted from the database. Cardiometabolic disorders included T2D, MetS, hypertension, hypercholesterolemia, obesity, and CVD. Cardiometabolic traits included markers of glucose regulation (fasting blood glucose, glycated haemoglobin [HbA1c] and skin autofluorescence), blood pressure (systolic and diastolic blood pressure [SBP and DBP]), inflammation (leukocyte count and c-reactive protein [CRP]), obesity (body mass index [BMI] and waist circumference), and lipid levels (total cholesterol, high density lipoprotein cholesterol [HDL], low density lipoprotein cholesterol [LDL], apolipoprotein A, apolipoprotein B, and triglycerides).

Signed informed consent was provided by all participants. Lifelines was conducted according to the principles of the Declaration of Helsinki and following the research code of the University Medical Center Groningen. Meanwhile, the study was approved by the medical ethical committee of the University Medical Center Groningen.

Measurement

Cardiometabolic disorders

Cardiometabolic disorders were defined in adult participants and included T2D, MetS, hypertension, hypercholesterolemia, obesity, and CVD. Individuals with at least one of these six conditions were considered prevalent cases of any cardiometabolic disorder.

Type 2 diabetes

We defined T2D as the combination of self-reported T2D and any supporting data (i.e., use of self-reported glucose-lowering medication, fasting plasma glucose level ≥ 7.0 mmol/L, and/or HbA1c ≥ 6.5%). The following participants were not considered to have T2D: (1) self-reported diabetes but without support by laboratory and medication data; (2) self-reported diabetes but age < 30 years at the time of visit; (3) self-reported diabetes but reported age of onset < 30 years old with self-reported insulin medication.

Metabolic syndrome

According to National cholesterol Education Program Adult Treatment Panel III (NCEP ATP III) definition, participants with three or more of the five following criteria were defined as having MetS: (1) SBP ≥ 130 mmHg, and/or DBP ≥ 85 mmHg, and/or use of antihypertensive medication (Anatomical Therapeutic Chemical [ATC] Classification Codes C02, C03, C07, C08, and/or C09; a key to the ATC codes can be found in the Additional file 2: Table S9); (2) fasting blood glucose ≥ 5.6 mmol/L, and/or use of blood glucose-lowering medication, and/or diagnosis of T2D; (3) HDL cholesterol levels < 1.03 mmol/L in men, and < 1.30 mmol/L in women, and/or use of lipid-lowering medication (ATC codes C10A and/or C10B); (4) triglyceride levels ≥ 1.70 mmol/L and/or use of lipid-lowering medication (ATC codes C10A and/or C10B); and (5) waist circumference ≥ 102 cm in men and ≥ 88 cm in women [26].

Hypertension

Hypertension was defined by SBP ≥ 140 mmHg, and/or DBP ≥ 90 mmHg, and/or the use of antihypertensive medication (ATC codes C02, C03, C07, C08, C09, and/or G04CA03).

Hypercholesterolemia

Total cholesterol ≥ 6.5 mmol/L and/or the use of lipid-lowering medication (ATC codes C10A and/or C10B) were used to define hypercholesterolemia. For participants with self-reported myocardial infarction, total cholesterol ≥ 5.0 mmol/L was used as a cut off to define hypercholesterolemia [27].

Obesity

BMI was calculated as weight (kg)/height squared (m2). BMI ≥ 30 was used to define obesity in adults [28].

Cardiovascular diseases

Four types of CVD were used to define CVD cases at baseline, including myocardial infarction with drugs (platelet aggregation inhibitors/antithrombotic drugs) or ECG abnormalities, self-reported heart failure with drug use (angiotensin converting enzyme inhibitors/angiotensin-II receptor antagonists/aldosterone antagonists) or therapy (pacemaker, ICD implantation or heart transplant), self-reported stroke, and self-reported cardiac surgery (i.e., Coronary Arterial By-pass Graft, Percutaneous Transluminal Coronary Angioplasty, and stent positioning). Participants with at least one of these four CVD were considered prevalent cases. To investigate a broader range of CVDs, we also used an extended CVD definition previously reported by van der Ende et al. [27]. The extended CVD definition additionally included self-reported heart valve problems, self-reported atherosclerosis, self-reported thrombosis, self-reported aneurysm, narrowing carotids, atrial fibrillation with drugs or CHA2DS2-VASc (congestive heart failure, hypertension, age ≥ 75 (doubled), diabetes mellitus, prior stroke or transient ischemic attack (doubled), vascular disease, age 65–74, female) score > 2, and self-reported arrhythmia.

Cardiometabolic traits

Biomarkers were measured from fasting blood samples at the laboratory centre of the University Medical Center Groningen. Skin autofluorescence was measured in adults during baseline visits to quantify the accumulation of Advanced Glycation End products in the skin. Anthropometrics, including body weight, body height, and waist circumference, were performed by a trained research nurse, following the Lifelines protocol [29]. SBP and DBP were measured repetitively 10 times within 10 min, but only the average of the last 3 measurements was used for analysis. Details of the cardiometabolic trait measurements are explained in Additional file 1: Supplementary Methods.

Statistical analysis

To adjust for treatment effects, 15 mmHg and 10 mmHg were added to the SBP and DBP values, respectively, in individuals taking antihypertensive medication. This method has been shown to reduce bias and improve statistical power [30]. Also, for individuals taking lipid-lowering medication, total cholesterol and LDL-C were adjusted by dividing their value by 0.8 and 0.7, respectively [31]. Data was presented depending on the type of variables and its distribution. Continuous variables were presented as mean ± SD when normally distributed and as median and interquartile range when non-normally distributed. Binary cardiometabolic disorders at baseline were described by the number of cases and their prevalences. All analyses were conducted in R version 4.2.2.

Familial aggregation and co-aggregation

Familial aggregation and co-aggregation of cardiometabolic disorders were quantified by the recurrence risk ratio (λR) introduced by Risch [32], and previously applied in Lifelines [13, 33]. The recurrence risk ratio is defined as a ratio between the risk in those with an affected first-degree relative (FDR) and the risk of the total Lifelines population, with λR > 1 indicating positive familial aggregation (i.e. elevated risk in those with positive family history). This λR was estimated using a modified, conditional Cox proportional hazards model applied to cross-sectional data. The modification of this model was performed by applying equal time-to-event for all participants [34]. This model was adjusted for age, age2 (to account for non-linear age effects), sex, and accounted for within-family correlations. We repeated this analysis on spouses as a negative control, as spouses are unlikely to be genetically related and estimates of λR therefore represent the general influence of shared environmental factors and/or assortative mating on cardiometabolic disorders.

As exploratory analyses, we performed familial aggregation analyses of the six cardiometabolic disorders stratified by age, sex, and familial relationship.

Heritability, genetic correlation and phenotypic correlation

Heritability (h2) was defined as the ratio of additive genetic variance to the total phenotypic variance (h2 = VA/VP), using family pedigree as a proxy. The contribution of household (or shared/common environment, c2) was defined as c2 = VC/VP. We estimated h2 and c2 of cardiometabolic traits using residual maximum likelihood-based (REML) variance decomposition under linear mixed models implemented in the ASReml-R package (ASReml, UK) [35]. Significance level of this estimate was derived from likelihood ratio tests, comparing the heritability model to a model in which additive genetic variances was constrained to zero. A detailed description of testing the significance of genetic correlations and phenotypic correlations can be found in the Supplementary Methods (Additional file 1).

Bivariate REML-based variance decomposition was used to estimate the genetic and phenotypic correlation between pairs of cardiometabolic traits. Genetic correlations between traits were estimated as: \(r_G=\frac{{\sigma}_{A1A2}}{\sqrt{{\sigma }_{A1}^{2}{\sigma }_{A2}^{2}}}\), where σA1A2 is the estimated additive genetic covariance between trait x and trait y, and σ2A1 and σ2A2 are the estimated additive genetic variances for traits x and y, respectively. Phenotypic correlations between pairs of traits were obtained from estimated phenotypic covariance and variance as: \(r_P=\frac{{\sigma }_{P1P2}}{\sqrt{{\sigma }_{P1}^{2}{\sigma }_{P2}^{2}}}\), where σP1P2 is the phenotypic covariance between trait x and trait y, and σ2P1 and σ2P2 is the phenotypic variance for traits x and y, respectively.

A rank-based inverse-normal transformation was applied to all variables prior to analysis. In addition, we performed log-transformation on non-normally distributed variables as a sensitivity analysis to further validate our results. All models for heritability, genetic correlations, and phenotypic correlations were adjusted for age, age2 and sex.

Results

Basic characteristics

Of the total of 152,723 included adult participants, 66,695 (43.7%) had at least one cardiometabolic disorder, of which approximately half had multiple (≥ 2) morbidities, shown in Table 1. The prevalence of co-occurring cardiometabolic disorders within individuals is illustrated in Fig. 1. The most common comorbidities in our study population involved a combination of MetS and hypertension with a prevalence of 2.47% among the total participants. Given the definition of MetS, it is possible that MetS overlaps with other cardiometabolic disorders. However, this overlap is incomplete, given that 12.28% of MetS cases only had a single morbidity without any other investigated morbidities (i.e. obesity, hypertension, CVD, T2D, hypercholesterolemia), and less than 1% of MetS cases had all morbidities (details in Additional file 2: Table S1). Meanwhile, 93.4% of participants with T2D had at least one comorbidity, as shown in Additional file 1: Fig. S1). In addition, the prevalence of all types of cardiometabolic disorders was expectedly higher in the older age group (Additional file 1: Fig. S2).

Table 1 Characteristics of the Lifelines participants (N = 162,416)
Fig. 1
figure 1

Upset plot showing the overlap in the prevalence of cardiometabolic disorders. The bars in the lower left panel show the total number of cases of each cardiometabolic disorder. The connected black dots in the matrix panel below represent the combination of multiple cardiometabolic disorders in participants, while a single dot without connection to others indicates a single morbidity in participants. The upper bars, representing the intersection size, show the number of individuals with cardiometabolic disorders and its comorbidities highlighted by the connected black dots in the matrix panel below

Familial aggregation

Positive familial aggregation of cardiometabolic disorders is shown in Fig. 2A; having an FDR or a spouse affected with a certain cardiometabolic disorder associated with a higher risk of the same disorder (ranging from λFDR: 1.23 (95% CI 1.20–1.25) for hypertension, to λFDR: 2.48 (95% CI 2.15–2.86) for T2D; and ranging from λspouse: 1.04 (95% CI 1.00–1.08) for hypercholesterolemia, to λspouse: 1.92 (95% CI 1.83–2.01) for obesity; details in Additional file 2: Table S2). In general, the recurrence risk ratios in individuals with an affected FDR were substantially higher than the recurrence risk ratios in individuals with an affected spouse (λFDR > λspouse). The exception was obesity, for which there was a modestly lower recurrence risk ratio among FDR than spouses (λFDR: 1.85 (95% CI 1.79–1.91) < λspouse: 1.92 (95% CI 1.83–2.01)). In addition to these analyses, we also calculated the familial aggregation of extended CVD and found slightly lower familial aggregation of this extended CVD compared to the original CVD (λFDR: 1.22 (95% CI 1.16–1.28) and λspouse: 1.08 (95% CI 1.01–1.15) for extended CVD, λFDR: 1.53 (95% CI 1.26–1.84) and λspouse: 1.17 (95% CI 0.97–1.41) for the original CVD; details in Additional file 2: Table S3).

Fig. 2
figure 2

Familial aggregation (A) and co-aggregation of cardiometabolic disorders among first-degree relatives (B) and spouses (C), adjusted for age, age2, and sex. Intensity of the colours in figure B and C indicates the magnitude of recurrence risk ratio for cardiometabolic disorders indicated in the X-axis (outcome) by having family member with cardiometabolic disorders indicated in Y-axis (exposure). CVD: any cardiovascular disease defined by at least one of four types of cardiovascular diseases; FDR: first-degree relative; MetS: metabolic syndrome; T2D: type 2 diabetes

Additionally, we performed exploratory stratified familial aggregation analysis to examine possible differential effects of kinship type (affected sibling vs. affected offspring vs. affected parent), sex (men vs. women), and age (age < 40 vs. 40–60 vs. > 60 years). The group with ≥ 1 affected offspring showed the highest recurrence risks, followed by the group with ≥ 1 affected sibling. The group with ≥ 1 affected parent showed the lowest recurrence risks (details in Additional file 2: Table S6). Sex-stratified analyses did not show differences in recurrence risks between men and women (details in Additional file 2: Table S7). Age stratification showed generally similar recurrence risks across age categories. Only recurrence risk of T2D in the group age < 40 seemed higher than in older age categories, although this estimate is relatively imprecise and possibly inflated due to low T2D prevalence in this category (details in Additional file 2: Table S8).

Familial co-aggregation

All pairs of cardiometabolic disorders showed positive familial co-aggregation among FDR (Fig. 2B), as well as spouses (Fig. 2C); either having FDR or a spouse affected with a certain cardiometabolic disorder conferred a higher risk of other cardiometabolic disorders. Significantly elevated recurrence risks of MetS, obesity, and T2D were observed among individuals with a first-degree relative affected by MetS, obesity, or T2D. Although the co-aggregation of the disorders was positive in both FDR and spouses, most of the co-aggregation of the disorders between FDR were higher than spouses (λFDR > λspouse). In spouses, obesity showed a higher co-aggregation with other disorders than in FDR. This higher co-aggregation in spouses than in FDR was also observed in pairs of CVD-MetS and CVD-T2D. In addition, using the extended CVD definition resulted in positive familial co-aggregation, yet with a lower lambda compared to the original CVD definition (from λFDR: 1.06 (95% CI 1.02–1.10) for hypercholesterolemia to λFDR: 1.12 (95% CI 1.08–1.16) for MetS, and from λspouse: 1.05 for hypertension to λspouse: 1.16 for obesity and T2D; details in Additional file 2: Table S3). Unexpectedly, for extended CVD we observed no familial co-aggregation of risk of T2D (λFDR: 1.00 (95% CI 0.92–1.08)).

We estimated λR of cardiometabolic disorders in FDR and spouses affected with one or more cardiometabolic disorders (i.e. ‘any cardiometabolic disorders’). We found a similar co-aggregation pattern, that is, those with a FDR with any cardiometabolic disorders had a higher recurrence risk (1.13 for both obesity and hypertension to 1.17 for T2D; see Additional file 2: Table S2), which were higher than for spouses (λspouse 1.08 for hypertension and T2D).

Heritability

Cardiometabolic traits had moderate levels of heritability (h2CRP: 0.26 to h2HDL: 0.50). Heritability explained much more variance than the shared environment in each cardiometabolic trait (c2Apolipoprotein B: 0.02 to c2 Skin autofluorescence: 0.18; see Fig. 3, details in Additional file 2: Table S4). Together with a substantial heritable component, BMI and waist circumference had the largest estimates for the shared environmental variance component (i.e., c2: 0.17) compared to other cardiometabolic traits. Furthermore, we performed sensitivity analysis using traditionally log-transformed cardiometabolic traits, yielding similar results (Additional file 1: Fig. S3).

Fig. 3
figure 3

Variance component of cardiometabolic biomarkers, adjusted for age, age2, and sex. Y-axis indicates the estimates of variance components of the cardiometabolic traits indicated in the X-axis. h2: variance of cardiometabolic traits due to genetic (Va/Vp), also known as heritability, c2: variance of cardiometabolic traits due to shared (or common) environment (Vc/Vp), e2: variance of cardiometabolic traits due to unique environment (Ve/Vp)

Phenotypic and genetic correlation

Cardiometabolic traits showed a wide range of phenotypic (range from rP HDL-Triglyceride: − 0.47 to rP Apolipoprotein B-LDL and rP total cholesterol-LDL: 0.91) and genetic correlations (range from rG HDL-Triglyceride: − 0.53 to rG Apolipoprotein B-LDL: 0.94), as shown in Fig. 4. Additional information can be found in Table S5 in Additional file 2. When we applied unsupervised hierarchical clustering, these correlations showed strong clustering of traits: (1) HDL cholesterol and apolipoprotein A, and (2) LDL cholesterol, apolipoprotein B, and total cholesterol. A third cluster included the remaining traits. Within this cluster, increased clustering was seen between obesity traits (i.e., BMI and waist circumference) and blood pressure traits (i.e., SBP and DBP). Although included in the third cluster together with other non-lipid traits, phenotypic correlations between skin autofluorescence and other traits were practically absent. In addition, weak-to-moderate genetic correlations were identified between triglycerides with glucose markers, blood pressure, and obesity markers (from rG Triglyceride-HbA1c: 0.16 to rG Triglyceride -waist circumference: 0.30). We also explored the environmental correlations across various cardiometabolic traits (Additional file 1: Fig. S4, details in Additional file 2: Table S5).

Fig. 4
figure 4

Phenotypic (A) and genetic (B) correlations between pairs of wide range of cardiometabolic traits. Dendrogram on the left side represents hierarchical distance between cardiometabolic traits. Clustering of cardiometabolic traits was performed using k-means clustering method, which generated three clusters indicated by the highlighted triangles. Phenotype and genetic correlations were adjusted for age, age2, sex, and shared common environment

Discussion

In this study, we aimed to estimate the genetic and environmental contribution to the co-occurrence of cardiometabolic disorders and traits within families, using objective measurements in a large multi-generational family study. We quantified the familial (co-)aggregation of six cardiometabolic disorders in first-degree relatives and spouses. Individuals with a first-degree relative affected with one of the cardiometabolic disorders had a higher risk of having the same or related disorders. Similarly, individuals with a spouse affected with cardiometabolic disorders had a higher risk of having the same or related disorders, suggesting an effect of shared environmental factors and/or assortative mating. Also, we estimated the heritability of fifteen cardiometabolic traits. These cardiometabolic traits had moderate heritability, indicating a role for genetics underlying the recurrence of cardiometabolic disorders within a family. Finally, we found moderate genetic correlations between cardiometabolic traits, suggesting genetics as an important but not exclusive underlying mechanism of the interrelation between cardiometabolic disorders.

Familial aggregation is evidence for a role of shared genetics and shared environment within a family in the occurrence of complex disorders. Positive familial aggregation of cardiometabolic disorders has been suggested by previous studies [10,11,12,13,14,15,16,17, 21]. Consistent with the literature, our study also found positive familial aggregation between first-degree relatives although of somewhat lower magnitude. For example, in The Framingham Offspring Study, risks of CVD was approximately 2 times higher in middle-aged adults with at least one parent with CVD [10] and 1.5 time if they had at least one sibling with CVD [11]. For T2D, a large Danish study also identified an up to 3.4 times higher risk in first-degree relatives than the general population [36], even higher than our estimate of 2.48 higher risk in first-degree relatives. Other studies also found higher familial aggregation of MetS compared to our recurrence risk estimate λFDR = 1.43 (95% CI 1.39–1.48) in individuals with affected first-degree relatives. A large population-based study in China identified a two to three times higher risk of MetS in younger siblings if their eldest sibling was affected by MetS [15]. Also, in the Tehran Lipid and Glucose Study, the risk of MetS was higher among offspring with affected parents (OR: 2.29–4.53) [16]. A possible explanation for the varying aggregation estimates is the heterogeneity between studies due to different family relationship included, age diversity between studies, and differences in ethnicity, lifestyle, and health behaviors between the Netherlands and other countries. Another possible explanation may arise from the diverse definitions of different disease phenotypes. In our study, we utilized both an extended CVD phenotype and a narrow CVD phenotype, the latter encompassing four CVD types: myocardial infarction, heart failure, stroke, and cardiac surgery. We observed a slightly lower familial aggregation in the extended CVD phenotype. This difference could be attributed to varying levels of heterogeneity in the extended CVD phenotype compared to the narrow CVD phenotype, with the latter demonstrating greater homogeneity. Despite these differences, the evidence converges on a major role of shared genetics in determining cardiometabolic risk. In exploratory analysis, we observed a higher recurrence risk in individuals with affected siblings and affected offspring when compared to individuals with affected parents, although possibly this is driven by age differences between offspring and parents. However age-stratified exploratory analysis showed that the recurrence risk estimates were stable across both age and sex.

We considered recurrence risk estimates in spouses a negative control to those in family: given that spouses are unlikely to be genetically related, estimates of spousal recurrence reflect the effects of shared environment and/or assortative mating on cardiometabolic risk. A Danish study observed around 1.5 times higher risk of T2D among individuals with a spouse affected with T2D [36], and in previous work in both Japanese and Dutch, we corroborated higher spousal cardiometabolic risk, for T2D (OR: 1.20 vs 1.59), hypertension (OR: 1.34 vs 1.45), and MetS (OR: 1.77 vs 1.77) in Japanese vs Dutch [37]. Although using the same cohort study, the estimates for spousal concordance in the Dutch Lifelines Population [37] are slightly higher than our current familial aggregation estimates, which may be attributed to differences in the statistical approaches used to estimate the risk of disorders among spouses. The spousal concordance was related to concordance in lifestyle factors, such as physical activity, smoking, and alcohol drinking, indicative of potential cohabitation effects, and assortative mating [37,38,39].

We found that cardiometabolic traits had moderate genetic extent in which the heritability estimates are consistent with those in the literature although some variation can occur due to differences in age, ethnicity, study design [40], and type of measurement used [41]. For example, slightly different estimates of heritability were found in a family study in a 1564 Chinese individuals from 494 families [20], reporting heritability of fasting glucose (h2: 0.17), waist circumference (h2: 0.26), SBP (h2: 0.24), DBP (h2: 0.17), triglycerides (h2: 0.41), HDL-C (h2: 0.49), LDL-C (h2: 0.47), total cholesterol (h2: 0.46), CRP (h2: 0.38), and BMI (h2 = 0.38). Also, a Dutch twin study presented moderate to high heritability for a range of cardiometabolic traits from 0.47 for insulin level to 0.78 for BMI [42], with estimates typically being higher than family studies including the present study.

Genetic similarity between first-degree relatives is likely to contribute to familial co-aggregation of related cardiometabolic disorders [21,22,23,24]. A previous study found that parental history of one or more CVD (i.e., myocardial infarction, stroke, and angina) at a younger age < 55 year in the father and < 65 year in the mother significantly increased the risk of MetS in women, with ORs from 1.62 to 1.84 [21]. Another family study identified increased risk of having multiple cardiometabolic disorders in relation to parental history of diabetes (OR: 1.54, 95% CI 1.01–2.33) and parental history of hypertension (OR 1.42, 95% CI 1.06–1.91). The risk was even higher when both parents were affected with hypertension or diabetes, suggesting an additive genetic effect on the risk of cardiometabolic disease co-occurrence [24]. Similarly, a population-based study in US found a higher risk of co-occurring cardiometabolic disorders when individuals had a family history of diabetes or hypertension, and only a slightly increased risk with family history of obesity [23]. These studies mostly found evidence for co-aggregation between CVD, T2D, and hypertension. Compared to these studies, our estimates of co-aggregation between obesity, T2D, and MetS were larger, while hypertension, CVD, and hypercholesterolemia showed only modest familial co-aggregation. A possible reason for this difference is that previous studies used mostly self-reported family history without actual validation with objective laboratory measures in family members.

Cardiometabolic disorders likely share pathophysiological mechanisms such as inflammation and insulin resistance [4,5,6]. The high insulin resistance in obesity and diabetes is thought to induce inflammation, causing vascular damage and endothelial dysfunction. Such vascular damage and dysfunction lead to increased production of vasoconstrictors, and subsequently to an increase in vascular resistance, a major contributor to CVD and hypertension [5, 6]. Consistent with this, we found significant phenotypic and genetic correlations between blood pressure, obesity traits, inflammatory markers, and fasting glucose, although these correlations were of modest strength. In the present study, we also found modest genetic correlations between HDL cholesterol and triglyceride with glucose markers, blood pressure, and obesity, which are the traits used for the definition of MetS. The accumulation of Advanced Glycation End products has been considered a potential cross-link between diabetes and cardiovascular events, by increasing inflammation and causing endothelial dysfunction [43, 44]. We found little evidence of this in the present study: we observed only weak correlations between skin autofluorescence and other cardiometabolic traits, although their genetic correlations are slightly higher than the phenotypic correlation. Furthermore, the association of skin autofluorescence, CVD, and T2D in the previous studies in this population were independent of glucose markers [43, 45].

Strengths of this study are that it is the largest family and comprehensive study of cardiometabolic outcomes to date, investigating six interrelated cardiometabolic disorders and fifteen intermediate cardiometabolic traits. It was conducted in a large-scale population-based, multi-generational cohort, representative of the general Dutch population [46]. We used a combination of objective laboratory measurements, questionnaire data, and medication data, resulting in precise outcome definitions and thus precise estimates of recurrence risk and heritability. Our study provides insights into the genetic and/or environmental mechanisms that underly of co-existing cardiometabolic disorders within individuals and families. Our findings highlight the role of shared genetics and environmental factors on the risk of cardiometabolic disorders and suggests overlapping genetic structure between disorders. Furthermore, our estimates of recurrence risk may inform clinicians and health services in diagnosis, patient communication, and potential screening efforts based on family history. However, several study limitations need to be addressed as well. Firstly, only data were available on family members participating in Lifelines. Missing data on non-participating family members may have caused underestimation of recurrence risk and heritability. Secondly, although largely representative of the Dutch population Lifelines predominantly consists of Dutch participants; generalization to other ancestries and maybe even to other European populations is therefore uncertain.

Conclusion

To conclude, in our large multi-generational family study, cardiometabolic disorders show positive (co-)aggregation within families and to a lesser extent between spouses. We found moderate heritability for a wide variety of intermediate cardiometabolic traits and moderate genetic correlations between traits. These results suggest that genetic factors are an important but moderate contributor to the co-occurrence of cardiometabolic traits. We find evidence for a strong contribution of shared environmental factors, especially for obesity. To further elucidate potential mechanisms for co-aggregation, future studies may focus on identifying the shared genetic factors, the specific shared environmental risk factors, and potential gene-environment interaction.

Availability of data and materials

Lifelines data may be obtained from a third party and are not publicly available. Researchers can apply to use the Lifelines data used in this study. More information about how to request Lifelines data and the conditions of use can be found on their website (https://www.lifelines.nl/researcher/how-to-apply).

Abbreviations

BMI:

Body mass index

CVD:

Cardiovascular diseases

CRP:

C-reactive protein

c2 :

Contribution of household

DBP:

Diastolic blood pressure

FDR:

First-degree relatives

HbA1c:

Glycated haemoglobin

HDL:

High density lipoprotein

h2 :

Heritability

LDL:

Low density lipoprotein

MetS:

Metabolic syndrome

REML:

Residual maximum likelihood

rG :

Genetic correlation

rP :

Phenotypic correlation

SBP:

Systolic blood pressure

T2D:

Type 2 diabetes

λR :

Recurrence risk ratios

References

  1. Roth GA, Mensah GA, Johnson CO, Addolorato G, Ammirati E, Baddour LM, Barengo NC, Beaton AZ, Benjamin EJ, Benziger CP, et al. Global burden of cardiovascular diseases and risk factors, 1990–2019: update from the GBD 2019 study. J Am Coll Cardiol. 2020;76(25):2982–3021.

    PubMed  PubMed Central  Google Scholar 

  2. International Diabetes Federation. IDF Diabetes Atlas, 10th edn. Brussels, Belgium: 2021 [https://www.diabetesatlas.org]

  3. Grundy SM. Pre-diabetes, metabolic syndrome, and cardiovascular risk. J Am Coll Cardiol. 2012;59(7):635–43.

    CAS  PubMed  Google Scholar 

  4. Reaven GM. Banting lecture 1988. Role of insulin resistance in human disease. Diabetes. 1988;37(12):1595–607.

    CAS  PubMed  Google Scholar 

  5. Aroor AR, McKarns S, DeMarco VG, Jia G, Sowers JR. Maladaptive immune and inflammatory pathways lead to cardiovascular insulin resistance. Metabolism. 2013;62(11):1543–52.

    CAS  PubMed  Google Scholar 

  6. Hill MA, Yang Y, Zhang L, Sun Z, Jia G, Parrish AR, Sowers JR. Insulin resistance, cardiovascular stiffening and cardiovascular disease. Metabolism. 2021;119: 154766.

    CAS  PubMed  Google Scholar 

  7. Yang G-R, Dye TD, Li D. Association between diabetes, metabolic syndrome and heart attack in US adults: a cross-sectional analysis using the Behavioral Risk Factor Surveillance System 2015. BMJ Open. 2019;9(9): e022990.

    PubMed  PubMed Central  Google Scholar 

  8. Aboonabi A, Meyer RR, Singh I. The association between metabolic syndrome components and the development of atherosclerosis. J Hum Hypertens. 2019;33(12):844–55.

    PubMed  Google Scholar 

  9. Pollack LM, Wang M, Leung MYM, Colditz G, Herrick C, Chang S-H. Obesity-related multimorbidity and risk of cardiovascular disease in the middle-aged population in the United States. Prev Med. 2020;139: 106225.

    PubMed  PubMed Central  Google Scholar 

  10. Lloyd-Jones DM, Nam BH, D’Agostino RB Sr, Levy D, Murabito JM, Wang TJ, Wilson PW, O’Donnell CJ. Parental cardiovascular disease as a risk factor for cardiovascular disease in middle-aged adults: a prospective study of parents and offspring. JAMA. 2004;291(18):2204–11.

    CAS  PubMed  Google Scholar 

  11. Murabito JM, Pencina MJ, Nam BH, D’Agostino RB Sr, Wang TJ, Lloyd-Jones D, Wilson PW, O’Donnell CJ. Sibling cardiovascular disease as a risk factor for cardiovascular disease in middle-aged adults. JAMA. 2005;294(24):3117–23.

    CAS  PubMed  Google Scholar 

  12. Hemminki K, Li X, Sundquist K, Sundquist J. Familial risks for type 2 diabetes in Sweden. Diabetes Care. 2010;33(2):293–7.

    PubMed  Google Scholar 

  13. Wang R, Snieder H, Hartman CA. Familial co-aggregation and shared heritability between depression, anxiety, obesity and substance use. Transl Psychiatry. 2022;12(1):108.

    PubMed  PubMed Central  Google Scholar 

  14. Fuentes RM, Notkola IL, Shemeikka S, Tuomilehto J, Nissinen A. Familial aggregation of body mass index: a population-based family study in eastern Finland. Horm Metab Res. 2002;34(7):406–10.

    CAS  PubMed  Google Scholar 

  15. Feng Y, Zang T, Xu X, Xu X. Familial aggregation of metabolic syndrome and its components in a large Chinese population. Obesity (Silver Spring). 2008;16(1):125–9.

    CAS  PubMed  Google Scholar 

  16. Azizi F, Farahani ZK, Ghanbarian A, Sheikholeslami F, Mirmiran P, Momenan AA, Asl SZ, Hadaegh F, Eskandari F. Familial aggregation of the metabolic syndrome: Tehran Lipid and Glucose Study. Ann Nutr Metab. 2009;54(3):189–96.

    CAS  PubMed  Google Scholar 

  17. Pollex RL, Hegele RA. Genetic determinants of the metabolic syndrome. Nat Clin Pract Cardiovasc Med. 2006;3(9):482–9.

    CAS  PubMed  Google Scholar 

  18. Poveda A, Chen Y, Brandstrom A, Engberg E, Hallmans G, Johansson I, Renstrom F, Kurbasic A, Franks PW. The heritable basis of gene-environment interactions in cardiometabolic traits. Diabetologia. 2017;60(3):442–52.

    CAS  PubMed  Google Scholar 

  19. Xia C, Amador C, Huffman J, Trochet H, Campbell A, Porteous D, Generation S, Hastie ND, Hayward C, Vitart V, et al. Pedigree- and SNP-associated genetics and recent environment are the major contributors to anthropometric and cardiometabolic trait variation. PLoS Genet. 2016;12(2): e1005804.

    PubMed  PubMed Central  Google Scholar 

  20. Lin CC, Peyser PA, Kardia SL, Li CI, Liu CS, Chu JS, Lin WY, Li TC. Heritability of cardiovascular risk factors in a Chinese population–Taichung Community Health Study and Family Cohort. Atherosclerosis. 2014;235(2):488–95.

    CAS  PubMed  Google Scholar 

  21. Dallongeville J, Grupposo MC, Cottel D, Ferrieres J, Arveiler D, Bingham A, Ruidavets JB, Haas B, Ducimetiere P, Amouyel P. Association between the metabolic syndrome and parental history of premature cardiovascular disease. Eur Heart J. 2006;27(6):722–8.

    PubMed  Google Scholar 

  22. Berntorp K, Lindgarde F. Familial aggregation of type 2 diabetes mellitus as an etiological factor in hypertension. Diabetes Res Clin Pract. 1985;1(6):307–13.

    PubMed  Google Scholar 

  23. Hunt KJ, Heiss G, Sholinsky PD, Province MA. Familial history of metabolic disorders and the multiple metabolic syndrome: the NHLBI family heart study. Genet Epidemiol. 2000;19(4):395–409.

    CAS  PubMed  Google Scholar 

  24. Liese AD, Mayer-Davis EJ, Tyroler HA, Davis CE, Keil U, Schmidt MI, Brancati FL, Heiss G. Familial components of the multiple metabolic syndrome: the ARIC study. Diabetologia. 1997;40(8):963–70.

    CAS  PubMed  Google Scholar 

  25. Sijtsma A, Rienks J, van der Harst P, Navis G, Rosmalen JGM, Dotinga A. Cohort profile update: lifelines, a three-generation cohort study and biobank. Int J Epidemiol. 2021;51: e295.

    PubMed Central  Google Scholar 

  26. Alberti KGMM, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, Fruchart J-C, James WPT, Loria CM, Smith SC. Harmonizing the metabolic syndrome. Circulation. 2009;120(16):1640–5.

    CAS  PubMed  Google Scholar 

  27. van der Ende MY, Hartman MH, Hagemeijer Y, Meems LM, de Vries HS, Stolk RP, de Boer RA, Sijtsma A, van der Meer P, Rienstra M, et al. The lifelines cohort study: prevalence and treatment of cardiovascular disease and risk factors. Int J Cardiol. 2017;228:495–500.

    PubMed  Google Scholar 

  28. Weir CB, Jan A. BMI classification percentile and cut off points. Treasure Island: StatPearls; 2022.

    Google Scholar 

  29. Scholtens S, Smidt N, Swertz MA, Bakker SJ, Dotinga A, Vonk JM, van Dijk F, van Zon SK, Wijmenga C, Wolffenbuttel BH, et al. Cohort profile: lifelines, a three-generation cohort study and biobank. Int J Epidemiol. 2015;44(4):1172–80.

    PubMed  Google Scholar 

  30. Tobin MD, Sheehan NA, Scurrah KJ, Burton PR. Adjusting for treatment effects in studies of quantitative traits: antihypertensive therapy and systolic blood pressure. Stat Med. 2005;24(19):2911–35.

    PubMed  Google Scholar 

  31. Baigent C, Keech A, Kearney PM, Blackwell L, Buck G, Pollicino C, Kirby A, Sourjina T, Peto R, Collins R, et al. Efficacy and safety of cholesterol-lowering treatment: prospective meta-analysis of data from 90,056 participants in 14 randomised trials of statins. Lancet. 2005;366(9493):1267–78.

    CAS  PubMed  Google Scholar 

  32. Risch N. Linkage strategies for genetically complex traits. I. Multilocus models. Am J Hum Genet. 1990;46(2):222–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Zhang J, Thio CHL, Gansevoort RT, Snieder H. Familial aggregation of CKD and heritability of kidney biomarkers in the general population: the lifelines cohort study. Am J Kidney Dis. 2021;77(6):869–78.

    CAS  PubMed  Google Scholar 

  34. Breslow N. Covariance analysis of censored survival data. Biometrics. 1974;30(1):89–99.

    CAS  PubMed  Google Scholar 

  35. Butler DG, Cullis BR, Gilmour AR, Gogel BG, Thompson R. ASReml-R reference manual version 4. Hemel Hempstead: VSN International Ltd; 2017.

    Google Scholar 

  36. Aasbjerg K, Norgaard CH, Vestergaard N, Sogaard P, Kober L, Weeke P, Gislason G, Torp-Pedersen C. Risk of diabetes among related and unrelated family members. Diabetes Res Clin Pract. 2020;160: 107997.

    CAS  PubMed  Google Scholar 

  37. Nakaya N, Xie T, Scheerder B, Tsuchiya N, Narita A, Nakamura T, Metoki H, Obara T, Ishikuro M, Hozawa A, et al. Spousal similarities in cardiometabolic risk factors: a cross-sectional comparison between Dutch and Japanese data from two large biobank studies. Atherosclerosis. 2021;334:85–92.

    CAS  PubMed  Google Scholar 

  38. Roberts ME, Banse R, Ebbeler C, Ferketich AK. Spousal concordance in the use of alternative tobacco products: a multi-country investigation. Drug Alcohol Depend. 2017;171:16–9.

    PubMed  Google Scholar 

  39. Cobb LK, Godino JG, Selvin E, Kucharska-Newton A, Coresh J, Koton S. Spousal influence on physical activity in middle-aged and older adults: the ARIC study. Am J Epidemiol. 2016;183(5):444–51.

    PubMed  Google Scholar 

  40. Asefa NG, Neustaeter A, Jansonius NM, Snieder H. Heritability of glaucoma and glaucoma-related endophenotypes: systematic review and meta-analysis. Surv Ophthalmol. 2019;64(6):835–51.

    PubMed  Google Scholar 

  41. Macgregor S, Cornes BK, Martin NG, Visscher PM. Bias, precision and heritability of self-reported and clinically measured height in Australian twins. Hum Genet. 2006;120(4):571–80.

    PubMed  Google Scholar 

  42. van Dongen J, Willemsen G, Chen WM, de Geus EJ, Boomsma DI. Heritability of metabolic syndrome traits in a large population-based sample. J Lipid Res. 2013;54(10):2914–23.

    PubMed  PubMed Central  Google Scholar 

  43. Boersma HE, van Waateringe RP, van der Klauw MM, Graaff R, Paterson AD, Smit AJ, Wolffenbuttel BHR. Skin autofluorescence predicts new cardiovascular disease and mortality in people with type 2 diabetes. BMC Endocr Disord. 2021;21(1):14.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Nedosugova LV, Markina YV, Bochkareva LA, Kuzina IA, Petunina NA, Yudina IY, Kirichenko TV. Inflammatory mechanisms of diabetes and its vascular complications. Biomedicines. 2022;10(5):1168.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. van Waateringe RP, Fokkens BT, Slagter SN, van der Klauw MM, van Vliet-Ostaptchouk JV, Graaff R, Paterson AD, Smit AJ, Lutgers HL, Wolffenbuttel BHR. Skin autofluorescence predicts incident type 2 diabetes, cardiovascular disease and mortality in the general population. Diabetologia. 2019;62(2):269–80.

    PubMed  Google Scholar 

  46. Klijs B, Scholtens S, Mandemakers JJ, Snieder H, Stolk RP, Smidt N. Representativeness of the lifelines cohort study. PLoS ONE. 2015;10(9): e0137203.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We acknowledge the services of the Lifelines Cohort Study, the contributing research centres delivering data to Lifelines, and all the study participants. We also appreciate the technical support from Arthur Gilmour for estimating heritability, phenotypic, and genetic correlations using ASReml-R.

Funding

The Lifelines initiative has been made possible by subsidy from the Dutch Ministry of Health, Welfare and Sport, the Dutch Ministry of Economic Affairs, the University Medical Center Groningen (UMCG), Groningen University and the Provinces in the North of the Netherlands (Drenthe, Friesland, Groningen). ZC and RDT acknowledge the support from the China Scholarship Council (no. 202008440344) and the Indonesia Endowment Funds for Education (LPDP), respectively.

Author information

Authors and Affiliations

Authors

Contributions

RDT, ZC, IMN, CHLT, HS contributed to the study conception and design. AA, RW, CAH contributed to the data preparation and analysis. RDT and ZC performed the data analysis and drafted the manuscript equally. All authors participated in revising it critically for important intellectual content.

Corresponding authors

Correspondence to Chris H. L. Thio or Harold Snieder.

Ethics declarations

Ethics approval and consent to participate

Lifelines was conducted according to the principles of the Declaration of Helsinki and following the research code of the University Medical Center Groningen. Meanwhile, the study was approved by the medical ethical committee of the University Medical Center Groningen. All the participants agreed to take part in the study and provided their written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no potential conflict of interest.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1:

Supplementary Methods. Figure S1. Morbidities in cardiometabolic disorders. Figure S2 Prevalence of cardiometabolic disorders per age category. Figure S3 Heritability estimates of log-transformed cardiometabolic traits. Figure S4 Correlations between cardiometabolic traits due to shared environment.

Additional file 2:

Table S1. Overlapping cardiometabolic disorders in adults participants (N=152,723). Table S2. Familial aggregation and coaggregation of cardiometabolic disorders (N=152,723). Table S3. Familial aggregation and coaggregation of cardiometabolic disorders using extended definition of cardiovascular diseases (N=152,723). Table S4. Variance component estimates (N=162,416). Table S5. Phenotypic, genetic and environmental correlation of cardiometabolic traits (N=162,416) Table S6. Familial aggregation of cardiometabolic disorders in first-degree relatives stratified by familial kinship (N=152,723). Table S7. Familial aggregation of cardiometabolic disorders in first-degree relatives stratified by sex (N=152,723). Table S8. Familial aggregation of cardiometabolic disorders in first-degree relatives stratified by age (N=152,723). Table S9. Explanation of different ATC codes used in the study.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Triatin, R.D., Chen, Z., Ani, A. et al. Familial co-aggregation and shared genetics of cardiometabolic disorders and traits: data from the multi-generational Lifelines Cohort Study. Cardiovasc Diabetol 22, 282 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12933-023-02017-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12933-023-02017-w

Keywords