ORIGINAL RESEARCH article

Front. Oncol., 18 June 2024

Sec. Thoracic Oncology

Volume 14 - 2024 | https://doi.org/10.3389/fonc.2024.1373380

Causal association of circulating cytokines with the risk of lung cancer: a Mendelian randomization study

  • 1. Department of Respiratory and Critical Care Medicine, Affiliated Hospital of North Sichuan Medical College, Nanchong, China

  • 2. Department of Pulmonary and Critical Care Medicine, Center of Respiratory Medicine, National Center for Respiratory Medicine, China-Japan Friendship Hospital, Beijing, China

Abstract

Background:

Lung cancer is the deadliest and most prevalent malignancy worldwide. While smoking is an established cause, evidence to identify other causal factors remains lacking. Current research indicates chronic inflammation is involved in tumorigenesis and cancer development, though the specific mechanisms underlying the role of inflammatory cytokines in lung cancer pathogenesis remain unclear. This study implemented Mendelian randomization (MR) analysis to investigate the causal effects of circulating cytokines on lung cancer development.

Methods:

We performed a two-sample MR analysis in Europeans utilizing publicly available genome-wide association study summary statistics. Single nucleotide polymorphisms significantly associated with cytokine were selected as genetic instrumental variables.

Results:

Genetically predicted levels of the chemokine interleukin-18 (IL-18) (OR = 0.942, 95% CI: 0.897–0.990, P = 0.018) exerted significant negative causal effects on overall lung cancer risk in this analysis. Examining specific histologic subtypes revealed further evidence of genetic associations. Stem cell factor (SCF) (OR = 1.150, 95% CI: 1.021–1.296, P = 0.021) and interleukin-1beta (IL-1β) (OR = 1.152, 95% CI: 1.003–1.325, P = 0.046) were positively associated with lung adenocarcinoma risk, though no inflammatory factors showed causal links to squamous cell lung cancer risk. Stratified by smoking status, interferon gamma-induced protein 10 (IP-10) (OR = 0.861, 95% CI: 0.781–0.950, P = 0.003) was inversely associated while IL-1β (OR = 1.190, 95% CI: 1.023–1.384, P = 0.024) was positively associated with lung cancer risk in ever smokers. Among never smokers, a positive association was observed between lung cancer risk and SCF (OR = 1.474, 95% CI: 1.105–1.964, P = 0.008). Importantly, these causal inferences remained robust across multiple complementary MR approaches, including MR-Egger, weighted median, weighted mode and simple mode regressions. Sensitivity analyses also excluded potential bias stemming from pleiotropy.

Conclusion:

This MR study found preliminary evidence that genetically predicted levels of four inflammatory cytokines—SCF, IL-1β, IL-18, and IP-10—may causally influence lung cancer risk in an overall and subtype-specific manner, as well as stratified by smoking status. Identifying these cytokine pathways that may promote lung carcinogenesis represents potential new targets for the prevention, early detection, and treatment of this deadly malignancy.

Introduction

As the deadliest cancer worldwide, lung cancer was responsible for 1,817,000 deaths in 2022, accounting for 18.7% of all cancer-related fatalities (). Due to typically late-stage diagnosis once symptoms arise, the 5-year survival rate for lung cancer patients is dismally low at around 20% (). Smoking is a well-known predominant risk factor, yet evidence to elucidate other potential causes remains scarce (, ). Identifying additional contributory risk factors to promote early detection and treatment thus represents a critical need in combating this devastating disease.

Over the past two decades, a sizable body of evidence has firmly established chronic infection and inflammation as key promoters of carcinogenesis (–). The inflammatory tumor microenvironment, comprised of leukocytes releasing cytokines, chemokines, reactive oxygen species and other cytotoxic mediators, can drive tumor progression through processes such as genotoxicity, aberrant tissue repair, heightened cellular proliferation, invasion and metastasis (, ). Critical transcription factors including STAT3 and NF-κB have been implicated in inflammation-fueled carcinogenesis (, ). Tumors can further manipulate the inflammatory milieu to suppress anti-tumor T cell responses (). However, despite substantial links between cytokines and cancer, the precise underlying mechanisms, particularly in lung cancer, remain to be fully deciphered.

Mendelian randomization (MR) analysis utilizes genetic variation as an instrumental variable to infer the potential causal relationship between modifiable exposures and disease outcomes (). As genetic variation is randomly assigned and not susceptible to reverse causation, this ingenious approach circumvents confounding and reverse causation bias that plague traditional observational studies (). In this study, we implemented MR to investigate the causal role of circulating cytokines in lung cancer pathogenesis. Using genetic instruments as surrogates for cytokine levels, we evaluated whether inflammation actively drives lung cancer development or simply represents an epiphenomenon. Elucidating these causal pathways will provide novel mechanistic insights and could reveal previously unrecognized therapeutic targets for this pernicious lung disease.

Method

Study design

In this two-sample MR analysis, we utilized single nucleotide polymorphisms (SNPs) as instrumental variables (IVs). To ensure validity, SNPs were chosen based on three critical assumptions: (1) significant associations with exposure factors, satisfying relevance; (2) effects on outcomes solely via exposures, not through alternative pathways, fulfilling exclusivity; and (3) independence from confounding factors ().

Data resource

SNPs associated with circulating cytokines were obtained from the latest genome-wide association studies (GWAS), as listed in Supplementary Table 1. Summary data from a large-scale cytokine GWAS meta-analysis were used to generate genetic instruments for cytokines (). This meta-analysis measured cytokine levels in plasma and blood samples from 8,293 Finnish individuals across three population-based cohorts.

To investigate causal effects of circulating cytokines on lung cancer risk, we obtained lung cancer GWAS summary statistics from the IEU open database. Lung cancer data were obtained from the International Lung Cancer Consortium (ILCCO) (https://ilcco.iarc.fr/), including 29,836 cases and 55,586 controls (). The study also provided associations between instrumental SNPs and different histologic subtypes of lung cancer, including lung adenocarcinoma (case: 11,245, control: 54,619), squamous cell carcinoma (case: 7,704, control: 54,763), and small cell lung cancer (case: 2,791, control: 20,580). Subgroup analyses were performed according to smoking status, including smokers (23,223 cases and 16,964 controls) and never-smokers (2,355 cases and 7,504 controls), and all populations were restricted to European ethnic groups.

SNPs selection

To identify valid SNPs for MR, we implemented several filtering steps: First, we selected independent SNPs strongly associated with different cytokine levels (P <5×10−6) to maximize instrument availability (, ). Second, we used clumping to prune correlated SNPs in linkage disequilibrium (r2<0.001, 10,000 kb) and avoid biased results. Third, we excluded pleiotropic SNPs associated with potential confounders including smoking, diabetes and anxiety using PhenoScanner (). Fourth, we retained only concordant SNPs between exposure and outcome datasets as valid instruments. Finally, we excluded weak instruments with F-statistics <10 (calculated as F = R2×(N-2)/1-R2) to minimize bias ().

Statistical analysis

After SNP filtering, our primary MR analysis utilized inverse-variance weighted (IVW) estimation to evaluate the overall causal effects, given its accuracy when all instruments are valid (). Complementary approaches including weighted median, MR Egger, weighted mode and simple mode were also implemented (, ). To probe potential horizontal pleiotropy, we performed MR Egger regression and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) testing (, ). Heterogeneity was assessed via Cochran’s Q and MR Egger regression, and result robustness was verified through leave-one-out analysis. Furthermore, we conducted Steiger testing to evaluate possible reverse causation ().

To correct for multiple hypothesis testing, we applied Bonferroni correction and set statistical significance at P < 0.0012 based on the number of cytokines analyzed. P-values ranging from 0.0012 to 0.05 were considered suggestive evidence for potential causality (, ). All Mendelian randomization analyses were implemented in R utilizing the TwoSampleMR and MR-PRESSO packages. A P-value below 0.05 was deemed statistically significant.

Results

Causal effects of circulating cytokines on lung cancer and different subtypes

After implementing the described filtering procedures, 4–16 SNPs were retained as instruments for the circulating cytokines (Supplementary Table 1). High F-statistics confirmed all selected SNPs were robust instruments (all F values >10). We then leveraged these SNPs to infer causality between cytokines and lung cancer. MR estimates from the various analytical methods are visualized in Figure 1 (Supplementary Table S2). Specifically, IVW analysis revealed interleukin-18 (IL-18) exerted significant negative causal effects on lung cancer risk (OR = 0.942, 95% CI: 0.897–0.990, P = 0.018), while eotaxin had significant positive causal effects (OR = 1.061, 95% CI: 1.002–1.123, P = 0.043) (Figures 1A, B). Eotaxin, however, has a statistical efficacy of only 67%, thus it remains to be verified. Consistent estimates were yielded by MR Egger and weighted median methods. Scatter plots visualized the effects for each method across datasets (Supplementary Figure 1).

Figure 1

Among the different histologic lung cancer types, we found further evidence of genetic associations (Figures 2-4). No inflammatory factor was found to be causally associated with the risk of squamous cell lung cancer (Figure 2) (Supplementary Table S3). Stem cell factor (SCF) (OR = 1.150, 95% CI: 1.021–1.296, P = 0.021) and interleukin-1beta (IL-1β) (OR = 1.152, 95% CI: 1.003–1.325, P = 0.046) were positively associated with the risk of lung adenocarcinoma (Figures 3A, B) (Supplementary Table S4). Vascular endothelial growth factor (VEGF) (OR = 1.117, 95% CI: 1.008–1.237, P = 0.035) was positively associated with the risk of small cell lung carcinoma (Figure 4A, B) (Supplementary Table S5), but the statistical efficacy of VEGF was limited to 21%, which makes the causal hypothesis unreliable. Scatter plots visualized the effects for each method across datasets (Supplementary Figures 2, 3).

Figure 2

Figure 3

Figure 4

MR Egger regression and MR-PRESSO global testing showed no evidence of horizontal pleiotropy (Table 1 and Supplementary Tables 6-9). Importantly, Cochran’s Q statistics and MR Egger regression revealed no significant heterogeneity among the individual SNP instruments (P > 0.05). MR-PRESSO global test p-values exceeding 0.05 for all cytokine-lung cancer associations, including different pathological types, further ruled out pleiotropy. Leave-one-out analysis confirmed the robustness of the causal estimates (Supplementary Figures 4-6).

Table 1

MethodBetaSEOR95% CIP -valueLOWUP
CTACK levelsInverse variance weighted-0.0590.0480.9420.858-1.0350.2140.8581.035
beta-nerve growth factor levelsInverse variance weighted0.0760.0761.0790.930-1.2530.3170.9301.253
Vascular endothelial growth factor levelsInverse variance weighted-0.0110.0330.9890.927-1.0560.7490.9271.056
Macrophage Migration Inhibitory Factor levelsInverse variance weighted-0.0700.0630.9320.824-1.0550.2650.8241.055
TRAIL levelsInverse variance weighted0.0180.0331.0180.955-1.0860.5800.9551.086
Tumor necrosis factor beta levelsInverse variance weighted-0.0450.0400.9560.885-1.0340.2610.8851.034
Tumor necrosis factor alpha levelsInverse variance weighted-0.0040.0580.9960.889-1.1160.9410.8891.116
Stromal-cell-derived factor 1 alpha levelsInverse variance weighted0.0700.0891.0720.901-1.2770.4330.9011.277
Stem cell growth factor beta levelsInverse variance weighted0.0080.0381.0080.936-1.0860.8350.9361.086
Stem cell factor levelsInverse variance weighted0.1090.0691.1150.975-1.2760.1120.9751.276
Interleukin-16 levelsInverse variance weighted0.0100.0351.0100.943-1.0820.7720.9431.082
RANTES levelsInverse variance weighted0.0100.0611.0100.897-1.1380.8630.8971.138
Platelet-derived growth factor BB levelsInverse variance weighted-0.0110.0500.9890.897-1.0920.8330.8971.092
Macrophage inflammatory protein 1b levelsInverse variance weighted-0.0320.0310.9690.911-1.0300.3070.9111.030
Macrophage inflammatory protein 1a levelsInverse variance weighted-0.0240.0560.9760.875-1.0890.6660.8751.089
Monokine induced by gamma interferon levelsInverse variance weighted-0.0220.0420.9780.902-1.0610.5940.9021.061
Macrophage colony stimulating factor levelsInverse variance weighted-0.0090.0520.9920.895-1.0990.8710.8951.099
Monocyte chemoattractant protein-3 levelsInverse variance weighted0.0710.0601.0740.955-1.2070.2360.9551.207
Monocyte chemoattractant protein-1 levelsInverse variance weighted0.1070.0541.1131.000-1.2380.0501.0001.238
Interleukin-12p70 levelsInverse variance weighted0.0240.0411.0250.945-1.1110.5550.9451.111
Interferon gamma-induced protein 10 levelsInverse variance weighted0.0260.0551.0260.922-1.1430.6360.9221.143
Interleukin-18 levelsInverse variance weighted-0.0070.0380.9930.922-1.0700.8570.9221.070
Interleukin-17 levelsInverse variance weighted-0.0850.0700.9190.802-1.0530.2240.8021.053
Interleukin-13 levelsInverse variance weighted0.0260.0341.0270.960-1.0970.4390.9601.097
Interleukin-10 levelsInverse variance weighted0.0040.0611.0040.891-1.1310.9470.8911.131
Interleukin-8 levelsInverse variance weighted0.0270.0851.0280.869-1.2150.7510.8691.215
Interleukin-6 levelsInverse variance weighted0.0180.0981.0180.840-1.2330.8570.8401.233
Interleukin-1-receptor antagonist levelsInverse variance weighted0.0230.0611.0230.907-1.1540.7130.9071.154
Interleukin-1-beta levelsInverse variance weighted-0.0170.0840.9830.834-1.1600.8420.8341.160
Hepatocyte growth factor levelsInverse variance weighted-0.0650.0710.9370.815-1.0770.3590.8151.077
Interleukin-9 levelsInverse variance weighted-0.0410.0650.9600.845-1.0910.5310.8451.091
Interleukin-7 levelsInverse variance weighted0.0380.0381.0390.964-1.1190.3190.9641.119
Interleukin-5 levelsInverse variance weighted0.0250.0671.0250.900-1.1690.7080.9001.169
Interleukin-4 levelsInverse variance weighted0.0400.0801.0410.890-1.2180.6160.8901.218
Interleukin-2 receptor antagonist levelsInverse variance weighted0.0370.0391.0380.961-1.1200.3450.9611.120
Interleukin-2 levelsInverse variance weighted-0.0140.0560.9860.883-1.1010.8050.8831.101
Interferon gamma levelsInverse variance weighted-0.1400.0830.8690.739-1.0220.0900.7391.022
Growth-regulated protein alpha levelsInverse variance weighted0.0260.0391.0270.952-1.1080.4980.9521.108
Granulocyte-colony stimulating factor levelsInverse variance weighted-0.0040.0980.9960.821-1.2070.9640.8211.207
Fibroblast growth factor basic levelsInverse variance weighted0.1260.1201.1340.896-1.4360.2940.8961.436
Eotaxin levelsInverse variance weighted0.0460.0461.0470.957-1.1450.3180.9571.145

Primary results of MR analysis on squamous cell lung cancer.

Among the other 35 examined cytokines, none demonstrated significant correlation with lung cancer risk or pathological types in IVW or secondary MR analyses (Tables 1–4). Heterogeneity testing showed no significant heterogeneity for any of the cytokines (Tables 5–8). Across all cytokines, MR Egger regression consistently revealed no evidence of pleiotropy (Tables 5–8). MR-PRESSO outlier testing validated the significant MR findings, except for monocyte chemoattractant protein-3 (MCP-3) in lung cancer (including different pathologic types) (Supplementary Tables 6-9), and macrophage migration inhibitory factor (MIF) and IL-5 in small-cell lung carcinoma, where limited SNPs were available (Supplementary Table 9). Furthermore, Steiger p-values <0.05 verified the detected causal direction was correct for all cytokines (Table 9–12).

Table 2

MethodBetaSEOR95% CIP -valueLOWUP
CTACK levelsInverse variance weighted-0.0560.0330.9460.887-1.0090.0900.8871.009
beta-nerve growth factor levelsInverse variance weighted0.0070.0391.0070.933-1.0860.8570.9331.086
Vascular endothelial growth factor levelsInverse variance weighted0.0190.0211.0190.978-1.0620.3690.9781.062
Macrophage Migration Inhibitory Factor levelsInverse variance weighted-0.0270.0410.9740.899-1.0550.5140.8991.055
TRAIL levelsInverse variance weighted0.0380.0231.0390.992-1.0880.1010.9921.088
Tumor necrosis factor beta levelsInverse variance weighted0.0170.0251.0170.969-1.0680.4840.9691.068
Tumor necrosis factor alpha levelsInverse variance weighted-0.0220.0360.9780.911-1.0500.5380.9111.050
Stromal-cell-derived factor 1 alpha levelsInverse variance weighted0.0500.0571.0510.941-1.1750.3790.9411.175
Stem cell growth factor beta levelsInverse variance weighted0.0330.0291.0330.976-1.0940.2600.9761.094
Stem cell factor levelsInverse variance weighted0.0860.0481.0900.992-1.1980.0730.9921.198
Interleukin-16 levelsInverse variance weighted0.0230.0221.0230.980-1.0680.2990.9801.068
RANTES levelsInverse variance weighted-0.0330.0330.9670.907-1.0320.3170.9071.032
Platelet-derived growth factor BB levelsInverse variance weighted-0.0620.0350.9400.878-1.0060.0720.8781.006
Macrophage inflammatory protein 1b levelsInverse variance weighted-0.0060.0200.9940.956-1.0340.7740.9561.034
Macrophage inflammatory protein 1a levelsInverse variance weighted-0.0280.0350.9720.907-1.0420.4220.9071.042
Monokine induced by gamma interferon levelsInverse variance weighted-0.0060.0270.9940.943-1.0480.8150.9431.048
Macrophage colony stimulating factor levelsInverse variance weighted0.0120.0281.0120.959-1.0680.6650.9591.068
Monocyte chemoattractant protein-3 levelsInverse variance weighted0.0030.0381.0030.931-1.0790.9470.9311.079
Monocyte chemoattractant protein-1 levelsInverse variance weighted0.0480.0311.0490.988-1.1150.1180.9881.115
Interleukin-12p70 levelsInverse variance weighted0.0400.0261.0410.989-1.0960.1260.9891.096
Interferon gamma-induced protein 10 levelsInverse variance weighted-0.0240.0360.9770.909-1.0490.5170.9091.049
Interleukin-18 levels*Inverse variance weighted-0.0590.0250.9420.897-0.9900.0180.8970.990
Interleukin-17 levelsInverse variance weighted-0.0080.0440.9920.911-1.0810.8610.9111.081
Interleukin-13 levelsInverse variance weighted0.0320.0221.0320.989-1.0770.1490.9891.077
Interleukin-10 levelsInverse variance weighted0.0170.0371.0180.946-1.0940.6390.9461.094
Interleukin-8 levelsInverse variance weighted0.0040.0591.0040.894-1.1280.9470.8941.128
Interleukin-6 levelsInverse variance weighted0.0610.0611.0630.944-1.1970.3130.9441.197
Interleukin-1-receptor antagonist levelsInverse variance weighted0.0240.0391.0240.949-1.1040.5430.9491.104
Interleukin-1-beta levelsInverse variance weighted0.0950.0541.1000.990-1.2220.0750.9901.222
Hepatocyte growth factor levelsInverse variance weighted0.0090.0451.0090.924-1.1010.8500.9241.101
Interleukin-9 levelsInverse variance weighted-0.0110.0410.9890.913-1.0720.7890.9131.072
Interleukin-7 levelsInverse variance weighted0.0250.0241.0260.978-1.0750.2940.9781.075
Interleukin-5 levelsInverse variance weighted0.0460.0421.0470.964-1.1380.2750.9641.138
Interleukin-4 levelsInverse variance weighted-0.0020.0550.9980.897-1.1120.9760.8971.112
Interleukin-2 receptor antagonist levelsInverse variance weighted-0.0270.0250.9730.927-1.0210.2680.9271.021
Interleukin-2 levelsInverse variance weighted-0.0190.0330.9810.921-1.0470.5680.9211.047
Interferon gamma levelsInverse variance weighted0.0100.0561.0100.905-1.1270.8590.9051.127
Growth-regulated protein alpha levelsInverse variance weighted0.0340.0201.0350.994-1.0770.0910.9941.077
Granulocyte-colony stimulating factor levelsInverse variance weighted0.0350.0481.0360.943-1.1370.4650.9431.137
Fibroblast growth factor basic levelsInverse variance weighted0.0850.0761.0880.937-1.2640.2680.9371.264
Eotaxin levels*Inverse variance weighted0.0590.0291.0611.002-1.1230.0431.0021.123

Primary results of MR analysis on lung cancer.

Table 3

MethodBetaSEOR95% CIP -valueLOWUP
CTACK levelsInverse variance weighted-0.0100.0480.9900.902-1.0870.8380.9021.087
beta-nerve growth factor levelsInverse variance weighted-0.0290.0530.9720.875-1.0790.5930.8751.079
Vascular endothelial growth factor levelsInverse variance weighted0.0290.0291.0290.972-1.0900.3190.9721.090
Macrophage Migration Inhibitory Factor levelsInverse variance weighted-0.0180.0550.9820.881-1.0930.7370.8811.093
TRAIL levelsInverse variance weighted0.0310.0281.0320.976-1.0900.2700.9761.090
Tumor necrosis factor beta levelsInverse variance weighted0.0180.0341.0180.953-1.0870.5920.9531.087
Tumor necrosis factor alpha levelsInverse variance weighted0.0260.0501.0270.932-1.1320.5950.9321.132
Stromal-cell-derived factor 1 alpha levelsInverse variance weighted-0.0020.0890.9980.838-1.1880.9790.8381.188
Stem cell growth factor beta levelsInverse variance weighted-0.0220.0370.9780.910-1.0520.5540.9101.052
Stem cell factor levels*Inverse variance weighted0.1400.0611.1501.021-1.2960.0211.0211.296
Interleukin-16 levelsInverse variance weighted0.0000.0311.0000.942-1.0630.9920.9421.063
RANTES levelsInverse variance weighted-0.0200.0570.9800.877-1.0950.7190.8771.095
Platelet-derived growth factor BB levelsInverse variance weighted-0.0060.0490.9940.903-1.0950.9020.9031.095
Macrophage inflammatory protein 1b levelsInverse variance weighted0.0440.0291.0450.988-1.1060.1200.9881.106
Macrophage inflammatory protein 1a levelsInverse variance weighted-0.0370.0490.9640.876-1.0600.4470.8761.060
Monokine induced by gamma interferon levelsInverse variance weighted0.0330.0341.0330.966-1.1050.3370.9661.105
Macrophage colony stimulating factor levelsInverse variance weighted0.0170.0381.0170.945-1.0950.6480.9451.095
Monocyte chemoattractant protein-3 levelsInverse variance weighted-0.0070.0590.9930.885-1.1140.9050.8851.114
Monocyte chemoattractant protein-1 levelsInverse variance weighted0.0140.0421.0140.933-1.1020.7440.9331.102
Interleukin-12p70 levelsInverse variance weighted0.0470.0371.0480.975-1.1260.2010.9751.126
Interferon gamma-induced protein 10 levelsInverse variance weighted-0.0080.0450.9920.908-1.0830.8590.9081.083
Interleukin-18 levelsInverse variance weighted-0.0540.0310.9470.891-1.0070.0830.8911.007
Interleukin-17 levelsInverse variance weighted-0.0310.0650.9700.853-1.1020.6350.8531.102
Interleukin-13 levelsInverse variance weighted0.0280.0311.0280.967-1.0930.3710.9671.093
Interleukin-10 levelsInverse variance weighted0.0520.0441.0530.966-1.1470.2390.9661.147
Interleukin-8 levelsInverse variance weighted0.0120.0691.0120.885-1.1580.8580.8851.158
Interleukin-6 levelsInverse variance weighted0.0200.1001.0200.839-1.2410.8390.8391.241
Interleukin-1-receptor antagonist levelsInverse variance weighted0.0720.0531.0750.969-1.1930.1730.9691.193
Interleukin-1-beta levels*Inverse variance weighted0.1420.0711.1521.003-1.3250.0461.0031.325
Hepatocyte growth factor levelsInverse variance weighted0.0500.0671.0520.922-1.2000.4540.9221.200
Interleukin-9 levelsInverse variance weighted-0.0140.0570.9870.882-1.1030.8110.8821.103
Interleukin-7 levelsInverse variance weighted0.0050.0341.0050.941-1.0750.8740.9411.075
Interleukin-5 levelsInverse variance weighted0.0790.0591.0820.964-1.2140.1790.9641.214
Interleukin-4 levelsInverse variance weighted0.0300.0601.0310.916-1.1600.6130.9161.160
Interleukin-2 receptor antagonist levelsInverse variance weighted-0.0360.0390.9650.895-1.0410.3580.8951.041
Interleukin-2 levelsInverse variance weighted-0.0180.0430.9820.903-1.0680.6690.9031.068
Interferon gamma levelsInverse variance weighted0.0860.0771.0900.938-1.2670.2610.9381.267
Growth-regulated protein alpha levelsInverse variance weighted0.0080.0281.0080.954-1.0650.7740.9541.065
Granulocyte-colony stimulating factor levelsInverse variance weighted0.1150.0651.1220.988-1.2730.0750.9881.273
Fibroblast growth factor basic levelsInverse variance weighted-0.0390.1430.9620.727-1.2730.7850.7271.273
Eotaxin levelsInverse variance weighted0.0550.0461.0560.966-1.1550.2320.9661.155

Primary results of MR analysis on lung adenocarcinoma.

Table 4

MethodBetaSEOR95% CIP -valueLOWUP
CTACK levelsInverse variance weighted-0.0170.0800.9840.841-1.1500.8360.8411.150
beta-nerve growth factor levelsInverse variance weighted-0.0960.1220.9080.715-1.1540.4310.7151.154
Vascular endothelial growth factor levels*Inverse variance weighted0.1100.0521.1171.008-1.2370.0351.0081.237
Macrophage Migration Inhibitory Factor levelsInverse variance weighted-0.1300.1370.8780.671-1.1500.3450.6711.150
TRAIL levelsInverse variance weighted0.0120.0551.0120.908-1.1280.8250.9081.128
Tumor necrosis factor beta levelsInverse variance weighted0.0960.0621.1000.975-1.2420.1220.9751.242
Tumor necrosis factor alpha levelsInverse variance weighted0.0310.1091.0310.833-1.2760.7770.8331.276
Stromal-cell-derived factor 1 alpha levelsInverse variance weighted0.0730.1441.0750.811-1.4250.6130.8111.425
Stem cell growth factor beta levelsInverse variance weighted-0.1270.0710.8810.767-1.0120.0740.7671.012
Stem cell factor levelsInverse variance weighted0.0380.1171.0390.826-1.3060.7460.8261.306
Interleukin-16 levelsInverse variance weighted-0.0640.0640.9380.827-1.0640.3210.8271.064
RANTES levels*Inverse variance weighted-0.0830.1010.9210.755-1.1220.4120.7551.122
Platelet-derived growth factor BB levelsInverse variance weighted0.1300.0861.1390.963-1.3480.1300.9631.348
Macrophage inflammatory protein 1b levelsInverse variance weighted0.0590.0501.0600.961-1.1700.2430.9611.170
Macrophage inflammatory protein 1a levelsInverse variance weighted-0.0740.0920.9290.776-1.1110.4180.7761.111
Monokine induced by gamma interferon levelsInverse variance weighted-0.1260.0770.8810.758-1.0250.1010.7581.025
Macrophage colony stimulating factor levelsInverse variance weighted-0.0650.1070.9370.760-1.1560.5450.7601.156
Monocyte chemoattractant protein-3 levelsInverse variance weighted-0.0490.1050.9520.775-1.1690.6400.7751.169
Monocyte chemoattractant protein-1 levelsInverse variance weighted0.0520.0811.0540.900-1.2340.5150.9001.234
Interleukin-12p70 levelsInverse variance weighted-0.1720.1270.8420.657-1.0800.1750.6571.080
Interferon gamma-induced protein 10 levelsInverse variance weighted0.0260.1131.0260.823-1.2800.8190.8231.280
Interleukin-18 levelsInverse variance weighted0.0200.0951.0210.847-1.2300.8290.8471.230
Interleukin-17 levelsInverse variance weighted0.1340.1251.1440.895-1.4620.2840.8951.462
Interleukin-13 levelsInverse variance weighted0.0450.0611.0460.929-1.1780.4570.9291.178
Interleukin-10 levelsInverse variance weighted0.0880.0911.0920.914-1.3050.3310.9141.305
Interleukin-8 levelsInverse variance weighted-0.0260.1010.9740.798-1.1880.7960.7981.188
Interleukin-6 levelsInverse variance weighted-0.2620.1520.7700.571-1.0370.0850.5711.037
Interleukin-1-receptor antagonist levelsInverse variance weighted-0.0030.0950.9970.828-1.2000.9740.8281.200
Interleukin-1-beta levelsInverse variance weighted0.1460.1501.1570.863-1.5520.3290.8631.552
Hepatocyte growth factor levelsInverse variance weighted0.1000.1121.1050.886-1.3770.3750.8861.377
Interleukin-9 levelsInverse variance weighted-0.1500.1340.8600.662-1.1180.2610.6621.118
Interleukin-7 levelsInverse variance weighted0.0980.0811.1020.941-1.2920.2270.9411.292
Interleukin-5 levelsInverse variance weighted-0.2920.1710.7470.534-1.0440.0880.5341.044
Interleukin-4 levelsInverse variance weighted-0.1570.1650.8540.619-1.1790.3390.6191.179
Interleukin-2 receptor antagonist levelsInverse variance weighted0.0250.0621.0260.908-1.1590.6850.9081.159
Interleukin-2 levelsInverse variance weighted-0.0310.0810.9700.828-1.1360.7020.8281.136
Interferon gamma levelsInverse variance weighted0.0630.1291.0650.827-1.3720.6270.8271.372
Growth-regulated protein alpha levelsInverse variance weighted-0.0210.0500.9790.887-1.0800.6720.8871.080
Granulocyte-colony stimulating factor levelsInverse variance weighted-0.0670.1240.9360.733-1.1940.5920.7331.194
Fibroblast growth factor basic levelsInverse variance weighted0.2270.2211.2550.814-1.9360.3050.8141.936
Eotaxin levelsInverse variance weighted0.0930.0741.0980.950-1.2680.2050.9501.268

Primary results of MR analysis on small cell lung carcinoma.

Table 5

HeterogenityMR-Egger intercept
QQ_P -valueEgger_interceptSEP -value
CTACK levels9.4670.2210.0150.0170.437
beta-nerve growth factor levels1.9720.922-0.0100.0250.701
Vascular endothelial growth factor levels4.6330.865-0.0020.0080.837
Macrophage Migration Inhibitory Factor levels5.1690.396-0.0100.0220.680
TRAIL levels18.4760.1860.0170.0080.058
Tumor necrosis factor beta levels1.3680.7130.0010.0120.958
Tumor necrosis factor alpha levels1.5200.8230.0080.0130.601
Stromal-cell-derived factor 1 alpha levels3.2510.8610.0000.0120.969
Stem cell growth factor beta levels4.4100.818-0.0110.0130.455
Stem cell factor levels10.4680.2340.0150.0130.297
Interleukin-16 levels3.8880.867-0.0050.0120.702
RANTES levels5.0890.748-0.0230.0190.259
Platelet-derived growth factor BB levels4.7140.9090.0010.0090.881
Macrophage inflammatory protein 1b levels10.3450.8480.0100.0070.175
Macrophage inflammatory protein 1a levels7.0540.531-0.0160.0170.387
Monokine induced by gamma interferon levels6.7540.8190.0010.0130.940
Macrophage colony stimulating factor levels3.1570.870-0.0020.0160.898
Monocyte chemoattractant protein-3 levels0.6500.7220.0240.0380.643
Monocyte chemoattractant protein-1 levels9.7020.7180.0110.0090.272
Interleukin-12p70 levels8.7630.459-0.0060.0080.513
Interferon gamma-induced protein 10 levels4.6490.703-0.0070.0140.611
Interleukin-18 levels17.7010.169-0.0090.0110.436
Interleukin-17 levels7.2670.6090.0050.0150.753
Interleukin-13 levels8.1970.414-0.0010.0110.952
Interleukin-10 levels12.2210.201-0.0100.0100.350
Interleukin-8 levels6.7480.080-0.0010.0220.971
Interleukin-6 levels1.8840.5970.0260.0300.481
Interleukin-1-receptor antagonist levels3.4930.745-0.0010.0170.975
Interleukin-1-beta levels1.9150.7510.0110.0170.555
Hepatocyte growth factor levels5.2290.5150.0110.0160.534
Interleukin-9 levels2.3130.8040.0030.0230.912
Interleukin-7 levels9.0320.4340.0340.0180.088
Interleukin-5 levels3.3900.4950.0000.0200.986
Interleukin-4 levels10.6380.1550.0080.0170.662
Interleukin-2 receptor antagonist levels3.9910.6780.0020.0130.868
Interleukin-2 levels10.1790.336-0.0130.0110.244
Interferon gamma levels12.3650.1930.0040.0160.821
Growth-regulated protein alpha levels5.7110.6800.0060.0150.711
Granulocyte-colony stimulating factor levels1.6060.952-0.0040.0100.740
Fibroblast growth factor basic levels0.9500.9170.0150.0320.667
Eotaxin levels13.9440.5300.0050.0100.592

Heterogeneity and pleiotropy analyses for lung cancer.

Table 6

HeterogenityMR-Egger intercept
QQ_P -valueEgger_interceptSEP -value
CTACK levels5.2940.6240.0140.0220.557
beta-nerve growth factor levels11.8420.1060.0090.0520.876
Vascular endothelial growth factor levels5.9240.748-0.0150.0130.286
Macrophage Migration Inhibitory Factor levels2.9790.7030.0180.0310.605
TRAIL levels8.3240.7590.0140.0140.358
Tumor necrosis factor beta levels0.5120.916-0.0020.0190.908
Tumor necrosis factor alpha levels3.7730.438-0.0240.0210.336
Stromal-cell-derived factor 1 alpha levels5.9940.5400.0180.0190.395
Stem cell growth factor beta levels11.6840.554-0.0150.0150.343
Stem cell factor levels8.5460.3820.0120.0200.579
Interleukin-16 levels6.3370.6090.0050.0180.801
RANTES levels3.3560.7630.0230.0370.568
Platelet-derived growth factor BB levels7.8310.6450.0150.0140.319
Macrophage inflammatory protein 1b levels13.8300.6790.0050.0110.667
Macrophage inflammatory protein 1a levels6.4140.601-0.0120.0270.670
Monokine induced by gamma interferon levels8.5900.737-0.0240.0200.261
Macrophage colony stimulating factor levels9.9630.1910.0020.0330.943
Monocyte chemoattractant protein-3 levels1.8100.405-0.0050.0830.960
Monocyte chemoattractant protein-1 levels9.8480.5440.0190.0170.302
Interleukin-12p70 levels6.2450.7150.0040.0130.740
Interferon gamma-induced protein 10 levels6.8120.557-0.0170.0190.412
Interleukin-18 levels18.2260.197-0.0180.0160.286
Interleukin-17 levels2.0110.9910.0080.0230.739
Interleukin-13 levels7.2600.5090.0080.0160.625
Interleukin-10 levels13.2180.153-0.0130.0170.480
Interleukin-8 levels5.5430.1360.0260.0260.435
Interleukin-6 levels2.7060.4390.0540.0480.382
Interleukin-1-receptor antagonist levels1.2670.9730.0020.0270.950
Interleukin-1-beta levels2.7660.5980.0380.0270.252
Hepatocyte growth factor levels4.9320.553-0.0100.0260.723
Interleukin-9 levels3.1160.6820.0120.0370.767
Interleukin-7 levels6.3270.7070.0510.0280.105
Interleukin-5 levels2.3250.6760.0020.0300.955
Interleukin-4 levels6.4430.375-0.0170.0240.510
Interleukin-2 receptor antagonist levels3.3670.7620.0100.0210.649
Interleukin-2 levels9.6270.292-0.0230.0190.254
Interferon gamma levels10.7360.2940.0290.0210.206
Growth-regulated protein alpha levels14.3790.1090.0430.0270.150
Granulocyte-colony stimulating factor levels11.9280.103-0.0050.0230.826
Fibroblast growth factor basic levels1.8140.770-0.0090.0490.865
Eotaxin levels11.0770.7470.0130.0160.437

Heterogeneity and pleiotropy analyses for squamous cell lung cancer.

Table 7

HeterogenityMR-Egger intercept
QQ_P -valueEgger_interceptSEP -value
CTACK levels7.2510.2980.0000.0260.991
beta-nerve growth factor levels2.2630.894-0.0140.0350.710
Vascular endothelial growth factor levels7.8090.5530.0040.0120.737
Macrophage Migration Inhibitory Factor levels1.2140.944-0.0140.0270.635
TRAIL levels10.8480.6980.0160.0110.158
Tumor necrosis factor beta levels0.9390.816-0.0100.0160.597
Tumor necrosis factor alpha levels2.3480.6720.0060.0190.760
Stromal-cell-derived factor 1 alpha levels9.3720.227-0.0190.0190.351
Stem cell growth factor beta levels9.8630.453-0.0220.0150.193
Stem cell factor levels8.8670.3540.0140.0170.444
Interleukin-16 levels9.7750.369-0.0090.0150.560
RANTES levels12.2600.140-0.0670.0260.037
Platelet-derived growth factor BB levels16.5570.167-0.0040.0150.809
Macrophage inflammatory protein 1b levels18.7550.3430.0180.0100.081
Macrophage inflammatory protein 1a levels8.0290.431-0.0380.0230.147
Monokine induced by gamma interferon levels8.3940.817-0.0100.0180.579
Macrophage colony stimulating factor levels2.2210.947-0.0240.0230.323
Monocyte chemoattractant protein-3 levels2.4880.2880.0590.0630.519
Monocyte chemoattractant protein-1 levels12.5100.4860.0010.0130.969
Interleukin-12p70 levels9.3230.4080.0040.0120.753
Interferon gamma-induced protein 10 levels3.9880.912-0.0030.0160.865
Interleukin-18 levels11.4330.575-0.0020.0120.866
Interleukin-17 levels10.6390.3010.0100.0230.690
Interleukin-13 levels1.1020.954-0.0140.0170.456
Interleukin-10 levels7.2090.6150.0010.0120.915
Interleukin-8 levels4.8140.1860.0090.0250.746
Interleukin-6 levels6.7640.149-0.0200.0310.566
Interleukin-1-receptor antagonist levels5.7110.4560.0080.0250.759
Interleukin-1-beta levels1.9170.751-0.0190.0230.470
Hepatocyte growth factor levels7.1940.3030.0060.0270.829
Interleukin-9 levels3.7620.584-0.0080.0320.821
Interleukin-7 levels9.5040.3920.0340.0240.199
Interleukin-5 levels2.8660.5800.0070.0260.818
Interleukin-4 levels6.9220.4370.0160.0180.412
Interleukin-2 receptor antagonist levels7.8410.2500.0070.0220.782
Interleukin-2 levels9.3850.4030.0080.0150.585
Interferon gamma levels12.2720.198-0.0060.0220.795
Growth-regulated protein alpha levels10.0260.3480.0030.0230.914
Granulocyte-colony stimulating factor levels4.9510.666-0.0120.0140.451
Fibroblast growth factor basic levels7.3580.1180.0440.0640.545
Eotaxin levels19.7140.1830.0000.0160.991

Heterogeneity and pleiotropy analyses for lung adenocarcinoma.

Table 8

HeterogenityMR-Egger intercept
QQ_P -valueEgger_interceptSEP -value
CTACK levels9.0170.2510.0070.0410.878
beta-nerve growth factor levels9.2570.160-0.0440.0800.608
Vascular endothelial growth factor levels7.0420.5320.0140.0210.515
Macrophage Migration Inhibitory Factor levels0.3690.8310.0450.1120.759
TRAIL levels11.4740.404-0.0290.0200.170
Tumor necrosis factor beta levels2.8200.4200.0470.0300.252
Tumor necrosis factor alpha levels1.0580.7870.0320.0410.521
Stromal-cell-derived factor 1 alpha levels4.9720.547-0.0450.0300.202
Stem cell growth factor beta levels4.1800.939-0.0080.0300.803
Stem cell factor levels6.0040.5390.0040.0390.929
Interleukin-16 levels11.2300.189-0.0420.0280.176
RANTES levels9.8070.200-0.0380.0590.540
Platelet-derived growth factor BB levels15.5460.213-0.0330.0240.203
Macrophage inflammatory protein 1b levels10.7770.768-0.0060.0180.750
Macrophage inflammatory protein 1a levels6.6050.4710.0350.0440.455
Monokine induced by gamma interferon levels10.4740.313-0.0250.0370.510
Macrophage colony stimulating factor levels6.8620.1430.0400.0610.555
Monocyte chemoattractant protein-3 levels0.1790.672Not ApplicableNot ApplicableNot Applicable
Monocyte chemoattractant protein-1 levels11.5120.4860.0030.0270.921
Interleukin-12p70 levels3.9950.677-0.0310.0290.327
Interferon gamma-induced protein 10 levels8.6250.196-0.0180.0490.734
Interleukin-18 levels13.7700.088-0.0780.0370.076
Interleukin-17 levels6.1630.521-0.0140.0480.780
Interleukin-13 levels8.8580.263-0.0170.0300.593
Interleukin-10 levels10.3910.239-0.0410.0230.114
Interleukin-8 levels3.0660.382-0.0350.0310.378
Interleukin-6 levels0.4650.9270.0000.0400.999
Interleukin-1-receptor antagonist levels3.9970.677-0.0120.0420.778
Interleukin-1-beta levels5.5250.2370.0160.0560.796
Hepatocyte growth factor levels2.2410.8960.0470.0410.307
Interleukin-9 levels5.7820.216-0.0230.0830.799
Interleukin-7 levels11.5030.118-0.0140.0610.831
Interleukin-5 levels0.9320.334Not ApplicableNot ApplicableNot Applicable
Interleukin-4 levels2.0890.719-0.0150.0720.848
Interleukin-2 receptor antagonist levels4.0420.543-0.0200.0320.563
Interleukin-2 levels5.1600.640-0.0120.0260.663
Interferon gamma levels2.4160.878-0.0280.0320.427
Growth-regulated protein alpha levels3.7330.880-0.0430.0380.298
Granulocyte-colony stimulating factor levels2.6880.748-0.0300.0260.318
Fibroblast growth factor basic levels3.9150.2710.0750.1140.579
Eotaxin levels5.3650.966-0.0060.0270.821

Heterogeneity and pleiotropy analyses for small cell lung carcinoma.

Table 9

ExposureOutcomeDirectionSteiger P -value
CTACK levelsLung cancerTRUE4.98E-84
beta-nerve growth factor levelsLung cancerTRUE5.01E-46
Vascular endothelial growth factor levelsLung cancerTRUE1.36E-259
Macrophage Migration Inhibitory Factor levelsLung cancerTRUE8.47E-51
TRAIL levelsLung cancerTRUE0
Tumor necrosis factor beta levelsLung cancerTRUE5.71E-44
Tumor necrosis factor alpha levelsLung cancerTRUE1.59E-25
Stromal-cell-derived factor 1 alpha levelsLung cancerTRUE5.03E-34
Stem cell growth factor beta levelsLung cancerTRUE2.72E-68
Stem cell factor levelsLung cancerTRUE6.56E-58
Interleukin-16 levelsLung cancerTRUE1.81E-83
RANTES levelsLung cancerTRUE3.50E-60
Platelet-derived growth factor BB levelsLung cancerTRUE2.61E-114
Macrophage inflammatory protein 1b levelsLung cancerTRUE0
Macrophage inflammatory protein 1a levelsLung cancerTRUE2.20E-41
Monokine induced by gamma interferon levelsLung cancerTRUE8.51E-85
Macrophage colony stimulating factor levelsLung cancerTRUE2.97E-60
Monocyte chemoattractant protein-3 levelsLung cancerTRUE3.21E-28
Monocyte chemoattractant protein-1 levelsLung cancerTRUE1.13E-123
Interleukin-12p70 levelsLung cancerTRUE9.05E-301
Interferon gamma-induced protein 10 levelsLung cancerTRUE5.14E-49
Interleukin-18 levelsLung cancerTRUE2.64E-168
Interleukin-17 levelsLung cancerTRUE1.32E-55
Interleukin-13 levelsLung cancerTRUE2.99E-117
Interleukin-10 levelsLung cancerTRUE1.09E-171
Interleukin-8 levelsLung cancerTRUE6.92E-20
Interleukin-6 levelsLung cancerTRUE1.71E-36
Interleukin-1-receptor antagonist levelsLung cancerTRUE4.62E-43
Interleukin-1-beta levelsLung cancerTRUE2.34E-23
Hepatocyte growth factor levelsLung cancerTRUE6.94E-49
Interleukin-9 levelsLung cancerTRUE1.66E-37
Interleukin-7 levelsLung cancerTRUE6.40E-137
Interleukin-5 levelsLung cancerTRUE8.87E-32
Interleukin-4 levelsLung cancerTRUE7.50E-55
Interleukin-2 receptor antagonist levelsLung cancerTRUE2.15E-73
Interleukin-2 levelsLung cancerTRUE6.91E-75
Interferon gamma levelsLung cancerTRUE2.44E-62
Growth-regulated protein alpha levelsLung cancerTRUE2.94E-134
Granulocyte-colony stimulating factor levelsLung cancerTRUE1.63E-56
Fibroblast growth factor basic levelsLung cancerTRUE7.28E-36
Eotaxin levelsLung cancerTRUE7.78E-136

Direction test for lung cancer.

Table 10

ExposureOutcomeDirectionSteiger P -value
CTACK levelsSquamous cell lung cancerTRUE7.19E-88
beta-nerve growth factor levelsSquamous cell lung cancerTRUE5.27E-50
Vascular endothelial growth factor levelsSquamous cell lung cancerTRUE1.35E-249
Macrophage Migration Inhibitory Factor levelsSquamous cell lung cancerTRUE3.75E-51
TRAIL levelsSquamous cell lung cancerTRUE2.88E-305
Tumor necrosis factor beta levelsSquamous cell lung cancerTRUE1.40E-43
Tumor necrosis factor alpha levelsSquamous cell lung cancerTRUE2.16E-24
Stromal-cell-derived factor 1 alpha levelsSquamous cell lung cancerTRUE7.87E-32
Stem cell growth factor beta levelsSquamous cell lung cancerTRUE3.92E-117
Stem cell factor levelsSquamous cell lung cancerTRUE3.81E-56
Interleukin-16 levelsSquamous cell lung cancerTRUE1.30E-81
RANTES levelsSquamous cell lung cancerTRUE4.19E-49
Platelet-derived growth factor BB levelsSquamous cell lung cancerTRUE8.27E-126
Macrophage inflammatory protein 1b levelsSquamous cell lung cancerTRUE0
Macrophage inflammatory protein 1a levelsSquamous cell lung cancerTRUE1.52E-40
Monokine induced by gamma interferon levelsSquamous cell lung cancerTRUE1.61E-86
Macrophage colony stimulating factor levelsSquamous cell lung cancerTRUE6.56E-59
Monocyte chemoattractant protein-3 levelsSquamous cell lung cancerTRUE1.84E-27
Monocyte chemoattractant protein-1 levelsSquamous cell lung cancerTRUE5.41E-92
Interleukin-12p70 levelsSquamous cell lung cancerTRUE3.46E-292
Interferon gamma-induced protein 10 levelsSquamous cell lung cancerTRUE1.23E-56
Interleukin-18 levelsSquamous cell lung cancerTRUE8.72E-175
Interleukin-17 levelsSquamous cell lung cancerTRUE1.23E-54
Interleukin-13 levelsSquamous cell lung cancerTRUE6.12E-117
Interleukin-10 levelsSquamous cell lung cancerTRUE3.66E-164
Interleukin-8 levelsSquamous cell lung cancerTRUE1.82E-19
Interleukin-6 levelsSquamous cell lung cancerTRUE3.56E-34
Interleukin-1-receptor antagonist levelsSquamous cell lung cancerTRUE1.81E-42
Interleukin-1-beta levelsSquamous cell lung cancerTRUE5.46E-23
Hepatocyte growth factor levelsSquamous cell lung cancerTRUE6.73E-46
Interleukin-9 levelsSquamous cell lung cancerTRUE8.94E-36
Interleukin-7 levelsSquamous cell lung cancerTRUE1.19E-134
Interleukin-5 levelsSquamous cell lung cancerTRUE6.16E-32
Interleukin-4 levelsSquamous cell lung cancerTRUE9.47E-48
Interleukin-2 receptor antagonist levelsSquamous cell lung cancerTRUE1.16E-70
Interleukin-2 levelsSquamous cell lung cancerTRUE1.03E-63
Interferon gamma levelsSquamous cell lung cancerTRUE2.43E-60
Growth-regulated protein alpha levelsSquamous cell lung cancerTRUE2.47E-134
Granulocyte-colony stimulating factor levelsSquamous cell lung cancerTRUE1.86E-53
Fibroblast growth factor basic levelsSquamous cell lung cancerTRUE2.52E-34
Eotaxin levelsSquamous cell lung cancerTRUE2.44E-132

Direction test for squamous cell lung cancer.

Table 11

ExposureOutcomeDirectionSteiger P -value
CTACK levelsLung adenocarcinomaTRUE5.14E-77
beta-nerve growth factor levelsLung adenocarcinomaTRUE4.85E-45
Vascular endothelial growth factor levelsLung adenocarcinomaTRUE3.62E-252
Macrophage Migration Inhibitory Factor levelsLung adenocarcinomaTRUE5.14E-51
TRAIL levelsLung adenocarcinomaTRUE0
Tumor necrosis factor beta levelsLung adenocarcinomaTRUE8.71E-44
Tumor necrosis factor alpha levelsLung adenocarcinomaTRUE7.78E-25
Stromal-cell-derived factor 1 alpha levelsLung adenocarcinomaTRUE2.21E-31
Stem cell growth factor beta levelsLung adenocarcinomaTRUE2.98E-80
Stem cell factor levelsLung adenocarcinomaTRUE8.86E-56
Interleukin-16 levelsLung adenocarcinomaTRUE1.02E-88
RANTES levelsLung adenocarcinomaTRUE1.95E-57
Platelet-derived growth factor BB levelsLung adenocarcinomaTRUE1.29E-131
Macrophage inflammatory protein 1b levelsLung adenocarcinomaTRUE0
Macrophage inflammatory protein 1a levelsLung adenocarcinomaTRUE2.90E-40
Monokine induced by gamma interferon levelsLung adenocarcinomaTRUE6.90E-95
Macrophage colony stimulating factor levelsLung adenocarcinomaTRUE6.96E-60
Monocyte chemoattractant protein-3 levelsLung adenocarcinomaTRUE1.38E-27
Monocyte chemoattractant protein-1 levelsLung adenocarcinomaTRUE1.09E-118
Interleukin-12p70 levelsLung adenocarcinomaTRUE4.36E-292
Interferon gamma-induced protein 10 levelsLung adenocarcinomaTRUE7.43E-64
Interleukin-18 levelsLung adenocarcinomaTRUE7.81E-148
Interleukin-17 levelsLung adenocarcinomaTRUE5.41E-52
Interleukin-13 levelsLung adenocarcinomaTRUE2.06E-100
Interleukin-10 levelsLung adenocarcinomaTRUE1.53E-168
Interleukin-8 levelsLung adenocarcinomaTRUE9.66E-20
Interleukin-6 levelsLung adenocarcinomaTRUE6.65E-41
Interleukin-1-receptor antagonist levelsLung adenocarcinomaTRUE3.88E-41
Interleukin-1-beta levelsLung adenocarcinomaTRUE7.61E-23
Hepatocyte growth factor levelsLung adenocarcinomaTRUE1.80E-45
Interleukin-9 levelsLung adenocarcinomaTRUE1.12E-36
Interleukin-7 levelsLung adenocarcinomaTRUE1.17E-134
Interleukin-5 levelsLung adenocarcinomaTRUE1.99E-31
Interleukin-4 levelsLung adenocarcinomaTRUE1.87E-54
Interleukin-2 receptor antagonist levelsLung adenocarcinomaTRUE1.61E-70
Interleukin-2 levelsLung adenocarcinomaTRUE6.07E-68
Interferon gamma levelsLung adenocarcinomaTRUE1.79E-60
Growth-regulated protein alpha levelsLung adenocarcinomaTRUE8.76E-137
Granulocyte-colony stimulating factor levelsLung adenocarcinomaTRUE2.31E-54
Fibroblast growth factor basic levelsLung adenocarcinomaTRUE5.71E-32
Eotaxin levelsLung adenocarcinomaTRUE7.40E-130

Direction test for lung adenocarcinoma.

Table 12

ExposureOutcomeDirectionSteiger P -value
CTACK levelsSmall cell lung carcinomaTRUE7.50E-69
beta-nerve growth factor levelsSmall cell lung carcinomaTRUE3.14E-36
Vascular endothelial growth factor levelsSmall cell lung carcinomaTRUE1.02E-192
Macrophage Migration Inhibitory Factor levelsSmall cell lung carcinomaTRUE2.27E-32
TRAIL levelsSmall cell lung carcinomaTRUE4.2579E-238
Tumor necrosis factor beta levelsSmall cell lung carcinomaTRUE1.74E-34
Tumor necrosis factor alpha levelsSmall cell lung carcinomaTRUE4.40E-18
Stromal-cell-derived factor 1 alpha levelsSmall cell lung carcinomaTRUE1.97E-24
Stem cell growth factor beta levelsSmall cell lung carcinomaTRUE4.40E-66
Stem cell factor levelsSmall cell lung carcinomaTRUE1.77E-40
Interleukin-16 levelsSmall cell lung carcinomaTRUE3.26E-66
RANTES levelsSmall cell lung carcinomaTRUE2.70E-45
Platelet-derived growth factor BB levelsSmall cell lung carcinomaTRUE5.20E-94
Macrophage inflammatory protein 1b levelsSmall cell lung carcinomaTRUE0
Macrophage inflammatory protein 1a levelsSmall cell lung carcinomaTRUE3.25E-31
Monokine induced by gamma interferon levelsSmall cell lung carcinomaTRUE1.34E-58
Macrophage colony stimulating factor levelsSmall cell lung carcinomaTRUE3.30E-36
Monocyte chemoattractant protein-3 levelsSmall cell lung carcinomaTRUE1.28E-21
Monocyte chemoattractant protein-1 levelsSmall cell lung carcinomaTRUE5.72E-87
Interleukin-12p70 levelsSmall cell lung carcinomaTRUE7.14E-124
Interferon gamma-induced protein 10 levelsSmall cell lung carcinomaTRUE6.96E-41
Interleukin-18 levelsSmall cell lung carcinomaTRUE2.86E-72
Interleukin-17 levelsSmall cell lung carcinomaTRUE1.74E-33
Interleukin-13 levelsSmall cell lung carcinomaTRUE1.07E-86
Interleukin-10 levelsSmall cell lung carcinomaTRUE6.85E-126
Interleukin-8 levelsSmall cell lung carcinomaTRUE1.40E-17
Interleukin-6 levelsSmall cell lung carcinomaTRUE1.19E-29
Interleukin-1-receptor antagonist levelsSmall cell lung carcinomaTRUE2.03E-28
Interleukin-1-beta levelsSmall cell lung carcinomaTRUE1.65E-19
Hepatocyte growth factor levelsSmall cell lung carcinomaTRUE4.00E-35
Interleukin-9 levelsSmall cell lung carcinomaTRUE4.27E-27
Interleukin-7 levelsSmall cell lung carcinomaTRUE2.99E-99
Interleukin-5 levelsSmall cell lung carcinomaTRUE6.73E-12
Interleukin-4 levelsSmall cell lung carcinomaTRUE1.19E-29
Interleukin-2 receptor antagonist levelsSmall cell lung carcinomaTRUE1.61E-59
Interleukin-2 levelsSmall cell lung carcinomaTRUE1.52E-54
Interferon gamma levelsSmall cell lung carcinomaTRUE1.09E-39
Growth-regulated protein alpha levelsSmall cell lung carcinomaTRUE3.86E-120
Granulocyte-colony stimulating factor levelsSmall cell lung carcinomaTRUE1.10E-24
Fibroblast growth factor basic levelsSmall cell lung carcinomaTRUE1.75E-18
Eotaxin levelsSmall cell lung carcinomaTRUE2.07E-98

Direction test for small cell lung carcinoma.

Subgroup analyses

Additional subgroup analyses stratified by smoking status were conducted to investigate whether the causal effect of cytokines on lung cancer risk was modified by smoking. In ever smokers, we found that interferon gamma-induced protein 10 (IP-10) (OR = 0.861, 95% CI: 0.781-0.950, P = 0.003) was inversely associated with lung cancer risk, while IL-1β (OR = 1.190, 95% CI: 1.023-1.384, P = 0.024) was positively associated with lung cancer risk (Supplementary Tables 10, 11). Scatter plots visualized the effects for each method across datasets (Supplementary Figures 7, 8). Heterogeneity tests showed no significant heterogeneity across cytokines, except for tumor necrosis factor (TNF)-α (Supplementary Table 12). MR-Egger regression revealed no evidence of horizontal pleiotropy for any cytokine (Supplementary Table 12). The MR-PRESSO outlier test confirmed significant MR findings, except for TNF-β, TNF-α, and MCP-3, due to limited SNP availability (Supplementary Table 13). Leave-one-out analysis confirmed the robustness of the causal estimates (Supplementary Figures 9, 10). Furthermore, Steiger p-values <0.05 verified the correct causal direction for all detected cytokines (Supplementary Table 14).

Among never smokers, positive association was observed between lung cancer risk and SCF (OR = 1.474, 95% CI: 1.105-1.964, P = 0.008) (Supplementary Tables 15, 16). Scatter plots visualized the effects for each method across datasets (Supplementary Figure 11). Heterogeneity was non-significant across all cytokines (Supplementary Table 17). MR-Egger and MR-PRESSO global tests showed no evidence of pleiotropy (Supplementary Tables 17, 18). The MR-PRESSO outlier test confirmed significant MR findings, except for MCP-3, IL-8, and fibroblast growth factor basic (FGF-basic), due to limited SNPs (Supplementary Table 18). Leave-one-out analysis confirmed the robustness of the causal estimates (Supplementary Figures 12). Additionally, Steiger p-values <0.05 verified the expected causal direction for all detected cytokines (Supplementary Table 19).

Discussion

In this MR study, we investigated potential causal relationships between 41 circulating cytokines and lung cancer risk. Our analysis provides preliminary evidence that genetically predicted levels of certain cytokines may influence cancer susceptibility. We identified four inflammatory mediators - SCF, IL-1β, IL-18, and IP-10 - involved in genetic susceptibility to lung cancer overall and in specific histologic subtypes, as well as differences based on smoking status. These results illuminate cytokine pathways that may promote cancer development, representing potential targets for prevention, early detection, and treatment. Our findings reveal putative causal effects of circulating cytokines in lung cancer pathogenesis and shed light on cytokine-mediated mechanisms influencing susceptibility across lung cancer subtypes and smoking strata.

Eotaxin, a chemokine, plays a pivotal role in managing a spectrum of inflammatory and immune-responsive conditions, chiefly by recruiting and activating eosinophils (, ). However, its specific function in lung cancer remains to be elucidated. Multiple studies have shown that eotaxin expression levels correlate with occurrence and prognosis in several cancers. Specifically, Yamaguchi et al. found lower eotaxin levels in healthy controls versus colon cancer patients, while Melisi et al. showed decreased eotaxin-2 levels in pancreatic cancer patients after treatment with a TGF-β receptor inhibitor and gemcitabine (, ). Additionally, Siva et al. reported reduced serum eotaxin-1 in non-small cell lung cancer patients following radiotherapy compared to radiotherapy alone, and Tsao et al. revealed an association between low serum eotaxin-1 and shorter progression-free survival in non-small cell lung cancer patients on vandetanib (, ). The statistical power of eotaxin was calculated to be only 67% and had a non-significant p-value after Bonferroni correction. In summary, although eotaxin is closely linked to pathogenesis and prognosis in some cancers, further research is warranted to elucidate its precise role in lung cancer. SCF activates the c-Kit signaling pathway to stimulate lung cancer cell proliferation, migration and invasion (). In vitro experiments show SCF promotes lung cancer cell proliferation, migration and metastasis (–). Moreover, clinical studies find elevated SCF and c-Kit expression correlates with lung cancer progression and metastasis (, ). IL-1β advances lung cancer progression through numerous mechanisms including inducing angiogenic factors like VEGF, activating oncogenic signaling, and promoting immunosuppression and inflammation (–). Mechanistically, IL-1β spurs tumor growth by activating MAPK/NF-κB pathways, recruiting immunosuppressive cells, promoting inflammation, angiogenesis, invasion and metastasis, downregulating tumor suppressors, upregulating oncogenes, and enabling immune evasion and inhibition of apoptosis (). Preclinical studies also demonstrate beneficial effects of IL-1β knockdown and inhibition. No inflammation or angiogenesis occurred in IL-1β-deficient mice, and IL-1β-deficient mice showed no local tumor growth or lung metastasis compared to wild-type mice (). Moreover, IL-1β antibody treatment inhibited tumor progression and boosted antitumor immunity in mice by reducing inflammation and promoting M1 macrophage maturation (). VEGF is a pivotal driver of angiogenesis in lung cancer, promoting tumor growth and metastasis (). Furthermore, IL-1β induces VEGF secretion, which also spurs tumor expansion and spread (). In lung cancer, VEGF overexpression critically supports angiogenesis and correlates with disease progression and prognosis (). Multiple studies demonstrate VEGF inhibition slows lung cancer growth and improves chemotherapy efficacy (, ). Anti-VEGF monoclonal antibody therapies including bevacizumab are now widely used in lung cancer treatment (53). Similarly, the calculated statistical power of VEGF is only 21%.

In contrast, increased genetic predisposition to higher interleukin-18 (IL-18) and interferon gamma-induced protein 10 (IP-10) levels are negatively associated with lung cancer risk. IL-18, a pro-inflammatory cytokine, plays a significant role in anti-tumor immunity through several mechanisms. Firstly, IL-18 activates natural killer (NK) cells and certain subsets of T cells, enhancing their cytotoxic activity against tumor cells by upregulating cytotoxic molecules such as perforin and granzymes, which are crucial for direct tumor cell killing (54). Secondly, IL-18 stimulates the production of key anti-tumor cytokines, notably interferon-gamma (IFN-γ), which has potent anti-tumor effects, including the inhibition of tumor cell proliferation, induction of tumor cell apoptosis, and enhancement of antigen presentation to help the immune system better recognize and target tumor cells (55). Thirdly, IL-18 facilitates the infiltration of immune cells, particularly T cells and NK cells, into the tumor microenvironment, thereby enhancing the overall immune response against the tumor (56). Additionally, IL-18 promotes the differentiation of T cells towards a Th1 phenotype; Th1 cells produce IFN-γ and activate macrophages and other immune cells that contribute to anti-tumor immunity, making this shift towards a Th1-dominated response crucial for effective anti-tumor activity (57, 58). Furthermore, IL-18 often works synergistically with other cytokines and immune signals to amplify the immune response, significantly boosting IFN-γ production and enhancing the cytotoxic activity of NK cells and T cells when combined with IL-12 (59, 60). IL-18 also exhibits anti-angiogenic and pro-lymphangiogenic properties that contribute to its anti-tumor activity (57). However, Jiang et al. found that IL-18 may promote metastasis by inhibiting E-cadherin expression (58). Conversely, experimental studies by Xiong et al. and Chen et al. demonstrated that IL-18 inhibited tumor proliferation and growth, enhanced apoptosis, and normalized the Th1/Th2 imbalance (59, 60).

In contrast, increased genetic predisposition to higher levels of interleukin-18 (IL-18) and interferon gamma-induced protein 10 (IP-10) is negatively associated with lung cancer risk. IL-18, a pro-inflammatory cytokine, plays a significant role in anti-tumor immunity through several mechanisms. Firstly, IL-18 activates natural killer (NK) cells and certain subsets of T cells, enhancing their cytotoxic activity against tumor cells by upregulating cytotoxic molecules such as perforin and granzymes, which are crucial for direct tumor cell killing (54). Secondly, IL-18 stimulates the production of key anti-tumor cytokines, notably interferon-gamma (IFN-γ), which has potent anti-tumor effects, including the inhibition of tumor cell proliferation, induction of tumor cell apoptosis, and enhancement of antigen presentation to help the immune system better recognize and target tumor cells (55). Thirdly, IL-18 facilitates the infiltration of immune cells, particularly T cells and NK cells, into the tumor microenvironment, thereby enhancing the overall immune response against the tumor (56). Additionally, IL-18 promotes the differentiation of T cells towards a Th1 phenotype; Th1 cells produce IFN-γ and activate macrophages and other immune cells that contribute to anti-tumor immunity, making this shift towards a Th1-dominated response crucial for effective anti-tumor activity (57, 58). Finally, IL-18 often works synergistically with other cytokines and immune signals to amplify the immune response, significantly boosting IFN-γ production and enhancing the cytotoxic activity of NK cells and T cells when combined with IL-12 (59, 60). IL-18 also has anti-angiogenic and pro-lymphangiogenic properties contributing to its anti-tumor activity (61). However, Jiang et al. found that IL-18 may promote metastasis by inhibiting E-cadherin expression (62). Conversely, experimental studies by Xiong et al. and Chen et al. demonstrated that IL-18 inhibited tumor proliferation and growth, enhanced apoptosis, and normalized the Th1/Th2 imbalance (63, 64). IP-10 attracts immune cells to tumors, inhibits angiogenesis, and reduces tumor burden (65). Previous studies confirm IP-10 participates in anti-tumor immunity by promoting immune cell migration to tumors, inducing apoptosis, and suppressing angiogenesis (66–69). In xenograft models of lymphoma, squamous cell carcinoma, and lung adenocarcinoma, CXCL10 production negatively correlated with tumor growth and significantly reduced tumor-associated angiogenesis (70). In advanced endometrial cancer, CXCL10 was shown to antagonize fibroblast growth factor action, thereby inhibiting angiogenesis. In estrogen receptor-positive breast tumors, CXCL10 inhibits vascular endothelial growth factor levels to reduce tumor burden (71).

We utilized Mendelian randomization to evaluate associations between circulating cytokines and lung cancer risk. Contrary to a prior study by Bouras et al. showing positive CTACK-nonsmoking lung cancer and negative IL-18-lung cancer/adenocarcinoma relationships (72), we found no evidence for CTACK-nonsmoking lung cancer or IL-18-adenocarcinoma associations, potentially attributable to divergent instrument variable selection and GWAS data pooling. Notably, our analysis revealed novel causal links between SCF, IL-1β, IL-18, and IP-10 in overall lung cancer as well as specific histological subtypes, highlighting important etiological roles for these cytokines. However, some limitations should be considered. First, the relaxed IV significance threshold of P < 5×10-6 introduces possible false positives and bias, although the consistent F-statistics >10 suggest weak instrument bias is less likely. Second, the single Finnish ethnicity limits generalizability to other populations. Third, no cytokines were statistically significantly associated with cancer risk or subtypes after Bonferroni correction, including six inflammatory factors with suggestive correlations - Eotaxin, SCF, IL-1β, VEGF, IL-18, and IP-10. However, excluding Eotaxin and VEGF, the other factors (SCF, IL-1β, IL-18, and IP-10) had statistical power over 80% but still require validation of these potential associations in larger cohorts and GWAS. Fourth, while efforts were made to mitigate confounding, pleiotropy cannot be completely ruled out. Finally, while not addressed here, inflammatory factors may influence lung cancer progression and survival rather than development. Therefore, further studies should analyze the role of inflammatory factors in lung cancer aggressiveness.

Conclusion

This MR study found preliminary evidence that genetically predicted levels of four inflammatory cytokines—SCF, IL-1β, IL-18, and IP-10—may causally influence lung cancer risk overall, in specific histologic subtypes, and stratified by smoking status. The identification of these cytokine pathways, which may promote lung carcinogenesis, represents potential new targets for the prevention, early detection, and treatment of lung cancer. Overall, these findings reveal putative causal effects of circulating cytokines in lung cancer pathogenesis, illuminating cytokine-mediated immunological mechanisms affecting susceptibility across lung cancer subtypes and smoking exposure groups.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author contributions

DL: Data curation, Formal analysis, Writing – original draft. ZG: Data curation, Formal analysis, Writing – original draft. QZ: Data curation, Formal analysis, Writing – original draft, Writing – review & editing. SL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Key Project of the Affiliated Hospital of North Sichuan Medical College (2023ZD008).

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/fonc.2024.1373380/full#supplementary-material

Supplementary Figure 1

Scatter plot of IL-18 levels on lung cancer.

Supplementary Figure 2

Scatter plot of SCF levels on lung lung adenocarcinoma.

Supplementary Figure 3

Scatter plot of IL-1β levels on lung lung adenocarcinoma.

Supplementary Figure 4

Leave-one-out plot of IL-1β levels on lung cancer.

Supplementary Figure 5

Leave-one-out plot of SCF levels on lung adenocarcinoma.

Supplementary Figure 6

Leave-one-out plot of IL-1β levels on lung adenocarcinoma.

Supplementary Figure 7

Scatter plot of IP-10 levels on lung cancer in ever smokers.

Supplementary Figure 8

Scatter plot of IL-1β levels on lung cancer in ever smokers.

Supplementary Figure 9

Leave-one-out plot of IP-10 levels on lung cancer in ever smokers.

Supplementary Figure 10

Leave-one-out plot of IL-1β levels on lung cancer in ever smokers.

Supplementary Figure 11

Scatter plot of SCF levels on lung cancer in never smokers.

Supplementary Figure 12

Leave-one-out plot of SCF levels on lung cancer in never smokers.

Abbreviations

GWASgenome-wide association studies
SNPssingle nucleotide polymorphisms
IVsinstrumental variables
ORodds ratio
CIconfidence interval
IVWinverse variance weighting
MR-PRESSOMR pleiotropy residual sum and outlier
ICDInternational Classification of Diseases
MRmendelian randomization
beta-NGFbeta nerve growth factor
CTACKcutaneous T-cell attracting (CCL27)
FGF-basicfibroblast growth factor basic
G-CSFgranulocyte colony-stimulating factor
GRO-alphagrowth regulated oncogene-alpha
HGFhepatocyte growth factor
IFN-gammainterferon-gamma
IL-1rainterleukin-1 receptor antagonist
IL-1 betainterleukin-1 beta
IL-2interleukin-2
IL-2rainterleukin-2 receptor antagonist
IL-4interleukin-4
IL-5interleukin-5
IL-6interleukin-6
IL-7interleukin-7
IL-8interleukin-8
IL-9interleukin-9
IL-10interleukin-10
IL-12p70interleukin-12p70
IL-13interleukin-13
IL-16interleukin-16
IL-17interleukin-17
IL-18interleukin-18
IP-10interferon gamma-induced protein 10
MCP-1monocyte chemoattractant protein-1
MCP-3monocyte chemoattractant protein-3
M-CSFmacrophage colony-stimulating factor
MIFmacrophage migration inhibitory factor
MIGmonokine induced by gamma interferon
MIP-1amacrophage inflammatory protein 1a
MIP-1bmacrophage inflammatory protein 1b
PDGF-bbplatelet-derived growth factor BB
RANTESregulated on activation, normal T-cell expressed and secreted (CCL5)
SCFstem cell factor
SCGF-betastem cell growth factor beta
SDF-1 alphastromal-cell-derived factor 1 alpha
TNF-alphatumor necrosis factor-alpha
TNF-betatumor necrosis factor-beta
TRAILTNF-related apoptosis inducing ligand
VEGFvascular endothelial growth factor

References

  • 1

    BrayFLaversanneMSungHFerlayJSiegelRLSoerjomataramIet al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2024) 74:229–63. doi: 10.3322/caac.21834

  • 2

    AllemaniCMatsudaTDi CarloVHarewoodRMatzMNikšićMet al. Global surveillance of trends in cancer survival 2000–14 (CONCORD-3): analysis of individual records for 37 513 025 patients diagnosed with one of 18 cancers from 322 population-based registries in 71 countries. Lancet. (2018) 391:1023–75. doi: 10.1016/S0140-6736(17)33326-3

  • 3

    Cannon-AlbrightLACarrSRAkerleyW. Population-based relative risks for lung cancer based on complete family history of lung cancer. J Thorac Oncol. (2019) 14:1184–91. doi: 10.1016/j.jtho.2019.04.019

  • 4

    LarssonSCCarterPKarSVithayathilMMasonAMMichaëlssonKet al. Smoking, alcohol consumption, and cancer: A mendelian randomisation study in UK Biobank and international genetic consortia participants. PloS Med. (2020) 17:e1003178. doi: 10.1371/journal.pmed.1003178

  • 5

    ElinavENowarskiRThaissCAHuBJinCFlavellRA. Inflammation-induced cancer: crosstalk between tumours, immune cells and microorganisms. Nat Rev Cancer. (2013) 13:759–71. doi: 10.1038/nrc3611

  • 6

    AggarwalBBShishodiaSSandurSKPandeyMKSethiG. Inflammation and cancer: how hot is the link? Biochem Pharmacol. (2006) 72:1605–21. doi: 10.1016/j.bcp.2006.06.029

  • 7

    CandidoJHagemannT. Cancer-related inflammation. J Clin Immunol. (2013) 33 Suppl 1:S79–84. doi: 10.1007/s10875-012-9847-0

  • 8

    QuXTangYHuaS. Immunological approaches towards cancer and inflammation: A cross talk. Front Immunol. (2018) 9:563. doi: 10.3389/fimmu.2018.00563

  • 9

    MantovaniAAllavenaPSicaABalkwillF. Cancer-related inflammation. Nature. (2008) 454:436–44. doi: 10.1038/nature07205

  • 10

    AggarwalBBVijayalekshmiRVSungB. Targeting inflammatory pathways for prevention and therapy of cancer: short-term friend, long-term foe. Clin Cancer Res. (2009) 15:425–30. doi: 10.1158/1078-0432.CCR-08-0149

  • 11

    YuHKortylewskiMPardollD. Crosstalk between cancer and immune cells: role of STAT3 in the tumour microenvironment. Nat Rev Immunol. (2007) 7:41–51. doi: 10.1038/nri1995

  • 12

    ZouW. Immunosuppressive networks in the tumour environment and their therapeutic relevance. Nat Rev Cancer. (2005) 5:263–74. doi: 10.1038/nrc1586

  • 13

    BurgessSButterworthAThompsonSG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. (2013) 37:658–65. doi: 10.1002/gepi.21758

  • 14

    EmdinCAKheraAVKathiresanS. Mendelian randomization. JAMA. (2017) 318:1925–6. doi: 10.1001/jama.2017.17219

  • 15

    Davey SmithGHolmesMVDaviesNMEbrahimS. Mendel's laws, Mendelian randomization and causal inference in observational data: substantive and nomenclatural issues. Eur J Epidemiol. (2020) 35:99–111. doi: 10.1007/s10654-020-00622-7

  • 16

    Ahola-OlliAVWürtzPHavulinnaASAaltoKPitkänenNLehtimäkiTet al. Genome-wide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors. Am J Hum Genet. (2017) 100:40–50. doi: 10.1016/j.ajhg.2016.11.007

  • 17

    McKayJDHungRJHanYZongXCarreras-TorresRChristianiDCet al. Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat Genet. (2017) 49:1126–32. doi: 10.1038/ng.3892

  • 18

    LuoJle CessieSBlauwGJFranceschiCNoordamRvan HeemstD. Systemic inflammatory markers in relation to cognitive function and measures of brain atrophy: a Mendelian randomization study. Geroscience. (2022) 44:2259–70. doi: 10.1007/s11357-022-00602-7

  • 19

    LinSMaoXHeW. Causal association of circulating cytokines with sepsis: a Mendelian randomization study. Front Immunol. (2023) 14:1281845. doi: 10.3389/fimmu.2023.1281845

  • 20

    KamatMABlackshawJAYoungRSurendranPBurgessSDaneshJet al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. (2019) 35:4851–3. doi: 10.1093/bioinformatics/btz469

  • 21

    BurgessSThompsonSGCRP CHD Genetics Collaboration. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. (2011) 40:755–64. doi: 10.1093/ije/dyr036

  • 22

    BowdenJDel GrecoMFMinelliCZhaoQLawlorDASheehanNAet al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. (2019) 48:728–42. doi: 10.1093/ije/dyy258

  • 23

    BowdenJDavey SmithGHaycockPCBurgessS. 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

  • 24

    BowdenJDavey SmithGBurgessS. 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

  • 25

    VerbanckMChenCYNealeBDoR. 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

  • 26

    BurgessSThompsonSG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. (2017) 32:377–89. doi: 10.1007/s10654-017-0255-x

  • 27

    HemaniGTillingKDavey SmithG. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PloS Genet. (2017) 13:e1007081. doi: 10.1371/journal.pgen.1007081

  • 28

    LarssonSCTraylorMMalikRDichgansMBurgessSMarkusHSet al. Modifiable pathways in Alzheimer's disease: Mendelian randomisation analysis. BMJ. (2017) 359:j5375. doi: 10.1136/bmj.j5375

  • 29

    LiuBLyuLZhouWSongJYeDMaoYet al. Associations of the circulating levels of cytokines with risk of amyotrophic lateral sclerosis: a Mendelian randomization study. BMC Med. (2023) 21:39. doi: 10.1186/s12916-023-02736-7

  • 30

    DaughertyBLSicilianoSJDeMartinoJAMalkowitzLSirotinaASpringerMS. Cloning, expression, and characterization of the human eosinophil eotaxin receptor. J Exp Med. (1996) 183:2349–54. doi: 10.1084/jem.183.5.2349

  • 31

    WilliamsTJ. Eotaxin-1 (CCL11). Front Immunol. (2015) 6:84. doi: 10.3389/fimmu.2015.00084

  • 32

    YamaguchiMOkamuraSYamajiTIwasakiMTsuganeSShettyVet al. Plasma cytokine levels and the presence of colorectal cancer. PloS One. (2019) 14:e0213602. doi: 10.1371/journal.pone.0213602

  • 33

    MelisiDGarcia-CarboneroRMacarullaTPezetDDeplanqueGFuchsMet al. TGFβ receptor inhibitor galunisertib is linked to inflammation- and remodeling-related proteins in patients with pancreatic cancer. Cancer Chemother Pharmacol. (2019) 83:975–91. doi: 10.1007/s00280-019-03807-4

  • 34

    SivaSMacManusMKronTBestNSmithJLobachevskyPet al. A pattern of early radiation-induced inflammatory cytokine expression is associated with lung toxicity in patients with non-small cell lung cancer. PloS One. (2014) 9:e109560. doi: 10.1371/journal.pone.0109560

  • 35

    TsaoASLiuSLeeJJAldenCMBlumenscheinGRJrHerbstRet al. Clinical and biomarker outcomes of the phase II vandetanib study from the BATTLE trial. J Thorac Oncol. (2013) 8:658–61. doi: 10.1097/JTO.0b013e31828d08ae

  • 36

    MartinFHSuggsSVLangleyKELuHSTingJOkinoKHet al. Primary structure and functional expression of rat and human stem cell factor DNAs. Cell. (1990) 63:203–11. doi: 10.1016/0092-8674(90)90301-T

  • 37

    MatsuiJYamamotoYFunahashiYTsuruokaAWatanabeTWakabayashiTet al. E7080, a novel inhibitor that targets multiple kinases, has potent antitumor activities against stem cell factor producing human small cell lung cancer H146, based on angiogenesis inhibition. Int J Cancer. (2008) 122:664–71. doi: 10.1002/ijc.23131

  • 38

    LevinaVMarrangoniAWangTParikhSSuYHerbermanRet al. Elimination of human lung cancer stem cells through targeting of the stem cell factor-c-kit autocrine signaling loop. Cancer Res. (2010) 70:338–46. doi: 10.1158/0008-5472.CAN-09-1102

  • 39

    MaulikGBhartiAKhanEBroderickRJKijimaTSalgiaR. Modulation of c-Kit/SCF pathway leads to alterations in topoisomerase-I activity in small cell lung cancer. J Environ Pathol Toxicol Oncol. (2004) 23:237–51. doi: 10.1615/JEnvPathToxOncol.v23.i4

  • 40

    PietschTNicotraMRFraioliRWolfHKMottoleseMNataliPG. Expression of the c-Kit receptor and its ligand SCF in non-small-cell lung carcinomas. Int J Cancer. (1998) 75:171–5. doi: 10.1002/(ISSN)1097-0215

  • 41

    MickePBasraiMFaldumABittingerFRönnstrandLBlaukatAet al. Characterization of c-kit expression in small cell lung cancer: prognostic and therapeutic implications. Clin Cancer Res. (2003) 9:188–94.

  • 42

    SharmaSYangSCZhuLReckampKGardnerBBaratelliFet al. Tumor cyclooxygenase-2/prostaglandin E2-dependent promotion of FOXP3 expression and CD4+ CD25+ T regulatory cell activities in lung cancer. Cancer Res. (2005) 65:5211–20. doi: 10.1158/0008-5472.CAN-05-0141

  • 43

    GaronEBChih-Hsin YangJDubinettSM. The role of interleukin 1β in the pathogenesis of lung cancer. JTO Clin Res Rep. (2020) 1:100001. doi: 10.1016/j.jtocrr.2020.100001

  • 44

    GreenDR. The coming decade of cell death research: five riddles. Cell. (2019) 177:1094–107. doi: 10.1016/j.cell.2019.04.024

  • 45

    WongJMagunBEWoodLJ. Lung inflammation caused by inhaled toxicants: a review. Int J Chron Obstruct Pulmon Dis. (2016) 11:1391–401. doi: 10.2147/COPD

  • 46

    DinarelloCA. Why not treat human cancer with interleukin-1 blockade? Cancer Metastasis Rev. (2010) 29:317–29. doi: 10.1007/s10555-010-9229-0

  • 47

    CarmiYDotanSRiderPKaplanovIWhiteMRBaronRet al. The role of IL-1β in the early tumor cell-induced angiogenic response. J Immunol. (2013) 190:3500–9. doi: 10.4049/jimmunol.1202769

  • 48

    HanahanDWeinbergRA. Hallmarks of cancer: the next generation. Cell. (2011) 144:646–74. doi: 10.1016/j.cell.2011.02.013

  • 49

    VoronovECarmiYApteRN. The role IL-1 in tumor-mediated angiogenesis. Front Physiol. (2014) 5:114. doi: 10.3389/fphys.2014.00114

  • 50

    KajdaniukDMarekBFoltynWKos-KudłaB. Vascular endothelial growth factor (VEGF) - part 2: in endocrinology and oncology. Endokrynol Pol. (2011) 62:456–64.

  • 51

    MaJWaxmanDJ. Combination of antiangiogenesis with chemotherapy for more effective cancer treatment. Mol Cancer Ther. (2008) 7:3670–84. doi: 10.1158/1535-7163.MCT-08-0715

  • 52

    RosenLS. VEGF-targeted therapy: therapeutic potential and recent advances. Oncologist. (2005) 10:382–91. doi: 10.1634/theoncologist.10-6-382

  • 53

    FrezzettiDGalloMMaielloMRD'AlessioAEspositoCChicchinelliNet al. VEGF as a potential target in lung cancer. Expert Opin Ther Targets. (2017) 21:959–66. doi: 10.1080/14728222.2017.1371137

  • 54

    DinarelloCA. Interleukin-18 and the pathogenesis of inflammatory diseases. Semin Nephrol. (2007) 27:98–114. doi: 10.1016/j.semnephrol.2006.09.013

  • 55

    CastroFCardosoAPGonçalvesRMSerreKOliveiraMJ. Interferon-gamma at the crossroads of tumor immune surveillance or evasion. Front Immunol. (2018) 9:847. doi: 10.3389/fimmu.2018.00847

  • 56

    WongJLBerkEEdwardsRPKalinskiP. IL-18-primed helper NK cells collaborate with dendritic cells to promote recruitment of effector CD8+ T cells to the tumor microenvironment. Cancer Res. (2013) 73:4653–62. doi: 10.1158/0008-5472.CAN-12-4366

  • 57

    EsmailbeigMGhaderiA. Interleukin-18: a regulator of cancer and autoimmune diseases. Eur Cytokine Netw. (2017) 28:127–40. doi: 10.1684/ecn.2018.0401

  • 58

    LiuYGuYHanYZhangQJiangZZhangXet al. Tumor exosomal RNAs promote lung pre-metastatic niche formation by activating alveolar epithelial TLR3 to recruit neutrophils. Cancer Cell. (2016) 30:243–56. doi: 10.1016/j.ccell.2016.06.021

  • 59

    YoshimotoTOkamuraHTagawaYIIwakuraYNakanishiK. Interleukin 18 together with interleukin 12 inhibits IgE production by induction of interferon-gamma production from activated B cells. Proc Natl Acad Sci U.S.A. (1997) 94:3948–53. doi: 10.1073/pnas.94.8.3948

  • 60

    WeinstockJVBlumAMetwaliAElliottDArsenescuR. IL-18 and IL-12 signal through the NF-kappa B pathway to induce NK-1R expression on T cells. J Immunol. (2003) 170:5003–7. doi: 10.4049/jimmunol.170.10.5003

  • 61

    ShenJChoyDFYoshidaTet al. Interleukin-18 has antipermeablity and antiangiogenic activities in the eye: reciprocal suppression with VEGF. J Cell Physiol. (2014) 229:974–83. doi: 10.1002/jcp.24575

  • 62

    JiangDYingWLuYet al. Identification of metastasis-associated proteins by proteomic analysis and functional exploration of interleukin-18 in metastasis. Proteomics. (2003) 3:724–37. doi: 10.1002/pmic.200300411

  • 63

    XiongDFengRYangSLinTChenX. Antitumor activity of interleukin-18 on A549 human lung cancer cell line. J Cancer Res Ther. (2019) 15:1635–41. doi: 10.4103/jcrt.JCRT_226_19

  • 64

    ChenXFengRXiongDYangSLinT. Effect of lentiviral vector-packaged interleukin-18 gene on the Malignant behavior of lung cancer. Exp Ther Med. (2020) 19:319–26. doi: 10.3892/etm.2019.8204

  • 65

    DufourJHDziejmanMLiuMTLeungJHLaneTELusterAD. IFN-gamma-inducible protein 10 (IP-10; CXCL10)-deficient mice reveal a role for IP-10 in effector T cell generation and trafficking. J Immunol. (2002) 168:3195–204. doi: 10.4049/jimmunol.168.7.3195

  • 66

    LuXLJiangXBLiuREZhangSM. The enhanced anti-angiogenic and antitumor effects of combining flk1-based DNA vaccine and IP-10. Vaccine. (2008) 26:5352–7. doi: 10.1016/j.vaccine.2008.08.012

  • 67

    YangXChuYWangYZhangRXiongS. Targeted in vivo expression of IFN-gamma-inducible protein 10 induces specific antitumor activity. J Leukoc Biol. (2006) 80:1434–44. doi: 10.1189/jlb.0306212

  • 68

    PersanoLCrescenziMIndraccoloS. Anti-angiogenic gene therapy of cancer: current status and future prospects. Mol Aspects Med. (2007) 28:87–114. doi: 10.1016/j.mam.2006.12.005

  • 69

    TokunagaRZhangWNaseemMet al. CXCL9, CXCL10, CXCL11/CXCR3 axis for immune activation - A target for novel cancer therapy. Cancer Treat Rev. (2018) 63:40–7. doi: 10.1016/j.ctrv.2017.11.007

  • 70

    LiuMGuoSHibbertJMet al. CXCL10/IP-10 in infectious diseases pathogenesis and potential therapeutic implications. Cytokine Growth Factor Rev. (2011) 22:121–30. doi: 10.1016/j.cytogfr.2011.06.001

  • 71

    AronicaSMRaiberLHanzlyMKiselaC. Antitumor/antiestrogenic effect of the chemokine interferon inducible protein 10 (IP-10) involves suppression of VEGF expression in mammary tissue. J Interferon Cytokine Res. (2009) 29:83–92. doi: 10.1089/jir.2008.0034

  • 72

    BourasEKarhunenVGillDet al. Circulating inflammatory cytokines and risk of five cancers: a Mendelian randomization analysis. BMC Med. (2022) 20:3. doi: 10.1186/s12916-021-02193-0

Summary

Keywords

cytokines, lung cancer, Mendelian randomization, causality, genome-wide association study

Citation

Luo D, Gong Z, Zhan Q and Lin S (2024) Causal association of circulating cytokines with the risk of lung cancer: a Mendelian randomization study. Front. Oncol. 14:1373380. doi: 10.3389/fonc.2024.1373380

Received

19 January 2024

Accepted

03 June 2024

Published

18 June 2024

Volume

14 - 2024

Edited by

Amancio Carnero, Sevilla University, Spain

Reviewed by

James K. Fields, Johns Hopkins University, United States

Maria V. Guijarro, University of Florida, United States

Updates

Copyright

*Correspondence: Shan Lin, ; Qingyuan Zhan,

†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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics