Yuchao Ma1
Changjiang Meng
Liang Weng- 1Department of Cardiothoracic Surgery, Third Xiangya Hospital, Central South University, Changsha, China
- 2Clinical Research Center, Third Xiangya Hospital, Central South University, Changsha, China
- 3Xiangya Cancer Center, Xiangya Hospital, Central South University, Changsha, China
- 4Key Laboratory of Molecular Radiation Oncology Hunan Province, Xiangya Hospital, Central South University, Changsha, China
- 5Hunan International Science and Technology Collaboration Base of Precision Medicine for Cancer, Xiangya Hospital, Central South University, Changsha, China
- 6Institute of Gerontological Cancer Research, National Clinical Research Center for Gerontology, Changsha, China
- 7Xiangya Lung Cancer Center, Xiangya Hospital, Central South University, Changsha, China
Background: Trauma is a well-known risk factor for many disease, but the effect of trauma on respiratory disease is unclarified. In the present study, we aimed to evaluate the association between trauma and respiratory disease.
Method: Using both United Kingdom biobank and Finnish biobank genome-wide association study data (GWAS), we performed a two-sample Mendelian randomization (MR) analysis to evaluate the relationship between trauma and respiratory disease. We used four methods including inverse-variance weighted (IVW), weighted median, Maximum likelihood, and MR-Egger in this MR analysis. The IVW MR was selected as the main method. We also performed multivariable Mendelian randomization (MVMR) to simultaneously assess the independent impact of trauma exposure on respiratory disease.
Results: In the main two-sample MR analysis, trauma exposure was significantly associated with increased risk of respiratory disease (OR 1.15, 95%CI: 1.05-1.25). Besides, there was no heterogeneity and horizontal pleiotropy observed in the sensitivity analysis. After adjusting for pack years of smoking and body mass index (BMI), trauma exposure retained its association with respiratory disease (OR, 1.13, 95%CI, 1.04-1.23 adjusted by pack years of smoking; and OR, 1.11, 95%CI, 1.04-1.18 adjusted by BMI).
Conclusion: Our study discovered the association between trauma exposure and the increased risk of respiratory disease, suggesting the prevention and treatment with trauma to reduce the risk of respiratory disease.
Introduction
Trauma is a public health problem worldwide and a leading cause of mortality. According to the report from World Health Organization (WHO) (1), the number of deaths due to trauma is about 10% of the world’s deaths each year, which is 5.8 million people. Thus, the socio-economic and health burden of trauma needs to be given more attention.
Previous studies have confirmed that trauma is a risk factor for many diseases including cardiovascular disease, respiratory disease, infectious diseases, sepsis, shock, and so on (2–5). Notably, the respiratory disease caused by trauma is very severe, such as acute respiratory distress syndrome (ARDS). A study conducted by Rubenfeld Gordon D. et al. found that 7% of ARDS cases are attributed to trauma, and the mortality rate for trauma-related ARDS is assessed at 24% (6). It is generally considered that trauma can promote a series of inflammatory responses and oxidative stress, leading to diverse diseases (7).
While a few cross-sectional or cohort studies have investigated the association between trauma and respiratory disease (8), these observational studies are susceptible to a variety of confounding factors. Mendelian randomization (MR) is a genetic instrumental variable analysis, that assesses the underlying association between exposure and outcome (9), which is less susceptible to confusion and reverses causal bias (10). Therefore, in the present study, we performed the MR analysis to evaluate the effect of trauma exposure on the risk of respiratory diseases.
Method
The overview design of the current MR study is shown in Figure 1. The causal relationships between trauma exposure and diseases of the respiratory system risk were explored by a two-sample MR analysis. The valid MR analysis needs to satisfy thise three assumptions: first, the instrument Single-nucleotide polymorphisms (SNPs) are associated with the exposure (trauma exposure); second, the instrument SNPs are not associated with confounding; and third, the instrument SNP influences the outcome via the exposure.
Figure 1 An overview of the MR assumptions. SNPs, single nucleotide polymorphisms; MR, Mendelian randomization.
Instrumental variable selection
Genetic instruments for “trauma exposure” were identified from the United Kingdom Biobank (UKB) (11). The study included 98,720 subjects, with 35,269 cases and 63,451 controls. The instrument SNP was filtered based on the threshold of P < 5 × 10-8. F statistic represents the strength of the relationship between SNPs and exposures. F statistic related to the explained variance for exposure (R2), sample size (n) and number of SNPs (k) by the formula F = [(n−k−1)/k]/[R2/(1−R2)]. Generally, F>10 indicates that selected SNPs strongly predict trauma.
Pack years of smoking data were also extracted from the UKB. Pack years calculated for individuals who have ever smoked. A pack year was defined as the number of cigarettes smoked per day, divided by twenty, multiplied by the number of years of smoking: pack years=Number of cigarettes per day/20*(Age stopped smoking - Age start smoking). If the individuals who gave up smoking for more than 6 months, the figure will be adjusted, and calculated as follow: pack years = Number of cigarettes per day/20*(Age stopped smoking - Age start smoking - 0.5). The BMI summary data were extracted from a large meta-analysis of GWASs conducted by the GIANT consortium, enrolling 125 GWASs in more than 339,224 individuals (12). In this study, we limited our selection of SNPs to European ancestry (n = 322,154) (Table S1).
Outcome data sources
We extracted GWAS summary data of diseases of the respiratory system (case/control=107,261/111,531; overall=218,792) from MR-base website (http://app.mrbase.org/).
Power calculations
We used a non-centrality parameter-based approach to evaluate the power of our study on a publicly available mRnd web tool (http://cnsgenomics.com/shiny/mRnd/). For the binary outcome (diseases of the respiratory system) in this MR study, after we imputed the required parameters in mRnd (α=0.05 in this study), the minimum detectable OR was roughly estimated. The results of power calculations are shown in Table S2.
Statistical analysis
The associations between exposure (trauma exposure) and outcome (diseases of the respiratory system) were calculated with a two-sample MR analysis using the inverse variance weighted (IVW) method as the primary causal effect estimates to explore the causal associations. The results were shown as odds ratios (OR) and 95% confidence intervals (CI).
Sensitivity analysis
We performed some sensitivity analysis using other methods including MR Egger, Weighted median, and Maximum likelihood to consider the consistency across all MR methods. To test the robustness of the results, we also performed the leave-one-out analysis using the IVW method (Figure 2B). Moreover, we used the heterogeneity test and Egger regression intercept to examine the reliability of MR estimates. Considering the effects of common risk factors such as smoking and body mass index (BMI) on the respiratory system, multivariable Mendelian randomization (MVMR) was also conducted to estimate the effect of trauma after adjusting for pack years of smoking and BMI. In this study, MR-PRESSO global test were conducted to examine the reliability of the results based on instrumental SNP extracted from trauma exposure.
Figure 2 Forest plot (A), leave-one-out analysis (B), scatter plot (C), and funnel plot (D) of the causal effect of trauma exposure on diseases of the respiratory system risk. MR, Mendelian randomization; SNPs, single nucleotide polymorphisms.
P value < 0.05 was thought to be statistically significant. All analyses were performed in R (version 4.0.3) with R packages “TwoSampleMR”, “MVMR” and “MR-PRESSO”.
Result
Character of SNP for analysis
Table 1 shows the source of GWAS summary data for exposure and outcome. Each SNP extracted from trauma exposure with P-value < 5×10-8 and its F statistic and R2 were presented in Table 2. There were 5 instrument SNPs selected for trauma with F statistics ranging from 30 to 42, representing that the MR estimates calculated based on these SNPs are reliable and do not lead to bias from weak instruments.
Table 1 Characteristics of trauma exposure and diseases of the respiratory system datasets.
Table 2 Single-nucleotide polymorphisms associated with trauma exposure.
Association between trauma exposure and diseases of the respiratory system
In the main analysis, we found there was evidence for a strong causal association between trauma exposure and disease of the respiratory system (OR, 1.15, 95%CI, 1.05-1.25 for the IVW method, Table 3). The results of weighted median and Maximum likelihood were consistent with IVW (OR, 1.14, 95%CI, 1.02-1.27 for the weighted median method; OR, 1.15, 95%CI, 1.05-1.25 for the maximum likelihood method, Table 3). However, the results of MR Egger showed that there was no statistically significant between trauma exposure and diseases of the respiratory system (OR, 2.94, 95%CI, 0.82-10.47 for MR Egger method, Table 3), but the direction was consistent with the results of the main analysis. In the MVMR analysis (Table 3), trauma exposure still showed an adverse effect on respiratory system disease after adjusting for pack years of smoking (OR, 1.13, 95%CI, 1.04-1.23) and BMI (OR, 1.11, 95%CI, 1.04-1.18), respectively. According to the results of the statistical power at the given OR, the association between trauma exposure and diseases of the respiratory system was robust. The forest plot of a single SNP was shown in Figure 2A.
Table 3 Two-sample Mendelian randomization estimations showing the effect of trauma exposure on the risk of diseases of the respiratory system.
The results of heterogeneity were not observed in overall MR results (all Ph ≥ 0.05, Table 3). Furthermore, horizontal pleiotropy was not found in overall MR results (all Pegger_intercept ≥ 0.05, Table 3). The scatter plot and funnel plot were presented in Figure 2C, D, respectively, showing the consistency of the results. The above results all suggested that the results in our analysis were convincing. The detailed results for the heterogeneity test and pleiotropy test were shown in Table 3. Moreover, the global test of MR-PRESSO analysis proved that there was no horizontal pleiotropy in our analysis between trauma exposure and diseases of the respiratory system (P = 0.360), and no outliers were found (Table S3).
Discussion
In this two-sample MR analysis using large sample GWAS data, we observed that the trauma expsoure will contribute to the risk of disease of respiratory systems, which indicated the care and treatment for the respiratory systems after trauma exposure.
Trauma usually refers to injury in four components, including physical, sexual, emotional abuse and general trauma (2). Systemic responses, such as hemostasis, inflammation, endocrine, and neurological interactions, are involved in severe trauma (13). Numerous modifications to the body’s immunological, endocrine, and metabolic pathways take place after trauma. The complement system is triggered by trauma-induced tissue injury, which causes neutrophils and macrophages to become active and release inflammatory mediators such as interleukin-1, tumor necrosis factor (TNF), and platelet-activating factor (14). The excessive infiltration and activation of neutrophils in respiratory system cause tissue injury, which may be caused by neutrophil extracellular traps (NETs) (15). Clinical research has demonstrated an association between the onset of chronic respiratory disease and local and systemic neutrophil increase. Additionally, research on animals has demonstrated that airway-infiltrating neutrophils have an impoertant role in both the beginning and continuation of airway inflammation (16).
Interleukin 1 receptor 1 (IL-1R1) signaling is critical for lung fibrosis and inflammation of mice (17). Transforming growth factor β (TGFβ) is a key pro-fibrotic cytokine, which is augmented by IL-1. Elevated levels of TGFβ activate the proliferation of epithelial cells and fibroblasts and their conversion to myofibroblasts, ultimately leading to pulmonary fibrosis (18, 19). TNFα levels are elevated in a variety of inflammatory lung diseases including asthma, chronic obstructive pulmonary disease, acute lung injury, acute respiratory distress syndrome, nodal disease, and interstitial lung disease pulmonary fibrosis. TNFα induces the formation of inflammatory cells, promotes the generation of inflammatory mediators, and causes oxidative and nitrosative stress, airway hyperreactivity, and tissue remodeling, among other pathological effects (20). Patients with non-small cell lung cancer were shown to have considerably higher expression levels of IL-1 and TNFα (21). In terms of immune response pathways, these findings help to partially explain the higher risk of respiratory disease following trauma exposure.
The most frequent respiratory diseases brought on by trauma are pneumothorax and ARDS. Pneumothorax can occur in up to 50% of people who have had severe multiple chest injuries (22). Traumatic pneumothorax is the presence of air between the visceral and mural pleura as a result of penetrating and non-penetrating (blunt) traumatic events, with ongoing lung retraction from the chest wall (23). In addition to pneumothorax, studies have shown that major trauma, including blunt force trauma, penetrating injuries, and burns predispose to ARDS (24). The pathophysiology of trauma-induced ARDS is complicated by several causes, of which cellular damage, the breakdown of cellular junctions, and disruption of cellular barriers lead to increased vascular permeability and greater air-blood distances are key players. First, trauma induces direct harm to epithelial and endothelial lung cells and prompts the release of endogenous danger-associated molecular patterns (DAMP). DAMP is capable of inducing pulmonary and systemic inflammatory responses, which in turn induce epithelial and endothelial lung injury via excessive immune responses. Additionally, the accompanying pro-inflammatory environment can have a deleterious impact on cell adhesion junctions (25). Involved in this process are both the activation of the complement system pathway and the synthesis of C5a by lipopolysaccharide (26).
Using data from a large estabished sample of GWAS data, our study systematically evaluated the association between trauma exposure and the chance of acquiring the respiratory disease. To our knowledge, this study is the first study to use MR analysis to investigate the causal association between trauma and respiratory disease. Our findings imply that suitable steps can be taken to prevent the development of respiratory diseases following trauma exposure. There are several limitations to this study. First, the sample relies too heavily on individual European pedigrees, limiting generalizability to other pedigrees. Second, the trauma exposure in this study was determined according to self-reported questionnaires. Thus, it should be cautious when considering our results. Third, our study explored only whether the liability to trauma exposure is associated with disease of respiratory system, so the causal pathways leading to trauma exposure that causes disease of respiratory system could not be ruled out.
Conclusion
In conclusion, this study corroborates the data that trauma exposure is related to an increased risk of respiratory disease, highlighting the significance of respiratory care following the onset of trauma.
Data availability statement
Publicly available datasets were analyzed in this study. The datasets were derived from sources in the public domain: UK biobank (https://www.ukbiobank.ac.uk/) and MR-Base (https://www.mrbase.org/).
Ethics statement
This study is a secondary analysis conducted through existing GWAS data and UK biobank. The specific ethics and consent statements reviewed in this study can be accessed in the original publication.
Author contributions
All authors participated in the field survey and data collection. YM and LW drafted the manuscript. YM and CM analyzed the data. LW designed the study. YM and LW obtained the funding. All authors participated in the field survey and data collection, critically revised the manuscript, and gave final approval to the version submitted for publication.
Funding
National Multidisciplinary Cooperative Diagnosis and Treatment Capacity Building Project for Major Diseases (Lung Cancer, Z027002), National Natural Science Foundations of China (81974465 and 81900199), the recruitment program for Huxiang Talents (2019RS1009) received by L. And Scientific Research Project of Hunan Provincial Health Commission (202204024507) received by YM.
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/fendo.2022.1001223/full#supplementary-material
Supplementary Table 1 | Description of public GWAS data used for in MVMR.
Supplementary Table 2 | Statistical power calculation for Mendelian randomization analyses.
Supplementary Table 3 | Application of MR-PRESSO to detect horizontal pleiotropy in Mendelian randomization analysis.
Abbreviations
WHO, World Health Organization; MR, Mendelian randomization; GWAS, genome-wide association study; MVMR, multivariable Mendelian randomization; ARDS, acute respiratory distress syndrome; BMI, body mass index; IVW, inverse-variance weighted; MR-PRESSO, MR-Pleiotropy Residual Sum and Outlier; OR, odds ratio; CI, confidence intervals; IV, instrumental variable; SNP, Single-nucleotide polymorphism; LD, linkage disequilibrium; TGFβ, Transforming growth factor β; NETs, neutrophil extracellular traps; DAMP, danger-associated molecular patterns.
References
1. WHO, Department of Violence, Injury Prevention and Disability. Injuries and Violence. The Facts (Geneva:World Health Organization)(2010) 1–6.
2. Islam SJ, Hwan Kim J, Joseph E, Topel M, Baltrus P, Liu C, et al. Association between early trauma and ideal cardiovascular health among black americans: Results from the Morehouse-Emory cardiovascular (MECA) center for health equity. Circ Cardiovasc Qual Outcomes (2021) 14(9):e007904. doi: 10.1161/circoutcomes.121.007904
3. Xu W, Song Y. Biomarkers for patients with trauma associated acute respiratory distress syndrome. Mil Med Res (2017) 4:25. doi: 10.1186/s40779-017-0134-5
4. Jin H, Liu Z, Xiao Y, Fan X, Yan J, Liang H. Prediction of sepsis in trauma patients. Burns Trauma (2014) 2(3):106–13. doi: 10.4103/2321-3868.135479
5. Hardaway RM. Traumatic and septic shock alias post-trauma critical illness. Br J Surg (1998) 85(11):1473–9. doi: 10.1046/j.1365-2168.1998.00911.x
6. Rubenfeld GD, Caldwell E, Peabody E, Weaver J, Martin DP, Neff M, et al. Incidence and outcomes of acute lung injury. N Engl J Med (2005) 353(16):1685–93. doi: 10.1056/NEJMoa050333
7. Messerer DAC, Halbgebauer R, Nilsson B, Pavenstädt H, Radermacher P, Huber-Lang M. Immunopathophysiology of trauma-related acute kidney injury. Nat Rev Nephrol (2021) 17(2):91–111. doi: 10.1038/s41581-020-00344-9
8. Tignanelli CJ, Hemmila MR, Rogers MAM, Raghavendran K. Nationwide cohort study of independent risk factors for acute respiratory distress syndrome after trauma. Trauma Surg Acute Care Open (2019) 4(1):e000249. doi: 10.1136/tsaco-2018-000249
9. Burgess S, Thompson SG. Multivariable mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol (2015) 181(4):251–60. doi: 10.1093/aje/kwu283
10. Smith GD, Ebrahim S. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol (2003) 32(1):1–22. doi: 10.1093/ije/dyg070
11. Coleman JRI, Peyrot WJ, Purves KL, Davis KAS, Rayner C, Choi SW, et al. Genome-wide gene-environment analyses of major depressive disorder and reported lifetime traumatic experiences in UK biobank. Mol Psych (2020) 25(7):1430–46. doi: 10.1038/s41380-019-0546-6
12. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature (2015) 518(7538):197–206. doi: 10.1038/nature14177
13. Lord JM, Midwinter MJ, Chen YF, Belli A, Brohi K, Kovacs EJ, et al. The systemic immune response to trauma: an overview of pathophysiology and treatment. Lancet (2014) 384(9952):1455–65. doi: 10.1016/s0140-6736(14)60687-5
14. Hastings J, Krepska A, Roodenburg O. The metabolic and endocrine response to trauma. Anaesthesia Intensive Care Med (2014) 15(9):432–5. doi: 10.1016/j.mpaic.2014.06.012
15. Schioppa T, Nguyen HO, Salvi V, Maugeri N, Facchinetti F, Villetti G, et al. The PDE4 inhibitor tanimilast restrains the tissue-damaging properties of human neutrophils. Int J Mol Sci (2022) 23(9):4982. doi: 10.3390/ijms23094982
16. Ham J, Kim J, Ko YG, Kim HY. The dynamic contribution of neutrophils in the chronic respiratory diseases. Allergy Asthma Immunol Res (2022) 14(4):361–78. doi: 10.4168/aair.2022.14.4.361
17. Gasse P, Mary C, Guenon I, Noulin N, Charron S, Schnyder-Candrian S, et al. IL-1R1/MyD88 signaling and the inflammasome are essential in pulmonary inflammation and fibrosis in mice. J Clin Invest (2007) 117(12):3786–99. doi: 10.1172/jci32285
18. Liu RM. Oxidative stress, plasminogen activator inhibitor 1, and lung fibrosis. Antioxid Redox Signal (2008) 10(2):303–19. doi: 10.1089/ars.2007.1903
19. Zhang S, Fan Y, Qin L, Fang X, Zhang C, Yue J, et al. IL-1β augments TGF-β inducing epithelial-mesenchymal transition of epithelial cells and associates with poor pulmonary function improvement in neutrophilic asthmatics. Respir Res (2021) 22(1):216. doi: 10.1186/s12931-021-01808-7
20. Malaviya R, Laskin JD, Laskin DL. Anti-TNFα therapy in inflammatory lung diseases. Pharmacol Ther (2017) 180:90–8. doi: 10.1016/j.pharmthera.2017.06.008
21. Liu Y, Gao Y, Lin T. Expression of interleukin-1 (IL-1), IL-6, and tumor necrosis factor-α (TNF-α) in non-small cell lung cancer and its relationship with the occurrence and prognosis of cancer pain. Ann Palliat Med (2021) 10(12):12759–66. doi: 10.21037/apm-21-3471
22. Shah K, Tran J, Schmidt L. Traumatic pneumothorax: updates in diagnosis and management in the emergency department. Emerg Med Pract (2022) 25(5, Suppl 1):1–28.
23. Swierzy M, Helmig M, Ismail M, Rückert J, Walles T, Neudecker J. [Pneumothorax]. Zentralbl Chir (2014) 139 Suppl 1:S69–86;quizS87. doi: 10.1055/s-0034-1383029
24. Matthay MA, Zemans RL, Zimmerman GA, Arabi YM, Beitler JR, Mercat A, et al. Acute respiratory distress syndrome. Nat Rev Dis Primers (2019) 5(1):18. doi: 10.1038/s41572-019-0069-0
25. Lupu L, Palmer A, Huber-Lang M. Inflammation, thrombosis, and destruction: The three-headed Cerberus of trauma- and SARS-CoV-2-Induced ARDS. Front Immunol (2020) 11:584514. doi: 10.3389/fimmu.2020.584514
Keywords: trauma, respiratory disease, Mendelian randomization, MVMR, instrumental variable (IV)
Citation: Ma Y, Meng C and Weng L (2022) Association between trauma exposure and respiratory disease-A Mendelian randomization study. Front. Endocrinol. 13:1001223. doi: 10.3389/fendo.2022.1001223
Received: 23 July 2022; Accepted: 16 August 2022;
Published: 05 September 2022.
Edited by:
Yimin Zhu, Hunan Provincial People’s Hospital, ChinaReviewed by:
Xianning Wu, USTC Life Sciences and Medicine, ChinaQi Li, First Affiliated Hospital of Hainan Medical University, China
Copyright © 2022 Ma, Meng and Weng. 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: Liang Weng, wengliang@csu.edu.cn