A multiethnic bidirectional Mendelian randomization study negates causal effects of C-reactive protein concentrations on lung cancer.

BACKGROUND: Significantly rising plasma circulating C-reaction protein (CRP) concentrations are pervasive in lung cancer (LC) development, demonstrating a bidirectional relation. However, it remains uncertain whether the causation between them exists, and the degree to which the effect varies across different ethnic ancestries remains unknown. Therefore, we attempted to investigate the causal relationship between these two phenotypes.
METHODS: With summary statistics of CRP-related single nucleotide polymorphisms (SNPs) identified by several large-scale genome-wide association studies (GWAS) datasets based on five ethnic ancestries coverage worldwide, we implemented bidirectional two-sample Mendelian randomization (MR) analyses. Genetic summary data of 11,348 LC cases and 15,861 controls from the International Lung Cancer Consortium (ILCCO) were applied. The inverse-variance weighted (IVW) approach was utilized as the principal analysis, supplemented by various complementary methods.
RESULTS: MR study did not reveal the causal relationship shared across genetically predisposed CRP blood concentrations and LC risk (OR =1.022, 95% CI: 0.965-1.083, P=0.455) including pathological subtypes (OR =1.026, 95% CI: 0.947-1.112, P=0.534 for lung adenocarcinoma; OR =1.060, 95% CI: 0.970-1.158, P=0.201 for squamous cell lung cancer). Further analyses among East Asian, Hispanic/Latin American, European, African American/Afro-Caribbean, and South Asian populations revealed consistent null causation. Additionally, the causal effects of LC on CRP concentrations were not statically significant (OR =0.999, 95% CI: 0.977-1.021, P=0.923).
CONCLUSIONS: We did not observe a bidirectional causal association between CRP blood concentrations on LC among East Asian, Hispanic/Latin American, European, African American/Afro-Caribbean, and South Asian ancestry individuals. 2021 Translational Lung Cancer Research. All rights reserved.

Introduction

Cancer is a leading global public health issue and the second major cause of death in the United States. Particularly, lung cancer (LC) is currently the second most common cancer and the major cause of cancer-related death globally (1). The American Cancer Society estimates that 235,760 new LC cases and 131,880 LC deaths are projected to occur in the United States in 2021 (2). Thus, early recognition of potentially modifiable risk factors contributes to better prevention of LC.

Growing studies have demonstrated that chronic inflammation, especially within the respiratory system, might be a risk factor for LC (3-6). Several inflammatory conditions, such as chronic pulmonary infection (7) and chronic obstructive pulmonary disease (8), were correlated with an increased risk of LC. It has been proposed that LC could be induced by pulmonary inflammation through several pathways, such as the increase in angiogenesis during the repair of the damaged tissue and the generation of reactive oxygen (9). Simultaneously, circulating C-reactive protein (CRP) was considered as a systemic marker of chronic inflammation (10), and elevated CRP blood concentrations have been widely evaluated for their association with LC (3,11-20). Up till now, two assumptions have been presented to explain the mechanisms between elevated CRP concentrations and the risk of LC (10). On the one hand, elevated CRP concentrations and chronic inflammation are causally related to the carcinogenesis of the lung indeed. On the other hand, CRP concentrations may also indicate a systemic response to ongoing disease or comorbidity, such as an occult cancer or a premalignant state. The majority of previous cohort studies and a meta-analysis (21) denoted that elevated CRP concentrations were correlated with LC occurrence. However, given that cigarette smoking is by far the leading risk factor for LC, CRP is widely seen as a predictor of LC in smokers and former smokers (22), while cigarette smoking itself has been reported to increase circulating levels of CRP directly (23). Hence, conventional observational studies are prone to bias by potential confounding factors or reverse causation, and the causality was underpowered for definitive conclusions from previous studies.

Using genetic variations as instrumental variables (IVs), Mendelian randomization (MR) analysis is a novel epidemiology method to estimate the causation between an exposure and an outcome, with less impressionability to reverse causation and unmeasured confounders. In settings for which the IV assumptions are well justified, the findings could help optimize drug development or clinical trials and inform clinical or public health decision-making (24). Single nucleotide polymorphism (SNP) is a single base-pair difference in the DNA sequence of individuals within a species, which is the most common type of genetic variation in humans (25). Given that genetic variants are allocated at conception randomly, they are generally independent of environmental risk factors and precede the risk factors and the diseases’ onset (26). As for CRP, it has been estimated that the heritability of CRP blood concentrations was estimated from 25% to 40% (27), suggesting that genetic architecture might contribute to modulating CRP concentrations. Meanwhile, based on the summary data from genome-wide association studies (GWAS), our two-sample MR analysis could assess their causal relationships more comprehensively with robust statistical power (28,29). Previous MR studies (5,30-32) have illustrated that no causation was detected between genetically regulated elevated CRP concentrations and risk of LC. However, the sample size of LC patients in previous MR studies was relatively limited, insufficient to provide adequate statistical power to evaluate their causal nexus. Then, due to the limitation of the GWASs at the time, causal inferences provided by previous studies solely applied to populations of European ancestry; therefore, the findings were lack of generalizable value. Moreover, since genetically higher LC risk can influence CRP concentrations, no bidirectional MR study has been conducted, and therefore, previous findings were likely to be biased by reverse causality.

Generally, using the recent largescale meta-analysis of the GWASs updating CRP-specific SNPs among diverse populations, the present study could offer the latest comprehensive evidence for assessing the bidirectional causal effect between genetically regulated CRP concentrations and risk of LC with the two-sample MR method.

Methods

Study participants of lung cancer

In general, 11,348 LC cases and 15,861 controls from the International Lung Cancer Consortium (ILCCO) (33), an international research group that established and implemented LC research from different geographical areas and ethnicities, were used as epidemiological summary-level data. Subgroup analyses were stratified according to the pathology classification of non-small cell lung cancer (NSCLC). Specifically, we divided types of cancer into lung adenocarcinoma (LUAD) (3442 cases, 14,894 controls) as well as squamous cell lung cancer (LUSC) (3275 cases, 15,038 controls). The study was conducted following the Declaration of Helsinki (as revised in 2013) (34).

Genetic variants associated with elevated CRP concentrations

The latest large-scale GWAS datasets were retrieved from MR-Base, an open GWAS project developed at the MRC Integrative Epidemiology Unit (IEU) at the University of Bristol (35). Through Bonferroni correction to control the family-wise error rate (FWER), the P=5×10−8 threshold was widely accepted for association identification between a common genetic variant and a trait of interest, given the linkage disequilibrium (LD) structure of the genome (36). Using the MR-Base platform, SNPs associated with elevated CRP concentrations were initially selected from Neale Lab, Pan-UK Biobank (UKB) team, RIKEN Center for Integrative Medical Sciences, and the European Bioinformatics Institute (EBI) database of complete GWAS summary data at the genome-wide significance threshold (P<5×10-8). Utilizing LD analysis, we attempted to exclude SNPs once mutual LD surpassed the limited value (R2<0.001) with a larger P value conjugately. Eventually, the final IV set was established, including 204 SNPs (available online: https://cdn.amegroups.cn/static/public/tlcr-21-750-01.xlsx). Among them, 174 SNPs were of European ancestry, 13 SNPs were of Hispanic or Latin American origin, 7 SNPs were of South Asian ancestry, 6 SNPs were of East Asian ancestry, and 4 SNPs were of African American or Afro-Caribbean ancestry. Risk alleles and baseline alleles were encoded separately in accordance with the association with the rate of CRP concentrations. As the sample size of 204 instruments and 310,305 individual samples in our study, the F-statistic was 3,135.39 as estimated (37) given the type I error rate =0.05, which suggested noticeable correlativity for the MR analysis.

Statistical analysis

MR method was implied as the statistical analysis tool, which possesses three suppositions as foundations (38): (I) the IVs are robustly correlating with elevated CRP concentrations; (II) the IVs affect LC merely via their effect on elevated CRP concentrations (i.e., the IVs are independent of the outcome given the exposure), and (III) the IVs are independent of any confounders. As demonstrated by the previous studies (35) and the screening procedure of SNPs described in the previous paragraph, the first assumption was met.

Two-sample MR analyses were conducted to evaluate the potential bidirectional associations between elevated CRP concentrations and LC risk (). In addition to the aggregate effect, causal inferences were obtained from five independent ancestries for specific evaluation among different ethnic groups. Five well-established MR analytical methods were applied, including inverse-variance weighted (IVW), weighted median, MR-Egger, weighted mode, and simple mode. IVW was applied to the combination of multiple IVs as sole estimation of genetic variants by weighted score, which has the most statistical power among the five methods (39). The rest of the methods were performed for sensitivity analysis and to indirectly test the second assumption. Global pleiotropic effects were obtained from the MR-Egger analyses based on the intercept. Cochran's Q statistic and I2 were inspected to estimate heterogeneity. Furthermore, we conducted leave-one-out analyses to investigate whether the estimation of MR was determined or biased by an individual SNP by omitting a single SNP successively. MR analyses were performed in R (version 3.6.2) using the package TwoSampleMR (version 0.5.0) (35).

Figure 1

The overview of our study design aiming at revealing the bidirectional relations between C-reactive protein (CRP) and lung cancer. Two hundred and four CRP-specific single nucleotide polymorphisms (SNPs) (originated from East Asian, Hispanic/Latin American, European, African American/Afro-Caribbean, and South Asian ancestries), and 20 lung cancer-specific SNPs (including lung adenocarcinoma and squamous cell lung cancer) were applied as instrument variables respectively to evaluate the bidirectional causal relation.

Results

The basic characteristics and their F statistics for selected summary level GWASs applied in our MR study are given in , and the detailed information of each SNP used in our study is listed in table available online: https://cdn.amegroups.cn/static/public/tlcr-21-750-01.xlsx.

Table 1

The basic characteristics and their individual F statistics for selected summary level GWASs applied in our MR study.

Trait	Sample size	Sex	First author/consortium	F statistics
Overall C-reactive protein (CRP)	310,305	Males and females	–	3,135.39
East Asian-specific CRP variants	75,391	Males and females	Ishigaki K	762.53
Hispanic/Latin American-specific CRP variants	15,912	Males and females	Wojcik GL	161.73
European-specific CRP variants	204,402	Males and females	Neale lab	2,065.67
African American/Afro-Caribbean-specific CRP variants	6,203	Males and females	Pan-UKB team	63.66
South Asian-specific CRP variants	8,397	Males and females	Pan-UKB team	85.82
Lung cancer	27,209	Males and females	ILCCO	275.84
Lung adenocarcinoma	18,336	Males and females	ILCCO	186.21
Squamous cell lung cancer	18,313	Males and females	ILCCO	185.98

ILCCO, the International Lung Cancer Consortium.

In the MR analysis of mixed ancestries using 204 SNPs, the outcome of IVW method demonstrated that although genetically regulated one-unit increased in the log-transformed CRP concentrations were associated with a relatively 2.2% increased risk of LC, no statistically meaningful conclusions could be drawn (OR =1.022, 95% CI: 0.965–1.083, P=0.455) (), which was also consistent regarding pathological subtypes including LUAD (OR =1.026, 95% CI: 0.947–1.112, P=0.534) and LUSC (OR =1.060, 95% CI: 0.970–1.158, P=0.201). Findings were consistent among sensitivity analyses, demonstrating null causal relationship between both phenotypes (available online: https://cdn.amegroups.cn/static/public/tlcr-21-750-02.xlsx). Meanwhile, as for the causal effect among different ethnic populations, no association reached statistical significance among East Asian (OR =1.075, 95% CI: 0.754–1.532, P=0.690), Hispanic/Latin American (OR =0.998, 95% CI: 0.884–1.127, P=0.980), European (OR =1.044, 95% CI: 0.945–1.155, P=0.397), African American/Afro-Caribbean (OR =1.055, 95% CI: 0.884–1.260, P=0.551), and South Asian population (OR =0.999, 95% CI: 0.911–1.096, P=0.984). Detailed results on LC subtypes of ethnic-specific analyses are presented in .

Figure 2

Inverse-variance weighted results of MR study investigating the causal association between genetically regulated CRP levels on lung cancer among different ethnic populations. MR, Mendelian randomization; CRP, C-reactive protein.

Figure 3

Inverse-variance weighted results of bidirectional MR study investigating the causal association between lung cancer and CRP levels among different ethnic populations. MR, Mendelian randomization; CRP, C-reactive protein.

Regarding the impact of LC upon CRP concentrations, similar modest and nonsignificant estimates of causal effect were observed (OR =0.999, 95% CI: 0.977–1.021, P=0.923), conforming the absence of bidirectional effects. Subtypes including LUSC (OR =0.996, 95% CI: 0.970–1.023, P=0.790) and LUAD (OR =1.011, 95% CI: 0.991–1.031, P=0.300) supported the same conclusion. Lack of causal association also predominated among population subgroups ().

The results of sensitivity analyses among all study populations are shown in table available online: https://cdn.amegroups.cn/static/public/tlcr-21-750-02.xlsx. Using MR-Egger regression, our findings did not support the existence of global pleiotropic assumptions in most of the study outcomes (Table S1), while evidence of heterogeneity was found in LC effect on European-specific CRP variants (Table S2). Leave-one-out studies of the LC overall and subgroup analyses showed no evidence that a single SNP had an impact upon the overall effect of CRP variants on LC risks (available online: https://cdn.amegroups.cn/static/public/tlcr-21-750-03.xlsx).

Discussion

This is the first study to conduct bidirectional MR analyses between CRP concentrations and LC among different ethnic backgrounds. No bidirectional causation between genetically regulated elevated concentrations of CRP and LC risk was observed in our study among East Asian (P=0.690, nSNP =6, n=75,391), Hispanic/Latin American (P=0.980, nSNP =13, n=15,912), European (P=0.397, nSNP =174, n=204,402), African American/Afro-Caribbean (P=0.551, nSNP =4, n=6,203), and South Asian populations (P=0.984, nSNP =7, n=8,397). Additionally, no causality was found in subgroup analyses concerning pathologic types as well.

The results are inconsistent with most previous prospective cohort studies (3,11,12,14,16-19) using measured plasma concentrations of circulating CRP to explore the risk of LC. Shiels et al. (16) demonstrated that elevated measured CRP concentrations were related to the risk of LC in two independent studies from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Allin et al. (11) also illustrated that elevated plasma concentrations of CRP in cancer-free individuals were relevant to an increased risk of LC. Besides, the latest meta-analysis (21) denoted that elevated plasma concentrations of CRP were also related to an increased risk of LC. Interestingly, after a stratified analysis (15) based on smoking status, the positive association between circulating high sensitivity CRP concentrations and the risk of LC disappeared when considering non-smokers. Thus, rather than a causal factor, the authors (15) concluded that circulating high sensitivity CRP concentration could be a pre-diagnostic marker of LC. Since traditional observational studies are vulnerable to potential confounding factors or reverse causation, we supposed that previous results might be influenced by underlying confounders (e.g., cigarette smoking status and older age stratifications).

Our study showed similar results in effect directions compared with previous MR studies (5,30-32), whereas these studies came to the conclusions from a limited LC sample size. Among them, the largest sample size of LC patients was only 416 (30). Therefore, they might not be precise enough to evaluate causation between genetically elevated CRP concentrations and LC risk due to the deficiency of statistical power. In comparison, we specifically selected a total of 27,209 participants from ILCCO with 11,348 LC cases and 15,861 controls as our study participants, which were much more extensive than previous ones and offered adequate statistical power. Meanwhile, two (5,30) of them didn’t identify SNPs through GWAS, and the majority of chosen SNPs haven’t been confirmed by subsequent GWAS (27,40,41), indicating those SNPs might not be authentically associated with elevated CRP concentrations. Even for the other two studies (31,32) with GWAS-identified SNPs, some of their chosen SNPs were also proposed to have nothing to do with elevated CRP concentrations subsequently (27,40,41), such as rs4903031, rs6901250, and rs4705952. Additionally, genetic pleiotropy could occur without sensitivity analyses, indicating previously chosen SNPs might relate to other inflammatory processes and result in a spurious association. Consequently, biased estimation for the causation could happen due to those irrelevant SNPs and potential global pleiotropic effects in previous MR studies. Furthermore, more relevant SNPs have been updated in recent years. Our study included 204 SNPs from the latest GWAS as the IV set, which would explain the variance of elevated CRP concentrations more comprehensively and precisely. Generally, our updated MR study provided the latest and much more substantial evidence for evaluating the causality between genetically regulated elevated CRP concentrations and LC risk.

Regarding the etiology of LC, increasing evidence have illustrated that chronic inflammation may play a significant role amid the carcinogenesis process (4-7,16,42,43). Inflammatory conditions, such as chronic obstructive pulmonary disease (3,44-47), chlamydia pneumonia (4,48), tuberculosis (6), and human immunodeficiency virus (HIV) (49) seem to correlate with an increased LC risk. Through the nuclear factor kappa B (NF-κB) pathway, chronic inflammation may actively promote the malignant process, the increase in angiogenesis and proliferation of cells during tissue repair, and the up-regulation of antiapoptotic genes through reactive oxygen and nitrogen species production (9). Considered a systemic marker of chronic inflammation (5), elevated plasma CRP concentrations were associated with an increased LC risk in many conventional observational studies. However, our MR analysis denoted that elevated CRP concentrations were not causally associated with LC risk, which could be explained from two aspects. For one thing, chronic inflammation promotes the carcinogenesis of the lung, possibly by means of other inflammatory mediators except for CRP. It has been reported that soluble tumor necrosis factor receptor-2 (sTNFRII), serum amyloid A (SAA), and monokine induced by gamma interferon (CXCL9/MIG) were also related to LC risk (16). For another, due to the natural defects of observational studies, their results could be influenced by unmeasured confounders, such as smoking status and elder age. In addition, inflammation and elevated CRP concentrations could be induced by tumor growth, and CRP could also be an indicator of an immune response to tumor antigens (10). Consequently, elevated CRP concentrations might also respond to occult or very early-stage cancer (16). However, CRP may not always act as a cancer promoter but sometimes play a role in treating cancer. Sasaki et al. observed that CRP can inhibit lymph node metastasis and lymphatic angiogenesis of squamous cell carcinoma through injecting CRP into mice model (50). As for the prognostic effects of CRP, Shinohara and colleagues found high serum circulating CRP levels on postoperative days were associated with enhanced 5-year overall survival (OS) as well as recurrence-free survival (RFS) in NSCLC patients (51). Instead, Okada et al. reported the high serum CRP levels during perioperative period were a poor prognostic factor for OS and RFS in NSCLC patients (52). Consequently, the impact of CRP levels on the incidence and prognosis of LC was nuanced and complex, and further studies with larger sample sizes and better designs are in an urgent need.

Several merits deserve explicit mention. First, this is the first study to assess the bidirectional causation between elevated CRP concentrations and LC risk to date. Second, we used the most comprehensive GWAS-identified SNPs as an IV set. Compared with the previous studies, our study with 204 SNPs explained more variance of elevated CRP concentrations. Thanks to the robustly associated IVs (F-statistics =3,135.39) and the great sample size (n=310,305), our study was capable of providing a relatively accurate evaluation of a causal inference. Third, to the best of our knowledge, this is the first MR study that conducted stratification analyses according to different ethnic populations. Thus, we were able to investigate whether their association could vary under different ancestry backgrounds.

Notwithstanding the advantages, several limitations in our study cannot be neglected. First, even though the most comprehensive set of multiethnic genetic variants were applied for our research so far, only a part of the variance of elevated CRP concentrations could be explained in the population. Some unknown CRP-related SNPs, which deserve further investigation in future GWASs, could play an essential role in LC development. Besides, due to the methodology limitations currently, MR assumptions could not be thoroughly tested, and therefore potential violations against the assumptions may occur. To overcome this difficulty, instead of directivity evaluated in the second assumption, we implemented additional sensitivity analyses, which showed the absence of pleiotropic effects in our study, suggesting no violation of the second MR assumption. Moreover, depending on GWAS summary statistics, the two-sample MR methods assume a linear relationship between the exposure and the outcome. Hence, it is possible that a few genetic variants were also related to confounding factors of CRP and LC.

Overall, our study negates bidirectional causal effects of CRP concentrations on LC among East Asian, Hispanic/Latin American, European, African American/Afro-Caribbean, and South Asian populations. Still, it opens a new concept for the current research orientation. Further investigation of the profound relationship between both phenotypes is required to be unveiled in pathologic and biochemistry aspects.

The article’s supplementary files as