Gang Li
Yang Liu- Department of Sports Medicine, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, Xinjiang, China
Background: This study aimed to clarify the relationship between tea consumption and osteoarthritis (OA).
Methods: Common single-nucleotide polymorphisms (SNPs) from the Open Genome-wide Association Studies database were obtained. Summary statistics on OA were retrieved from the second dataset that enrolled 50,508 participants (10,083 OA cases) of European ancestry. The causal association between tea intake and OA was tested using two-sample Mendelian randomization (MR) analysis.
Results: Tea consumption has adverse effects on OA. (inverse-variance weighted method: OR = 1.19, 95% CI = 1.08–1.30; weighted median method: OR = 1.22, 95% CI = 1.07–1.40). The MR–Egger regression intercept (MR intercept = −0.002; p = 0.73) showed no evidence of directional pleiotropy. Moreover, no evidence of underlying heterogeneity in MR analysis was found according to Cochran’s Q test and funnel and forest analyses.
Conclusion: A genetically predicted high daily tea intake can increase the risk of OA.
Background
Osteoarthritis (OA) is the most prevalent type of arthritis and is characterized by articular cartilage degeneration, remodeling of underlying bone, and synovitis (Loeser et al., 2012; Martel-Pelletier et al., 2016). According to previous studies of 291 conditions, hip and knee OA is the 11th highest disease in terms of global disability and 38th highest in disability-adjusted life years (Cross et al., 2014). Epidemiological studies from Prieto-Alhambra et al. (2014) showed that OA is prevalent among aged individuals and severely affects the quality of life. Patients with OA suffer higher risks of disability and all-cause mortality than their healthy counterparts. These burdens motivate researchers to explore underlying mechanisms and develop effective methods against OA (Palazzo et al., 2016).
Accumulating lines of evidence show that OA is a multifactorial disorder of articular cartilage. Aging, obesity, inflammation, trauma, and overloading are closely associated with OA (Kulkarni et al., 2016; Rahmati et al., 2017; Nickien et al., 2018; Berenbaum and Walker, 2020; Wang et al., 2020). Therapies that suppress inflammatory mediators, such as cyclooxygenase-2, and oxidative stress and enhance autophagy, have shown efficacy in animal models and clinical settings to some extent. The elderly may benefit from complementary and alternative medicine, such as tea (Lapane et al., 2013). However, the association and underlying mechanism between tea intake and OA remains inconsistent and needs further investigation.
Tea is a popular drink worldwide. Tea contains chemicals, such as caffeic acid, caffeine, catechin, coumaric acid, and gallic acid, which have effects on OA (Luk et al., 2020). However, few studies have explored the relationship between tea intake and OA in humans, and the findings are inconsistent. In a 5-years cohort study, Takiguchi et al. (2019) reported that low intake of green tea is associated with a high incidence of knee OA in males but not in females. A randomized clinical study revealed that green tea intake could not decrease the pain and stiffness of patients with OA as assessed by the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) questionnaire (Hashempur et al., 2018). However, in animal models, tea intake can decrease cartilage degeneration, inhibit inflammation, and function as a protective factor. Further studies are needed to explore their relationship (Luk et al., 2020).
Summary-level datasets were used to verify the relationship between tea intake and OA. Mendelian randomization (MR) was adopted to test their causal links.
Materials and methods
Data source for genetic variants
The MRC IEU Open GWAS Project (https://gwas.mrcieu.ac.uk/) is an online database for searching GWAS datasets and traits. SNP information of tea intake was extracted from Neale Lab (http://www.nealelab.is/uk-biobank, GWAS round 2), which consists of more than 3,49,376 samples from European ancestry (Wang et al., 2021). Based on a dietary questionnaire, daily tea consumption was determined at baseline. The touchscreen question that needs to be answered is “How many cups of tea do you drink every day?” (including black and green tea). A total of 2,672 unique single-nucleotide polymorphisms (SNPs) were chosen based on minor frequency >1% of SNPs with strong correlation as p < 5 × 108. We clumped the 2,672 SNPs at a 10,000 kb window and linkage disequilibrium with r2 < 0.001 by using the 1,000 Genomes Project reference for the European samples. The results confirmed the genetic independence among exposure variants. A total of 45 significant SNPs related to exposure were identified. Large-scale GWAS was carried out from the UK Biobank, which included 50,508 European ancestry, and yielded summary-level data with a clinical diagnosis for OA (10,083 cases and 40,425 controls). Protocols related to the data have been released and described in previous studies (Wu et al., 2013; Okada et al., 2014; Zengini et al., 2018). Summary data can be obtained from the UK MRC IEU Open GWAS Project database (http://gwas.mrcieu.au.uk). Informed consent was obtained from all participants. All GWAS-related current analyses were approved by the relevant ethics committees.
Statistical analysis
MR analysis requires genetically predicted SNP exposure that is not associated with potential confounders (Burgess et al., 2013). In the first step, we selected 45 SNPs that were independently associated with tea consumption. In the second step, we identified the association of each SNP with OA risk. In the third step, MR was used to estimate the uncompounded causal relationship between exposure and outcome risk. Based on summary statistics from two different GWASs, we performed a two-sample MR to estimate the causal effect of exposure on outcomes. We then assessed the causal relation by using GWAS datasets from tea intake and OA with the independent 45 SNPs as instruments (Hartwig et al., 2016). An inverse-variance weighted (IVW) approach based on combining Wald ratio estimates from different SNPs was employed to estimate the causal effect and provide consistent estimates of the effect of genetic variants on outcome when the instrument’s assumptions are verified (Pierce and Burgess, 2013). Although multiple variants in MR analysis increased the statistical power, they might not be valid instrumental variants due to their pleiotropic nature (Hartwig et al., 2016). Therefore, MR–Egger regression and weighted median estimators were used to avoid pleiotropy (multiple genetic variants associated with multiple variables). MR and Egger are weighted regression models that introduce an intercept to account for directional pleiotropy. A horizontal pleiotropy occurs when intercept values differ from zero (Bowden et al., 2015). Weighted linear regression was conducted using MR–Egger to estimate the coefficients of gene exposure and outcome (Bowden et al., 2015). In this method, unbiased estimates are used when pleiotropic instruments are used, in which the pleiotropic effects are unrelated to instrument size (Burgess and Thompson, 2017). Weighted median estimation supports a robust evaluation of causal power, even when genetic variants are invalid instrument factors contributing up to 50% of the explanation (Bowden et al., 2016). In comparison with MR–Egger, weighted median approach maintains greater precision in the estimates. Additionally, mode-based causal estimation consistently estimates true causal effects when most instruments generate consistent MR estimates.
Heterogeneities between SNP exposure were examined by Cochran’s Q statistics and funnel plots, respectively (Egger et al., 1997). Leave-one-out was used to visualize whether the causal relation is influenced by one outlier SNP alone. We also performed Wald ratio estimates on each SNP associated with OA.
All MR analyses were performed using Two-sample MR packages, and statistical significance was tested at p < 0.05 (Hemani et al., 2018).
Results
Instrumental variables for Mendelian randomization
Forty-five independent SNPs related to genetic tea intake were determined using R2 statistic. F-statistics were achieved to assess the strength of exposure-related instruments (Pierce et al., 2011; Shim et al., 2015). Table 1 provides detailed information about SNPs and exposures associated with selected SNPs.
TABLE 1. Characteristics of tea consumption-associated SNPs.
Mendelian randomization between tea intake and OA
Tea consumption was positively associated with OA according to the IVW approach [odds ratio (OR) = 1.22; 95% confidence interval (CI): 1.08, 1.30; Figures 1, 2]. An intercept represents the average pleiotropic effect across genetic variants (that is, how a variant affects the outcome on average). MR–Egger regression found that directional pleiotropy was unlikely to bias the results. If the intercept is different from zero (MR–Egger test), then directional pleiotropy exists (MR–Egger intercept = −0.002; p = 0.73). Tea consumption was not related to OA via MR–Egger (OR = 1.22; 95% CI: 0.98, 1.53; Figures 1, 2). However, weighted median (OR = 1.22; 95% CI: 1.07, 1.40) and weighted mode (OR = 1.22; 95% CI: 1.02, 1.45) suggested a causal relation between them (Table 2; Figure 2). As a result, IVW, weighted median, modes, and MR–Egger methods did not present consistent results regarding tea consumption and OA. Compared with the MR–Egger analysis, the weighted median estimator maintains greater precision in estimates, which may suggest a causal association between high tea intake and increased risk of OA (Bowden et al., 2015).
FIGURE 1. Forest plots of MR study using genetically predicted tea consumption with OA. Inverse-variance weighted methods, MR–Egger analysis, weighted median, and simple and weighted modes were used in this study.
FIGURE 2. Scatter plot of the effect size for each SNP on tea consumption and OA.
TABLE 2. MR estimates from each method of the causal effect of tea consumption on osteoarthritis risk.
Heterogeneity and sensitivity tests
Cochran’s Q test and funnel plot did not indicate heterogeneity between SNP estimates based on single variants (Table 2; Figure 3). In the leave-one-out analysis, no SNP was found to be a potential outlier that could influence point estimates (Figure 4). In addition, the forest plot of each SNP-associated tea intake on OA showed that the majority of SNPs were in accordance with the same direction (Figure 5).
FIGURE 3. Funnel plot for the causal estimates of each tea SNP with OA.
FIGURE 4. Leave-one-out sensitivity analysis for the causal association between tea consumption and OA.
FIGURE 5. Forest plot of the effect size for each SNP on tea consumption and OA.
Discussion
In recent years, aging has become increasingly pronounced globally. As a condition closely related to aging, the prevalence of OA faces an upward trend, placing heavy pressure on the healthcare system and patients’ quality of life. Tea intake may be associated with the incidence of OA, but the conclusions from animal models and clinical studies vary (Luk et al., 2020). Our study is the first to adopt MR to examine the causal link between tea intake and OA to provide new evidence in clinical settings.
Our results showed that tea intake increases the OA risk, which is inconsistent with previous reports. Epigallocatechin gallate (EGCG) is the major bioactive component of polyphenolic fraction in green tea. In cellular studies, Huang et al. (2015) used EGCG to treat primary rabbit articular chondrocytes and found that the EGCG-treated group had higher amounts of cartilage extracellular matrix and collagen II synthesized than the control group. Thus, EGCG may protect the cartilage from degradation and prevent OA. Similar results were reported by Elder et al. (2017) and Bae et al. (2010). This protective phenomenon is also replicated in animal models. Haqqi et al. (1999) reported that tea extracts and polyphenolic fraction could alleviate collagen-induced arthritis in mice. However, few studies were performed in humans, and their conclusions vary. In an open-label clinical study with 40 knee OA cases, patients treated with green tea extract for 4 weeks showed alleviated pain as assessed by the Visual Analog Scale but had the same extent of stiffness as assessed by WOMAC (Hashempur et al., 2018). Among the 40 participants, 80% were males. Considering the small number and uneven baseline population, the robustness of results may be careful. As a covariate in another study that enrolled 11,091 participants, the protective effects of green tea intake were observed in males but not in females (Takiguchi et al., 2019). In a further subgroup analysis, the protective effects disappeared in males aged more than 60 years, who are the high-risk population of OA. This preliminary finding was deduced based on observational study, which could not overcome endogeneity. By contrast, Chet et al. determined no causal association between tea intake and OA (knee OA and hip OA) by using MR method (Chen et al., 2022). However, in their studies, the odds ratios were greater than 1.00 with marginal significance (1.11 for knee OA; 1.20 for hip OA), indicating the likelihood of less statistical power. Therefore, further studies with more statistical power are required.
The mechanisms of the relation of tea to OA remain unknown. In previous studies, EGCG is used to imitate the function of tea and further explore the specific molecular mechanisms in cellular studies and animal models (Leong et al., 2014; Rasheed et al., 2016; Jin et al., 2020). Leong et al. (Leong et al., 2014) found that EGCG-treated mice (25 mg kg−1 per day) had reduced expression of pro-inflammatory cytokines, namely, IL-1β and TNF-α. Pain symptoms were also alleviated compared with those in the control group. However, EGCG cannot be used as a substitute of tea in human study. In addition to EGCG, other major components, such as caffeine, may influence OA. In a recent review, Guillán comprehensively summarized the evidence regarding the relationship between caffeine and OA and concluded that caffeine increases the OA risk (Guillan-Fresco et al., 2020). The influence of other compounds may lead to different conclusion between our study and others.
MR is an effective method used to evaluate causal links because it can overcome endogeneity in observational studies. It can also be used to test the relationships between other risk factors, such as caffeine and OA (Lee, 2018). To the best of our knowledge, this study is the first to use MR to investigate the relationship between tea intake and OA. However, this study has some limitations. First, the SNPs used are from the European population, which may lead to bias. Further validations in different races should be made. Moreover, this study used green tea and black tea, whose chemical components vary. Consequently, their effects also differ. The exposure and outcome examples were obtained from UK Biobank; as such, sample overlap could lead to substantial bias and inflated Type 1 error (Burgess, Davies, Thompson). Therefore, sensitivity analyses should be conducted using fewer genetic variants with stronger power, and two-sample method can be performed in further causal research.
Currently, tea intake is prevalent among aged individuals (Tang et al., 2019). Given that increased risks were reported in our study, tea should be carefully used as a common supplement and alternative medicine.
Conclusion
New causal evidence is provided for the association of genetically predicted tea consumption with increased risks of OA.
Data availability statement
Publicly available datasets were analyzed in this study. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Author contributions
GL writes the manuscript; ZZ and YL revise and review the manuscript; YL drafts the study design and supervises all process.
Acknowledgments
Thanks to the Open GWAS team for collecting and sharing GAWS datasets.
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.
References
Bae, J. Y., Han, D. W., Wakitani, S., Nawata, M., and Hyon, S. H. (2010). Biological and biomechanical evaluations of osteochondral allografts preserved in cold storage solution containing epigallocatechin gallate. Cell Transpl. 19, 681–689. doi:10.3727/096368910X508771
Berenbaum, F., and Walker, C. (2020). Osteoarthritis and inflammation: A serious disease with overlapping phenotypic patterns. Postgrad. Med. 132, 377–384. doi:10.1080/00325481.2020.1730669
Bowden, J., Smith, G. D., and Burgess, S. (2015). Mendelian randomization with invalid instruments: Effect estimation and bias detection through egger regression. Int. J. Epidemiol. 44, 512–525. doi:10.1093/ije/dyv080
Bowden, J., Smith, G. D., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314. doi:10.1002/gepi.21965
Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37, 658–665. doi:10.1002/gepi.21758
Burgess, S., Davies, N. M., and Thompson, S. G. (201). Bias due to participant overlap in two-sample Mendelian randomization. Genet. Epidemiol. 40 (7), 597–608. doi:10.1002/gepi.21998
Burgess, S., and Thompson, S. G. (2017). Interpreting findings from mendelian randomization using the mr-egger method. Eur. J. Epidemiol. 32, 377–389. doi:10.1007/s10654-017-0255-x
Chen, S., Chen, T., Chen, Y., Huang, D., Pan, Y., and Chen, S. (2022). Causal association between tea consumption and bone health: A mendelian randomization study. Front. Nutr. 9, 872451. doi:10.3389/fnut.2022.872451
Cross, M., Smith, E., Hoy, D., Nolte, S., Ackerman, I., Fransen, M., et al. (2014). The global burden of hip and knee osteoarthritis: Estimates from the global burden of disease 2010 study. Ann. Rheum. Dis. 73, 1323–1330. doi:10.1136/annrheumdis-2013-204763
Egger, M., Smith, G. D., and Phillips, A. N. (1997). Meta-analysis: Principles and procedures. Bmj-British Med. J. 315, 1533–1537. doi:10.1136/bmj.315.7121.1533
Elder, S., Clune, J., Walker, J., and Gloth, P. (2017). Suitability of egcg as a means of stabilizing a porcine osteochondral xenograft. J. Funct. Biomater. 8, E43. doi:10.3390/jfb8040043
Guillan-Fresco, M., Franco-Trepat, E., Alonso-Perez, A., Jorge-Mora, A., Lopez-Fagundez, M., Pazos-Perez, A., et al. (2020). Caffeine, a risk factor for osteoarthritis and longitudinal bone growth inhibition. J. Clin. Med. 9, E1163. doi:10.3390/jcm9041163
Haqqi, T. M., Anthony, D. D., Gupta, S., Ahmad, N., Lee, M. S., Kumar, G. K., et al. (1999). Prevention of collagen-induced arthritis in mice by a polyphenolic fraction from green tea. Proc. Natl. Acad. Sci. U. S. A. 96, 4524–4529. doi:10.1073/pnas.96.8.4524
Hartwig, F. P., Davies, N. M., Hemani, G., and Davey Smith, G. (2016). Two-sample mendelian randomization: Avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int. J. Epidemiol. 45, 1717–1726. doi:10.1093/ije/dyx028
Hashempur, M. H., Sadrneshin, S., Mosavat, S. H., and Ashraf, A. (2018). Green tea (camellia sinensis) for patients with knee osteoarthritis: A randomized open-label active-controlled clinical trial. Clin. Nutr. 37, 85–90. doi:10.1016/j.clnu.2016.12.004
Hemani, G., Zheng, J., Wade, K. H., Laurin, C., Elsworth, B., Burgess, S., et al. (2018). Mr-base: A platform for systematic causal inference across the phenome using billions of genetic associations. Eur. J. Hum. Genet. 26, 702–703.
Huang, H. J., Liu, Q., Liu, L., Wu, H. Y., and Zheng, L. (2015). Effect of epigallocatechin-3-gallate on proliferation and phenotype maintenance in rabbit articular chondrocytes in vitro. Exp. Ther. Med. 9, 213–218. doi:10.3892/etm.2014.2057
Jin, Y. J., Koh, R. H., Kim, S. H., Kim, K. M., Park, G. K., and Hwang, N. S. (2020). Injectable anti-inflammatory hyaluronic acid hydrogel for osteoarthritic cartilage repair. Mat. Sci. Eng. C Mat. Biol. Appl., 111096. doi:10.1016/j.msec.2020.111096
Kulkarni, K., Karssiens, T., Kumar, V., and Pandit, H. (2016). Obesity and osteoarthritis. Maturitas 89, 22–28. doi:10.1016/j.maturitas.2016.04.006
Lapane, K. L., Yang, S., Jawahar, R., McAlindon, T., and Eaton, C. B. (2013). CAM use among overweight and obese persons with radiographic knee osteoarthritis. BMC Complement. Altern. Med. 13, 241. doi:10.1186/1472-6882-13-241
Lee, Y. H. (2018). Investigating the possible causal association of coffee consumption with osteoarthritis risk using a mendelian randomization analysis. Clin. Rheumatol. 37, 3133–3139. doi:10.1007/s10067-018-4252-6
Leong, D. J., Choudhury, M., Hanstein, R., Hirsh, D. M., Kim, S. J., Majeska, R. J., et al. (2014). Green tea polyphenol treatment is chondroprotective, anti-inflammatory and palliative in a mouse posttraumatic osteoarthritis model. Arthritis Res. Ther. 16, 508. doi:10.1186/s13075-014-0508-y
Loeser, R. F., Goldring, S. R., Scanzello, C. R., and Goldring, M. B. (2012). Osteoarthritis: A disease of the joint as an organ. Arthritis Rheum. 64, 1697–1707. doi:10.1002/art.34453
Luk, H. Y., Appell, C., Chyu, M. C., Chen, C. H., Wang, C. Y., Yang, R. S., et al. (2020). Impacts of green tea on joint and skeletal muscle health: Prospects of translational nutrition. Antioxidants 9, E1050. doi:10.3390/antiox9111050
Martel-Pelletier, J., Barr, A. J., Cicuttini, F. M., Conaghan, P. G., Cooper, C., Goldring, M. B., et al. (2016). Osteoarthritis. Nat. Rev. Dis. Prim. 2, 16072. doi:10.1038/nrdp.2016.72
Nickien, M., Heuijerjans, A., Ito, K., and Van Donkelaar, C. C. (2018). Comparison between in vitro and in vivo cartilage overloading studies based on a systematic literature review. J. Orthop. Res. 36, 2076–2086. doi:10.1002/jor.23910
Okada, Y., Wu, D., Trynka, G., Raj, T., Terao, C., Ikari, K., et al. (2014). Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature 506, 376–381. doi:10.1038/nature12873
Palazzo, C., Nguyen, C., Lefevre-Colau, M. M., Rannou, F., and Poiraudeau, S. (2016). Risk factors and burden of osteoarthritis. Ann. Phys. Rehabil. Med. 59, 134–138. doi:10.1016/j.rehab.2016.01.006
Pierce, B. L., Ahsan, H., and Vanderweele, T. J. (2011). Power and instrument strength requirements for mendelian randomization studies using multiple genetic variants. Int. J. Epidemiol. 40, 740–752. doi:10.1093/ije/dyq151
Pierce, B. L., and Burgess, S. (2013). Efficient design for mendelian randomization studies: Subsample and 2-sample instrumental variable estimators. Am. J. Epidemiol. 178, 1177–1184. doi:10.1093/aje/kwt084
Prieto-Alhambra, D., Judge, A., Javaid, M. K., Cooper, C., Diez-Perez, A., and Arden, N. K. (2014). Incidence and risk factors for clinically diagnosed knee, hip and hand osteoarthritis: Influences of age, gender and osteoarthritis affecting other joints. Ann. Rheum. Dis. 73, 1659–1664. doi:10.1136/annrheumdis-2013-203355
Rahmati, M., Nalesso, G., Mobasheri, A., and Mozafari, M. (2017). Aging and osteoarthritis: Central role of the extracellular matrix. Ageing Res. Rev. 40, 20–30. doi:10.1016/j.arr.2017.07.004
Rasheed, Z., Rasheed, N., and Al-Shobaili, H. A. (2016). Epigallocatechin-3-o-gallate up-regulates microrna-199a-3p expression by down-regulating the expression of cyclooxygenase-2 in stimulated human osteoarthritis chondrocytes. J. Cell. Mol. Med. 20, 2241–2248. doi:10.1111/jcmm.12897
Shim, H., Chasman, D. I., Smith, J. D., Mora, S., Ridker, P. M., Nickerson, D. A., et al. (2015). A multivariate genome-wide association analysis of 10 ldl subfractions, and their response to statin treatment, in 1868 caucasians. Plos One 10, e0120758. doi:10.1371/journal.pone.0120758
Takiguchi, R., Komatsu, R., Kitamura, K., Watanabe, Y., Takahashi, A., Kobayashi, R., et al. (2019). Modifiable factors associated with symptomatic knee osteoarthritis: The murakami cohort study. Maturitas 128, 53–59. doi:10.1016/j.maturitas.2019.06.013
Tang, G. Y., Meng, X., Gan, R. Y., Zhao, C. N., Liu, Q., Feng, Y. B., et al. (2019). Health functions and related molecular mechanisms of tea components: An update review. Int. J. Mol. Sci. 20, E6196. doi:10.3390/ijms20246196
Wang, L. J., Zeng, N., Yan, Z. P., Li, J. T., and Ni, G. X. (2020). Post-traumatic osteoarthritis following acl injury. Arthritis Res. Ther. 22, 57. doi:10.1186/s13075-020-02156-5
Wang, M. M., Bai, Y. J., Wang, Z., Zhang, Z. Z., Liu, D. D., and Lian, X. G. (2021). Higher tea consumption is associated with decreased risk of small vessel stroke. Clin. Nutr. 40, 1430–1435. doi:10.1016/j.clnu.2020.08.039
Wu, J. H. Y., Lemaitre, R. N., Manichaikul, A., Guan, W. H., Tanaka, T., Foy, M., et al. (2013). Genome-wide association study identifies novel loci associated with concentrations of four plasma phospholipid fatty acids in the de novo lipogenesis pathway results from the cohorts for heart and aging research in genomic epidemiology (charge) consortium. Circ. Cardiovasc. Genet. 6, 171–183. doi:10.1161/CIRCGENETICS.112.964619
Keywords: tea, osteoarthritis, Mendelian randomization, risk factor, causal association
Citation: Li G, Zhang Z and Liu Y (2022) Genetically predicted tea intake increases the risk of osteoarthritis: A Mendelian randomization study. Front. Genet. 13:1004392. doi: 10.3389/fgene.2022.1004392
Received: 27 July 2022; Accepted: 21 September 2022;
Published: 04 October 2022.
Edited by:
Xi Wang, Xi’an Jiaotong University, ChinaReviewed by:
Chengping Wen, Zhejiang Chinese Medical University, ChinaMasahiro Yoshikawa, Nihon University School of Medicine, Japan
Copyright © 2022 Li, Zhang and Liu. 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: Yang Liu, 8861450@qq.com