Identification and Validation of a Ferroptosis-Related Long Non-coding RNA Signature for Predicting the Outcome of Lung Adenocarcinoma

BackgroundFerroptosis is a recently recognized type of programmed cell death that is involved in the biological processes of various cancers. However, the mechanism of ferroptosis in lung adenocarcinoma (LUAD) remains unclear. This study aimed to determine the role of ferroptosis-associated long non-coding RNAs (lncRNAs) in LUAD and to establish a prognostic model.MethodsWe downloaded ferroptosis-related genes from the FerrDb database and RNA sequencing data and clinicopathological characteristics from The Cancer Genome Atlas. We randomly divided the data into training and validation sets. Ferroptosis-associated lncRNA signatures with the lowest Akaike information criteria were determined using COX regression analysis and the least absolute shrinkage and selection operator. The risk scores of ferroptosis-related lncRNAs were calculated, and patients with LUAD were assigned to high- and low-risk groups based on the median risk score. The prognostic value of the risk scores was evaluated using Kaplan–Meier curves, Cox regression analyses, and nomograms. We then explored relationships between ferroptosis-related lncRNAs and the immune response using gene set enrichment analysis (GSEA).ResultsTen ferroptosis-related lncRNA signatures were identified in the training group, and Kaplan–Meier and Cox regression analyses confirmed that the risk scores were independent predictors of LUAD outcome in the training and validation sets (all P < 0.05). The area under the curve confirmed that the signatures could determine the prognosis of LUAD. The predictive accuracy of the established nomogram model was verified using the concordance index and calibration curve. The GSEA showed that the 10 ferroptosis-related lncRNAs might be associated with tumor immune response.ConclusionWe established a novel signature involving 10 ferroptosis-related lncRNAs (LINC01843, MIR193BHG, AC091185.1, AC027031.2, AL021707.2, AL031667.3, AL606834.1, AC026355.1, AC124045.1, and AC025048.4) that can accurately predict the outcome of LUAD and are associated with the immune response. This will provide new insights into the development of new therapies for LUAD.

Download Full-text

A model of twenty-three metabolic-related genes predicting overall survival for lung adenocarcinoma

PeerJ ◽

10.7717/peerj.10008 ◽

2020 ◽

Vol 8 ◽

pp. e10008

Author(s):

Zhenyu Zhao ◽

Boxue He ◽

Qidong Cai ◽

Pengfei Zhang ◽

Xiong Peng ◽

...

Keyword(s):

Prognostic Factor ◽

Lung Adenocarcinoma ◽

Cox Regression ◽

Enrichment Analysis ◽

Independent Prognostic Factor ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Risk Scores ◽

Mechanism Analysis ◽

Cox Regression Analysis

Background The highest rate of cancer-related deaths worldwide is from lung adenocarcinoma (LUAD) annually. Metabolism was associated with tumorigenesis and cancer development. Metabolic-related genes may be important biomarkers and metabolic therapeutic targets for LUAD. Materials and Methods In this study, the gleaned cohort included LUAD RNA-SEQ data from the Cancer Genome Atlas (TCGA) and corresponding clinical data (n = 445). The training cohort was utilized to model construction, and data from the Gene Expression Omnibus (GEO, GSE30219 cohort, n = 83; GEO, GSE72094, n = 393) were regarded as a testing cohort and utilized for validation. First, we used a lasso-penalized Cox regression analysis to build a new metabolic-related signature for predicting the prognosis of LUAD patients. Next, we verified the metabolic gene model by survival analysis, C-index, receiver operating characteristic (ROC) analysis. Univariate and multivariate Cox regression analyses were utilized to verify the gene signature as an independent prognostic factor. Finally, we constructed a nomogram and performed gene set enrichment analysis to facilitate subsequent clinical applications and molecular mechanism analysis. Result Patients with higher risk scores showed significantly associated with poorer survival. We also verified the signature can work as an independent prognostic factor for LUAD survival. The nomogram showed better clinical application performance for LUAD patient prognostic prediction. Finally, KEGG and GO pathways enrichment analyses suggested several especially enriched pathways, which may be helpful for us investigative the underlying mechanisms.

Download Full-text

Mining novel cell glycolysis related gene markers that can predict the survival of colon adenocarcinoma patients

Bioscience Reports ◽

10.1042/bsr20201427 ◽

2020 ◽

Vol 40 (8) ◽

Author(s):

Sihan Chen ◽

Guodong Cao ◽

Wei Wu ◽

Yida Lu ◽

Xiaobo He ◽

...

Keyword(s):

Cox Regression ◽

Mechanistic Model ◽

Single Gene ◽

Colon Adenocarcinoma ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Risk Scores ◽

Gene Markers ◽

Cox Regression Analysis

Abstract Colon adenocarcinoma (COAD) is a malignant gastrointestinal tumor, often occurring in the left colon, which is regulated by glycolysis-related processes. In past studies, multiple genes that influence the prognosis for survival have been discovered through bioinformatics analysis. However, the prediction of disease prognosis using a single gene is not an accurate method. In the present study, a mechanistic model was established to achieve better prediction for the prognosis of COAD. COAD-related data downloaded from The Cancer Genome Atlas (TCGA) were correlated with the glycolysis process using gene set enrichment analysis (GSEA) to determine the glycolysis-related genes that regulate COAD. Using COX regression analysis, glycolysis-related genes associated with the prognosis of COAD were identified, and the genes screened to establish a predictive model. The risk scores of this model were correlated with relevant clinical data to obtain a connection diagram between the model and survival rate, tumor characteristic data, etc. Finally, genes in the model were correlated with cells in the tumor microenvironment, finding that they affected specific immune cells in the model. Seven genes related to glycolysis were identified (PPARGC1A, DLAT, 6PC2, P4HA1, STC2, ANKZF1, and GPC1), which affect the prognosis of patients with COAD and constitute the model for prediction of survival of COAD patients.

Download Full-text

Upregulation of GNPNAT1 Predicts Poor Prognosis and Correlates With Immune Infiltration in Lung Adenocarcinoma

Frontiers in Molecular Biosciences ◽

10.3389/fmolb.2021.605754 ◽

2021 ◽

Vol 8 ◽

Author(s):

Wenting Liu ◽

Kaiting Jiang ◽

Jingya Wang ◽

Ting Mei ◽

Min Zhao ◽

...

Keyword(s):

Lung Adenocarcinoma ◽

Signaling Pathways ◽

Poor Prognosis ◽

Cox Regression ◽

Clinical Stage ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Clinical Value ◽

Cox Regression Analysis ◽

Normal Tissues

BackgroundGlucosamine 6-phosphate N-acetyltransferase (GNPNAT1) is a key enzyme in the hexosamine biosynthetic pathway (HBP), which functions as promoting proliferation in some tumors, yet its potential biological function and mechanism in lung adenocarcinoma (LUAD) have not been explored.MethodsThe mRNA differential expression of GNPNAT1 in LUAD and normal tissues was analyzed using the Cancer Genome Atlas (TCGA) database and validated by real-time PCR. The clinical value of GNPNAT1 in LUAD was investigated based on the data from the TCGA database. Then, immunohistochemistry (IHC) of GNPNAT1 was applied to verify the expression and clinical significance in LUAD from the protein level. The relationship between GNPNAT1 and epigenetics was explored using the cBioPortal database, and the miRNAs regulating GNPNAT1 were found using the miRNA database. The association between GNPNAT1 expression and tumor-infiltrating immune cells in LUAD was observed through the Tumor IMmune Estimation Resource (TIMER). Finally, Gene set enrichment analysis (GSEA) was used to explore the biological signaling pathways involved in GNPNAT1 in LUAD.ResultsGNPNAT1 was upregulated in LUAD compared with normal tissues, which was verified through qRT-PCR in different cell lines (P < 0.05), and associated with patients’ clinical stage, tumor size, and lymphatic metastasis status (all P < 0.01). Kaplan–Meier (KM) analysis suggested that patients with upregulated GNPNAT1 had a relatively poor prognosis (P < 0.0001). Furthermore, multivariate Cox regression analysis indicated that GNPNAT1 was an independent prognostic factor for LUAD (OS, TCGA dataset: HR = 1.028, 95% CI: 1.013–1.044, P < 0.001; OS, validation set: HR = 1.313, 95% CI: 1.130–1.526, P < 0.001). GNPNAT1 overexpression was correlated with DNA copy amplification (P < 0.0001), low DNA methylation (R = −0.52, P < 0.0001), and downregulation of hsa-miR-30d-3p (R = −0.17, P < 0.001). GNPNAT1 expression was linked to B cells (R = −0.304, P < 0.0001), CD4+T cells (R = −0.218, P < 0.0001), and dendritic cells (R = −0.137, P = 0.002). Eventually, GSEA showed that the signaling pathways of the cell cycle, ubiquitin-mediated proteolysis, mismatch repair and p53 were enriched in the GNPNAT1 overexpression group.ConclusionGNPNAT1 may be a potential prognostic biomarker and novel target for intervention in LUAD.

Download Full-text

BACH1: A Potential Predictor of Survival in Early-Stage Lung Adenocarcinoma

Oxidative Medicine and Cellular Longevity ◽

10.1155/2022/3921095 ◽

2022 ◽

Vol 2022 ◽

pp. 1-16

Author(s):

Jin Zhou ◽

Zheming Liu ◽

Huibo Zhang ◽

Tianyu Lei ◽

Jiahui Liu ◽

...

Keyword(s):

Lung Adenocarcinoma ◽

Prognostic Model ◽

Cox Regression ◽

Early Stage ◽

External Validation ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Functional Enrichment ◽

Cox Regression Analysis

Purpose. Recent researches showed the vital role of BACH1 in promoting the metastasis of lung cancer. We aimed to explore the value of BACH1 in predicting the overall survival (OS) of early-stage (stages I-II) lung adenocarcinoma. Patients and Methods. Lung adenocarcinoma cases were screened from the Cancer Genome Atlas (TCGA) database. Functional enrichment analysis was performed to obtain the biological mechanisms of BACH1. Gene set enrichment analysis (GSEA) was performed to identify the difference of biological pathways between high- and low-BACH1 groups. Univariate and multivariate COX regression analysis had been used to screen prognostic factors, which were used to establish the BACH1 expression-based prognostic model in the TCGA dataset. The C-index and time-dependent AUC curve were used to evaluate predictive power of the model. External validation of prognostic value was performed in two independent datasets from Gene Expression Omnibus (GEO). Decision analysis curve was finally used to evaluate clinical usefulness of the BACH1-based model beyond pathologic stage alone. Results. BACH1 was an independent prognostic factor for lung adenocarcinoma. High-expression BACH1 cases had worse OS. BACH1-based prognostic model showed an ideal C-index and t -AUC and validated by two GEO datasets, independently. More importantly, the BACH1-based model indicated positive clinical applicability by DCA curves. Conclusion. Our research confirmed that BACH1 was an important predictor of prognosis in early-stage lung adenocarcinoma. The higher the expression of BACH1, the worse OS of the patients.

Download Full-text

Identification of a novel glycolysis-related gene signature for predicting metastasis and survival in patients with lung adenocarcinoma

Journal of Translational Medicine ◽

10.1186/s12967-019-02173-2 ◽

2019 ◽

Vol 17 (1) ◽

Cited By ~ 19

Author(s):

Lei Zhang ◽

Zhe Zhang ◽

Zhenglun Yu

Keyword(s):

Lung Cancer ◽

Lung Adenocarcinoma ◽

Cox Regression ◽

Single Gene ◽

Gene Signature ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Test Series ◽

Cox Regression Analysis ◽

Incidence And Mortality

Abstract Background Lung cancer (LC) is one of the most lethal and most prevalent malignant tumors, and its incidence and mortality are increasing annually. Lung adenocarcinoma (LUAD) is the most common pathological type of lung cancer. Several biomarkers have been confirmed by data excavation to be related to metastasis, prognosis and survival. However, the moderate predictive effect of a single gene biomarker is not sufficient. Thus, we aimed to identify new gene signatures to better predict the possibility of LUAD. Methods Using an mRNA-mining approach, we performed mRNA expression profiling in large LUAD cohorts (n = 522) from The Cancer Genome Atlas (TCGA) database. Gene Set Enrichment Analysis (GSEA) was performed, and connections between genes and glycolysis were found in the Cox proportional regression model. Results We confirmed a set of nine genes (HMMR, B4GALT1, SLC16A3, ANGPTL4, EXT1, GPC1, RBCK1, SOD1, and AGRN) that were significantly associated with metastasis and overall survival (OS) in the test series. Based on this nine-gene signature, the patients in the test series could be divided into high-risk and low-risk groups. Additionally, multivariate Cox regression analysis revealed that the prognostic power of the nine-gene signature is independent of clinical factors. Conclusion Our study reveals a connection between the nine-gene signature and glycolysis. This research also provides novel insights into the mechanisms underlying glycolysis and offers a novel biomarker of a poor prognosis and metastasis for LUAD patients.

Download Full-text

A Comprehensive Prognostic and Immunological Analysis of a Six-Gene Signature Associated With Glycolysis and Immune Response in Uveal Melanoma

Frontiers in Immunology ◽

10.3389/fimmu.2021.738068 ◽

2021 ◽

Vol 12 ◽

Author(s):

Jun Liu ◽

Jianjun Lu ◽

Wenli Li

Keyword(s):

Immune Response ◽

Risk Stratification ◽

Cox Regression ◽

Predictive Ability ◽

Gene Signature ◽

Risk Scores ◽

Survival Prediction ◽

Regression Analyses ◽

Training Set ◽

Kaplan Meier

Uveal melanoma (UM) is a subtype of melanoma with poor prognosis. This study aimed to construct a new prognostic gene signature that can be used for survival prediction and risk stratification of UM patients. In this work, transcriptome data from the Molecular Signatures Database were used to identify the cancer hallmarks most relevant to the prognosis of UM patients. Weighted gene co-expression network, univariate least absolute contraction and selection operator (LASSO), and multivariate Cox regression analyses were used to construct the prognostic gene characteristics. Kaplan–Meier and receiver operating characteristic (ROC) curves were used to evaluate the survival predictive ability of the gene signature. The results showed that glycolysis and immune response were the main risk factors for overall survival (OS) in UM patients. Using univariate Cox regression analysis, 238 candidates related to the prognosis of UM patients were identified (p < 0.05). Using LASSO and multivariate Cox regression analyses, a six-gene signature including ARPC1B, BTBD6, GUSB, KRTCAP2, RHBDD3, and SLC39A4 was constructed. Kaplan–Meier analysis of the UM cohort in the training set showed that patients with higher risk scores had worse OS (HR = 2.61, p < 0.001). The time-dependent ROC (t-ROC) curve showed that the risk score had good predictive efficiency for UM patients in the training set (AUC > 0.9). Besides, t-ROC analysis showed that the predictive ability of risk scores was significantly higher than that of other clinicopathological characteristics. Univariate and multivariate Cox regression analyses showed that risk score was an independent risk factor for OS in UM patients. The prognostic value of risk scores was further verified in two external UM cohorts (GSE22138 and GSE84976). Two-factor survival analysis showed that UM patients with high hypoxia or immune response scores and high risk scores had the worst prognosis. Moreover, a nomogram based on the six-gene signature was established for clinical practice. In addition, risk scores were related to the immune infiltration profiles. Taken together, this study identified a new prognostic six-gene signature related to glycolysis and immune response. This six-gene signature can not only be used for survival prediction and risk stratification but also may be a potential therapeutic target for UM patients.

Download Full-text

Overexpression of HOXA10 is associated with unfavorable prognosis of acute myeloid leukemia

10.21203/rs.3.rs-16941/v1 ◽

2020 ◽

Author(s):

Chao Guo ◽

Qian-qian Ju ◽

Chun-xia Zhang ◽

Ming Gong ◽

Zhen-ling Li ◽

...

Keyword(s):

Cox Regression ◽

Cox Model ◽

Prognostic Significance ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Risk Scores ◽

High Group ◽

Cox Regression Analysis ◽

Kaplan Meier ◽

The Impact

Abstract Background HOXA family genes were crucial transcription factors involving cell proliferation and apoptosis. While few studies have focused on HOXA10 in AML. We aimed to investigate the prognostic significance of HOXA10. Methods We downloaded datasets from GEO and BeatAML database, to compare HOXA expression level between AML patients and controls. Kaplan-Meier curves were used to estimate the impact of HOXA10 expression on AML survival. The differentially expressed genes, miRNAs, lncRNAs and methylated regions between HOXA10-high and -low groups were obtained using R (version 3.6.0). Accordingly, the gene set enrichment analysis (GSEA) was accomplished using MSigDB database. Moreover, the regulatory TFs/microRNAs/lncRNAs of HOXA10 were identified. A LASSO-Cox model fitted OS to clinical and HOXA10-associated genetic variables by glmnet package. Results HOXA10 was overexpressed in AML patients than that in controls. The HOXA10-high group is significantly associated with shorter OS and DFS. A total of 1219 DEGs, 131 DEmiRs, 282 DElncRs were identified to be associated with HOXA10. GSEA revealed that 12 suppressed and 3 activated pathways in HOXA10-high group. Furthermore, the integrated regulatory network targeting HOXA10 was established. The LASSO-Cox model fitted OS to AML-survival risk scores, which included age, race, molecular risk, expression of IKZF2/LINC00649/LINC00839/FENDRR and has-miR-424-5p. The time dependent ROC indicated a satisfying AUC (1-year AUC 0.839, 3-year AUC 0.871 and 5-year AUC 0.813). Conclusions Our study identified HOXA10 overexpression as an adverse prognostic factor for AML. The LASSO-COX regression analysis revealed novel prediction model of OS with superior diagnostic utility.

Download Full-text

Identification and Validation of a Novel Immune-Related Four-lncRNA Signature for Lung Adenocarcinoma

Frontiers in Genetics ◽

10.3389/fgene.2021.639254 ◽

2021 ◽

Vol 12 ◽

Author(s):

Jixin Wang ◽

Xiangjun Yin ◽

Yin-Qiang Zhang ◽

Xuming Ji

Keyword(s):

Gene Expression ◽

Lung Adenocarcinoma ◽

Immune Cell ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Expression Omnibus ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Lasso Regression ◽

Kaplan Meier

Lung adenocarcinoma (LUAD) is a major subtype of lung cancer, the prognosis of patients with which is associated with both lncRNAs and cancer immunity. In this study, we collected gene expression data of 585 LUAD patients from The Cancer Genome Atlas (TCGA) database and 605 subjects from the Gene Expression Omnibus (GEO) database. LUAD patients were divided into high and low immune-cell-infiltrated groups according to the single sample gene set enrichment analysis (ssGSEA) algorithm to identify differentially expressed genes (DEGs). Based on the 49 immune-related DE lncRNAs, a four-lncRNA prognostic signature was constructed by applying least absolute shrinkage and selection operator (LASSO) regression, univariate Cox regression, and stepwise multivariate Cox regression in sequence. Kaplan–Meier curve, ROC analysis, and the testing GEO datasets verified the effectiveness of the signature in predicting overall survival (OS). Univariate Cox regression and multivariate Cox regression suggested that the signature was an independent prognostic factor. The correlation analysis revealed that the infiltration immune cell subtypes were related to these lncRNAs.

Download Full-text

Significance of m 6A Methyltransferase‑Related lncRNAs in the Prognosis and Immunotherapy of Hepatocellular Carcinoma

10.21203/rs.3.rs-165071/v1 ◽

2021 ◽

Author(s):

Lili Li ◽

Rongrong Xie ◽

Qichun Wei

Keyword(s):

Hepatocellular Carcinoma ◽

Correlation Analysis ◽

Immune Cells ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Immune Checkpoints ◽

Regression Analyses ◽

Cox Regression Analysis

Abstract Background: Hepatocellular carcinoma (HCC) is one of the leading causes of mortality worldwide. N6-methyladenosine (m6A) methyltransferase, has been proved to act as an oncogene in several human cancers. However, little is known about its relationship with the long non-coding RNAs (lncRNAs) that remains elusive in HCC.Methods: We comprehensively integrated gene expression data acquired from 371 HCC and 50 normal tissues in The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression & Kaplan-Meier (K-M) analysis was performed to identify m6A methyltransferase‑related lncRNAs that were related to overall survival (OS). m6A methyltransferase‑related lncRNA signature was constructed using the Least absolute shrinkage and selection operator (LASSO) Cox regression analyses. Furthermore, Cox regression analysis was applied to identify independent prognostic factors in HCC. The signature was validated in an internal validation set. Finally, the correlation analysis between gene signature and immune cells infiltration was also investigated via single-sample Gene Set Enrichment Analysis (ssGSEA) and immunotherapy response was calculated through Tumor Immune Dysfunction and Exclusion (TIDE) algorithm.Results: A total of 21 m6A methyltransferase-related lncRNAs were screened out according to Spearman correlation analysis with the immune score (|R| > 0.3, P < 0.05). We selected 3 prognostic lncRNAs to construct m6A methyltransferase-related lncRNA signature through univariate and LASSO Cox regression analyses. The univariate and multivariate Cox regression analyses demonstrated that the lncRNAs signature was a robust independent prognostic factor in OS prediction with high accuracy. The GSEA also suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes and pathways which were very well-known in the context of HCC tumorigenesis. Besides, we found that the lncRNAs signature was strikingly correlated with the tumor microenvironment (TME) immune cells infiltration and expression of critical immune checkpoints. Finally, results from the TIDE analysis revealed that the m6A methyltransferase-related lncRNAs could efficiently predict the clinical response of immunotherapy in HCC.Conclusion: Together, our study screened potential prognostic m6A methyltransferase related lncRNAs and established a novel m6A methyltransferase-based prognostic model of HCC, which not only provides new potential prognostic biomarkers and therapeutic targets but also deepens our understanding of tumor immune microenvironment status and lays a theoretical foundation for immunotherapy.

Download Full-text

Identification of m6A methyltransferase-related lncRNA signature for predicting immunotherapy and prognosis in patients with hepatocellular carcinoma

Bioscience Reports ◽

10.1042/bsr20210760 ◽

2021 ◽

Author(s):

Lili Li ◽

Rongrong Xie ◽

Guangrong Lu

Keyword(s):

Hepatocellular Carcinoma ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Immune Checkpoints ◽

Regression Analyses ◽

Protein Coding ◽

Normal Tissues ◽

Kaplan Meier

N6-methyladenosine (m6A) methyltransferase has been shown to be an oncogene in a variety of cancers. Nevertheless, the relationship between the long non-coding RNAs (lncRNAs) and hepatocellular carcinoma (HCC) remains elusive. We integrated the gene expression data of 371 HCC and 50 normal tissues from The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression & Kaplan-Meier (K-M) analysis were performed to identify m6A methyltransferase‑related lncRNAs. Three prognostic lncRNAs were selected by univariate and LASSO Cox regression analyses to construct the m6A methyltransferase-related lncRNA signature. Multivariate Cox regression analyses illustrated that this signature was an independent prognostic factor for overall survival (OS) prediction. The Gene Set Enrichment Analysis (GSEA) suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes and pathways. Besides, we discovered that the lncRNAs signature was correlated with the tumor microenvironment (TME) and the expression of critical immune checkpoints. Tumor Immune Dysfunction and Exclusion (TIDE) analysis revealed that the lncRNAs could predict the clinical response to immunotherapy. Our study had originated a prognostic signature for HCC based on the potential prognostic m6A methyltransferase-related lncRNAs. This study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.

Download Full-text