ORIGINAL RESEARCH article

Front. Genet., 03 May 2019

Sec. Applied Genetic Epidemiology

Volume 10 - 2019 | https://doi.org/10.3389/fgene.2019.00415

Causal Effects of Genetically Predicted Cardiovascular Risk Factors on Chronic Kidney Disease: A Two-Sample Mendelian Randomization Study

  • 1. Center of System Biology and Data Information, School of Basic Medical Science, Central South University, Changsha, China

  • 2. Center of Reproductive Health, School of Basic Medical Science, Central South University, Changsha, China

  • 3. Kangda College of Nanjing Medical University, Nanjing, China

  • 4. College of Public Health, Zhengzhou University, Zhengzhou, China

  • 5. Institute of Biomedical Engineering, West China School of Basic Medical Sciences and Forensic Medicine, Sichuan University, Chengdu, China

  • 6. Tulane Center of Bioinformatics and Genomics, Department of Biostatistics and Data Science, Tulane University School of Public Health and Tropical Medicine, New Orleans, LA, United States

Abstract

Observational studies have demonstrated that cardiovascular risk factors are associated with chronic kidney disease (CKD). However, these observational associations are potentially influenced by the residual confounding, including some unmeasured lifestyle factors and interaction risk factors. Two-sample mendelian randomization analysis was conducted in this study to evaluate whether genetically predicted cardiovascular risk factors have a causal effect on the risk of CKD. We selected genetic variants associated with cardiovascular risk factors and extracted the corresponding effect sizes from the largest GWAS summary-level dataset of CKD. Cardiovascular risk factors contain high density lipoprotein (HDL) cholesterol, low density lipoprotein (LDL) cholesterol, total cholesterol (TC), triglyceride (TG), glycated hemoglobin (HbA1c), fasting glucose, systolic blood pressure (SBP) and diastolic blood pressure (DBP). A Bonferroni corrected threshold of P = 0.006 was considered as significant, and 0.006 < P < 0.05 was considered suggestive of evidence for a potential association. Genetically predicted DBP was significantly associated with CKD [odds ratio (OR) was 1.35 (95% confidence interval (CI) (1.10, 1.65); P = 0.004)]. There was suggestive evidence for potential associations between genetically predicted higher HDL cholesterol [OR: 0.88, 95%CI (0.80, 0.98), P = 0.025] and lower adds of CKD, and between higher SBP [OR: 1.36, 95%CI (1.07, 1.73), P = 0.013] and higher adds of CKD. However, genetically predicted LDL cholesterol, TC, TG, HbA1c, and fasting glucose did not show any causal association with CKD.

Introduction

Chronic kidney disease is a global-health challenge affecting at least 225.7 million people in the world, especially in developing countries, in 2010 (). The prevalence rates of CKD have dramatically increased in the past several decades, which threaten human health seriously (). Observational studies have detected a close relationship between increased risk of cardiovascular diseases and CKD, showing that major cardiovascular diseases account for approximately 50% of the causes of death in CKD patients (). Multiple cardiovascular risk factors, including lipids (), glycemic traits (), and blood pressure (; ), have been reported associated with CKD in observational studies. However, a review reported a reverse effect of CKD on hypertension and dyslipidemia (), demonstrating that CKD can promote hypertension and dyslipidemia, which in turn can lead to the progression of renal failure (). The reverse association between cardiovascular risk factors and CKD might be influenced by potential confounders of unmeasured lifestyle factors.

Mendelian randomization analysis, analogous to randomized controlled trials (RCTs), has been widely performed to investigate the potential causality between genetically predicted environmental factors and diseases (; ; ). As the genotype is randomly assigned during the meiotic process, MR analysis results will not be distorted by confounders (), a major limitation of traditional observational studies. Compared with one-sample MR, which extracts the effect estimates for IV-exposure association and IV-outcome association from the same sample, TSMR estimates that these associations in different samples and the estimates are then combined to infer the potential exposure-outcome causal association ().

Recently, a TSMR analysis demonstrated that higher HDL cholesterol concentration was causally associated with better kidney function, while LDL cholesterol or triglyceride concentration showed no association with kidney function (). Another MR analysis detected that genetically determined type 2 diabetes (T2D) was causally associated with decreased eGFR in populations of Chinese aged above 40 (). However, no studies have been conducted utilizing TSMR to summarize the causal relationship between cardiovascular risk factors and CKD. Therefore, in this study, we performed TSMR analysis to examine the causal effect of cardiovascular risk factors, including lipids, glycemic traits, and blood pressure on CKD, based on GWAS summary-level data.

Materials and Methods

Data Sources

We selected genetic variants associated with cardiovascular risk factors, including lipids, glycemic traits and blood pressure, and then extracted the corresponding effect sizes for CKD using the largest GWAS summary-level dataset (; ; ; ; ). No ethical approval was conducted in our study due to this being a re-analysis based on previous collected and published data. Lipids data of HDL cholesterol (n = 187, 167), LDL cholesterol (n = 173, 082), TC (n = 187, 365), and triglyceride (TG, n = 177, 861) were extracted from the GLGC consortium (). We obtained glycemic traits data of fasting glucose () (n = 133, 010) and HbA1c (n = 46, 368) () from the MAGIC consortium. And blood pressure data () of SBP (n = 317, 754) and DBP (n = 317, 756) were extracted from UK Biobank imputed genotype data1. Additionally, people with eGFR based on serum creatinine (eGFRcrea) < 60 mL/min/1.73 m2 were defined as CKD, and the dataset were acquired from CKDGen consortium (n = 117, 165) (). All of the dataset were conducted on European populations from RCTs and population-based cohorts. Genomic control to each sample was applied to correct for inflated test statistics due to potential population stratification in our datasets. Age, sex, and body mass index were all adjusted in regression models of GLGC, MAGIC and UK Biobank, and age and sex were also corrected in CKDGen (; ; ; ; ) (Supplementary Table S1).

Study Design

The genetic variants used as IVs in TSMR analysis must satisfy three assumptions as follows (Figure 1): (1) IVs are strongly associated with cardiovascular risk factors, including lipids, glycemic traits, and blood pressure. (2) The IVs are independent of any known confounders. (3) The selected IVs are conditionally independent of CKD, given cardiovascular risk factors and confounders. The second and third assumptions are known as independence from pleiotropy (). In this study, we used cardiovascular risk factors as the exposures, including HDL cholesterol, LDL cholesterol, TC and TG for lipids, HbA1c and fasting glucose for glycemic traits, SBP and DBP for blood pressure, and CKD as the outcome to perform the TSMR analysis.

FIGURE 1

IVs Selection and Validation

Instrumental variables must be associated with cardiovascular risk factors, including HDL cholesterol, LDL cholesterol, TC, TG, HbA1c, fasting glucose, SBP, and DBP. To ensure the close relationship between IVs and cardiovascular risk factors, we selected variants with P < 5 × 10−8 in the corresponding GWAS summary-level dataset. In addition, pairwise-linkage disequilibrium (LD) was calculated by PLINK 1.90 () to ensure the independence among selected IVs and SNPs with r2> 0.001 will be removed from our analysis.

Selected IVs should be independent of any known confounders and conditionally independent of CKD, given the related traits of cardiovascular risk factors. The two assumptions indicate that IVs must influence CKD only through cardiovascular risk factors rather than another pathway (). Firstly, we obtained the corresponding effect estimates of these variables on CKD. For the SNPs that were not available in the CKD, we used proxy SNPs that were highly correlated (r2 > 0.8) based on the SNP Annotation and Proxy (SNAP) search system2 (). Secondly, MR-Egger regression was performed to assess the horizontal pleiotropic (), which will be introduced in the following statistic analysis. Additionally, we excluded any palindromic SNPs that have minor allele frequency above 0.42 to ensure that the effects of the SNPs on the exposures correspond to the same allele as their effects on CKD (). Moreover, to adjust for potential confounding, GWAS Catalog was used to check for the associations between selected IVs, smoking and type II diabetes (H1bAc and fasting glucose were excluded). Additionally, inspired by the idea of , we used the F statistic to investigate the association of selected IVs with the exposure on a web application3.

Pleiotropy Assessment

Mendelian randomization-Egger regression was performed to assess the horizontal pleiotropic pathway between IVs and CKD, independent of cardiovascular risk factors (). MR-Egger regression was developed from Egger regression, which has been used to examine the publication bias in meta-analysis (ref needed here). This approach is expressed as αi = βγi0, αi represents the estimated effect between IVs and CKD; γi represents the estimated effect between IVs and cardiovascular risk factors, including HDL cholesterol, LDL cholesterol, TC, TG, HbA1c, fasting glucose, SBP, and DBP; slope β represents the estimated causal effect of cardiovascular risk factors on CKD; intercept β0 could be explained as the estimated average value of horizontal pleiotropic. Intercept with P > 0.05 indicates no horizontal pleiotropic exists. Additionally, the slope estimate provides the pleiotropy-corrected causal effect. However, this estimate may be underpowered if the selected SNPs collectively fail to explain a large proportion of the variance in the exposure ().

TSMR Analysis

In this study, IVW method was used for TSMR analysis to estimate the causal effect between cardiovascular risk factors and CKD (). The causal effect β was estimated as wiii), where i refers to the ith IV, αi defines as the association effect of IVs on CKD, γi represents the association effect of IVs on cardiovascular risk factors, and wi means the weights of the causal effect of cardiovascular risk factors on CKD. MR Steiger test was also performed to infer the causal direction between exposures and CKD. It calculates the variance explained in the exposures and outcome by the instrumenting SNPs, and tests if the variance in the outcome is less than the exposures. Given the multiple testing situation, we used a conservative approach and applied a Bonferroni corrected significance level of 0.006 (0.05/8). 0.006 < P < 0.05 was considered as suggestive evidence for a potential association.

Robust Adjusted Profile Score

The TSMR might fail if the selected SNPs are weak instruments. Therefore, we carried out a recently proposed method called robust adjusted profile score (RAPS) () which considers the measurement error in SNP-exposure effects and is unbiased even when there are many (e.g., hundreds of) weak instruments, and is robust to systematic and idiosyncratic pleiotropy. Detailed information about this method please refer to the original paper ().

Positive Control and Negative Control

Previous studies already reported established causal relationship between coronary artery disease (CAD) and CKD (), while there is little evidence indicating that cardiovascular risk factors are associated with myopia. To further demonstrate the validity of the selected IVs, we included coronary heart disease and myopia as positive and negative controls in our analysis. The summary statistics of them were respectively derived from CARDIoGRAMplusC4D consortium () and UK Biobank imputed genotype data1, including 184, 305 and 335, 700 individuals from European population.

Sensitivity Analysis of TSMR

In the current study, weighted median and simple median methods were also applied as the follow-up sensitivity analysis (). Compared with IVW, the weighted median and simple median methods have greater robustness to individual genetics with strongly outlying causal estimates and would generate a consistent estimate of the causal effect when valid IVs exceed 50% (; ). Furthermore, leave-one-out sensitivity analysis was performed to identify if the association was disproportionately influenced by a single SNP. The TSMR analysis is performed again but leaving out each SNP in turn and the overall analysis including all SNPs was shown for comparison (). All of the analysis was implemented by the “TwoSampleMR” package in R software environment.

Results

Selection and Validation of IVs

We obtained 84, 66, 78, 54, 10, 31, 103, 111 LD-independent (r2 < 0.001) IVs in total that achieved genome-wide significance (P < 5 × 10−8) from HDL cholesterol, LDL cholesterol, TC, TG, HbA1c, fasting glucose, SBP and DBP datasets, respectively. Not all of the SNPs were directly found in the CKD dataset; detailed information of all independent IVs in this TSMR analysis were shown in Supplementary Table S2. In addition, the intercept term, estimated for the exposures from MR-Egger regression, demonstrated that no horizontal pleiotropic exists in our TSMR analysis (Table 1). Moreover, we did not detect a direct association between selected IVs and smoking or type II diabetes (H1bAc and fasting glucose were excluded) in GWAS Catalog.

Table 1

ExposuresOutcomeMR-Egger regression
intercepts (95% CI)P-value
HDL cholesterolCKD−0.005 (−0.014, 0.004)0.282
LDL cholesterolCKD0.007 (−0.001, 0.016)0.080
TCCKD0.002 (−0.007, 0.011)0.687
TGCKD0.006 (−0.004, 0.016)0.241
HbA1cCKD−0.006 (−0.034, 0.023)0.707
Fasting glucoseCKD−0.006 (−0.021, 0.01)0.463
DBPCKD−0.003 (−0.018, 0.013)0.746
SBPCKD−0.006 (−0.022, 0.011)0.503

Mendelian randomization (MR)-Egger regression intercepts.

A significant result (P > 0.05) indicates that the y-intercept of the MR-Egger regression line is not significantly different from zero and thus no pleiotropy exists. MR, mendelian randomization; HDL, high density lipoprotein; CKD, chronic kidney disease; LDL, low density lipoprotein; TC, total cholesterol; TG, triglyceride; SBP, systolic blood pressure; DBP, diastolic blood pressure; CI, confidence interval.

F statistics were presented to demonstrate the strength of relationship between IVs and exposures, and F statistics greater than 10 are often considered as strong enough to mitigate against any bias of the causal IV estimate. Our selected IVs showed strong strength with F statistics ranging between 802 and 9447 (Supplementary Table S3).

TSMR Analysis

According to the IVW analysis results, the odds ratio (OR) and 95% confidence interval (CI) per 1-SD increase of DBP within CKD was 1.35 (1.10, 1.65) (P = 0.004). There was suggestive evidence for potential associations between genetically predicted higher HDL cholesterol [OR: 0.88, 95%CI (0.80, 0.98), P = 0.025] and lower adds of CKD, and between higher SBP [OR: 1.36, 95%CI (1.07, 1.73); P = 0.013] and higher adds of CKD. However, genetically predicted LDL cholesterol, TC, TG, HbA1c and fasting glucose were not associated with CKD (Figure 2).

FIGURE 2

Robust Adjusted Profile Score

The results turned out to be consistent with the TSMR results that increased HDL cholesterol decreases the risk of CKD [OR = 0.89, 95% CI (0.79, 0.98), P = 0.025], but increased SBP [OR = 1.37, 95% CI (1.17, 1.56), P = 0.002] and DBP [OR = 1.33, 95% CI (1.14, 1.51), P = 0.003] increases the risk of CKD (Figure 2).

MR Steiger Directionality Test

The inferred causal direction between exposures (HDL cholesterol, LDL cholesterol, TC, TG, HbA1c, fasting glucose, DBP, SBP) and CKD were “TRUE” in our MR Steiger test (Supplementary Table S4).

Sensitivity Analysis of TSMR

In sensitivity analysis, there was suggestive evidence for potential associations between higher genetically predicted DBP [Weighted median: OR = 1.37, P = 0.013, 95% CI (1.07, 1.76); Simple median: OR = 1.32, P = 0.028, 95% CI (1.03, 1.70)] and SBP [Weighted median: OR = 1.30, P = 0.050, 95% CI (1.10, 1.70); Simple median: OR = 1.45, P = 0.007, 95% CI (1.11, 1.89)] and higher odds of CKD. However, genetically predicted HDL cholesterol, LDL cholesterol, TC, TG, HbA1c and fasting glucose were not associated with CKD (Figure 2). Furthermore, leave-one-out analysis showed a consistent significant causal effect of cardiovascular risk factors on CKD, supporting the robustness of our IVW analysis findings (Supplementary Table S5).

Positive and Negative Control

Consistent with the conclusion of previous published paper (), our analysis demonstrated a causal relationship between all exposures and the positive control (coronary artery disease). As for negative control (myopia), our analysis showed no causal relationship between all exposures and myopia (Supplementary Table S6). Independent IVs of CAD and myopia were listed in Supplementary Tables S7, S8, respectively. Both negative and positive control analysis results suggest that the selected IVs of CKD are appropriate.

Discussion

In the current study, by using genetic variants associated with cardiovascular risk factors as proxies, our TSMR analysis confirmed the causal association between DBP and increased risk of CKD. Genetically predicted SBP showed suggestive evidence for a possible causal association with CKD. We also detected suggestive evidence for an inverse causal association between genetically predicted HDL cholesterol and CKD. However, no evidence was found to support the associations between genetically predicted LDL cholesterol, TC, triglycerides, HbA1c, fasting glucose, and CKD.

Our conclusion of the causal relationship between LDL cholesterol, TC and CKD was consistent with the results of . They also reported MR evidence for the association between genetically predicted higher HDL cholesterol and lower odds of CKD, but we did not draw the same conclusion in our sensitivity analysis (). This inconsistency may be caused by a different selection of IVs and different analysis methods applied. Previous MR analysis, based on the GLGC and CKDGen dataset, reported a causal effect of HDL cholesterol on kidney function (). However, our study focused on CKD (eGFR < 60 mL/min/1.73 m2) rather than kidney function, which may be more directional compared with eGFR, due to the CKD’s clearly stages definition in clinical diagnosis(). Another MR analysis performed in a Chinese population demonstrated a causation between T2D and CKD, but we did not detect a significant causal relationship between glycemic traits (HbA1c and fasting glucose) and CKD in a European population (), which demonstrates the importance of ethnicity in MR analysis.

Observational studies demonstrated that hypertension is almost invariably present in patients with renal failure, which is a powerful risk factor for cardiovascular diseases and CKD (). The renin–angiotensin system and sympathetic nervous system have been considered as important mechanisms involved in the elevation of blood pressure in subjects with CKD (). The elevation of plasma catecholamine might increase nerve sympathetic traffic in CKD. Serving as a regulator of blood pressure (; ), renalase could metabolize catecholamines through dopamine, epinephrine, and norepinephrine (). Investigators have discovered that renalase is an oxidase expressed mainly in glomeruli and proximal tubules of the kidney and cardiomyocytes (). All the detected potential mechanisms suggest the relationship between blood pressure, cardiovascular disease and CKD, but the mechanism remains unclear (). Our current study successfully detected a causal relationship between genetically predicted blood pressure and CKD, which may provide novel evidence to further explain the mechanism of hypertension in CKD.

There are several strengths in our current study. First, we included multiple cardiovascular risk factors (8 in total) as the exposures, hence we were able to include a relatively large number of IVs. Then, to provide relatively consistent causal effect estimates, we also performed sensitivity analysis using several different approaches. In addition, we used a web application3 to investigate the magnitude of bias arising from sample overlap with a conservative value of concentration parameter, and it would not be substantial due to sample overlap in our study. Furthermore, CAD and myopia were used as positive and negative controls, respectively, to demonstrate the validity of selected IVs. However, this study also has some limitations. Firstly, as we only used summary statistics and had no access to the original individual clinical outcome measures, we could not conduct analyses stratified by subtypes of CKD. Secondly, different standards of quality control in individual-level GWAS may affect our results. Therefore, the results cannot be easily generalized.

Using a genetic approach, we found DBP is causally associated with CKD risk. Furthermore, we provided suggestive evidence that SBP is causally associated with CKD risk and HDL cholesterol is inversely causally associated with CKD. However, additional human and animal studies are still needed to further confirm our TSMR results.

Statements

Author contributions

H-ML and H-WD: research idea and study design. H-ML: data acquisition, statistical analysis, and writing the manuscript. H-ML, H-WD, and QH: data analysis and interpretation. G-YS and H-MX: supervision or mentorship. H-ML, QZ, B-YL, W-DS, XQ, and W-QL: new manuscript revision and point to point response to the reviewers. All authors read and approved the final manuscript.

Funding

H-WD was partially supported by the National Institutes of Health (Grant Nos. P50AR055081, R01AR057049, R01AR059781, D43TW009107, P20GM109036, R01MH107354, R01MH104680, R01GM109068, and R01AR069055), the Edward G. Schlieder Endowment fund to Tulane University. We appreciate the support from Central South University in providing necessary support for this collaborative project.

Acknowledgments

The authors would like to thank MR-base dataset for providing relevant publicly available summary statistics. Thanks to all the consortiums that provide public data, which have been mentioned in the data sources in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00415/full#supplementary-material

Abbreviations

  • CI

    confidence interval

  • CKD

    chronic kidney disease

  • DBP

    diastolic blood pressure

  • eGFR

    estimated glomerular filtration rate

  • GWAS

    genome wide association study

  • HbA1c

    glycated hemoglobin

  • HDL

    high density lipoprotein

  • IVs

    instrumental variables

  • IVW

    Inverse variance weighted

  • LD

    linkage disequilibrium

  • LDL

    low density lipoprotein

  • MR

    mendelian randomization

  • OR

    odds ratio

  • SBP

    systolic blood pressure

  • TC

    total cholesterol

  • TG

    triglyceride

  • TSMR

    two-sample mendelian randomization.

References

Summary

Keywords

two-sample mendelian randomization, genome-wide association study, cardiovascular risk factors, chronic kidney disease, causation

Citation

Liu H-M, Hu Q, Zhang Q, Su G-Y, Xiao H-M, Li B-Y, Shen W-D, Qiu X, Lv W-Q and Deng H-W (2019) Causal Effects of Genetically Predicted Cardiovascular Risk Factors on Chronic Kidney Disease: A Two-Sample Mendelian Randomization Study. Front. Genet. 10:415. doi: 10.3389/fgene.2019.00415

Received

12 January 2019

Accepted

16 April 2019

Published

03 May 2019

Volume

10 - 2019

Edited by

Robert Klein, Icahn School of Medicine at Mount Sinai, United States

Reviewed by

David Fardo, University of Kentucky, United States; Jing Hua Zhao, University of Cambridge, United Kingdom

Updates

Copyright

*Correspondence: Hong-Wen Deng,

This article was submitted to Applied Genetic Epidemiology, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics