Translational Research Based on Pathogen-Derived Molecules Inducing Innate And Adaptive Immune ResponsesView this Special Issue
Research Article | Open Access
Yanni Xu, Qiongchao Jiang, Hejun Liu, Xiaoyun Xiao, Dinghong Yang, Phei Er Saw, Baoming Luo, "DHX37 Impacts Prognosis of Hepatocellular Carcinoma and Lung Adenocarcinoma through Immune Infiltration", Journal of Immunology Research, vol. 2020, Article ID 8835393, 20 pages, 2020. https://doi.org/10.1155/2020/8835393
DHX37 Impacts Prognosis of Hepatocellular Carcinoma and Lung Adenocarcinoma through Immune Infiltration
Background. RNA helicases have various essential functions in basically all aspects of RNA metabolism, not only unwinding RNA but also disturbing the interaction of RNA with proteins. Recently, RNA helicases have been considered potential targets in cancers. So far, there has been no detailed investigation of the biological functions of RNA helicase DHX37 in cancers. Objective. We aim to identify the prognostic value of DHX37 associated with tumor microenvironments in cancers. Methods. DHX37 expression was examined via the Oncomine database and Tumor Immune Estimation Resource (TIMER). We explored the prognostic role of DHX37 in cancers across various databases. Coexpression genes, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), and fundamental regulators were performed via LinkedOmics. Confirming the prognostic value of DHX37 in liver hepatocellular carcinoma (LIHC) and lung adenocarcinoma (LUAD), we explored the role of DHX37 in infiltrated lymphocytes in cancers using the Gene Expression Profiling Interactive Analysis (GEPIA) and TIMER databases. Results. Through GO and KEGG analyses, expression of DHX37 was also correlated with complex function-specific networks involving the ribosome and RNA metabolic signaling pathways. In LIHC and LUAD, DHX37 expression showed significant positive correlations with markers of Tregs, myeloid-derived suppressor cells (MDSCs), and T cell exhaustion, contributing to immune tolerance. Conclusion. These results indicate that DHX37 can serve as a prognostic biomarker in LIHC and LUAD while having an important role in immune tolerance by activating the function of Tregs, MDSC, and T cell exhaustion.
Cancer is now known to be a disease which involves multiple players, including its relationship with micro- and macroenvironments. Tumor microenvironment (TME) is crucial in cancer progression and therapeutic responses [1, 2] and consists of complex components, among which inflammatory cells make up for the majority proportion and are valuable for diagnostic and prognostic assessment of tumors. The widely used therapy of immune checkpoint modulators, including programmed cell death protein 1 (PD-1), cytotoxic T lymphocyte-associated antigen-4 (CTLA4), lymphocyte activation gene-3 (LAG-3), and T cell immunoglobulin mucin 3 (TIM-3) [3, 4], has become a promising anticancer therapy by reversing the suppressed immune status in tumors [5, 6]. However, immunotherapy failed to have satisfactory effects in that only 20% of cancer patients benefited with significantly increased survival rates [7, 8], due to a combination of different clinical and biological behaviors .
In recent years, cumulative evidence has revealed that RNA helicases can modulate physiological processes like innate immune reactions, carcinogenic disorders, and inflammatory disorders [10–13]. For example, RNA helicase A/DHX9, as a potential therapeutic target, is associated with cancer risk and inflammation . The abnormally high expression of RNA helicases has been tested in a variety of cancers . DHX15 is significantly upregulated in HCC, and its high expression was correlated with poor prognosis. DDX3X drives posttranscriptional programming that dictates melanoma phenotype and poor disease prognosis. RNA helicases p68 and p72 show increased expression during colon carcinogenesis.
DHX37 is a highly conserved DEAH box RNA helicase essential for development of the ribosome . It conjugates to the U3 small nucleolar RNA and is significant for remodeling the U3 snoRNA-pre-18S rRNA structure during 40S maturation ; this ensures the formation of the central pseudoknot structure. The production of eukaryotic ribosomes is a massively complicated and energy-draining intracellular activity responsible for the translation of mRNAs into proteins .
Previous immunological characterization of CD8 T cells indicated that DHX37 suppressed effector functions, cytokine production, and T cell activation by modulating NF-κB , which first showed that DHX37 was related to the progression of cancers. This suggests that DHX37 may be a potential biomarker for cancers and is involved in the immune response.
However, the intrinsic mechanisms of effects of DHX37 on malignant tumor development and immune regulation have not been investigated. In this research, we conducted a comprehensive bioinformatic analysis of DHX37 expression profiles in multiple cancers. This is the first finding that reveals the association between DHX37 and Pan-Cancer in multidimensional biological functions. Our study has identified DHX37 as a potential marker of lung and liver cancers, which may guide the development of novel anticancer therapies.
2.1. Oncomine Database Analysis
We examined the gene expression of DHX37 in Pan-Cancer using the Oncomine database . The gene was assessed for differential expression with -statistics using Total Access Statistics 2002 (FMS Inc., Vienna, VA). -tests were conducted both as two-sided for differential expression analysis and one-sided for specific overexpression analysis. A value of 0.001 and a fold change of 1.5 were set as significance thresholds.
2.2. PrognoScan Database Analysis
We utilized PrognoScan database analysis  to assess the correlations between DHX37 expression and survival time in multiple cancers. Survival analysis in PrognoScan employs the minimum value approach to find the cut point in continuous gene expression measurement for grouping patients. An adjusted Cox value < 0.05 was considered statistically significant.
2.3. Kaplan-Meier Plotter Database Analysis
The correlation between the gene level of DHX37 and overall survival as well as relapse-free survival in 21 types of tumors was performed by Kaplan-Meier Plotter . The Kaplan-Meier Plotter was set up by searching the GEO, EGA, TCGA, and PubMed repositories to identify datasets with published gene expression and clinical data. The hazard ratio (HR) with 95% confidence intervals (CI) was analyzed. The correlation was significant at the 0.05 level.
2.4. UALCAN Analysis
UALCAN  was used to obtain the relative expression of DHX37 both in normal and neoplastic tissues and in several clinicopathological subgroups. UALCAN analysis use the TCGA-Assembler to download the TCGA level 3 RNA-seq data related to 31 cancer types and used TPM as the measure of expression. The value cutoff was 0.05.
Gene Expression Profiling Interactive Analysis (GEPIA) was applied to confirm the correlational analyses from TIMER . GEPIA analyzed the RNA sequencing expression data of both tumor samples and normal samples from TCGA and the Genotype Tissue Expression projects. To solve the imbalance between the tumor and normal data which can cause inefficiency in various differential analyses, GEPIA downloaded the TCGA and GTEx gene expression data that are recomputed from raw RNA-seq data by the UCSC Xena project based on a uniform pipeline. The value cutoff was 0.05.
2.6. TIMER Database Analysis
The Tumor Immune Estimation Resource aims to evaluate relative proportions of various immune cell subsets based on data from TCGA with a deconvolution approach . We assessed DHX37 expression and the correlation of DHX37 expression with the density of 6 types of infiltrating immune cells in diverse cancers by the TIMER algorithm database. In addition, the “differential expression module” was used to evaluate clinical prognosis affected by overexpressed DHX37.
2.7. TISIDB Database Analysis
Tumor and immune system interactions  were performed to explore the abundance of 28 types of tumor-infiltrating lymphocytes (TILs) and MHC markers by precalculating for 30 TCGA cancer types. The relational coefficient between DHX37 and TILs as well as MHC markers was measured by Spearman’s test.
2.8. LinkedOmics Analysis
The LinkedOmics database is a publicly available web server for analyzing multidimensional datasets based on TCGA . DHX37 coexpression was presented in heat maps with analysis of Pearson’s correlation coefficient. GO (CC (cellular component), BP (biological process), and MF (molecular function)) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways and enrichment of cellular regulators including kinase targets and miRNA targets and transcription factor target were analyzed by gene set enrichment analysis (GSEA). The rank criterion was an , a minimum number of genes of 3, and a simulation of 500, using the LIHC and LUAD compared datasets.
3.1. DHX37 Expression in Various Cancers
DHX37 gene expression was retrieved using the Oncomine database to determine differences between tumor and normal tissues over a cancer-wide range. As depicted in Figure 1(a) and Table S1, the DHX37 expression was elevated in breast (1.5, ), colorectal (3.526, ), gastric (3.866, ), kidney (3.023, ), and lung cancers (2.357, ) as well as lymphoma (2.032, ), whereas DHX37 was only observed significantly reduced in the sarcoma dataset (-1.632, ).
To validate DHX37 high expression in other databases, the RNA-seq data in TCGA from TIMER were investigated. The aberrant expression of tumor masses compared with adjacent normal tissues in Pan-Cancer is shown in Figure 1(b). DHX37 expression was significantly upregulated in cancer groups, including BLCA (bladder urothelial carcinoma), BRCA (breast invasive carcinoma), CHOL (cholangiocarcinoma), COAD (colon adenocarcinoma), ESCA (esophageal carcinoma), HNSC (head and neck squamous cell carcinoma), KIRC (head and neck squamous cell carcinoma), KIRP (kidney renal papillary cell carcinoma), LIHC (liver hepatocellular carcinoma), LUAD (lung adenocarcinoma), LUSC (lung squamous cell carcinoma), PRAD (prostate adenocarcinoma), READ (rectum adenocarcinoma), STAD (stomach adenocarcinoma), and UCEC (uterine corpus endometrial carcinoma) as compared with the surrounding tissues.
3.2. Prognostic Value of DHX37 in Various Cancers
The association between DHX37 and survival time from PrognoScan based on the Gene Expression Omnibus (GEO) database is summarized in Table S2. DHX37 expression significantly impacted prognosis in 5 types of cancers, including breast, colorectal, skin, blood, and lung cancers (Figures 2(a)–2(h)). Two cohorts (GSE31210 and GSE11117) comprising 204 samples and 41 samples, respectively, of lung cancer revealed that upregulated expression of DHX37 was related to poorer final outcome (RFS (relapse-free survival) , to 7.66, Cox ; OS , , Cox ) (Figures 2(a) and (b)). Therefore, increased DHX37predicted poor outcomes in lung cancer.
Unlike the findings from PrognoScan, we found a high expression of DHX37 reduced survival in BRCA (breast invasive carcinoma) (Figures 2(k) and 2(l)) by using the Kaplan-Meier Plotter database. Similarly, the poor outcome  in liver hepatocellular carcinoma (OS (overall survival) , to 2.42, ; RFS , to 1.73, ) and lung adenocarcinoma (OS , to 2.64, ; RFS , to 1.96, ) was shown to be relevant to higher DHX37 expression (Figures 2(i), 2(j), 2(m), and 2(n)). However, DHX37 expression has less influence on head and neck squamous cell carcinoma (Figure S1C, D). For lung squamous cell carcinoma (Figures 2(o) and 2(p)), thyroid carcinoma (Figure S1E, F), rectum adenocarcinoma (Figure S1I, J), stomach adenocarcinoma (Figure S1M, N), and uterine corpus endometrial carcinoma (Figure S1O, P), DHX37 plays a protective role in their OS but not RFS. For esophageal adenocarcinoma, DHX37 was found to have a favorable effect on relapse-free survival while it worsened overall survival (Figure S1A, B). In addition, DHX37 only had significant correlation with RFS for pancreatic ductal adenocarcinoma and ovarian cancers (Figure S1G, H, S, T).
To further verify the value of differentially expressed DHX37 in the progression of different cancers, the RNA-seq data in TCGA were also exploited to confirm the prognostic implications of DHX37 in each cancer type via GEPIA (Gene Expression Profiling Interactive Analysis). We assessed correlations between DHX37 expression and clinical outcomes in 33 types of cancer (Figure S2). DHX37 overexpression was related to poor outcomes of OS (overall survival) and DFS (disease-free survival) in ACC (adrenocortical carcinoma), LGG (brain lower grade glioma), and LIHC (liver hepatocellular carcinoma); OS in LUAD (lung adenocarcinoma), MESO (mesothelioma), and THCA (thyroid carcinoma); and DFS (disease-free survival) in SKCM (Skin Cutaneous Melanoma). These results validated the predictive value of DHX37 in particular types of cancer, such as LIHC and LUAD.
3.3. DHX37 Expression Is Associated with Advanced Clinicopathological Characteristics in LIHC and LUAD
To further reveal the potential relevance of DHX37 expression in cancers, we explored the relationship between DHX37 expression and several clinical features of LIHC and LUAD patients in TCGA cohorts by UALCAN. Subgroup analysis of several clinicopathological characteristics of 421 LIHC samples and 574 LUAD samples consistently showed significantly elevated DHX37 mRNA expression. As shown in Figure 3, the transcription level of DHX37 was significantly upregulated in LIHC and LUAD patients compared to the healthy group (with subgroup analysis based on gender, disease stages, pathological grade, cancer status, and TP53 mutation). Upregulated expression of DHX37 correlates with advanced stage and poor differentiation, particularly in older men and smokers.
Overall, these data analyses indicate that patients of LIHC and LUAD with high levels of DHX37 expression tend to have tumors with advanced clinicopathological parameters.
3.4. DHX37 Coexpression Networks in Patients with LIHC and LUAD
Given the above prognostic findings from multiple databases, we chose LIHC and LUAD as representative cancers for further research. To elucidate the biological function of DHX37, LinkedOmics was applied to analyze DHX37 coexpression genes by comparing LIHC and LUAD cohorts. We found that 3682 overlap genes were positively correlated with DHX37, whereas 2002 overlap genes were negatively correlated (Table S4). This result suggests an extensive influence of DHX37 on the transcriptome. A heat map (Figures 4(a) and 4(b)) provides details about the top 50 significant (positively and negatively correlated) genes with DHX37. DHX37 expression had a strong positive association with expression of DDX54, GCN1L1, and PUS1, which reflect changes in mRNA modifications, transcriptional regulation, and DNA repair [29–31]. In line with the fact that DHX37 suppresses the immune system by modulating NF-κB, DHX37 expression is positively correlated with the expression of PDCD11, which is known as the NF-κB binding protein (Figure 4(a)). Notably, 43/50 and 18/50 genes in LIHC and LUAD, respectively, had high HR () in the top 50 significantly positive genes. Instead, there were 12/50 and 13/50 genes with low HR () in negatively significant ones (Figure 4(c)).
The functions of DHX37 were predicted by analyzing GO and KEGG by GSEA. The most highly enriched signaling pathway was determined by their normalized enrichment score (NES). As illustrated in Table 1 and Figure S3B, the biological processes, cellular components, and molecular functions strongly associated with DHX37 were cell cycle regulation, DNA replication, mRNA processing, and respiratory activity. Interestingly, GO analysis also uncovered that MHC (major histocompatibility complex), which plays a crucial role in antigen presenting in cancers, was one of the negatively correlated categories. KEGG analysis defined enrichment in DNA replication, cell cycle, ribosome biogenesis in eukaryotes, homologous recombination, Fanconi anemia pathway, notch signaling pathway, and microRNAs in cancers, while the activities like fatty acid degradation, drug metabolism, metabolism of xenobiotics by cytochrome, chemical carcinogenesis, complement and coagulation cascades, and activity of amino acid including histidine, arginine, and tyrosine were inhibited (Figure S3A). These results reveal that the functions involving cell circle modulation, amino acid metabolism, and immune activity were highly correlated with DHX37 expression.
3.5. Regulators of DHX37 in LIHC and LUAD
To delve further into the regulators of DHX37 in LIHC and LUAD, we analyzed the kinase, miRNA, and transcription factor derived from positively correlated gene sets. The top 5 most important target networks were the kinase-target ones related mainly to the polo-like kinase 1, checkpoint kinase 2, cyclin dependent kinase 2, ATR serine/threonine kinase, and cyclin dependent kinase 1 (Table 2). The miRNA-target network was associated with (TCCGTCC) MIR-184, (TGCACGA) MIR-517, (GTGGTGA) MIR-197, (CCAGGGG) MIR-331, and (CAGCAGG) MIR-370. The transcription factor-target network was related primarily to the E2F Transcription Factor family, including E2F1DP2_01, E2F1_Q6, E2F1DP1_01, E2F1DP2_01, and E2F4DP2_01. The gene set enriched for kinase is responsible mainly for regulating stability and integrality of the genome.
LeadingEdgeNum: the number of leading-edge genes; FDR: false discovery rate from Benjamini and Hochberg from gene set enrichment analysis (GSEA); V$: the annotation found in the Molecular Signatures Database for transcription factors.
3.6. DHX37 Expression Impacts Immune Infiltration Level
DHX37 is expressed in immune cells and TILs (tumor infiltrating lymphocytes). These cells serve as independent prognostic factors of clinicopathological parameters and outcome in cancers . Hence, we analyzed the correlation between DHX37 expression and immunophenotypic characteristics in Pan-Cancer from TIMER and TISIDB (tumor and immune system interaction). The comprehensive analysis indicated that the DHX37 expression had significant correlations with tumor purity in 16 types of cancers. Moreover, DHX37 expression also significantly correlated with CD8 T cells, CD4 T cells, B cells, macrophages, neutrophils, and dendritic cells (DCs) in 16, 19, 12, 18, 16, and 14 types of cancer, respectively (Table S3 and Figure S4). Furthermore, we also found that DHX37 expression weakly to moderately negatively correlated with 28 types of TILs and MHC expression across all human heterogeneous cancers, except in KIRP (kidney renal papillary cell carcinoma), THCA (thyroid carcinoma), and LGG (brain lower grade glioma) (Figures 5(a) and 5(b)).
On the basis of the findings in immune infiltration landscapes, we further found that DHX37 expression level correlates with poorer clinical outcomes and specific immune cell infiltration in LIHC and LUAD (Figures 5(c) and 5(d)). The DHX37 expression level has significant positive correlations with infiltrating levels of CD8 T cells (, ), CD4T cells (, ), macrophages (, ), neutrophils (, ), B cells (, ), and DCs (dendritic cells) (, ) in LIHC. In LUAD, except for CD8 T cells (, ), DHX37 expression positively correlated with infiltration levels of CD4 T cells (, ) and neutrophils (, ). Although DHX37 expression has no significant correlations with tumor purity in both LIHC (, ) and LUAD (, ), these findings strongly suggest that DHX37 may participate in immune response to affect patient survival in cancers like LIHC and LUAD.
3.7. Correlation Analysis between DHX37 Expression and Immune Markers
To broaden our understanding of DHX37 crosstalk with immune signatures, we assessed the correlations between DHX37 expression and immune marker genes of CD8+ T cells, T cells (general), B cells, monocytes, tumor associated macrophages (TAMs), M1 and M2 macrophages, neutrophils, NK cells, myeloid-derived suppressor cells (MDSCs), cancer-associated fibroblasts (CAFs), and DCs in LIHC and LUAD (Table 3). We also analyzed the different functional T cells, including Th1 cells, Th2 cells, Tregs, and exhausted T cells. In TIMER, after adjustments for tumor purity, the DHX37 expression level was significantly correlated with 59 out of 72 immune cell markers in LIHC and 38 out of 72 in LUAD. In LIHC, B cells, macrophages, and various types of T cells were strongly correlated with DHX37 expression (Table 3) and less significant in LUAD. Moreover, we found that expression of DHX37 positively correlates with markers of CAFs including ACTA2 (, ), FAP (, ), PDGFR (, ), and S100A4 (, ) in LIHC.
, , , and .
The expression levels of most marker set of Tregs and exhausted T cells, such as forkhead box P3 (FOXP3), C-C chemokine receptor type 8 (CCR8), signal transducer and activator of transcription 5B (STAT5B), PD-1, CTLA4, and LAG3, have strong positive correlations with DHX37 expression in LIHC and LUAD (Figure 6). FOXP3 regulates the immune suppression and is a strong prognostic factor for distant metastases . PD-1, a widely known marker related to T lymphocyte function, showed strong positive correlation with DHX37 expression, indicating that high DHX37 expression itself may be a novel predictor for immunotherapy response. In addition, DHX37 expression showed strong negative correlations with complement and strong positive correlations with markers of myeloid-derived suppressor cells. We further evaluated the correlation between DHX37 expression and the above strong significant markers in GEPIA (Table 3). Correlation results between DHX37 and the above markers are similar to this in TIMER. These findings further confirm that DHX37 plays a vital role in cancer immune escape.
DHX37, a member of the DEAH box family of RNA helicases , plays indispensable roles in many aspects of gene expression . Differential expression and dysfunction of RNA helicases have been reported in various cancers [36–39]. Although DHX37 functions have not been extensively elucidated, it was once discovered that it suppressed T cell activation in breast cancer . As elucidated in other research, LAYN was identified as a prognostic biomarker and is highly correlated with immune infiltrates in gastric and colon cancers . Therefore, it might reasonably be assumed that DHX37 expression may influence patients’ clinical outcomes through immune infiltration. Here, we report that aberrant expression level of DHX37 correlated to prognosis in multiple cancers. Overexpression of DHX37 predicted higher rates of recurrence and shorter survival times in LIHC and LUAD. Interestingly, increased levels of DHX37 expression were related to advanced clinicopathological characteristics, indicating that DHX37 can be used as a prognostic indicator of cancer stages and metastasis. In addition, this study found that in liver hepatocellular carcinoma and lung adenocarcinoma, immune infiltration levels and various immune marker sets were correlated with the level of DHX37 expression. Therefore, this study revealed the potential usage of DHX37 expression as a novel predictive biomarker.
Notably, DHX37 expression is related to various immune infiltration levels in cancers, especially in liver cancer and lung adenocarcinoma (Tables 3 and 4 and Figures 5 and 6). Surprisingly, the expression level of DHX37 in LIHC and LUAD is not correlated with tumor purity, indicating that the expression of DHX37 in cancer cells and the tumor microenvironment is equally important. Our results suggest that there is weak to moderate positive correlations between DHX37 expression and infiltration levels of CD4+ T cells and neutrophils in LIHC and LUAD (Figures 5 and 6). Moreover, gene markers of Tregs, T cell exhaustion, and MDSCs (myeloid-derived suppressor cells) showed significantly strong positive correlation with DHX37 expression, which gives clues about DHX37 function in modulating tumor immunology in LIHC and LUAD. Furthermore, significantly negative correlations can be found between DHX37 expression and several markers of complement. Immunotherapy based on checkpoint inhibitors yields significant clinical benefit for multiple cancers. However, PD-1 inhibition meets resistance partially resulting from the accumulation of MDSCs  or Tregs . MDSCs, comprising macrophages, granulocytes, dendritic cells, and immature myeloid cells , were initially identified leading to tumor persistence and metastasis . A set of microRNAs was associated with MDSCs and resistance to immunotherapy , which was consistent with our results that show expression of DHX37 induces microRNA in cancer (Figure S3A). C3 is an important component of the complement system . Recent studies find that the complement activates and functions within cells , and that this takes effect in the induction of key metabolic pathways  and the regulation of cell death. Whether the complement fights cancer  or promotes the development of cancer  or both is yet undetermined raising the possibility that the function of complement depends on the type of cancer. A recent study  showed that C3d induced the antitumor immunological effect by increasing infiltrating CD8+ T cells, by decreasing Tregs and by suppressing expression of PD-1. Interestingly, other complement components, like anaphylatoxin C5a, contributed to cancer progression by promoting an immunosuppressive microenvironment in which MDSCs are involved [52, 53]. After adjustment for tumor purity, DHX37 in LIHC and LUAD showed negative correlations with C3 and other complement regulators but not significantly correlated with C5. Our results substantiate that DHX37 overexpression has far-reaching effects in the immune response and complement systems in cancers, which may ultimately affect patient clinical outcomes.
, , , .
To probe regulators potentially responsible for DHX37, we performed enrichment analyses of target gene sets, which illustrated that DHX37 participated primarily in the spliceosome, ribosome, DNA replication, and cell cycle (Table 2). The aberrant expression of cell cycle regulatory factors in tumor cells leads to rapid multiplication and decreased apoptosis. The E2F family, always related to the progression of liver cancer , served as the main transcription factors for DHX37 dysregulation. In research into cancer biology and molecular pathways, we also found DHX37 expression was correlated with metabolic changes including inhibiting fatty acid degradation, amino acid (arginine, histidine valine, leucine, and isoleucine), metabolism, and mitochondrial function. These findings were in line with the molecular pathways illustrated in liver cancer oncogenes . In addition, DHX37 is positively correlated with notch signaling that regulates cell proliferation, differentiation, and survival . The notch pathway is an important target for many types of solid cancers . The functional effects of DHX37 on cancer cell proliferation and survival of liver cancer and lung adenocarcinoma are probably partially modulated by notch signaling. Future studies may develop compounds targeting DHX37 for precision medicine in cancers.
However, even though we utilized online tools based on widely used bioinformatic theories from public databases, this study still had one major limitation. We only performed a bioinformatic analysis of DHX37 expression without further confirmation using in vivo/in vitro experiments. Future prospective studies focusing on these aspects in a comprehensive manner could help identify the function of DHX37 in cancers.
In conclusion, DHX37 can impact cancer prognosis by not only playing direct regulatory roles in cancer cells but also affecting the immune microenvironment. Based on multilevel evidence, DHX37 plays an oncogenic role and induces a suppressive tumor microenvironment in LIHC and LUAD. These findings for the first time offer evidence that DHX37 serves as an immunobased potential therapeutic target for cancer treatment.
|MDSCs:||Myeloid-derived suppressor cells|
|TCGA:||The Cancer Genome Atlas|
|TIMER:||Tumor Immune Estimation Resource|
|GEPIA:||Gene Expression Profiling Interactive Analysis|
|TISIDB:||Tumor and immune system interaction|
|KEGG:||Kyoto Encyclopedia of Genes and Genomes|
|GSEA:||Gene set enrichment analysis|
|FDR:||False discovery rate|
|BLCA:||Bladder urothelial carcinoma|
|BRCA:||Breast invasive carcinoma|
|CESC:||Cervical squamous cell carcinoma|
|HNSC:||Head and neck squamous cell carcinoma|
|KIRC:||Kidney renal clear cell carcinoma|
|KIRP:||Kidney renal papillary cell carcinoma|
|AML:||Acute myeloid leukemia|
|LGG:||Brain lower grade glioma|
|LIHC:||Liver hepatocellular carcinoma|
|LUSC:||Lung squamous cell carcinoma|
|OV:||Ovarian serous cystadenocarcinoma|
|PCPG:||Pheochromocytoma and paraganglioma|
|SKCM:||Skin cutaneous melanoma|
|STES:||Stomach and esophageal carcinoma|
|UCEC:||Uterine corpus endometrial carcinoma|
The datasets analyzed for this study can be found in the Oncomine, PrognoScan, GEPIA, Kaplan-Meier Plotter, UALCAN, TIMER, TISIDB, and LinkedOmics web resources.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
PS and BL designed the project and wrote and revised the manuscript. YX and QJ completed most of the data and analysis; XX, HL, XG, and DY participated in the data preparation and analysis. All authors have read and approved the final manuscript. Yanni Xu and Qiongchao Jiang contributed equally to this work.
This study was supported by the National Natural Science Foundation of China (81701715, 81873899), the Natural Science Foundation of Guangdong Province (2017A030310200, 2018A030313097), and the People’s Livelihood and Technology Planning Project of Guangzhou City of China (201803010035).
Per the publisher’s request, the details of the supplementary materials are added here. However, the editable version of each table and figure is given in a separate folder, as requested. Table S1: DHX37 expression in cancers versus normal tissue in the Oncomine database. The DHX37 expression was elevated in breast, colorectal, gastric, kidney, lung cancers as well as lymphoma, whereas DHX37 was only observed significantly reduced in the sarcoma dataset. Table S2: relationship between DHX37 expression and patient prognosis of different cancers in the PrognoScan database. Table S3: DHX37 cooccurrence genes shown in the PDF file. We found that 3682 overlap genes were positively correlated with DHX37, whereas 2002 overlap genes were negatively correlated. Table S4: summarization of correlation between DHX37 expression and immune infiltration level in diverse type cancers via the TIMER database. N: not significantly; , , , and . DHX37 expression showed significantly correlated with CD8 T cells, CD4 T cells, B cells, macrophages, neutrophils, and dendritic cells in 16, 19, 12, 18, 16, and 14 types of cancer, respectively. Figure S1: correlation of DHX37 expression with diverse types of cancer via Kaplan-Meier Plot. For esophageal adenocarcinoma, DHX37 was found to have a favorable effect on relapse-free survival while worsening overall survival. For head and neck squamous cell carcinoma, DHX37 expression has less influence. For thyroid carcinoma, rectum adenocarcinoma, stomach adenocarcinoma, and uterine corpus endometrial carcinoma, DHX37 plays a protective role in their OS but not RFS. DHX37 only had significant correlation with RFS for pancreatic ductal adenocarcinoma and ovarian cancer. Figure S2: correlation of DHX37 expression with diverse types of cancer via GEPIA. Overall survival and disease-free survival comparing the high and low expression of DHX37 in various cancers. DHX37 overexpression was related to worsening outcome of OS and DFS in ACC, LGG, and LIHC; OS in LUAD, MESO, and THCA; and DFS in SKCM. These results validated the predictive value of DHX37 in particular types of cancer, such as LIHC and LUAD. Figure S3: KEGG pathways and GO of DHX37 in the LIHC and LUAD cohort. These results reveal that the functions involving cell circle modulation, amino acid metabolism, and immune activity were highly correlated with DHX37 expression. Figure S4: correlation of the DHX37 expression with immune infiltration level in various cancers via the TIMER database. Diverse type of cancers via the TIMER database. DHX37 expression was shown as significantly correlated with CD8 T cells, CD4 T cells, B cells, macrophages, neutrophils, and dendritic cells in 16, 19, 12, 18, 16, and 14 types of cancer, respectively. (Supplementary Materials)
- D. Zeng, R. Zhou, Y. Yu et al., “Gene expression profiles for a prognostic immunoscore in gastric cancer,” The British Journal of Surgery, vol. 105, no. 10, pp. 1338–1348, 2018.
- Y. Jiang, Q. Zhang, Y. Hu et al., “Immunoscore signature: a prognostic and predictive tool in gastric cancer,” Annals of Surgery, vol. 267, no. 3, pp. 504–513, 2018.
- S. L. Topalian, C. G. Drake, and D. M. Pardoll, “Immune checkpoint blockade: a common denominator approach to cancer therapy,” Cancer Cell, vol. 27, no. 4, pp. 450–461, 2015.
- B. Huang, H. Zhang, L. Gu et al., “Advances in immunotherapy for glioblastoma multiforme,” Journal of Immunology Research, vol. 2017, Article ID 3597613, 11 pages, 2017.
- G. Bindea, B. Mlecnik, M. Tosolini et al., “Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer,” Immunity, vol. 39, no. 4, pp. 782–795, 2013.
- J. Liu, M. Xu, and Z. Yuan, “Immunoscore guided cold tumors to acquire “temperature” through integrating physicochemical and biological methods,” BIO Integration, vol. 1, no. 1, pp. 6–14, 2020.
- S. Bayraktar, S. Batoo, S. Okuno, and S. Gluck, “Immunotherapy in breast cancer,” Journal of Carcinogenesis, vol. 18, no. 1, p. 2, 2019.
- T. W. Hulett, B. A. Fox, D. J. Messenheimer et al., “Future research goals in immunotherapy,” Surgical Oncology Clinics of North America, vol. 28, no. 3, pp. 505–518, 2019.
- D. F. Quail and J. A. Joyce, “Microenvironmental regulation of tumor progression and metastasis,” Nature Medicine, vol. 19, no. 11, pp. 1423–1437, 2013.
- S. Ahmad and S. Hur, “Helicases in antiviral immunity: dual properties as sensors and effectors,” Trends in Biochemical Sciences, vol. 40, no. 10, pp. 576–585, 2015.
- H. Kato and T. Fujita, “RIG-I-like receptors and autoimmune diseases,” Current Opinion in Immunology, vol. 37, pp. 40–45, 2015.
- M. M. Rahman, E. Bagdassarian, M. A. M. Ali, and G. McFadden, “Identification of host DEAD-box RNA helicases that regulate cellular tropism of oncolytic myxoma virus in human cancer cells,” Scientific Reports, vol. 7, no. 1, pp. 1–14, 2017.
- M. Lewinsohn, A. L. Brown, L. M. Weinel et al., “Novel germ line DDX41 mutations define families with a lower age of MDS/AML onset and lymphoid malignancies,” Blood, vol. 127, no. 8, pp. 1017–1023, 2016.
- T. Lee and J. Pelletier, “The biology of DHX9 and its potential as a therapeutic target,” Oncotarget, vol. 7, no. 27, pp. 42716–42739, 2016.
- F. Raza, J. A. Waldron, and J. L. Quesne, “Translational dysregulation in cancer: eIF4A isoforms and sequence determinants of eIF4A dependence,” Biochemical Society Transactions, vol. 43, no. 6, pp. 1227–1233, 2015.
- S. Ozgur, G. Buchwald, S. Falk, S. Chakrabarti, J. R. Prabu, and E. Conti, “The conformational plasticity of eukaryotic RNA-dependent ATPases,” The FEBS Journal, vol. 282, no. 5, pp. 850–863, 2015.
- A. Roychowdhury, C. Joret, G. Bourgeois, V. Heurgue-Hamard, D. L. J. Lafontaine, and M. Graille, “The DEAH-box RNA helicase Dhr1 contains a remarkable carboxyl terminal domain essential for small ribosomal subunit biogenesis,” Nucleic Acids Research, vol. 47, no. 14, pp. 7548–7563, 2019.
- A. M. Anger, J. P. Armache, O. Berninghausen et al., “Structures of the human and Drosophila 80S ribosome,” Nature, vol. 497, no. 7447, pp. 80–85, 2013.
- M. B. Dong, G. C. Wang, R. D. Chow et al., “Systematic immunotherapy target discovery using genome-scale in vivo CRISPR screens in CD8 T cells,” Cell, vol. 178, no. 5, pp. 1189–1204.e23, 2019.
- D. R. Rhodes, S. Kalyana-Sundaram, V. Mahavisno et al., “Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles,” Neoplasia, vol. 9, no. 2, pp. 166–180, 2007.
- H. Mizuno, K. Kitada, K. Nakai, and A. Sarai, “PrognoScan: a new database for meta-analysis of the prognostic value of genes,” BMC Medical Genomics, vol. 2, no. 1, p. 18, 2009.
- A. Lánczky, A. Nagy, G. Bottai et al., “miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients,” Breast Cancer Research and Treatment, vol. 160, no. 3, pp. 439–446, 2016.
- D. S. Chandrashekar, B. Bashel, S. A. H. Balasubramanya et al., “UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses,” Neoplasia, vol. 19, no. 8, pp. 649–658, 2017.
- Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–W102, 2017.
- T. Li, J. Fan, B. Wang et al., “TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells,” Cancer Research, vol. 77, no. 21, pp. e108–e110, 2017.
- B. Ru, C. N. Wong, Y. Tong et al., “TISIDB: an integrated repository portal for tumor-immune system interactions,” Bioinformatics, vol. 35, no. 20, pp. 4200–4202, 2019.
- S. V. Vasaikar, P. Straub, J. Wang, and B. Zhang, “LinkedOmics: analyzing multi-omics data within and across 32 cancer types,” Nucleic Acids Research, vol. 46, no. D1, pp. D956–D963, 2018.
- L. J. van 't Veer, H. Dai, M. J. van de Vijver et al., “Gene expression profiling predicts clinical outcome of breast cancer,” Nature, vol. 415, no. 6871, pp. 530–536, 2002.
- M. Milek, K. Imami, N. Mukherjee et al., “DDX54 regulates transcriptome dynamics during DNA damage response,” Genome Research, vol. 27, no. 8, pp. 1344–1359, 2017.
- T. D. Cambiaghi, C. M. Pereira, R. Shanmugam et al., “Evolutionarily conserved IMPACT impairs various stress responses that require GCN1 for activating the eIF2 kinase GCN2,” Biochemical and Biophysical Research Communications, vol. 443, no. 2, pp. 592–597, 2014.
- T. M. Carlile, N. M. Martinez, C. Schaening et al., “mRNA structure determines modification by pseudouridine synthase 1,” Nature Chemical Biology, vol. 15, no. 10, pp. 966–974, 2019.
- K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, 2013.
- A. Facciabene, G. T. Motz, and G. Coukos, “T-regulatory cells: key players in tumor immune escape and angiogenesis,” Cancer Research, vol. 72, no. 9, pp. 2162–2171, 2012.
- A. Colley, J. D. Beggs, D. Tollervey, and D. L. J. Lafontaine, “Dhr1p, a putative DEAH-box RNA helicase, is associated with the box C+D snoRNP U3,” Molecular and Cellular Biology, vol. 20, no. 19, pp. 7238–7246, 2000.
- A. Luking, U. Stahl, and U. Schmidt, “The protein family of RNA helicases,” Critical Reviews in Biochemistry and Molecular Biology, vol. 33, pp. 259–296, 2008.
- K. V. Rosinski, N. Fujii, J. K. Mito et al., “DDX3Y encodes a class I MHC-restricted H-Y antigen that is expressed in leukemic stem cells,” Blood, vol. 111, no. 9, pp. 4817–4826, 2008.
- S. Lin, L. Tian, H. Shen et al., “DDX5 is a positive regulator of oncogenic NOTCH1 signaling in T cell acute lymphoblastic leukemia,” Oncogene, vol. 32, no. 40, pp. 4845–4853, 2013.
- J. Tang, H. Chen, C. C. Wong et al., “DEAD-box helicase 27 promotes colorectal cancer growth and metastasis and predicts poor survival in CRC patients,” Oncogene, vol. 37, no. 22, pp. 3006–3021, 2018.
- Y. Xue, X. Jia, C. Li et al., “DDX17 promotes hepatocellular carcinoma progression via inhibiting Klf4 transcriptional activity,” Cell Death & Disease, vol. 10, no. 11, p. 814, 2019.
- J. H. Pan, H. Zhou, L. Cooper et al., “LAYN is a prognostic biomarker and correlated with immune infiltrates in gastric and colon cancers,” Frontiers in Immunology, vol. 10, p. 6, 2019.
- J. I. Youn, S. Nagaraj, M. Collazo, and D. I. Gabrilovich, “Subsets of myeloid-derived suppressor cells in tumor-bearing mice,” Journal of Immunology, vol. 181, no. 8, pp. 5791–5802, 2008.
- M. Mkrtichyan, Y. G. Najjar, E. C. Raulfs et al., “Anti-PD-1 synergizes with cyclophosphamide to induce potent anti-tumor vaccine effects through novel mechanisms,” European Journal of Immunology, vol. 41, no. 10, pp. 2977–2986, 2011.
- D. I. Gabrilovich and S. Nagaraj, “Myeloid-derived suppressor cells as regulators of the immune system,” Nature Reviews. Immunology, vol. 9, no. 3, pp. 162–174, 2009.
- S. Nagaraj and D. I. Gabrilovich, “Tumor escape mechanism governed by myeloid-derived suppressor cells,” Cancer Research, vol. 68, no. 8, pp. 2561–2563, 2008.
- V. Huber, V. Vallacchi, V. Fleming et al., “Tumor-derived microRNAs induce myeloid suppressor cells and predict immunotherapy resistance in melanoma,” The Journal of Clinical Investigation, vol. 128, no. 12, pp. 5505–5516, 2018.
- Q. Peng, K. Li, K. Anderson et al., “Local production and activation of complement up-regulates the allostimulatory function of dendritic cells through C3a-C3aR interaction,” Blood, vol. 111, no. 4, pp. 2452–2461, 2008.
- M. K. Liszewski, M. Kolev, G. le Friec et al., “Intracellular complement activation sustains T cell homeostasis and mediates effector differentiation,” Immunity, vol. 39, no. 6, pp. 1143–1157, 2013.
- M. Kolev, S. Dimeloe, G. le Friec et al., “Complement regulates nutrient influx and metabolic reprogramming during Th1 cell responses,” Immunity, vol. 42, no. 6, pp. 1033–1047, 2015.
- R. Pio, L. Corrales, and J. D. Lambris, “The role of complement in tumor growth,” Advances in Experimental Medicine and Biology, vol. 772, pp. 229–262, 2014.
- D. Ricklin, S. K. Ricklin-Lichtsteiner, M. M. Markiewski, B. V. Geisbrecht, and J. D. Lambris, “Cutting edge: members of the Staphylococcus aureus extracellular fibrinogen-binding protein family inhibit the interaction of C3d with complement receptor 2,” Journal of Immunology, vol. 181, no. 11, pp. 7463–7467, 2008.
- J. L. Platt, I. Silva, S. J. Balin et al., “C3d regulates immune checkpoint blockade and enhances antitumor immunity,” JCI Insight, vol. 2, no. 9, 2017.
- L. Corrales, D. Ajona, S. Rafail et al., “Anaphylatoxin C5a creates a favorable microenvironment for lung cancer progression,” Journal of Immunology, vol. 189, no. 9, pp. 4674–4683, 2012.
- P. Berraondo, L. Minute, D. Ajona, L. Corrales, I. Melero, and R. Pio, “Innate immune mediators in cancer: between defense and resistance,” Immunological Reviews, vol. 274, no. 1, pp. 290–306, 2016.
- Y. L. Chen, Y. H. Uen, C. F. Li et al., “The E2F transcription factor 1 transactives stathmin 1 in hepatocellular carcinoma,” Annals of Surgical Oncology, vol. 20, no. 12, pp. 4041–4054, 2013.
- T. Couri and A. Pillai, “Goals and targets for personalized therapy for HCC,” Hepatology International, vol. 13, no. 2, pp. 125–137, 2019.
- S. Artavanis-Tsakonas and M. A. Muskavitch, “Notch: the past, the present, and the future,” Current Topics in Developmental Biology, vol. 92, pp. 1–29, 2010.
- P. Ranganathan, K. L. Weaver, and A. J. Capobianco, “Notch signalling in solid tumours: a little bit of everything but not all the time,” Nature Reviews Cancer, vol. 11, no. 5, pp. 338–351, 2011.
Copyright © 2020 Yanni Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.