Tea consumption and risk of lower respiratory tract infections: a two-sample mendelian randomization study.

BACKGROUND: Observational studies have reported the association between tea consumption and the risk of lower respiratory tract infections (LRTIs). However, a consensus has yet to be reached, and whether the observed association is driven by confounding factors or reverse causality remains unclear.
METHOD: A two-sample Mendelian randomization (MR) analysis was conducted to determine whether genetically predicted tea intake is causally associated with the risk of common LRTI subtypes. Genome-wide association study (GWAS) from UK Biobank was used to identify single-nucleotide polymorphisms (SNPs) associated with an extra cup of tea intake each day. The summary statistics for acute bronchitis, acute bronchiolitis, bronchiectasis, pneumonia, and influenza and pneumonia were derived from the FinnGen project.
RESULTS: We found that genetically predicted an extra daily cup of tea intake was causally associated with the decreased risk of bronchiectasis [odds ratio (OR) = 0.61, 95% confidence interval (CI) = 0.47-0.78, P < 0.001], pneumonia (OR = 0.90, 95% CI = 0.85-0.96, P = 0.002), influenza and pneumonia (OR = 0.91, 95% CI = 0.85-0.97, P = 0.002), but not with acute bronchitis (OR = 0.91, 95% CI = 0.82-1.01, P = 0.067) and acute bronchiolitis (OR = 0.79, 95% CI = 0.60-1.05, P = 0.100). Sensitivity analyses showed that no heterogeneity and pleiotropy could bias the results.
CONCLUSIONS: Our findings provided new evidence that genetically predicted an extra daily cup of tea intake may causally associated with a decreased risk of bronchiectasis, pneumonia, and influenza and pneumonia.

Introduction

Lower respiratory tract infections (LRTIs), which include acute bronchitis, acute bronchiolitis, bronchiectasis, pneumonia and others, remain to be the leading cause of morbidity and mortality around the world, particularly affecting adults over 70 years and children under 5 years old [1, 2]. LRTIs are caused by various pathogens, such as bacteria (Streptococcus pneumoniae, Haemophilus influenzae, Staphylococcus aureus, etc.), viruses (respiratory syncytial virus [RSV], influenza virus, human coronavirus [hCoV], etc.), mycoplasma, chlamydia, and legionella [2]. The systematic analysis for the Global Burden of Disease Study 2016 estimated that there were 2,377,697 deaths due to LRTIs globally in 2016 [3]. Over the past decades, many countries have taken a range of actions to prevent and treat LRTIs, including the improvement of nutrition and hygiene, distribution and usage of influenza antibiotics and antivirals, and global promotion of Haemophilus influenzae type b and pneumococcal conjugate vaccines [4, 5]. Despite remarkable progress has been made with these measures, LRTIs still bring huge burden to individuals and society. As a result, sustained efforts and preventive strategies remain significant for managing LRTIs.

For thousands of years, plants have exerted an important role in maintaining health and improving the quality of life for humans [6]. Tea is the most popular beverages after water, consumed by more than two-thirds of the general population worldwide [7]. Two polyphenols, catechins and theaflavins, were found in abundance in green and black tea, and were reported to have multiple physiological activities, such as antibacterial, antiviral, antioxidative, and anticancer effects, etc. [8, 9]. Accumulating evidence has linked tea intake to various health benefits, including the prevention of tumors, cardiovascular diseases, and kidney stones [10, 11]. The association between tea consumption and LRTIs risk has also been investigated in many epidemiologic researches. However, these studies failed to reach a consensus. Several studies observed that tea consumption has preventive effects against LRTIs [12, 13], but others did not find similar link [14, 15]. Traditional epidemiologic studies are susceptible to underlying confounding factors and reverse causality, which can overestimate or underestimate the causal association between determinants and outcomes. Thus, it is uncertain whether the observed relationship between tea intake and LRTIs risk is causal.

Mendelian randomization (MR), which integrates summary data from genome-wide association studies (GWAS) to infer the causality in the putative exposure–outcome pathway, is an emerging genetic epidemiologic method [16, 17]. Since the genetic variations are randomly arranged in meiosis and fixed after fertilization, MR approach enables to overcome the shortcomings of traditional epidemiologic studies, including confounding factors, reverse causation and selection biases [18]. The method has never been used to explore the association between tea intake and LRTIs risk. Hence, we conducted a two-sample MR analysis to determinate whether tea consumption is causally correlated with the risk of the common subtypes of LRTIs.

Materials and methods

All analyses performed in the current study were based on publicly available GWAS summary statistics. No additional ethical statements or informed consent were required.

Data source

Single-nucleotide polymorphisms (SNPs) were selected as instrumental variables (IVs) for exposure (tea intake) and outcome (LRTIs) from the GWAS databases. SNPs associated with tea intake were derived from a large GWAS conducted by Neale Lab, which involves 349,376 samples of European ancestry. The GWAS adjusted for sex, age, age2, sex × age, sex × age2, and first twenty ancestry principal components. Information about tea intake were acquired by a dietary questionnaire. Participants were asked the following question at baseline: “How many cups of tea do you drink every day? (Include black and green tea)”. Details regarding the study design, data analyses, and ethical approval were described at https://github.com/Nealelab/UK_Biobank_GWAS.

SNPs associated with common LRTI subtypes were extracted from the data of FinnGen Release 5 (released to public on May 11, 2021), including acute bronchitis, acute bronchiolitis, bronchiectasis, pneumonia, and influenza and pneumonia. We limited LRTI data to samples of European ancestry to avoid potential bias aroused by population stratification. The dataset consists of 7338 acute bronchitis cases (208,689 controls), 976 acute bronchiolitis cases (208,689 controls), 1107 bronchiectasis cases (186,723 controls), 27,370 pneumonia cases (191,422 controls), and 29,924 influenza and pneumonia cases (188,868 controls). FinnGen project is a public–private partnership that integrates disease endpoint genetic data from the Finnish biobanks and the Finnish health registries. Detailed documentations were recorded at FinnGen website (https://finngen.gitbook.io/documentation/).

Instrumental SNPs selection

To ensure the correctness and robustness of the conclusion on the causal link between tea intake and the risk of LRTIs, a series of quality control steps were adopted to select valid IVs associated with tea consumption. First, SNPs had a genome-wide significant association (P < 5E-08) with tea intake were selected as IVs. Second, the minor allele frequency (MAF) threshold of selected SNPs should more than 0.01. Third, due to the presence of strong linkage disequilibrium (LD) among the selected SNPs may bias the results, the clumping process (r2 < 0.001, clumping distance = 10,000 kb) were carried out to eliminate the LD between the included IVs. Fourth, to guarantee that the impact of SNPs on exposure corresponds to the same allele as the impact on the outcome, palindromic SNPs with intermediate allele frequencies were excluded. Fifth, when the exposure-related SNPs were not available in the outcome dataset, the proxy SNPs significantly correlated with the variants of interest were used (r2 > 0.8).

MR assumptions

To minimize the bias on the results, the MR approach must obey three key assumptions. First, the selected IVs must be closely associated with the exposure. Herein, the strength of the correlation between IVs and exposure was evaluated by F statistic, which is presented as R2(n−k−1)/k(1−R2). In this formula, R2 represents the cumulative interpreted variance of selected SNPs on exposure, k refers to the number of selected SNPs and n is the sample size. If F value is more than 10, the correlation is considered strong enough to avoid bias caused by weak IVs. Second, IVs are independent of confounders that influence exposure and outcome. Third, IVs affect the outcome merely via exposure, that is, there is no horizontal pleiotropy effect between IVs and outcome. Figure 1 presents an overview of the design for the current study.

Fig. 1

An overview of this Mendelian randomization study design

Effect size estimate

A two-sample MR analysis was conducted to evaluate the causal relationship between exposure (tea intake) and outcome (LRTIs). The inverse variance weighted (IVW) method was perceived as the primary analysis. IVW is essentially a meta-analysis approach, which translates into a weighted regression of outcome effects of IVs on exposure effects to acquire an overall estimate of causal effect, with a cut-off limit of zero. In the absence of horizontal pleiotropic effect, IVW can avoid the influence of confounding factors and acquire unbiased estimation [19]. In addition, three complementary methods were adopted to correct horizontal pleiotropy, despite at the expense of reduced statistical power. First, the MR-Egger regression was conducted. MR-Egger may be strongly affected by outlier genetic instruments, resulting in inaccurate estimates. However, even though all the selected IVs are invalid, MR-Egger method can still calculate unbiased estimates [20]. Second, the weighted median method was conducted, which combines data from multiple genetic instruments for consistency analysis by calculating a single weighted median estimator [21]. Third, the weighted mode method was performed, which is flexible for genetic variables that violate the pleiotropy hypothesis [22].

Sensitivity analysis

To ensure the reliability of the conclusion, six sensitivity analyses were adopted to verify whether heterogeneity and pleiotropy within the genetic variables could bias the MR results. First, the MR-Egger regression was applied to detect and adjust for the underlying horizontal pleiotropic effects among the selected IVs through the assessment of the intercept. Second, MR pleiotropy residual sum and outlier (MR-PRESSO) test was used to detect potential horizontal pleiotropy, and MR-PRESSO global test of heterogeneity was conducted to identify the underlying horizontal pleiotropy. Third, the asymmetry of funnel plots was also used as an indicator of horizontal pleiotropy. Fourth, the Cochran’s Q statistic was applied to quantify the heterogeneity across the selected genetic instruments. Fifth, the leave-one-out sensitivity analysis was implemented by deleting each SNP in turn, which can ensure that the MR estimates are not driven by certain strong SNPs. Finally, the PhenoScanner database was searched to check whether any of the SNPs were significantly associated (P < 5E-08) with other phenotypes at risk of affecting the five LRTIs subtypes independent of the tea consumption [23]. We further assessed the effect of removing these SNPs from the MR estimates to exclude potential pleiotropic effects.

A Bonferroni corrected P-value less than 0.01 (correcting for five outcomes) was defined as statistically significant. All data analyses were carried out using the “TwoSampleMR” and “MRPRESSO” packages of the R 4.0.3 software (www.r-project.org).

Results

Instrumental variables selection

A total of 45 SNPs (all P < 5E-08, r2 < 0.001) associated with tea consumption were identified, explaining 0.76% of the variance for tea intake. Supplementary Table S1 provided the detailed information of the 45 SNPs, including effect allele, other allele, beta, SE, P-value, F statistic, etc. The F statistic of these SNPs ranged from 30 to 408 with the median value of 40, which was greater than the conventional threshold of 10, indicating that the instrument bias was weak and cannot significantly affect the estimation of causal effects. Three SNPs (rs12901092, rs140775622, rs28676340) were excluded because they were not available in the outcome datasets, and four SNPs (rs11487328, rs2315024, rs397074, rs79413667) were excluded due to palindrome. Ultimately, a total of 38 SNPs with two proxy SNPs were selected as IVs for the following MR analysis.

MR analysis and sensitivity analysis

The results of IVW analyses revealed that genetically predicted an extra daily cup of tea intake was causally associated with the decreased risk of bronchiectasis (odds ratio (OR) = 0.61, 95% confidence interval (CI) = 0.47–0.78, P < 0.001), pneumonia (OR = 0.90, 95% CI = 0.85–0.96, P = 0.002), influenza and pneumonia (OR = 0.91, 95% CI = 0.85–0.97, P = 0.002) (Table 1, Fig. 2). Moreover, MR Egger, weighted median and weighted mode also obtained similar effect estimates, although wider CIs led to lower statistical power (Figs. 2, 3). The horizontal pleiotropy between IVs and outcomes was evaluated by MR-Egger regression, and the results indicated that there was weak evidence of horizontal pleiotropy (Table 1). Similarly, the MR-PRESSO test also identified no outlier SNPs or a horizontal pleiotropic effect of tea intake on the risk of bronchiectasis (P = 0.095), pneumonia (P = 0.121), and influenza and pneumonia (P = 0.117). The funnel plots showed symmetric results, by which the directional and horizontal pleiotropy was not significant (Supplementary Figure S1). The Cochran’s Q test suggested that there was no significant heterogeneity across the selected genetic instruments (P > 0.05; Table 1). The leave-one-out and forest plots showed that there was no single SNP driving the MR estimates, indicating that the results of the current MR analysis were robust (Supplementary Figure S2 and S3). Using the PhenoScanner tool, we found that two SNPs (rs199621380 and rs79217743) were significantly associated hematological traits (eg., monocyte count, eosinophil count and neutrophil percentage of white cells). It is known that changes in hematological parameters in patients with LRTIs are very common [24]. To avoid horizontal pleiotropy due to possible causal mechanistic associations between hematological traits and LRTIs subtypes, we removed these SNPs. Notably, the significant associations between tea intake and bronchiectasis, pneumonia and influenza and pneumonia did not change (Supplementary Figure S4).

Table 1

MR estimates from each method of assessing the causal effects of tea intake on LRTIs risk

Outcome	Nsnp	Methods	Beta	SE	OR (95% CI)	P-value	Horizontal pleiotropy			Heterogeneity
							Egger intercept	SE	P-value	Cochran’s Q	P-value
Bronchitis	38	MR Egger	− 0.07	0.13	0.93 (0.73–1.19)	0.572	− 0.0016	0.01	0.852	32.57	0.633
		IVW	− 0.09	0.05	0.91 (0.82–1.01)	0.067
		Weighted median	− 0.06	0.08	0.94 (0.81–1.10)	0.432
		Weighted mode	− 0.08	0.11	0.93 (0.75–1.14)	0.486
Bronchiolitis	38	MR Egger	− 0.24	0.348	0.79 (0.40–1.56)	0.504	0.0002	0.02	0.994	22.16	0.966
		IVW	− 0.23	0.141	0.79 (0.60–1.05)	0.100
		Weighted median	− 0.29	0.203	0.75 (0.50–1.11)	0.149
		Weighted mode	− 0.30	0.256	0.74 (0.45–1.23)	0.251
Bronchiectasis	38	MR Egger	− 0.86	0.32	0.42 (0.23–0.79)	0.011	0.0259	0.02	0.229	35.87	0.475
		IVW	− 0.50	0.13	0.61 (0.47–0.78)	< 0.001
		Weighted median	− 0.56	0.19	0.57 (0.39–0.84)	0.004
		Weighted mode	− 0.51	0.23	0.60 (0.38–0.95)	0.035
Pneumonia	38	MR Egger	− 0.10	0.08	0.90 (0.77–1.06)	0.233	− 0.0001	0.01	0.984	48.12	0.085
		IVW	− 0.10	0.03	0.90 (0.85–0.96)	0.002
		Weighted median	− 0.14	0.04	0.87 (0.79–0.94)	0.001
		Weighted mode	− 0.11	0.05	0.89 (0.80–1.00)	0.050
Influenza and pneumonia	38	MR Egger	− 0.09	0.08	0.92 (0.78–1.07)	0.289	− 0.0009	0.01	0.868	48.30	0.083
		IVW	− 0.10	0.03	0.91 (0.85–0.97)	0.002
		Weighted median	− 0.13	0.04	0.88 (0.81–0.95)	0.002
		Weighted mode	− 0.11	0.06	0.90 (0.80–1.00)	0.057

P < 0.01 was considered statistically significant

CI confidence interval, IVW inverse variance weighted, LRTI lower respiratory tract infection, MR Mendelian randomization, OR odds ratio, SNP single nucleotide polymorphism, SE standard error

Fig. 2

Mendelian randomization estimates of tea consumption with lower respiratory tract infections risk

Fig. 3

Scatterplot of the effect size for each SNP on tea consumption and the risk of bronchiectasis, pneumonia, and influenza and pneumonia

On the contrary, the results of IVW method revealed that genetically predicted tea intake was not causally associated with the risk of acute bronchitis (OR = 0.91, 95% CI = 0.82–1.01, P = 0.067) and acute bronchiolitis (OR = 0.79, 95% CI = 0.60–1.05, P = 0.100) (Table 1, Fig. 2). The methods of MR Egger, weighted median and weighted mode also provided consistent results (Table 1, Figs. 2, 4). The MR-Egger regression analyses indicated that there was weak evidence of horizontal pleiotropy (Table 1). Additionally, no outlier SNPs or a horizontal pleiotropic effect of tea intake on the risk of acute bronchitis (P = 0.629) and acute bronchiolitis (P = 0.986) were identified by the MR-PRESSO approaches. The funnel plots were also visually symmetric (Supplementary Figure S5). The Cochran’s Q test found no significant heterogeneity across the selected IVs (P > 0.05; Table 1). The leave-one-out and forest plots further proved the robustness of the MR results (Supplementary Figure S6 and S7).

Fig. 4

Scatterplot of the effect size for each SNP on tea consumption and the risk of acute bronchitis and acute bronchiolitis

Discussion

To the best of our knowledge, this is the first study to determinate the causal relationship between tea consumption and LRTIs risk. The main finding was that genetically predicted higher tea intake was causally associated with the decreased risk of bronchiectasis, pneumonia, and influenza and pneumonia. However, we did not find the association between tea intake and the risk of acute bronchitis and acute bronchiolitis.

The association between tea consumption and LRTIs risk remains controversial in previous epidemiologic studies. In a large cohort of 19,079 males and 21,493 males, Watanabe et al. found that tea consumption was associated with a lower risk of death from pneumonia in Japanese women [13]. In the CoPanFlu-France cohort of 1121 participants, tea consumption at least twice a week was a protective factor for influenza infection [25]. In an observational study, Park et al. reported that the consumption of 1–5 cups per day of tea may prevent influenza infection in children [12]. However, a case–control study conducted by Kondo et al. observed no association between pneumonia and tea consumption [15]. The discrepancy between these studies may be caused by the residual confounders, such as selection bias, recall bias, and detection bias.

In the current MR analysis, genetically predicted intake of one extra cup of tea per day was associated with 39% lower risk of bronchiectasis, 10% lower risk of pneumonia, and 9% lower risk of influenza and pneumonia. Several potential biological mechanisms may account for the inverse association between tea intake and these LRTIs. LRTIs are usually caused by various pathogenic microorganisms. Experimental studies have indicated that tea may be effective against many Gram-positive and gram-negative microorganisms as well as some viruses, fungi and parasites [26-28]. Polyphenols are the most important components in tea, which exert antimicrobial effects by damaging the bacterial cell membrane, inhibiting fatty acid synthesis, enzyme activity and inflammation [26]. Moreover, polyphenols can also bind to virus’s hemagglutinin molecule, thereby inhibiting the virus’s adsorption to host cells and preventing its assembly or maturation lysis [29]. Epigallocatechin gallate (EGCG), a major and highly bioactive catechin, was verified to reduce the infectivity of influenza A virus and influenza B virus in Madin–Darby canine kidney cells [27, 30]. ECGG can not only inhibit the migration of neutrophils through endothelial monolayer, thus decreasing vascular permeability [31], but also reduce neutrophil elastase, which is a proteolytic molecule involved in the body’s inflammatory response that has been associated with increased alveolar epithelial permeability [32]. The above provides a theoretical basis for the role of tea in preventing bronchiectasis, pneumonia and influenza, but the specific mechanisms remain to be further elucidated.

The strengths of this study are that it has large sample size, and it is the first MR analysis to evaluate the causal association between tea consumption and LRTIs risk. The application of MR approach reduces the interference of confounding factors and reverse causality on the results and is more convincing than observational studies. However, several shortcomings worth further consideration. First, the data of this study were derived from large-scale GWAS analyses, and the detailed information of the subjects and tea type (black or green tea) was not available, which makes it impossible to conduct stratification analyses or analyses adjusted for other covariables. Second, all the included subjects were of European ancestry, so the extrapolation of the findings to other ethnic groups may be limited. Third, there might be overlapping participants in the exposure and outcome studies, but the degree of sample overlap was difficult to estimate. Fortunately, the strong IVs (F statistic much larger than 10) used in this study minimizes the bias due to sample overlap. Finally, standard MR method assumes a linear relationship between exposure and outcome, so the non-linear association and threshold effect between tea intake and LRTIs cannot be detected.

In summary, our study found that genetically predicted an extra daily cup of tea intake may causally associated with the reduced risk of bronchiectasis, pneumonia and influenza and pneumonia. The finding of this study provided new evidence for the prevention and treatment of these LRTIs.

Below is the link to the electronic supplementary material.

Supplementary file1 (DOCX 1329 KB)