1.1 Material
Based on the transcriptomic data and clinical information of 9661 cases of tumor tissues without COVID-19 infection collected from the TGCA database from the establishment of the database to August 6, 2020, including 77 cases of adrenal cortical carcinoma, 404 cases of bladder urothelial cell carcinoma, 1085 cases of breast invasive carcinoma, 306 cases of cervical cancer, 36 cases of gallbladder cancer, 275 cases of colon cancer, 47 cases of diffuse large B-cell lymphoma, 182 cases of esophageal cancer, 163 cases of glioma, 518 cases of head and neck squamous cell carcinoma, 66 cases of renal chromophobe cell carcinoma, 522 cases of renal clear cell carcinoma, 286 cases of renal papillary cell carcinoma, 173 cases of acute myeloid leukemia, 518 cases of low-grade brain glioma, 369 cases of hepatocellular carcinoma, 483 cases of lung adenocarcinoma, 486 cases of squamous cell carcinoma of the lung, 86 cases of mesothelioma, 426 cases of ovarian cancer, 179 cases of pancreatic cancer, 182 cases of pheochromocytoma and paraganglioma, 492 cases of prostate cancer, 92 cases of rectal cancer, 262 cases of sarcoma, 461 cases of cutaneous melanoma, 408 cases of gastric cancer, 137 cases of testicular germ cell tumor, 512 cases of thyroid cancer, 118 cases of thymic carcinoma, 174 cases of endometrial carcinoma, 57 cases of uterine sarcoma and 79 cases of uveal melanoma. The effects of IL6, TNF, IL2RA and IL10 in CRS caused by COVID-19 on the prognosis of patients with above tumors were analyzed.
Using the clinical data and laboratory test results of 232 patients with confirmed tumor diagnosis and co-infected with COVID-19 collected from 9 hospitals including Wuhan Tongji Hospital, Union Hospital and Wuhan First Hospital in China from January 13 to March 18, 2020 [9], namely 33 cases of bladder cancer, 31 cases of breast cancer, 27 cases of colorectal cancer, 25 cases of thyroid cancer, 23 cases of lung cancer, 18 cases of gastric cancer, 11 cases of cervical cancer, 8 cases of liver cancer, and 8 cases of prostate cancer, 6 cases of esophageal cancer and 42 other cases, to analyze the effect of tumor-expressed inflammatory factors on the outcome of COVID-19.
1.2 Methods
TCGA Visualization data platform GEPIA[10] (Gene Expression Profiling Interactive Analysis, http://gepia.cancer-pku.cn) is an interactive web server of cancer expression profile data developed by Peking University. It contains the RNA sequencing expression data and corresponding clinical information of 33 kinds of cancer tissues from TCGA database and 32 kinds of normal tissues from GTEx database. Gene expression and its correlation with prognosis of cancer patients can be analyzed online. Select the Survival Plots, enter IL6, TNF, IL2RA, and IL10 in Gene, and then select each cancer type in Datasets Selection, setting 95% confidence intervals. The median expression of IL6, TNF, IL2RA and IL10 in each cancer species was used as the boundary for survival analysis to evaluate the impact of inflammatory factors produced by COVID-19 on the prognosis of cancer patients.
Based on the data of COVID-19 retrospective clinical studies jointly published by 9 hospitals in Wuhan, China, cancer types with no less than 6 co-infected cases of COVID-19 were selected, namely bladder cancer, breast cancer, colorectal cancer, thyroid cancer, lung cancer, gastric cancer, cervical cancer, liver cancer, prostate cancer and esophageal cancer. According to the median mRNA expression of IL6, TNF, IL2RA and IL10 in the 10 types of cancer, TCGA data and tumor patients co-infected with COVID-19 were divided into high and low expression groups, that is, The tumor tissues with low expression of IL6 included liver cancer, thyroid cancer, breast cancer, prostate cancer and cervical cancer. while the IL6 high expression group included colorectal cancer, bladder cancer, esophageal cancer, gastric cancer and lung cancer. The tumor tissues with low expression of TNF included liver cancer, prostate cancer, thyroid cancer, bladder cancer and colorectal cancer, while the TNF high expression group included gastric cancer, esophageal cancer, lung cancer, breast cancer and cervical cancer. The IL2RA low expression group included liver cancer, prostate cancer, thyroid cancer, breast cancer and bladder cancer, while the IL2RA high expression group included cervical cancer, colorectal cancer, esophageal cancer, gastric cancer and lung cancer. The IL10 low expression group included liver cancer, prostate cancer, thyroid cancer, colorectal cancer and esophageal cancer, while the IL10 high expression group included bladder cancer, gastric cancer, cervical cancer, breast cancer and lung cancer. Using the Graphpad Prism8.0 analysis software, select the "Grouped" option to create a new stacked graph, and analyze the effect of tumors with high expression of IL6, TNF, IL2RA and IL10 on the outcome of co-infection with COVID-19 by chi-square test.
TIMER2.0 (Tumor Immune Estimation Resource [11–13], http://timer.cistrome.org) is a comprehensive analysis of the Tumor and Immune cell interactions data platform, It covers the information of 32 Cancer species in TCGA database except acute myeloid leukemia and consists of three major online analysis modules: Immune Association, Cancer Exploration and Immune Estimation. Select the "Gene" module in the "Immune Association" section, enter IL6 in "Gene Expression" after entering the page, Then, six kinds of common intratumor infiltrating Immune cells, CD4 + T, CD8 + T, B cells, neutrophils, macrophages, and dendritic cells, were selected respectively at "Immune Infiltrates". After clicking "Submit", the platform automatically calculated and generated the Pearson correlation between IL6 and immune cell infiltration in each tumor tissue.
1.3 Statistical methods
Statistical analysis was performed using Graphpad Prism 8.0 and R4.0.0 software. Measurement data were expressed as mean ± standard deviation, using t test; Enumeration data were expressed as percentage (%), using chi-square test; Kaplan-Meier method was used to draw survival curve, log-rank test was used; Pearson correlation coefficient was used for IL6 Correlation analysis with immune cell infiltration; P < 0.05 was considered statistically significant.