Causal effects from non-alcoholic fatty liver disease on kidney function: A Mendelian randomization study

Ethical considerations

This study was approved by the Institutional Review Boards of Seoul National University Hospital (E-2010-007-1160). The requirement for informed consent was waived as we investigated anonymous public databases. The investigation of the UK Biobank was approved by the consortium (No. 53799).

Concept of MR analysis

Unlike direct assessment of the effect of a SNP, MR tests the association between genetically predicted exposure and an outcome, which would be mediated by phenotypical exposure (Figure 1).¹¹ As not all individuals with a genetic predisposition develop the phenotype of interest, MR requires a large study sample to demonstrate causal estimates. The randomization process occurs before birth, and MR stratifies the population according to their genetic risk for exposure. For instance, as the rs738409 G allele has been identified to be the most significant genetic variant for NAFLD, those with the G/G rs738409 genotype would have a higher probability of developing NAFLD than those with the C/C genotype. Thus, if individuals with the rs738409 G allele have worse kidney function than those carrying other alleles, the genetic effect of the G allele on NAFLD would suggest a causal effect of NAFLD on kidney function. As one’s genotype precedes the development of clinical confounders or reverse causation, a causal inference of the identified association is possible. To validly demonstrate a causal effect, MR requires attaining three key assumptions: the relevance assumption, independence assumption, and exclusion-restriction assumption. The relevance assumption is that the genetic instrument should be strongly associated with the exposure of interest. The independence assumption and exclusion-restriction assumption indicate the absence of a horizontal pleiotropic pathway, and the assumptions are violated when the genetic effect from an instrument is not through the exposure of interest but through a confounder or another mechanism.

• Download figure
• Open in new tab

Figure 1. Study flow diagram.

The study included individual-level Mendelian randomization analysis with the UK Biobank and summary-levle Mendelian randomization analysis with the CKDGen data. To deomnstrate causal estiamtes from genetically predicted exposures, Mendelian randomization analysis requires to attain three key assumptions. The relevance assumption is that the genetic instrument should be strongly associated with the exposure of interest, and this assumption was considered attained as the study instrumented two most significant SNPs from previous genome-wide assiciation studies. The independence assumption is that the genetic effect should not be from a confounder. The exclusion-restriction assumptions means that there is no horizontal pleiotropic pathway other than a confounder (e.g. direct effect from a SNP to outcome).

Genetic instruments for NAFLD

The abovementioned SNP rs738409 is a confirmed causal variant for NAFLD and is also associated with disease severity, including steatosis, fibrosis, and NAFLD activity score.⁵ As the SNP genetically explained NAFLD to the largest degree with confirmed direct causality, we first performed a single-variant MR with the SNP. Although including a large number of variants to genetically predict a phenotype in MR has value for increasing the explained variance by genotype, the issue of horizontal pleiotropy increases according to the method.

This issue was of particular concern because other SNPs reported in a previous GWAS were associated with metabolism-related phenotypes rather than with NAFLD.⁵ Considering that kidney function and NAFLD are both closely related to metabolism-related phenotypes, we first performed MR analysis with single-variant MR, which has the advantage of demonstrating specific causal effects and minimizing the effects of the horizontal pleiotropic pathway. This was also because of the particular interest in rs738409, considering previous reports for its’ association with kidney function.

Next, to include other relevant variants for NAFLD and to implement an analysis with a recently popularized method involving pleiotropy-robust MR sensitivity analysis of a group of instrumented SNPs, which has strength reducing the possibility of weak-instrument bias, we used the results from the current largest-scale GWAS for clinically diagnosed NAFLD in the Million Veteran Program by Vojkovic M et al.^16,17 We also considered another GWAS with histologically confirmed NAFLD, but the reported 6 genetic loci were mostly associated with obesity or metabolism-related traits, raising concerns of a pleiotropic effect; thus, we implemented the larger-scale GWAS with the Million Veteran Program cohort.^5,17 To minimize the effects of different ethnic distributions, we screened 55 SNPs from the European ancestry-specific analysis, including 68,725 NAFLD cases and 95,472 controls (Supplemental Table 1). From Steiger filtering, we identified a SNP (rs12500824) with a reverse-direction effect (association with eGFR outcome, P = 2.61E-61 and with NAFLD exposure, P = 1.37E-9);¹⁸ thus, the SNP was disregarded. Regarding a horizontal pleiotropic effect from shared metabolic/circulatory risk factors between kidney function and NAFLD, we disregarded 14 SNPs that had a strong association with any phenotype of the circulatory system (e.g., hypertension) or endocrine/metabolism category (e.g., body mass index or dyslipidemia) in the phenome-wide association study with the UK Biobank data reported in a previous study (Supplemental Table 2).⁵ Thus, 40 SNPs remained as the genetic instrument with a group of SNPs for NAFLD. The genetic instrument developed from the Million Veteran Program data had particular strength, as the cohort was independent from the outcome data (UK Biobank and CKDGen) in this study, enabling us to perform two-sample MR analysis, which has strength due to its low possibility of false-positive findings.

Considerations for MR assumptions

The relevance assumption was considered attained, as the single variant MR implemented the genetic variants with strongest association with NAFLD or liver steatohepatitis. Also, the alternate genetic instrument for NAFLD consisting a group of variants was also reported from a large-scale GWAS study with association reaching genome-wide significant level.

The independence assumption was considered attained, particularly in the single-SNP MR analysis because considering the direct association with the variants with NAFLD or related phenotypes, the effect from the SNP would minimally from a confounding factor. In the genetic instrument including multiple SNPs, we carefully considered the confounding effects from metabolic disease-related phenotypes and disregarded the variants related to such confounders as those described above. In addition, we performed several pleiotropy-robust MR sensitivity analyses to support the attainment of the independence assumption.

Along with the consideration for exclusion-restriction assumption, in the current MR analysis, we paid additional attention to the possibility of a horizontal pleiotropic pathway, as some reports, despite the inclusion of only a few hundred individuals, suggest that rs738409 may directly affect kidney function.^6,8,9 If the effect from rs738409 on kidney function is direct, than as the genetic effect is not through NAFLD, the causal effect from genetically predicted NAFLD by rs738409 on kidney function could not be studied by MR. Currently, there is no direct statistical test assessing the exclusion-restriction assumption. Therefore, ws suggested in the literature, to test whether there is an independent effect of rs738409 on kidney function, we implemented the “negative control approach”, which has been suggested to be useful for testing the absence of a horizontal pleiotropic pathway when individual-level data is available.^11,19 The negative control approach re-evaluates the association between genetic instruments and outcome in those with a lower probability of developing the exposure phenotype. For instance, if an estimate of the causal effect of genetically predicted alcohol drinking behavior on a health outcome is significant in a female population, this would be less likely to indicate a true causal effect than the presence of a horizontal pleiotropic pathway. To perform this analysis, we constructed two negative control subgroups with a low probability of having the NAFLD phenotype with the available individual-level data. First, we divided subgroups according to whether the participants had a normal range body mass index (< 25 kg/m²), which is not in a range or overweight or obese and would have a small probability of having NAFLD. Second, we divided subgroups according to an alanine aminotransferase level < 20 IU/mL. Although the absolute diagnostic criterion with alanine aminotransferase level has not been determined, it would be reasonable that individuals with a low level would less like to have NAFLD. The association between genetical predisposition for NAFLD and eGFR was reassessed in the negative control group, with minimal risk of NAFLD, and their counterparts, with remaining possibility of NAFLD. If a significant association is identified only in the population with a higher probability of NAFLD, this would support that phenotypical NAFLD is the driving factor of the genetic effects of the instruments on kidney function. Conversely, if the association is similarly identified or more prominent in the negative controls, this would discourage the interpretation that phenotypical NAFLD is the mediating factor from the genetically proxied NAFLD on kidney function. Lastly, the median-based MR sensitivity analysis, which was performed in our summary-level MR analysis with multiple SNPs, relax the assumption in up to half of the genetic variants, thus has been considered as a valid sensitivity analysis for this assumption.

UK Biobank data source

The UK Biobank is a prospective cohort of > 500,000 individuals aged 50-65 years in the United Kingdom.^20,21 The database includes invaluable resources of various phenotypical data, and with deeply genotyped information, the UK Biobank data have been implemented in many genetic studies in the current literature. Regarding this study, the database has a particular strength, as other phenotypical data are available, and direct adjustment for kidney function-related metabolic comorbidities was possible. Moreover, the individual-level data enabled us to perform a negative control analysis, which was crucial for supporting the absence of a horizontal pleiotropic pathway.

MR analysis with the individual-level UK Biobank data

In the analysis, we included 337,138 white British UK Biobank participants passing the sample quality control.¹² We calculated the eGFR based on cystatin C levels, which has some superiority over creatinine-based values as they are less affected by body shape or diet,²² and the values were the outcome variable in the individual-level MR analysis. There were 321,405 individuals with available cystatin C-based eGFR levels. The first exposure was an addition of the G allele in rs738409 as a single-variant MR. The second exposure was the polygenic core (PGS) for NAFLD, calculated from the genetic instrument including 40 NAFLD-associated SNPs by PLINK 2.0. We constructed a regression model including age, sex, body mass index, waist circumference, and the first 10 genetic principal components. In addition, a covariate-adjusted model was constructed and also included diabetes mellitus, medication history for hypertension, systolic blood pressure, diastolic blood pressure, medication history for dyslipidemia, serum triglyceride level, low-density lipoprotein level, high-density lipoprotein level, income grade, and smoking history (Supplemental Methods).

The analysis in the negative control group and the counterparts was also performed. A two-sided P value < 0.05 was considered indicative of a statistically significant finding.¹²

MR analysis with the CKDGen data

To replicate the findings on a large scale, we demonstrated the causal estimates with the largest GWAS meta-analysis for eGFR performed by the CKDGen consortium.²³ The outcome summary statistics in the database, which are available in the public domain (URL: https://ckdgen.imbi.uni-freiburg.de/), have been utilized for various MR analyses in the current literature regarding kidney function.^12,13 The CKDGen GWAS summary result is independent from the UK Biobank data; thus, a similar finding in the data would be an independent replication of the main analysis.

To minimize the effects of different ethnic distributions, we limited the analysis to individuals of European ancestry (N = 567,460). The summary statistics provided in the data were for log-transformed eGFR traits, which were rescaled to a percent change unit.

In the single-variant MR, the effect from addition of the minor alleles of rs738409 on eGFR was tested by the Wald ratio method, as the analysis instrumented a single SNP, with the exposure effect size for NAFLD introduced from a previous GWAS.⁵

In the summary-level MR with multiple instrumented SNPs, we first performed the conventional inverse variance weighted method. Next, pleiotropy-robust MR sensitivity analyses were performed. First, MR-Egger regression with bootstrapped standard error, which provides pleiotropy-corrected causal estimates, was performed with MR-Egger intercept to detect the presence of a pleiotropic effect.²⁴ Second, median-based methods, including weighted median and simple median methods, were performed, and such methods have strength because they allow the presence of invalid instruments up to 50%.²⁵ Third, the contamination-mixture method was implemented, and this method has strength, as it detects groups of SNPs with similar causal estimates and performs robust MR analysis even in the presence of invalid instruments.²⁶ Fourth, MR-PRESSO analysis, which detects and corrects the effects from outliers, was performed.²⁷ The analysis was performed with the TwoSampleMR and MendelianRandomization packages in R (version 3.6.2).²⁸ The negative control approach analysis was not performed with the CKDGen data, as individual-level data to define subgroups were unavailable.

Characteristics of the data sources

The included UK Biobank participants had a median age of 58 [interquartile range 50, 63] years old, and 48% of them were male. The prevalence of obesity and diabetes was 24% and 5%, respectively. There were 20% and 17% of individuals taking medications for blood pressure control and dyslipidemia, respectively. While 5% of the individuals had the highest income grade (> 100000 ₤ per month), 21%, 26%, 26%, and 22% of individuals reported being in the 52000 to 1000000 ₤, 31000 to 51999 ₤, 18000 to 30999 ₤, and < 18000 ₤ income categories, respectively.

Individuals of European ancestry in the CKDGen dataset had a median age of 54 years old, and 50% of them were male, with a median eGFR of 91.4 mL/min/1.73 m² and a prevalence of CKD of 9%.

MR results in the UK Biobank data

The G allele of rs738409 was significantly associated with reduced eGFR values in the UK Biobank dataset (Table 1). In addition, the PGS for NAFLD showed a significant association with a lower eGFR. When the analysis was restricted to those with normal weight (body mass index < 25 kg/m²) and those with a decreased probability of NAFLD development, the association was attenuated and became nonsignificant. However, when the analysis was performed including those with body mass index ≥ 25 kg/m², genetically predicted NAFLD was significantly associated with lower eGFR levels. In another negative control group including those with an alanine aminotransferase level < 20 IU/mL, genetically predicted NAFLD was nonsignificantly associated with eGFR. On the other hand, in those with alanine aminotransferase ≥ 20 IU/mL, who would have an increased probability of NAFLD, the G allele of rs738409 or PGS for NAFLD was significantly associated with reduced eGFR values as the main findings.

MR results in the CKDGen data

In the summary-level MR from the independent CKDGen dataset, NAFLD that was genetically predicted by rs738409 was significantly associated with reduced eGFR values (Table 2). When summary-level MR was performed with multiple SNPs to predict NAFLD, all implemented MR methods, namely, the inverse variance weighted, MR-Egger regression, weighted median, simple median, contamination mixture, and MR-PRESSO methods, all provided significant causal estimates, indicating that genetically predicted NAFLD was significantly associated with reduced eGFR values.