- Traditional Chinese Medicine and Inflammation Regulation Research Group, School of Basic Medical Sciences, Chengdu University of Traditional Chinese Medicine, Chengdu, China
Background: Observational studies have shown that anti-tumor necrosis factor (TNF) therapy may be beneficial for patients with coronavirus disease 2019 (COVID-19). Nevertheless, because of the methodological restrictions of traditional observational studies, it is a challenge to make causal inferences. This study involved a two-sample Mendelian randomization analysis to investigate the causal link between nine TNFs and COVID-19 severity using publicly released genome-wide association study summary statistics.
Methods: Summary statistics for nine TNFs (21,758 cases) were obtained from a large-scale genome-wide association study. Correlation data between single-nucleotide polymorphisms and severe COVID-19 (18,152 cases vs. 1,145,546 controls) were collected from the COVID-19 host genetics initiative. The causal estimate was calculated by inverse variance-weighted (IVW), MR–Egger, and weighted median methods. Sensitivity tests were conducted to assess the validity of the causal relationship.
Results: Genetically predicted TNF receptor superfamily member 6 (FAS) positively correlated with the severity of COVID-19 (IVW, odds ratio = 1.10, 95% confidence interval = 1.01–1.19, p = 0.026), whereas TNF receptor superfamily member 5 (CD40) was protective against severe COVID-19 (IVW, odds ratio = 0.92, 95% confidence interval = 0.87–0.97, p = 0.002).
Conclusion: Genetic evidence from this study supports that the increased expression of FAS is associated with the risk of severe COVID-19 and that CD40 may have a potential protective effect against COVID-19.
1 Introduction
Coronavirus disease 2019 (COVID-19) is caused by the infection of the respiratory system with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a member of the betacoronavirus family. According to a survey from the Center for Systems Science and Engineering at Johns Hopkins University, as of 3 March 2023, there were over 676 million confirmed cases of COVID-19 worldwide, including over 6.88 million fatalities (CSSE, 2023). Several treatments have been investigated during the pandemic, including antiviral drugs, anti-inflammatory drugs, cellular therapy, plasma therapy, monoclonal antibodies, and intravenous immunoglobulin therapy (Moeinafshar et al., 2022; Niknam et al., 2022; Lin et al., 2023; Yang et al., 2023). However, some of them were revealed to have no therapeutic effect on COVID-19 (Abdool Karim and Devnarain, 2022; Marcec et al., 2022).
Many studies have demonstrated that an excessive inflammatory response may be a key contributor to an adverse outcome of COVID-19 (Heidari Nia et al., 2022; Zheng et al., 2022; Falck-Jones et al., 2023; Li et al., 2023) and that tumor necrosis factors (TNFs) play a vital role in initiating the inflammatory cascade response. Furthermore, observational studies have indicated that TNF therapy may be beneficial for patients with COVID-19 (Evangelatos et al., 2022; Tripathi et al., 2022). Additionally, a retrospective study by Stallmach et al. (2020) found a decline in pro-inflammatory cytokines in seven patients with severe COVID-19 but no potential immune-mediated inflammatory disease who received a single dose of infliximab 5 mg/kg within 0–3 days of admission. Only one of the seven patients, with a median age of 60 years, died. Nonetheless, whether TNFs are genetically causally associated with severe COVID-19 remains largely unclear because of possible confounding factors and reverse causality in traditional observational studies. Therefore, this study extracted large genome-wide association study (GWAS) data for TNFs for two-sample Mendelian randomization (MR) analysis to resolve this issue.
MR utilizes genetic variation to identify whether the observed correlation between risk factors and outcomes is concordant with a causal effect (Smith and Ebrahim, 2003). In general, the use of MR analysis minimizes confounding bias and avoids reverse causality. In the absence of large randomized controlled trials, MR methods are an important alternative strategy to provide credible evidence for a causal relationship between exposure and disease risk (Zuccolo and Holmes, 2017). Currently, MR analysis has an essential role in predicting disease risk factors, especially during pandemics (Baranova et al., 2022; Chen et al., 2022; Au Yeung et al., 2023). For the work described herein, MR was applied to assess the potential causal contribution of nine TNFs to severe COVID-19.
2 Methods
A two-sample MR analysis was performed to examine the causal effects of nine TNFs on COVID-19 severity using GWAS pooled statistics. This MR study was based on three important assumptions: i) genetic instruments predict the exposure of interest (P < 5E-06); ii) genetic instruments are independent of potential confounders; and iii) genetic instruments affect the outcome through risk factors only (Davies et al., 2018).
2.1 Data source
A recent study by Folkersen et al. (2020) mapped and replicated protein quantitative trait loci (pQTLs) for 90 cardiovascular proteins in more than 30,000 participants. A total of 90 proteins from 21,758 individuals within 13 European pedigree cohorts passed quality control criteria for use in GWAS meta-analysis. This study determined and duplicated 315 primary and 136 secondary pQTLs for 85 of these cyclins, providing new insights for translational research and drug development. The pQTL-based framework was designed to solve several critical tasks related to drug development, which included (a) mapping of protein regulatory pathways, (b) identification of new target candidates, (c) drug repositioning, (d) target-related safety, and (e) matching target mechanisms to patients through protein biomarkers or genetic polygenic risk scores (Folkersen et al., 2020). The findings provided a comprehensive toolbox for evaluating and developing therapeutic hypotheses and precision medicine methods for complex diseases. Supplementary Table S1 lists the details of the nine TNF-related phenotypes involved in our study. These nine TNFs are as follows: TNF-related activation-induced cytokine (TRANCE); TNF-related apoptosis-inducing ligand (TRAIL); TNF-related apoptosis-inducing ligand receptor 2 (TRAIL-R2); tumor necrosis factor ligand superfamily member 14 (TNFSF14); tumor necrosis factor receptor 1 (TNF-R1); tumor necrosis factor receptor 2 (TNF-R2); tumor necrosis factor receptor superfamily member 5 (CD40); tumor necrosis factor receptor superfamily member 6 (FAS); and CD40 ligand (CD40L). The genome-wide significance of TNF-related single-nucleotide polymorphisms (SNPs) (P < 5E-06) was used as an instrumental variable.
The COVID-19 host genetics initiative was an international genetics collaboration that intended to elucidate the genetic determinants of COVID-19 susceptibility and severity-related outcomes (The COVID-19 Host Genetics Initiative, 2020). To achieve this goal, researchers worldwide collected individual-level clinical and genetic data and carried out GWAS analyses. This study selected summary statistics from GWAS meta-analysis Round 7 involving 18,152 patients with confirmed very severe respiratory COVID-19 released on 3 April 2022 (https://www.covid19hg.org/results/r7/). Details of the GWAS dataset for severe COVID-19 included within our work are shown in Figure 1 and Supplementary Table S1.
TRANCE, TNF-related activation-induced cytokine; TRAIL, TNF-related apoptosis-inducing ligand; TRAIL-R2, TNF-related apoptosis-inducing ligand receptor 2; TNFSF14; tumor necrosis factor ligand superfamily member 14; TNF-R1, tumor necrosis factor receptor 1; TNF-R2, tumor necrosis factor receptor 2; CD40, tumor necrosis factor receptor superfamily member 5; FAS, tumor necrosis factor receptor superfamily member 6; and CD40L, CD40 ligand.
2.2 Statistical analyses
This study first harmonized effect alleles for exposures and outcomes and then used inverse variance-weighted (IVW), weighted median, and MR‒Egger analyses to determine MR estimates of exposure versus outcome. IVW is most widely used in MR studies and was considered the primary outcome of our study because of its ability to provide reliable causal estimates without directional pleiotropy (Bowden et al., 2015). Weighted median and MR‒Egger methods are utilized to complement IVW estimates because they present more reliable estimates over a wider range of situations but are less efficient. The linkage disequilibrium threshold was set to r2 < 0.01 within a clumping window of 1,000 kb.
The listed SNPs were excluded in the following cases. First, no specific requested SNPs existed during the extraction of specific SNPs from the resulting GWAS, and it was not possible to search for proxies with requested SNPs in LD from the resulting GWAS. Second, the effect of ambiguous SNPs with incongruent alleles or palindromic SNPs with ambiguous strands could not be corrected. In addition, to assess the strength of each instrumental variable, the F-statistic (F-statistic≥10 suggested that the selected SNPs were valid with sufficient strength) was calculated for each SNP (Burgess et al., 2011). R2 and F-statistics were computed as follows:
β, genetic estimate of exposure per SNP; N, sample size (number of cases plus controls); k, number of SNPs; EAF, effect allele frequency; and SE, standard error.
Sensitivity tests were conducted to ensure the stability of the results. First, Cochran’s Q test was conducted to assess heterogeneity in individual causal effects (a value of p > 0.05 was considered to indicate no heterogeneity). Then, MR‒Egger regression was used to evaluate the horizontal pleiotropy of the instruments (a value of p > 0.05 was considered to indicate no pleiotropy). If horizontal pleiotropy was found, the MR-pleiotropy residual sum and outlier methods (MR-PRESSO) were performed to correct for horizontal pleiotropy. Finally, a leave-one-out test was performed to evaluate the effect of individual variation on the observed associations.
The results are expressed as odds ratios (ORs) and their 95% confidence intervals (CIs). Statistical significance was accepted for p < 0.05. All analyses were carried out with R statistical software (version 4.2.1) with the R packages “TwoSampleMR” and “MRPRESSO.”
3 Results
3.1 Identification of genetic instruments
For each of the 9 TNFs (TRANCE, TRAIL, TRAIL-R2, TNFSF14, TNF-R1, TNF-R2, CD40, IRF, FAS, and CD40L), SNPs exceeding the genome-wide significant p-value (P < 5E-06) were selected as genetic instruments. The number of genetic instruments varied from 14 for TRAIL-R2 to 26 for TNF-R2 (Supplementary Table S2). The F-statistics of each SNP varied from 15.54 to 1864.79, indicating a low probability of weak instruments (Pierce et al., 2011).
3.2 Associations of TNFs with severe COVID-19
MR estimates for different methods to assess the causal impact of TNFs on severe COVID-19 are presented in Supplementary Table S3. As illustrated in Figure 2, genetically predicted CD40 levels were significantly negatively correlated with severe COVID-19 (IVW, OR = 0.92, 95% confidence interval (CI) = 0.87–0.97, p = 0.002). This result was also supported by MR‒Egger (OR = 0.90, 95% CI = 0.83–0.97, p = 0.012) and weighted median (OR = 0.91, 95% CI = 0.85–0.97, p = 0.002) analyses. Using the IVW approach under a multiplicative random-effects model, a positive correlation was detected between FAS (OR = 1.10, 95% CI = 1.01–1.19, p = 0.026) and confirmed very severe respiratory COVID-19, which was consistent with MR‒Egger (OR = 1.16, 95% CI = 1.01–1.33, p = 0.049) and weighted median (OR = 1.14, 95% CI = 1.03–1.26, p = 0.014) analyses. The estimated effect sizes of the SNPs on both exposures (CD40 and FAS) and outcomes (severe COVID-19) are presented in scatter plots in Figure 3.
FIGURE 3. Scatter plot of SNPs associated with CD40 levels (A) and FAS levels (B) and severe COVID-19 after outlier removal with MR-PRESSO. Analyses were performed through inverse variance-weighted, weighted median, and MR‒Egger methods. The slope of each line corresponds to the estimated MR effect per method.
Sensitivity analysis did not reveal heterogeneity (MR‒Egger, Q_P-val = 0.483; IVW, Q_P-val = 0.497) or significant horizontal pleiotropy (intercept = 0.005, se = 0.006, p = 0.388) for CD40 with severe COVID-19. Similarly, there was no clear evidence of heterogeneity (MR‒Egger, Q_P-val = 0.796; IVW, Q_P-val = 0.782) or horizontal pleiotropy (intercept = −0.009, se = 0.009, p = 0.322) between FAS and severe COVID-19 (Supplementary Table S4). Funnel plots, indicating directional horizontal pleiotropy for an outcome, are illustrated in Supplementary Figure S1, indicating that the results were robust and convincing. The results of the leave-one-out method showed that the causal relationship between CD40 and severe COVID-19 was robust, while FAS required caution (Supplementary Figure S2).
4 Discussion
Proteins circulating in the blood are critical for predicting human disease and are often used as biomarkers for clinical decision making (Folkersen et al., 2020). Previous studies have revealed the role of TNF superfamily members in cancer (Flieswasser et al., 2022; Curtis et al., 2023), autoimmune diseases (Tian et al., 2020; Chen et al., 2023), and respiratory diseases (Cho et al., 2020; Cross et al., 2023), which have been proven to be involved with severe COVID-19 (Zhang et al., 2021). Our genetic evidence suggests that FAS is positively associated with the severity of COVID-19 and that CD40 is protective against severe COVID-19. Given the global prevalence of COVID-19, this study may provide new evidence to explore biomarkers and drug targets for cases of severe disease.
The outcome of COVID-19 is thought to be associated with an excess of pro-inflammatory cytokines, leading to a “cytokine storm” and acute respiratory distress syndrome (Cron et al., 2023; Zebardast et al., 2023). FAS is a death receptor that mediates programmed apoptosis to maintain immune homeostasis. Supporting this hypothesis, patients with COVID-19 show elevated FAS expression on the T cell in plasma, which is known to be related to higher levels of caspase activation and PS exposure on the T cell and ultimately leads to T cell apoptosis, which is a feature of severe COVID-19 (Andre et al., 2022). Banerjee et al. (2022) showed that nedagliptin may significantly block TNF-induced inflammation/Fas-induced apoptotic death signaling compared to the positive control drug oseltamivir. This study may provide genetic evidence for the identification of COVID-19 biomarkers. Moreover, although the leave-one-out method showed that the causal relationship between FAS and severe cases of COVID-19 is not robust, this would not lead to a large bias in the results because FAS and severe COVID-19 are driven by potentially influential SNPs (F > 10).
The role of CD40 in COVID-19 is mainly applicable to vaccines. The mammalian innate immune response has a memory capacity called “trained immunity” (Netea et al., 2020). Murphy et al. (2023) used an adenoviral vector for the ChAdOx1 nCoV-19 vaccine (AZD1222) to observe whether it induces training immunity in humans and found that this vaccine induces metabolic reprogramming of monocytes and leads to enhanced expression of HLA-DR, CD40, and CD80 on monocytes, with many of the characteristics of trained immunity.Marlin et al. (2021) proposed a viral antigen vaccine targeting CD40-expressing antigen-presenting cells. They targeted the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein toward CD40 (αCD40.RBD) to induce significant levels of specific T and B cells with a long-term memory phenotype. This vaccine significantly enhanced protection against a new high-dose toxicity challenge compared to unvaccinated recovery animals. These findings, combined with our results, exemplify the strong potential of CD40 in vaccine development. However, more work is needed to determine the role of CD40 in vaccine development for SARS-CoV-2.
This work adopted a solid quasi-experimental method grounded in international consortia and utilized large, high-quality GWAS data. The genetic instrument for TNFs consists of multiple SNPs that are robustly associated with each TNF, thus providing a powerful genetic tool. In addition, there are several limitations. First, the probability that TNF-associated SNPs impact severe COVID-19 outcomes via causal pathways apart from TNF levels cannot be completely excluded. Second, genetic variants related to severe COVID-19 were obtained from comparisons between severe COVID-19 and general cohorts. Information on the SARS-CoV-2 infection status of control participants was limited, which might influence the MR analyses in this study (Yao et al., 2022) (https://www.covid19hg.org/results/r7/). Third, the majority of the sample involved was European, which minimized the bias in population stratification. However, we must acknowledge that Europe is made up of very diverse ethnicities, and the applicability of the results obtained so far to different ethnicities needs further study. Nevertheless, as genetic variation is quite strongly related to exposure, the bias due to sample overlap is quite small.
5 Conclusion
In conclusion, this MR study provides genetic evidence supporting that increased expression of FAS may be associated with the risk of severe COVID-19, with a possible protective effect of CD40 against COVID-19.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
All authors had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. JM, YS, and JL designed the study. HS, NL, and LZ acquired the study data. XL and CJ analyzed and interpreted the data. QF provided advice and material support for conducting the study. JM and HS wrote the first draft of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Interdisciplinary Innovation Group of Traditional Chinese Medicine—“Innovative Team for Clinical Efficacy Evaluation and Mechanism Research of Traditional Chinese Medicine for COVID-19” (ZYYCXTD-D-202201).
Acknowledgments
The authors would like to thank Folkersen et al. (2020) and the COVID-19 host genetics initiative (The COVID-19 Host Genetics Initiative, 2020) for conducting the GWASs and making the summary statistics publicly available.
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/fphar.2023.1171404/full#supplementary-material
References
Abdool Karim, S. S., and Devnarain, N. (2022). Time to stop using ineffective covid-19 drugs. N. Engl. J. Med. 387 (7), 654–655. doi:10.1056/NEJMe2209017
Andre, S., Picard, M., Cezar, R., Roux-Dalvai, F., Alleaume-Butaux, A., Soundaramourty, C., et al. (2022). T cell apoptosis characterizes severe Covid-19 disease. Cell Death Differ. 29 (8), 1486–1499. doi:10.1038/s41418-022-00936-x
Au Yeung, S. L., Wong, T. H. T., He, B., Luo, S., and Kwok, K. O. (2023). Does ACE2 mediate the detrimental effect of exposures related to COVID-19 risk: A mendelian randomization investigation. J. Med. Virol. 95 (1), e28205. doi:10.1002/jmv.28205
Banerjee, A., Kanwar, M., Das Mohapatra, P. K., Saso, L., Nicoletti, M., and Maiti, S. (2022). Nigellidine (Nigella sativa, black-cumin seed) docking to SARS CoV-2 nsp3 and host inflammatory proteins may inhibit viral replication/transcription and FAS-TNF death signal via TNFR 1/2 blocking. Nat. Prod. Res. 36 (22), 5817–5822. doi:10.1080/14786419.2021.2018430
Baranova, A., Cao, H., and Zhang, F. (2022). Severe COVID-19 increases the risk of schizophrenia. Psychiatry Res. 317, 114809. doi:10.1016/j.psychres.2022.114809
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: Effect estimation and bias detection through egger regression. Int. J. Epidemiol. 44 (2), 512–525. doi:10.1093/ije/dyv080
Burgess, S., Thompson, S. G., and Collaboration, C. C. G. (2011). Avoiding bias from weak instruments in Mendelian randomization studies. Int. J. Epidemiol. 40 (3), 755–764. doi:10.1093/ije/dyr036
Chen, L., Sun, X., Han, D., Zhong, J., Zhang, H., and Zheng, L. (2022). Visceral adipose tissue and risk of COVID-19 susceptibility, hospitalization, and severity: A mendelian randomization study. Front. Public Health 10, 1023935. doi:10.3389/fpubh.2022.1023935
Chen, X., Zhang, Y., Zhang, Y., Zhang, Y., Wang, S., Yu, Z., et al. (2023). Increased IFN-γ+ and TNF-α+ mucosal-associated invariant T cells in patients with aplastic anemia. Cytom. B Clin. Cytom. 104, 253–262. doi:10.1002/cyto.b.22115
Cho, J. W., Hong, M. H., Ha, S. J., Kim, Y. J., Cho, B. C., Lee, I., et al. (2020). Genome-wide identification of differentially methylated promoters and enhancers associated with response to anti-PD-1 therapy in non-small cell lung cancer. Exp. Mol. Med. 52 (9), 1550–1563. doi:10.1038/s12276-020-00493-8
Cron, R. Q., Goyal, G., and Chatham, W. W. (2023). Cytokine storm syndrome. Annu. Rev. Med. 74, 321–337. doi:10.1146/annurev-med-042921-112837
Cross, A. R., de Andrea, C. E., Villalba-Esparza, M., Landecho, M. F., Cerundolo, L., Weeratunga, P., et al. (2023). Spatial transcriptomic characterization of COVID-19 pneumonitis identifies immune circuits related to tissue injury. JCI Insight 8 (2), e157837. doi:10.1172/jci.insight.157837
CSSE (2023). COVID-19 dashboard by the center for systems science and engineering (CSSE) at Johns Hopkins university. [Online]. Available at: https://coronavirus.jhu.edu/map.html (Accessed March 26, 2023).
Curtis, J. R., Yamaoka, K., Chen, Y. H., Bhatt, D. L., Gunay, L. M., Sugiyama, N., et al. (2023). Malignancy risk with tofacitinib versus TNF inhibitors in rheumatoid arthritis: Results from the open-label, randomised controlled ORAL surveillance trial. Ann. Rheum. Dis. 82 (3), 331–343. doi:10.1136/ard-2022-222543
Davies, N. M., Holmes, M. V., and Davey Smith, G. (2018). Reading mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ 362, k601. doi:10.1136/bmj.k601
Evangelatos, G., Bamias, G., Kitas, G. D., Kollias, G., and Sfikakis, P. P. (2022). The second decade of anti-TNF-a therapy in clinical practice: New lessons and future directions in the COVID-19 era. Rheumatol. Int. 42 (9), 1493–1511. doi:10.1007/s00296-022-05136-x
Falck-Jones, S., Osterberg, B., and Smed-Sorensen, A. (2023). Respiratory and systemic monocytes, dendritic cells, and myeloid-derived suppressor cells in COVID-19: Implications for disease severity. J. Intern Med. 293 (2), 130–143. doi:10.1111/joim.13559
Flieswasser, T., Van den Eynde, A., Van Audenaerde, J., De Waele, J., Lardon, F., Riether, C., et al. (2022). The CD70-CD27 axis in oncology: The new kids on the block. J. Exp. Clin. Cancer Res. 41 (1), 12. doi:10.1186/s13046-021-02215-y
Folkersen, L., Gustafsson, S., Wang, Q., Hansen, D. H., Hedman, A. K., Schork, A., et al. (2020). Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals. Nat. Metab. 2 (10), 1135–1148. doi:10.1038/s42255-020-00287-2
Heidari Nia, M., Rokni, M., Mirinejad, S., Kargar, M., Rahdar, S., Sargazi, S., et al. (2022). Association of polymorphisms in tumor necrosis factors with SARS-CoV-2 infection and mortality rate: A case-control study and in silico analyses. J. Med. Virol. 94 (4), 1502–1512. doi:10.1002/jmv.27477
Li, S., Zhang, Y., Guan, Z., Ye, M., Li, H., You, M., et al. (2023). SARS-CoV-2 Z-RNA activates the ZBP1-RIPK3 pathway to promote virus-induced inflammatory responses. Cell Res. 33, 201–214. doi:10.1038/s41422-022-00775-y
Lin, Y., Yue, S., Yang, Y., Yang, S., Pan, Z., Yang, X., et al. (2023). Nasal spray of neutralizing monoclonal antibody 35B5 confers potential prophylaxis against severe acute respiratory syndrome coronavirus 2 variants of concern: A small-scale clinical trial. Clin. Infect. Dis. 76 (3), e336–e341. doi:10.1093/cid/ciac448
Marcec, R., Dodig, V. M., Radanovic, I., and Likic, R. (2022). Intravenous immunoglobulin (IVIg) therapy in hospitalised adult COVID-19 patients: A systematic review and meta-analysis. Rev. Med. Virol. 32 (6), e2397. doi:10.1002/rmv.2397
Marlin, R., Godot, V., Cardinaud, S., Galhaut, M., Coleon, S., Zurawski, S., et al. (2021). Targeting SARS-CoV-2 receptor-binding domain to cells expressing CD40 improves protection to infection in convalescent macaques. Nat. Commun. 12 (1), 5215. doi:10.1038/s41467-021-25382-0
Moeinafshar, A., Yazdanpanah, N., and Rezaei, N. (2022). Immune-based therapeutic approaches in COVID-19. Biomed. Pharmacother. 151, 113107. doi:10.1016/j.biopha.2022.113107
Murphy, D. M., Cox, D. J., Connolly, S. A., Breen, E. P., Brugman, A. A., Phelan, J. J., et al. (2023). Trained immunity is induced in humans after immunization with an adenoviral vector COVID-19 vaccine. J. Clin. Invest. 133 (2), e162581. doi:10.1172/JCI162581
Netea, M. G., Dominguez-Andres, J., Barreiro, L. B., Chavakis, T., Divangahi, M., Fuchs, E., et al. (2020). Defining trained immunity and its role in health and disease. Nat. Rev. Immunol. 20 (6), 375–388. doi:10.1038/s41577-020-0285-6
Niknam, Z., Jafari, A., Golchin, A., Danesh Pouya, F., Nemati, M., Rezaei-Tavirani, M., et al. (2022). Potential therapeutic options for COVID-19: An update on current evidence. Eur. J. Med. Res. 27 (1), 6. doi:10.1186/s40001-021-00626-3
Pierce, B. L., Ahsan, H., and Vanderweele, T. J. (2011). Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int. J. Epidemiol. 40 (3), 740–752. doi:10.1093/ije/dyq151
Smith, G. D., and Ebrahim, S. (2003). Mendelian randomization': Can genetic epidemiology contribute to understanding environmental determinants of disease? Int. J. Epidemiol. 32 (1), 1–22. doi:10.1093/ije/dyg070
Stallmach, A., Kortgen, A., Gonnert, F., Coldewey, S. M., Reuken, P., and Bauer, M. (2020). Infliximab against severe COVID-19-induced cytokine storm syndrome with organ failure-a cautionary case series. Crit. Care 24 (1), 444. doi:10.1186/s13054-020-03158-0
The COVID-19 Host Genetics Initiative (2020). The COVID-19 Host Genetics Initiative, a global initiative to elucidate the role of host genetic factors in susceptibility and severity of the SARS-CoV-2 virus pandemic. Eur. J. Hum. Genet. 28 (6), 715–718. doi:10.1038/s41431-020-0636-6
Tian, J., Zhang, B., Rui, K., and Wang, S. (2020). The role of GITR/GITRL interaction in autoimmune diseases. Front. Immunol. 11, 588682. doi:10.3389/fimmu.2020.588682
Tripathi, K., Godoy Brewer, G., Thu Nguyen, M., Singh, Y., Saleh Ismail, M., Sauk, J. S., et al. (2022). COVID-19 and outcomes in patients with inflammatory bowel disease: Systematic review and meta-analysis. Inflamm. Bowel Dis. 28 (8), 1265–1279. doi:10.1093/ibd/izab236
Yang, W., Cao, J., Cheng, H., Chen, L., Yu, M., Chen, Y., et al. (2023). Nanoformulations targeting immune cells for cancer therapy: mRNA therapeutics. Bioact. Mater 23, 438–470. doi:10.1016/j.bioactmat.2022.11.014
Yao, Y., Song, H., Zhang, F., Liu, J., Wang, D., Feng, Q., et al. (2022). Genetic predisposition to blood cell indices in relation to severe COVID-19. J. Med. Virol. 95, e28104. doi:10.1002/jmv.28104
Zebardast, A., Hasanzadeh, A., Ebrahimian Shiadeh, S. A., Tourani, M., and Yahyapour, Y. (2023). COVID-19: A trigger of autoimmune diseases. Cell Biol. Int. 47, 848–858. doi:10.1002/cbin.11997
Zhang, F., Mears, J. R., Shakib, L., Beynor, J. I., Shanaj, S., Korsunsky, I., et al. (2021). IFN-γ and TNF-α drive a CXCL10+ CCL2+ macrophage phenotype expanded in severe COVID-19 lungs and inflammatory diseases with tissue inflammation. Genome Med. 13 (1), 64. doi:10.1186/s13073-021-00881-3
Zheng, Q., Lin, R., Chen, Y., Lv, Q., Zhang, J., Zhai, J., et al. (2022). SARS-CoV-2 induces "cytokine storm" hyperinflammatory responses in RA patients through pyroptosis. Front. Immunol. 13, 1058884. doi:10.3389/fimmu.2022.1058884
Keywords: COVID-19, genome-wide association study, tumor necrosis factor, Mendelian randomization analyses, causal associations
Citation: Song H, Lei N, Zeng L, Li X, Jiang C, Feng Q, Su Y, Liu J and Mu J (2023) Mendelian randomization analysis identified tumor necrosis factor as being associated with severe COVID-19. Front. Pharmacol. 14:1171404. doi: 10.3389/fphar.2023.1171404
Received: 23 February 2023; Accepted: 02 June 2023;
Published: 16 June 2023.
Edited by:
Antonio Javier Carcas Sansuán, Universidad Autónoma de Madrid, SpainReviewed by:
Ali Razaghi, Karolinska Institutet (KI), SwedenIngrid Fricke-Galindo, Instituto Nacional de Enfermedades Respiratorias-México (INER), Mexico
Copyright © 2023 Song, Lei, Zeng, Li, Jiang, Feng, Su, Liu and Mu. 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: Jibin Liu, ZWFlYXMxMkAxNjMuY29t; Yue Su, c3V5dWVAY2R1dGNtLmVkdS5jbg==; Jie Mu, MTA0MTQxNTU2MEBxcS5jb20=
†These authors have contributed equally to this work