Association of Lipoprotein (a) variants with risk of cardiovascular disease: a Mendelian randomization study

Abstract

Background

There is a well-documented empirical relationship between lipoprotein (a) [Lp(a)] and cardiovascular disease (CVD); however, causal evidence, especially from the Chinese population, is lacking. Therefore, this study aims to estimate the causal association between variants in genes affecting Lp(a) concentrations and CVD in people of Han Chinese ethnicity.

Methods

Two-sample Mendelian randomization analysis was used to assess the causal effect of Lp(a) concentrations on the risk of CVD. Summary statistics for Lp(a) variants were obtained from 1256 individuals in the Cohort Study on Chronic Disease of Communities Natural Population in Beijing, Tianjin and Hebei. Data on associations between single-nucleotide polymorphisms (SNPs) and CVD were obtained from recently published genome-wide association studies.

Results

Thirteen SNPs associated with Lp(a) levels in the Han Chinese population were used as instrumental variables. Genetically elevated Lp(a) was inversely associated with the risk of atrial fibrillation [odds ratio (OR), 0.94; 95% confidence interval (95%CI), 0.901–0.987; P = 0.012)], the risk of arrhythmia (OR, 0.96; 95%CI, 0.941–0.990; P = 0.005), the left ventricular mass index (OR, 0.97; 95%CI, 0.949–1.000; P = 0.048), and the left ventricular internal dimension in diastole (OR, 0.97; 95%CI, 0.950–0.997; P = 0.028) according to the inverse-variance weighted method. No significant association was observed for congestive heart failure (OR, 0.99; 95% CI, 0.950–1.038; P = 0.766), ischemic stroke (OR, 1.01; 95%CI, 0.981–1.046; P = 0.422), and left ventricular internal dimension in systole (OR, 0.98; 95%CI, 0.960–1.009; P = 0.214).

Conclusions

This study provided evidence that genetically elevated Lp(a) was inversely associated with atrial fibrillation, arrhythmia, the left ventricular mass index and the left ventricular internal dimension in diastole, but not with congestive heart failure, ischemic stroke, and the left ventricular internal dimension in systole in the Han Chinese population. Further research is needed to identify the mechanism underlying these results and determine whether genetically elevated Lp(a) increases the risk of coronary heart disease or other CVD subtypes.

Background

Cardiovascular disease (CVD) is a class of diseases that include coronary artery diseases, stroke, heart failure, hypertensive heart disease, rheumatic heart disease, cardiomyopathy, abnormal heart rhythms, congenital heart disease, valvular heart disease, carditis, aortic aneurysms, peripheral artery disease, thromboembolic disease, and venous thrombosis [1]. CVD is the leading cause of death in China and worldwide [2, 3], accounting for two out of every five deaths in China, according to the China Cardiovascular Disease Report 2019 [4]. It has been estimated that the number of patients with CVD in China is over 330 million, and the prevalence of CVD continues to rise [4]. Recently, a study by Ooi EM and colleagues showed that the concentration of lipoprotein (a) [Lp(a)] is independently associated with angiographic extent and severity of coronary artery disease [5], a study by Rallidis LS and colleagues showed that high levels of Lp(a) are continuously and independently associated with an increased risk of acute coronary syndrome [6]; and a study by Verbeek R and colleagues showed that Lp(a) and low-density lipoprotein cholesterol are independently associated with CVD risk [7]. There is a well-documented empirical relationship between Lp(a) and CVD [8]; however, causal evidence, especially from the Chinese population, is lacking.

Confounding or reverse causation are two major problems in observational studies that render conclusions about causality uncertain. The Mendelian randomization (MR) approach has been developed to assess whether an association represents a causality in epidemiologic settings [9]. For example, the results of an observational study showed that Lp(a) is an independent risk factor associated with the progression of intima-media thickness [10], while the results from MR did not support an effect of Lp(a) variants on carotid intima-media thickness [11]. Considering the limitations of observational studies, several studies have attempted to use MR to assess the relationship between Lp(a) and CVD: Yuesong Pan and colleagues found that Lp(a) concentrations were causally associated with an increased risk of large artery stroke but a decreased risk of small vessel stroke among European individuals [12]; Helgadottir S and colleagues showed that high Lp(a) concentrations were a cause of myocardial infarction and ischemic heart disease [13]; available evidence from genetic studies supports a causal role of Lp(a) and CVD, mainly including coronary heart disease, peripheral arterial disease, aortic valve stenosis and ischemic stroke [14]. Nevertheless, these MR studies were performed in populations such as Europeans, African-Americans, and residents of Copenhagen, while limited results are available for the Chinese population. Given that Lp(a) levels are highly variable and ethnically specific [15] and that the level of Lp(a) was lower in the Han Chinese population than in Western populations [16, 17], there is a critical need for a precise characterization of the Lp(a)-associated metabolic disease risk in the Chinese population. The present study aimed to use two-sample MR to comprehensively evaluate the causal association between variants in genes affecting Lp(a) concentrations and CVD in East Asian populations.

Methods

Study design and data sources

Two-sample MR analysis was used to assess the causal effect of Lp(a) concentrations on the risk of CVD (Fig. 1). Summary statistics for Lp(a) variants were obtained from the Cohort Study on Chronic Disease of Communities Natural Population in Beijing, Tianjin and Hebei (CHCN-BTH), which is a large, prospective, population-based chronic disease study of the Communities Natural Population initiated in 2017 to explore risk factors for major chronic diseases; the design and rationale of the study have been described previously in detail [18, 19]. Individuals were excluded if there had missing data on the key demographic variables (including age and gender), Lp(a), or CVD, or if they were of non-Han ethnicity. After screening and matching, a total of 1256 individuals with coronary heart diseases and frequency-matched controls were ultimately included in the CHCN-BTH, and each was genotyped. Patients were eligible if they had been diagnosed with CVD in the hospital, and controls were residents of the same communities who did not have CVD. Data on associations between single-nucleotide polymorphisms (SNPs) and CVD were obtained from a recently published online report on the platform of genome-wide association studies (GWAS, https://www.ebi.ac.uk/gwas/) [20]. All data in this analysis were based on subjects of Han Chinese ethnicity.

Fig. 1

Conceptual framework for the Mendelian randomization analysis of Lp(a) and the risk of cardiovascular diseases

The CHCN-BTH study was approved by the Ethics Committee of the Center of Disease Control (IRB2017–003, CYCDPCIRB-20170830-1) and Capital Medical University (2018SY81), and written informed consent was obtained from all individuals prior to the baseline survey. The data on associations between SNPs and outcomes were based on publicly shared databases, and no additional participant ethical consent is required.

Laboratory procedures

Serum Lp(a) levels were determined by using a particle-enhanced turbidimetric immunoassay, Lp(a) Latex [DAIICHI] (Sekisui Diagnostic Ltd., Japan). This assay uses a highly specific monoclonal antibody and has high sensitivity. Measurements are linear in the range of 1–100 mg/dL. Three independent calibrators and two quality control materials were used to control for apo(a) size heterogeneity. The intra-assay coefficients of variation were lower than 5%, and inter-assay coefficients of variation were lower than 10%.

Selection of genetic variants

Lp(a)-related candidate genes were selected based on the National Human Genome Research Institute (NHGRI) and published GWAS data. The present study used the SNP function prediction (https://snpinfo.niehs.nih.gov/) to assess the predicted functional effects of the selected SNPs. Then, linkage disequilibrium (LD) analysis with an R² threshold of 0.8 was performed with Haploview software (https://www.broad.mit.edu/haploview/haploview) for tagging SNP selection. Under these criteria, a total of 29 SNP loci were selected in the present study for further evaluation: four SNPs (rs1018234, rs2048327, rs2457574, rs520829, and rs641990) in SLC22A3, 2 SNPs (rs14224 and rs783147) in PLG, 3 SNPs (rs1510224, rs2504921, and rs3127599) in LPAL2, 1 SNP (rs2140650) in MAP3K4, 4 SNPs (rs3798220, rs6415084, rs7765781, and rs7770628) in LPA, 2 SNPs (rs429358 and rs7412) in APOE, 2 SNPs (rs505151 and rs662145) in PCSK9, 5 SNPs (rs5925, rs5927, rs5929, rs5930, and rs688) in LDLR, and 5 SNPs (rs117052562, rs1406888, rs41269133, rs56393506 and rs9457778) in intergenic regions.

SNP genotyping

A total of 1256 blood samples were collected and centrifuged at 4000 r/min for 5 min. After centrifugation, the samples were stored at − 80 °C. Genomic DNA was extracted from 300 μL blood samples using a CWE9600 Blood DNA Kit (Kangwei Shiji Biotechnology Co., Ltd.). The Sequenom MassARRAY platform (San Diego, Calif) was applied for the candidate SNP genotyping.

Outcomes

Summary statistics for the associations of each SNP with any cardiovascular diseases in the East Asian population were obtained from a previously published GWAS. Three large observational studies of the disease traits in the Japanese population were included, with sample sizes of 212,453 [21], 128,184 [22] and 162,255 [23], respectively. This study calculated the power using Mendelian Randomization with a Binary Outcome by PASS software (https://www.ncss.com/software/pass/). The results showed that these sample size can make sure this study’s statistical power was at least 0.83. Finally, the summarized data for the associations of the 13 SNPs for Lp(a) level with any congestive heart failure (CHF) [21], ischemic stroke (IS) [21], atrial fibrillation (AF) [22], arrhythmia [21], left ventricular mass index (LVMI) [23], left ventricular internal dimension in diastole (LVDd) [23], and left ventricular internal dimension in systole (LVDs) [23] were derived; these data are presented in Additional File 1 Table S1. The clinical information was obtained from medical records, and all the clinical end points were diagnosed by physicians.

Statistical analyses

Effect sizes of the SNP and Lp(a) levels were estimated in the present study population using multiple linear regression under an additive model with adjustment for age and sex. Two-sample MR analyses were performed to compute estimates of Lp(a)-outcome (CHF, IS, AF, arrhythmia, LVMI, LVDd, and LVDs) associations using summarized data on the SNP-Lp(a) level and SNP-outcome associations. The MR-Egger regression method was used to evaluate the potential pleiotropic effects, and also a leave-one-out analysis was performed to investigate the influence of outlying or pleiotropic genetic variants. Heterogeneity across all instrumental SNPs was estimated by Cochran’s Q statistic. All analyses were carried out using R software (version 4.0.2, https://www.R-project.org/).

Results

Participant characteristics

Table 1 presents a summary of participant characteristics. The median Lp(a) concentration was 28.04 mg/dL. There were more female (n = 1028) than male participants (n = 228), and the median age of the participants was 67 for males and 62 for females. The participants were slightly overweight [median body mass index (BMI) 25.37 kg/m²]. Of the 1256 individuals in the sample, 22.7% were diagnosed with diabetes, 56.37% with hypertension, and 53.42% had hyperlipidemia. The median values for glucose, triglyceride, total cholesterol, high-density lipoprotein cholesterol and low-density lipoprotein cholesterol were 5.60 mmol/L, 1.40 mmol/L, 5.22 mmol/L, 1.41 mmol/L, and 3.02 mmol/L, respectively.

Table 1 Characteristics of the study population

Identification of Lp(a)-associated genetic variants

Only one of the 29 SNPs showed significant deviation from the Hardy-Weinberg equilibrium (HWE) (rs662145, P = 0.03321). After adjusting for age and sex, 13 SNPs were significantly associated with Lp(a) levels (P < 0.05): four SNPs in SLC22A3, 3 SNPs in LPA, 1 SNP in APOE, and 5 SNPs in undefined pathways (Table 2).

Table 2 GWAS identified SNPs Associated with Lp(a) Level

Causal association of Lp(a) with CVD

Figure 2 summarizes the main MR results exploring the causal association of Lp(a) with CVD. The results support a causal relationship between Lp(a) levels and atrial fibrillation [odds ratio (OR), 0.94; 95% confidence interval (95%CI), 0.901–0.987; P = 0.012], arrhythmia (OR, 0.96; 95%CI, 0.941–0.990; P = 0.005), the left ventricular mass index (OR, 0.97; 95%CI, 0.949–1.000; P = 0.048), and the left ventricular internal dimension in diastole (OR, 0.97; 95%CI, 0.950–0.997; P = 0.028) in the inverse-variance weighted (IVW) analysis. In contrast, no significant association of genetically elevated Lp(a) levels with congestive heart failure (OR, 0.99; 95% CI, 0.950–1.038; P = 0.766), ischemic stroke (OR, 1.01; 95%CI, 0.981–1.046; P = 0.422), and left ventricular internal dimension in systole (OR, 0.98; 95%CI, 0.960–1.009; P = 0.214) was observed. Three additional MR methods yielded consistent results (Fig. 2, Fig. 3 -1 through Fig. 3–7).

Fig. 2

Causal estimates of genetically predicted Lp(a) level in cardiovascular diseases

Fig. 3

1 Associations of Lp(a) variants with Congestive Heart Failure in different methods. Fig. 3–2. Associations of Lp(a) variants with Ischemic Stroke in different methods. Fig. 3–3. Associations of Lp(a) variants with Atrial Fibrillation in different methods. Fig. 3–4. Associations of Lp(a) variants with Arrhythmia in different methods. Fig. 3–5. Associations of Lp(a) variants with Left Ventricular Mass Index in different methods. Fig. 3–6. Associations of Lp(a) variants with Left Ventricular Internal Dimension in Diastole in different methods. Fig. 3–7. Associations of Lp(a) variants with Left Ventricular Internal Dimension in Systole in different methods

MR-Egger regression suggested no evidence of horizontal pleiotropy and no substantial heterogeneity among individual SNPs (Table 3). The results of leave-one-out sensitivity analysis showed that the association between Lp(a) levels and CHF, IS, AF, arrhythmia, LVMI, LVDd, and LVDs was not substantially driven by any individual SNP (Additional File 1 Fig. 1–1 through 1–7). Single SNP analysis for Lp(a) on CVD subtypes is also displayed in Additional File 1 Table S2 and Fig. 2–1 through Fig. 2–7.

Table 3 Horizontal pleiotropy

Discussion

Using a two-sample MR analysis design to evaluate the causal association between variants in genes encoding Lp(a) concentrations and CVD in East Asian populations, this study observed an association between genetically elevated Lp(a) and decreased risk of atrial fibrillation and arrhythmia, as well as the values of the left ventricular mass index and the left ventricular internal dimension in diastole. In contrast, no significant association was found between genetically elevated Lp(a) and the risk of congestive heart failure, the risk of ischemic stroke, and the value of the left ventricular internal dimension in systole. The results were robust in sensitivity analysis with different instruments and statistical models.

Consistent with Steffen BT and colleagues’ findings [24], no significant association between Lp(a) and heart failure was observed in the Han Chinese population. However, the results concerning heart failure were less consistent with previous observational studies [25, 26]. Elevated Lp(a) levels were found to be associated with an increased risk of incident heart failure hospitalization in the Atherosclerosis Risk in Communities (ARIC) study, with 14,154 participants, whereas after excluding the prevalent and incident myocardial infarction, the association was no longer significant [26]. Gudbjartsson DF and colleagues’ results also indicated that the odds ratio of 1.05 (95%CI: 1.02–1.07) for heart failure with a 50 nmol/L increase in Lp(a) lost significance after controlling for Lp(a) concentration (OR: 1.01, 95%CI: 0.98–1.04) [27]. Using instrumental variable analysis, Kamstrup PR and colleagues indicated that each 10-fold increase in Lp(a) levels was associated with a genetic relative risk for heart failure of 1.18 (95%CI: 1.04–1.34) [25]. Racial differences may be one of the reasons for the inconsistent findings, as it is well known that Lp(a) concentrations are highly variable with race/ethnicity, with the lowest Lp(a) levels occurring in Chinese individuals [15, 28]. When Lp(a) levels were below the 66th percentile (8 to 19 mg/dL), their effect on cardiovascular risk was not statistically significant. However, based on the data of the Cohort Study on Chronic Diseases of the General Community Population in the Beijing-Tianjin-Hebei Region (CHCN-BTH Study, Registration number: ChiCTR1900024725) project [18], the median level of Lp(a) was 11.9 mg/dL (IQR 5.9 to 23.7 mg/dL), which was in the statistically nonsignificant range in the study by Kamstrup PR and colleagues.

The findings revealed a null association between genetically predicted Lp(a) levels and ischemic stroke, which is similar to the previous results in a Mendelian randomization study published by Pan et al. [12] However, such associations were less consistent in some previous observational and MR studies [29]. In the Copenhagen City Heart Study and Copenhagen General Population Study population [30], subjects with high Lp(a) levels were associated with an increased risk of ischemic stroke both observationally (hazard ratio: 1.20; 95%CI: 1.13–1.28) and causally (hazard ratio: 1.20; 95%CI: 1.02–1.43). Gudbjartsson DF and colleagues found that the odds ratio of 1.03 (95%CI: 1.00–1.07) for stroke for a 50 nmol/L increase in Lp(a) lost significance when KIV-2 repeats were used (OR: 1.01; 95%CI: 0.98–1.04) [27]. A recent retrospective case-control study indicated that elevated Lp(a) was significantly associated with increased ischemic stroke risk (OR = 2.03, 2.36, and 2.03 for quartiles 2, 3 and 4, respectively, vs. quartile 1) in Han Chinese individuals, especially among men and younger patients [31]. Some of the factors may have contributed to the differential findings among the present study and previous studies. First, Lp(a) levels had differential associations by sex and race/ethnicity [32]; differences in participant characteristics might be one of the main reasons for the discrepancy. Second, the results of the cross-sectional study might be biased by some confounders, such as myocardial infarction [26]. Furthermore, most previous studies included all stroke subtypes [13, 33,34,35,36], while the association between Lp(a) levels and the risk of different subtypes of stroke might be different; for example, elevated Lp(a) concentrations were reported to have a positive association with large artery stroke, and an inverse association with small vessel stroke. All of these may cause variable results.

Fewer studies have reported the associations of Lp(a) with AF, arrhythmia, and the left ventricular mass index. In a prospective study of 23,738 healthy middle-aged and older women, the results indicated that Lp(a) was not associated with AF in multivariable models and that there was no significant linear trend (p for direction = 0.63) for the risk of AF across Lp(a) quintiles [37]. Aronis KN and colleagues also reported no significant association between Lp(a) and incident AF (hazard ratio: 0.98, 95%CI: 0.82–1.17) [38]. However, in the Multi-Ethnic Study of Atherosclerosis (MESA) study, the results showed that an Lp(a) level ≥ 30 mg/dL was inversely associated with the risk of AF (hazard ratio: 0.84, 95%CI: 0.71–0.99) [39], which was consistent with the results of the present MR analysis. The reasons for the disparity of the above results may lie in the demographic differences between the samples. For example, the study by Mora S and colleagues included only women [37], and the study by Aronis KN included White and Blacks participants [38]. However, in Chinese individuals, the results showed a significant difference in both cross-sectional and MR studies. Consistent with the AF results, significant associations of Lp(a) with arrhythmia, the left ventricular mass index and the left ventricular internal dimension in diastole were also observed.

Although the mechanism inversely linking high Lp(a) levels with the risk of AF and arrhythmia remains quite unclear, a similar relationship has been reported for blood lipids, in which elevated total cholesterol, low-density lipoprotein cholesterol, and high-density lipoprotein levels were found to be associated with a low risk of AF in a meta-analysis of cohort studies [40]. Inverse associations of AF with small cholesterol-poor low-density lipoprotein cholesterol, the total number of low-density lipoprotein particles, and the total number of very-low-density lipoprotein particles were driven by the numbers of cholesterol-poor small low-density lipoprotein particles and small very-low-density lipoprotein particles rather than larger cholesterol-rich LDL particles, suggesting that there might be mechanisms beyond the direct effects of cholesterol [37]. These findings could be partially explained the consideration that Lp(a) consists of a low-density lipoprotein-like particle and contains approximately 45% cholesterol [41, 42]. Further research is needed to explore the specific mechanisms in the future.

The present study has potential implications for clinical practice. Some previous studies showed that very high Lp(a) levels (above the threshold of 60 mg/dL) had a significant association with recurrent cardiovascular events in patients with premature coronary artery disease and atherosclerotic cardiovascular disease [43]. Therefore, therapeutic strategies that specifically address the reduction of Lp(a) levels are receiving a great deal of attention in clinical practice [44]. The present study results can help clinicians manage patients with CVD by identifying the risks associated with very low Lp(a) levels in patients taking Lp(a)-lowering treatment and improving the prognostic assessment. These findings could be used to provide evidence for the design of further research towards the treatment and prevention of CVD in the Chinese population.

Strengths and limitations

An advantage of this study was that the present study evaluated the associations of genetically predicted Lp(a) levels with CVD in the Han Chinese population, and the results have been limited thus far. Since the results vary among ethnic groups, this study provides evidence to enable the precise treatment of the Chinese population. Through the use of MR, the potential bias can be greatly reduced, and reverse causation can also be avoided [45]. However, the results were limited to the East Asian population, which may limit the generalizability of these findings to other populations. Some of the individuals were in a chronic disease registry, considering that Lp(a) is mainly determined by the LPA gene, there might be a slight influence on the results. In addition, due to financial constraints, additional Lp(a) related loci could not be detected; thus, the present study could not validate the relationship of Lp(a) with coronary heart disease and some other CVD subtypes. This will be pursued further in future research.

Conclusions

In conclusion, the present study provided evidence that genetically elevated Lp(a) was inversely associated with atrial fibrillation, arrhythmia, the left ventricular mass index and the left ventricular internal dimension in diastole but not with congestive heart failure, ischemic stroke, or the left ventricular internal dimension in systole. These results can help clinicians manage patients with CVD by identifying the risks associated with very low Lp(a) levels in patients taking Lp(a)-lowering treatment and improving the prognostic assessment. Further research should address the mechanism underlying these results and determine whether genetically elevated Lp(a) increases the risk of coronary heart disease or other CVD subtypes.