- Research
- Open access
- Published:
Causal associations between insulin and Lp(a) levels in Caucasian population: a Mendelian randomization study
Cardiovascular Diabetology volume 23, Article number: 316 (2024)
Abstract
Background
Numerous observational studies have demonstrated that circulating lipoprotein(a) [Lp(a)] might be inversely related to the risk of type 2 diabetes (T2D). However, recent Mendelian randomization (MR) studies do not consistently support this association. The results of in vitro research suggest that high insulin concentrations can suppress Lp(a) levels by affecting apolipoprotein(a) [apo(a)] synthesis. This study aimed to identify the relationship between genetically predicted insulin concentrations and Lp(a) levels, which may partly explain the associations between low Lp(a) levels and increased risk of T2D.
Methods
Independent genetic variants strongly associated with fasting insulin levels were identified from meta-analyses of genome-wide association studies in European populations (GWASs) (N = 151,013). Summary level data for Lp(a) in the population of European ancestry were acquired from a GWAS in the UK Biobank (N = 361,194). The inverse-variance weighted (IVW) method approach was applied to perform two-sample summary-level MR. Robust methods for sensitivity analysis were utilized, such as MR‒Egger, the weighted median (WME) method, MR pleiotropy residual sum and outlier (MR-PRESSO), leave-one-out analysis, and MR Steiger.
Results
Genetically predicted fasting insulin levels were negatively associated with Lp(a) levels (β = − 0.15, SE = 0.05, P = 0.003). The sensitivity analysis revealed that WME (β = − 0.26, SE = 0.07, P = 0.0002), but not MR‒Egger (β = − 0.22, SE = 0.13, P = 0.11), supported a causal relationship between genetically predisposed insulin levels and Lp(a).
Conclusion
Our MR study provides robust evidence supporting the association between genetically predicted increased insulin concentrations and decreased concentrations of Lp(a). These findings suggest that hyperinsulinaemia, which typically accompanies T2D, can partially explain the inverse relationship between low Lp(a) concentrations and an increased risk of T2D.
Background
Lipoprotein(a) [Lp(a)] is classified as a highly atherogenic lipoprotein that consists of a lipid-rich domain, apolipoprotein B100, and apolipoprotein(a) [apo(a), not to be confused with apolipoprotein A] encoded by the LPA gene on chromosome 6 (q25.3–q26) [1]. The serum Lp(a) concentration ranges from 70 to 90%, as determined by LPA gene variants [2]. The variability in Lp(a) concentrations is primarily due to “copy number variations” in the number of repeats of Kringle-IV 2 (KIV-2) domains and, to a lesser extent, by single-nucleotide polymorphisms (SNPs) within and in the proximity of this gene, as well as loss-of-function mutations [3]. In addition, sex, diet, physical activity, hormonal balance, and kidney and liver diseases may influence Lp(a) levels [4, 5]. It is worth noting that in the interpretation of available results, there is a significant drawback related to the strongly skewed distribution of Lp(a) concentrations in a population and differences between patients from different ethnic groups [6].
The physiological role of Lp(a) in humans has not yet been fully clarified. However, it has been suggested that it is responsible for the transport of oxidized phospholipids and affects the processes of carcinogenesis, inflammation, coagulation, and wound healing [7]. The impact of elevated Lp(a) concentrations on the development of cardiovascular disease (CVD) has been confirmed by epidemiological studies, including genome-wide association (GWASs) and Mendelian randomization (MR) studies [8,9,10,11].
Several studies have suggested that there is an inverse association between Lp(a) concentration and the risk of type 2 diabetes (T2D), whereby people with a low Lp(a) concentration have a greater risk of T2D [12,13,14]. However, it is unclear whether this inverse association between Lp(a) concentration and T2D is causal. MR analysis is a method that can explore causal associations between exposures and outcomes arising from observational epidemiological studies [15]. Nevertheless, evidence from MR studies has not confirmed the causal relationship between low Lp(a) concentrations and a greater risk of T2D [16,17,18]. Indeed, some data indicate that not the Lp(a) concentration but the high number of KIV-2 repeats may be associated with an increased risk of T2D [19, 20]. The difficulty in interpreting the results of MR studies can be partly explained by differences in analysis strategies and the heterogeneity of instrumental variables, as individual SNPs may have different functional significance (only some are loss-of-function mutations) and differ in their relationships with KIV-2 repeats [21].
Previously published in vitro studies have shown that insulin can regulate Lp(a) synthesis in the liver [22, 23]. Specifically, increased insulin levels may reduce Lp(a) synthesis. This relationship may partially explain the association of Lp(a) with T2D observed in epidemiological studies because increasing tissue insulin resistance, leading to T2D, is accompanied by hyperinsulinaemia. Unfortunately, there is insufficient in vivo data to properly validate these results.
Several MR analyses focusing on the impact of insulin on lipid traits were conducted in the past, however, the impact of insulin on Lp(a) was not extensively studied [24, 25]. We have hypothesized that the reason for the apparent correlation of low Lp(a) with T2D without clear confirmation from MR studies may be explained by reverse causation where the increased fasting insulin affects Lp(a) levels.
This study applied a two-sample MR approach using summary-level GWAS data to examine this hypothesis to identify whether genetically predicted insulin levels influence Lp(a) concentrations in the Caucasian population.
Materials and methods
Study design
In this study, we used a two-sample MR approach utilising publicly available GWAS summary statistics to investigate the causal relationship between insulin and Lp(a). This approach eliminates the need for additional ethical approval or informed consent [26]. This study adhered to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines, which are mainly designed for MR studies (STROBE-MR) (Supplementary Table 1, Additional file 1) [27].
Exposure data
IVs (instrument variables) for BMI-adjusted fasting insulin were derived from meta-analysis of aggregated GWASs for N = 151,013 individuals of European ancestry without diabetes (a subset of the total number of 281,416 individuals included in the full study) [28]. Access to the GWAS summary statistics was obtained through the OpenGWAS project developed at the MRC Integrative Epidemiology Unit (IEU) [29]. The GWAS ID corresponding to insulin was “ebi-a-GCST90002238”.
Outcome data
Summary-level data for Lp(a), Apolipoprotein B (ApoB) and LDL-cholesterol (LDL-C) from UK Biobank, which were inverse-rank normalized by the Neale laboratory (N = 361,194), were obtained through the OpenGwas project [29]. The GWAS IDs corresponding to Lp(a) ApoB and LDL-C are “ukb-d-30790_irnt”, “ukb-d-30640_irnt” and “ukb-d-30780_irnt” respectively. There were no overlapping populations between the GWASs exposure and outcome.
Selection of instrumental variables
To be valid IVs in MR, SNPs must meet three core assumptions: (1) relevance, (2) independence or exchangeability and (3) exclusion restriction. Whether a SNP meets those criteria cannot be confirmed but can be falsified. Therefore, SNPs were considered appropriate IVs if they were associated with insulin levels of genome-wide significance (p-value < 5.00e−08) and not in linkage disequilibrium (LD) with each other (r2 < 0.001) within a clumping distance of 10,000 kb. For each IV, the variance (R2) was calculated using the following formula: R2 = 2×beta2×EAF×(1 − EAF), where EAF is the effect allele frequency and beta is the association estimate. Additionally, F-statistics were calculated based on the formula F=(beta/se)2, where ‘beta’ and ‘se’ refer to the genetic association of SNPs with the exposure and its standard error, respectively [15]. SNPs with F-statistics < 10 were considered weak IVs and excluded from the analysis [30]. Full lists of all SNPs at different stages of selection are presented in Supplementary Tables 2–3 (Additional file 2). To prevent potential pleiotropy, associations of the selected SNPs with confounding factors influencing Lp(a) concentration were evaluated using PhenoScanner V2 [31]. For further analysis, two sets of IVs were used: all identified SNPs except pleiotropic variants (conservative analysis) and all identified SNPs with pleiotropic variants (liberal analysis).
Statistical analyses
The MR approach was used to estimate the association between the genetically predicted concentration of fasting insulin and Lp(a). In the first step, variant harmonization between datasets was conducted to confirm that the association between SNPs and the exposure and between SNPs and the outcome reflected the same alleles. If the selected SNPs were unavailable in the outcome dataset, they were replaced with proxy SNPs with an LD of r2 > 0.8 or excluded from further MR analysis. Furthermore, outlier pleiotropic variants were excluded via MR-PRESSO [32]. After that, two-sample MR analyses were performed with the inverse-variance weighted (IVW) method with multiplicative random effects as the primary analysis for evaluating the causal effect estimates in our study [33].
In the next step, sensitivity analyses were performed using the following methods: MR‒Egger, weighted median (WME), and MR pleiotropy residual sum and outlier (MR-PRESSO) [32, 34, 35]. Additionally, the presence of horizontal pleiotropy was assessed by calculating the MR‒Egger intercept deviation [36]. The Cochran Q statistic was also calculated to measure heterogeneity between variant-specific causal estimates [37]. Furthermore, to examine the possibility that reverse causality exists in these studies, the MR Steiger test was performed and a bidirectional two-sample MR analysis between genetically predicted Lp(a) and insulin was performed using the same approach described earlier [38]. Finally, leave-one-out analyses were performed to determine the possible effect of SNPs on the causal estimates. As an additional sensitivity analysis, we used variants employed as insulin IVs in Buchmann and colleagues’ paper to validate our findings on a different IVs set [39]. Moreover, due to the complex structure of Lp(a), we conducted additional MR analyses on the primary components of Lp(a), such as LDL-C and ApoB, to measure the mediating effects of insulin on each component separately. In these analyses, we used only the IVs set from the liberal analysis. A p-value (P) of < 0.05 in this study was considered to indicate statistical significance for all tests. All these analyses were conducted with R software (version 4.1.1) [40] using “TwoSampleMR” (version 0.5.6) [41] and the “MR-PRESSO” package (version 1.0) [32].
Results
According to the criteria set, 38 IVs associated with fasting insulin levels were obtained (Supplementary Table 3, Additional file 2). After searching in PhenoScanner, 22 SNPs were removed from the conservative analysis because they are related to known or likely confounders (kidney diseases, lipid parameters, psychosocial and lifestyle outcomes) for Lp(a) level (Supplementary Table 4, Additional file 2, records with exclusion reasons for SNPs are presented in bold) to remove pleiotropic SNPs that may affect Lp(a) levels through mechanisms other than affecting apo(a) expression via insulin [42]. No outliers were detected by the MR-PRESSO. The genetic variants used in conservative, liberal, and additional sensitivity analyses are listed in Supplementary Tables 5–7 (Additional file 3).
As shown in Table 1; Fig. 1, the primary conservative analysis revealed higher genetically predicted fasting insulin levels to be associated with decreased Lp(a) levels (β = − 0.15, SE = 0.05, P = 0.003). Moreover, these findings were supported by the results of the liberal analysis (Supplementary Tables 8 and Supplementary Fig. 1, Additional file 3), in which genetically determined fasting insulin levels also influenced Lp(a) levels (βIVM = −0.13, SEIVM = 0.04, PIVM = 0.0002; βWME = −0.14, SEWME = 0.05, PWME = 0.005). In the sensitivity analysis, the results of the WME (β = −0.26, SE = 0.07; P = 0.0002) method supported this association. The MR‒Egger results were not significant (β = −0.22, SE = 0.13, P = 0.11) but yielded a consistent estimate. This is possibly due to the reduced power of MR‒Egger compared to IVM due to the additional intercept term. Based on Cochran’s Q test (PIVM = 0.53, PMR-Egger = 0.49) and MR-PRESSO analysis (P = 0.60), heterogeneity or outliers were not detected between the exposure and outcome genetic variants. Additionally, the MR‒Egger intercept test did not show evidence of horizontal pleiotropy (Pintercept = 0.54). The leave-one-out analysis indicated no significant change in the MR results after removing one SNP at a time (Fig. 2).
Genetic associations between fasting insulin and Lp(a). Each genetic variant included in the analysis is represented as a point + 95% CI. Location on the horizontal axis represents the correlation of the variant with exposure (plasma fasting insulin, inverse variance normal-transformed values); location on the vertical axis represents the correlation of the variant with the outcome (Lp(a), inverse-rank normalized). Lines represent estimates of different MR methods
Results of leave-one-out analysis. Horizontal axis—inverse-variance weighted-mean estimate with 95% confidence intervals of genetic associations between insulin and Lp(a). Each dot and confidence interval represents an inverse-variance weighted mean obtained when a variant listed on the left side is removed from the analysis.
Furthermore, no evidence of reverse causality was found between the exposure and outcome, as determined by the Steiger test (P = 5.67e−78, Table 1) and by conducting bidirectional MR (Table 2; Fig. 3).
Genetic associations between Lp(a) levels and fasting insulin.Each genetic variant included in the analysis is represented as a point + 95% CI. Localization on the horizontal axis represents the correlation of the variant with exposure (Lp(a), inverse-rank normalized). Localization on the vertical axis represents the correlation of the variant with the outcome (plasma fasting insulin, inverse variance normal transformed values). Lines represent estimates of different MR methods.
Additional analysis, which included IVs from Buchmann and colleagues’ paper [39], supported an association between genetically elevated fasting insulin and decreased Lp(a) by IVM (β = −0.19, SE = 0.07, P = 0.003) and WME (β = − 0.17, SE = 0.08, P = 0.04), but the results were not statistically significant by MR‒Egger (β = − 0.33, SE = 0.16, P = 0.07) method (Table 3; Fig. 4).
Genetic associations between fasting insulin levels and Lp(a) in additional sensitivity analysis. Each genetic variant included in the analysis is represented as a point + 95% CI. Localization on the horizontal axis represents the correlation of the variant with exposure (plasma fasting insulin, inverse variance normal transformed values). Localization on the vertical axis represents the correlation of the variant with the outcome (Lp(a), inverse-rank normalized). Lines represent estimates of different MR methods.
We have also analysed the impact of fasting insulin on ApoB and LDL-C. Among 38 IV SNPs identified for fasting insulin, MR-PRESSO analysis has identified some outlier SNPs: 8 for LDL-C and 10 for ApoB. Full lists of SNPs excluded from the analysis for each parameter are presented in Supplementary Tables 9–10, Additional file 3. In the analysis of secondary Lp(a) components we have found that ApoB levels increased with fasting insulin levels (IVW: β = 0.24, SE = 0.004, P = 5.78e−08). The results of the WME (β = 0.27, SE = 0.05; P = 4.17e−07) and MR-Egger (β = 0.30, SE = 0.13; P = 2.53e−02) method supported this association (Supplementary Tables 11 and Supplementary Figs. 2–3, Additional file 3). LDL-C was found to increase due to an increase in genetic fasting insulin concentrations but only when the estimate was computed with IVW method (β = 0.09, SE = 0.04, P = 0.03) (Supplementary Tables 12 and Supplementary Figs. 4–5, Additional file 3). Notably, all analyses demonstrated significant heterogeneity in the IVs (Cochran Q test P < 0.05).
Discussion
Our analysis revealed that an increase in the fasting insulin concentration causes a decrease in the Lp(a) concentration. It is known that hyperinsulinaemia is usually associated with the development of T2D. First, the emerging insulin resistance of peripheral tissues leads to excessive insulin synthesis, and both factors lead to prediabetes and then to T2D. In patients with T2D, insulin levels decrease as pancreatic beta cells become depleted [43]. Many studies have shown an inverse association between lower Lp(a) and insulin resistance [44, 45].
In vitro studies on the effect of insulin on Lp(a) synthesis have shown that high insulin concentrations reduce apo(a) synthesis. Such a relationship was indicated in a study conducted by Neele et al. based on cell culture of monkey hepatocytes. The study was motivated by the fact that elevated Lp(a) levels decreased after insulin therapy in patients with type 1 diabetes (T1D). The experiment showed that high insulin concentrations inhibit apo(a) synthesis at the posttranscriptional level, with the peak of inhibition occurring after 72 h of incubation of cells with insulin [23]. In another study published by Suzuki H et al. utilizing the HepG2 cell line, the researchers detected reduced apo(a) promoter activity due to the influence of insulin [22]. Importantly, MR tests chronic rather than short term exposures and this differentiates it from conventional in vitro studies. Unfortunately, there is a paucity of in vitro studies exploring the relationship between insulin and Lp(a) and further research in this area is necessary.
We would like to emphasize that, to the best of our knowledge, our MR study is the second analysis of this type published to date aiming to verify whether there are actual indications of regulation of Lp(a) concentration by insulin. Previous MR analyses typically concentrated on the causal impact of Lp(a) on T2D and did not produce fully consistent results [16,17,18].
The first previously published MR study was based on a meta-analysis of cross-sectional data for patients from three independent cohorts, i.e., Berlin Ageing Study II (BASE-II; N = 2012), LIFE-Adult (N = 3281) and LIFE-Heart (N = 2816), and used 2 or 12 IVs. That previous study comprised a one-sample MR analysis in which the genetic IVs risk factor associations and genetic IVs outcome associations were derived from the same sample, and individual-level data were used to derive MR estimates. In contrast to the results of our investigation, that study yielded no significant evidence of insulin-dependent regulation of Lp(a) concentration, although some analyses showed an insignificant tendency for insulin to reduce Lp(a) levels [39]. In our analysis, we performed two-sample MR on aggregated data, i.e., MR analyses in which the genetic IVs-risk factor associations and genetic IVs-outcome associations were (ideally) obtained from different (nonoverlapping) samples. Two-sample MR reduces the severity of “weak instrument bias”, “sample overlap”, and the appearance of “Winner’s curse” compared to one-sample MR, and the bias direction tends to approach null. The reliability of our results is primarily due to the type of analysis employed, the significantly larger group of subjects used, the larger number of IVs involved, and the more restrictive cut-off point for genome-wide significance applied (p-value < 5.00e−08). In addition, to validate our findings on a different IV set, we utilized variants used as insulin IVs in Buchmann and colleagues’ paper and performed MR, and we also obtained a result that partially supported our findings.
Interestingly, studies on acute exposure to exogenous insulin differ in their findings when compared to our MR analysis. An experimental study published by Sidhu et al. verified the relationship between serum Lp(a) concentration and insulin metabolism in healthy men without obesity undergoing an intravenous glucose tolerance test. Their study revealed a significant inverse relationship between the insulin level and Lp(a) concentration in the first, but not second phase of the plasma insulin response [46]. Another study conducted by Sutherland et al. evaluated the effect of changes in Lp(a) concentrations in T2D patients and healthy controls receiving intravenous insulin infusion. Plasma Lp(a) levels showed a small but significant increase in mean concentrations in T2D patients, but mean Lp(a) concentrations did not change significantly in healthy subjects [47]. In addition, two studies conducted by Riemens et al. are worth mentioning. In the first study, no Lp(a) concentration changes were observed when assessing the effects of 3-h hyperinsulinaemia induced by hyperglycaemia in healthy subjects [48]. In the second study, Riemes et al. verified the impact of 24-h intravenous insulin infusion on the Lp(a) concentration in plasma, and no significant changes in Lp(a) levels [49]. The results described above contrast with our observations, but it should be noted that MR analyses take into account long-term minor differences in exposure and not short-term, high-intensity interventions, as in the case of the results described by the researchers. The half-life of Lp(a) in human plasma is approximately 3 days [50]. Therefore, even several hours of insulin infusion in those studies may not have been long enough to detect significant differences in Lp(a) synthesis. This seems particularly likely considering the previously discussed results published by Neele et al., in which monkey hepatocytes showed a peak inhibition of Lp(a) synthesis under the influence of insulin only after 72 h [23]. It should be noted that, in comparison with the results of other research, our results are based on the MR, which is a method that enables estimating causal relationships in the observational study using genetic variants as instrumental variables. Our results are based on a large sample and the exclusion of diabetes as a confounding variable because these individuals are a heterogeneous group and may show different degrees of diabetes development.
As mentioned, MR analysis examines the effects of long-term changes in exposure over the lifespan rather than short-term effects. Studies determining the effects of chronic insulin therapy on Lp(a) levels yielded results more consistent with those of our MR analysis. In the first such study, Wolffenbuttel et al. examined the effect of improving blood glucose control with insulin therapy on Lp(a) and other lipoproteins in patients with T2D with poor disease control. It was concluded that insulin therapy does not alter Lp(a) levels in T2D patients after 6 months of insulin treatment [51]. Another study by Kuusi et al. aimed to verify the effect of insulin treatment on Lp(a) levels in T2D patients. Lp(a) concentration increased in patients with Lp(a) concentrations less than 300 mg/L but not in patients with Lp(a) concentrations greater than 300 mg/L. Interestingly, metabolic control was not related to the described effect [52]. Another study examined the effects of improved glycaemic control following 3 months of insulin therapy on Lp(a) concentration in patients with poorly controlled T2D. Lp(a) levels did not change, regardless of glycaemic control [53]. Although not a study of insulin administration but also an interventional study, a trial on 140 obese children (mean age 12.5 ± 1.6 years), demonstrated that a 3-week diet has significantly reduced both serum insulin and serum Lp(a) levels [54]. Once again, we would like to emphasize that in our analyses, the exposure dataset excluded patients with diabetes due to potential confounding due to disease treatment and progression. Our model reflects the dependence of Lp(a) on insulin in a physiological state and not during the progression and treatment of diabetes, when beta cells may eventually fail, insulin be of exogenous origin and insulin-sensitivity of tissues be modified. Additionally, in the case of individuals using insulin therapy, there may also be differences between insulin concentrations compared to the profile resulting from endogenous insulin secretion. These situations may impact the change in Lp(a) secretion. Therefore, in our opinion, caution should be taken when comparing the results of our analysis, which is based on the effect of insulin on Lp(a) concentration in non-diabetic individuals, with the results based on the assessment of insulin therapy in diabetes.
The available results of cross-sectional studies aimed at establishing the relationship between insulin and Lp(a) concentrations also show divergent conclusions, but they are more in line with the results of our analysis. In a study conducted by Zamora-González et al. in the Mexican population, the fasting insulin concentration was found to be directly and independently inversely related to the Lp(a) level in patients of both sexes [55]. Another important study also concerning the Mexican population was conducted by Posadas-Romero et al., including included lean and obese subjects with normoglycaemia and normal blood pressure. Lp(a) concentrations were similar in women with obesity and high and normal concentrations of insulin (19.9 mg/dl vs. 18.6 mg/dl), but men with obesity and hyperinsulinaemia had significantly lower Lp(a) levels than did men with obesity and normo-insulinaemia (7.9 mg/dl vs. 29.4 mg/dl). These observations were not found in the female group [56]. On the other hand, in Jichi Medical Cohort Study, which included 1121 men and 1480 women, Inoue et al. showed a weak, statistically insignificant inverse correlation between insulin and Lp(a) in both men and women [57]. Another study by Habib et al. assessed the relationship between basal insulin concentration and Lp(a) concentration in a group of patients with and without T2D. The analysed groups were additionally divided into two subgroups based on normal (< 10 µU/mL) or elevated fasting insulin concentrations (≥ 10 µU/mL). The Lp(a) concentration was significantly lower in individuals with T2D with increased fasting insulin concentrations than in individuals with normal fasting insulin concentrations. Regression analysis confirmed that the Lp(a) concentration correlated inversely with the insulin concentration in individuals with T2D (P < 0.05) [58]. Another interesting study on Lp(a) and insulin concentrations in different types of diabetes was conducted by Heller et al. The study group consisted of patients with T1D and patients with T2D treated with diet, diet in combination with oral hypoglycaemic drugs, or requiring insulin therapy. The results of the analyses indicated that the concentration of Lp(a) in serum was not significantly related to the duration of diabetes in the individual or to the degree of glycaemic control. T2D patients receiving insulin therapy were observed to have increased serum Lp(a) concentrations [59]. A study conducted by Vaverkova et al. has found that in dyslipidemic patients there is a negative correlation between serum Lp(a) and insulin [60]. Finally, in a study of 7633 Chinese men, patients with insulin resistance were found to exhibit both higher insulin levels and lower Lp(a) levels [44]. On the other hand, the positive correlation between fasting insulin and Lp(a) was noted in a sample of 131 normal-weight Mexican children with an average age of 7.7 ± 0.8 years. It was reported that children with hyperinsulinemia had significantly increased Lp(a) levels [61]. It may be possible that factors affecting Lp(a) and insulin relationship are different in pre-pubertal children. In summary, several observational studies may suggest that endogenous fasting insulin is indeed negatively associated with Lp(a).
Overall, the results of studies published thus far indicate that endogenous insulin has a slight tendency to reduce Lp(a) concentrations. In contrast, acute administration of exogenous insulin may have the opposite effect or have no significant impact on Lp(a) concentration. In this respect, the previously discussed in vitro study by Neele et al. may be difficult to classify as either acute or chronic exposure, but in that case, the exposure of liver cells to insulin was also more long-lasting and constant than in the case of bolus insulin therapy or an infusion of this hormone lasting a mere few hours and the study explored apo(a) expression rather than simple Lp(a) levels. As already mentioned, our MR study examined long-term changes in exposure to endogenous insulin. Hence, the results of the observational studies seem consistent with MR results, suggesting a negative albeit weak effect of insulin concentration on plasma Lp(a) concentrations.
Lp(a) levels appear to be related to sex and sex hormones. The literature is not fully consistent in terms of the correlation between gender and Lp(a) levels. An increased serum Lp(a) levels were found in diabetic females compared to males and non-diabetic female, regardless of menopausal status [62], in white women but not in black women [63], while in a different study no significant differences between males and females in regards to Lp(a) category were detected [64]. Hormone replacement therapy was found to lower Lp(a) concentration in post-menopausal women [65, 66]. Interestingly, a decrease in Lp(a) concentrations was also found in healthy males after testosterone administration [67, 68] but not in post-menopausal women [69]. On the other hand, no correlation between Lp(a) and sex hormone concentrations was found in a study of Finnish and American men [70], while in the study of Italian males, only a negative correlation between Lp(a) and DHEA-S was found [71]. In one study, Lp(a) was reported to be positively correlated with estradiol concentrations but only in women with self-reported heart disease [72].
In our study, we have also analysed the impact of fasting insulin levels on components of Lp(a) different than apo(a), namely ApoB and LDL-cholesterol, which was chosen because Lp(a) lipid core mostly resembles that of LDL particles. Lipoprotein(a) assembly mechanism does not appear to be clearly elucidated, but some have suggested ApoB to be a limiting factor based on studies in subjects with abetalipoproteinemia [73, 74]. However, we have found ApoB levels to increase with insulin levels which is the exact opposite of the expected effect if fasting insulin decreased Lp(a) through reducing ApoB availability. Interestingly, other MR analyses have found ApoB to be either unaffected by fasting insulin levels or to be increased by BMI-adjusted insulin levels, like ones used in our study [24, 25]. This finding contrasts with the results of in vitro and in vivo studies [75,76,77,78][ It may be possible that in our study the increase in ApoB due to elevated fasting insulin reflects metabolic dysfunction and insulin resistance rather than simple direct action of insulin. The results regarding the influence of fasting insulin on LDL-C were not fully consistent. Some have found that fasting insulin only affects some LDL fractions in terms of their composition [24], which may explain this phenomenon.
Our analysis concluded that increased serum insulin concentration may contribute to a decrease in Lp(a) concentration. Nevertheless, it is important to acknowledge some limitations of this study. This study relied on publicly accessible summary-level data, and the available data precluded us from performing additional subgroup analyses to explore the influence of nongenetic factors affecting Lp(a) levels, such as sex, hormones, or diet. Moreover, the exposure data excluded patients with diabetes, which raises the possibility of selection bias and collider bias [79], although insulin levels in patients with diabetes can also be affected by treatment and disease progression [80]. Additionally, the study was limited to individuals of European descent and did not distinguish between males and females, raising questions about the generalizability of our findings to populations of different ethnicities and genders. Finally, insulin and glycemia-related traits are highly complex and interconnected. Further analyses are required to determine if there exist indirect effects of insulin on Lp(a) levels due to other mechanisms related to insulin resistance rather than direct action of insulin.
There are two clinical implications of this study. First, it is clear that the relationship between Lp(a) and diabetes is not trivial and that Lp(a) is likely not a risk factor for the development of diabetes independent of preexisting hyperinsulinaemia and insulin resistance. Second, the observational relationship between low Lp(a) and diabetes risk may not translate to adverse effects of therapies that reduce Lp(a) levels due to causal relationships. Although therapies aimed at reducing insulin resistance, insulinaemia and glycaemia may increase Lp(a), it is almost certain that their cardiometabolic benefits outweigh the increased cardiovascular risk caused by an increase in Lp(a) as demonstrated by the fact that good glycemic control improves patient survival [81].
Conclusion
The results of our MR study support evidence that elevated insulin concentrations may be causally related to decreased serum Lp(a) levels. Our findings suggest that hyperinsulinaemia, which typically accompanies T2D [82], can partially explain the inverse relationship between low Lp(a) concentrations and an increased risk of T2D. Further studies are needed to validate and establish the exact mechanism behind this relationship.
Availability of data and materials
The statistical code used in this work is available upon reasonable request to the corresponding author. Access to the GWAS summary statistics was obtained through the OpenGWAS project developed at the MRC Integrative Epidemiology Unit (IEU) [26]. The GWAS ID corresponding to insulin was “ebi-a-GCST90002238”. The GWAS ID corresponding to Lp(a) was “ukb-d-30790_irnt”. No original raw data were generated in this study.
Abbreviations
- apo(a):
-
Apolipoprotein(a)
- apoB:
-
Apolipoprotein B100
- CVD:
-
Cardiovascular diseases
- EAF:
-
Effect allele frequency
- GWASs:
-
Genome-wide association studies
- IVs:
-
Instrument variables
- IVW:
-
Inverse-variance weighted
- KIV-2:
-
Kringle-IV 2
- LD:
-
Linkage disequilibrium
- Lp(a):
-
Lipoprotein(a)
- MR:
-
Mendelian randomization
- T1D:
-
Type 1 diabetes
- T2D:
-
Type 2 diabetes
- WME:
-
Weighted median
- MR-PRESSO:
-
MR pleiotropy residual sum and outlier
- SNP:
-
Single-nucleotide polymorphism
- STROBE:
-
Strengthening the reporting of observational studies in epidemiology
References
Kronenberg F. Lipoprotein(a). Handb Exp Pharmacol. 2022;270:201–32.
Reyes-Soffer G, Ginsberg HN, Berglund L, Duell PB, Heffron SP, Kamstrup PR, et al. Lipoprotein(a): a genetically determined, causal, and prevalent risk factor for atherosclerotic cardiovascular disease: a scientific statement from the American Heart Association. Arterioscler Thromb Vasc Biol. 2022;42:e48-60.
Coassin S, Kronenberg F. Lipoprotein(a) beyond the kringle IV repeat polymorphism: the complexity of genetic variation in the LPA gene. Atherosclerosis. 2022;349:17–35.
Enkhmaa B, Berglund L. Non-genetic influences on lipoprotein(a) concentrations. Atherosclerosis. 2022;349:53–62.
Ruscica M, Sirtori CR, Corsini A, Watts GF, Sahebkar A. Lipoprotein(a): knowns, unknowns and uncertainties. Pharmacol Res. 2021;173:105812.
Enkhmaa B, Anuurad E, Berglund L. Lipoprotein (a): impact by ethnicity and environmental and medical conditions. J Lipid Res. 2016;57:1111–25.
Orsó E, Schmitz G. Lipoprotein(a) and its role in inflammation, atherosclerosis and malignancies. Clin Res Cardiol Suppl. 2017;12:31–7.
Duarte Lau F, Giugliano RP. Lipoprotein(a) and its significance in cardiovascular disease: a review. JAMA Cardiol. 2022;7:760–9.
Kronenberg F. Lipoprotein(a): from causality to treatment. Curr Atheroscler Rep. 2024;26:75–82.
Welsh P, Al Zabiby A, Byrne H, Benbow HR, Itani T, Farries G, et al. Elevated lipoprotein(a) increases risk of subsequent major adverse cardiovascular events (MACE) and coronary revascularisation in incident ASCVD patients: a cohort study from the UK Biobank. Atherosclerosis. 2024;389:117437.
Banach M. Lipoprotein(a): the enemy that we still don’t know how to defeat. Eur Heart J Open. 2023;3:oead080.
Paige E, Masconi KL, Tsimikas S, Kronenberg F, Santer P, Weger S, et al. Lipoprotein(a) and incident type–2 diabetes: results from the prospective Bruneck study and a meta-analysis of published literature. Cardiovasc Diabetol. 2017;16:38.
Muhanhali D, Zhai T, Cai Z, Ling Y. Lipoprotein(a) concentration is associated with risk of type 2 diabetes and cardiovascular events in a Chinese population with very high cardiovascular risk. Endocrine. 2020;69:63–72.
Mora S, Kamstrup PR, Rifai N, Nordestgaard BG, Buring JE, Ridker PM. Lipoprotein(a) and risk of type 2 diabetes. Clin Chem. 2010;56:1252–60.
Burgess S, Thompson SG. Mendelian randomization: methods for causal inference using genetic variants. 2nd ed. Boca Raton London New York: CRC Press; 2021.
Ye Z, Haycock PC, Gurdasani D, Pomilla C, Boekholdt SM, Tsimikas S, et al. The association between circulating lipoprotein(a) and type 2 diabetes: is it causal? Diabetes. 2014;63:332–42.
Kamstrup PR, Nordestgaard BG. Lipoprotein(a) concentrations, isoform size, and risk of type 2 diabetes: a mendelian randomisation study. Lancet Diabetes Endocrinol. 2013;1:220–7.
Goławski M, Lejawa M, Banach M, Jóźwiak J, Gierlotka M, Osadnik T. Association between Lp(a) and T2D: a Mendelian randomization study. Arch Med Sci. 2024;20:1002–5.
Zhai TY, Ling Y, Gao X. LPA kringle IV type 2 is associated with type 2 diabetes in a Chinese population with very high cardiovascular risk. J Lipid Res. 2018;59:884–91.
Tolbus A, Mortensen MB, Nielsen SF, Kamstrup PR, Bojesen SE, Nordestgaard BG, Kringle IV. Type 2, not low Lipoprotein(a), as a cause of diabetes: a novel genetic approach using SNPs associated selectively with lipoprotein(a) concentrations or with Kringle IV type 2 repeats. Clin Chem. 2017;63:1866–76.
Lamina C, Ward NC. Lipoprotein (a) and diabetes mellitus. Atherosclerosis. 2022;349:63–71.
Suzuki H, Matsushima T, Yamashita K. Transcriptional regulation of apolipoprotein(a) synthesis by insulin. Atherosclerosis. 1997;134:150–1.
Neele DM, de Wit EC, Princen HM. Insulin suppresses apolipoprotein(a) synthesis by primary cultures of cynomolgus monkey hepatocytes. Diabetologia. 1999;42:41–4.
Wong THT, Mo JMY, Zhou M, Zhao JV, Schooling CM, He B, et al. A two-sample Mendelian randomization study explores metabolic profiling of different glycemic traits. Commun Biol. 2024;7:293.
Zhao JV, Luo S, Schooling CM. Sex-specific Mendelian randomization study of genetically predicted insulin and cardiovascular events in the UK Biobank. Commun Biol. 2019;2:332.
Lawlor DA. Commentary two-sample Mendelian randomization: opportunities and challenges. Int J Epidemiol. 2016;45:908–15.
Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomization: the STROBE-MR statement. JAMA. 2021;326:1614–21.
Chen J, Spracklen CN, Marenne G, Varshney A, Corbin LJ, Luan J, et al. The trans-ancestral genomic architecture of glycemic traits. Nat Genet. 2021;53:840–60.
Elsworth B, Lyon M, Alexander T, Liu Y, Matthews P, Hallett J, et al. The MRC IEU OpenGWAS data infrastructure. Genetics. 2020. https://doi.org/10.1101/2020.08.10.244293.
Burgess S, Thompson SG, CRP CHD Genetics Collaboration. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol. 2011;40:755–64.
Phenoscanner [Internet]. Phenoscanner. http://www.phenoscanner.medschl.cam.ac.uk/about/
Verbanck M, Chen C-Y, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet. 2018;50:693–8.
Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. 2017;36:1783–802. https://doi.org/10.1002/sim.7221.
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40:304–14. https://doi.org/10.1002/gepi.21965.
Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28:30–42.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.
Greco MFD, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in mendelian randomisation studies with summary data and a continuous outcome. Stat Med. 2015;34:2926–40.
Hemani G, Tilling K, Davey Smith G, Correction. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13:e1007149.
Buchmann N, Scholz M, Lill CM, Burkhardt R, Eckardt R, Norman K, et al. Association between lipoprotein(a) level and type 2 diabetes: no evidence for a causal role of lipoprotein(a) and insulin. Acta Diabetol. 2017;54:1031–8.
R Core Team. R: A language and environment for statistical computing. [Internet]. Vienna, Austria: R Foundation for Statistical Computing. 2021. https://www.R-project.org/
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408.
Enkhmaa B, Berglund L. Non-genetic influences on lipoprotein(a) concentrations. Atherosclerosis. 2022;349:53–62.
Thomas DD, Corkey BE, Istfan NW, Apovian CM. Hyperinsulinemia: an early indicator of metabolic dysfunction. J Endocr Soc. 2019;3:1727–47.
Wang H, Fan J-L. Association between lipoprotein(a) and insulin resistance in Chinese adults: results from the China health and nutrition survey. Front Endocrinol (Lausanne). 2023;14:1234140.
Meroni M, Longo M, Lombardi R, Paolini E, Macchi C, Corsini A, et al. Low Lipoprotein(a) levels predict hepatic fibrosis in patients with nonalcoholic fatty liver disease. Hepatol Commun. 2022;6:535–49.
Sidhu M, Crook D, Godsland IF, Walton C, Wynn V, Oliver MF. Inverse relationship between serum Lp(a) levels and first-phase insulin secretion. Diabetes. 1992;41:1341–5.
Sutherland WH, Lewis-Barned NJ, Pratt MC, Tillman HC, Walker RJ. The response of plasma lipoprotein (a) concentration to acute hyperinsulinemia in patients with non-insulin dependent diabetes and in healthy subjects. Res Commun Chem Pathol Pharmacol. 1994;83:77–86.
Riemens SC, Ligtenberg JJ, Dullaart RP. Hyperglycemia-induced hyperinsulinemia acutely lowers plasma apolipoprotein B but not lipoprotein (a) in man. Clin Chim Acta. 1997;261:149–58.
Riemens SC, Dullaart RP. Twenty-four hours of insulin infusion does not lower plasma lipoprotein(a) in healthy men. Scand J Clin Lab Invest. 2000;60:119–23.
Krempler F, Kostner GM, Bolzano K, Sandhofer F. Turnover of Lipoprotein (a) in man. J Clin Investig. 1980;65:1483–90.
Wolffenbuttel BH, Leurs PB, Sels JP, Rondas-Colbers GJ, Menheere PP, Nieuwenhuijzen Kruseman AC. Improved blood glucose control by insulin therapy in type 2 diabetic patients has no effect on lipoprotein(a) levels. Diabet Med. 1993;10:427–30.
Kuusi T, Yki-Järvinen H, Kauppinen-Mäkelin R, Jauhiainen M, Ehnholm C, Kauppila M, et al. Effect of insulin treatment on serum lipoprotein(a) in non-insulin-dependent diabetes. Eur J Clin Invest. 1995;25:194–200.
Caixas A, Pérez A, Qrdóñez-Llanos J, Bonet R, Rigla M, Castellví A, et al. Lack of change of lipoprotein(a) levels by the optimization of glycemic control with insulin therapy in NIDDM patients. Diabetes Care. 1997;20:1459–61.
Brandstätter A, Lingenhel A, Zwiauer K, Strobl W, Kronenberg F. Decrease of Lp(a) during weight reduction in obese children is modified by the apo(a) kringle-IV copy number variation. Int J Obes. 2009;33:1136–42.
Zamora-González J, Yamamoto-Kimura L, Lerman-Garber I, Cardoso-Saldaña G, Fajardo-Gutierrez A, Posadas-Romero C. Clustering of metabolic disorders and hyperinsulinemia in Mexico City. Int J Obes Relat Metab Disord. 1996;20:311–8.
Posadas-Romero C, Hernández-Ono A, Zamora-González J, Cardoso-Saldaña G, Yamamoto-Kimura L, Brito-Zurita OR. Differential effects of obesity with and without hyperinsulinemia on plasma lipoprotein(a) concentrations in men. Metabolism. 2001;50:178–83.
Inoue K, Nago N, Matsuo H, Goto T, Miyamoto T, Saegusa T, et al. Serum insulin and lipoprotein(a) concentrations. The Jichi Medical School Cohort Study. Diabetes Care. 1997;20:1242–7.
Habib SS, Aslam M, Shah SFA, Naveed AK. Lipoprotein (a) is associated with basal insulin levels in patients with type 2 diabetes Mellitus. Arq Bras Cardiol. 2009;93:28–33.
Heller FR, Jamart J, Honore P, Derue G, Novik V, Galanti L, et al. Serum lipoprotein(a) in patients with diabetes mellitus. Diabetes Care. 1993;16:819–23.
Vaverková H, Karásek D, Halenka M, Cibíčková L, Kubíčková V. Inverse association of Lipoprotein (a) with markers of insulin resistance in dyslipidemic subjects. Physiol Res. 2017;66:113–20.
Rodríguez-Moran M, Gamboa-Gómez CI, Preza-Rodríguez L, Guerrero-Romero F. Lipoprotein(a) and hyperinsulinemia in healthy normal-weight, prepubertal Mexican Children. Endocr Res. 2021;46:87–91.
Nakhjavani M, Morteza A, Esteghamati A, Khalilzadeh O, Zandieh A, Safari R. Serum Lipoprotein(a) levels are greater in female than male patients with type–2 diabetes. Lipids. 2011;46:349–56. https://doi.org/10.1007/s11745-010-3513-1.
Bovet P, Rickenbach M, Wietlisbach V, Riesen W, Shamlaye C, Darioli R, et al. Comparison of serum Lipoprotein(a) distribution and its correlates among black and white populations. Int J Epidemiol. 1994;23:20–7. https://doi.org/10.1093/ije/23.1.20.
Bigazzi F, Minichilli F, Sbrana F, Pino BD, Corsini A, Watts GF, et al. Gender difference in lipoprotein(a) concentration as a predictor of coronary revascularization in patients with known coronary artery disease. Biochimica et Biophysica Acta (BBA) Mol Cell Biol Lipids. 2021;1866:158869.
Mijatovic V, Kenemans P, Netelenbos JC, Peters-Muller ERA, Van Kamp GJ, Voetberg GA, et al. Oral 17β-estradiol continuously combined with dydrogesterone lowers serum Lipoprotein(a) concentrations in healthy postmenopausal women. J Clin Endocrinol Metab. 1997;82:3543–7.
Kim CJ, Jang HC, Cho DH, Min YK. Effects of hormone replacement therapy on lipoprotein(a) and lipids in postmenopausal women. Arterioscler Thromb. 1994;14:275–81.
Zmuda JM, Thompson PD, Dickenson R, Bausserman LL. Testosterone decreases lipoprotein(a) in men. Am J Cardiol. 1996;77:1244–7.
Anderson RA, Wallace EM, Wu FCW. Effect of testosterone enanthate on serum lipoproteins in man. Contraception. 1995;52:115–9.
Zang H, Carlström K, Arner P, Hirschberg AL. Effects of treatment with testosterone alone or in combination with estrogen on insulin sensitivity in postmenopausal women. Fertility and Sterility. 2006;86:136–44.
Haffner SM, Mykkanen L, Gruber KK, Rainwater DL, Laakso M. Lack of association between sex hormones and Lp(a) concentrations in American and Finnish men. Arterioscler Thromb. 1994;14:19–24.
Denti L, Pasolini G, Ablondi F, Valenti G. Correlation between plasma Lipoprotein Lp(a) and sex hormone concentrations: a cross-sectional study in healthy males. Horm Metab Res. 1994;26:602–8.
Sowers M, Crawford SL, Cauley JA, Stein E. Association of Lipoprotein(a), insulin resistance, and reproductive hormones in a multiethnic cohort of pre- and perimenopausal women (The SWAN Study). Am J Cardiol. 2003;92:533–7.
Dieplinger H. The Kidney is the Heart of the organs: its role in Lp(a) physiology and pathophysiology. In: Kostner K, Kostner GM, Toth PP, editors. Lipoprotein(a). Cham: Springer International Publishing; 2023. p. 207–30.
Menzel HJ, Dieplinger H, Lackner C, Hoppichler F, Lloyd JK, Muller DR, et al. Abetalipoproteinemia with an ApoB–100-lipoprotein(a) glycoprotein complex in plasma. Indication for an assembly defect. J Biol Chem. 1990;265:981–6.
Haas ME, Attie AD, Biddinger SB. The regulation of ApoB metabolism by insulin. Trends Endocrinol Metab. 2013;24:391–7.
Salhanick AI, Schwartz SI, Amatruda JM. Insulin inhibits apolipoprotein B secretion in isolated human hepatocytes. Metabolism. 1991;40:275–9.
Mondola P, Patti L, Santangelo F, Coraggio S. Changes in serum apolipoprotein B concentration after oral glucose or intramuscular insulin. Acta Diabet Lat. 1984;21:235–40. https://doi.org/10.1007/BF02642896.
Karimian Pour N, Adeli K. Insulin silences apolipoprotein B. mRNA translation by inducing intracellular traffic into cytoplasmic RNA granules. Biochemistry. 2011;50:6942–50. https://doi.org/10.1021/bi200711v.
Schneider EB. Collider bias in economic history research. Explor Econ History. 2020;78:101356.
Fonseca VA. Defining and characterizing the progression of type 2 diabetes. Diabetes Care. 2009;32:S151-6.
Gaster B, Hirsch IB. The effects of improved glycemic control on complications in type 2 diabetes. Arch Intern Med. 1998;158:134.
Weyer C, Hanson RL, Tataranni PA, Bogardus C, Pratley RE. A high fasting plasma insulin concentration predicts type 2 diabetes independent of insulin resistance: evidence for a pathogenic role of relative hyperinsulinemia. Diabetes. 2000;49:2094–101.
Acknowledgements
Not applicable.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
MG and MLe were responsible for the conception and design. MLe performed the statistical analysis. MLe, MG and MF were involved in drafting the manuscript. TO, MB, FP, MR, NP, and MLi were involved in reviewing and revising the manuscript. All the authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study only used published or publicly available data; no original data were collected. Ethical approval by the corresponding institutional review board or ethics committee
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary material
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/.
About this article
Cite this article
Lejawa, M., Goławski, M., Fronczek, M. et al. Causal associations between insulin and Lp(a) levels in Caucasian population: a Mendelian randomization study. Cardiovasc Diabetol 23, 316 (2024). https://doi.org/10.1186/s12933-024-02389-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12933-024-02389-7