Skip to main content

ORIGINAL RESEARCH article

Front. Nutr., 29 March 2022
Sec. Nutritional Epidemiology

Causal Association Between Tea Consumption and Kidney Function: A Mendelian Randomization Study

\r\nYangchang Zhang,,,&#x;Yangchang Zhang1,2,3,4†Yang Xiong&#x;Yang Xiong5†Shisi ShenShisi Shen6Jialu Yang,,Jialu Yang1,2,3Wei Wang,,Wei Wang1,2,3Tingting WuTingting Wu7Li ChenLi Chen8Qiuhua YuQiuhua Yu9Hangjia ZuoHangjia Zuo6Xu Wang,,Xu Wang1,2,3Xun Lei,,,*\r\nXun Lei1,2,3,4*
  • 1School of Public Health and Management, Chongqing Medical University, Chongqing, China
  • 2Research Center for Medicine and Social Development, Chongqing Medical University, Chongqing, China
  • 3The Innovation Center for Social Risk Governance in Health, Chongqing Medical University, Chongqing, China
  • 4Research Center for Public Health Security, Chongqing Medical University, Chongqing, China
  • 5The West China Hospital, Sichuan University, Chengdu, China
  • 6The First School of Clinical Medicine, Chongqing Medical University, Chongqing, China
  • 7Chongqing Collaborative Innovation Center for Functional Food, Chongqing University of Education, Chongqing, China
  • 8School of Public Health & Institute of Child and Adolescent Health, Peking University, Beijing, China
  • 9Department of Critical Care Medicine, Peking Union Medical College Hospital, Peking Union Medical College and Chinese Academy of Medical Sciences, Beijing, China

Background: Causal research concerning the consumption of tea and the risk of chronic kidney disease (CKD) is limited. This study identified the potential causal effects of tea intake on CKD, the estimated glomerular filtration rate (eGFR), and albuminuria.

Methods: Genome-wide association studies (GWASs) from UK Biobank were able to identify single-nucleotide polymorphisms (SNPs) associated with an extra cup of tea each day. The summary statistics for the kidney function from the CKDGen consortium include 11,765 participants (12,385 cases of CKD) and 54,116 participants for the urinary albumin-to-creatinine ratio who were mostly of European descent. A two-sample Mendelian randomization (MR) analysis was performed to test the relationship between the selected SNPs and the risk of CKD.

Results: A total of 2,672 SNPs associated with tea consumption (p < 5 × 10–8) were found, 45 of which were independent and usable in CKDGen. Drinking more cups of tea per day indicates a protective effect for CKD G3-G5 [odds ratio (OR) = 0.803; p = 0.004] and increases eGFR (β = 0.019 log ml/min/1.73 m2 per cup per day; p = 2.21 × 10–5). Excluding two SNPs responsible for directional heterogeneity (Cochran Q p = 0.02), a high consumption of tea was also negatively correlated with a lower risk of albuminuria (OR = 0.758; p = 0.002).

Conclusion: From the perspective of genes, causal relationships exist between daily extra cup of tea and the reduced risk of CKD and albuminuria and increased eGFR.

Introduction

Chronic kidney disease (CKD) is an increasingly serious public health issue, which affects 9–16% of the world’s population (1). Globally, the number of deaths attributed to CKD has increased to approximately 1.2 million, and the mortality rate in all ages has increased by 41.5% from 1990 to 2017 (2). The time-series model from the global burden of disease (GBD) group indicated an increasing trend in the number of years of life lost, from around 26 million annually in 2016 to 52.5 million in 2040 (3). Although treatment has been shown to slow the progression, CDK may proceed at different paces over time (4). CKD can progress to kidney failure and early cardiovascular disease in the end stage, namely, end-stage renal disease (ESRD). Dialysis or kidney transplant is necessary for survival on account of dysfunctional kidneys (4). CKD is associated with an increased risk of many other conditions, including cognitive impairment, renal bone disease, chronic anemia, and death by sepsis and cardiovascular disease (5, 6). Empirical studies have found that lifestyle factors, such as smoking, alcohol consumption, and obesity, are related to a higher risk of the disease (79). Aside from these factors, an increasing interest is being paid to the important role of diet (10). As an essential part of the diet, beverages can influence the maintenance of general health and renal function and inhibit/abet the high-risk factors that may lead to CKD (e.g., hypertension, obesity, and diabetes) (10, 11).

Tea is one of the most popular beverages in the world (12). Many longitudinal and cross-sectional studies have investigated the association between tea consumption and the risk of CKD and the estimated glomerular filtration rate (eGFR), and studies on the association between tea intake and albuminuria are limited (1216). Nevertheless, these studies failed to reach a robust and consistent conclusion. This discrepancy may be due to diverse tea types (e.g., green, black, and oolong tea) and demographic heterogeneity (e.g., sex) (16). For example, no evidence showed that extra tea intake could improve renal function or delay the stage of nephropathy in Iran, Singapore, and Netherlands (12, 14, 16). However, another follow-up survey reported that tea intake can exert positive effects on patients with metabolic syndrome (MetS) (13). Additionally, a study from China found that oolong tea could obviously promote the efficiency of eGFR compared to black tea and green tea (15). In view of the inconsistent results and different indicators of previous studies, the effects of tea consumption on the different indicators of renal function need to be investigated. Additionally, considering the limitations of observational studies and the potential influence of confounding factors or reverse causality, a causal relationship between tea consumption and the risk of abnormal CKD indicators cannot be concluded. Therefore, whether a causal relationship exists between tea intake and kidney function remains unclear.

Mendelian randomization (MR) is a genetic epidemiology design, which improves the power of causal inference by applying proxy germline genetic variants as instrumental variables for exposure (e.g., tea intake) on an outcome (e.g., CKD) (17). Single-nucleotide polymorphism (SNP) sites are randomly assigned in the beginning, so the bias of reverse causation and residual confounding is avoided (18). The two-sample MR design has not yet been used to determine the causal association between tea intake and the risk of abnormal renal function. Thus, a conventional MR method was performed in this study to estimate whether tea intake is causally associated with the risk of abnormal renal function (e.g., CKD, eGFR, and albuminuria).

Materials and Methods

Genetic Instrument Selection

The genome-wide association study (GWAS) summary data set for tea consumption (Phenotype Code:1488_raw) based on UK Biobank was downloaded from Neale Lab1 (GWAS round 2), including over 349,376 samples of European ancestry. GWAS was adjusted for sex, age, age2, sex × age, sex × age2, and first 20 ancestry principal components. The data of habitual tea consumption were obtained as a baseline from a dietary questionnaire. In the questionnaire, the following question was asked: “How many cups of tea do you drink each day (including black and green tea)?” GWAS summary statistics for unconverted daily tea consumption were used to identify SNPs associated with tea drinking. Detailed information regarding the phenotype and the process of quality control in UK Biobank is available on the Neale Lab website.2 We selected autosomal biallelic SNPs with p < 5 × 10–8 and performed further quality control based on minor frequency >1%, ending up with 2,672 unique SNPs. Furthermore, based on the European sample reference data from the 1000 Genomes Project (19), we clumped these 2,672 SNPs with linkage disequilibrium r2 < 0.001 at a 10,000 kb window, confirming the independence of the selected genetic variants. Finally, 45 independent SNPs were associated with tea consumption. The proportion of variance in tea consumption explained by each SNP was estimated using the R2 value (20), and the instrumental strength of each SNP was assessed using the F-statistic (21). Detailed information on the relationship between the selected SNPs and exposures is shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Characteristics of SNPs associated with tea consumption.

Genetic Summary Data of Kidney Function

We extracted summary statistics for CKD and eGFR from the CKDGen consortium. The meta-analyses of GWAS of the kidney function comprised 43 studies of European ancestry (n = 117,165, 12,385 cases) (22). Participants were diagnosed with CKD GFR categories 3–5 (G3–G5) based on eGFR < 60 ml/min/1.73 m2 (23). Except the two studies that reported data, all definitions of CKD G3-G5 came from a single assessment of eGFR (22). eGFR is calculated by the Schwartz formula (<18 years) (24) and the Chronic Kidney Disease Epidemiology Collaboration equation for adults (>18 years) (25). Wuttke et al. reported the characteristics of the CKDGen alliance study (26). Albuminuria was also extracted from the CKDGen consortium, which enrolled a total of 54,116 participants of European ancestry (27). Albuminuria is defined as a condition where the urinary albumin–creatinine ratio is >17 mg/g (>1.92 mg/mmol) in men and >25 mg/g (>2.83 mg/mmol) in women (27). Sex differences in albuminuria were obtained from the study of Warram et al. (28). All summary data can be obtained from the UK Medical Research Council Integrative Epidemiology Unit Open GWAS Project database.3

All research analyses were based on publicly available GWAS summary statistics, and no additional ethics approval and informed consent were required.

Statistical Analyses

The conventional MR method was applied in this study4 (29). The fixed-effect inverse-variance-weighted (IVW) model was used to examine a causal association, and it was considered as the main analytical method (30, 31). The IVW method included individual MR effects of SNPs to derive overall weighted effects of the potential causal association. Multiplicative random IVW was used as a supplement if necessary. Furthermore, the forest plots showed the MR-derived odds ratio (OR) of CKD and albuminuria, and log odds of eGFR were predicted by genes when an additional cup of tea is consumed daily. “Leave-one-out” analyses were performed to determine whether a causal association was reliant on any single SNP. Causal Analysis Using Summary Effect estimates (CAUSE) account for uncorrelated and correlated pleiotropy simultaneously (32).

The IVW method assumes that all genetic variants satisfy the three assumptions of the instrumental variables: (1) closely associated with tea intake; (2) not associated with confounders of the association between tea intake, CKD, eGFR, and albuminuria; and (3) risk factors associated with the risks of CKD, eGFR, and albuminuria, which could only be induced via tea consumption (33, 34). F-statistics were used to test weak instrumental variables, F-statistics = (N − 2) × R2/(1 − R2), where R2 is the variance in tea consumption explained by the genetic instrument, k is the number of genetic variants, and n is the sample size. F > 10 proved to be a strong genetic instrument in the MR study. Cochran’s Q test was used to quantify the size of heterogeneity effect between the genetic instruments (35), which might indicate that the potential horizontal pleiotropy violated the third MR assumption. Potential violation of the second and third MR assumptions was tested using several approaches, such as the MR-Egger regression (36) and the weighted median (37) methods. The intercept from MR-Egger also provides a formal test for directional pleiotropy. The association between the selected SNPs and exposures was validated in the PhenoScanner database5 (Supplementary Table 1). SNPs associated with the traits other than tea intake were recorded at the significance level (p < 5 × 10–8).

MR-Egger is an adaption of Egger regression, which realizes directional pleiotropy by introducing an intercept in the weighted regression model. When the value of the intercept term is away from 0, it indicates horizontal pleiotropy (38). Based on this approach, unbiased estimates are performed in the presence of pleiotropic instruments, assuming that the magnitude of pleiotropic effects is independent of the size of the instrumental variables—SNPs associated with tea intake (38).

The weighted median method uses each instrumental variable to weigh the estimated value of the reciprocal of its variance to rank them in MR. The median result is selected, and the single MR estimated value and CIs based on the bootstrapping technique are displayed (39). The weighted median requires and assumes that at least half of the instruments are effective (40).

Steiger-MR was used to test whether the SNPs explained significantly more variance in exposure than the outcome (the opposite may indicate reverse causation) (41).

Bonferroni correction (p = 0.05/3 outcomes/3 methods) was applied to adjust multiple testing (p = 0.006) in univariable MR. The “TwoSampleMR” package (version 0.5.5), “CAUSE” package (version 1.2.0), and R software version 3.6.1 were used for all statistical analyses.

Results

Causal Association Between Tea Consumption and CKD G3-G5

The estimation of the causal effect of tea consumption on CKD G3-G5 from the MR analyses is shown in Figure 1. In the IVW MR analyses, the OR of CKD for an additional daily cup of tea was 0.803 (p = 0.004). The MR-Egger test detected no directional pleiotropy (p = 0.732). In the leave-one-out analyses, the estimates ranged from 0.779 to 0.830, suggesting that the observed result was not the effect of a single SNP (Supplementary Figure 1). The weighted median analysis (OR = 0.824; p = 0.031) of the estimated values was concordant and similar in size, but it could not confirm the protective effect of tea against CKD G3-G5 in the adjusted value of p (Figure 1). The estimates from the MR-Egger (OR = 0.758; p = 0.187) were not statistically significant, suggesting limited effectiveness. Additionally, there was no indication of pleiotropy when the intercept was derived from the MR-Egger regression (Egger intercept: 0.006, the value of p: 0.732). The value of p for the MR-Egger method was 0.146, which also suggested no significant sign of heterogeneity. The causal effect of tea intake on the CKD estimated by CAUSE was unsupported (p = 0.132). The scatter plot of the SNP—CKD associations against SNP—tea associations is shown in Figure 2. The forest plots of tea—CKD estimates in each SNP are presented in Supplementary Figure 2.

FIGURE 1
www.frontiersin.org

Figure 1. Forest plots of Mendelian randomization (MR) study using genetically predicted tea consumption with chronic kidney disease (CKD), glomerular filtration rate (GFR), and albuminuria. Inverse-variance-weighted (IVW), MR-Egger, weighted median, simple, and weighted mode were used in this study. *Denotes the removal of two single-nucleotide polymorphism (SNPs; rs1030510 and rs4410790) that might give rise to significant heterogeneity (MR-Egger intercept p = 0.02).

FIGURE 2
www.frontiersin.org

Figure 2. Scatter plot of the effect size for each SNP on tea consumption and CKD G3–G5.

Causal Association Between Tea Consumption and eGFR

Inverse-variance-weighted analyses of tea consumption and eGFR provided correlative evidence for an association (β = 0.019 log ml/min/1.73 m2 per cup per day; p = 2.21 × 10–5) (Figure 1). In the leave-one-out analyses, β ranged from 0.015 to 0.020 (Supplementary Figure 3). There was no evidence of directional pleiotropy (MR-Egger intercept p = 0.183) and horizontal pleiotropy (heterogeneity p = 0.521). This was consistent with the estimates from the weighted median (β = 0.020; p = 1.00 × 10–4) and MR-Egger (β = 0.029; p = 1.34 × 10–2) analyses (original p < 0.05; adjusted p > threshold value), which were more robust to pleiotropy (Figure 1). The causal effect of tea intake on eGFR estimated by CAUSE was insignificant (p = 0.141). The scatter plot of the SNP—eGFR associations against the SNP—tea associations is shown in Figure 3. The forest plot of tea—eGFR estimates in each SNP is presented in Supplementary Figure 4.

FIGURE 3
www.frontiersin.org

Figure 3. Scatter plot of the effect size for each SNP on tea consumption and the estimated glomerular filtration rate (eGFR).

Causal Association Between Tea Consumption and Albuminuria

Initial IVW analyses of tea consumption and albuminuria did not provide strong evidence for their correlation (OR = 0.749; p = 0.009) (Figure 1). There was evidence of directional pleiotropy (MR-Egger intercept p = 0.019) and no horizontal pleiotropy (heterogeneity p = 0.575). Based on the leave-one-out method, these outlying SNPs (rs1030510 and rs4410790), which mainly caused heterogeneity, were removed. Then, through the random-effect IVW (OR = 0.758; p = 0.002) and weighted median (OR = 0.701; p = 0.027) methods (original p < 0.05; adjusted p > threshold value), the estimate of the causal effect of tea consumption on albuminuria was similar in direction and magnitude to CKD G3-G5 (Figure 1). In the leave-one-out analyses, ORs ranged from 0.723 to 0.894, showing a consistency in the entire estimate (Supplementary Figure 5). There was neither significant horizontal pleiotropy (heterogeneity p = 0.893) nor directional pleiotropy (MR-Egger test p = 0.383). The scatter plot of the SNP—albuminuria associations against the SNP—tea associations before and after removing two SNPs is shown in Figures 4A,B. The causal effect of tea intake on albuminuria estimated by CAUSE was insignificant (p = 0.620). The forest plots of tea—eGFR estimates in each SNP after removing outlying SNPs are presented in Supplementary Figure 6.

FIGURE 4
www.frontiersin.org

Figure 4. Scatter plot of the effect size for each SNP on tea consumption and albuminuria. (A) No SNP was excluded. (B) Two SNPs were removed (rs1030510 and rs4410790) due to potential horizontal pleiotropy.

Sensitivity Analyses

A GWAS involving CKD diagnosis in European men and women from the NHGRI-EBI GWAS Catalog6 was used to further validate the reliability of results. The MR analysis was conducted using SNP sites associated with CKD exposure in EBI, indicating the similar causal associations of the diagnosed CKD with IVW (OR = 0.803, p = 0.004). However, weighted median (OR = 0.824, p = 0.039) and MR Egger (OR = 0.758, p = 0.187) could not confirm a casual association (Supplementary Figure 7).

The Steiger-MR analysis was used to identify the presence of horizontal pleiotropy and the robustness of the causal effect estimates. There was no sign of heterogeneity, and Steiger-MR indicated that the SNPs explained more variance in exposure than the outcome (p > 0.05) (e.g., CKD, eGFR, and albuminuria).

Discussion

This two-sample MR study used summary-level data from UK Biobank and the CKDGen consortium to estimate the potential causal association of tea consumption with CKD, eGFR, and albuminuria. A total of 2,672 SNPs associated with tea consumption were found in UK Biobank, and 45 independent and available SNPs were found in CKDGen. The MR analyses showed that an increase in tea consumption appeared to be protective against CKD, eGFR, and albuminuria. The effects were generally similar in magnitude across diverse sensitivity analyses. Therefore, tea intake may be beneficial to renal function.

The association between tea intake and the risk of CKD remains inconsistent in previous longitudinal and cross-sectional studies. In an Iranian longitudinal study covering 1,780 adult Iranians (2006–2008 to 2012–2014), high tea intake was not associated with the risk of CKD (12). In the Doetinchem Cohort Study, including 4,722 individuals aged 26–65 years, tea consumption was not significantly associated with the changes in eGFR (14). However, in another Spanish cohort, including 5,851 overweight/obese elderly with MetS, overweight/obese adults with Mets, and high tea intake (>at least 1 cup/day) had a greater decline in eGFR in 1 year (13). Apart from this, the relationship between the intake of different types of tea (black tea, oolong tea, or green tea) and kidney function is also heterogenous. A Guangzhou Biobank cohort study including 12,428 elderly people did not find a significant association of green or black tea with eGFR. Nevertheless, a negative association between an intake of oolong tea and eGFR was found (15). The Singapore Chinese Health Study, a prospective cohort study of 63,257 participants aged 45–74 years, identified that tea intake is not associated with the risk of ESRD over an average follow-up of 17 years (16). Compared with previous studies, our study found that genetically predicted tea intake was causally associated with the decrease of CKD and albuminuria and the increase of eGFR. This MR study, in which confounding bias and reverse causality were avoided and a large number of individuals were from two-sample designs, is different from previous observational studies.

Several potential biological mechanisms may underlie an inverse association between tea intake and the risk of renal impairment. Growing evidence has shown that the generation of reactive oxygen species (ROS) plays an important role in the pathogenesis of kidney disease (42). Several interventional studies in humans have reported that the consumption of black and green tea improves vasodilator effects and decrease ROS concentrations in patients with renal failure (43, 44). Ardalan et al. reported that short-term intake of black tea could improve endothelial function and endothelium-dependent arterial vasodilation in renal transplant recipients (43). Hsu et al. revealed that the decaffeinated green tea extract (catechins) reduced hemodialysis-induced production of hydrogen peroxide and hypochlorous acid, the risk factors of atherosclerotic disease, and proinflammatory substances (44).

Previous studies have shown that tea catechins improve the metabolic mechanism of endothelial function. For example, the most abundant catechins in tea are epigallocatechin-3-gallate (EGCG), epigallocatechin, and epicatechin (45). These substances have been identified to be fairly strong antioxidants and free radical scavengers (46). Additionally, some murine and in vitro studies found the potential effects of EGCG on renal function. Yamabe et al. suggested that EGCG may reduce urinary protein excretion and serum glucose in mice with diabetic nephropathy (47). In rodent models with induced renal failure, tea catechins (including EGCG) have shown uniquely beneficial effects in decreasing nephrotoxicity (48, 49). The potential signal pathway mechanisms demonstrated that tea catechins could increase antioxidant activity, prostaglandin levels in the kidney, changes in transforming growth factor-β1 expression, and the regulation of nuclear factor-κβ (47, 48, 50, 51). Furthermore, compared with black tea, green tea has shown positive associations with renal function. The intake of green tea has been reported to reduce urinary oxalate excretion and deposition in rat models (52). In addition, green tea is very rich in catechins, including EGCG, which exerts an in vivo antiproliferative effect on renal cell carcinoma cell lines (53). Chinese green tea (CGT) has also been demonstrated for its protection against ROS, which causes apoptosis, inflammation, and damage in the lung tissue, among rats exposed to tobacco smoke for long periods; at the same time, CGT could also relieve nicotine toxicity, exerting antioxidant and anti-inflammatory properties (54, 55). Caffeine intake has been identified to exacerbate renal failure in a rat model (56). A cup of pure green tea usually contains around 25 mg of caffeine per 8 oz serving, which is considered to be a low amount (57). It is roughly one-fourth of the amount of caffeine in a typical cup of coffee and roughly half of the amount of caffeine in a typical cup of black tea (58).

The main advantage of this study is the two-sample method, which has a large summary-level genetic data that is able to avoid potential confounding factors and reverse causation in observational studies. The second advantage is that genetically predicted consumption of tea has been identified in a large GWAS of 349,376 European individuals, to diminish weak instrument bias (F-statistic > 10). However, this study has some limitations. First, the results did not apply to other populations on account of deviations from the data limited in European populations. Second, there were no non-linear relationships or stratification effects due to the summary-level data. Lastly, the type of tea and the amount of intake are important for exploring their causal association with the exposed genes. Further studies should take into consideration the overall impact of tea consumption.

Conclusion

It can be genetically predicted that there is a causal relationship of an extra cup of tea a day with the reduced risk of CKD and albuminuria and the increased level of eGFR. This MR study conducted a complete and detailed analysis of tea intake and renal function, which provided new evidence for the prevention and treatment of CKD.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

Ethical review and approval were waived for this study, all the data from Mendelian randomization is publicly accessible. Informed consent was obtained from all subjects in the original genome-wide association studies. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

YZ, YX, WW, TW, SS, JY, XW, LC, HZ, QY, and XL designed the study. YX analyzed the data. YZ wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Natural Science Foundation General Project of Chongqing Science and Technology Bureau (Grant No. cstc2020jcyj-msxmX0279) and the Science and Technology Research Project of Chongqing Education Commission (KJQN202100432).

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.

Publisher’s Note

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.

Supplementary Material

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

Footnotes

  1. ^ http://www.nealelab.is/uk-biobank
  2. ^ https://github.com/Nealelab/UK_Biobank_GWAS
  3. ^ http://gwas.mrcieu.au.uk
  4. ^ https://mrcieu.github.io/TwoSampleMR/articles/perform_mr.html#mr-methods-1
  5. ^ http://www.phenoscanner.medschl.cam.ac.uk/
  6. ^ https://www.ebi.ac.uk/gwas/downloads/summary-statistics

References

1. Shi Z, Taylor AW, Riley M, Byles J, Liu J, Noakes M. Association between dietary patterns, cadmium intake and chronic kidney disease among adults. Clin Nutr. (2018) 37:276–84. doi: 10.1016/j.clnu.2016.12.025

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Gbd Chronic Kidney Disease Collaboration. Global, regional, and national burden of chronic kidney disease, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. (2020) 395:709–33. doi: 10.1016/S0140-6736(20)30045-3

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Foreman KJ, Marquez N, Dolgert A, Fukutaki K, Fullman N, McGaughey M, et al. Forecasting life expectancy, years of life lost, and all-cause and cause-specific mortality for 250 causes of death: reference and alternative scenarios for 2016-40 for 195 countries and territories. Lancet. (2018) 392:2052–90. doi: 10.1016/S0140-6736(18)31694-5

CrossRef Full Text | Google Scholar

4. Webster AC, Nagler EV, Morton RL, Masson P. Chronic kidney disease. Lancet. (2017) 389:1238–52. doi: 10.1016/S0140-6736(16)32064-5

CrossRef Full Text | Google Scholar

5. Astor BC, Matsushita K, Gansevoort RT, van der Velde M, Woodward M, Levey AS, et al. Lower estimated glomerular filtration rate and higher albuminuria are associated with mortality and end-stage renal disease. A collaborative meta-analysis of kidney disease population cohorts. Kidney Int. (2011) 79:1331–40. doi: 10.1038/ki.2010.550

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Chronic Kidney Disease Prognosis Consortium, Matsushita K, van der Velde M, Astor BC, Woodward M, Levey AS, et al. Association of estimated glomerular filtration rate and albuminuria with all-cause and cardiovascular mortality in general population cohorts: a collaborative meta-analysis. Lancet. (2010) 375:2073–81. doi: 10.1016/S0140-6736(10)60674-5

CrossRef Full Text | Google Scholar

7. Xia J, Wang L, Ma Z, Zhong L, Wang Y, Gao Y, et al. Cigarette smoking and chronic kidney disease in the general population: a systematic review and meta-analysis of prospective cohort studies. Nephrol Dial Transplant. (2017) 32:475–87. doi: 10.1093/ndt/gfw452

PubMed Abstract | CrossRef Full Text | Google Scholar

8. White SL, Polkinghorne KR, Cass A, Shaw JE, Atkins RC, Chadban SJ. Alcohol consumption and 5-year onset of chronic kidney disease: the AusDiab study. Nephrol Dial Transplant. (2009) 24:2464–72. doi: 10.1093/ndt/gfp114

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zelle DM, Klaassen G, van Adrichem E, Bakker SJ, Corpeleijn E, Navis G. Physical inactivity: a risk factor and target for intervention in renal care. Nat Rev Nephrol. (2017) 13:152–68. doi: 10.1038/nrneph.2016.187

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Cosola C, Sabatino A, di Bari I, Fiaccadori E, Gesualdo L. Nutrients, nutraceuticals, and xenobiotics affecting renal health. Nutrients. (2018) 10:808. doi: 10.3390/nu10070808

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Cirillo M, Cavallo P, Palladino R, Terradura-Vagnarelli O, Zulli E, Villa R, et al. Relationship of the intake of water and other beverages with renal endpoints: cross-sectional and longitudinal data-observational, population-based study. J Ren Nutr. (2022) 32:68–77. doi: 10.1053/j.jrn.2021.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Gaeini Z, Bahadoran Z, Mirmiran P, Azizi F. Tea, coffee, caffeine intake and the risk of cardio-metabolic outcomes: findings from a population with low coffee and high tea consumption. Nutr Metab (Lond). (2019) 16:28. doi: 10.1186/s12986-019-0355-6

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Díaz-López A, Paz-Graniel I, Ruiz V, Toledo E, Becerra-Tomás N, Corella D, et al. Consumption of caffeinated beverages and kidney function decline in an elderly Mediterranean population with metabolic syndrome. Sci Rep. (2021) 11:8719. doi: 10.1038/s41598-021-88028-7

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Herber-Gast GC, van Essen H, Verschuren WM, Stehouwer CD, Gansevoort RT, Bakker SJ, et al. Coffee and tea consumption in relation to estimated glomerular filtration rate: results from the population-based longitudinal Doetinchem Cohort Study. Am J Clin Nutr. (2016) 103:1370–7. doi: 10.3945/ajcn.115.112755

PubMed Abstract | CrossRef Full Text | Google Scholar

15. van Hasselt TJ, Pickles O, Midgley-Hunt A, Jiang CQ, Zhang WS, Cheng KK, et al. Effects of tea consumption on renal function in a metropolitan Chinese population: the Guangzhou biobank cohort study. J Ren Nutr. (2014) 24:26–31. doi: 10.1053/j.jrn.2013.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lew QJ, Jafar TH, Jin A, Yuan JM, Koh WP. Consumption of coffee but not of other caffeine-containing beverages reduces the risk of end-stage renal disease in the singapore chinese health study. J Nutr. (2018) 148:1315–22. doi: 10.1093/jn/nxy075

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Stephen B, Thompson SG. Mendelian Randomization: Methods for Using Genetic Variants in Causal Estimation. London: Chapman and Hall/CRC (2015).

Google Scholar

18. Larsson SC. Mendelian randomization as a tool for causal inference in human nutrition and metabolism. Curr Opin Lipidol. (2021) 32:1–8. doi: 10.1097/MOL.0000000000000721

PubMed Abstract | CrossRef Full Text | Google Scholar

19. 1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, et al. A global reference for human genetic variation. Nature. (2015) 526:68–74. doi: 10.1038/nature15393

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Shim H, Chasman DI, Smith JD, Mora S, Ridker PM, Nickerson DA, et al. A multivariate genome-wide association analysis of 10 LDL subfractions, and their response to statin treatment, in 1868 Caucasians. PLoS One. (2015) 10:e0120758. doi: 10.1371/journal.pone.0120758

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol. (2011) 40:740–52. doi: 10.1093/ije/dyq151

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Pattaro C, Teumer A, Gorski M, Chu AY, Li M, Mijatovic V, et al. Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun. (2016) 7:10023. doi: 10.1038/ncomms10023

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Chen X, Kong J, Pan J, Huang K, Zhou W, Diao X, et al. Kidney damage causally affects the brain cortical structure: a Mendelian randomization study. EBioMedicine. (2021) 72:103592. doi: 10.1016/j.ebiom.2021.103592

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Schwartz GJ, Schneider MF, Maier PS, Moxey-Mims M, Dharnidharka VR, Warady BA, et al. Improved equations estimating GFR in children with chronic kidney disease using an immunonephelometric determination of cystatin C. Kidney Int. (2012) 82:445–53. doi: 10.1038/ki.2012.169

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF III, Feldman HI, et al. CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration). A new equation to estimate glomerular filtration rate. Ann Intern Med. (2009) 150:604–12. doi: 10.7326/0003-4819-150-9-200905050-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Wuttke M, Li Y, Li M, Sieber KB, Feitosa MF, Gorski M, et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet. (2019) 51:957–72. doi: 10.1038/s41588-019-0407-x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Teumer A, Tin A, Sorice R, Gorski M, Yeo NC, Chu AY, et al. Genome-wide association studies identify genetic loci associated with albuminuria in diabetes. Diabetes. (2016) 65:803–17. doi: 10.2337/db15-1313

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Warram JH, Gearin G, Laffel L, Krolewski AS. Effect of duration of type I diabetes on the prevalence of stages of diabetic nephropathy defined by urinary albumin/creatinine ratio. J Am Soc Nephrol. (1996) 7:930–7. doi: 10.1681/ASN.V76930

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Guillán-Fresco M, Franco-Trepat E, Alonso-Pérez A, Jorge-Mora A, López-Fagúndez M, Pazos-Pérez A, et al. Caffeine, a risk factor for osteoarthritis and longitudinal bone growth inhibition. J Clin Med. (2020) 9:1163. doi: 10.3390/jcm9041163

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Burgess STS. Mendelian Randomization: Methods for Using Genetic Variants in Causal Estimation. 1st ed. Boca Raton, FL: CRC Press (2015).

Google Scholar

31. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Morrison J, Knoblauch N, Marcus JH, Stephens M, He X. Mendelian randomization accounting for correlated and uncorrelated pleiotropic effects using genome-wide summary statistics [published correction appears in Nat Genet. 2020 May 29;:]. Nat Genet. (2020) 52:740–7. doi: 10.1038/s41588-020-0631-4

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Didelez V, Sheehan N. Mendelian randomization as an instrumental variable approach to causal inference. Stat Methods Med Res. (2007) 16:309–30. doi: 10.1177/0962280206077743

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Greenland S. An introduction to instrumental variables for epidemiologists. Int J Epidemiol. (2018) 47:358. doi: 10.1093/ije/dyx275

PubMed Abstract | CrossRef Full Text | Google Scholar

35. 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. doi: 10.1002/sim.6522

PubMed Abstract | CrossRef Full Text | Google Scholar

36. 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. doi: 10.1093/ije/dyv080

PubMed Abstract | CrossRef Full Text | Google Scholar

37. 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. doi: 10.1002/gepi.21965

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. (2017) 46:1985–98. doi: 10.1093/ije/dyx102

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Verbanck M, Chen CY, 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. doi: 10.1038/s41588-018-0099-7

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan NA, Thompson JR. Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic. Int J Epidemiol. (2016) 45:1961–74. doi: 10.1093/ije/dyw220

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. (2017) 13:e1007081. doi: 10.1371/journal.pgen.1007081

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Agharazii M, St-Louis R, Gautier-Bastien A, Ung RV, Mokas S, Larivière R, et al. Inflammatory cytokines and reactive oxygen species as mediators of chronic kidney disease-related vascular calcification. Am J Hypertens. (2015) 28:746–55. doi: 10.1093/ajh/hpu225

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Ardalan MR, Tarzamni MK, Shoja MM, Tubbs RS, Rahimi-Ardabili B, Ghabili K, et al. Black tea improves endothelial function in renal transplant recipients. Transplant Proc. (2007) 39:1139–42. doi: 10.1016/j.transproceed.2007.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Hsu SP, Wu MS, Yang CC, Huang KC, Liou SY, Hsu SM, et al. Chronic green tea extract supplementation reduces hemodialysis-enhanced production of hydrogen peroxide and hypochlorous acid, atherosclerotic factors, and proinflammatory cytokines. Am J Clin Nutr. (2007) 86:1539–47. doi: 10.1093/ajcn/86.5.1539

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lee MJ, Maliakal P, Chen L, Meng X, Bondoc FY, Prabhu S, et al. Pharmacokinetics of tea catechins after ingestion of green tea and (-)-epigallocatechin-3-gallate by humans: formation of different metabolites and individual variability. Cancer Epidemiol Biomarkers Prev. (2002) 11(10Pt 1):1025–32.

Google Scholar

46. Henning SM, Niu Y, Liu Y, Lee NH, Hara Y, Thames GD, et al. Bioavailability and antioxidant effect of epigallocatechin gallate administered in purified form versus as green tea extract in healthy individuals. J Nutr Biochem. (2005) 16:610–6. doi: 10.1016/j.jnutbio.2005.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Yamabe N, Yokozawa T, Oya T, Kim M. Therapeutic potential of (-)-epigallocatechin 3-O-gallate on renal damage in diabetic nephropathy model rats. J Pharmacol Exp Ther. (2006) 319:228–36. doi: 10.1124/jpet.106.107029

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Sahin K, Tuzcu M, Gencoglu H, Dogukan A, Timurkan M, Sahin N, et al. Epigallocatechin-3-gallate activates Nrf2/HO-1 signaling pathway in cisplatin-induced nephrotoxicity in rats. Life Sci. (2010) 87:240–5. doi: 10.1016/j.lfs.2010.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Yokozawa T, Chung HY, He LQ, Oura H. Effectiveness of green tea tannin on rats with chronic renal failure. Biosci Biotechnol Biochem. (1996) 60:1000–5. doi: 10.1271/bbb.60.1000

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Rhee SJ, Kim MJ, Kwag OG. Effects of green tea catechin on prostaglandin synthesis of renal glomerular and renal dysfunction in streptozotocin-induced diabetic rats. Asia Pac J Clin Nutr. (2002) 11:232–6. doi: 10.1046/j.1440-6047.2002.00312.x

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Shi SH, Zheng SS, Jia CK, Zhu YF, Xie HY. Inhibitory effect of tea polyphenols on transforming growth factor-beta1 expression in rat with cyclosporine A-induced chronic nephrotoxicity. Acta Pharmacol Sin. (2004) 25:98–103.

Google Scholar

52. Itoh Y, Yasui T, Okada A, Tozawa K, Hayashi Y, Kohri K. Preventive effects of green tea on renal stone formation and the role of oxidative stress in nephrolithiasis. J Urol. (2005) 173:271–5. doi: 10.1097/01.ju.0000141311.51003.87

CrossRef Full Text | Google Scholar

53. Carvalho M, Jeronimo C, Valentão P, Andrade P, Silva BM. Green tea: a promising anticancer agent for renal cell carcinoma. Food Chem. (2010) 122:49–54. doi: 10.1016/j.foodchem.2010.02.014

CrossRef Full Text | Google Scholar

54. Al-Awaida W, Akash M, Aburubaiha Z, Talib WH, Shehadeh H. Chinese green tea consumption reduces oxidative stress, inflammation and tissues damage in smoke exposed rats. Iran J Basic Med Sci. (2014) 17:740–6.

Google Scholar

55. Al-Awaida WJ, Zihlif MA, Al-Ameer HJ, Sharab A, Akash M, Aburubaiha ZA, et al. The effect of green tea consumption on the expression of antioxidant- and inflammation-related genes induced by nicotine. J Food Biochem. (2019) 43:e12874. doi: 10.1111/jfbc.1287

CrossRef Full Text | Google Scholar

56. Tofovic SP, Kost CK Jr., Jackson EK, Bastacky SI. Long-term caffeine consumption exacerbates renal failure in obese, diabetic, ZSF1 (fa-fa(cp)) rats. Kidney Int. (2002) 61:1433–44. doi: 10.1046/j.1523-1755.2002.00278.x

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Unno K, Noda S, Kawasaki Y, Yamada H, Morita A, Iguchi K, et al. Ingestion of green tea with lowered caffeine improves sleep quality of the elderly via suppression of stress. J Clin Biochem Nutr. (2017) 61:210–6. doi: 10.3164/jcbn.17-6

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Kole J, Barnhill A. Caffeine content labeling: a missed opportunity for promoting personal and public health. J Caffeine Res. (2013) 3:108–13. doi: 10.1089/jcr.2013.0017

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: tea consumption, chronic kidney disease (CKD), SNPs, albuminuria, Mendelian randomization

Citation: Zhang Y, Xiong Y, Shen S, Yang J, Wang W, Wu T, Chen L, Yu Q, Zuo H, Wang X and Lei X (2022) Causal Association Between Tea Consumption and Kidney Function: A Mendelian Randomization Study. Front. Nutr. 9:801591. doi: 10.3389/fnut.2022.801591

Received: 01 November 2021; Accepted: 07 February 2022;
Published: 29 March 2022.

Edited by:

Raul Zamora-Ros, Institut d’Investigacio Biomedica de Bellvitge (IDIBELL), Spain

Reviewed by:

Youxin Wang, Capital Medical University, China
Wajdy Awaida, American University of Madaba, Jordan

Copyright © 2022 Zhang, Xiong, Shen, Yang, Wang, Wu, Chen, Yu, Zuo, Wang and Lei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xun Lei, leixun521@163.com

These authors have contributed equally to this work

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.