Jiejun Zhou1†
Youqian Zhang
Kun Zhang
Meng Li
Mingwei Chen- 1Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
- 2Health Science Center, Yangtze University, Jingzhou, Hubei, China
Background: Former research has emphasized a correlation between lung cancer (LC) and sepsis, but the causative link remains unclear.
Method: This study used univariate Mendelian Randomization (MR) to explore the causal relationship between LC, its subtypes, and sepsis. Linkage Disequilibrium Score (LDSC) regression was used to calculate genetic correlations. Multivariate MR was applied to investigate the role of seven confounding factors. The primary method utilized was inverse-variance-weighted (IVW), supplemented by sensitivity analyses to assess directionality, heterogeneity, and result robustness.
Results: LDSC analysis revealed a significant genetic correlation between LC and sepsis (genetic correlation = 0.325, p = 0.014). Following false discovery rate (FDR) correction, strong evidence suggested that genetically predicted LC (OR = 1.172, 95% CI 1.083–1.269, p = 8.29 × 10−5, Pfdr = 2.49 × 10−4), squamous cell lung carcinoma (OR = 1.098, 95% CI 1.021–1.181, p = 0.012, Pfdr = 0.012), and lung adenocarcinoma (OR = 1.098, 95% CI 1.024–1.178, p = 0.009, Pfdr = 0.012) are linked to an increased incidence of sepsis. Suggestive evidence was also found for small cell lung carcinoma (Wald ratio: OR = 1.156, 95% CI 1.047–1.277, p = 0.004) in relation to sepsis. The multivariate MR suggested that the partial impact of all LC subtypes on sepsis might be mediated through body mass index. Reverse analysis did not find a causal relationship (p > 0.05 and Pfdr > 0.05).
Conclusion: The study suggests a causative link between LC and increased sepsis risk, underscoring the need for integrated sepsis management in LC patients.
1 Introduction
Sepsis is described as “a life-threatening organ dysfunction caused by dysregulated host systemic inflammatory and immune response to infection”. (Singer et al., 2016).Sepsis remains a serious global health challenge (Tiru et al., 2015; Zhu et al., 2022). Sepsis impacted close to 50 million individuals globally and accounted for about 20% of global deaths before the COVID-19 pandemic (Rudd et al., 2020). Key factors contributing to the development of sepsis include the pathogen’s virulence, the site and type of infection, and host factors such as age, genetic predisposition, and comorbidities (Mao et al., 2013; Pan et al., 2017; Moon et al., 2023). Despite advancements in understanding its mechanisms, sepsis remains a significant challenge in healthcare due to its rapid progression and high mortality rate (Bauer et al., 2020). Early recognition and prompt management are crucial in improving patient outcomes.
Lung cancer (LC) is responsible for 2.2 million new cases each year, ranking as the world’s second most prevalent cancer. Furthermore, it is the primary cause of cancer-related mortality, resulting in approximately 1.79 million deaths annually (Sung et al., 2021; Thai et al., 2021; Huang et al., 2022). From a pathological classification perspective, LC can be roughly divided into two types: non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC). NSCLC primarily includes histological subtypes such as adenocarcinoma (LUAD) and squamous cell carcinoma (SqCLC). The relationship between LC and sepsis is complex and intricate (Pavon et al., 2013; Hensley et al., 2019; Mirouse et al., 2020; Moore et al., 2020; Shvetsov et al., 2021; Xia et al., 2022). Most research evidence supports the idea that LC poses a risk for sepsis (Pavon et al., 2013; Hensley et al., 2019; Moore et al., 2020; Xia et al., 2022). The immunosuppressive state induced by the tumor itself or by cancer treatments can make these patients more susceptible to infections (Pavon et al., 2013; Hensley et al., 2019; Xia et al., 2022). The results of a prospective study on the likelihood of sepsis following cancer showed a heightened risk of sepsis in cancer survivors, these cancers include lung cancer, breast cancer, prostate cancer, and other solid tumors, as well as hematological tumors (Moore et al., 2020). Additionally, studies have shown that the incidence and mortality rates of sepsis among cancer patients are higher, severe sepsis is associated with 8.5% of all cancer deaths, costing 3.4 billion dollars per year (Williams et al., 2004). However, some study results differ from this view; for example, a study suggests both complementary and antagonistic relationships between cancer and sepsis (Mirouse et al., 2020). Yurii B. Shvetsov and colleagues, in a multiethnic cohort study concerning the association between sepsis mortality and specific cancer sites and treatment types, found that lung cancer was associated with a significantly lower increase in sepsis mortality compared to non-sepsis mortality (Shvetsov et al., 2021).
Conversely, there is currently no consensus on whether sepsis increases the risk of cancer incidence. A multicenter observational study suggests that the incidence of sepsis does not alter the oncological and prognostic results in patients with epithelial ovarian cancer (Said et al., 2023). However, another study confirms that sepsis was significantly linked to a higher risk of nine types of cancer within 5 years after sepsis diagnosis, including LC (Liu Z. et al., 2019).
Given these inconsistent academic findings and the constraints of observational research in establishing cause-and-effect relationships, Mendelian randomization (MR) can provide insights into causality that observational studies lack. MR uses genetic variations as instrumental variables (IVs) derived from genomic-wide linkage analyses for causality inference. MR functions like a natural randomized controlled trial (RCT), offering more substantial evidence and less vulnerability to confounding factors than observational studies. MR is extensively used in cancer and disease research (Li and Wang, 2023; Xu et al., 2023). Hence, performing a bidirectional MR study could be critical in deciphering causal links between sepsis and LC, paving the way for better prevention and therapies.
2 Materials and methods
2.1 Study design
This study explores the causal relationship between LC and sepsis using summary-level data from the largest publicly accessible genome-wide association study (GWAS) currently available on these conditions. A suite of sophisticated analyses was conducted, incorporating bidirectional univariate MR, complementary multivariable MR (MVMR) analysis, and in-depth genetic correlation evaluations. IVs for exposure were established based on stringent criteria: (i) strong association of the genetic instrument with the exposure; (ii) independence of the instrument from confounding variables; (iii) the exclusive pathway of the genetic variants’ impact on the outcome is through the exposure (Lawlor et al., 2008). The methodological intricacies of the MR framework are presented in Figure 1, while the comprehensive summary data are systematically detailed in Table 1. This study is reported following the Strengthening the Reporting of Observational Studies in Epidemiology Using Mendelian Randomization guidelines (STROBE-MR) (Skrivankova et al., 2021).
Figure 1. Overview of research design and analysis strategy. Overview of the research design. The MR framework is based on three fundamental MR assumptions. IVs, instrumental variables; MR, Mendelian randomization; MR-PRESSO, MR Pleiotropy Residual Sum and Outlier; SNP, single nucleotide polymorphism.
Table 1. Detailed information of data sources.
2.2 Selection of genetic instrumental variables
The MR analysis operationalized rigorous selection parameters for Single nucleotide polymorphism (SNP) identification: (i) SNP, to serve as instrumental variables, showcased genome-wide significant associations with the exposure (p < 5 × 10−8). In the reverse analysis, due to the inability to obtain SNPs at the genome-wide significance level for the sepsis phenotype, we adjusted to a more relaxed threshold (p < 5 × 10−6) based on previous MR analysis experience to acquire a sufficient number of SNPs for the analysis (Liu et al., 2024; Xu et al., 2024; Yang et al., 2024). (ii) The selection of SNPs underwent rigorous scrutiny to exclude confounding variable associations and to confirm independence, thus preventing biases due to linkage disequilibrium (r2 < 0.001, clumping distance = 10,000 kb). (iii) SNP validity as instrumental variables was gauged by F-statistics (F = R2 ×(N− 2)/(1 − R2)), where R2 denotes the percentage of variance in the exposure explained by the SNPs, and N is the sample size of the GWAS from which the exposure is drawn). This criterion helped to eliminate weak instruments, with an F-statistic threshold of >10, ensuring robust instrument strength (Teumer, 2018). (iv) MR-Steiger filtering was applied to remove variants demonstrating stronger associations with outcomes than with exposures (Hemani et al., 2017). (v) In instances of SNP unavailability in the outcome dataset, the SNiPa web interface (http://snipa.helmholtz-muenchen.de/snipa3/) was employed, leveraging genotype data from the European cohort of the 1000 Genomes Project Phase 3, to locate a proxy SNP in strong linkage disequilibrium with the primary SNP (r2 > 0.8). (vi) Consistency was essential, with the SNP’s effects on exposure and outcome required to be in the same allelic direction.
2.3 Source of lung cancer phenotype
This study utilizes the most extensive dataset to date, drawn from a meta-analysis by McKay JD et al., encompassing European ancestry GWAS for LC with 29,266 cases and 56,450 controls, SqCLC with 7,426 cases and 55,627 controls, LUAD with 11,273 cases and 55,483 controls, and SCLC with 2,664 cases and 21,444 controls (McKay et al., 2017). The research integrates novel data from the OncoArray genotyping platform with existing data from previous LC GWAS, conducting a large-scale association analysis on over 29,000 patients and 56,000 controls of European descent.
2.4 Source of sepsis phenotype
The latest and most exhaustive aggregate GWAS analysis about sepsis is derived from the UK Biobank (Sudlow et al., 2015). Methodological adjustments in this study included age, sex, ten principal genetic components, and genotyping batch effects. It comprised 11,643 sepsis cases juxtaposed against 474,841 controls, all of European descent. Case identification hinged on the presence of ICD-10 codes A02, A39, A40, and A41.
2.5 MVMR models
Lung cancer and sepsis share common risk factors, including smoking, chronic obstructive pulmonary disease (COPD), body mass index (BMI), chronic conditions like type 2 diabetes mellitus (T2DM), and other factors such as education level (Biesalski et al., 1998; Pallis and Syrigos, 2013; Henriksen et al., 2015; Bohl et al., 2016; Malhotra et al., 2016; Locham et al., 2021; Bladon et al., 2024). Considering potential confounding factors, seven major confounding factors were selected. BMI data was sourced from the Genetic Investigation of Anthropometric Traits (GIANT) consortium (Yengo et al., 2018). Educational attainment (EA) data came from the Social Science Genetic Association Consortium (SSGAC) (Lee et al., 2018). T2DM data was acquired from the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) consortium (Mahajan et al., 2022). Data on cigarettes per day (CPD), age of smoking initiation (ASI), and alcoholic drinks per week (ADPW) were derived from the GWAS and Sequencing Consortium of Alcohol and Nicotine Use (GSCAN) (Liu M. et al., 2019). COPD data came from the Global Biobank Meta-analysis Initiative (GBMI) (Zhou et al., 2022).
2.6 Statistical analyses
2.6.1 MR analysis
The univariate MR framework evaluated individual IVs using the Wald ratio, which calculates the causal effect by dividing SNP-outcome association (β_Y) by the SNP-exposure association (β_X). This method provides a causal estimate for each genetic variant, assuming IVs significantly influence the exposure, are independent of confounders, and affect the outcome solely through the exposure (Burgess et al., 2013). Concurrently, to elucidate the causal associations involving multiple IVs (two or more), use the multiplicative random-effects inverse-variance-weighted (IVW) method (Burgess et al., 2013). It is critical to note that when the heterogeneity index I2 is below 50%, outcomes derived from the fixed-effects model are considered robust. This statistical strategy was additionally refined by incorporating the MR-Egger and weighted median methodologies. The IVW technique’s weighting schema is coherent with the Wald ratio estimates for each SNP, inversely correlating with its variance (Hemani et al., 2018). By integrating the full complement of genetic variants, the IVW approach ensures systematic and reliable results. Contrariwise, the weighted median method gains prominence when more than half of the genetic variants are presumed invalid; concurrently, the MR-Egger method operates under the premise that all such variants are invalid (Bowden et al., 2016). Furthermore, the constrained maximum likelihood (CML) process was employed, allowing for collective analysis over an expansive array of genetic variants while adjusting for possible confounders and intrinsic genetic heterogeneity. Particularly when addressing a comprehensive array of genetic variants and confounders, the CML method is indispensable for obtaining accurate and robust results (Zhang et al., 2008).
Additional MVMR analyses were carried out to delineate the direct causal pathways from exposure to outcome (Burgess and Thompson, 2015). These analyses were to define the direct causal connections precisely, thereby differentiating them from the univariate MR model. Contrary to UVMR, which concentrates on a singular exposure, MVMR considers genetic variations linked to multiple exposures. The initial stage involved generating Mendelian Randomization effect estimates for the exposure-to-outcome relationships using the IVW method. Subsequently, an MVMR assessment was carried out to assess the impact of six mediators on the outcome, considering the specific attributes of the exposure.
This study conducted its analysis using the “MendelianRandomization,” “TwoSampleMR,” “MR-PRESSO,” and “MRMR” packages in R version 4.3.0 software.
2.6.2 LDSC regression analysis
The linkage disequilibrium score (LDSC) regression, designed for analyzing GWAS summary data, constitutes an effective instrument for discerning genetic correlations among complex diseases or traits. This method facilitates the separation of authentic polygenic influences from confounding factors, which include subtle family structures and population stratification (Bulik-Sullivan et al., 2015). A significant genetic correlation, characterized by statistical solidity and substantial effect size, implies that the correlation between phenotypes is not merely attributable to environmental influences. The LDSC, available at (https://github.com/bulik/ldsc), provides a direct path for investigating the genetic foundations linking exposure and outcome traits.
2.6.3 Sensitivity analysis
Diversity among chosen genetic variants was measured using Cochran’s Q test, with a p-value less than 0.05, denoting notable differences within the SNPs under study (Kulinskaya et al., 2020). MR-Egger regression was utilized to investigate directional pleiotropy within the MR context (Burgess and Thompson, 2017). An MR-Egger intercept with a p-value below 0.05 indicates notable directional pleiotropy despite the acknowledged limitations of this method (Wu et al., 2020). The MR Pleiotropy Residual Sum and Outlier (MR-PRESSO) approach was used to pinpoint outliers and evaluate horizontal pleiotropy, with a global p-value below 0.05 confirming its presence (Verbanck et al., 2018). Outliers were rigorously removed to refine analytical accuracy. This was followed by a leave-one-out sensitivity analysis to appraise the effect of single SNPs on the collective outcomes (Cheng et al., 2017). The false discovery rate (FDR) method was employed to correct multiple comparisons rigorously. Post-correction, p-values less than 0.05 denoted significant causal associations. Conversely, outcomes with raw p-values under 0.05 that did not maintain this significance after FDR adjustment were classified as suggestive rather than definitive.
For the calculation of R2 values, equation 2×MAF×(1-MAF)×beta2 was used, where MAF represents the minor allele frequency for each SNP. These calculated R2 values were then gathered to establish the combined parameter for power calculation (Guan et al., 2014). The mRnd platform (Brion et al., 2013) (https://shiny.cnsgenomics.com/mRnd/) provided the means for statistical power assessment.
3 Results
3.1 Genetic instrument selection and genetic correlation between phenotypes
The research findings suggest that the F-statistics for all instrumental variables surpassed 500, suggesting a significant decrease in bias due to weak instrument variation. The quantity of SNPs selected as instrumental variables ranged from one to 14, with the explained variance in genetic variation ranging from 1.75% to 12.36% (Supplementary Table S1). The scatter plot (Figure 2) provides an intuitive representation of the direction of causal associations, while the forest plot (Supplementary Figure S1) displays the effects contributed by all IVs. Detailed SNP information can be found in Supplementary Tables S2–S8.
Figure 2. Scatterplot summary of all MR analyses. The vertical and horizontal lines denote the 95% confidence intervals for the effect size, while the slopes of the fitted lines indicate the estimated Mendelian randomization effect per method. (A) Lung adenocarcinoma on sepsis (B) Lung cancer on sepsis (C) Squamous cell lung carcinoma on sepsis (D) Sepsis on lung adenocarcinoma (E) Sepsis on small cell lung carcinoma (F) Sepsis on lung cancer (G) Sepsis on squamous cell lung carcinoma. MR, Mendelian randomization; SNP, single nucleotide polymorphism.
The LDSC study uncovered a substantial genetic link between LC and sepsis (rg = 0.325, p = 0.014); however, no genetic correlation was detected between the subtypes of LC: SqCLC (rg = 0.095, p = 0.605), LUAD (rg = 0.295, p = 0.052), and SCLC (rg = 0.151, p = 0.305), and sepsis. The SNP-based liability-scale heritability (h2) ranged from 0.24% to 10.13% (Supplementary Table S3).
3.2 Association of genetically predicted lung cancer with sepsis
In the progressive MR analysis (Figure 3), a significant causal relationship was identified between all subtypes of LC and sepsis. Expressly, after FDR correction for multiple comparisons, the primary analysis method, IVW, indicated that each standard deviation (SD) increase in genetically predicted LC was linked to a 17.2% increase in the risk of sepsis, with both random (OR = 1.172, 95% CI 1.083–1.269, p = 8.29 × 10−5, Pfdr = 2.49 × 10−4) and fixed (OR = 1.172, 95% CI 1.098–1.252, p = 2.24 × 10−6, Pfdr = 6.72 × 10−6, I2 = 30%) effects models showing concordance. Supplementary methods provided consistent evidence of this causal relationship. With an OR of 1.172, we achieved 100% statistical power to detect the association between LC and sepsis. Among LC subtypes, for each SD increase in genetically predicted SqCLC, the risk of sepsis increased by 9.8%, with IVW random (OR = 1.098, 95% CI 1.021–1.181, p = 0.012, Pfdr = 0.012) and fixed (OR = 1.098, 95% CI 1.021–1.181, p = 0.012, Pfdr = 0.012, I2 = 0%) effects models in agreement. Supplementary methods, including the Weighted median (OR = 1.148, 95% CI 1.043–1.264, p = 0.005) and CML (OR = 1.099, 95% CI 1.022–1.181, p = 0.011), also yielded consistent causal evidence. At an OR of 1.098, we possessed 78% statistical strength to identify the relationship between SqCLC and sepsis. Similarly, for each SD increase in genetically predicted LUAD, there was a 9.8% increase in sepsis risk, with concordance across random (OR = 1.098, 95% CI 1.024–1.178, p = 0.009, Pfdr = 0.012) and fixed (OR = 1.098, 95% CI 1.040–1.159, p = 7.27 × 10−4, Pfdr = 0.001, I2 = 39%) effects models in IVW and consistency in supplementary methods excluding the simple mode (OR = 1.150, 95% CI 0.998–1.326, p = 0.076). With an OR of 1.098, we had a 96% statistical strength to determine the relationship between LUAD and sepsis. Finally, in the case of SCLC and sepsis, due to the limited IVs available for analysis, only Wald ratio (OR = 1.156, 95% CI 1.047–1.277, p = 0.004) and CML (OR = 1.159, 95% CI 1.048–1.283, p = 0.004) were used, which also showed a 100% statistical strength to identify the association. However, due to the inability to perform multidimensional validation, the evidence indicated a suggestive risk.
Figure 3. Summary of results from MR analysis of genetically predicted lung cancer phenotypes for sepsis. CML, constrained maximum likelihood; FDR, false discovery rate; IVW, inverse-variance-weighted; LC, lung cancer; LUAD, lung adenocarcinoma; MR, Mendelian randomization; MR-PRESSO, MR Pleiotropy Residual Sum and Outlier; OR, odd ratio; P-val, p-value; SCLC, small cell lung carcinoma; SqCLC, squamous cell lung carcinoma; SNP, single nucleotide polymorphism.
In the sensitivity analysis (Table 2), all IVs passed the MR-Steiger filtering, and Cochran’s Q statistic indicated no significant heterogeneity (p > 0.05). Likewise, MR-Egger and MR-PRESSO tests revealed no pleiotropy (p > 0.05). Leave-one-out analysis confirmed that the causal inference was not driven by any single SNP (Supplementary Figure S2), and the funnel plot exhibited a symmetrical distribution (Supplementary Figure S3).
Table 2. Summary of sensitivity results.
3.3 Association of genetically predicted sepsis with lung cancer
In the reverse MR analysis (Figure 4), the primary method, IVW, indicated no reverse causal effects between sepsis and LC (OR = 0.983, 95% CI 0.878–1.102, p = 0.771, Pfdr = 0.771), SqCLC (OR = 1.087, 95% CI 0.875–1.350, p = 0.451, Pfdr = 0.771), LUAD (OR = 1.039, 95% CI 0.887–1.216, p = 0.638, Pfdr = 0.771), or SCLC (OR = 0.938, 95% CI 0.684–1.286, p = 0.691, Pfdr = 0.771). Additional supplementary methods provided consistent evidence supporting these associations (p > 0.05). Comprehensive sensitivity analyses confirmed the robustness of the findings (p > 0.05), as detailed in Table 2. Corresponding figures are provided in Supplementary Figures S2–S3.
Figure 4. Summary of results from MR analysis of genetically predicted sepsis for lung cancer phenotypes. CML, constrained maximum likelihood; FDR, false discovery rate; IVW, inverse-variance-weighted; LC, lung cancer; LUAD, lung adenocarcinoma; MR, Mendelian randomization; MR-PRESSO, MR Pleiotropy Residual Sum and Outlier; OR, odds ratio; P-val, p-value; SCLC, small cell lung carcinoma; SqCLC, squamous cell lung carcinoma; SNP, single nucleotide polymorphism.
3.4 MVMR analysis
In the univariate MR analysis, evidence supported a causal relationship between LC and its subtypes and the risk of sepsis, reaching statistical significance (p < 0.05 and Pfdr < 0.05). Since only one IV related to SCLC was obtained, we excluded the SCLC phenotype from further analysis to avoid biased results in MVMR adjustments. In the MVMR analysis (Table 3), the causal link between LC and sepsis was no longer notable when adjusting for BMI (OR = 1.023, 95% CI 0.944–1.108, p = 0.581). Similarly, SqCLC adjusted for BMI (OR = 1.031, 95% CI 0.975–1.091, p = 0.279), EA (OR = 1.052, 95% CI 0.990–1.118, p = 0.099), T2DM (OR = 1.070, 95% CI 0.992–1.154, p = 0.081), and LUAD adjusted for BMI (OR = 0.989, 95% CI 0.923–1.060, p = 0.752) and CPD (OR = 1.060, 95% CI 0.989–1.136, p = 0.101) also showed no significant causal relationship with sepsis. This suggests that these confounding factors may partially mediate the causal relationship between LC, subtypes, and sepsis.
Table 3. Summary of analytical results for MVMR.
4 Discussion
This research undertook a comprehensive MR analysis to explore the link between genetic predisposition to sepsis and LC. The results of the MR support earlier epidemiological research (Pavon et al., 2013; Rhee et al., 2017; Xia et al., 2022), confirming a causal link between sepsis and LC. Furthermore, no reverse causal link was found between LC and sepsis. Additional MVMR analysis suggested that factors such as body mass index, level of education, type 2 diabetes, and information on daily cigarette consumption might play a role in mediating part of this causative link.
Prior studies have indicated a link between sepsis and LC, with results showing a correlation between LC and an increased likelihood of sepsis (Pavon et al., 2013; Hensley et al., 2019; Xia et al., 2022). The immunosuppressive state induced by the tumor itself or by cancer treatments can make these patients more susceptible to infections (Pavon et al., 2013). In a group of more than one million hospital admissions for sepsis in the U.S., over 20% had a connection to cancer (Hensley et al., 2019). The results of a prospective study on the risk of sepsis after cancer showed an increased risk of sepsis in cancer survivors (Moore et al., 2020). Nonetheless, previous research has yielded inconsistent results concerning the relationship between general sepsis and LC. Notably, one study demonstrated that the relationship between cancer and sepsis is complementary and antagonistic (Mirouse et al., 2020). Our results are similar to the results of previous mainstream studies. Our MR results support the idea that LC contributes to sepsis.
Alternatively, there is currently no consensus on whether sepsis increases the risk of cancer incidence. A multicenter observational study suggests that the occurrence of sepsis does not affect the oncological and survival outcomes in patients with epithelial ovarian cancer (Said et al., 2023). However, another study confirms that sepsis was notably linked with a heightened risk of nine different cancer types in 5 years after a sepsis diagnosis, including LC (20). In our study, sepsis was not found to act as a genetic predisposing element for LC.
Based on previous disparities in research and incorporating the findings of this study, we posit that, given cancer and sepsis are not singularly unique diseases, it is evident that the risk associated with breast cancer differs from that of pancreatic cancer. Biological distinctions also exist between solid malignancies and malignancies of the hematopoietic system. Consequently, categorizing all cancer or sepsis patients uniformly is erroneous, and treatments conducted in this context are likely to yield suboptimal outcomes. Sepsis, akin to cancer, exhibits intricate diversity. Gaining a deeper understanding of the distinctive physiological states induced by sepsis and cancer is complex and crucial.
Observational studies often face limitations due to unobserved confounding factors and reverse causality, focusing more on correlation than causation. While data indicates a connection, the causal link between sepsis and LC has not yet been conclusively proven. We used MR analysis to investigate the genetic underpinnings of the causative link between sepsis and LC to counteract biases and confounders. Our study indicates that we must pay attention to monitoring patients with LC infection and inflammation factors and prevention and intervention in promptly treating sepsis.
Further MVMR analyses underscored the significance of BMI, educational attainment, T2DM, and data on cigarettes per day. Firstly, high BMI, especially obesity, may increase the risk of infection because obesity may affect the function of the immune system and may be associated with chronic inflammation. An MR study also showed that obesity was linked to a heightened likelihood of developing sepsis (Hu et al., 2023). Secondly, educational attainment may indirectly affect an individual’s risk of sepsis by influencing their lifestyle, health behaviors, access to medical resources, and the environment in which they live and work. For example, previous observational research has shown that lower educational attainment (EA) levels are connected to a heightened likelihood of COVID-19 (Jian et al., 2021). Additionally, individuals with diabetes have a higher propensity to develop wounds and ulcers that do not heal and can become infected, resulting in sepsis. In addition, diabetes alters the immune system, leading to an increased risk of sepsis (Schuetz et al., 2011). Lastly, smoking may elevate infection risk by increasing proinflammatory cytokines, damaging endothelial cells, and correlates with poor health habits (Alroumi et al., 2018; Zhang et al., 2022). A causative link between smoking and infectious disease risk was also shown in an MR study (Zhu et al., 2023). Therefore, this implies a comprehensive strategy for managing sepsis, considering these factors combined.
Our research has several advantages. This MR study represents the inaugural exploration of the causative link between sepsis and LC at the genetic level. All the SNPs selected as instrumental variables (IVs) originated from the European demographic, thus diminishing the probability of population stratification bias and bolstering the credibility of the bidirectional MR hypothesis. Our robust tools in this research (such as an F statistic significantly exceeding 10) should mitigate potential bias from sample overlap. Nonetheless, our investigation has its limitations. Several initial exposures were sourced from the UKB cohort, and the absence of additional GWAS hindered the execution of a confirmatory control analysis. Additionally, the exclusive access to summary-level GWAS data impeded the conduct of more detailed subgroup analyses.
In conclusion, our study used a comprehensive approach to investigate the association between lung cancer and sepsis, providing novel insights. Our results suggest that LC is a significant risk factor for sepsis, however, sepsis was not found to act as a genetic predisposing factor for LC. Further in-depth research is warranted to unravel the additional intricacies of this relationship. These efforts underscore the need for integrated sepsis management in LC patients.
5 Conclusion
To summarize, our study establishes a causal relationship between LC and increased risk of sepsis, with no evidence for a reverse association. Comprehensive prevention and treatment of sepsis should be carried out in LC patients, especially those with high BMI, low educational attainment, T2DM, and smoking.
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 author.
Author contributions
JZ: Conceptualization, Data curation, Investigation, Methodology, Writing–original draft, Writing–review and editing. YZ: Conceptualization, Data curation, Investigation, Methodology, Writing–original draft, Writing–review and editing. TY: Funding acquisition, Validation, Visualization, Writing–review and editing. KZ: Methodology, Validation, Visualization, Writing–review and editing. AL: Validation, Visualization, Writing–review and editing. ML: Validation, Visualization, Writing–review and editing. XP: Methodology, Visualization, Writing–review and editing. MC: Conceptualization, Funding acquisition, Investigation, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Key R&D Program of Shaanxi Province (grant number 2021ZDLSF02-05), the Science and Technology Project of Shaanxi Province (grant number 2022SF-580), and Wu Jie-ping Medical Foundation (grant number HX202293).
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/fgene.2024.1381303/full#supplementary-material
Abbreviations
ADPW, alcoholic drinks per week; ASI, age of smoking initiation; BMI, body mass index; CI, confidence interval; CML, constrained maximum likelihood; COPD, chronic obstructive pulmonary disease; CPD, cigarettes per day; df, degree of freedom; DIAGRAM, DIAbetes Genetics Replication And Meta-analysis; EA, education attainment; GBMI, Global Biobank Meta-analysis Initiative; GIANT, Genetic Investigation of Anthropometric Traits; GSCAN, GWAS and Sequencing Consortium of Alcohol and Nicotine use; IVs, instrumental variables; LC, lung cancer; LUAD, lung adenocarcinoma; MR, Mendelian Randomization; MR-PRESSO, MR Pleiotropy Residual Sum and Outlier; MVMR, Multivariable Mendelian randomization; OR, odds ratio; P-val, p-value; RSSobs, residual sum of squares observation; SCLC, small cell lung carcinoma; SNP, single nucleotide polymorphism; SqCLC, squamous cell lung carcinoma; SSGAC, Social Science Genetic Association Consortium; T2DM, type 2 diabetes mellitus.
References
Alroumi, F., Abdul Azim, A., Kergo, R., Lei, Y., and Dargin, J. (2018). The impact of smoking on patient outcomes in severe sepsis and septic shock. J. Intensive Care 6, 42. doi:10.1186/s40560-018-0312-x
Bauer, M., Gerlach, H., Vogelmann, T., Preissing, F., Stiefel, J., and Adam, D. (2020). Mortality in sepsis and septic shock in Europe, North America and Australia between 2009 and 2019- results from a systematic review and meta-analysis. Crit. Care 24 (1), 239. doi:10.1186/s13054-020-02950-2
Biesalski, H. K., Bueno de Mesquita, B., Chesson, A., Chytil, F., Grimble, R., Hermus, R. J., et al. (1998). European consensus statement on lung cancer: risk factors and prevention. Lung cancer panel. CA Cancer J. Clin. 48 (3), 167–176. doi:10.3322/canjclin.48.3.167
Bladon, S., Ashiru-Oredope, D., Cunningham, N., Pate, A., Martin, G. P., Zhong, X., et al. (2024). Rapid systematic review on risks and outcomes of sepsis: the influence of risk factors associated with health inequalities. Int. J. Equity Health 23 (1), 34. doi:10.1186/s12939-024-02114-6
Bohl, D. D., Sershon, R. A., Fillingham, Y. A., and Della Valle, C. J. (2016). Incidence, risk factors, and sources of sepsis following total joint arthroplasty. J. Arthroplasty 31 (12), 2875–2879. doi:10.1016/j.arth.2016.05.031
Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40 (4), 304–314. doi:10.1002/gepi.21965
Brion, M. J., Shakhbazov, K., and Visscher, P. M. (2013). Calculating statistical power in Mendelian randomization studies. Int. J. Epidemiol. 42 (5), 1497–1501. doi:10.1093/ije/dyt179
Bulik-Sullivan, B. K., Loh, P. R., Finucane, H. K., Ripke, S., Yang, J., Patterson, N., et al. (2015). LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47 (3), 291–295. doi:10.1038/ng.3211
Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37 (7), 658–665. doi:10.1002/gepi.21758
Burgess, S., and Thompson, S. G. (2015). Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am. J. Epidemiol. 181 (4), 251–260. doi:10.1093/aje/kwu283
Burgess, S., and Thompson, S. G. (2017). Interpreting findings from Mendelian randomization using the MR-Egger method. Eur. J. Epidemiol. 32 (5), 377–389. doi:10.1007/s10654-017-0255-x
Cheng, H., Garrick, D. J., and Fernando, R. L. (2017). Efficient strategies for leave-one-out cross validation for genomic best linear unbiased prediction. J. Anim. Sci. Biotechnol. 8, 38. doi:10.1186/s40104-017-0164-6
Guan, W., Steffen, B. T., Lemaitre, R. N., Wu, J. H. Y., Tanaka, T., Manichaikul, A., et al. (2014). Genome-wide association study of plasma N6 polyunsaturated fatty acids within the cohorts for heart and aging research in genomic epidemiology consortium. Circ. Cardiovasc Genet. 7 (3), 321–331. doi:10.1161/CIRCGENETICS.113.000208
Hemani, G., Tilling, K., and Davey Smith, G. (2017). Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 13 (11), e1007081. doi:10.1371/journal.pgen.1007081
Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base platform supports systematic causal inference across the human phenome. Elife 7, e34408. doi:10.7554/eLife.34408
Henriksen, D. P., Pottegård, A., Laursen, C. B., Jensen, T. G., Hallas, J., Pedersen, C., et al. (2015). Risk factors for hospitalization due to community-acquired sepsis - a population-based case-control study. PLoS One 10 (4), e0124838. doi:10.1371/journal.pone.0124838
Hensley, M. K., Donnelly, J. P., Carlton, E. F., and Prescott, H. C. (2019). Epidemiology and outcomes of cancer-related versus non-cancer-related sepsis hospitalizations. Crit. Care Med. 47 (10), 1310–1316. doi:10.1097/CCM.0000000000003896
Hu, J., Gan, Q., Zhou, D., Xia, X., Xiang, W., Xiao, R., et al. (2023). Evaluating the risk of sepsis attributing to obesity: a two-sample Mendelian randomization study. Postgrad. Med. J. 99 (1178), 1266–1271. doi:10.1093/postmj/qgad072
Huang, J., Deng, Y., Tin, M. S., Lok, V., Ngai, C. H., Zhang, L., et al. (2022). Distribution, risk factors, and temporal trends for lung cancer incidence and mortality: a global analysis. Chest 161 (4), 1101–1111. doi:10.1016/j.chest.2021.12.655
Jian, Z., Wang, M., Jin, X., and Wei, X. (2021). Genetically predicted higher educational attainment decreases the risk of COVID-19 susceptibility and severity: a mendelian randomization study. Front. Public Health 9, 731962. doi:10.3389/fpubh.2021.731962
Kulinskaya, E., Dollinger, M. B., and Bjørkestøl, K. (2020). On the moments of Cochran's Q statistic under the null hypothesis, with application to the meta-analysis of risk difference. Res. Synth. Methods 11 (6), 920. doi:10.1002/jrsm.1446
Lawlor, D. A., Harbord, R. M., Sterne, J. A., Timpson, N., and Davey Smith, G. (2008). Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat. Med. 27 (8), 1133–1163. doi:10.1002/sim.3034
Lee, J. J., Wedow, R., Okbay, A., Kong, E., Maghzian, O., Zacher, M., et al. (2018). Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat. Genet. 50 (8), 1112–1121. doi:10.1038/s41588-018-0147-3
Li, W., and Wang, W. (2023). Causal effects of exposure to ambient air pollution on cancer risk: insights from genetic evidence. Sci. Total Environ. 912, 168843. doi:10.1016/j.scitotenv.2023.168843
Liu, M., Jiang, Y., Wedow, R., Li, Y., Brazel, D. M., Chen, F., et al. (2019b). Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat. Genet. 51 (2), 237–244. doi:10.1038/s41588-018-0307-5
Liu, W., Wang, X., Zhao, S., Yang, S., Zheng, X., Gong, F., et al. (2024). Unraveling the immunological landscape in acute pancreatitis progression to sepsis: insights from a Mendelian randomization study on immune cell traits. Front. Immunol. 15, 1374787. doi:10.3389/fimmu.2024.1374787
Liu, Z., Mahale, P., and Engels, E. A. (2019a). Sepsis and risk of cancer among elderly adults in the United States. Clin. Infect. Dis. 68 (5), 717–724. doi:10.1093/cid/ciy530
Locham, S., Naazie, I., Canner, J., Siracuse, J., Al-Nouri, O., and Malas, M. (2021). Incidence and risk factors of sepsis in hemodialysis patients in the United States. J. Vasc. Surg. 73 (3), 1016–1021.e3. doi:10.1016/j.jvs.2020.06.126
Mahajan, A., Spracklen, C. N., Zhang, W., Ng, M. C. Y., Petty, L. E., Kitajima, H., et al. (2022). Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation. Nat. Genet. 54 (5), 560–572. doi:10.1038/s41588-022-01058-3
Malhotra, J., Malvezzi, M., Negri, E., La Vecchia, C., and Boffetta, P. (2016). Risk factors for lung cancer worldwide. Eur. Respir. J. 48 (3), 889–902. doi:10.1183/13993003.00359-2016
Mao, K., Chen, S., Chen, M., Ma, Y., Wang, Y., Huang, B., et al. (2013). Nitric oxide suppresses NLRP3 inflammasome activation and protects against LPS-induced septic shock. Cell Res. 23 (2), 201–212. doi:10.1038/cr.2013.6
McKay, J. D., Hung, R. J., Han, Y., Zong, X., Carreras-Torres, R., Christiani, D. C., et al. (2017). Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat. Genet. 49 (7), 1126–1132. doi:10.1038/ng.3892
Mirouse, A., Vigneron, C., Llitjos, J. F., Chiche, J. D., Mira, J. P., Mokart, D., et al. (2020). Sepsis and cancer: an interplay of friends and foes. Am. J. Respir. Crit. Care Med. 202 (12), 1625–1635. doi:10.1164/rccm.202004-1116TR
Moon, J. S., Lee, S., Park, M. A., Siempos, I. I., Haslip, M., Lee, P. J., et al. (2023). UCP2-induced fatty acid synthase promotes NLRP3 inflammasome activation during sepsis. J. Clin. Invest 133 (6), e169986. doi:10.1172/JCI169986
Moore, J. X., Akinyemiju, T., Bartolucci, A., Wang, H. E., Waterbor, J., and Griffin, R. (2020). A prospective study of community mediators on the risk of sepsis after cancer. J. Intensive Care Med. 35 (12), 1546–1555. doi:10.1177/0885066619881122
Pallis, A. G., and Syrigos, K. N. (2013). Lung cancer in never smokers: disease characteristics and risk factors. Crit. Rev. Oncol. Hematol. 88 (3), 494–503. doi:10.1016/j.critrevonc.2013.06.011
Pan, B., Wang, X., Kojima, S., Nishioka, C., Yokoyama, A., Honda, G., et al. (2017). The fifth epidermal growth factor like region of thrombomodulin alleviates LPS-induced sepsis through interacting with GPR15. Thromb. Haemost. 117 (3), 570–579. doi:10.1160/TH16-10-0762
Pavon, A., Binquet, C., Kara, F., Martinet, O., Ganster, F., Navellou, J. C., et al. (2013). Profile of the risk of death after septic shock in the present era: an epidemiologic study. Crit. Care Med. 41 (11), 2600–2609. doi:10.1097/CCM.0b013e31829a6e89
Rhee, C., Dantes, R., Epstein, L., Murphy, D. J., Seymour, C. W., Iwashyna, T. J., et al. (2017). Incidence and trends of sepsis in US hospitals using clinical vs claims data, 2009-2014. Jama 318 (13), 1241–1249. doi:10.1001/jama.2017.13836
Rudd, K. E., Johnson, S. C., Agesa, K. M., Shackelford, K. A., Tsoi, D., Kievlan, D. R., et al. (2020). Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the Global Burden of Disease Study. Lancet 395 (10219), 200–211. doi:10.1016/S0140-6736(19)32989-7
Said, S. A., de Hullu, J. A., van der Aa, M. A., Walraven, J. E. W., Bekkers, R. L. M., Slangen, B. F. M., et al. (2023). Impact of sepsis on the oncologic outcomes of advanced epithelial ovarian cancer patients: a multicenter observational study. Cancers (Basel) 15 (18), 4642. doi:10.3390/cancers15184642
Schuetz, P., Castro, P., and Shapiro, N. I. (2011). Diabetes and sepsis: preclinical findings and clinical relevance. Diabetes Care 34 (3), 771–778. doi:10.2337/dc10-1185
Shvetsov, Y. B., Ogino, M. H., Glibetic, N., Asato, C. B., Wilkens, L. R., Le Marchand, L., et al. (2021). Association of sepsis mortality with specific cancer sites and treatment type: the multiethnic cohort study. J. Pers. Med. 11 (2), 146. doi:10.3390/jpm11020146
Singer, M., Deutschman, C. S., Seymour, C. W., Shankar-Hari, M., Annane, D., Bauer, M., et al. (2016). The third international consensus definitions for sepsis and septic shock (Sepsis-3). Jama 315 (8), 801–810. doi:10.1001/jama.2016.0287
Skrivankova, V. W., Richmond, R. C., Woolf, B. A. R., Yarmolinsky, J., Davies, N. M., Swanson, S. A., et al. (2021). Strengthening the reporting of observational studies in epidemiology using mendelian randomization: the STROBE-MR statement. Jama 326 (16), 1614–1621. doi:10.1001/jama.2021.18236
Sudlow, C., Gallacher, J., Allen, N., Beral, V., Burton, P., Danesh, J., et al. (2015). UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 12 (3), e1001779. doi:10.1371/journal.pmed.1001779
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Teumer, A. (2018). Common methods for performing mendelian randomization. Front. Cardiovasc Med. 5, 51. doi:10.3389/fcvm.2018.00051
Thai, A. A., Solomon, B. J., Sequist, L. V., Gainor, J. F., and Heist, R. S. (2021). Lung cancer. Lancet 398 (10299), 535–554. doi:10.1016/S0140-6736(21)00312-3
Tiru, B., DiNino, E. K., Orenstein, A., Mailloux, P. T., Pesaturo, A., Gupta, A., et al. (2015). The economic and humanistic burden of severe sepsis. Pharmacoeconomics 33 (9), 925–937. doi:10.1007/s40273-015-0282-y
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50 (5), 693–698. doi:10.1038/s41588-018-0099-7
Williams, M. D., Braun, L. A., Cooper, L. M., Johnston, J., Weiss, R. V., Qualy, R. L., et al. (2004). Hospitalized cancer patients with severe sepsis: analysis of incidence, mortality, and associated costs of care. Crit. Care 8 (5), R291–R298. doi:10.1186/cc2893
Wu, F., Huang, Y., Hu, J., and Shao, Z. (2020). Mendelian randomization study of inflammatory bowel disease and bone mineral density. BMC Med. 18 (1), 312. doi:10.1186/s12916-020-01778-5
Xia, S., Zhao, Y. C., Guo, L., Gong, H., Wang, Y. K., Ma, R., et al. (2022). Do antibody-drug conjugates increase the risk of sepsis in cancer patients? A pharmacovigilance study. Front. Pharmacol. 13, 967017. doi:10.3389/fphar.2022.967017
Xu, H., Wang, X., Xu, X., Liu, L., Zhang, Y., Yan, X., et al. (2023). Association of plasma branched-chain amino acid with multiple cancers: a mendelian randomization analysis. Clin. Nutr. 42 (12), 2493–2502. doi:10.1016/j.clnu.2023.10.019
Xu, J., Zhu, G., and Zhang, H. (2024). Causal relationship between telomere length and sepsis: a bidirectional Mendelian randomization study. Sci. Rep. 14 (1), 5397. doi:10.1038/s41598-024-56205-z
Yang, R., Xiang, H., and Zheng, T. (2024). Causal associations between severe mental illness and sepsis: a Mendelian randomization study. Front. Psychiatry 15, 1341559. doi:10.3389/fpsyt.2024.1341559
Yengo, L., Sidorenko, J., Kemper, K. E., Zheng, Z., Wood, A. R., Weedon, M. N., et al. (2018). Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry. Hum. Mol. Genet. 27 (20), 3641–3649. doi:10.1093/hmg/ddy271
Zhang, L., Mukherjee, B., Ghosh, M., Gruber, S., and Moreno, V. (2008). Accounting for error due to misclassification of exposures in case-control studies of gene-environment interaction. Stat. Med. 27 (15), 2756–2783. doi:10.1002/sim.3044
Zhang, W., Lin, H., Zou, M., Yuan, Q., Huang, Z., Pan, X., et al. (2022). Nicotine in inflammatory diseases: anti-inflammatory and pro-inflammatory effects. Front. Immunol. 13, 826889. doi:10.3389/fimmu.2022.826889
Zhou, W., Kanai, M., Wu, K. H., Rasheed, H., Tsuo, K., Hirbo, J. B., et al. (2022). Global biobank meta-analysis initiative: powering genetic discovery across human disease. Cell Genom 2 (10), 100192. doi:10.1016/j.xgen.2022.100192
Zhu, H., Zhan, X., Wang, C., Deng, Y., Li, X., Song, L., et al. (2023). Causal associations between tobacco, alcohol use and risk of infectious diseases: a mendelian randomization study. Infect. Dis. Ther. 12 (3), 965–977. doi:10.1007/s40121-023-00775-4
Keywords: causal association, genome-wide association study, lung cancer, mendelian randomization, sepsis
Citation: Zhou J, Zhang Y, Yang T, Zhang K, Li A, Li M, Peng X and Chen M (2024) Causal relationships between lung cancer and sepsis: a genetic correlation and multivariate mendelian randomization analysis. Front. Genet. 15:1381303. doi: 10.3389/fgene.2024.1381303
Received: 13 February 2024; Accepted: 10 June 2024;
Published: 28 June 2024.
Edited by:
Alexandre Bureau, Laval University, CanadaReviewed by:
Akira Sugawara, Tohoku University, JapanGeorgia Damoraki, National and Kapodistrian University of Athens, Greece
Copyright © 2024 Zhou, Zhang, Yang, Zhang, Li, Li, Peng and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mingwei Chen, Y2hlbm13MzZAMTYzLmNvbQ==
†These authors have contributed equally to this work